Soil aggregation status under conventional and conservation management
Several studies on agricultural production in different countries (Brazil, China, India, Canada, Argentina, and the United States of America) have shown that conservation management, in particular, NT practices, may positively influence soil quality. Factors such as time since system adoption, continuous input of organic matter (plant residues), beneficial action of plant root residues, and protection provided by cover crops are associated with the increase in soil quality (Andrade et al. 2018).
In Paraná State, Brazil, Loss et al. (2014) observed higher MWD values in Rhodic Paleudalfs (Red Nitosol, Nitossolo Vermelho) under NT and pasture than in soil under CT. The authors attributed the results to the high organic carbon contents stemming from the well-developed root systems of Axonopus compressus in pasture and ryegrass in NT. Sales et al. (2016), in evaluating the structural quality of a Hapludox (Red-Yellow Latosol, Latossolo Vermelho-Amarelo) in Minas Gerais, Brazil, also observed higher MWD values and macroaggregate proportions at NT areas than at CT areas, both of which were planted with maize, sorghum, and sunflower. The authors argued that NT contributes to soil aggregation, especially in fields planted with grasses.
In a Typic Hapludoll in Northeast China, Zheng et al. (2018) found that long-term adoption of conservation management significantly increased macroaggregate proportion and MWD values compared with conventional management. According to the authors, aggregate stability is associated with the ability of soil to resist exogenous action and remain stable when exposed to changes in the external environment. In a Ferralsol (Yellow Latosol, Latossolo Amarelo) in Piauí, Brazil, Silva et al. (2018) found that soil aggregation was highest in an old NT system, followed by pasture and forest (eucalyptus). However, at a CT site, aggregate stability was low.
The results of macroaggregate proportion and MWD in the 0.00–0.10 m layer of 6NT and 18NT might be related to the lower soil mobilization associated with organic carbon content and fractions. The results also suggest that the use of cover crops (millet, Brachiaria, and Crotalaria) in rotation with maize/soybean or common bean (annual crops) favored the formation of larger aggregates; macroaggregate proportion and MWD values were higher in NT systems than in CT (Table 2). The higher proportions of meso- (A2, A3, and A4) and microaggregates (A5 and A6) in CT further confirmed that conventional practices decreased the proportion of macroaggregates (A1) compared with NT (Table 2).
In Central India, Somasundaram et al. (2018) found higher proportions of meso- and microaggregates and lower MWD values in CT than in NT and reduced tillage systems. Similar results were reported by Comin et al. (2018) in a Humic Dystrudept (Humic Cambisol, Cambissolo Húmico) in Santa Catarina State, Brazil. According to the authors, the negative effects of management on soil quality were associated with reductions in aggregate stability, especially with the decrease in macroaggregate proportion.
The higher proportions of meso- and microaggregates in the 0.00–0.10 m layer of CT may be attributed to intense tillage practices (plowing and leveling) and lack of soil cover and plant diversity. Such factors contribute to reducing macroaggregate formation and stabilization and increasing the proportion of meso- and microaggregates in surface/subsurface layers. Thus, it can be inferred that macro-, meso-, and microaggregate proportions were directly related to the low MWD values in CT (Table 2). Another negative effect of conventional management on soil aggregation status is the decrease in hierarchization process. Higher proportions of soil meso- and microaggregates in detriment to macroaggregates indicate that hierarchization is at the early-to-mid stage (Somasundaram et al. 2018). Thus, the high macroaggregate proportions and MWD values in NT indicate an advanced stage of hierarchization process (Table 2).
Effect of management systems on SOM compartments
The effects of management practices on soil carbon dynamics may vary depending on the quantity and quality of organic waste deposited on the soil surface. During soil preparation, SOM may be redistributed in the environment. Given this, small variations in organic carbon content and fractions can significantly affect aggregate stability (Somasundaram et al. 2018). SOM exerts a decisive influence on the formation and stabilization of aggregates because of its diverse molecular structure and the differences in decomposition rates of SOM fractions in soil (Schiller et al. 2018).
Regarding total SOM contents, Loss et al. (2014) observed higher carbon contents in the 0.00–0.05 m layer of NT soil and lower carbon contents in the 0.00–0.10 m layer of CT soil. Sales et al. (2016) found that NT accumulated more soil carbon than CT in the 0.00–0.05 m layer in all cropping systems. Silva et al. (2018) found that carbon contents were lower in an 8-year-old CT system and a 6-year-old forest than in the other systems evaluated. Somasundaram et al. (2018) observed that management practices had a significant effect on the contents of carbon associated with soil aggregates, especially macroaggregates. The authors argued that higher carbon values favored the formation of macro- and mesoaggregates in conservation systems compared with conventional ones.
In Santa Catarina, Brazil, Comin et al. (2018) found that the use of grasses in succession or rotation with other crops under NT increased carbon contents in the surface layer, particularly when in rotation with maize and winter grasses. In Northeast China, Zheng et al. (2018) observed that the contents of carbon associated with soil aggregates in conservation management systems were higher than in conventional systems. The findings were related to the large amount of SOM and good soil structure of the conservation system. In a Hapludox (Red-Yellow Latosol, Latossolo Vermelho-Amarelo) in Mato Grosso, Brazil, 9 years of NT led to higher carbon contents at 0.00–0.05 m depth than CT (Souza et al. 2018). The authors attributed this effect to the increased protection of SOM against rapid microbial decomposition in NT soil. Through the interaction of the organic fraction with soil structure, SOM is protected from oxidation, maximizing the formation of more recalcitrant and stable carbon fractions (Blanco-Canqui and Lal 2004).
Crop rotation, used in NT systems, also contributes to carbon sequestration and maintenance. This practice increases the diversity of substrates and exudates released by plant roots and during above-ground organic matter decomposition. Crop residues and root biomass stimulate microbial communities, especially that of mycorrhizal fungi in soil aggregates, positively influencing carbon and aggregation dynamics (Zhang et al. 2012). Ramalho et al. (2019), in studying an Umbric Ferralsol (Latossolo Bruno) in Paraná State, Brazil, found that a 9-year-old NT system promoted greater carbon accumulation in the 0.00–0.20 m layer than CT. This finding underscored the potential of NT in improving soil quality and increasing SOM. In the same region, Assunção et al. (2019) found that a Ferritic Ferralsol (Red Latosol, Latossolo Vermelho) under NT showed higher carbon contents at 0.00–0.10 m depth than soil under pasture or CT. According to the authors, NT is one of the most efficient agricultural practices for sequestering and accumulating carbon in soil.
From the above, it is possible to infer that the potential of NT in maintaining or increasing carbon contents in tropical soils depends on the quantity and quality of organic residues added to the soil surface as well as the magnitude of anthropogenic disturbances (Souza et al. 2018). The results of Sales et al. (2016), Silva et al. (2018), Somasundaram et al. (2018), Comin et al. (2018), Zheng et al. (2018), Souza et al. (2018), Ramalho et al. (2019), and Assunção et al. (2019) agree with our findings for TC (Table 3).
The lower TC contents in the 0.00–0.05 m layer of CT soil (Table 3) can be attributed to disaggregation caused by soil tillage (Balesdent et al. 2000; Comin et al. 2018). Aggregate fragmentation is intensified with the use of agricultural implements and irrigation/drying cycles (Comin et al. 2018). Soil disaggregation was evidenced by the higher proportion of meso- and microaggregates and the lower MWD values in CT (Table 2).
Previous studies have also investigated more labile SOM fractions. Assunção et al. (2019) observed higher values of POXC in soil (0.00–0.05 m depth) under native forest (1130 mg kg−1) as a function of increasing organic matter input, although values did not differ significantly from those in NT and pasture systems (980 and 900 mg kg−1, respectively). POXC contents are mainly influenced by management practices that contribute to SOM accumulation or stabilization (long-term carbon sequestration) (Hurisson et al. 2016). In quantifying POXC and LOM contents, Santos et al. (2020) did not observe differences in POXC contents but found significant variations in LOM between sites with Inceptisol (Haplic Cambisol, Cambissolo Háplico) under sugarcane in southeastern Brazil, in agreement with the results of the current study. The results indicate that the more labile chemical fraction (POXC) was less sensitive to management practices than the more labile physical fraction (LOM) of SOM (Table 3).
Loss et al. (2014) observed higher LOM contents at NT, pasture, and forest sites and associated LOM results with MWD values. In assessing a 6-year-old NT site and a CT site of a Typic Hapludalf in Canada, Ramnarine et al. (2015) observed that light fraction weight and organic carbon content were higher at the NT area. The authors argued that NT provides better protection of this SOM fraction, at least in the short term. As noted by Briedis et al. (2018), such an increase in the carbon contents of labile and stabilized fractions of SOM at NT areas may result in carbon sequestration in the long term.
Briedis et al. (2018) investigated SOM dynamics at NT, CT, and native vegetation areas in tropical and subtropical regions of Brazil. The authors found higher contents of coarse particulate organic carbon (cPOC) at NT and native vegetation areas. In a Haplic Kastanozem in Argentina, Duval et al. (2018) found that cPOC contents were 58 % higher at NT areas than at CT areas. cPOC and LOM are labile SOM fractions with similar diameters (2.0 > Ø ≥ 0.25 mm); these parameters are efficient indicators of soil quality. According to Duval et al. (2018), cPOC can be used as an initial indicator of changes in soil quality resulting from soil management practices.
Following physical protection via encapsulation within macroaggregates, labile SOM fractions (e.g, LOM, cPOC, or light fraction) become a substrate for soil microfauna. These fractions also act as a nucleation site for the formation of microaggregates within macroaggregates, contributing to hierarchization process (Six et al. 2000; Briedis et al. 2018). Also noteworthy is the intrinsic relationship between litterfall/crop residues and LOM. Incorporation of plant material in soil can be carried out directly, using soil preparation techniques, or indirectly, through biological homogenization. Biological homogenization is more efficient in NT than in CT systems (Buurman and Roscoe 2011). The environmental conditions of NT sites favor biological activity, especially that of soil invertebrates. These individuals are responsible for various ecosystem functions and services, including SOM fragmentation and decomposition, aggregate formation, and stabilization and accumulation of organic carbon (Lavelle et al. 2006; Brown et al. 2018).
The higher LOM contents in 6NT and 18NT soils (Table 3) are likely associated with their higher macroaggregate proportion, MWD, and TC (Tables 2 and 3). In CT soil, the low LOM contents (Table 3) can be attributed to the destructive effect of tillage on soil macroaggregates, as tillage exposes physically protected LOM to microbial attack in a more oxidizing environment (Cambardella and Elliott 1993; Briedis et al. 2018). This hypothesis is supported by the aggregation status of CT soil (Table 2).
Soil quality indices under native vegetation
Soil quality indices were lowest in CA, except for LILOM at 0.05–0.10 m depth, which was lower in 6NT (Table 4). The increase in soil carbon content is a slow process because of the complexity of stable organic fractions. Such fractions depend on the quantity and quality of deposited organic waste and on the prevailing climatic conditions, which directly affect decomposing microorganisms (Torres et al. 2019). These factors help explain the results of carbon indices (CPITC, LIPOXC, and LILOM) in NT systems compared with the reference area.
6NT and 18NT clearly influenced soil quality indices, particularly SIMWD, CPITC, and LILOM (Table 4). In the 0.00–0.05 m layer, Sales et al. (2016) found that the aggregation status of NT was similar to that of native forest (SIMWD ≥ 1.0); CT, by contrast, had lower soil aggregation (SIMWD < 1.0) and carbon accumulation. Rosset et al. (2019) observed an increase in CPITC as a function of time since NT implementation, with values of 0.64, 0.67, and 0.76 at NT areas aged 6, 14, and 22 years, respectively; significant differences were observed between NT areas aged 6 and 22 years. According to the authors, this pattern is related to the slow and gradual increase of TC contents in soil under NT.
Ramalho et al. (2019) reported that NT sites had higher LI values than CT, reinforcing the hypothesis that conservation management can enhance soil quality. High LI values were also recorded by Rosset et al. (2019) at NT areas, especially after 22 years of system implementation. In the 0.00–0.20 m layer, Duval et al. (2018) found that CPI and LI were higher in NT than in CT systems. According to the authors, soil carbon indices were more sensitive to management practices than total organic carbon.
Effect of management practices on CO2-C contents
Changes in SOM quantity and quality in response to management practices directly influence microbial activity and, consequently, CO2-C emission from aggregates during incubation. This fraction is best associated with environmental practices and conditions that promote rapid SOM mineralization (Hurisson et al. 2016; Wade et al. 2018). The greater release of CO2-C by NT aggregates (Fig. 1) can be attributed to practices such as crop rotation (maize/soybean or common bean), use of cover crops (millet, Brachiaria, and Crotalaria), and minimal soil disturbance.
The lower CO2-C values of CT aggregates throughout incubation (Fig. 1) can be related to the negative effects of conventional management practices on soil structure, as discussed above. A similar finding was obtained by Loss et al. (2014). The authors attributed their results to the detrimental effect of CT on soil aggregates, which increase SOM mineralization rates, culminating in low aggregation and low LOM contents. Low availability of the most labile physical fraction (e.g., LOM) for soil biota might have resulted in the low CO2-C in CT samples.
Variations in the evolution of CO2-C during the first 15 days of incubation (Fig. 1) might have been caused by the consumption of SOM by soil microbiota. These microorganisms release CO2-C when decomposing available SOM, resulting in peaks of CO2-C emission. CO2-C emission occurs until the substrate ends, after which the microorganisms die and CO2-C values decrease. Dead microfauna then become a source of energy for resistant microorganisms, which multiply and generate new peaks of CO2-C emission (Carvalho et al. 2008; Pinto et al. 2018). Such a pattern explains the peaks in CO2-C observed between days 7 and 13 of incubation, as also reported by Loss et al. (2013; 2014) and Rosset et al. (2019). In the present study, microbial activity stabilized from the 15th day onward (Fig. 1).
Total CO2-C values were highest at areas with the highest macroaggregate proportion, MWD values (Table 2), and TC and LOM contents (Table 3). Such findings agree with those of Loss et al. (2013; 2014), Wade et al. (2018), and Rosset et al. (2019). Total CO2-C values of CA, 6NT, and 18NT samples (Fig. 1) suggest that there is a direct relationship between LOM contents and CO2-C accumulation. According to Loss et al. (2014), greater LOM availability in aggregates contributes to microbial activity, leading to greater CO2-C accumulation.
Several studies have investigated the relationship between CO2-C and POXC under different soil and environmental conditions (Hurisson et al. 2016; Wade et al. 2018). Culman et al. (2013) analyzed a Typic Hapludalf cropped with maize and managed under different systems in the United States and observed that POXC was mostly influenced by amendment with organic matter (manure). On the other hand, CO2-C was more influenced by rotation of annual crops with cover crops. These results help explain the lack of differences in POXC contents as well as the significant variations in total CO2-C between management systems (Tables 3 and 5).
Dissimilarity between evaluated areas
Multivariate techniques were used to identify similarity patterns for the analyzed data. Such techniques have been applied in previous research (Silva et al. 2018; Rosset et al. 2019; Assunção et al. 2019; Briedis et al. 2020), contributing to understanding the patterns of variations between areas in association with statistical tests. Duval et al. (2018) used PCA to investigate associations between management systems and fertilizer applications; the authors observed that soil attributes had low sensitivity to management practices compared with soil quality indices.
PCA showed that aggregation status, more labile SOM compartments (LOM and total CO2-C), and SIMWD and LI indices were directly related to the management practices adopted in CT and NT. Total SOM compartment (TC), POXC values, and their respective indices were more associated with CA (more stable and balanced environment). The variables that contributed most to discriminating areas in the 0.00–0.10 m layer were aggregate classes (A1–A6), MWD, TC, C/N ratio, SIMWD, and CPITC (Figs. 2 and 3). Such findings highlight the importance of assessing aggregation and SOM compartmentalization to determine the quality of anthropized environments.
The reduction in dissimilarity between CT and NT areas with depth (from ± 56 to ± 40 %), as shown by hierarchical cluster analysis (Fig. 4), suggests that CT management practices exert stronger effects on the surface layer (0.00–0.05 m), resulting in greater dissimilarity to NT systems, mainly in surface.