Next Article in Journal
Genetic Parameters, Prediction of Gains and Intraspecific Hybrid Selection of Paspalum notatum Flügge for Forage Using REML/BLUP
Next Article in Special Issue
Impact of the Cultivation System and Pre-Sprouting of Selected Potato Cultivars on the Physico-Chemical Properties and Enzymatic Activity of Soil in the Conditions of Central-Eastern Poland
Previous Article in Journal
Prediction of Cultivation Areas for the Commercial and an Early Flowering Wild Accession of Salvia hispanica L. in the United States
Previous Article in Special Issue
Application of Conservation Tillage in China: A Method to Improve Climate Resilience
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Organic Matter and Aggregate Stability in Soybean, Maize and Urochloa Production Systems in a Very Clayey Soil of the Brazilian Savanna

by
Juslei Figueiredo da Silva
1,
Miguel Marques Gontijo Neto
2,
Gustavo Ferreira da Silva
1,*,
Emerson Borghi
2 and
Juliano Carlos Calonego
1
1
Department of Crop Science, School of Agriculture, São Paulo State University (UNESP), Botucatu 18610-034, Brazil
2
Brazilian Agricultural Research Corporation (EMBRAPA), EMBRAPA Maize and Sorghum, Sete Lagoas 35701-970, Brazil
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(7), 1652; https://doi.org/10.3390/agronomy12071652
Submission received: 2 June 2022 / Revised: 29 June 2022 / Accepted: 29 June 2022 / Published: 11 July 2022
(This article belongs to the Special Issue Effects of Tillage, Cover Crop and Crop Rotation on Soil)

Abstract

:
Forage grasses cultivation in production system with soybean and maize is an alternative to improve tropical weathered soils quality in Brazil. The aim of the study was to evaluate the effects in the production systems involving cultivation of Urochloa brizantha cv. Piatã, in monoculture or in succession with soybean and maize crops, on organic matter and structuring of soil in Brazilian savanna. The experiment was implemented in the 2010/2011 season. The treatments consisted of nine production systems and a native forest (savanna) as a reference area. In March 2017, soil sampling was carried out for C and N analysis, physical and chemical fractionation of SOM and aggregate stability. Production systems influenced total organic carbon (TOC) and aggregate stability, mainly in the surface layers, leading to changes in SOM quality. TOC was 31% lower in monoculture soybean production system, when compared to native savanna area, in the 0.00–0.20 m layer. The agricultural production systems influence organic matter quality and soil aggregates stability. For the Brazilian savanna conditions, grain cultivation systems under no-tillage that integrate Urochloa brizantha cv. Piatã contribute to the soil quality improvement. Soybean monoculture generally provides worse soil quality indices compared to other agricultural production systems.

1. Introduction

Conservation agriculture has been based on the search for soil quality improvement aiming at the sustainability of the production system, maintaining and/or increasing grain and forage yield [1,2]. A no-tillage system (NT) is the main agricultural management technique related to soil conservation, and its introduction into agriculture was one of the greatest advances in the Brazilian production process [3]. With the NT, it was possible to expand the cultivation to the tropical region in Brazil, and it was possible to increase the efficiency of the second crop production (fall–winter crop). In the Brazilian savanna, the main challenge for the establishment of NT is the difficulty in maintaining plant residues on the soil surface due to high temperatures and the dry period, which makes it difficult to implement an off-season or winter crop, compromising crop rotation and increasing soil organic matter in these regions [4].
Covering species, especially grasses, play an important role in the search for technical and economic viability of NT in different regions of Brazil, especially in tropical regions, due to the difficulty of production in the off-season and the rapid decomposition of the accumulated biomass [5,6], leaving the soil unprotected [7]. Intercropping of perennial forages with maize is an alternative for establishing ground cover crops, with the objective of producing maize and soybean grains in NT, keeping the soil permanently covered [8,9,10].
In addition to the high production and persistence of straw, perennial forages have been considered to improve the soil environment. The high contribution of aerial and root residues results in an increase in the soil organic matter content (SOM), which contributes to the formation of stable aggregates, due to the cementing and binding action that SOM exerts on the soil mineral particles [7,11,12,13]. Salton et al. [14] observed that systems with the use of forages for straw production, compared to those exclusively with crops, showed higher levels of C in the soil, which must be associated with the high input of plant material commonly provided by pastures.
SOM has an important role in the physical, chemical, and biological soil properties and, as a result, the study of its labile and recalcitrant fractions has been used as a soil quality indicator [15,16]. SOM analysis does not always detect the effects promoted by changes in management systems. One of the ways to improve the understanding of this issue is through the study of total organic carbon compartments, which are more sensitive to soil management and is characterized one of the better indicators of this dynamics [17].
Thus, the aim of this research was to evaluate the aggregates stability and soil organic carbon compartments in different production systems conducted in a long-term experiment under the no-tillage (NT) in the Brazilian savanna, involving the cultivation of Urochloa brizantha cv. Piatã in monoculture or in succession with soybean and maize crops.

2. Materials and Methods

2.1. Site Description, History and Experimental Area Characterizationubsection

The experiment was carried out in an experimental area belonging to Embrapa Maize and Sorghum in the Sete Lagoas, Minas Gerais, Brazil (19°28′ S; 44°15′ W and 732 m). Soil of the experimental area is a Typic Haplorthox [18], very clayey texture (Table 1), subdeciduous cerrado phase and smooth wavy relief [19]. Before the implementation of this experiment, the area had been already conducted in NT since 1995 with interspersed sowing of maize and soybeans. Table 1 shows the granulometric and chemical characteristics of the experimental area in a collection carried out in the area adjacent to the experiment (native savanna) in March 2017.

2.2. Treatments and Experimental Design

Experimental design was carried out in completely randomized blocks, with three replications. Treatments consisted of nine cropping systems, plus a control with native forest (Table 2 and Figure 1). Each plot was 12 m wide and 14 m long (168 m2) with lanes between plots and between blocks of 2.5 m. Production systems were implemented in the 2010/11 crop year and conducted for seven years. Before the implementation of this experiment, the area had been already managed in NT since 1995, with intercrop sowing of maize and soybeans in the summer and fallow in the off-season.

2.3. Conducting the Experimental Area

Sowing and management of the treatments, with the exception of system 4, Piatã pasture (Urochloa brizantha cv BRS Piatã), were carried out mechanically from the first year of implantation, with soybean and maize being sown with 0.5 m spacing between rows. For maize and soybean crops, an average of 3.5 and 15 seeds m−1 were used over the years, respectively, with 360 kg ha−1 of 08-28-16 NPK for maize and 315 kg ha−1 of formulated 06-30-16 NPK for soybeans.
In the V2 stage of maize development, topdressing fertilization was carried out in maize with 112.5 kg ha−1 of N via urea. In the maize plots intercropped with grass, 4.0 kg ha−1 of viable pure seeds of U. brizantha cv. Piatã were treated with fipronil and mixed with fertilizer at the time of maize sowing. Management was always carried out according to the need and current recommendations.
In the plots with Piatã grass, throughout the years, the management was carried out simulating hay production; that is, when the plants, for the most part, reached 0.5 m in height, the cut was carried out at 0.2 m in height, and the material was then removed from the area.

2.4. Soil Sampling

Samples were sampled in March 2017. For analysis of total organic carbon (TOC) and total nitrogen (TN) and the physical and chemical fractionation of carbon, the deformed soil samples were collected using a soil probe at three random points in each plot at depths of 0.00–0.10, 0.10–0.20 and 0.20–0.40 m. In the laboratory, the samples were air-dried, crushed and passed through a 2 mm sieve to obtain air-dried soil.
For the aggregate stability analysis, trenches measuring 0.30 m × 0.30 m × 0.30 m were opened in each experimental plot. In each trench, at a depth of 0.00–0.10 m, two monoliths were removed (one on each side of the trench) with approximate dimensions of 0.10 m in height by 0.20 m in length and 0.10 m in width, which were then placed in plastic containers to avoid crumbling.

2.5. Soil Analysis

2.5.1. Total Organic Carbon (TOC), Total Nitrogen (TN) and Soil Carbon Fractionation

