INTRODUCTION
The predominant soils in the central-southern zone of Chile are of volcanic origin. The forestry, agricultural, and livestock activities carried out in this area, along with the natural processes occurring in the soil, modify the chemical, physical, and biological properties of the soil. The Andisol and Ultisol orders are the most extensively used for agricultural purposes in this area (Pizarro, 2000). Andisols cover approximately 124 million ha worldwide and account for approximately 0.8% of the total land area (Tsai et al., 2010).
Soils derived from volcanic materials represent between 50% and 60% of the total arable land in Chile. These soils cover large areas under agricultural (cereal production) and livestock use. They also support a considerable part of the forests (Besoain, 1985).
Volcanic ash soils have a number of unique properties that are directly attributable to the properties of the parent material (Shoji et al., 1993). They are characterized by a mineralogy dominated by non-crystalline components, such as clay minerals, organic matter (OM), Fe, and aluminum oxides. This provides these soils with distinctive characteristics, such as pH-dependent charges and electrolyte concentration present in the solution (Moustakas and Georgoulias, 2005). These distinctive properties largely favor the formation of amorphous materials and organic C (OC) accumulation, which are the main pedogenic processes in soils derived from volcanic ash. The clay fraction of these soils is dominated by the presence of allophane, para-allophanic components, imogolite, and several phyllosilicates (smectites, chlorites) in different amounts; these are products of the weathering of pyroclastic materials from recent volcanic deposits (Besoain, 1985).
Iron plays an important role within the mineralogy of these soils and can explain part of the pedogenic processes. It is used as an indicator of the environmental conditions in which soil is formed and developed (Sánchez and Rubiano, 2015). Varying amounts of Fe are found in the soil, but levels are usually high. It is the most abundant microelement in the soil, but only a small fraction corresponds to assimilable forms in most cases (Solís, 1998). Therefore, Fe is scarcely available because it can easily change its oxidation state and form insoluble hydroxides (Harrington and Crumbliss, 2009). It can form Fe-humus complexes, oxides (crystalline and non-crystalline), and ferromagnesium silicates (Aguado-Santacruz et al., 2012). Magnetite, goethite, hematite, and ferrihydrite are highly crystalline Fe minerals, whereas minerals with little or no crystalline structure are called amorphous oxyhydroxides (Torrent et al., 1980). The weathering of these primary minerals releases Fe to the soil solution where it can be used by organisms, form part of different organic complexes, or form secondary minerals, such as sulfides, carbonates, clay minerals, but mainly oxides and hydroxides. These minerals have different compositions and degree of crystallization, which control Fe solubility in the soil (Murad and Fischer, 1988). The Fe oxides, which include hydroxides, oxides, and oxyhydroxides, promote aggregation, adsorption of nutrients and contaminants, and also serve as electron receptors (Stoppe et al., 2015). Therefore, it is important to understand the composition of the soil and sedimentary Fe pools.
The dynamics and availability of Fe are governed by the amount and type of clay present, pH, texture, and aeration of the soil. Organic matter is an important factor in the Fe availability because it forms stable compounds that generate the so-called chelates (soluble organic complexes of Fe) (Solís, 1998). Organic matter favors growth of microbial populations, which can consume oxygen under waterlogging conditions or decompose OM to release Fe from organic compounds in forms that can be assimilated by plants (Aguado-Santacruz et al., 2012).
The application of dairy slurry to the soil is a common fertilization practice; it can meet the total or partial fertilization requirements of permanent grasslands in southern Chile (Salazar et al., 2003) and provide microelements. Its use in soils with low OM content represents a valuable resource to increase this parameter (Salazar et al., 2007).
The Fe contents can be determined by physical and chemical methods. Physical methods include X-ray diffraction, differential thermal analysis, Mössbauer spectroscopy, and infrared spectroscopy (Hernández and Meurer, 1997). Chemical extraction approaches include the dithionite-citrate-bicarbonate (DCB) (Mehra and Jackson, 1960), ammonium oxalate (Feox), and sodium or potassium pyrophosphate methods. The DCB method (Mehra and Jackson, 1960) allows extracting available Fe (Fed), which represents most of the pedogenic Fe minerals (crystalline and amorphous), and determining labile Fe oxides (Vazques et al., 2014). The acid ammonium oxalate extraction method (Schwertmann, 1964) is used to estimate the amount of amorphous or poorly crystalline Fe in soils and sediments. The sodium or potassium pyrophosphate method is commonly used to extract Fe from soil organic complexes (Smith, 1994). Some Fe ratios are used to determine the presence of these minerals in the soil. For example, Feox:d is a measurement of the amount of total pedogenic Fe, while the difference, Fed-ox, is used to estimate crystallized Fe (goethite and hematite) (Acevedo- Sandoval et al., 2002).
The objective of this study was to determine the effect of applying dairy slurry over a 12-yr period on the levels of available, amorphous, and total Fe in an Andisol soil as indicators of pedogenic alteration.
MATERIALS AND METHODS
Study site
The study was conducted at the Livestock Experimental Station of the Universidad de Concepción (36º33’17’’ S, 71º52’45’’ W), Chillán, Ñuble Region, Chile. It is located in the central depression of the country and has a temperate Mediterranean climate with a mean annual temperature of 13.9 °C and a minimum of 3 °C during July and a maximum of 28.6 °C in January. Mean annual rainfall is 1025 mm.
The soil is a Humic Haploxerands (Arrayán Series) (Stolpe, 2006). The experimental site has a predominantly silt loam texture. The soil is well-structured on the surface with high porosity associated with good rooting in the whole pedon (Table 1). It has a flat or almost flat topography with moderate permeability and slow surface runoff.
Experimental design
A completely randomized design with repeated measurements was used, consisting of six treatments and four replicates per treatment. A control treatment (untreated) was also included. The experimental unit was a 1 ha (10 000 m2) plot. The treatments were T0 (control with no slurry application), T2 (2-yr slurry application), T4 (4-yr slurry application), T6 (6-yr slurry application), T8 (8-yr slurry application), T10 (10-yr slurry application), and T12 (12-yr slurry application).
Slurry was accumulated in an open air lagoon (10 m long, 2 m wide, and 1.5 m deep) with a 15 000 L capacity and applied at a mean rate of 15 m3 ha-1 yr-1 using a cart with a 5000 L capacity. The slurry used in this study complies with the provisions of “Emission standard for the regulation of pollutants associated with liquid waste discharges to marine and inland surface waters in Chile”, which ensures the composition of waste over time, originating from a homogeneous and constant flock, number of animals, breed, diet, and winter housing period (Table 2). On average, slurry was applied every 15 d. The experiment was conducted on permanent Lolium perenne L. and Trifolium repens L. grasslands. For maintenance fertilization, fertilizers were applied at rates of 46 kg N ha-1, 120 kg P ha-1, and 60 kg K ha-1.
Sampling
A test pit was dug on each hectare plot in each treatment to identify genetic horizons and determine the depth and morphological characteristics for pedogenic horizon differentiation (A and B) (Table 1). Twenty-five subsamples were randomly taken from the soil with a cylindrical soil sampler, thus forming a homogeneous sample for each hectare in each treatment and pedogenic horizon. Therefore, each treatment consisted of four combined samples, which were placed in polyethylene containers and sent for analysis. Soil samples were dried at room temperature, sieved through a 2 mm sieve according to the procedure described by Sandoval et al. (2012), and then analyzed.
Soil chemical analysis
Extraction of different forms of Fe
The DCB method proposed by Mehra and Jackson (1960) was used to selectively extract available Fe oxides. A 2 g homogenized and previously sieved soil sample (< 2 mm) was weighed on an analytical balance (model L5201, BEL Engineering, Monza, Italy). A 50 mL solution of (NH4)2C6H6O7 (0.26 M) - NaHCO3 (0.11 M) and 1 g of sieved soil were placed in an Erlenmeyer flask and agitated for 16 h in a shaking incubator (BJPX-Wichita, Biobase, Jinan, China) at 120 rpm and 25 °C . The mixture was filtered through filter paper (Whatman Nr 42, 125 mm) and stored in a polyethylene jar until the samples were read in an atomic absorption spectrophotometer (AAS). To extract amorphous Fe oxides, 2 g homogenized and previously sieved soil (< 2 mm) were weighed on an analytical balance (BEL Engineering) and mixed with 50 mL extractant (NH4)2C2O4 × H2O (0.2 M) at pH 3.0 in a 100 mL Erlenmeyer flask according to the method described by Schwertmann (1964). The solution was agitated in a shaking incubator (Biobase) in the dark at 120 rpm and 25 °C for 120 min. Once the agitation period was complete, the mixture was filtered through filter paper (Whatman Nr 125, 125 mm) and the extract was placed in a polyethylene container until the samples were read in an AAS. Calcination and acid digestion procedures were used to extract total Fe (Sadzawka et al., 2006).
Atomic absorption spectrophotometry (AAS)
Based on Fe levels measured in AAS, 1:100 dilutions were made for all the evaluated forms of Fe. The lamp to analyze Fe and the light wavelength were selected in the AAS (GE- 503115, SOLAAR 969, Thermo Fisher Scientific, Waltham, Massachusetts, USA). The Fe patterns were determined, the calibration curve was performed, and the corresponding samples were read for different forms of Fe. Results were expressed in g kg-1 of sample.
RESULTS AND DISCUSSION
Selective extractions of available, amorphous and total Fe
The Fed contents, which are labile Fe, varied between 4.55 and 9.85 g kg-1 in the A horizon and accounted for 14.63% to 29.09% of total Fe (Fet) content (Table 3). These levels were lower than those observed by Lu et al. (2014) in an Ultisol with values up to 35.38 g kg-1. However, a recent study conducted by Kilic et al. (2018) reported similar values to those obtained in our work in soils of volcanic origin (Andisol). Significant differences (p ≤ 0.05) were found between all the treatments studied in the A horizon; Fed levels increased in treatments T2, T4, and T6, but decreased significantly in treatments T8, T10, and T12. A similar behavior was observed in the B horizon (Table 3). However, values in the B horizon were lower compared with the surface A horizon in all treatments. Acevedo-Sandoval et al. (2011) described a similar behavior in the distribution of Fe oxides in three soil profiles of volcanic origin in the Northwest Region of the State of Mexico.
Amorphous Fe (Feox) are usually non-crystalline and poorly arranged compounds (Hernández and Meurer, 1997). The Feox values were higher in the A horizon. These results concur with findings by Pérez et al. (2018), who conducted selective Fe extractions in an Andisol from eastern Antioquia, Colombia, and reported that the Feox content decreases as depth in the profile increases. The Feox content for the surface horizon varied between 5.23 and 7.79 g kg-1 (Table 3). Similar values were reported by Chevallier et al. (2010) in an Andisol in Martinique, French West Indies. Nonsignificant differences were found (p ≤ 0.05) in Feox levels between treatments T0, T2, and T12. For the B horizon, Feox values varied between 3.03 and 7.34 g kg-1 and accounted for 9.78% to 23.72% of the Fet content in the soil. The Feox values increased in T2, T4, and T6 compared with the control. However, levels decreased significantly in T8 and T10, but increased again in T12. Acevedo et al. (2014) indicated that Feox is usually present in soils with pH-dependent charges and a high P fixation capacity; this coincides with the characteristics of the parental material of the soil in the present study. Acevedo-Sandoval et al. (2002) reported that Feox predominates in the upper horizons of Inceptisols and Alfisols. In our study, this parameter showed a higher coefficient of variation compared with Fed and Fet in both horizons. Van Dam et al. (2008) reported that Feox is considered to be unstable under certain environmental conditions that occur during initial soil formation or in the presence of silicates, which prevent the formation of crystalline Fe oxides.
Means in columns followed by the same letter are not different according to Tukey’s test (P ≤ 0.05).
Fed: Available Fe extracted by dithionite-citrate-bicarbonate (DCB) method; Feox: amorphous Fe extracted with ammonium oxalate; Fet: total Fe extracted by acid digestion; CV: coefficient of variation; LSD: least significant difference.
The Fet contents ranged from 37.80 to 31.10 g kg-1 and 38.35 to 29.70 g kg-1 for the A and B horizons, respectively (Table 3). Winkler et al. (2016) evaluated the Fe response in rice growing soils (Vertisols, Andisols, and Alfisols) and reported Fet values ranging from 59.3 to 62.4 g kg-1. These values were higher than those obtained in our study. A similar behavior of Fet levels was observed in both genetic horizons, and the values decreased as the number of years of slurry application increased.
Ratios of different evaluated forms of Fe
The content of crystalline Fe oxides obtained from the difference between Fed and Feox (Parfitt and Childs, 1988) was low, and values were < 3.4 g kg-1 in both studied horizons (Table 4). The low oxidization reduction potential in the soil could explain this situation (Smith, 1994). The Feox:d ratio (Table 3) has been widely used to determine the degree of crystallinity of pedogenic Fe forms and soil evolution (Kendrick and McFadden, 1996). In general, the progressive weathering of soils leads to a decreased Feox:d ratio as soil aging advances (Vacca et al., 2003). However, these values are typically high in Andisols (> 0.75), while more developed soils exhibit values < 0 (Kleber et al., 2004). In our study, the Feox:d ratio ranged from 0.66 to 2.12, suggesting that these soils are young and poorly evolved. The Feox:d values for treatments T2, T4, T6, T8, and T10 were significantly lower (p ≤ 0.05) than those in T0 and T12 (Table 4). This is explained by a greater pedological development in both horizons (A and B) as a result of slurry application. However, T2 did not show greater development in the A or B horizons, which might be due to the degradation of the grassland after 12 yr and consequently to decreased bioprocesses. Most of the Fe in the profile is part of the silicate structure (Fet -d) and accounts for 70.91% to 88.89% of the total Fe in both genetic horizons.
Chemical and microbiological properties (pH, electrical conductivity, and organic matter)
Soil acidity (pH) was significantly higher (p ≤ 0.05) in T0 in the A horizon compared with the different evaluated treatments (Table 5). As the time (years) of slurry application increased, pH became less acidic. Whalen et al. (2000) stated that there is evidence that animal slurries can increase pH in acid soils due to the buffer effect caused by the content of bicarbonates and organic acids. The pH values for the B horizon ranged from 6.41 to 6.81, varying from slightly acidic to almost neutral. It is important to note that there was a slight decrease in acidity in both profiles as depth increased.
Means in columns followed by the same letter are not different according to Tukey’s test (P ≤ 0.05). Fed-ox: Crystalline Fe oxide; Feox:d: activity index; Fet -d: silicate Fe; LSD: least significant difference.
Means in columns followed by the same letter are not different according to Tukey’s test (P ≤ 0.05).
EC: Electrical conductivity; OM: organic matter; MR: microbiological respiration; CV: coefficient of variation; LSD: least significant difference.
Nonsignificant differences were found (p ≤ 0.05) between treatments T4, T6, T8, and T10 (Table 5) for electrical conductivity (EC) in the A horizon. Values ranged from 0.13 to 0.32 dS m-1 with T0 and T2 recording the highest values. The EC is significantly influenced by water content, clay content, and the presence of interchangeable ions capable of conducting an electrical current, which also affects the nutritional characteristics of the soil. Wei et al. (2009) indicated that high variations in EC can be explained by unequal applications of amendments or fertilizers. In the B horizon, EC was less variable and values varied between 0.09 and 0.13 dS m-1.
Organic matter contents were higher in the A horizon treatments compared with those in the B horizon (Table 5). Organic C (OC) from grasslands and slurry applications would explain this situation. A similar behavior was observed by Venegas (2008) in an Andisol of the Osorno series in Chile. He indicated that OM content decreases as depth in the profile increases. The OM levels in the A horizon showed significant differences (p ≤ 0.05) in all the treatments; T6 and T8 recorded OM levels > 15%. These results are similar to those obtained by Alarcón et al. (2010), who reported similar OM values in a soil of volcanic origin in an artificial pasture. For the B horizon, OM contents varied between 7.10% and 10.35%; they were in the lower OM percentage and similar (7.2%) to those obtained by Fernández et al. (2017) in an Alfisol in the town of San Ignacio, Chile (36°46’51.01” S, 72°00’29.62” W). A similar trend was observed in both horizons. The T8 treatment recorded the highest OM levels, while T10 and T12 had decreased OM. This effect could be attributed to grassland degradation in T10 and T12. OM can be subjected to leachate losses and the humification processes diminished; therefore, the good use of these organic residues would be based on the quality of the grassland, which in turn would be related to microorganism activity. The lower microbial activity decreases OM mineralization, restricting N availability (NH4+ and NO3-) and negatively impacting crop growth (Robertson and Groffman, 2015).
Tables 6 and 7 show correlations between the studied variables. There is a correlation between the different forms of Fe and the physicochemical properties of the soil under study. Specifically, there is a positive and highly significant correlation (p ≤ 0.001) between OM and the Fed fraction in the A horizon. Conversely, there is a highly significant negative correlation (p ≤ 0.001) between OM and Fet. There is a negative correlation between pH and all the forms of evaluated Fe. This coincides with the results reported by Bronick and Lal (2005). In our study, the negative correlation between pH and Fet was highly significant (p ≤ 0.001) in both genetic horizons.
For the selected physicochemical properties, no correlation was found between EC and the forms of Fe evaluated in both horizons.
EC: Electrical conductivity; OM: organic matter; Fed: available Fe extracted by dithionite-citrate-bicarbonate (DCB) method; Feox: amorphous Fe extracted with ammonium oxalate; Fet: total Fe extracted by acid digestion.
*, **Significant at the 0.05 and 0.01 probability levels, respectively.
CONCLUSIONS
Slurry applications in the soil resulted in a significant increase in the available Fe (Fed) contents up to 8 yr in both the A and B horizons. However, total Fe (Fet) levels decreased as the number of years of application increased.
The amorphous Fe (Feox): Fed ratio (Feox:d) indicates that these soils are pedologically underdeveloped. It also confirms that slurry application has an anthropogenic effect on its pedogenesis.
There was a positive correlation between organic matter (OM) and Fed. This would indicate that the pedogenesis of this soil can be accelerated as a result of increased OM, but conditioned to the bioprocesses related to the state and quality of the grassland.