Next Article in Journal
Do Farming Scale and Mechanization Affect Moral Hazard in Rice Harvest Outsourcing Service in China?
Previous Article in Journal
Design and Evaluation of a Mechanical Floor-Standing Jujube Picker
Previous Article in Special Issue
Conservation Agricultural Practices Impact on Soil Organic Carbon, Soil Aggregation and Greenhouse Gas Emission in a Vertisol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

No-Till Mitigates SOC Losses after Grassland Renovation and Conversion to Silage Maize

1
Institute of Crop Sciences and Plant Breeding, Grass and Forage Sciences/Organic Agriculture, Christian-Albrechts-University, Hermann-Rodewald-Straße9, D-24118 Kiel, Germany
2
Grass-Based Dairy Systems, Animal Production Systems Group, Wageningen University (WUR), 6700 HB Wageningen, The Netherlands
*
Author to whom correspondence should be addressed.
Agriculture 2022, 12(8), 1204; https://doi.org/10.3390/agriculture12081204
Submission received: 2 July 2022 / Revised: 4 August 2022 / Accepted: 8 August 2022 / Published: 12 August 2022
(This article belongs to the Special Issue Tillage Effects on Soil Carbon Sequestration and Soil Organic Matter)

Abstract

:
Many studies recommend no-till (NT) to increase soil organic carbon (SOC) in the topsoil (<30 cm) of arable land to counterbalance greenhouse gas emissions. Its potential use to mitigate SOC losses during conversion and renovation of grassland ecosystems in the top meter soil is yet to be determined. The SOC dynamics of a 10-year-old grassland converted to silage maize (CM) and renovated and seeded (GR) using either conventional tillage (CT) or NT were compared to an undisturbed grassland control (GC) for 7 years, across three fixed soil depth increments (0–30, 30–60, 60–90 cm). The annual C inputs (Cinput) from crop residues were further analyzed. The systems were either non-fertilized (N0) or fertilized with mineral N (N1) according to a demand of 180 and 380 kg N ha−1 year−1 in the silage maize and grassland systems, respectively. For the 7-year period, the renovated grassland using NT ensured maintenance of the initial SOC in the topsoil, while a conversion toward arable cropping resulted in SOC losses, regardless of the tillage method. The use of NT during conversion significantly reduced these losses from 2.5 Mg ha−1 year−1 to 1.8 Mg ha−1 year−1, for a 28% reduction compared to CT. In the subsoil (30–90 cm), SOC remained stable and was not affected by the cropping systems nor by the tillage method. Reduced annual Cinput was found as the main factor affecting SOC losses after grassland removal, regardless of the tillage method. Our findings highlight the potential of NT to mitigate annual SOC losses after grassland conversion if annual Cinput remains high.

1. Introduction

The conversion and renovation of grasslands using conventional tillage (CT) often result in high soil organic carbon (SOC) losses in agricultural ecosystems [1]. Conventional tillage promotes both disruption of soil aggregates, exposing aggregate-protected SOC to microbial activity, and enhanced conditions for faster SOC mineralization and erosion [2]. These losses can be up to two-thirds of the initial SOC levels within a few decades after conversion to arable land use and up to a quarter within the first years after a single renovation event, particularly in the topsoil [3,4,5]. With these losses, soil fertility is threatened, and high amounts of carbon dioxide (CO2) and other greenhouse gases (GHGs) of greater warming potential are released, contributing to global warming [6,7].
To counter this problem, ceasing the use of CT and replacing it with less intensive tillage methods, such as reduced (RT) and no-till (NT) during conversion and renovation, has been recommended as a suitable method not only to reduce losses but also to conserve SOC and soil fertility [8,9]. Despite these claims, animated debates are arising concerning the efficacy of NT over CT in promoting SOC sequestration in arable systems. Several studies suggest that replacing CT with NT could increase SOC stocks by at least 0.3 Mg C ha−1 year−1 in the topsoil (<30 cm) or by 4 to 10% across temperate and tropical climates [10,11]. Other studies argue that this increase is limited to the top 10 cm soil, and thus, these apparent benefits disappear when larger soil depths are taken into account, resulting in a simple redistribution or stratification of SOC and soil mass without real gains in SOC [12,13,14,15].
While contrasting evidence exists for arable land, research investigating the use of NT for SOC conservation during grassland conversion and renovation for the top and subsoil is still very limited and non-conclusive. Research conducted under drought conditions or low SOC levels found that NT can conserve SOC after grassland conversion at the same level as an unplowed adjacent grassland site, at 30 cm soil depth [9,16]. Nonetheless, other research works combining results from different sites reported similar SOC losses using both NT and CT after conversion of natural sites to arable land when the subsoil was considered at 60 cm soil depth [12,13].
Subsoil SOC can contribute one-third to more than half of the top-meter SOC stocks [17,18]. It can be affected by grassland conversion through changes in vegetation, e.g., when perennials are displaced by annual crops with lower C inputs (Cinput), particularly from belowground [19,20,21]. Perennials can mobilize high Cinput through the soil profile by growing deep taproots and by high rhizodeposition production, favoring an increase in SOC [22,23,24,25,26]. Increasing subsoil SOC has been observed in the conversion of arable land to grassland, whereas the opposite has been reported in grassland conversions to arable or to other land uses [5,20,27]. Subsoil SOC can also be affected by the tillage method, as this affects both crop residue incorporation and root growth patterns. While CT incorporates crop residues deeper in the soil and supports root penetration, favoring Cinput in the subsoil, NT accumulates these on or near the soil surface, limiting Cinput only to the topsoil [12,28,29,30]. Many studies have been limited to determining SOC changes in the topsoil where plowing takes place [15,31]. This approach neglects the SOC changes occurring in the subsoil, creating high uncertainties regarding the efficacy of replacing CT with NT as a method to mitigate SOC losses following grassland removal [8,15,32]. It is thus necessary to assess the effects of grassland conversion and renovation using different tillage methods on SOC across the top-meter profile, as both the top and subsoil SOC can turn into important sink or source of CO2 in different ways depending on the management applied [12,19,22].
In this study, we investigated the effects of the tillage method on SOC dynamics after conversion to silage maize and after renovation of a 10-year-old permanent grassland for the 0–90 cm soil depth over a 7-year period. We hypothesized that (1) NT can reduce SOC losses after grassland conversion to arable land; (2) the conversion and renovation of permanent grassland negatively affects the subsoil SOC; and (3) SOC changes are negatively affected by the reduction in annual soil Cinput in the unplowed systems.

2. Materials and Methods

2.1. Study Site