The TOC and NT analyses were performed in an automatic elemental analyzer (LECO-TruSpec1 CHN, Leco Corp., St. Joseph, MI, USA) using 0.2 g of soil.
Soil organic matter physical fractionation was performed according to the method proposed by [20]. A 20 g soil sample (air-dried and sieved through a 2 mm mesh) was homogenized in an 80 mL sodium hexametaphosphate solution (5 g L−1). The samples were horizontally shaken for 15 h and sieved through 0.053 mm mesh sieves using deionized water. The retained material was transferred to aluminum dishes and dried at 45 C in a forced-air oven to a constant mass, after which it was ground with a porcelain mortar and pestle and homogenized. The particulate organic carbon (POC) was measured via dry combustion with an elemental analyzer, which is defined as the N and C content in the 0.053–2 mm soil fraction. The mineral-associated organic carbon (MOC) was calculated by the difference between TOC and POC. From the results of the physical fractionation, the carbon management index (CMI) was calculated [21]. This index is a relative measure of alterations caused by soil management, comparing it with a situation considered original or ideal. In this research, the native savanna (NS) was taken as the original situation (CMI = 100).
To obtain the CMI, the carbon pool index (CPI) and the TOC lability index (LI) are needed (Equations (1)–(3)):
CMI = CPI × LI × 100
where CPI is the carbon pool index and LI is the lability index.
The CPI and the LI were calculated as follows:
CPI = TOCtreatment/TOCreference
LI = Ltreatment/Lreference
where L refers to the C lability, calculated according Equation (4):
L = POC/MOC
The chemical fractionation of soil organic matter was determined according to the method suggested by the International Humic Substances Society and described by [22]. Fulvic acids (FA), humic acids (HA) and humin (HUM) were obtained based on the solubility in acids and alkalis. After the separation of each fraction, the levels of C were determined by obtaining CHum, CHA and CFA.

2.5.2. Aggregates Stability in Water

The soil monoliths were initially sieved using 8 and 4 mm sieves. Above of the 8 mm sieve, the monoliths were manually fragmented by the planes of weakness, using moderate force to avoid grinding all the soil. The material retained in the 4 mm sieve was used for the aggregate stability test [23]. Water stability of the aggregates was estimated by wet sieving of 4–8 mm aggregates by a set of 2.00, 1.00, 0.50, 0.25 and 0.105 mm sieves [23] using a vertical oscillation device [24]. Based on the distribution of aggregates in each sieve, the percentage of aggregates retained on the 2 mm sieve (%Aggregation > 2mm), mean weight-diameter (MWD), geometric mean diameter (GMD) and aggregate stability index were calculated. (ASI) according to [23].

2.6. Data Analysis

Data were submitted to analysis of variance (ANOVA) using the F test and the means compared using the t test (LSD-Student) at 5% error probability.

3. Results

3.1. Total Organic Carbon (TOC), Total Nitrogen (NT) and Soil Carbon Physical Fractionation

In NS, the highest TN levels and TOC (Table 3) occurred in the layers from 0.00 to 0.10 m and 0.10 to 0.20 m (p < 0.05); however, for TN, there were no differences (p > 0.05) between NS and M and MP/P/P in the 0.00–0.10 m layer and between NS and MP in the 0.10–0.20 m layer (Table 3). For TOC, in the 0.10–0.20 m layer, there was no difference (p > 0.05) between NS with the treatments MP and MP/S (Table 3), and thus they present higher levels of TOC but without differing from P and MP/P.
Among the systems studied, MP/P, MP/S, S/MP/P/P and M presented higher TOC content in relation to the S system in the 0.00–0.10 m layer. For the 0.10–0.20 m layer, the MP and MP/S systems did not differ from the NS and thus showed higher TOC contents but without differing from P and MP/P. In the 0.20–0.40 m layer, there were no differences between treatments for both TOC and TN (Table 3).
In the TOC physical fractionation, the main differences occurred in the POC (Table 3). In the 0.00–0.10 m layer, the highest content occurred in the NS (17.4 g kg−1), and the lowest levels occurred in the S and MP/P/S treatments (7.2 and 8.4 g kg−1, respectively). In the 0.10–0.20 m layer, the highest POC content occurred in the MP treatment and in the 0.20–0.40 m layer in the P, MP/S, M/P/P and MP/P/S treatments. As for MOC (Table 3), there was no significant difference (p > 0.05) between treatments in the surface layer and at 0.20–0.40 m (Table 3). However, in the 0.10–0.20 m layer, the highest content was found in the NS, without differing (p > 0.05) from the MP and MP/S treatments. On the other hand, the S/MP/P/P system promoted the lowest MOC content, differing from the MP and MP/S systems, in addition to NS (p < 0.05).

3.2. Soil C/N Ratio, Carbon Pool Index (CPI), Lability (L), Lability Index (LI) and Carbon Management Index (CMI)

In the 0.0–0.10 and 0.20–0.40 m layers, there was no difference (p > 0.05) between treatments for the soil C/N ratio, with the average of these layers being 16.73 and 17.69, respectively. In the 0.10–0.20 m layer, areas with P, MP/S and MP/P showed a higher C/N ratio than NS, which did not differ (p > 0.05) from the other treatments (Table 4).
In general, an increased CPI was observed at the 0.20–0.40 m depth for the cropping systems. At a 0.0–0.10 m depth, all treatments were below 1.00, the standard value of the reference area (NS); however, the MP/P, MP/S and S/MP/P/P systems obtained indices closer to the EC, differing (p < 0.05) from S and MP/P/S with lower CPI. For the 0.10–0.20 m layer, it was found that MP and MP/S did not differ (p > 0.05) from the NS, and in the 0.20–0.40 m layer, there was no difference (p > 0.05) between the evaluated treatments (Table 4).
As for the C lability (L), in the 0.0–0.10 m layer, the L index was higher in NS, without differing from the MP/P/P system (Table 4). Treatments P, MP/P and MP/P/P resulted in higher L compared to the S system, which showed lower lability due to its lower POC content (Table 3). In the 0.10–0.20 m layer, in treatments with MP and P, the greatest soil carbon lability was observed in relation to the other cropping systems, which did not differ from the NS. At 0.20–0.40 m, it was observed that L in treatments with M, MP and S did not differ from NS; however, higher L was observed in P, MP/S, MP/P/P and MP/P/S (Table 4).
The LI increased in depth, due to the reduction in L values in NS and in the treatments in general. In the 0.0–0.10 m layer, the LI in the MP/P/P system did not differ from the NS, being, therefore, closer to 1.00. It was also verified that the LI in the P and MP/P treatments were superior to area with S, where the lowest LI was observed, as well as verified for L. In the 0.10–0.20 m layer, the MP and P treatments presented the highest LI. In this layer, the lowest LI was obtained in the M system, which did not differ (p > 0.05) from S, MP/S, S/MP/P/P and NS. At 0.20–0.40 m, the LI results followed the same behavior observed for L described previously (Table 4).
The CMI values were higher with the increasing soil depth for the nine cropping systems evaluated, adopting 100% for the reference area (NS) (Table 4). In the 0–0.10 m layer, the CMI with cropping systems were lower than in the reference area. In this layer, treatment S resulted in the lowest CMI (34%) without differing from MP/P/S (46%), and the closest management to NS was MP/P (66%). The CMIs at 0.10–0.20 m with the MP and P treatments were higher than NS, 142% and 125%, respectively. The CMI with MP/P and MP/P/S did not differ from NS. In this layer, the lowest CMI occurred in areas with M (58%) and S (65%) systems, without differing from MP/S and S/MP/P/P. In the 0.20–0.40 m layer, all production systems showed a CMI greater than 100 (NS), with the highest CMI observed in areas cultivated with MP/P/S (235%), MP/S (230%), P (222%) and MP/P/P (212%) (Table 4).

3.3. Soil Organic Matter Chemical Fractionation

