DOI: http://dx.doi.org/10.22267/rcia.173402.70

Relationship between physical properties and the magnetic susceptibility in two soils of Valle del Cauca

Relación entre propiedades físicas y la susceptibilidad magnética en dos suelos del Valle del Cauca

Cristian Jiménez A.¹; Jimmy Jhony Benavides B.²; Daniel I. Ospina-Salazar³; Orlando Zúñiga E.⁴; Oscar Ochoa B.⁵; Carlos Mosquera G.⁶

Received: December 12 2016.				Accepted: November 1 2017.

ABSTRACT

Magnetic susceptibility (MS) is a property that determines the degree of magnetization of a material according to its composition; therefore, it has potential in the assessment of agricultural soils. This study aimed to the application of this attribute to the physical analysis of soils, by determining its correlation with some physical properties in soils of Valle del Cauca. Samples were taken in two lots of sugar cane (Chondular and Santa Rosa), of 55 and 98 hectares, respectively. The lots were analyzed by descriptive statistics and spatial and Pearson correlation between MS and the physical properties of the soil, through a geographic information system software. High spatial correlations were found between MS and the properties analyzed, particularly sand and clay content (0.9 and -0.88, respectively, P<0.001) in Chondular lot, although correlations above 50% were also found with field capacity, thermal conductivity, bulk density and total pore space. Moreover, the correlations of MS were lower in Santa Rosa lot, being the most relevant in the order of 0.59 and -0.65 (P<0.005) also for sand and clay content respectively. These correlation differences were attributed to alterations due to soil tillage in Santa Rosa lot. The determination of MS and its spatial correlation with some physical properties of agricultural soils is a technique that could simplify its characterization, particularly the proportion of sand and clay.

Keywords: geospatial data, paramagnetism, pedotransfer functions, soil texture.

RESUMEN

La susceptibilidad magnética (SM) es una propiedad que determina el grado de magnetización de un material de acuerdo con su composición, dado su potencial en el diagnóstico de suelos agrícolas. Con el objetivo de evaluar la aplicabilidad de este atributo en el análisis físico de suelos, se determinó su correlación con algunas propiedades físicas en suelos del Valle del Cauca. Se realizaron muestreos en dos lotes de caña de azúcar (Chondular y Santa Rosa), de 55 y 98 hectáreas, respectivamente. Los datos se analizaron por estadísticos descriptivos, así como por matriz de correlación espacial y Pearson entre SM y las propiedades físicas del suelo, a través de un software de sistemas de información geográfica. Se encontraron correlaciones espaciales altas entre SM y las propiedades analizadas, particularmente con arenas y arcillas (0,9 y -0,88, respectivamente, P<0,001) en el lote Chondular, aunque también se encontraron correlaciones superiores al 50% con capacidad de campo, conductividad térmica, densidad aparente y espacio poroso total. Por otra parte, las correlaciones de SM fueron más bajas en el lote Santa Rosa, siendo las más relevantes del orden de 0,59 y -0,65 (P<0,005) para arenas y arcillas en su orden. Estas diferencias de correlación se atribuyeron a alteraciones por labranza del suelo en el lote Santa Rosa. De acuerdo con lo anterior, la determinación de SM y su correlación espacial con algunas propiedades físicas de suelos agrícolas, es una técnica que podría simplificar su caracterización, particularmente la proporción de arenas y arcillas.

Palabras clave: datos geoespaciales, funciones de edafotransferencia, paramagnetismo, textura del suelo..

INTRODUCTION

The survey of soils plays a key role in crop production, and mapping the spatial variability of soil properties is important for planning sustainable agricultural practices (Barbieri et al., 2008). However, due to the high costs and time required for soil sampling, testing and analysis, its enforcement is not always spread among farmers. In order to make more efficient this practice, one alternative is to study physical attributes of matter (e.g. conductivity, refraction index, radiation emission, etc.) that can be correlated with soil properties such as texture, bulk density or hydraulic conductivity. These alternative tools for indirect estimation in soil survey have been called pedotransfer functions. According to McBratney et al. (2002) pedotransfer functions can be defined as mathematical models used to estimate the soil properties from other parameters measured with greater ease and lower cost.

