Organic matter compartments in an Ultisol under integrated agricultural and livestock production systems in the Cerrado. Organic matter compartments in an Ultisol under integrated agricultural and livestock production systems in the Cerrado

: I ntegrated agricultural production systems have the potential to increase organic matter content, which is reflected in the soil carbon (C) and nitrogen (N) concentrations. Here, we evaluated the C and N stocks and its compartments in a typical distro cohesive yellow Ultisol under the no-till (NT) and crop-livestock integration (CLI) systems, in eastern Maranhão. Five areas with different management strategies were evaluated, more specifically, one area was managed under the NT system in succession for 14 years (soybean/millet), three areas had different CLI system adoption histories (i.e., CLI was adopted 2, 4, or 8 years prior to sample collection); and finally, one area consisted of native Cerrado (savannah) vegetation. Soil samples were collected at depths of 0–0.10, 0.10–0.20, 0.20–0.30, and 0.30–0.50 m,to analyze the content and total stocks of carbon (C) and its compartments (physical and chemical), N content, soil microbial biomass, and basal respiration. Results revealed higher content and stock of C and N in areas with CLI. The principal humic substance reported was humin, and higher concentrations of C associated with minerals were present in areas under more recent CLI adoption (2 and 4 years). We demonstrated that the adoption of CLI, even if recent, leads to immediate increases in the concentrations of C and its fractions as a result of using forage crops. NT for 14 years and CLI for 8 years exhibited higher levels of C management and higher soil biological activity due to the greater stability of these systems.

benefiting both agricultural production and livestock; in addition, the resultant increase in dry matter production per unit area can be used as soil cover for cultivation in the no-till (NT) system, or even as an alternative in the recovery of degraded pastures (CHIODEROLI et al., 2012;MAZZUCHELLI et al., 2020).
Due to its capacity to provide economic and environmental benefits, the NT system is widely considered as a conservation-orientated system. This perception, together with the ensuing increase in crop yields, has made the NT system popular in Brazil, the management practice is being studied in various parts of the country and in the 2017/18 crop year, it was implemented in 32.8 million hectares of land (FERREIRA et al., 2020).
IAPS can be modified to use the most suitable intercropping for the edaphoclimatic conditions inherent to different regions of Brazil. A main advantage of adopting these systems, which are widely accepted as contributing to the generation of more sustainable and efficient practices, is the improvement of soil quality (SKORUPA & MANZATTO, 2019). Compared to other conservation management models that do not make use of pastures, the CLI system fosters a greater increase in carbon (C) content in the soil due to its maintenance principles, i.e., permanent soil cover, minimum disturbance, intercropping and crop rotation, and the introduction of animals into the system (GUESMI et al., 2019). In fact, in the Brazilian Cerrado (savannah) reported in the state of Goiás, increases of up to 40% in soil C content have been found under CLI systems compared to areas under the NT system (GAZOLLA et al., 2015).In addition, studies in the Cerrado by SILvA et al. (2016), SOARES et al. (2019), andSOUSA et al. (2020) also report improvement in the content and compartments of soil C through the adoption of CLI.
However, several factors have compromised the successful development of these management systems, including inherent risks related to the complexity of IAPS management, regional intricacies, limited information, and the unavailability of suitable technologies (TAKAHASHI et al., 2019).
Sandy soils, reported in agricultural frontier regions, the major constraints related to management information and organic matter gains when adopting IAPS, mainly by improving soil attributes (DONAGEMMA et al., 2016).
Despite the sandy properties of some soils in agricultural frontier regions, the presence of the cohesive horizon can diminish the potential of agricultural systems designed for soils in the Cerrado that have developed in sediments of the Formação Barreiras. This is as a result of cohesion and the ensuing high soil density and resistance to penetration, which forces dramatic restrictions to the deepening of root systems, leading to are duction in the water and nutrient storage capacity of plants, development (VIEIRA et al., 2012;RIBEIRO et al., 2016), and ultimately, soil C stocks (GATTO et al., 2010).
In eastern Maranhão (MA), the position of the landscape determines the distinction of cohesive soils, for example, those that develop in concave pedoformex hibit the greatest expression of the cohesive character. Soil surveys carried out in the region indicated that the B textural horizon is below 60 cm depth (DANTAS et al., 2014;RESENDE et al., 2014). Thus, the adoption of systems that favor root production, such as NT and IAPS should be considered.
In this study, we hypothesized that the adoption of IAPS, such as the CLI, increases the soil C and nitrogen (N) stock compared to the NT systems. Thus, this study evaluated and compared the composition and accumulation of total C and N stocks and their compartments under the NT and CLI systems in a typical dystrocohesive yellow Ultisol in eastern Maranhão.