The NS area showed the highest CHum in the surface layer of the soil (18.12 g kg−1) (Table 5). Still in this layer, the MP/P system resulted in a higher CHum content compared to M, with the other treatments being statistically similar to both. In the 0.10–0.20 m layer, NS also showed higher CHum (14.67 g kg−1) but did not differ from MP and P. At 0.20–0.40 m, the MP/P/S treatment provided the highest CHum (13.01 g kg−1), however, without differing from NS and MP, MP/S and S/MP/P/P, and in the S system, there was the lowest CHum (8.14 g kg−1), without differing from M and MP/P/P (Table 5).
For C contents in humic acid fraction (CHA), there was only a difference between treatments in the 0.0–0.10 m and 0.10–0.20 m layers. In the 0.0–0.10 m layer, higher levels were observed in NS (7.15 g kg−1) but without differing from MP/P, MP/P/P and M (p > 0.05). In the 0.10–0.20 m layer, the highest CHA content occurred in MP/S (8.37 g kg−1) but without differing (p > 0.05) from MP/P (Table 5).
The CFA contents in the soil surface layer were higher in M and MP/S (7.85 and 8.04 g kg−1, respectively), without differing (p > 0.05) from MP and NS. In the 0.10–0.20 m layer, the highest CFA occurred in the MP/P/S treatment (7.36 g kg−1), being a similar value to that obtained in NS (6.48 g kg−1). The lowest levels in this fraction occurred in the S treatment (3.44 g kg−1), without differing from M, MP, MP/S, MP/P and MP/P/P. At 0.20–0.40 m, a lower CFA value occurred in NS (3.38 g kg−1) and in the S system (3.65 g kg−1), without differing (p > 0.05) from the M, MP/P/P and MP/P/S cultivation systems, while MP/P resulted in the highest CFA content (5.87 g kg−1) without differing from MP, P, MP/S and S/MP/P/P (Table 5).
The S and MP/P treatments resulted the highest CHA/CFA ratios in the 0.0–0.10 m layer (1.31 and 1.24, respectively) (Table 5). At 0.10–0.20 m, the highest ratios were obtained in the MP, S, MP/S and MP/P systems without differing from M and MP/P/P, and at 0.20–0.40 m, the NS had the highest CHA/CFA ratio but was not different from M, S and MP/P/S.

3.4. Aggregate Stability

The NS soil presented greater aggregates stability, with higher values of %Aggregation > 2 mm (98.1%), GMD (4.76 mm) and MWD (4.91 mm); however, it did not differ (p < 0.05) from the P, MP/P and MP/P/P systems (Table 6). On the other hand, lower soil structuring was observed in the S treatment, with values of %Aggregation > 2mm, GMD and MWD of 41, 53 and 35%, respectively. These values were lower than those found in NS, however, without differing (p < 0.05) from those obtained in MP/S and S/MP/P/P.
The ASI values did not differ between the cropping systems (p < 0.05), with an overall average of 96.52%.

4. Discussion

TOC, TN, Physical Fractionation and Carbon Management Index of Soil

The MP treatment provided higher levels of TN and TOC in the 0.10–0.20 m layer, similar to those obtained in NS (Table 3). This system is implemented annually in the summer, and after the maize harvest, the forage remains in the area; thus, the soil remains covered throughout the year, favoring the development of forage roots and organic material accumulation in subsurface, evidenced by the high POC values (Table 3), which are highly related to recent plant materials deposition [20]. In systems with one or two years of U. brizantha cv. Piatã, the forage remains in the area without receiving maintenance fertilization; thus, the nutrients supply via replacement is lower, compromising dry matter productivity. In the MP system, forage is desiccated annually for maize sowing, and these residues are important for increasing the N and C levels in the soil [25,26].
Considering the superficial layers average (0.0–0.20 m), where differences for TOC were significant, it is observed that in area with the P treatment there was a reduction of 22% in the TOC content in relation to the NS. However, in the MP, MP/S and MP/P systems, this average proportion is around 15%, suggesting that these treatments were more efficient in the maintenance and/or accumulation of TOC in the soil. This result may be due to the plant residues’ input and management with annual fertilization (MP and MP/S), or every two years in the case of MP/P. For the other treatments, the mean TOC was between 24 and 26%. The CPI establishes a relationship between the TOC of the cultivated areas and the TOC of the reference area (NS). The treatments MP/P, MP/S and S/MP/P/P in the 0–0.10 m layer and MP and MP/S in the 0.10 to 0.20 m layer showed CPIs closer to 1. In the 0.20–0.40 m layer, all treatments showed a CPI equal to or greater than 1. It can be inferred in these systems with higher CPI that there was an increase in, or at least maintenance of, C in the soil over the years with no-tillage [27]. Thus, the contribution of residues from crop and forage straw over the years provided increasing increases in TOC, mainly in depth.
In treatments referring to monocultures (M and S), there is a lower proportion of labile C being added compared to other systems with the presence of forage, in which the highest POC levels were observed in the 0.20–0.40 m layer. This result may be related to the greater presence of residues from roots that are still in the initial decomposition stage [28]. Forage roots have a well-developed root system and are able to reach greater depths in the soil, increasing POC levels [29,30]. Systems with the use of forages for straw production, compared to those exclusively with crops, showed higher C levels, which must be associated with the high plant material input commonly provided by pastures [14].
For MOC contents, there was no significant difference (p > 0.05) between treatments in the surface layer and at 0.20–0.40 m (Table 3). Changes in this carbon fraction occur more slowly and gradually, and it represents the stable carbon in the soil [31,32,33]. However, in the 0.10–0.20 m layer, high MOC levels in the MP and MP/S systems were similar to those observed in the NS, which can be explained by the constant C supply by cultivated species and also because they were treatments that received fertilization annually, or with N input via soybean in the case of the MP/S system. Increases in C levels in more stable SOM fractions depend on the balanced relationship between C and N inputs, in addition to the availability of other nutrients [16]. Biomass-associated N availability has also been reported as another strategy to increase stabilized SOM content, as it increases microbial activity [34]. It is understood that most stable SOM compounds have been transformed through N-accelerated biological activities [35].
The lower lability (L) in the S system, in the 0.0–0.10 m layer, occurs due to the lower labile C (POC) content obtained in this treatment. Soybean crop provides less plant residues, and these residues have a low C/N ratio; therefore, it decomposes faster than Piatã residues. As the soil collection took place in March, that is, approximately one year after the previous year’s soybean harvest, it is likely that a large part of labile C was lost. Salton et al. [36] studying SOM and soil aggregation in crop-pasture rotation, observed L lower than the values found in the present research, both for the area with native vegetation, ranging from 0.08 to 0.28, as in the area with NT farming, from about 0.08 to 0.26, in the 0.0–0.20 m layer. The divergence in the results is related to the lower TOC levels verified by Salton et al. [36] and shorter conduct time in the experimental area.
However, as observed by Salton et al. [36], there was greater L in the soil surface layer, which can be attributed to the greater amount of plant residues and POC concentration. The LI increased in depth, due to the reduction in L values in the reference area and, in general, in the treatments. This observation is also verified in studies by [36,37]. Although there is a difference in C levels between the treatments, through the LI, it is possible to observe whether the L of these is close to or below the natural condition without anthropic intervention, so that LI values close to or above 1.00 indicate that the production system employee is providing labile C addition similar to or greater than the original condition, respectively.
Therefore, the greatest contributions to the addition of labile C by production systems occurred at the deepest layer (0.20 to 0.40 m), where all cultivation systems that use Piatã grass presented LI > 1, with emphasis on the P, MP/P/S, MP/S and MP/P/P systems with ILs of 2.06, 2.01, 1.94 and 1.94, respectively, most likely due to the addition of labile C by grass roots. These values were lower than those obtained by [37] with rotation systems involving Urochloa ruziziensis, soybean and sorghum, in which the authors obtained LI between 1.1 and 4.9 in the 0.0–0.10 m layer and above 5.0 in the 0.10–0.20 m layer.
In general, the low CMI values in the 0.0 to 0.10 m layer indicate the need to improve cropping systems in order to match the original condition of the soil in a native forest condition. This is a difficult condition to be obtained in the region where the research was conducted due to the rapid degradation of labile organic matter in a tropical climate because of the higher decomposition rates favored by high temperatures characteristic of these regions [38]. However, at greater depths in the profile, it is already possible to obtain CMI > 100 due to the constant contribution of root biomass, lower micro-organisms activity in deeper layers and the greater aggregates preservation due to lower anthropic action. Production systems with higher CMI values show the ability to promote the sustainability in tropical regions by maintaining C in the agricultural system [39,40].
Thus, higher CMI values mean a positive effect of soil use and management practices on organic matter content [41,42], while the lowest values indicate that the C compounds are being degraded [21], reflecting the lower quantity and quality of residues added by the production systems, that is, lower labile C content [43]. The CMI represents a measure of sustainability of different production systems or land uses and can be used to compare the changes that occur in labile C and TOC contents as a result of agricultural practices [44]. Thus, low CMI values with monocultures S and M obtained in layers 0 to 0.10 and 0.10 to 0.20 m, respectively, indicate the lower capacity of these crops to promote sustainability to production systems.
Stable C presence in the soil is also an indication of the quality and sustainability of production systems, which generally suffer the greatest interference in the more superficial layers. In the layer from 0 to 0.10 m, there were the greatest differences in CHum between the production systems and the NS. Humin fraction makes up the most stable carbon and is very resistant to decomposition [45], generally representing most of humidified C in tropical soils [46], as verified in this research, since the CHum represented on average about 53% of humic substances in cultivated areas and 57% in the NS. The C reduction in this fraction of SOM indicates that soil management practices led to the loss of part of stable C over the cultivation years because of the native forest’s conversion into arable land.
The lower CFA levels in soybean monoculture, mainly in the layer from 0 to 0.10 m, is related to the lower POC levels with this production system, due to the lower biomass contributions. The POC detects the portion of SOM that has been newly deposited and is less humified, justifying the lower CFA levels in this treatment, since, among the humic substances evaluated, fuvic acid has a lower polymerization and humification degree [47].
In this study, it was found that the CHA/CFA ratio varied between soil layers and cropping systems (Table 5). Humic acid fractions were not always higher than those of fuvic acid, although, for the most part, the CHA/CFA values are close to or above 1.00, which indicate areas where there is a predominance of more stable and better-quality organic material [38,48]. However, it is noted that treatment S showed a CHA/CFA ratio greater than 1.00 in the entire profile, but this did not lead to higher CHA and CHum levels compared to the other treatments, indicating that the highest CHA/CFA ratio in S is actually due to a lower C labile proportion (Table 3), from which the humified SOM is formed, leading to lower CFA levels in this treatment.
As for soil structure, the area without anthropic interference (NS) maintains its structure unchanged, and therefore, it has better soil aggregation indices. Production systems that involved the exclusive Piatã (P) growed and without soybean crop succession (M, MP, MP/P, MP/P/P) obtained, in general, better aggregation indices. Effective surface coverage by forages reduces or even prevents raindrop impact [49], promotes hydraulic roughness and reduces surface runoff [50,51], favoring the preservation of soil moisture and contributing to a more favorable environment for aggregation [52,53].
In the research, treatments that presented higher CMI (Table 4) also had better aggregation, except for the MP/S treatment. A possible explanation is the period with uncovered soil from the previous crop (2015/16), after the soybean harvest. It is observed that the cultivation time in the S treatment has a negative effect on the aggregation, and to a lesser extent, this effect occurs in other systems in which the area is left without cultivation in the fall–winter season (when there is soybean cultivation in the summer). This shows that the fact that soil remains part of year without vegetation implies soil destabilization over time because, in addition to soil protection, the vegetation acts in forming aggregates through the roots’ mechanical action and/or through the excretion of substances with cementing action [54], which can serve as a substrate for microorganisms, stimulating their activity and leading to the production of new cementing agents [55].
Larger aggregates can be formed around residues recently added to the soil, which make up the particulate organic matter [56]. These residues allow the macroaggregates’ formation, as they are a source of labile C for microbial activity and for the production of binding agents [57]. Therefore, particulate carbon probably favored aggregates formation and MWD in NS, P, MP/P and MP/P/P, given that NS and MP/P/P, followed by P and MP/ P, showed higher carbon lability (L) in relation to the other systems (Table 4). In addition, in NS, MP/P and MP/P/P, there was a higher POC content in the soil surface layer (Table 3).
In addition to effect of labile C (POC), the aggregation in treatments NS, MP/P and MP/P/P may be related to the higher CHA values in these treatments (Table 5) and to CHum values in the MP/P system in the surface layer. Humified organic matter has the ability to associate with soil mineral particles, forming clay–metal–humic bonds, which contributes to soil aggregation [58].
The greater aggregation in NS and P areas can also be attributed, in part, to chemical processes, since these areas were not chemically corrected; thus, an acidic condition with high aluminum contents is expected, which may favor aggregation, as polyvalent cations have aggregating action in the soil [59].