The long-term experiment was conducted at the research farm “Hohenschulen” of the Christian-Albrechts-University of Kiel in northwest Germany (54°19′07′′ N; 9°59′32′′ E; 30 m a.s.l.). The climate is maritime, with a mean annual temperature of 9.4 °C and a mean annual precipitation of 759 mm (1991–2020). The dominant soil type is stagnic Luvisol (according to the FAO Guidelines for soil description, 2006) with a sandy loam to loam texture (Table 1).
Before the study, a 10-year-old permanent grassland experiment organized in a split-plot design was located on the field between 2005 and 2014, comparing different cutting frequencies (3 vs. 5 cuts year−1) and N rates (0 vs. 360 kg N ha−1 year−1) in 6 × 6 m plots. The seeded grassland consisted of a multi-species mixture of grasses and legumes, including Lolium perenne, Festuca pratensis, Poa pratensis, Phleum pratense, Dactylis glomerata, Trifolium repens and Medicago sativa. By 2014, the botanical composition of the swards was dominated by Medicago sativa and Trifolium repens from the legume species and by Dactilys glomerata and Lolium perenne from the grass species in the 3- and 5-cut treatments, respectively [23,34]. Despite the differences in botanical composition, no significant differences in SOC stocks were found between the treatments at different soil depths prior to initiation of the experiment (Supplementary Table S1).
In spring 2015, three cropping systems were introduced into the grassland sward randomly. Randomly, plots of 6 × 6 m were either converted to continuous silage maize Zea Mays L. (CM), using either CT or NT (named CT-CM and NT-CM, respectively), or renovated and re-seeded with Lolium perenne and Trifolium repens using NT (NT-GR). The undisturbed permanent grassland served as the control treatment (GC). Further, the systems were either non-fertilized (N0) or fertilized (N1) with calcium ammonium nitrate (CAN) according to an expected demand in the silage maize and grassland systems of 180 and 380 kg N ha−1 year−1, respectively. The number of replicates per treatment used after 2015 is provided in Appendix A (Table A1).
During conversion and renovation, all plots except for the GC were sprayed with glyphosate-containing herbicides (Round-Up Ultra and Round-Up Powerflex). Plots with the CT treatment were subsequently rotovated (10 cm depth), moldboard plowed (25–30 cm depth) and treated with rotary harrow for seedbed preparation. The plots with the NT treatment received no soil cultivation. Seeding in the CM systems was performed with the maize variety Ronaldinio (KWS), using a direct drilling seeder (Horsch, Germany), at a rate of 10.5 seeds m−2, a row distance of 0.75 m and a sowing depth of 5 cm. In the GR treatment, a 5:1 seed mixture of Lolium perenne to Trifolium repens was sown using a direct seeder (Vredo, Netherlands) at a rate of 40 kg ha−1, once, in spring 2015. Nutrients in addition to N (P, K, Mg and S) were supplied in sufficient amounts at the beginning of the growing season to avoid nutrient deficiencies. Further details about the management can be found in Struck et al. (2019) [35].

2.2. Measurements

2.2.1. Soil Sampling and Carbon Analysis

Soil sampling was performed annually from 2014 (the year prior to conversion) to 2021 at the end of each growing season in fixed soil-depth increments: 0–30, 30–60 and 60–90 cm. To do so, samples from each plot were collected using a hydraulic-driven soil auger (inner Ø 2 cm). Afterward, soil samples were oven dried at 30 °C and sieved to pass a 2 mm mesh size for further analyses. The C and N contents (%) were determined using a CN elemental analyzer (Vario Max CN, Elementar) (ISO). No carbonates were detected within the soil samples after testing with HCl. Thus, the estimated total C was equal to SOC. The soil sample results from 2018 had to be discarded due to technical problems.
Soil bulk density measurements were conducted for the years 2014, 2015, 2016 and 2020 in all the treatments in the topsoil (0–30 cm). To do so, four soil cores from randomly chosen plots were collected using a steel cylinder (Ø = 5.1 cm; L = 5 cm), according to DIN ISO 11272 (HBU, 1998). No significant differences in soil bulk densities were found between the treatments and the years. Thus, a mean soil bulk density of 1.6 g cm−3 was adopted in all treatments in the 0–30 cm soil depth. For the 30–60 and 60–90 cm soil depths, soil bulk densities of 1.77 and 1.70 g cm−3 were used as reported in Zink et al. (2010) [33] for the study site. These values were adopted, as no significant differences between tillage treatments were reported in the subsoil at the study site by Hartmann et al. (2012) [36]. The SOC stocks calculations for each soil depth were performed following Equation (1):
SOC   [ Mg   C   ha 1 ] =   C   content   [ % ] × Soil   Bulk   Density   [ g   cm 3 ] × Soil   Depth   [ cm ]

2.2.2. Plant biomass sampling and C input calculations

The annual crop yields (YP) in the CM systems were obtained using a two-row plot harvester (Kemper 1200, Haldrup, Logstor, Denmark) at the silage maturity stage, at a stubble height of 25 cm. After harvesting, the stubbles were left on the soil surface in the NT system, whereas in the CT system, these remained on the field until their incorporation via tillage during spring. The grassland systems were harvested with a cutting frequency of 4 cuts year−1 (May, July, August and October), using a forage plot harvester (Haldrup, Logstor, Denmark) at a cutting height of 5 cm. Dry matter yields were determined from fresh matter subsamples that were oven dried at 58 °C to estimate their humidity content. The mean YP measured for the period 2015–2021 are provided in Supplementary Table S2.
To obtain yield-based allocation coefficients, total above- (AGB) and belowground biomass (BGB) measurements conducted during and after the growing seasons (March–October and November–March, respectively) by Struck (2018) [37] in the years 2015 and 2016 were used. During the growing season, AGB samples were taken at each grassland cutting and maize harvest from an area of 0.25 and 1 m2 in 4 plots, at a cutting height of 5 and 25 cm, respectively, and separated into harvestable biomass and stubbles. In the same plots, BGB was measured using the ingrowth core method [38], starting from mid-March and following the grassland cutting intervals (every 6 weeks). Roots were washed from the soil using a hydro-pneumatic elutriation system [39] and collected using a mesh size of 0.63 mm. After drying the samples, the ash-corrected dry matter and carbon contents were determined for both AGB and BGB samples. Harvest (HI) and stubble indices (SI) were calculated by the ratio of YP and stubbles to total AGB, respectively. The R:S was determined by the ratio between the annual sums of the ash-corrected BGB and AGB. To determine the BGB produced after the growing season, the same method was applied as explained above, using the cumulative biomass of three sampling dates (November, February and March). A detailed description of the measurements is available in Struck (2018) [37].
The annual Cinput were estimated in each cropping system by using the measured YP and yield-based allocation coefficients to calculate the Cinput produced above- (AG) and belowground (BG) during the growing season, and by adding fixed inputs to account for the residues produced AG and BG in the grassland systems after the growing season (AW–AGinput and AW–BGinput, respectively). For the AG residues produced during the growing season, stubbles were calculated using the estimated SI in the CM systems, whereas in the grassland systems, the harvest residues and litter deposition were calculated by using a SI equal to 15% on a DM basis, as proposed by Bolinder et al. 2007 [40]. The BG residues produced in the same season were calculated using the mean ash-corrected R:S ratios measured in the first two years. For the fixed inputs in the grassland systems, the BGB measured after the growing season in the first two experimental years was used. Since AGB was not measured in this period, we assumed these inputs using the AGB measured after the growing season from a similar grassland experiment located at a nearby experimental farm [26,41]. To account for the extra inputs derived from rhizodeposition, an additional 50% of the total BGB produced during and after the growing season was added to the BG Cinput calculations, as suggested by Pausch and Kuzyakov (2018) [42]. Ash-corrected dry matter values were used for the Cinput calculations, as these are preferred to total dry matter [40]. Carbon content of 48% was calculated on the ash-corrected AG and BG plant material. The steps and equations used to calculate the annual soil Cinput in each cropping system using the YP, yield-based allocation coefficients and fixed inputs are provided in Supplementary Tables S3 and S4.

2.3. Statistical Analysis

To determine their effects on SOC dynamics, the factors cropping system, N rate and soil depth) were tested for their interaction with time (year) using analysis of covariance (ANCOVA) in a linear mixed-effects model. To do so, the nlme package [43] was used in R statistical software [44]. Normal distribution and heteroscedasticity of the data with respect to the different treatments were assumed in the statistical models after graphical residual analysis. Two separate ANCOVAs were performed: one considering soil depths and another summing the SOC stocks across the fixed depth increments (0–90 cm). When the study factors significantly interacted with time, separate regressions were estimated and tested for whether the ΔSOC statistically differed from zero. Additionally, a Tukey contrast was used for mean ΔSOC comparisons using the multcomp package [45]. All the tested effects were considered significant with a p-value < 0.05.
The relationships between the ΔSOC and the corresponding annual soil Cinput were fitted to a linear and a non-linear asymptotic regression model in the 0–30 cm soil depth increment, using the observations of the unplowed systems (GC, NT-GR and NT-CM). Both approaches estimate the SOC sequestration efficiency and the limits to soil C sequestration [46,47].

3. Results

3.1. Effects of Grassland Conversion and Renovation on Soil Organic Carbon Dynamics

A significant interaction between the factors soil depth, cropping system and time was found after grassland conversion and renovation in this study (p < 0.0001, Supplementary Table S5). Similarly, across soil depths (0–90 cm), cropping system significantly interacted with time (p < 0.0001, Supplementary Table S6). No interaction was found between N rate and time. Thus, the ΔSOC of the cropping systems were pooled over the N rates.
At the 0–30 cm depth, significant SOC losses over time were observed in CM for both tillage methods (p < 0.0001, Figure 1a). By using NT during conversion, SOC losses were reduced from −2.5 to −1.8 Mg C ha−1 year−1, for a 28% reduction in annual SOC losses in comparison to the CT treatment. By the end of the 7-year period, the SOC levels in NT-CM and CT-CM were 14 and 19 Mg ha−1 lower than the levels measured before conversion, for a reduction in 18 and 26%, respectively. In contrast, SOC remained unchanged in the topsoil of NT-GR, similar to the GC. In the subsoil (30–60 and 60–90 cm), no changes in SOC were observed after conversion or renovation by using any of the tillage methods. Across the soil depth increments (0–90 cm, Figure 1b), the ΔSOC were significantly different to zero in both NT-CM and CT-CM (p < 0.0001), with SOC losses of −1.58 and −2.63 Mg C ha−1 year−1, respectively. In contrast, in the grassland systems, SOC remained stable, with non-significant ΔSOC.
The mean comparisons of ΔSOC between the cropping systems revealed significant differences overall between the grassland and maize systems, and between tillage methods within the maize systems; the grassland systems did not differ in ΔSOC (Figure 2a,b). These differences occurred mainly in the topsoil, whereas in the subsoil, no significant differences between cropping systems were observed.

3.2. Effects of Soil C Inputs on Topsoil SOC Dynamics

The total annual Cinput decreased by 73% (from 6.7 to 1.8) after conversion to silage maize, whereas after grassland renovation, the total Cinput was comparable to the GC (6.7 compared to 6.3) (Figure 3). Lower Cinput after conversion occurred mainly due to a 70% reduction in BG Cinput (from 5.5 to 1.5) compared to a 33% reduction in AG Cinput (from 1.2 to 0.8) across the N rates (Figure 3a,b) due to the altered root-to-shoot ratios of the following crop. After renovation, BG Cinput decreased only by 7% (from 5.5 to 5.1), whereas AG Cinput remained unchanged (1.2 Mg C ha−1 year−1). Differences between the annual Cinput of the two tillage methods in the CM systems were marginal.
In all four systems, N fertilization increased AG Cinput, while it reduced BG Cinput except for CT-CM. The grassland systems decreased their annual Cinput from 8.1 to 5.2 (35%) and from 7.5 to 5 Mg ha−1 (33%) in GC and NT-GR, respectively. The CM systems, on the other hand, did not increase their annual Cinput with N fertilization in the NT-CM, whereas in the CT-CM, a 0.4 Mg ha−1 increase was observed.
The relationship between ΔSOC with annual soil Cinput for the unplowed systems (GC, NT-GR and NT-CM) showed both a significant linear and non-linear relationship in the topsoil (Figure 4). A better fit was obtained using a non-linear asymptotic function compared to the linear function due to a lower score in the Akaike information criteria (AIC). Annual SOC losses (−ΔSOC) were reduced by increasing the annual Cinput, up to ~5 Mg ha−1. Additional Cinput did not produce higher ΔSOC, resulting in SOC equilibrium at 1.6% C content. A reduction in the annual Cinput in the grasslands at N1 was not different to the ΔSOC at N0. In the maize systems, both N rates produced a similar annual Cinput and mean negative ΔSOC (SOC losses). Due to the absence of ΔSOC, any effects of the Cinput across the unplowed systems in the subsoil were discarded.

4. Discussion

4.1. No-Till (NT) Effects on Topsoil SOC Dynamics after Grassland Conversion and Renovation

In this study, it was hypothesized that NT can reduce SOC losses after grassland conversion to arable land and maintain SOC stocks when used for grassland renovation. In agreement with the initial hypothesis, annual SOC losses were significantly reduced by almost 30% when NT was used instead of CT, despite the substantial losses observed (Figure 1).
These results present both similarities and differences with those of previous research works in different site conditions. While our results show SOC losses using both tillage methods, a 3-year study by Franzluebbers and Stuedemann (2008) [16] found SOC losses in CT but not in NT during grassland conversion. Moreover, the annual SOC losses in the CT were higher in our study with −2.5 compared to −1.2 Mg ha−1 in their study, whereas the absolute difference between CT and NT in our study was only 0.7 Mg ha−1 after seven years of land use change to arable land. These results were obtained at a site with a similar soil texture as in the present study but with much higher precipitation (1250 vs. 847 mm year−1) and much lower initial SOC stocks (~35 vs. ~74 Mg ha−1). In another study, no SOC losses were found 6 years after grassland was converted to continuous maize using NT, at a site with high initial SOC stocks (~91 Mg ha−1) but with a different soil texture (silt loam vs. sandy loam) combined with an extended drought period [9].
Environmental differences can be crucial in this context, as the water potential in the soil affects the turnover rates of SOC. Periods of drought can bring microbial activity to its minimum, while wet conditions increase it [5], regardless of the tillage method used. Moreover, these rates can be proportional to the initial SOC content, the soil C balance and the amount of silt and clay protecting the SOC [30]. Soils high in SOC are often rich in labile carbon fractions and are thus more susceptible to high decay rates during the first years of grassland conversion compared to soils with low initial SOC and higher proportions of inert SOC. This could potentially explain the differences between the studies. Accordingly, the turnover of SOC in the soil depends on many factors, making mechanistic dynamic models necessary in order to predict potential losses after management changes. A previously applied soil carbon model calibrated with soil carbon and yield data from the presented experimental site confirmed the decay of SOC regardless of the tillage method used [28].
Contrary to the results observed for grassland conversion, using NT during grassland renovation (NT-GR) conserved SOC stocks and showed no differences in the 0–30 cm soil depth compared to the GC (Figure 1 and Figure 2). Research works using NT for SOC conservation during grassland renovation are scarce. The existing studies are inconclusive regarding the advantages of NT over CT during grassland renovation. For instance, Linsler et al. (2013) [48] reported no differences in SOC stocks between grassland that was renovated using CT and the adjacent undisturbed grassland system sampled 5 years after renovation. Similarly, no changes in SOC stocks over time in either NT or CT during grassland renovation in the 0–30 cm soil depth were observed, sampled 3 years after renovation [49]. In contrast, a 25% reduction in SOC stocks was observed after 2.5 years with a single use of CT for grassland renovation compared to the adjacent undisturbed grassland in a study located in a temperate wet region [3]. Our study did not include a grassland renovation treatment using CT; however, high decay rates in the first year after grassland renovation using CT were found close to our experimental site at similar soil conditions in a long-term grassland renovation trial [50]. Thus, we can assume that replacing CT with NT for renovation can significantly decrease the risks of high carbon losses in the initial phase of management and, consequently, the risk of high nitrogen emissions in the environment [6]. By comparing it with GC, it can be deduced that SOC balance was not affected by the grassland renovation measure using NT in the mid-term.

4.2. Land Use Change and Tillage Effects on the Subsoil SOC Dynamics

It was hypothesized that due to changes in vegetation and residue distribution caused by tillage treatments, the conversion and renovation of permanent grassland negatively affect subsoil SOC. Contrary to what was expected, in this study, subsoil SOC was not affected by the grassland conversion or renovation nor by the tillage method over the measured period (Figure 1 and Figure 2). The same results were obtained by Jarvis et al. (2017) [51], who observed significant effects of different ley durations, varying from 1 to 5 years, in 6-year rotations on topsoil SOC stocks but no effects on the subsoil SOC, after 54 years of establishment. Similarly, Syswerda et al. (2011) [52] found no differences in subsoil SOC comparing annual cropping systems using different tillage methods with perennial crops and successional systems, after 12 years of establishment. Meanwhile, Luo et al. (2010) observed similar SOC losses between NT and CT after conversion of natural to arable land in the subsoil [12].
Despite the higher Cinput usually observed in soils under both permanent grassland and the grassland phase of a ley-arable rotation compared to arable soils [20,26,53], the absence of SOC changes in this and other studies supports the existence of factors other than management and annual Cinput that exert a great influence on the subsoil SOC dynamics. The absence of SOC changes observed over time might be attributed to the influence of site-specific conditions associated with soil texture, structure and pedological processes that affect SOC in the soil profile [18,54,55]. The amount of soil mineral-binding particles that protect SOC from microbial degradation and the existence of anoxic conditions that inhibit microbial activity can increase the turnover time of SOC in the subsoil [56]. Likewise, the presence of high soil bulk densities impeding root penetration due to both anoxic conditions and high soil compaction might also limit the influence of Cinput on the subsoil. For instance, changes in subsoil SOC stocks for a loamy site but not for a clayey site were reported after using similar rotations and N fertilization rates at different sites [57]. Similarly, after observing the subsoil SOC dynamics of 24 sites, Poeplau and Don (2013) [55] found the processes affecting subsoil SOC varied between the sites. At our site, the presence of a stagnic subsoil (Table 1) could have impeded both the decomposition of SOC and the root growth of both the grassland and maize systems equally, keeping SOC relatively stable over time. In addition to this, the nature of C compounds in the subsoil is highly recalcitrant, dominated mostly by slow and passive pools of mean residence time that range from several hundreds to thousands of years [18]. This means that any change from vegetation and thereby the amount of fresh Cinput entering the soil will have a faster impact on the topsoil than on the subsoil, as the decoupling of Cinput with the SOC cycling increases with depth [58]. This makes the subsoil SOC less sensitive to land use change or to management than the topsoil SOC [55,58]. Assuming the introduction of deep rooting crop species into shallow-rooted systems or inversely regulating subsoil SOC differently through the different belowground Cinput sources, it would probably require a longer time to observe changes occurring due to grassland conversion than the 7-year period of this study. Based on these observations, the potential site-specific effects on subsoil SOC turnover time and on translocation of freshly added Cinput would require to be analyzed in addition to the treatment effects to draw definite conclusions [18,55,58].

4.3. Plant-Derived C Inputs Affect Annual SOC Changes (ΔSOC)

It was hypothesized that SOC changes are negatively affected by the reduction in annual soil Cinput in the unplowed systems. Both significant linear and non-linear asymptotic relationships between ΔSOC and annual Cinput were found across the unplowed systems, with higher SOC losses (negative ΔSOC) related to lower annual Cinput obtained in the maize systems after conversion and with SOC conservation in the grassland systems related to similar annual Cinput (Figure 3 and Figure 4).
The use of NT during conversion and renovation offers the opportunity to compare both the reduction and the maintenance of Cinput to the levels existing before a management change. According to studies, the equilibrium (i.e., steady state) levels of SOC are determined by the balance between the inputs and losses of organic C [59]. This means that a reduction in Cinput to soil could lead to SOC losses even when decomposition rates remain unaltered due to the absence of tillage [60]. In support of this, the SOC equilibrium maintained in the grassland systems was interrupted mainly as a result of a reduction in BG Cinput, which represented 70 to 90% of the total inputs in the grassland systems (Figure 3c). As opposed to conversion, SOC levels after grassland renovation were maintained due to similar Cinput. These results are consistent with those found in a 3-year study where the conversion of grasslands to arable land using NT led to a reduction in BG Cinput of 43% in the top meter and to the degradation of the labile C pools, which is used as an early indicator for the deterioration in soil quality [21]. Another study comparing NT continuous wheat and NT wheat–fallow systems with a native undisturbed grassland also observed 17 and 22% lower SOC levels at the top 30 cm in the annual systems compared to the undisturbed grassland, respectively, due to the 90% lower BG Cinput [61]. This suggests that, with the use of NT alone to minimize the turnover rates, it is unlikely that SOC can be conserved to the levels seen prior to grassland conversion unless the Cinput remains unchanged or minimally reduced. On the other hand, to increase the SOC stocks in cropland soils under NT, it would be necessary to increase the Cinput either by increasing the cropping frequency from single to double cropping [12,62] or by rotating annual crops with crop species high in belowground net primary productivity, such as perennial forages [28,30,63]. This confirms the expected importance of maintaining Cinput levels to conserve SOC stocks after grassland conversion.

5. Conclusions

In this study, the measurement of all relevant sources of C inputs to the soil allowed a deeper insight on SOC regarding the impact of land use change and tillage method during grassland conversion and renovation for the top-meter soil. While no-till (NT) could not conserve SOC after grassland conversion, it strongly reduced its losses by 28% compared to conventional tillage (CT) and conserved SOC during renovation under the conditions of the study site. Subsoil SOC was not affected by the tillage method, which is likely to be attributed to the lack of carbon inputs to the subsoils, as both grass and maize systems expand their rooting systems predominately in the topsoil. Overall, the much lower soil carbon inputs in the maize (1.8) system compared to the grassland (6.7 Mg C ha−1) system was the reason for the observed SOC losses under NT-CM. Since topsoil SOC followed an asymptotic relationship with increasing annual C inputs, a SOC equilibrium was indicated. This strengthens the need to maintain annual C input high in order to conserve SOC even when no-till approaches are implemented. Our results indicate that NT is a suitable option for mitigating SOC losses when a grassland conversion and renovation measure seems unavoidable.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture12081204/s1, Table S1: Analysis of variance (ANOVA) evaluating the effects of cutting frequency and N rate on the soil organic carbon (SOC) stocks at different soil depths (0–30, 30–60 and 60–90 cm) in the 10-year-old grass sward in year 2014. Table S2: Mean annual yields (YP) with standard deviations (S.D.) measured for the period 2015–2021. Table S3: Yield-based allocation coefficients and fixed inputs used to calculate the annual soil C inputs (Cinput) from plants during and after the growing season, respectively. Table S4: Steps and equations used to calculate annual soil C inputs (Cinput) in each cropping system, using the annual yields (YP) from Table S2 and the yield-based allocation coefficients and fixed inputs from Table S3. Table S5: Analysis of covariance (ANCOVA) evaluating the interaction between cropping system, N rate and soil depth with time on soil organic carbon (SOC) stocks after conversion and renovation of 10-year-old permanent grassland. Table S6: Analysis of covariance (ANCOVA) evaluating the interaction between cropping system and N rate with time across soil depth increments (0–90 cm) on soil organic carbon (SOC) stocks after conversion and renovation of 10-year-old permanent grassland. Table S7: Analysis of covariance (ANCOVA) evaluating the influence of the covariate annual soil Cinput on annual SOC changes (ΔSOC) at the 0–30 cm soil depth after conversion and renovation of 10-year-old permanent grassland in the unplowed system (GC, NT-GR and NT-CM).

Author Contributions

Conceptualization, J.D.L.R., A.P., T.R., R.L. and F.T.; Methodology, J.D.L.R., A.P., T.R. and R.L.; Software, J.D.L.R.; Formal analysis, J.D.L.R.; Investigation, J.D.L.R. and R.L.; Data curation, C.K., J.D.L.R. and R.L.; Writing—original draft preparation, J.D.L.R.; Writing—review and editing, A.P., T.R. and F.T; Visualization, J.D.L.R.; Supervision, A.P., T.R. and F.T.; Funding acquisition, T.R. and F.T. All authors have read and agreed to the published version of the manuscript.

Funding

J.D.L.R. is supported by the Evangelisches Studienwerk Villigst foundation under the research program: “Third Ways of Feeding the World” in the form of a doctoral scholarship. This research was financially supported by the Stiftung Schleswig-Holsteinische Landschaft. The Land Schleswig-Holstein provided financial support within the funding programme Open Access Publikationsfonds. Their financial support is gratefully acknowledged.

Acknowledgments

The authors thank M. Hasler for the valuable statistical consultation, HORSCH Maschinen GmbH for the technical assistance and R. Kopp, K. Helmich, M. Hansen, T. Ehmsen and R. Ströh for their expert technical support during the fieldwork and lab analyses.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Appendix A

Table A1. Number of replicates for the different cropping systems and N rates used in the current study, after their introduction into the 10-year-old permanent grassland in 2015. Treatment descriptions are provided in the text.
Table A1. Number of replicates for the different cropping systems and N rates used in the current study, after their introduction into the 10-year-old permanent grassland in 2015. Treatment descriptions are provided in the text.
Cropping Systems from 2015N RatePermanent Grassland (2005–2014)
Cutting Frequency × N RatesTreatment Replicates
3 × 05 × 03 × 3605 × 360
GCN0333312
N1333312
NT-GRN0333312
N1333312
NT-CMN0337720
N1226616
CT-CMN0336618
N1336618
Total number of replicates23233737120

References

  1. Conant, R.T.; Cerri, C.E.P.; Osborne, B.B.; Paustian, K. Grassland management impacts on soil carbon stocks: A new synthesis. Ecol. Appl. 2017, 27, 662–668. [Google Scholar] [CrossRef] [PubMed]
  2. Lal, R. World cropland soils as a source or sink for atmospheric carbon. In Advances in Agronomy; Sparks, D.L., Ed.; Academic Press: Cambridge, MA, USA, 2001; Volume 71, pp. 145–191. [Google Scholar]
  3. Necpálová, M.; Li, D.; Lanigan, G.; Casey, I.A.; Burchill, W.; Humphreys, J. Changes in soil organic carbon in a clay loam soil following ploughing and reseeding of permanent grassland under temperate moist climatic conditions. Grass Forage Sci. 2013, 69, 611–624. [Google Scholar] [CrossRef]
  4. Kayser, M.; Müller, J.; Isselstein, J. Grassland renovation has important consequences for C and N cycling and losses. Food Energy Secur. 2018, 7, e00146. [Google Scholar] [CrossRef]
  5. Guo, L.B.; Gifford, R.M. Soil carbon stocks and land use change: A meta analysis. Glob. Chang. Biol. 2002, 8, 345–360. [Google Scholar] [CrossRef]
  6. Reinsch, T.; Loges, R.; Kluß, C.; Taube, F. Renovation and conversion of permanent grass-clover swards to pasture or crops: Effects on annual N2O emissions in the year after ploughing. Soil Tillage Res. 2018, 175, 119–129. [Google Scholar] [CrossRef]
  7. Velthof, G.L.; Hoving, I.E.; Dolfing, J.; Smit, A.; Kuikman, P.J.; Oenema, O. Method and timing of grassland renovation affects herbage yield, nitrate leaching, and nitrous oxide emission in intensively managed grasslands. Nutr. Cycl. Agroecosyst. 2010, 86, 401–412. [Google Scholar] [CrossRef]
  8. Ogle, S.M.; Alsaker, C.; Baldock, J.; Bernoux, M.; Breidt, F.J.; McConkey, B.; Regina, K.; Vazquez-Amabile, G.G. Climate and Soil Characteristics Determine Where No-Till Management Can Store Carbon in Soils and Mitigate Greenhouse Gas Emissions. Sci. Rep. 2019, 9, 11665. [Google Scholar] [CrossRef] [PubMed]
  9. Follett, R.F.; Varvel, G.E.; Kimble, J.M.; Vogel, K.P. No-Till Corn after Bromegrass: Effect on Soil Carbon and Soil Aggregates. Agron. J. 2009, 101, 261–268. [Google Scholar] [CrossRef]
  10. West, T.O.; Post, W.M. Soil Organic Carbon Sequestration Rates by Tillage and Crop Rotation: A global data analysis. Soil Sci. Soc. Am. J. 2002, 66, 1930–1946. [Google Scholar] [CrossRef]
  11. Ogle, S.M.; Wakelin, S.J.; Buendia, L.; McConkey, B.; Baldock, J.; Akiyama, H.; Kishimoto, A.W.M.; Chirinda, N.; Bernoux, M.; Bhattacharya, S.N.; et al. Chapter 5 Cropland. In 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories; Buendia, E.C., Tanabe, K., Kranjc, A., Baasansuren, J., Fukuda, M., Ngarize, S., Osako, A., Pyrozhenko, Y., Shermanau, P., Federici, S., Eds.; IPCC: Geneva, Switzerland, 2019. [Google Scholar]
  12. Luo, Z.; Wang, E.; Sun, O.J. Can no-tillage stimulate carbon sequestration in agricultural soils? A meta-analysis of paired experiments. Agric. Ecosyst. Environ. 2010, 139, 224–231. [Google Scholar] [CrossRef]
  13. Christopher, S.F.; Lal, R.; Mishra, U. Regional Study of No-Till Effects on Carbon Sequestration in the Midwestern United States. Soil Sci. Soc. Am. J. 2009, 73, 207–216. [Google Scholar] [CrossRef]
  14. Powlson, D.S.; Stirling, C.M.; Jat, M.L.; Gerard, B.; Palm, C.A.; Sanchez, P.A.; Cassman, K.G. Limited potential of no-till agriculture for climate change mitigation. Nat. Clim. Chang. 2014, 4, 678–683. [Google Scholar] [CrossRef]
  15. Gál, A.; Vyn, T.J.; Michéli, E.; Kladivko, E.J.; McFee, W.W. Soil carbon and nitrogen accumulation with long-term no-till versus moldboard plowing overestimated with tilled-zone sampling depths. Soil Tillage Res. 2007, 96, 42–51. [Google Scholar] [CrossRef]
  16. Franzluebbers, A.J.; Stuedemann, J.A. Early Response of Soil Organic Fractions to Tillage and Integrated Crop–Livestock Production. Soil Sci. Soc. Am. J. 2008, 72, 613–625. [Google Scholar] [CrossRef]
  17. Poeplau, C.; Jacobs, A.; Don, A.; Vos, C.; Schneider, F.; Wittnebel, M.; Tiemeyer, B.; Heidkamp, A.; Prietz, R.; Flessa, H. Stocks of organic carbon in German agricultural soils—Key results of the first comprehensive inventory. J. Plant Nutr. Soil Sci. 2020, 183, 665–681. [Google Scholar] [CrossRef]
  18. Rumpel, C.; Kögel-Knabner, I. Deep soil organic matter—A key but poorly understood component of terrestrial C cycle. Plant Soil 2011, 338, 143–158. [Google Scholar] [CrossRef]
  19. VandenBygaart, A.J.; Bremer, E.; McConkey, B.G.; Ellert, B.H.; Janzen, H.H.; Angers, D.A.; Carter, M.R.; Drury, C.F.; Lafond, G.P.; McKenzie, R.H. Impact of Sampling Depth on Differences in Soil Carbon Stocks in Long-Term Agroecosystem Experiments. Soil Sci. Soc. Am. J. 2011, 75, 226–234. [Google Scholar] [CrossRef]
  20. Gentile, R.; Martino, D.; Entz, M. Influence of perennial forages on subsoil organic carbon in a long-term rotation study in Uruguay. Agric. Ecosyst. Environ. 2005, 105, 419–423. [Google Scholar] [CrossRef]
  21. Dupont, S.T.; Culman, S.; Ferris, H.; Buckley, D.H.; Glover, J.D. No-tillage conversion of harvested perennial grassland to annual cropland reduces root biomass, decreases active carbon stocks, and impacts soil biota. Agric. Ecosyst. Environ. 2010, 137, 25–32. [Google Scholar] [CrossRef]
  22. Lal, R.; Follett, R.F. Soil Carbon Sequestration and the Greenhouse Effect, 2nd ed.; Soil Science Society of America: Madison, WI, USA, 2009; ISBN 9780891188506. [Google Scholar]
  23. Nüsse, A.; Linsler, D.; Loges, R.; Reinsch, T.; Taube, F.; Ludwig, B. Effect of grassland harvesting frequency and N-fertilization on stocks and dynamics of soil organic matter in the temperate climate. Arch. Agron. Soil Sci. 2018, 64, 1925–1931. [Google Scholar] [CrossRef]
  24. Johnson, J.M.; Allmaras, R.R.; Reicosky, D.C. Estimating Source Carbon from Crop Residues, Roots and Rhizodeposits Using the National Grain-Yield Database. Agron. J. 2006, 98, 622–636. [Google Scholar] [CrossRef]
  25. Rasse, D.P.; Rumpel, C.; Dignac, M.-F. Is soil carbon mostly root carbon? Mechanisms for a specific stabilisation. Plant Soil 2005, 269, 341–356. [Google Scholar] [CrossRef]
  26. Loges, R.; Bunne, I.; Reinsch, T.; Malisch, C.; Kluß, C.; Herrmann, A.; Taube, F. Forage production in rotational systems generates similar yields compared to maize monocultures but improves soil carbon stocks. Eur. J. Agron. 2018, 97, 11–19. [Google Scholar] [CrossRef]
  27. Deng, L.; Zhu, G.-Y.; Tang, Z.-S.; Shangguan, Z.-P. Global patterns of the effects of land-use changes on soil carbon stocks. Glob. Ecol. Conserv. 2016, 5, 127–138. [Google Scholar] [CrossRef]
  28. Reinsch, T.; Struck, I.J.A.; Loges, R.; Kluß, C.; Taube, F. Soil carbon dynamics of no-till silage maize in ley systems. Soil Tillage Res. 2021, 209, 104957. [Google Scholar] [CrossRef]
  29. Angers, D.A.; Eriksen-Hamel, N.S. Full-Inversion Tillage and Organic Carbon Distribution in Soil Profiles: A Meta-Analysis. Soil Sci. Soc. Am. J. 2008, 72, 1370–1374. [Google Scholar] [CrossRef]
  30. Paul, E.A.; Elliott, E.T.; Paustian, K.; Cote, C.V. (Eds.) Soil Organic Matter in Temperate Agroecosystems: Long-Term Experiments in North America, 1st ed.; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  31. Baker, J.M.; Ochsner, T.; Venterea, R.T.; Griffis, T.J. Tillage and soil carbon sequestration—What do we really know? Agric. Ecosyst. Environ. 2007, 118, 1–5. [Google Scholar] [CrossRef]
  32. Don, A.; Scholten, T.; Schulze, E.-D. Conversion of cropland into grassland: Implications for soil organic-carbon stocks in two soils with different texture. J. Plant Nutr. Soil Sci. 2009, 172, 53–62. [Google Scholar] [CrossRef]
  33. Zink, A.; Fleige, H.; Horn, R. Load Risks of Subsoil Compaction and Depths of Stress Propagation in Arable Luvisols. Soil Sci. Soc. Am. J. 2010, 74, 1733–1742. [Google Scholar] [CrossRef]
  34. Schmeer, M.; Loges, R.; Dittert, K.; Senbayram, M.; Horn, R.; Taube, F. Legume-based forage production systems reduce nitrous oxide emissions. Soil Tillage Res. 2014, 143, 17–25. [Google Scholar] [CrossRef]
  35. Struck, I.J.; Reinsch, T.; Herrmann, A.; Kluß, C.; Loges, R.; Taube, F. Yield potential and nitrogen dynamics of no-till silage maize (Zea mays L.) under maritime climate conditions. Eur. J. Agron. 2019, 107, 30–42. [Google Scholar] [CrossRef]
  36. Hartmann, P.; Zink, A.; Fleige, H.; Horn, R. Effect of compaction, tillage and climate change on soil water balance of Arable Luvisols in Northwest Germany. Soil Tillage Res. 2012, 124, 211–218. [Google Scholar] [CrossRef]
  37. Struck, I. No-Tillage Silage Maize (Zea mays L.) in Ley-Arable Systems–Crop Performance and Environmental Effects under Maritime Climates. Ph.D. Thesis, Christian-Albrechts University of Kiel, Kiel, Germany, 2018. [Google Scholar]
  38. Steingrobe, B.; Schmid, H.; Claassen, N. The use of the ingrowth core method for measuring root production of arable crops–Influence of soil and root disturbance during installation of the bags on root ingrowth into the cores. Eur. J. Agron. 2001, 15, 143–151. [Google Scholar] [CrossRef]
  39. Smucker, A.J.M.; McBurney, S.L.; Srivastava, A.K. Quantitative Separation of Roots from Compacted Soil Profiles by the Hydropneumatic Elutriation System. Agron. J. 1982, 74, 500–503. [Google Scholar] [CrossRef]
  40. Bolinder, M.; Janzen, H.; Gregorich, E.; Angers, D.; VandenBygaart, A. An approach for estimating net primary productivity and annual carbon inputs to soil for common agricultural crops in Canada. Agric. Ecosyst. Environ. 2007, 118, 29–42. [Google Scholar] [CrossRef]
  41. De Los Rios, J.; Poyda, A.; Reinsch, T.; Kluß, C.; Taube, F.; Loges, R. Integrating Crop-Livestock System Practices in Forage and Grain-Based Rotations in Northern Germany: Potentials for Soil Carbon Sequestration. Agronomy 2022, 12, 338. [Google Scholar] [CrossRef]
  42. Pausch, J.; Kuzyakov, Y. Carbon input by roots into the soil: Quantification of rhizodeposition from root to ecosystem scale. Glob. Chang. Biol. 2017, 24, 1–12. [Google Scholar] [CrossRef]
  43. Bates, D.M.; Watts, D.G. Nonlinear Regression Analysis and Its Applications; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1988; ISBN 9780470316757. [Google Scholar]
  44. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  45. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous Inference in General Parametric Models. Biom. J. 2008, 50, 346–363. [Google Scholar] [CrossRef]
  46. Johnson, J.M.F.; Novak, J.; Varvel, G.E.; Stott, D.E.; Osborne, S.L.; Karlen, D.L.; Lamb, J.A.; Baker, J.; Adler, P.R. Crop Residue Mass Needed to Maintain Soil Organic Carbon Levels: Can It Be Determined? BioEnergy Res. 2014, 7, 481–490. [Google Scholar] [CrossRef]
  47. Stewart, C.E.; Paustian, K.; Conant, R.; Plante, A.; Six, J. Soil carbon saturation: Concept, evidence and evaluation. Biogeochemistry 2007, 86, 19–31. [Google Scholar] [CrossRef]
  48. Linsler, D.; Geisseler, D.; Loges, R.; Taube, F.; Ludwig, B. Temporal dynamics of soil organic matter composition and aggregate distribution in permanent grassland after a single tillage event in a temperate climate. Soil Tillage Res. 2013, 126, 90–99. [Google Scholar] [CrossRef]
  49. Gaweł, E.; Grzelak, M. The Influences of Different Methods of Grassland Renovation on the Weight of Post-Harvest Residues and the Abundance of Selected Soil Nutrients. Agronomy 2020, 10, 1590. [Google Scholar] [CrossRef]
  50. Reinsch, T.; Loges, R.; Kluß, C.; Taube, F. Effect of grassland ploughing and reseeding on CO2 emissions and soil carbon stocks. Agric. Ecosyst. Environ. 2018, 265, 374–383. [Google Scholar] [CrossRef]
  51. Jarvis, N.; Forkman, J.; Koestel, J.; Kätterer, T.; Larsbo, M.; Taylor, A. Long-term effects of grass-clover leys on the structure of a silt loam soil in a cold climate. Agric. Ecosyst. Environ. 2017, 247, 319–328. [Google Scholar] [CrossRef]
  52. Syswerda, S.; Corbin, A.; Mokma, D.; Kravchenko, A.; Robertson, G. Agricultural Management and Soil Carbon Storage in Surface vs. Deep Layers. Soil Sci. Soc. Am. J. 2011, 75, 92–101. [Google Scholar] [CrossRef]
  53. Haynes, R. Labile organic matter as an indicator of organic matter quality in arable and pastoral soils in New Zealand. Soil Biol. Biochem. 2000, 32, 211–219. [Google Scholar] [CrossRef]
  54. Mordhorst, A.; Fleige, H.; Zimmermann, I.; Burbaum, B.; Filipinski, M.; Cordsen, E.; Horn, R. Organische Kohlenstoffvorräte von Bodentypen in den Hauptnaturräumen Schleswig-Holsteins (Norddeutschland). Die Bodenkultur. J. Land Manag. Food Environ. 2018, 69, 85–95. [Google Scholar] [CrossRef]
  55. Poeplau, C.; Don, A. Sensitivity of soil organic carbon stocks and fractions to different land-use changes across Europe. Geoderma 2013, 192, 189–201. [Google Scholar] [CrossRef]
  56. Lorenz, K.; Lal, R. The depth distribution of soil organic carbon in relation to land use and management and the potential of carbon sequestration in subsoil horizons. In Advances in Agronomy; Sparks, D.L., Ed.; Elsevier: Amsterdam, Netherlands, 2005; Volume 88, pp. 35–66. [Google Scholar]
  57. Börjesson, G.; Bolinder, M.A.; Kirchmann, H.; Kätterer, T. Organic carbon stocks in topsoil and subsoil in long-term ley and cereal monoculture rotations. Biol. Fertil. Soils 2018, 54, 549–558. [Google Scholar] [CrossRef]
  58. Jobbagy, E.G.; Jackson, R.B. The Vertical Distribution of Soil Organic Carbon and Its Relation to Climate and Vegetation. Ecol. Appl. 2000, 10, 423–436. [Google Scholar] [CrossRef]
  59. Follett, R. Soil management concepts and carbon sequestration in cropland soils. Soil Tillage Res. 2001, 61, 77–92. [Google Scholar] [CrossRef]
  60. Poyda, A.; Wizemann, H.-D.; Ingwersen, J.; Eshonkulov, R.; Högy, P.; Demyan, M.S.; Kremer, P.; Wulfmeyer, V.; Streck, T. Carbon fluxes and budgets of intensive crop rotations in two regional climates of southwest Germany. Agric. Ecosyst. Environ. 2019, 276, 31–46. [Google Scholar] [CrossRef]
  61. Frank, A.; Liebig, M.; Tanaka, D. Management effects on soil CO2 efflux in northern semiarid grassland and cropland. Soil Tillage Res. 2006, 89, 78–85. [Google Scholar] [CrossRef]
  62. Conceição, P.C.; Dieckow, J.; Bayer, C. Combined role of no-tillage and cropping systems in soil carbon stocks and stabilization. Soil Tillage Res. 2013, 129, 40–47. [Google Scholar] [CrossRef]
  63. Johnston, A.E.; Poulton, P.R.; Coleman, K. Chapter 1 soil organic matter: Its importance in sustainable agriculture and carbon dioxide fluxes. Adv Agron. 2009, 101, 1–57. [Google Scholar] [CrossRef]
Figure 1. Mean temporal development of soil organic carbon (SOC) stocks (±S.D.) in the different cropping systems, (a) for the different soil depth increments (0–30, 30–60 and 60–90 cm) and (b) across the fixed soil depth increments (0–90 cm) for the period 2014–2021, pooled over the two N rates. Cropping systems with **** indicate ΔSOC significantly different from zero at p < 0.0001. Estimates represent the number of replicates for each treatment, as described in Table A1. GC is undisturbed grassland control, NT-GR is no-till grassland renovation, NT-CM is no-till continuous silage maize, CT-CM is conventionally tilled continuous silage maize; N0 is non-fertilized and N1 is fertilized with CAN according to an expected demand of 180 and 380 kg N ha−1 year−1 in the silage maize and grassland systems, respectively.
Figure 1. Mean temporal development of soil organic carbon (SOC) stocks (±S.D.) in the different cropping systems, (a) for the different soil depth increments (0–30, 30–60 and 60–90 cm) and (b) across the fixed soil depth increments (0–90 cm) for the period 2014–2021, pooled over the two N rates. Cropping systems with **** indicate ΔSOC significantly different from zero at p < 0.0001. Estimates represent the number of replicates for each treatment, as described in Table A1. GC is undisturbed grassland control, NT-GR is no-till grassland renovation, NT-CM is no-till continuous silage maize, CT-CM is conventionally tilled continuous silage maize; N0 is non-fertilized and N1 is fertilized with CAN according to an expected demand of 180 and 380 kg N ha−1 year−1 in the silage maize and grassland systems, respectively.
Agriculture 12 01204 g001
Figure 2. Mean annual SOC change rates (±95% confidence intervals) in the different cropping systems (a) for the different soil depth increments (0–30, 30–60 and 60–90 cm) and (b) across the soil depths increments (0–90 cm) during the period 2014–2021, pooled over the two N rates. Different letters indicate significantly different mean ΔSOC between cropping systems. Treatment descriptions are provided in Figure 1.
Figure 2. Mean annual SOC change rates (±95% confidence intervals) in the different cropping systems (a) for the different soil depth increments (0–30, 30–60 and 60–90 cm) and (b) across the soil depths increments (0–90 cm) during the period 2014–2021, pooled over the two N rates. Different letters indicate significantly different mean ΔSOC between cropping systems. Treatment descriptions are provided in Figure 1.
Agriculture 12 01204 g002
Figure 3. Mean annual C inputs from above-(AG) and belowground (BG) crop residues (±S.D.) (a,b) and relative values (c) estimated for the period 2015–2021 in the different cropping systems and N rates using annual yields and yield-based allocation coefficients, for the topsoil. Treatment descriptions are provided in Figure 1.
Figure 3. Mean annual C inputs from above-(AG) and belowground (BG) crop residues (±S.D.) (a,b) and relative values (c) estimated for the period 2015–2021 in the different cropping systems and N rates using annual yields and yield-based allocation coefficients, for the topsoil. Treatment descriptions are provided in Figure 1.
Agriculture 12 01204 g003
Figure 4. Soil organic carbon change rates (ΔSOC) expressed as a function of annual soil C input rates for the unplowed systems (GC, NT-GR and NT-CM). Observations indicate the replicates for the different cropping systems and N rates. Treatment descriptions are provided in Figure 1.
Figure 4. Soil organic carbon change rates (ΔSOC) expressed as a function of annual soil C input rates for the unplowed systems (GC, NT-GR and NT-CM). Observations indicate the replicates for the different cropping systems and N rates. Treatment descriptions are provided in Figure 1.
Agriculture 12 01204 g004
Table 1. Soil characteristics at the study site. Adopted from Zink et al. 2010 [33].
Table 1. Soil characteristics at the study site. Adopted from Zink et al. 2010 [33].
Layer (cm)Texture (%)pHSoil TypeHorizon
SandSiltClayin CaCl2
0–285430166.2Highly loamy sandAp
28–455031196.5Medium sandy loamBtg1
45–704933186.5Medium sandy loamBtg2
>70553114-Highly loamy sandBg
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De Los Rios, J.; Poyda, A.; Taube, F.; Kluß, C.; Loges, R.; Reinsch, T. No-Till Mitigates SOC Losses after Grassland Renovation and Conversion to Silage Maize. Agriculture 2022, 12, 1204. https://doi.org/10.3390/agriculture12081204

AMA Style

De Los Rios J, Poyda A, Taube F, Kluß C, Loges R, Reinsch T. No-Till Mitigates SOC Losses after Grassland Renovation and Conversion to Silage Maize. Agriculture. 2022; 12(8):1204. https://doi.org/10.3390/agriculture12081204

Chicago/Turabian Style

De Los Rios, Josue, Arne Poyda, Friedhelm Taube, Christof Kluß, Ralf Loges, and Thorsten Reinsch. 2022. "No-Till Mitigates SOC Losses after Grassland Renovation and Conversion to Silage Maize" Agriculture 12, no. 8: 1204. https://doi.org/10.3390/agriculture12081204

APA Style

De Los Rios, J., Poyda, A., Taube, F., Kluß, C., Loges, R., & Reinsch, T. (2022). No-Till Mitigates SOC Losses after Grassland Renovation and Conversion to Silage Maize. Agriculture, 12(8), 1204. https://doi.org/10.3390/agriculture12081204

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