Description of the area and management systems
The study was conducted on the Barbosa farm located in the municipality of Brejo, MA, Brazil (03°42ʹ0.93ʺS and 42°56ʹ25.57ʺW) (Figure 1).The Barbosa farm has an average altitude of 95 m, a flat to gently undulating relief, and a hot and humid tropical climate. The average annual temperature of the study area is 27 °C and rainfall, which typically occurs between the months of November and April, has an annual average of 1,835 mm ( Figure 1). The area falls under the Cerrado type biome, and the soil is classified as typical dystrocohesive yellow Ultisol (DANTAS et al., 2014).
Three areas were selected which were managed using NT systems that followed soybean/ millet succession and had different histories of CLI adoption; one area was managed exclusively under NT systems; and one area consisted of native Cerrado vegetation ( Figure 2). The areas are described as follows: (i) The NT system for 14 years (NT14); this area was first initiated (originally native Cerrado biome) in 2003 with plowing and harrowing. Thereafter, 2 t ha -1 of limestone was incorporated, and the first year of cultivation (soybean) was in 2004. The subsequent years followed the NT system (first crop, soybean; second crop, millet) until 2018, and fertilization was under taken based on the recommendation for soybean culture (SOUSA & LOBATO, 2004).
(ii) The CLI system was adopted 8 years ago (CLI8); this area, which has the same management history as NT14 until 2009; in 2010, the soil was turned over (via plowing and harrowing) and 3.7 t ha -1 of calcareous limestone was applied. Thereafter, the area was managed in a CLI system (maize + Urochloa ruziziensis with an input of 0,7 animal units [male/ female cattle -AU = 450 kg] per hectare in the offseason) until 2018, when the area returned to the NT system (soybean and millet) and fertilization was undertaken based on the recommendation for soybean and corn crops (SOUSA & LOBATO, 2004). (iii) CLI adopted 4 years ago (CLI4); this area, which has the same management history as NT14 until 2009, had its soil turned over (via plowing and harrowing) and 3.7 t ha -1 of limestone incorporated in 2010. In 2014, the area was managed using the CLI system (maze + Urochloa ruziziensis with an input of 0.8 animal units [male cattle -AU = 450 kg] per hectare in the off-season) until 2018, when the area returned to the NT system (soybean and millet) and fertilization was undertaken based on the recommendation for soybean and maze culture (SOUSA & LOBATO, 2004). (iv) CLI adopted two years ago but with two animal introductions (CLI2); this area, has the same management history as NT14 until 2011. Thereafter, in 2012, the area had its soil turned over (via plowing and harrowing), 3.8 t ha -1 of calcitic limestone incorporated, and CLI adopted (maize + Urochloa ruziziensis with an input of 1.5 animal units [male cattle -AU = 450 kg] per hectare in the off-season). From 2013 to 2015, the area was managed under the NT system (soybean and millet), and in 2016, after being subsoiled to 0.30 m, the area returned to being managed under the CLI system (corn + Urochloa ruziziensis with an input of 0.8 animal units [male cattle -AU = 450 kg] per hectare in the offseason). In the following years, the area returned to the NT system (soybeans and millet) and fertilization was undertaken based on the recommendation for soybean and maize culture (SOUSA&LOBATO, 2004). (v) Native Cerrado (FOREST); this area of native Cerrado biome vegetation was included as a reference point.