5. Conclusions

Soybean monoculture generally provides worse soil quality indices compared to other agricultural production systems.
The agricultural production systems influence organic matter quality and soil aggregates stability. Maize intercropped with U. brizantha cv. Piatã in annual succession and Maize + Piatã intercropped and Soybean in annual rotation are better options to improve the organic matter quality in Brazilian savanna conditions and very clayey soil. Agricultural production systems involving maize cultivation provide better soil structuring than those with soybean.

Author Contributions

Conceptualization M.M.G.N., E.B., J.F.d.S. and J.C.C.; methodology, J.F.d.S., J.C.C.; M.M.G.N. and E.B.; formal analysis, J.F.d.S.; investigation, J.F.d.S. and J.C.C.; writing—original draft preparation, J.F.d.S. and G.F.d.S.; writing—review and editing, J.C.C., M.M.G.N. and E.B.; supervision, J.C.C.; funding acquisition, J.F.d.S. and J.C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Coordination for the Improvement of Higher-Level Personnel (CAPES), grant number 001; the Agrisus Foundation, grant number 2038/17; and the National Council for Scientific and Technological Development (CNPq) for granting an award for excellence in research to the fifth author, grant number 301596/2019-0.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to acknowledge the College of Agricultural Sciences at São Paulo State University (UNESP) for supporting the development of this research, and to Brazilian Agricultural Research Corporation (EMBRAPA) for conducting the experimental area.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carducci, C.E.; de Oliveira, G.C.; Silva, É.A.; de Oliveira, L.M.; Barbosa, S.M. Retenção de água do solo sob sistema conservacionista de manejo com diferentes doses de gesso. Rev. Ciências Agrárias Amazon. J. Agric. Environ. Sci. 2015, 58, 284–291. [Google Scholar] [CrossRef]
  2. Da Silva, P.C.G.; Tiritan, C.S.; Echer, F.R.; dos Santos Cordeiro, C.F.; Rebonatti, M.D.; dos Santos, C.H. No-tillage and crop rotation increase crop yields and nitrogen stocks in sandy soils under agroclimatic risk. Field Crop. Res. 2020, 258, 107947. [Google Scholar] [CrossRef]
  3. Kassam, A.; Friedrich, T.; Derpsch, R. Global spread of Conservation Agriculture. Int. J. Environ. Stud. 2019, 76, 29–51. [Google Scholar] [CrossRef]
  4. Alvarenga, R.C.; Cabezas, W.A.L.; Cruz, J.C.; Santana, D.P. Plantas de cobertura de solo para sistema plantio direto. Inf. Agropecuário 2001, 22, 25–36. [Google Scholar]
  5. Grassmann, C.S.; Mariano, E.; Rocha, K.F.; Gilli, B.R.; Rosolem, C.A. Effect of tropical grass and nitrogen fertilization on nitrous oxide, methane, and ammonia emissions of maize-based rotation systems. Atmos. Environ. 2020, 234, 117571. [Google Scholar] [CrossRef]
  6. Rigon, J.P.G.; Calonego, J.C.; Capuani, S.; Franzluebbers, A.J. Soil organic C affected by dry-season management of no-till soybean crop rotations in the tropics. Plant Soil 2021, 462, 577–590. [Google Scholar] [CrossRef]
  7. Chioderoli, C.A.; de Mello, L.M.M.; Grigolli, P.J.; Furlani, C.E.A.; Silva, J.O.R.; Cesarin, A.L. Atributos físicos do solo e produtividade de soja em sistema de consórcio milho e braquiária. Rev. Bras. Eng. Agrícola Ambient. 2012, 16, 37–43. [Google Scholar] [CrossRef]
  8. Dalla Côrt, A.S.; Feitosa, P.B.; Pacheco, L.P.; Greco, T.M.; Silva, I.D.F.; Souza, E.D.; Santos, F.S.; Petter, F.A.; Crusciol, C.A.C. Accumulation and efficiency of nutrient use in crop systems in second crop under no-tillage. Pesqui. Agropecuária Bras. 2021, 56, e01879. [Google Scholar] [CrossRef]
  9. Lauriault, L.; Guldan, S.; Popiel-Powers, F.; Steiner, R.; Martin, C.; Heyduck, R.; Falk, C.; Petersen, M.; May, T. Relay Intercropping with Cover Crops Improved Autumn Forage Potential of Sweet Maize Stover. Agriculture 2018, 8, 103. [Google Scholar] [CrossRef] [Green Version]
  10. Bybee-Finley, K.; Ryan, M. Advancing Intercropping Research and Practices in Industrialized Agricultural Landscapes. Agriculture 2018, 8, 80. [Google Scholar] [CrossRef] [Green Version]
  11. Braida, J.A.; Bayer, C.; Albuquerque, J.A.; Reichert, J.M. Matéria orgânica e seu efeito na física no solo. In Tópicos em Ciência do Solo; Klauberg Filho, O., Mafra, A.L., Gatiboni, L.C., Eds.; Sociedade Brasileira de Ciência do Solo: Viçosa, Brazil, 2011; pp. 221–278. [Google Scholar]
  12. Ferrari Neto, J.; Franzluebbers, A.J.; Crusciol, C.A.C.; Rigon, J.P.G.; Calonego, J.C.; Rosolem, C.A.; do Nascimento, C.A.C.; Ribeiro, L.C. Soil carbon and nitrogen fractions and physical attributes affected by soil acidity amendments under no-till on Oxisol in Brazil. Geoderma Reg. 2021, 24, e00347. [Google Scholar] [CrossRef]
  13. Momesso, L.; Crusciol, C.A.C.; Leite, M.F.A.; Bossolani, J.W.; Kuramae, E.E. Forage Grasses Steer Soil Nitrogen Processes, Microbial Populations, and Microbiome Composition in A Long-term Tropical Agriculture System. Agric. Ecosyst. Environ. 2022, 323, 107688. [Google Scholar] [CrossRef]
  14. Salton, J.C.; Mielniczuk, J.; Bayer, C.; Fabrício, A.C.; Macedo, M.C.M.; Broch, D.L. Teor e dinâmica do carbono no solo em sistemas de integração lavoura-pecuária. Pesqui. Agropecuária Bras. 2011, 46, 1349–1356. [Google Scholar] [CrossRef] [Green Version]
  15. Dhaliwal, S.S.; Naresh, R.K.; Mandal, A.; Singh, R.; Dhaliwal, M.K. Dynamics and transformations of micronutrients in agricultural soils as influenced by organic matter build-up: A review. Environ. Sustain. Indic. 2019, 1, 100007. [Google Scholar] [CrossRef]
  16. Raphael, J.P.A.; Calonego, J.C.; Milori, D.M.B.P.; Rosolem, C.A. Soil organic matter in crop rotations under no-till. Soil Tillage Res. 2016, 155, 45–53. [Google Scholar] [CrossRef] [Green Version]
  17. Dou, F.; Wright, A.L.; Hons, F.M. Sensitivity of Labile Soil Organic Carbon to Tillage in Wheat-Based Cropping Systems. Soil Sci. Soc. Am. J. 2008, 72, 1445–1453. [Google Scholar] [CrossRef]
  18. Soil Survey Staff. Keys to Soil Taxonomy, 12th ed.; USDA-Natural Resources Conservation Service: Washington, DC, USA, 2014; Volume 12.
  19. Santos, H.G.; Jacomine, P.K.T.; Anjos, L.H.C.; Oliveira, V.A.; Lumbreras, J.F.; Coelho, M.R.; Almeida, J.A.; Araujo Filho, J.C.; Oliveira, J.B.; Cunha, T.J.F. Sistema Brasileiro de Classificação de Solos, 5th ed.; Embrapa: Brasília, Brazil, 2018. [Google Scholar]
  20. Cambardella, C.A.; Elliott, E.T. Particulate Soil Organic-Matter Changes across a Grassland Cultivation Sequence. Soil Sci. Soc. Am. J. 1992, 56, 777–783. [Google Scholar] [CrossRef]
  21. Blair, G.; Lefroy, R.; Lisle, L. Soil carbon fractions based on their degree of oxidation, and the development of a carbon management index for agricultural systems. Aust. J. Agric. Res. 1995, 46, 1459–1466. [Google Scholar] [CrossRef]
  22. Swift, R.S. Organic matter characterization. In Methods of Soil Analysis: Part 3. Chemical Methods, 5.3; Soil Science Society of America Book, Series; Sparks, D.L., Page, A.L., Helmke, P.A., Loeppert, R.H., Soltanpour, P.N., Tabatabai, M.A., Johnston, C.T., Sumner, M.E., Eds.; Soil Science Society of America: Madison, WI, USA; American Society of Agronomy: Madison, WI, USA, 1996. [Google Scholar]
  23. Kemper, W.D.; Chepil, W.S. Size Distribution of Aggregates. In Methods of Soil Analysis: Part 1 Physical and Mineralogical Properties, Including Statistics of Measurement and Sampling; American Society of Agronomy: Madison, WI, USA, 1965; pp. 499–510. [Google Scholar]
  24. Yoder, R.E. A Direct Method of Aggregate Analysis of Soils and a Study of the Physical Nature of Erosion Losses 1. Agron. J. 1936, 28, 337–351. [Google Scholar] [CrossRef]
  25. De Queiroz Barcelos, J.P.; de Souza, M.; do Nascimento, C.A.C.; Rosolem, C.A. Soil acidity amelioration improves N and C cycles in the short term in a system with soybean followed by maize-guinea grass intercropping. Geoderma 2022, 421, 115909. [Google Scholar] [CrossRef]
  26. Canisares, L.P.; Rosolem, C.A.; Momesso, L.; Crusciol, C.A.C.; Villegas, D.M.; Arango, J.; Ritz, K.; Cantarella, H. Maize-Brachiaria intercropping: A strategy to supply recycled N to maize and reduce soil N2O emissions? Agric. Ecosyst. Environ. 2021, 319, 107491. [Google Scholar] [CrossRef] [PubMed]
  27. Paul, B.K.; Vanlauwe, B.; Ayuke, F.; Gassner, A.; Hoogmoed, M.; Hurisso, T.T.; Koala, S.; Lelei, D.; Ndabamenye, T.; Six, J.; et al. Medium-term impact of tillage and residue management on soil aggregate stability, soil carbon and crop productivity. Agric. Ecosyst. Environ. 2013, 164, 14–22. [Google Scholar] [CrossRef]
  28. Garcia, R.A.; Rosolem, C.A. Agregados em um Latossolo sob sistema plantio direto e rotação de culturas. Pesqui. Agropecuária Bras. 2010, 45, 1489–1498. [Google Scholar] [CrossRef] [Green Version]
  29. Oduor, C.O.; Karanja, N.K.; Onwonga, R.N.; Mureithi, S.M.; Pelster, D.; Nyberg, G. Enhancing soil organic carbon, particulate organic carbon and microbial biomass in semi-arid rangeland using pasture enclosures. BMC Ecol. 2018, 18, 45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Lei, L.; Thompson, J.A.; McDonald, L.M. Soil Organic Carbon Pools and Indices in Surface Soil: Comparing a Cropland, Pasture, and Forest Soil in the Central Appalachian Region, West Virginia, U.S.A. Commun. Soil Sci. Plant Anal. 2022, 53, 17–29. [Google Scholar] [CrossRef]
  31. Poeplau, C.; Don, A.; Six, J.; Kaiser, M.; Benbi, D.; Chenu, C.; Cotrufo, M.F.; Derrien, D.; Gioacchini, P.; Grand, S.; et al. Isolating organic carbon fractions with varying turnover rates in temperate agricultural soils—A comprehensive method comparison. Soil Biol. Biochem. 2018, 125, 10–26. [Google Scholar] [CrossRef] [Green Version]
  32. Haddix, M.L.; Gregorich, E.G.; Helgason, B.L.; Janzen, H.; Ellert, B.H.; Francesca Cotrufo, M. Climate, carbon content, and soil texture control the independent formation and persistence of particulate and mineral-associated organic matter in soil. Geoderma 2020, 363, 114160. [Google Scholar] [CrossRef]
  33. Lavallee, J.M.; Soong, J.L.; Cotrufo, M.F. Conceptualizing soil organic matter into particulate and mineral-associated forms to address global change in the 21st century. Glob. Change Biol. 2020, 26, 261–273. [Google Scholar] [CrossRef] [Green Version]
  34. Sarker, J.R.; Singh, B.P.; Dougherty, W.J.; Fang, Y.; Badgery, W.; Hoyle, F.C.; Dalal, R.C.; Cowie, A.L. Impact of agricultural management practices on the nutrient supply potential of soil organic matter under long-term farming systems. Soil Tillage Res. 2018, 175, 71–81. [Google Scholar] [CrossRef]
  35. Grandy, A.S.; Neff, J.C. Molecular C dynamics downstream: The biochemical decomposition sequence and its impact on soil organic matter structure and function. Sci. Total Environ. 2008, 404, 297–307. [Google Scholar] [CrossRef]
  36. Salton, J.C.; Mielniczuk, J.; Bayer, C.; Fabrício, A.C.; Macedo, M.C.M.; Broch, D.L.; Boeni, M.; Conceição, P.C. Matéria Orgânica do Solo na Integração Lavoura-Pecuária em Mato Grosso do Sul; Embrapa Agropecuária Oeste: Dourados, Brazil, 2005; Volume 29, p. 58. [Google Scholar]
  37. Rossi, C.Q.; Pereira, M.G.; Giácomo, S.G.; Betta, M.; Polidoro, J.C. Frações orgânicas e índice de manejo de carbono do solo em Latossolo Vermelho sob plantio de soja no cerrado goiano. Rev. Bras. Ciências Agrárias 2012, 7, 233–241. [Google Scholar] [CrossRef] [Green Version]
  38. Gmach, M.-R.; Dias, B.O.; Silva, C.A.; Nóbrega, J.C.A.; Lustosa-Filho, J.F.; Siqueira-Neto, M. Soil organic matter dynamics and land-use change on Oxisols in the Cerrado, Brazil. Geoderma Reg. 2018, 14, e00178. [Google Scholar] [CrossRef]
  39. Duval, M.E.; Martinez, J.M.; Galantini, J.A. Assessing soil quality indices based on soil organic carbon fractions in different long-term wheat systems under semiarid conditions. Soil Use Manag. 2020, 36, 71–82. [Google Scholar] [CrossRef]
  40. Guimarães, D.V.; Silva, M.L.N.; Beniaich, A.; Pio, R.; Gonzaga, M.I.S.; Avanzi, J.C.; Bispo, D.F.A.; Curi, N. Dynamics and losses of soil organic matter and nutrients by water erosion in cover crop management systems in olive groves, in tropical regions. Soil Tillage Res. 2021, 209, 104863. [Google Scholar] [CrossRef]
  41. Torres, J.L.R.; e Silva, V.R.; de Almeida Costa, D.D.; Pereira, M.G.; Assunção, S.A.; Gonzalez, A.P.; Pinto, L.A.; Loss, A. Aggregation and dynamics of soil organic matter under different management systems in the Brazilian Cerrado. Soil Res. 2021, 59, 715. [Google Scholar] [CrossRef]
  42. Gopinath, K.A.; Rajanna, G.A.; Venkatesh, G.; Jayalakshmi, M.; Kumari, V.V.; Prabhakar, M.; Rajkumar, B.; Chary, G.R.; Singh, V.K. Influence of Crops and Different Production Systems on Soil Carbon Fractions and Carbon Sequestration in Rainfed Areas of Semiarid Tropics in India. Sustainability 2022, 14, 4207. [Google Scholar] [CrossRef]
  43. Barbosa, L.R.; Souza, H.A.D.; Teixeira Neto, M.L.; Leite, L.F.C. Organic matter compartments in an Ultisol under integrated agricultural and livestock production systems in the Cerrado. Ciência Rural 2022, 52, 20200845. [Google Scholar] [CrossRef]
  44. Meena, V.S.; Mondal, T.; Pandey, B.M.; Mukherjee, A.; Yadav, R.P.; Choudhary, M.; Singh, S.; Bisht, J.K.; Pattanayak, A. Land use changes: Strategies to improve soil carbon and nitrogen storage pattern in the mid-Himalaya ecosystem, India. Geoderma 2018, 321, 69–78. [Google Scholar] [CrossRef]
  45. Corbeels, M.; Marchão, R.L.; Neto, M.S.; Ferreira, E.G.; Madari, B.E.; Scopel, E.; Brito, O.R. Evidence of limited carbon sequestration in soils under no-tillage systems in the Cerrado of Brazil. Sci. Rep. 2016, 6, 21450. [Google Scholar] [CrossRef] [Green Version]
  46. Benites, V.M.; Machado, P.L.O.A.; Madari, B.E.; Fontana, A. Fracionamento químico da matéria orgânica. In Manual de Métodos de Análises de Solo; Teixeira, P.C., Donagemma, G.K., Fontana, A., Teixeira, W.G., Eds.; Embrapa: Brasilia, Brazil, 2017; pp. 401–414. [Google Scholar]
  47. Guerra, J.G.M.; Santos, G.A. Métodos químicos e físicos. In Fundamentos da Matéria Orgânica do Solo: Ecossistemas Tropicais e Subtropicais; Santos, G.A., Silva, L.S., Canellas, L.P., Camargo, F.A.O., Eds.; Genesis: Porto Alegre, Brazil, 2008; pp. 185–199. [Google Scholar]
  48. Gazolla, P.R.; Guareschi, R.F.; Perin, A.; Pereira, M.G.; Rossi, C.Q. Frações da matéria orgânica do solo sob pastagem, sistema plantio direto e integração lavoura-pecuária. Semin. Ciências Agrárias 2015, 36, 693. [Google Scholar] [CrossRef] [Green Version]
  49. Cavalcante, D.M.; de Castro, M.F.; Chaves, M.T.L.; da Silva, I.R.; de Oliveira, T.S. Effects of rehabilitation strategies on soil aggregation, C and N distribution and carbon management index in coffee cultivation in mined soil. Ecol. Indic. 2019, 107, 105668. [Google Scholar] [CrossRef]
  50. Cantalice, J.R.B.; Silveira, F.P.M.; Singh, V.P.; Silva, Y.J.A.B.; Cavalcante, D.M.; Gomes, C. Interrill erosion and roughness parameters of vegetation in rangelands. CATENA 2017, 148, 111–116. [Google Scholar] [CrossRef]
  51. Zhao, C.; Gao, J.; Huang, Y.; Wang, G.; Zhang, M. Effects of Vegetation Stems on Hydraulics of Overland Flow Under Varying Water Discharges. Land Degrad. Dev. 2016, 27, 748–757. [Google Scholar] [CrossRef]
  52. Costa, N.R.; Andreotti, M.; Mascarenhas Lopes, K.S.; Yokobatake, K.L.; Ferreira, J.P.; Pariz, C.M.; Bonini, B.; dos Santos, C.; Longhini, V.Z. Soil Properties and Carbon Accumulation in an Integrated Crop-Livestock System Under No-Tillage. Rev. Bras. Ciência Solo 2015, 39, 852–863. [Google Scholar] [CrossRef] [Green Version]
  53. Arantes, A.C.C.; Cotta, S.R.; Conceição, P.M.D.; Meneghin, S.P.; Martinelli, R.; Próspero, A.G.; Boaretto, R.M.; Andreote, F.D.; Mattos, D., Jr.; Azevedo, F.A.D. Implication of Urochloa spp. Intercropping and Conservation Agriculture on Soil Microbiological Quality and Yield of Tahiti Acid Lime in Long Term Orchard Experiment. Agriculture 2020, 10, 491. [Google Scholar] [CrossRef]
  54. Bronick, C.J.; Lal, R. Soil structure and management: A review. Geoderma 2005, 124, 3–22. [Google Scholar] [CrossRef]
  55. Silva, I.R.; Mendonça, E.S. Matéria Orgânica do Solo. In Fertilidade do Solo; Novais, R.F., Alvarez, V.V.H., Barros, N.F., Ferreira, R.L., Cantarutti, R.B., Neves, J.C.L., Eds.; Sociedade Brasileira de Ciência do Solo: Viçosa, Brazil, 2007; pp. 275–374. [Google Scholar]
  56. Six, J.; Elliott, E.; Paustian, K. Soil macroaggregate turnover and microaggregate formation: A mechanism for C sequestration under no-tillage agriculture. Soil Biol. Biochem. 2000, 32, 2099–2103. [Google Scholar] [CrossRef]
  57. dos Reis Ferreira, C.; da Silva Neto, E.C.; Pereira, M.G.; do Nascimento Guedes, J.; Rosset, J.S.; dos Anjos, L.H.C. Dynamics of soil aggregation and organic carbon fractions over 23 years of no-till management. Soil Tillage Res. 2020, 198, 104533. [Google Scholar] [CrossRef]
  58. Tavares Filho, J.; Melo, T.R.D.; Machado, W.; Maciel, B.V. Structural changes and degradation of Red Latosols under different management systems for 20 years. Rev. Bras. Ciência Solo 2014, 38, 1293–1303. [Google Scholar] [CrossRef] [Green Version]
  59. Auler, A.C.; Pires, L.F.; Caires, E.F. Surface and incorporated liming effects on clay dispersion, water availability, and aeration capacity of a Dystrudept soil. Bragantia 2017, 76, 433–446. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Crop sequence representation referring to the 9 treatments in the 7 agricultural years between the 2010/11 and 2016/2017 seasons. Year 1: 2010/11 season; Year 2: 2011/12 season; Year 3: 2012/13 season; Year 4: 2013/14 season; Year 5: 2014/15 season; Year 6: 2015/16 season; Year 7: 2016/17 season.