Some authors have suggested that there is great potential for using pedotransfer functions in countries where there is limited information or paucity of data from soils, in combination with remote sensing tools (Patil and Singh, 2016). Therefore, in situ and ex situ measurement of intensive physical properties, such as paramagnetism, is an option from which other soil attributes can be estimated.

Paramagnetism is a property of non-ferromagnetic elements (as aluminum and titanium) which is manifested in a transitory rearrangement of the electron spins in the presence of a magnetic field, which makes them attracted by it. In terms of magnetic susceptibility (MS), that is, the degree of magnetization of a given material under the influence of a magnetic field, ferromagnetic materials (iron, nickel, cobalt) exhibits very high values for this feature, in comparison with the paramagnetic ones, such as the above mentioned elements (Zúñiga et al., 2016). Magnetic susceptibility is the result of the translation and rotation of the electrons, which constitute the atoms of the minerals that are present in soils, sediments and rocks (Luque, 2008). Kanu et al. (2014) states that magnetic minerals present in soils may either be obtained from the parent rocks (lithogenic origin) during pedogenesis or because of anthropogenic activities, concluding that the content of magnetic minerals in the soil can broadly be expressed by its MS.

Numerous studies pointed out that the measurement of magnetic susceptibility can be very useful in the identification of properties related with factors and processes of soil formation, applied to soil mapping (Ramos et al., 2017). Magnetic susceptibility became popular in soil-related researches including paleoclimatic field studies (Huang et al., 2006; Torrent et al., 2010; Urcia et al., 2012), soil erosion (Royall, 2001; Rahimi et al., 2013; Jordanova et al., 2014), and soil contamination by heavy metals (Pedroso, 2013). In addition, this technique has been tested to assist in the indirect quantification of soil properties that have magnetic expression in agricultural surveys (Mathé and Lévêque, 2003; Blundell et al., 2009; Cortez et al., 2011; Marques et al., 2014). Whereas sampling and laboratory analysis are some of the more expensive steps and can hinder variability and application studies, magnetic susceptibility data can be used as a pedotransfer function to estimate clay and sand content with good accuracy (Siqueira et al., 2010). Studies such as those of Souza et al. (2010) and Torrent et al. (2007) have confirmed predictive potential of MS in relation to other soil attributes (physical, chemical and mineralogical) at different soil types.

Regarding the soil texture, for example, it has been observed that MS is associated with the nature of clay minerals and parent material (Camargo et al., 2014), and that the magnetic properties of a soil along its profile are associated with the clay content and organic carbon in the different horizons (Bautista et al., 2013). Marques et al. (2014) showed that the variability in soil properties can be represented by the spatial distribution of the MS, because this property (along with diffuse reflectance) express mineralogy, mainly clay, which in turn exhibits spatial variability correlated with spatial variability of soil attributes. Grimley et al. (2004) concluded that differences in critical MS values among their study site at Illinois (USA) seem to be controlled principally by surface soil texture and ambient pH. Dos Reis Barrios et al. (2017) studied the MS as an indicator of soil quality in sugarcane fields in southeastern Brazil, concluding that this feature can be used to identify changes in physical, chemical and mineralogical traits of Latosols under sugarcane harvested with and without previous burning, indicating the changes in soil quality. Ramos et al. (2017) reported that the limits indicated by a spatial boundary analysis based on MS values were effective to outline representative compartments of pedogenetic environments on the study area. Thus, the use of MS enables soil characterization to replace costly and time-consuming chemical and mineralogical soil analysis (Camargo et al., 2014).

The MS measurement for soil survey, complemented with computational tools such as software of geographic information, allow an alternative mapping of the study area in order to undertake corrective actions, optimize the tillage or spot any issue (e.g., slow hydraulic conductance associated with clay content). Hence, the aim of this work was to assess the spatial relationship of MS with several physical properties in two soils of southern Valle del Cauca region that could lead to an easier determination of these properties by correlation with MS.

MATERIALS AND METHODS

Location characteristics. The experiment was conducted in the south of Valle del Cauca department, Colombia, in the municipalities of Palmira and Candelaria. Measurements were done in two soils under cultivation of sugar cane, located specifically in two lots: Santa Rosa and Chondular. The Santa Rosa lot is located in Palmira, with geographical coordinates 3º35’05’’N-76º15’16’’W, altitude 1060 m.a.s.l., and average temperature of 22.8ºC. The Chondular lot is located in Candelaria, with geographical coordinates 3º26’35’’N-76º22’45’’W, at an altitude of 967 m.a.s.l., and average temperature of 23.2ºC (Figura 1). Total annual rainfall and relative humidity averages 1053.7mm and 81.9% in Santa Rosa, and 1023.1mm and 79.5% in Chondular, respectively, being both sites in the tropical dry forest ecological zone.

Figure 1. Location areas: Santa Rosa lot located in Palmira Municipality and Chondular lot located in Candelaria Municipality, in Valle del Cauca.

According to the study of soils conducted by IGAC - CVC (2004), in Santa Rosa lot dominates the Consociation of Pachic Haplustolls family, isohyperthermic, with fine texture, loamy PL symbol belonging to the Palmira complex soils. While in the Chondular lot, predominates the Consociation of Cumulic Haplustolls family, isohyperthermic, with fine texture, loamy MN symbol belonging to the Manuelita complex soils.

Soil sampling. Altered and unaltered soil samples were collected from 0 to 0.30m for physical analysis. Each sample was geo-referenced with a Garmin GPS, comprising an area of approximately 98.3 hectares in the lot Santa Rosa (subdivided in three irregular lots), and 55.3 hectares in the lot Chondular (subdivided in four irregular plots). Forty-two and 28 reference-sampling points were distributed for Santa Rosa and Chondular lots, respectively, with a density of one point per two hectares, as it is shown in (Figura 2).

Figure 2. Distribution of sampling points in Santa Rosa (left) and Chondular lots (right).

Disturbed samples were taken with a bore, extracting four sub-samples within a radius of ten meters around the georeferenced point for the representative sample mixture. For the undisturbed samples three replicates per sampling point were taken, using a core auger-hole. The samples were processed for physical properties and determination of MS in the laboratory of Environmental Physics at Universidad del Valle in Cali - Colombia. Physical analysis included the determination of sand (S), clay (C) and silt (M), bulk density (BD), field capacity moisture (FC), drainable porosity (DP), total pore space (TPS), and thermal conductivity (TC) (Table 1).

Table 1. Physical properties of the soils analyzed.

Magnetic susceptibility determination. For the determination of MS, it was used a paramagnetic soil meter (PCSM Pike Agri-lab Supplies Inc., USA). According to the method developed by the authors of this paper, the altered soil samples were dried in an oven (40GC, Quincy Lab Corporation, USA) at 50ºC for 12 hours, then passed through a sieve of 2mm size to separate the soil from rocks, roots and detritus. A sample of 25g was weighed and poured in a plastic cylindrical container of 4cm height and 2.5cm of diameter. Each container was placed in the paramagnetic soil meter device for one minute until obtain a steady value that shows the relative alignment of the electrons of the soil sample with the magnetic field generated.

Geospatial mapping and experimental design. The experimental design consisted in a randomized soil sampling across the two lots (Santa Rosa and Chondular), where the sampling points were evenly spread in the field (one point per two hectares) and georeferenced. The information obtained from the soil survey was processed with the software ArcMap of ArcGis 10.1 (ESRI, USA) in order to generate maps of each of the soil characteristics, including MS. The data set of the soil properties was analyzed by descriptive statistics, including a test of normality and Levene’s test of homogeneity of variances. The geostatistical analysis was carried out using the ArcGis tool “Band collection statistics”, a multivariate analysis from the covariance of each pair or layers, from which a correlation matrix is calculated. The correlation matrix shows the values of the correlation coefficients that represent the relationship between the two datasets, and is a measure of dependence between the layers. It is calculated from the covariance between the two layers divided by the product of their standard deviations. It is a number without units. In addition, a Pearson correlation of the data with corresponding significance t-test was done to complement the spatial correlation, using Minitab 17 program.

RESULTS AND DISCUSSION

The descriptive analysis (Table 2) showed high variability of MS in both experimental locations, with values of 44.9% and 42.8%, for Santa Rosa and Chondular lots, respectively. Properties such as total pore space, field capacity and bulk density showed less variability within the study areas.

Tabla 2. Descriptive statistics for the physical properties of soil.

Despite the high variability of the data of MS, when comparing these values with sand and clay content in Chondular lot, a consistent trend was observed, i.e., the areas of high MS match with the areas of high content of sand and the areas of low content of clay (Figure 3). Further relationships were also detected between MS and other properties such as field capacity, bulk density and total pore space, although its correspondence with the former variable was not particularly evident.

Figure 3. Spatial distribution of soil properties in Chondular lot.

Figure 4 shows the values of MS, clay and sand for Santa Rosa lot. In these maps it can be observed the spatial distribution of MS (A), sand (B) and clay content (C), in which the same apparent relation with sand and clay was detected, although unlike Chondular lot, the matching of MS with the areas of high content of sand and low content of clay was only partial.

Figure 4. Spatial distribution of soil properties in Santa Rosa lot. A: MS. B: sand. C: clay.

The relationship between the figures displayed in Figures 3 and 4 analyzed by correlation matrix in ArcGis and Pearson correlation are shown in Table 3. It is corroborated the high and significant correlation between sand and clay content with MS in both locations, being more consistent that of Chondular lot. However, MS did not exhibit a coherent trend with silt and total pore space, in view that it showed a positive and negative correlation in Chondular and Santa Rosa lots, respectively. Field capacity, thermal conductivity and bulk density had higher and significant (P<0.01) spatial correlations with MS in Chondular, whereas these correlations were less appreciable in Santa Rosa location. For drainable porosity, there was little correlation and significance with MS in both sites.

Table 3. Coefficients of spatial and Pearson correlation between MS and the physical properties of the soil evaluated.

Further analysis of the correlation between MS and sand/clay is depicted as dispersion diagrams in Figure 5, of each sampling point with its corresponding value of MS contrasted with that of sand and clay. The determination coefficient revealed the same pattern of direct and inverse relation of the latter properties with MS at both sites, although this relation was not as accurate in Santa Rosa as it was in Chondular lot.

Figure 5. Dispersion diagram of sand and clay content versus MS in Chondular (left) and Santa Rosa (right).

The expression levels of MS can be related to the texture composition (Williams and Cooper, 1990). Matias et al. (2014) found correlations between different physicochemical properties of oxisols and MS, but unlike this work, their results showed that the relation with sand and clay are negative and positive, respectively, and concluded that the spatial variability of MS can be used to infer on clay content and soil base saturation variability. Aversely, but supporting the results of the present investigation, Siqueira et al. (2010) found a positive relation of MS with clay content. In relation with this, it is important to consider that the nature of the clay minerals and MS is dependent on the variation of the soil parent material (Camargo et al., 2014). The presence of ferromagnetic iron in the molecular structure of clays seems to be determinant for the expression of MS. The results reported by Ramos et al. (2017), emphasized that magnetism is most evident in clayey soils. Among the ferromagnetic minerals, magnetite is more easily found in coarser soil fractions, such as sand and silt, whereas maghemite is mostly found in finer fractions, such as clay. Grimley et al. (2004) found MS values consistently higher in well-drained soils and lower in hydric soils, reflecting anaerobic deterioration of both detrital magnetite and soil-formed ferrimagnetics. The latter authors indicated that the higher magnetite content and slower dissolution rate in sandy soils might explain the difference. This in part corroborates the findings of this study, in which the highest MS values were observed in the areas of sandy texture. Hence, it is possible that the water retention capacity of the soil related with its texture and the degradation of magnetite were the mechanism that elicits the paramagnetic response. Other researchers mention that the magnetic behavior is in fact more evident in soils whose clay fraction is greater because in clay, magnetite is oxidized to maghemite, while in the sand fraction, magnetite is oxidized directly to hematite that has lower MS than magnetite (Fontes et al., 2000). Highly variable critical MS values in areas of South Africa, probably a result of differing magnetite contents in soil parent materials among the sites, show the need for calibration of critical MS values on a site specific basis (Grimley et al., 2004). However, these authors also informed that MS values of clayey soils are significantly lower than those of sandy loam soils in the same landscape position, which is consistent with the values here reported. For these reasons, the identification of high spatial correlation between MS and the contents of clay and sand could be an alternative for inferring other texture-related parameters, such as hydraulic conductance.

The relationship between MS and the other physical properties (BD, TPS, FC and DP) can be inferred from the content of clay and sand. Regarding clay, a high content lead to a greater formation of micropores, responsible for water retention in the soil; since microporosity predominates in clayey soils, it occupies most of the total pore space, whereby there is a direct relationship between clay and microporosity. In addition, soil micropores allow greater heat conductance (Gutiérrez et al., 2016), which explains its correlation with thermal conductivity. Moreover, a high content of sand lead to a higher bulk density in the soil, because it is a heavier and thicker material which in turn contributes to the formation of macropores, improving drainable porosity.

It is also important to determine the extent of tillage influence on MS in the selected sampling points. The correlation values (and the coefficient of variation of the properties that rely on the structure, like bulk density and pore space) observed in Santa Rosa were lower and less significant than in Chondular, which is attributed to a recent tillage in the area prior to the sampling. This is an issue previously mentioned by Mathé et al. (2003), who reported magnetic anomalies in marshland soils recently tilled. Uniform mineralogy is typical for both surface soils and samples from different depths in non-disturbed soil profile (Jordanova et al., 2011). This supports the hypothesis that, within the studied area, the observed magnetic differences comes from variations in soil movement and redistribution of the ferromagnetic carriers due to tillage operations. This statement agrees with previous results (Nazarok et al., 2014), who indicated that soil MS has a high degree of statistical relationship with erosion index and humus content, both parameters highly dependent on human intervention.

As pointed out by Grimley et al. (2004), more research is required to validate the MS method according to the different regions, land uses, soil types and parent materials. It is important to determine and compare its critical values in contrasting sites to develop an accurate database.

CONCLUSIONS

The significant correlation values (either by spatial matrix or Pearson) between MS and soil particle composition (particularly clay and sand percentage), are evidence of this property as a pedotransfer function for texture determination. The lack of correlation with silt content could be avoided by the proper determination of sand and clay through MS. However, this parameter might not be suitable as a pedotransfer function for other physical properties such as total pore space, drainable porosity and bulk density. Moreover, the physical disturbance of the soil by tillage, that alters the attributes dependent on the soil structure (like bulk density and pore space), also affects the expression of MS and therefore its correlation with these variables. Notwithstanding the above stated, MS could be a proper tool to the indirect quantification of soil texture with measurable paramagnetism. It must be considered the degree of physical disturbance of the soil that impairs the accuracy of MS as a function of other physical properties.

BIBLIOGRAPHIC REFERENCES

1. Barbieri, D.M.; Marques-Junior, J.; Pereira, G.T. 2008. Variabilidad espacial de atributos químicos de um argissolo para aplicação de insumos à taxa variável em diferentes formas de relevo. Eng. Agr. 28(4):645 - 653. doi: http://dx.doi.org/10.1590/S0100-69162008000400004.

2. Bautista F.; Cejudo, R.; Sánchez, A.; Aguilar, B.; Delgado, M.; Goguitchaichvili, A.; Marín, P.; Gil, J.; Díaz, E. 2013. Propiedades magnéticas y pedogénesis en un perfil de suelo con horizontes contrastantes. Latinmag Lett. 3:1 - 6.

3. Blundell, A.; Dearing, J.A.; Boyle, J.F.; Hannam, J.A. 2009. Controlling factors for the spatial variability of soil magnetic susceptibility across England and Wales. Earth Sci. Rev. 95(3-4):158 - 188. doi: http://dx.doi.org/10.1016/j.earscirev.2009.05.001.

4. Camargo, L.A.; Marques, J.; Pereira, G.T.; De Souza Bahia, A.S. 2014. Clay mineralogy and magnetic susceptibility of Oxisols in geomorphic surfaces. Sci. Agric. 71(3):244 – 256. doi: http://dx.doi.org/10.1590/S0103-90162014000300010.

5. Cortez, L.A.; Marques, J.R.; Peluco, R.G.; Teixeira, D.B.; Siqueira, D.S. 2011. Suscetibilidade magnética para identificação de áreas de manejo específico em citricultura. Energia Agr. 26(3):60 - 79.

6. Dos Reis Barrios, M.; Marques, J.; Rocha, S.; Panosso, A.R.; Silva, D.; Scala, N. 2017. Magnetic susceptibility as indicator of soil quality in sugarcane fields. Revista Caatinga 30(2):287-295. doi: http://dx.doi.org/10.1590/1983-21252017v30n203rc.

7. Fontes, M.P.F.; Oliveira, T.S.; Costa, L.M.; Campos, A.A.G. 2000. Magnetic separation and evaluation of magnetization of Brazilian soils from different parent materials. Geoderma. 96(1-2):81 - 99. doi: http://dx.doi.org/10.1016/S0016-7061(00)00005-7.

8. Grimley, D.A.; Arruda, N.K.; Bramstedt, M.W. 2004. Using magnetic susceptibility to facilitate more rapid, reproducible and precise delineation of hydric soils in the midwestern USA. Catena 58(2):183 - 213. doi: http://dx.doi.org/10.1016/j.catena.2004.03.001.

9. Gutiérrez C., M.A.; Zúñiga E., O.; Ospina-Salazar, D.I. 2016. Effect of three biowastes on the productivity potential of a sodic soil. Agron. colomb. 34(2):250 - 259. doi: http://dx.doi.org/10.15446/agron.colomb.v34n2.55044.

10. Huang, C.C.; Jia, Y.; Pang, J.; Zha, X.; Su, H. 2006. Holocene colluviation and its implications for tracing human-induced soil erosion and redeposition on the piedmont loess lands of the Qinling Mountains, northern China. Geoderma 136(3-4):838 - 851. doi: http://dx.doi.org/10.1016/j.geoderma.2006.06.006.

11. IGAC - Instituto Geográfico Agustín Codazzi; CVC - CORPORACIÓN Autónoma Regional Del Valle Del Cauca. 2004. Levantamiento de suelos y zonificación de tierras del Departamento del Valle del Cauca. IGAC. Bogotá D.C. 775p.

12. Jordanova, D.; Jordanova, N.; Petrov, P. 2014. Pattern of cumulative soil erosion and redistribution pinpointed through magnetic signature of Chernozem soils. Catena. 120:46 - 56. doi: http://dx.doi.org/10.1016/j.catena.2014.03.020.

13. Jordanova, D.; Jordanova, N.; Atanasova, A.; Tsacheva, T.; Petrov, P. 2011. Soil tillage erosion estimated by using magnetism of soils – a case study from Bulgaria. Environ Monit. Assess. 183:381 - 394. doi: http://dx.doi.org/10.1007/s10661-011-1927-8.

14. Kanu, M.O.; Meludu, O.C.; Oniku, S.A. 2014. Comparative study of top soil magnetic susceptibility variation based on some human activities. Geofí. Int. 53(4):411 - 423. doi: http://dx.doi.org/10.1016/S0016-7169(14)70075-3.

15. Luque, E.C.L. 2008. Propiedades magnéticas de los óxidos de hierro en suelos mediterráneos. En: https://dialnet.unirioja.es/servlet/tesis?codigo=54579; consulted: agosto, 2016.

16. Marqués, J.R.; Siqueira, D.S.; Camargo, L.A.; Teixeira, D.B.; Barrón, V.; Torrent, J. 2014. Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf. Geoderma. 219-220:63 - 71. doi: http://dx.doi.org/10.1016/j.geoderma.2013.12.007.

17. Mathé, V.; Lévêque, F. 2003. High resolution magnetic survey for soil monitoring: detection of drainage and soil tillage effects. Earth Planet Sc. Lett. 212(1-2):241 - 251. doi: http://dx.doi.org/10.1016/S0012-821X(03)00241-3.

18. Matias, S.S.R.; Marques, J.; Siqueira, D.S.; Pereira, G.T. 2014. Outlining precision boundaries among areas with different variability standards using magnetic susceptibility and geomorphic surfaces. Eng. Agríc. 34(4):695 - 706. doi: http://dx.doi.org/10.1590/S0100-69162014000400009.

19. McBratney, A.B.; Minasny, B.; Cattle, S.R.; Vervoort, R.W. 2002. From pedotranfer functions to soil inference systems. Geoderma. 109(1-2):41 - 73. doi: http://dx.doi.org/10.1016/S0016-7061(02)00139-8.

20. Nazarok, P.; Kruglov, O.; Menshov, O.; Kutsenko, M.; Sukhorada, A. 2014. Mapping soil erosion using magnetic susceptibility. A case study in Ukraine. Solid Earth Discuss. 6:831 - 848. doi: http://dx.doi.org/10.5194/sed-6-831-2014.

21. Patil, N.G.; Singh, S.K. 2016. Pedotransfer functions for estimating soil hydraulic properties: a review. Pedosphere 26(4):417 - 430. doi: http://dx.doi.org/10.1016/S1002-0160(15)60054-6.

22. Pedroso, I. 2013. Zonación de la contaminación por metales pesados en la cuenca del Almendares según mapeo de la susceptibilidad magnética. Miner. Geol. 29(3):1 - 17.

23. Rahimi, M.R.; Ayoubi, S.; Abdi, M.R. 2013. Magnetic susceptibility and Cs-137 inventory variability as influenced by land use change and slope positions in a hilly, semiarid region of west-central Iran. J. App. Geophy. 89:68 - 75. doi: http://dx.doi.org/10.1016/j.jappgeo.2012.11.009.

24. Ramos, P.V.; Dalmolin, R.S.; Marques, J.; Siqueira, D.S.; Almeida, J.A.; Moura-Bueno, J.M. 2017. Magnetic susceptibility of soil to differentiate soil environments in southern Brazil. Rev. Bras. Cienc. Solo. 41:e0160189. doi: http://dx.doi.org/10.1590/18069657rbcs20160189.

25. Royall, D. 2001. Use of mineral magnetic measurements to investigate soil erosion and sediment delivery in a small agricultural catchment in limestone terrain. Catena. 46(1):15 - 34. doi: http://dx.doi.org/10.1016/S0341-8162(01)00155-2.

26. Siqueira, D.S.; Marques, J.; Matias, S.R.; Barrón, V.; Torrent, J.; Baffa, O.; Oliveira, L.C. 2010. Correlation of properties of Brazilian Haplustalfs with magnetic susceptibility measurements. Soil Use and Management. 26(4):425 - 431. doi: http://dx.doi.org/10.1111/j.1475-2743.2010.00294.x.

27. Souza, I.G.; Costa, A.C.S.; Vilar, C.C.; Hoespers, A. 2010. Mineralogia e susceptibilidade magnética dos óxidos de ferro do horizonte B de solos do Estado do Paraná. Ciência Rural. 40(3):513 - 519.

28. Torrent, J.; Liu, Q.S.; Barrón, V. 2010. Magnetic minerals in Calcic Luvisols (Chromic) developed in a warm Mediterranean region of Spain: origin and paleoenvironmental significance. Geoderma. 154(3-4):465 - 472. doi: http://dx.doi.org/10.1016/j.geoderma.2008.06.020.

29. Torrent, J., Liu, Q.S., Bloemendal, J., Barro, N.V. 2007. Magnetic enhancement and iron oxides in the upper Luochuan loess – paleosol sequence, Chinese Loess Plateau. Soil Science Society of America Journal.71:1 - 9. doi: http://dx.doi.org/10.2136/sssaj2006.0328.

30. Urcia, O.; Larrasoaña, J.C.; Muñoz, A.; González, A.; Pérez, A.; Luzón, A.; Román, T.; Villalaín, J. 2012. La susceptibilidad magnética como marcador paleoambiental en un abanico aluvial del Pleistoceno superior: la cuenca de Añavieja, Cordillera Ibérica (NE de España). p. 742 - 745. En: VIII Congreso Geológico de España. Oviedo, España.

31. Williams, R.D.; Cooper, J.R. 1990. Locating soil boundaries using magnetic susceptibility. Soil Sc. 150:889 - 895.

32. Zúñiga, O.; Benavides, J.; Ospina-Salazar, D.I.; Jiménez, C.O.; Gutiérrez, M.A. 2016. Magnetic treatment of irrigation water and seeds in agriculture. R. Ing. Compet. 18(2):217 - 232.

33. Zúñiga Escobar, O.; Reyes Trujillo, A. Electro-thermal measurement device for e.g. evaluating compaction of agricultural soil has upper lid having orifice through which soil sample is introduced for evaluating soil's thermal conductivity. En: https://worldwide.espacenet.com/publicationDetails/biblio?FT=D&date=20061001&DB=&locale=en_EP&CC=ES&NR=2259498A1&KC=A1&ND=4; consulted: agosto, 2016.