Soil Sampling and Analysis
In each area, four trenches with a width, length, and depth of approximately 0.60 × 0.60 × 0.60 m, respectively, were dug in randomly distributed positions, 50 m apart. Four not deformed soil samples were taken from each trench with the aid of aluminum volumetric rings (98.17 cm 3 volume, 0.05m diameter, 0.05 m height). An additional nine deformed samples were collected using a Dutch auger; one was taken from the central point of the respective trench, and eight were taken from points around the trenches spaced 3 m apart. Soil samples were collected in June 2018 (the dry season), at 0-0.10, 0.10-0.20, 0.20-0.30, and 0.30-0.50 m soil depth. Ciência Rural, v.52, n.10, 2022. Barbosa et al. Part of the not deformed samples were used in determining the soil density (SD) via the volumetric ring method, and another part of the samples were used for chemical characterization and particle size analyses using the pipette method (TEIXEIRA et al., 2017) (Table 1).
Soil biomass carbon (SBM) was only analyzed in the superficial layer (0-0.20 m) using the irradiation extraction method (ISLAM & WEIL, 1998), and after the samples were collected, refrigerated, and taken to the laboratory, the basal respiration (C-CO 2 ) was obtained by incubating the samples for 7 days (ALEF et al., 1995). The microbial quotient (qMIC) obtained using the CBM (soil microbial biomass carbon)/TOC ratio (ANDERSON & DOMSCH, 1989) and the metabolic quotient (qCO 2 ) from the C-CO 2 /CBM ratio (ANDERSON & DOMSCH, 1990) were also calculated.
After obtaining the laboratory data, the C stock (CS) and N stock (NS) were determined using the layer method (BAYER et al., 2000) and equivalent soil mass (ELLERT & BETTANY, 1995), this was carried out using the calculations by layer and thereafter using the sum of stocks of all the layers, resulting in an accumulated stock (0-0.50 m). The carbon management index (CMI) and its components were determined (BLAIR; LEFROY & LISLE, 1995) while considering the POC as representative of the labile fraction of the TOC, the MAOC as the non-labile fraction, and the carbon lability (CL) as the ratio between POC/MAOC (DIEKOW et al., 2005;VIEIRA et al., 2007).

Data analysis
The limitations for extrapolation produced by pseudo replication were considered, and an analysis of variance (F test) was performed. Considering the different management strategies as treatments in an entirely randomized design and, depending on normality, a test of means (Tukey, 5% probability) was performed using the statistical program SISvAR 5.6 (FERREIRA, 2014). Correlation and principal component analysis (PCA) were performed using the statistical program R (R CORE TEAM, 2017) to assess the response variables as a function of depth. variables were standardized for the multivariate tests (MANLY, 2008). Correlation between the variables were interpreted according to the following

RESULTS
The highest TOC values in the 0-0.20 m layer were observed in the CLI systems (CLI2, CLI4 and CLI8), while the lowest values were in NT14 and FOREST (Table 2). In the subsurface layers, the highest concentrations of TOC were observed in CLI2, an area with recent CLI adoption; these values differed from the values observed in NT14 and FOREST (Table  2). In general, the TOC content observed in CLI2 was 15% higher than in NT14 ( Table 2).
The highest concentrations of TN were observed in area CLI4, while they did not differ from those of CLI8 in the 0-0.10, 0.20-0.30,and 0.30-0.50 m layers, nor from those of NT14 in the 0.20-0.30 m layer ( Table 2). The C/N ratio; however, showed an inverse signal to that of TN and a similar signal to that observed in TOC, with area CLI2 having the highest ratio regardless of soil layer ( Table 2).
The highest CS values were observed in the areas CLI2 and CLI4, and these values differed from those observed in areas NT14 and FOREST     until the 0-0.30 m layer (Table 2). Below this layer, the highest CS was observed in CLI4, being 9.5 and 10.2 Mg ha -1 higher than those observed in NT14 and Forest, respectively ( Table 2). The highest NS values were observed in area CLI4 for the 0-0.10, 0.20-0.30, and 0.30-0.50 m layers. These values; however, did not differ from CLI8 in the 0-0.10, 0.20-0.30, and 0.30-0.50 m layer, nor from FOREST in the 0-0.10 and 0.30-0.50 m layers (Table 2). Finally, in the 0.10-0.20 m layer, the highest NS was observed in CLI2 and CLI8 (Table 2).
Areas CLI4 and CLI2 displayed a higher CS accumulation than areas CLI8, NT14, and FOREST, the latter being 22.9% lower compared to CLI4 (Table 2). NS accumulation for the 0-0.50 m layer was highest in area CLI4 (Table 2).
high HA/FA ratios did not differ from all layers in a FOREST area, nor from the superficial layers (0-0.10 and 0.10-0.20 m) in area CLI8 (Table 3). Among the humic fractions, humin was reported in the highest concentrations. The humin concentrations observed in the 0-0.10, 0.10-0.20,and 0.20-0.30 m layers in area CLI2 did not differ from the same layers in area CLI4, nor from the superficial layer (0-0.10 m) in area CLI8, nor from the 0.10-0.20 m layer in area NT14 (Table 3). In the 0.30-0.50 m layer, the highest concentrations of humin were reported in area CLI4, and these concentrations did not differ from areas CLI2 and CLI8 (Table 3).
With regard to the physical fractionation of C, no differences in POC by area or depth were observed, however, changes in the MAOC were observed (Table 3). The areas with more recent adoption of CLI (i.e., CLI2 and CLI4) showed higher values of MAOC in the 0-0.10 and 0.10-0.20 m layers compared to areas NT14 and FOREST (Table 3). In the subsurface layers (0.20-0.30 and 0.30-0.50 m), area CLI4 had higher MAOC values than area NT14, and higher values in the 0.20-0.30 m layer compared to area FOREST (Table 3).
The highest CL values in the superficial layer (0-0.10 m) were observed in area NT14; in the 0.10-0.20 and 0.30-0.50 m layers, the highest values were observed in areas NT14 and CLI8; and in the 0.20-0.30 m layer, the highest values were observed in areas NT14 and FOREST (Table 3).     Considering the POC, MAOC, and CL data used when determining the CMI, it was found that for the 0-0.20 m layer, all areas under agricultural management had CMI values higher than the reference value (100) of the FOREST area, with the highest CMI values being observed in area NT14. However, in the subsurface layers (0.20-0.30 and 0.30-0.50 m), higher CMI values were observed in areas CLI8 and NT14, with values being 11% and 46% higher than the CMI of the reference area, respectively ( Table 3).
The CBM and qMIC values reflected differences according to the management type adopted in each area. The highest values were observed in areas NT14 and CLI8 and the lowest were observed in the FOREST area (Table 4). The highest concentration for basal respiration (C-CO 2 ) was observed in area NT14, while the highest metabolic quotient (qCO 2 ) value was observed in area CLI2.
The  Figure 3D). Joint analyses of the results (the soil biology data were excluded, as they were sampled in the superficial layer i.e., the 0-0.20 m layer), revealed relationships between the soil variables. More specifically, two principal components (CP1 and CP2) explained 84.40%, 88.50%, 82.79%, and 82.53% of the total data variance for the 0-0.10, 0.10-0.20, 0.20-0.30,and 0.30-0.50 m layers, respectively. Variables with weight coefficients > 0.3 were considered relevant in the PCA (Table 5).
All layers in area CLI4 were associated with positive CP1values, with the TN and NS variables displaying particularly high scores. However, area NT14 was associated with negative values for CP1, with the most notable variables being the C/N ratio for all layers and POC for the 0-0.10, 0.10-0.20, and 0.20-0.30 m layers ( Figures 4A, B, C and D).
Nonetheless, apart from the 0.20-0.30 m layer, area FOREST displayed some similar results; the FOREST area was associated with positive CP2 values, while conversely, area CLI2 was associated with negative CP2 values, largely as a result of a negative humin concentration ( Figures 4A, B, C and D).
In the surface layer, no attribute was associated with positive values for CP2, which was defined by the FOREST area. Thus, at the depth of 0-0.10 m, the highest values of N, C, its fractioning (chemical and physical), and stocks were all related to the agricultural management strategies ( Figure 4A). CLI2 = no-till soy/no-till twice with history of crop-livestock integration, CLI4 = no-till soy/no-till with history of crop-livestock integration system once for 4 years prior to soil collection, CLI8 = no-till soy/no-till with history of crop-livestock integration system once for 8 years prior to soil collection, NT14 = no-till soy/no-till for 14 years and no history of animal inputs, FOREST = area under native Cerrado vegetation, CBM = soil microbial biomass carbon, SBM = soil microbial biomass nitrogen, qMIC = microbial quotient, Respiration = basal respiration, qCO2 = metabolic quotient. 1 Means followed by the same letter in the column do not have differences between them based on the Tukey's test (5%).

DISCUSSION
The higher contents and stocks of C and N, and C/N ratio observed in areas under the CLI system compared to areas under exclusive soil management in NT (Table 2), is likely attributed to the intercropping of corn with Urochloa, which can lead to a greater amount of cultural residues in the soil (MATEUS et al., 2020). This affect is particularly prevalent in Urochloa pasture, a crop which causes a considerable addition of organic matter to the soil, ultimately increasing C and N contents (SIGNOR et al., 2018). The productivity of forage dry mass when Urochloa is intercropped with corn, depending on the cut, in fall/winter can enter the system at approximately 1.6-1.7 t/ha in the first cut and 3.9-4.2 t/ha in the second cut (MATEUS et al., 2020). Furthermore, the higher values of accumulated C in the CLI system, independent of the time of adoption, can be justified by the increment of C in depth. The land use change provided gains or stability for soil C at up to 0.5 m depth, a finding which can be justified by the co-dependence verified in the correlation results.
The correlations presented here allow several other inferences to be made. Firstly, there was a high correlation between TOC and MAOC in all layers, which also provided high correlations with the stock of C, this suggested that the stabilization of soil organic matter (SOM) in the long term occurs by increasing the levels of MAOC (LIMA et al., 2016). This result indicates that in the edaphoclimatic conditions of eastern MA, the accumulation of TOC depends on the increase in MAOC content. Therefore, it is possible to infer that the stabilization of SOM is fundamental for an increase in TOC content. The HA/FA ratio observed in the 0-0.20 m layer among all the studied areas was higher than the average value of 1.2 for the Cerrado biome (SANTOS et al., 2013). This finding is important, because an increase in the HA/FA ratio allows the chemical quality of the soil to improve, as a higher HA/FA ratio contributes to a greater cation exchange capacity, it allows the formation of complexes with various metal ions and acts as a buffer for chemical reactions (ROSA et al., 2017). In addition, the higher humin content found, especially in areas under the CLI, corroborate with what has been observed in tropical soils (ZHU et al., 2014). Furthermore, the higher humin and MAOC values observed in areas CLI2 and CLI4 (Table 3) likely result from the more recent intercropping of maize with pasture in area CLI8, these are plants that release large amounts of organic material via rhizo deposition (CARMO et al., 2012).
The similar CBM and qMIC values reported among the native vegetation and areas with different management strategies, indicates the efficiency of both the CLI system and the strategy applied in area NT14, in the immobilization of the available C in microbial cells in the soil (DADALTO et al., 2015). However, in area CLI2, the higher qCO 2 value is likely due to the quality of the organic matter, which has a high C/N ratio (Table 2) and lower C lability (Table 4). This may have generated microbial stress during its degradation, which along with the fact that the recent soil mobilization promotes exposure of organic matter to oxidation and mineralization, may have ultimately allowed a greater influx of energy (C) through the CBM (LACERDA et al., 2013).
TOC correlated best with the other variables and considering that the PCA showed a distinction between the systems studied, it suggests that the CLI system is differentiated in relation to NT14 and FOREST (Figure 1).
The recent introduction of pastures and animals in the off-season of soybean cultivation in areas CLI2 and CLI4 were likely the main factors that distinguished these areas from area NT14. Pastures have high dry biomass production capacity and root mass density (SARTO et al., 2020). Thus, the use of Urochloa in the CLI system promotes better soil Attributes chemical quality, which contributes to the C stock and fractionation, and the soil CMI (SOARES et al., 2019); in addition, the added material contributed by animal excrement also improved soil quality by raising C and N contents (GUESMI et al., 2019). In addition to the above, it is worth noting that the distance observed between the NT14 and CLI areas is due to the type and quality of the C present in each environment, while the CLI and the NT14 area was more associated with particulate organic C (Figure 3), that is, more easily degraded, it was also the area that showed higher C lability (Table 3).
Consideration of the effect of the management strategy adopted and the C reported in the soil over time reveals an interesting finding, 14 years after the removal of native vegetation for agricultural ex-ploration, area NT14 showed a 4.7 t/ha increase in C, which translates to an average gain of 0.33 t/ha/year; area CLI8, however, showed a 14.5 t/ha increase in C, which represents a gain of 1.03 t/ha/year. In other words, the exclusive use of soybean with millet under the NT strategy does not increase the C in the system in the same quantity as in the CLI system, which uses crop rotation, intercropping of corn with forage, and the input of animals, is adopted. Although, the 112.3 t/ha increase in C content observed in area CLI2 was higher than the increases observed in areas NT14 and CLI8, it was still 8.7 t/ha lower than the increase in C content observed in area CLI4 (121.0 t/ha), this is despite the sub soiling that occurred, with aggregate breakage exposing the physically protected C to microbial oxidation (MELO et al., 2016), which corroborates with the higher rate of qCO 2 (Table 4). These results reveal that the increase in soil C is easily influenced by management strategies, even in conservation systems.
Maize intercropped with Urochloa contribute residues to the surrounding system that have a high C/N ratio and high lignin quantity compared to the residues derived from soybean/ millet crops, which allows the straw greater perenniality (GARCIA et al., 2014;SOUSA et al., 2019). In this sense, when observing the trends of the results presented here, which were processed through uni and multivariate statistical analyses showed that the area without recent history of soil tillage which had the CLI system introduced 4 years ago, was most efficient in increasing soil C (121 t/ha) and N (3.4 t/ha) ( Table 2), since it is linked to the largest number of studied C and N response variables (Figure 4).

CONCLUSION
CLI adopted in typical distro cohesive yellow ultisol in the Cerrado of eastern Maranhão incorporates higher soil C and N contents and stocks compared to the NT system.
Areas under the CLI system with adoption times of 2, 4, and 8 years ago, increase the C stock in soil up to 0.5 m in depth by 14.5%, 23.4%, and 10.0%, respectively, in relation to the NT system with soybean and millet.
The NT system with soybean and millet has lower N stocks in soil up to 0.5 m in depth, in relation to the native Cerrado and areas managed by the CLI system. Humin is the predominant soil fraction in the studied areas, with higher concentrations being associated with areas managed by the CLI system. Higher concentrations of C fraction associated with minerals are related to the recent adoption of the CLI system. The NT and CLI conservationist agricultural management systems presented CMI percentages above 100% for soil up to 0.20 m in depth. In conservation management, tillage history alters C biomass and microbiological activity.