Figure 1. Crop sequence representation referring to the 9 treatments in the 7 agricultural years between the 2010/11 and 2016/2017 seasons. Year 1: 2010/11 season; Year 2: 2011/12 season; Year 3: 2012/13 season; Year 4: 2013/14 season; Year 5: 2014/15 season; Year 6: 2015/16 season; Year 7: 2016/17 season.
Agronomy 12 01652 g001
Table 1. Chemical and granulometric soil analysis of the experimental area in March 2017.
Table 1. Chemical and granulometric soil analysis of the experimental area in March 2017.
Soil DepthPHPresin 1H + AlCaMgKSBCTCVSand 2Silt 3Clay 4
mCaCl2mg dm−3mmolc dm−3%g kg−1
0.00–0.204.05.01144081.24816230117113770
0.20–0.403.62.01231550.92114414110116774
1 P determination using ion exchange resin; 2 particles > 0.05 mm and <2 mm; 3 particles > 0.002 mm and <0.05 mm; 4 particles < 0.002 mm.
Table 2. Treatments description.
Table 2. Treatments description.
TreatmentsProduction SystemCode
T1Single Maize in annual succession (monoculture)M
T2Maize intercropped with U. brizantha cv. Piatã in annual successionMP
T3Soybean in annual succession (monoculture)S
T4U. brizantha cv. Piatã pasture implanted in single cultivationP
T5Maize + Piatã intercropped and Soybean in annual rotationMP/S
T6Sequence of Maize + Piatã intercropping and Piatã pasture in annual rotationMP/P
T7Sequence of Maize + Piatã intercropping and Piatã pasture for two yearsMP/P/P
T8Sequence of Maize + Piatã intercropping, Piatã pasture and Soybean in annual rotationMP/P/S
T9Sequence of Soybean, Maize + Piatã intercropping and Piatã pasture for two years in rotationS/MP/P/P
T10Native savanna (reference for soil parameters)NS
Table 3. Content of total nitrogen (TN), total organic carbon (TOC), particulate organic carbon (POC) and mineral-associated organic carbon (MOC) in layers 0.0–0.10, 0.10–0.20 and 0.20–0.40 m.
Table 3. Content of total nitrogen (TN), total organic carbon (TOC), particulate organic carbon (POC) and mineral-associated organic carbon (MOC) in layers 0.0–0.10, 0.10–0.20 and 0.20–0.40 m.
Production SystemTNTOCPOCMOC
g kg−1
0.00–0.10 m
M 12.12 ab33.07 bcd9.58 cd23.48 a
MP1.98 b31.90 cde9.98 cd21.91 a
S1.78 b27.60 e7.18 e20.41 a
P1.75 b31.40 cde10.42 bcd20.97 a
MP/S1.97 b35.47 bc10.71 bcd24.75 a
MP/P2.02 b36.50 b12.19 b24.30 a
MP/P/P2.16 ab31.93 cde11.28 bc20.65 a
MP/P/S1.95 b29.07 de8.09 de20.16 a
S/MP/P/P1.89 b34.70 bc10.22 bcd24.48 a
NS2.53 a42.90 a17.09 a25.81 a
P Fcalc0.0506<0.0001<0.00010.1011
0.10–0.20 m
M1.56 c26.07 d4.70 g21.36 bc
MP2.00 a34.63 ab9.95 a24.68 abc
S1.59 c26.77 d5.16 fg21.60 bcd
P1.54 c29.83 bcd8.74 b21.09 cd
MP/S1.64 c31.90 abc5.90 ef26.00 ab
MP/P1.64 c30.50 bcd7.02 cde23.47 bcd
MP/P/P1.55 c27.77 cd6.14 def21.62 bcd
MP/P/S1.76 bc29.03 cd7.19 cd21.84 bcd
S/MP/P/P1.45 c25.47 d5.45 fg20.01 d
NS2.23 a35.93 a7.72 bc28.21 a
P Fcalc0.00410.0046<0.00010.0255
0.20–0.40 m
M1.45 a24.90 a4.38 bc20.52 a
MP1.58 a27.13 a5.08 b22.04 a
S1.40 a22.73 a4.02 bc18.70 a
P1.32 a25.13 a6.55 a18.57 a
MP/S1.47 a27.53 a6.86 a20.67 a
MP/P1.38 a25.90 a5.09 b20.80 a
MP/P/P1.50 a25.73 a6.43 a19.30 a
MP/P/S1.50 a27.20 a6.98 a20.22 a
S/MP/P/P1.34 a24.67 a5.00 b19.67 a
NS1.52 a23.30 a3.44 c19.86 a
P Fcalc0.27250.1317<0.00010.4648
1 Means followed by different letters differ from each other by the t test (LSD, p ≤ 0.05). M: Single Maize in annual succession (monoculture); MP: Maize intercropped with U. brizantha cv. Piatã in annual succession; S: Soybean in annual succession (monoculture); P: U. brizantha cv. Piatã pasture implanted in single cultivation; MP/S: Maize + Piatã intercropped and Soybean in annual rotation; MP/P: Sequence of Maize + Piatã intercropping and Piatã pasture in annual rotation; MP/P/P: Sequence of Maize + Piatã intercropping and Piatã pasture for two years; MP/P/S: Sequence of Maize + Piatã intercropping, Piatã pasture and Soybean in annual rotation; S/MP/P/P: Sequence of Soybean, Maize + Piatã intercropping and Piatã pasture for two years in rotation; NS: Native savanna.
Table 4. Soil C/N Ratio, Carbon Pool Index (CPI), Lability (L), Lability Index (LI) and Carbon Management Index (CMI) in layers 0.0–0.10, 0.10–0.20 and 0.20–0.40 m.
Table 4. Soil C/N Ratio, Carbon Pool Index (CPI), Lability (L), Lability Index (LI) and Carbon Management Index (CMI) in layers 0.0–0.10, 0.10–0.20 and 0.20–0.40 m.
Production SystemC/NCPILLICMI
0.00–0.10 m
M 115.70 a0.77b cd0.41 bc0.64 bcd49 cd
MP16.22 a0.74 cd0.45 bc0.70 bcd52 cd
S15.54 a0.64 e0.35 c0.54 d34 e
P17.93 a0.73 cde0.51 b0.80 bc57 bcd
MP/S17.99 a0.83 bc0.43 bc0.68 bcd56 bcd
MP/P18.08 a0.85 b0.50 b0.78 bc66 b
MP/P/P15.13 a0.74 cde0.54 ab0.83 ab62 bc
MP/P/S15.13 a0.67 de0.44 bc0.67 bcd46 de
S/MP/P/P18.61 a0.81 bc0.42 bc0.63 cd51 cd
NS16.95 a1.00 a0.67 a1.00 a100 a
P Fcalc0.0846<0.00010.00970.0076<0.0001
0.10–0.20 m
M16.74 bc0.72 d0.22 d0.80 e58 f
MP17.58 abc0.97 ab0.40 a1.48 a142 a
S16.83 bc0.75 cd0.24 cd0.88 cde65 f
P19.29 a0.83 bcd0.42 a1.53 a125 b
MP/S19.46 a0.89 abc0.22 cd0.82 de73 ef
MP/P18.55 ab0.85 bcd0.30 b1.09 bc92 cd
MP/P/P17.97 abc0.77 cd0.28 bc1.04 bcd81 de
MP/P/S16.62 bc0.81 cd0.33 b1.23 b98 c
S/MP/P/P17.55 abc0.71 d0.27 bcd1.00 cde71 ef
NS16.17 c1.00 a0.27 bcd1.00 cde100
P Fcalc0.00090.0046<0.0001<0.0001<0.0001
0.20–0.40 m
M17.20 a1.07 a0.21 bc1.25 bc135 bc
MP17.25 a1.16 a0.23 bc1.34 bc155 b
S16.34 a0.98 a0.22 bc1.26 bc121 bc
P19.00 a1.08 a0.35 a2.06 a222 a
MP/S18.69 a1.17 a0.33 a1.94 a230 a
MP/P18.76 a1.11 a0.24 b1.42 b157 b
MP/P/P17.37 a1.11 a0.33 a1.94 a212 a
MP/P/S18.19 a1.17 a0.34 a2.01 a235 a
S/MP/P/P18.63 a1.06 a0.25 b1.46 b155 b
NS15.40 a1.00 a0.17 c1.00 c100 c
P Fcalc0.11170.1453<0.0001<0.0001<0.0001
1 Means followed by different letters differ from each other by the t test (LSD, p ≤ 0.05). M: Single Maize in annual succession (monoculture); MP: Maize intercropped with U. brizantha cv. Piatã in annual succession; S: Soybean in annual succession (monoculture); P: U. brizantha cv. Piatã pasture implanted in single cultivation; MP/S: Maize + Piatã intercropped and Soybean in annual rotation; MP/P: Sequence of Maize + Piatã intercropping and Piatã pasture in annual rotation; MP/P/P: Sequence of Maize + Piatã intercropping and Piatã pasture for two years; MP/P/S: Sequence of Maize + Piatã intercropping, Piatã pasture and Soybean in annual rotation; S/MP/P/P: Sequence of Soybean, Maize + Piatã intercropping and Piatã pasture for two years in rotation; NS: Native savanna.
Table 5. Organic carbon content in the chemical fractions of soil organic matter: Humin (CHum), humic acid (CHA), fulvic acid (CFA) and CHA/ CFA ratio in layers 0.0–0.10, 0.10–0.20 and 0.20–0.40 m.
Table 5. Organic carbon content in the chemical fractions of soil organic matter: Humin (CHum), humic acid (CHA), fulvic acid (CFA) and CHA/ CFA ratio in layers 0.0–0.10, 0.10–0.20 and 0.20–0.40 m.
Production SystemCHumCHACFACHA/CFA
g kg−1
0.0–0.10 m
M 110.60 c5.87 abcd7.85 a0.74 cd
MP12.97 bc4.63 de6.70 abc0.70 d
S12.52 bc4.61 de3.75 f1.24 a
P12.47 bc3.86 e6.08 bcde0.64 d
MP/S11.82 bc5.00 bcde8.04 a0.62 d
MP/P13.89 b6.22 ab4.75 ef1.31 a
MP/P/P13.23 bc6.03 abc6.20 bcde1.00 b
MP/P/S13.21 bc4.75 cde5.15 def0.94 bc
S/MP/P/P12.87 bc5.04 bcde5.37 cde0.94 bc
NS18.12 a7.15 a7.04 ab1.01 b
P Fcalc0.00350.0028<0.0001<0.0001
0.10–0.20 m
M10.60 c5.45 de4.17 cde1.35 ab
MP13.33 ab6.92 bc4.44 cde1.57 a
S10.47 c5.78 cd3.44 e1.86 a
P13.33 ab5.01 de5.46 bc0.92 bc
MP/S11.06 bc8.37 a5.03 bcde1.73 a
MP/P11.22 bc7.20 ab3.80 de1.90 a
MP/P/P10.27 c5.17 de3.67 de1.42 ab
MP/P/S11.17 bc4.30 e7.36 a0.60 c
S/MP/P/P11.52 bc5.54 de5.15 bcd0.88 bc
NS14.67 a5.65 d6.48 ab0.88 bc
P Fcalc0.0163<0.00010.00140.0006
0.20–0.40 m
M8.21 de4.82 a4.55 bc1.07 abc
MP12.59 ab4.30 a5.27 ab0.82 c
S8.14 e4.30 a3.65 c1.21 ab
P10.56 bc4.06 a5.68 ab0.72 c
MP/S11.83 ab4.67 a5.41 ab0.86 bc
MP/P10.43 bcd4.56 a5.87 a0.78 c
MP/P/P9.50 cde3.36 a4.52 bc0.84 c
MP/P/S13.01 a4.45 a4.47 bc1.00 abc
S/MP/P/P10.97 abc4.07 a5.41 ab0.75 c
NS11.88 ab4.593.38 c1.37 a
P Fcalc0.00190.14950.00410.0218
1 Means followed by different letters differ from each other by the t test (LSD, p ≤ 0.05). M: Single Maize in annual succession (monoculture); MP: Maize intercropped with U. brizantha cv. Piatã in annual succession; S: Soybean in annual succession (monoculture); P: U. brizantha cv. Piatã pasture implanted in single cultivation; MP/S: Maize + Piatã intercropped and Soybean in annual rotation; MP/P: Sequence of Maize + Piatã intercropping and Piatã pasture in annual rotation; MP/P/P: Sequence of Maize + Piatã intercropping and Piatã pasture for two years; MP/P/S: Sequence of Maize + Piatã intercropping, Piatã pasture and Soybean in annual rotation; S/MP/P/P: Sequence of Soybean, Maize + Piatã intercropping and Piatã pasture for two years in rotation; NS: Native savanna.
Table 6. Percentage of water-stable aggregates larger than 2.0 mm (%Aggregation > 2 mm), geometric mean diameter (GMD), mean weight-diameter (MWD) and the aggregate stability index (ASI).
Table 6. Percentage of water-stable aggregates larger than 2.0 mm (%Aggregation > 2 mm), geometric mean diameter (GMD), mean weight-diameter (MWD) and the aggregate stability index (ASI).
Production System%Aggregation > 2 mmGMDMWDASI
%Mm%
M 178.6 bcd3.13 bcde4.08 bcd96.13 a
MP82.1 bcd3.40 bcd4.24 bcd97.22 a
S57.4 e2.21 e3.21 e94.61 a
P90.9 ab4.08 ab4.60 ab98.06 a
MP/S68.0 de2.54 de3.66 de95.27 a
MP/P85.1 ab3.62 bc4.37 ab97.04 a
MP/P/P84.4 abc3.53 bcd4.34 abc96.70 a
MP/P/S77.6 bcd3.14 bcde4.06 bcd96.36 a
S/MP/P/P69.4 cde2.69 cde3.70 cde94.49 a
NS98.1 a4.76 a4.91 a99.38 a
P Fcalc0.00150.00210.00170.1300
1 Means followed by different letters differ from each other by the t test (LSD, p ≤ 0.05). M: Single Maize in annual succession (monoculture); MP: Maize intercropped with U. brizantha cv. Piatã in annual succession; S: Soybean in annual succession (monoculture); P: U. brizantha cv. Piatã pasture implanted in single cultivation; MP/S: Maize + Piatã intercropped and Soybean in annual rotation; MP/P: Sequence of Maize + Piatã intercropping and Piatã pasture in annual rotation; MP/P/P: Sequence of Maize + Piatã intercropping and Piatã pasture for two years; MP/P/S: Sequence of Maize + Piatã intercropping, Piatã pasture and Soybean in annual rotation; S/MP/P/P: Sequence of Soybean, Maize + Piatã intercropping and Piatã pasture for two years in rotation; NS: Native savanna.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Silva, J.F.d.; Gontijo Neto, M.M.; Silva, G.F.d.; Borghi, E.; Calonego, J.C. Soil Organic Matter and Aggregate Stability in Soybean, Maize and Urochloa Production Systems in a Very Clayey Soil of the Brazilian Savanna. Agronomy 2022, 12, 1652. https://doi.org/10.3390/agronomy12071652

AMA Style

Silva JFd, Gontijo Neto MM, Silva GFd, Borghi E, Calonego JC. Soil Organic Matter and Aggregate Stability in Soybean, Maize and Urochloa Production Systems in a Very Clayey Soil of the Brazilian Savanna. Agronomy. 2022; 12(7):1652. https://doi.org/10.3390/agronomy12071652

Chicago/Turabian Style

Silva, Juslei Figueiredo da, Miguel Marques Gontijo Neto, Gustavo Ferreira da Silva, Emerson Borghi, and Juliano Carlos Calonego. 2022. "Soil Organic Matter and Aggregate Stability in Soybean, Maize and Urochloa Production Systems in a Very Clayey Soil of the Brazilian Savanna" Agronomy 12, no. 7: 1652. https://doi.org/10.3390/agronomy12071652

APA Style

Silva, J. F. d., Gontijo Neto, M. M., Silva, G. F. d., Borghi, E., & Calonego, J. C. (2022). Soil Organic Matter and Aggregate Stability in Soybean, Maize and Urochloa Production Systems in a Very Clayey Soil of the Brazilian Savanna. Agronomy, 12(7), 1652. https://doi.org/10.3390/agronomy12071652

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop