3.1. Plant and Soil Effects on SOC Formation
Our results indicate that crop rotation had a significant effect on the content and isotope composition of SOC, while there was no significant fertilization effect. The increase in SOC content between 2009 and 2017 was 500–600 kg ha
−1 yr
−1 higher in CR10 (0% maize, 75% leys) compared to CR6 (50% maize, 25% leys), providing strong evidence in support of hypothesis i. In a comparable field experiment on a sandy loam soil in northern Germany [
30], stable SOC contents were found for crop rotations including both grass-clover leys and silage maize in one out of three years and with cattle slurry application. Since SOC content decreased under continuous maize and increased under permanent grassland, the authors concluded that the ley proportion needed to be higher than 33% to significantly enhance SOC stocks, as was the case in our study. Although SOC content in 2020 was higher in CR10, particularly in the topsoil, this difference was not significant and the crop rotation effect was much stronger for δ
13C than for SOC content. This effect was also clearly visible in the subsoil with the difference in δ
13C between CR6 and CR10 being approximately half as large as the difference in the topsoil.
These findings demonstrate the significant contribution of C
4 carbon to SOC, even in the subsoil of CR6, where maize was grown in eight of the 16 years since the start of the experiment. Fertilization significantly increased the amounts of C
4 carbon in the soil, a strong indication that the fertilization with biogas digestate in this system not only increased yields, but also the belowground productivity of maize. This is in contrast to findings by Hirte et al. [
31], who observed no difference in belowground carbon inputs of maize, irrespective of aboveground productivity. Interestingly, the fertilized plots also tended to have higher fractions of maize carbon, particularly in the subsoil. Although this difference was not significant, it could be interpreted (with some caution) that digestate fertilization had a larger effect on carbon inputs to the subsoil from maize than from C
3 crops. Considering the fertilization effect on crop yields, this result is indeed obvious, since fertilization increased maize yields by 112% on average, while winter wheat and clover-grass yields in CR6 only increased by 66 and 18%, respectively. Note that clover-grass in CR6 did not receive digestate directly. Since maize was grown with undersown white clover, which has a shallower rooting system, fertilization likely increased the competitiveness of maize and thereby also led to the formation of a deeper rooting system, resulting in higher carbon inputs to the subsoil.
We found a significant difference in δ
13C of 1.0 (CR6) and 1.5‰ (CR10) between topsoil and subsoil SOC, which must be due to how the site was cultivated in the years before the start of the experiment, since these differences were even higher than the difference between CR6 and CR10 in the topsoil (0.9‰). In the model approach, the optimization of the initial δ
13C resulted in a difference of 1.2‰ at the beginning of the experiment. One main reason for this difference is probably related to the Suess effect [
32], which was also accounted for in our model. The Suess effect causes recent biomass δ
13C to be more negative than old biomass δ
13C under otherwise identical conditions (i.e., identical
13C fractionation due to identical photosynthetic pathways and growing conditions) and implies that the δ
13C of old SOC must be less negative than that of younger SOC. This agrees with our finding that initial δ
13C in the subsoil was higher, where SOC is older and less influenced by the Suess effect, than in the topsoil. The change in δ
13C in atmospheric CO
2 was about 2‰ during the last two centuries, with the rate of change increasing since the 1950s. The difference of 1.2‰ between topsoil and subsoil at the beginning of the experiment is thus not unlikely and very close to the mean difference of 1.1‰ between 0–10 cm (−27.4‰) and 70–100 cm (−26.3‰) in cropland soils without maize and groundwater tables deeper than 200 cm, as determined in the German Agricultural Soil Inventory [
33]. In addition to the Suess effect, several other processes potentially contributed to these differences between topsoil and subsoil:
(i) Maize cultivation began in this area in the early 1960s. As changes due to maize carbon would be faster in the topsoil than in the subsoil, this potentially cancelled out a part of the Suess effect. Due to lower carbon turnover in the subsoil, there is most likely a higher fraction of older maize carbon in the subsoil than in the topsoil, explaining the higher initial subsoil δ
13C. This difference can no longer be measured since the same crops and maize proportions were cultivated over the entire field in the past. (ii) Roots of C
3 plants are known to have higher δ
13C values than shoots [
25,
34]. Given that the subsoil receives carbon mainly from roots while the topsoil also receives a substantial amount of carbon from aboveground crop biomass (particularly in the years before the experiment started when crop residues were left on the field), less negative δ
13C values in the subsoil are not unlikely. (iii) During the decomposition of organic matter, microbial biomass, which forms only part of SOC, may become enriched by 2‰ compared to SOC, while CO
2 becomes depleted [
35]. However, there is usually negligible fractionation of
13C during the transformation of plant residues into SOC [
36,
37]. Due to these different processes, a difference in initial δ
13C between topsoil and subsoil is thus not unlikely. It is, however, impossible to disentangle these effects given that the magnitude of the Suess effect, the root effect, and the decomposition effect would be considerably less than 3‰.
The size of topsoil SOC pools increased in all treatments, despite the fact that all aboveground biomass, except for stubbles, was used for biogas production. These observations were recently described for this field trial by Levin et al. [
19], who analyzed the change in SOC between 2010 and 2017 for all 10 crop rotations. They found that SOC content increased in 96% of the plots, independent of the type of rotation and the fertilization level. The proportion of clover-grass leys and row crops with undersown clover or clover-grass within a crop rotation had positive and highly significant effects on SOC content. The ley effect was approximately twice as large as the undersowing effect, underlining the generally positive impact of increased perenniality from clover-grass leys on carbon sequestration [
18], and given that clover-grass was not grown on that field for several decades before the year 2000. The 0.22% higher SOC content in the 0–30 cm depth of CR10 (75% ley proportion) compared to CR6 (25% ley proportion) found in soil samples from spring 2020 corresponds well with the 0.04% increase in SOC content with every 10% increase in ley proportion identified by Levin et al. [
19]. SOC content typically decreases with increasing soil depth in mineral soils due to lower inputs of plant residues, and this was also observed in our study, with significantly lower SOC content in the subsoil. As documented by Schneider et al. [
33], climate- and soil-related factors become more important for SOC content than land use at depths below 10 cm.
The amounts of C
4 carbon found in the soil of CR6 correspond to increments of 34, 290, 70, and 353 kg C
4-C ha
−1 per maize cropping year in the subsoil and topsoil of the unfertilized and fertilized treatments, respectively. Despite different pedoclimatic conditions, the values for the topsoil are close to the accumulation rates of 200–300 kg C
4-C ha
−1 yr
−1 in the 0–20 cm depth reported for a long-term maize mono-cropping system in southern China [
38]. The resulting carbon turnover, i.e., the increase in the fraction of C
4 carbon in total SOC per maize year, was 0.8, 0.3, 0.9, and 0.5% in our study. The values for the topsoil fit well with the overall average of 0.8% carbon turnover per year in cropland soils with maize found in the German Agricultural Soil Inventory [
33]. In this inventory, however, carbon turnover in the subsoil could not be determined. Irrespective of the fertilization treatment, the enrichment of maize carbon in the topsoil contributed around 20% of the observed overall SOC pool changes of 1409 and 1800 kg ha
−1 yr
−1 in the unfertilized and fertilized treatments in CR6. This means that the C
3 crops contributed 80% of the increase in SOC pool sizes, despite their proportion of the crop rotation only being 50%, which also supports hypothesis ii. However, a C
3 crop was present on the field during maize years since the white clover cover crop was only mulched and regrew as living mulch in maize. Although the productivity of the white clover was probably low, particularly in the fertilized treatment as described above, there was at least a small input of C
3 carbon in maize years.
The amounts of 2.3 (unfertilized) and 2.8 Mg ha
−1 (fertilized) of C
4 carbon found after sixteen years with eight years of maize cropping in the CR6 plots were remarkably low compared to 2.5–3.2 Mg ha
−1 found after only two years of conventionally managed maize in a field trial in northern Germany, established after conversion from permanent grassland [
15]. This difference can be explained by the much higher productivity of maize crops (shoot biomass of 14–27 Mg dry matter ha
−1) on the conventionally managed plots in the study by Reinsch et al. [
15]. Nevertheless, the huge amounts of SOC these systems lost due to the break-up of the grassland sward could not be offset by the carbon inputs from maize.
Despite the positive effect of digestate fertilization on maize carbon in the soil, this effect was low (+31%) compared to the increase in yield (+112%), and there was no overall fertilization effect on the content or isotopic composition of SOC in the two crop rotations and soil depths. In contrast to these results, the analyses of all plots of the field trial by Levin et al. [
19] revealed significantly higher SOC content in fertilized plots. While they found that SOC content was 0.14% higher in fertilized compared to unfertilized plots in autumn 2017, this difference was only 0.06% in the plots sampled for this study in spring 2020. The reason is most likely that the carbon inputs from the two main crops in the crop rotations of this study, legume–grass mixtures and silage maize, had a low sensitivity to fertilization. While the mixtures only showed a small increase in overall productivity, most likely due to a high biological nitrogen fixation of the legume component in the unfertilized treatments [
39], maize reacted mainly with increased aboveground productivity, as indicated by the low fraction of belowground carbon allocation.
These results are in line with the frequently observed effect of increasing shoot-to-root ratios with the increasing availability of belowground resources, such as nutrients and water [
40,
41,
42]. Since almost all the aboveground biomass was removed from the field for biogas production in this cropping system, soil carbon inputs from harvest residues were only slightly increased by fertilization. However, in our model, we assumed that the fraction of belowground carbon allocation remained constant with increasing yield, which is mathematically described by the partitioning factor (pf) and the yield-specific index (ysi). The ysi was only reduced for silage maize in the optimization procedure, while a value of 1 was used for the other crops. However, this does not mean that there was no interaction between above- and belowground carbon allocation for these crops; the factor affecting this interaction is not the aboveground productivity directly, but rather the resource availability belowground. The inter-annual yield variability was much higher than the fertilization effect, particularly in the legume–grass mixtures, suggesting that additional important factors had an impact on productivity, and not just nutrient availability.
In addition to the limited effect of fertilization on root carbon inputs in the two crop rotations, the direct effect of digestate carbon on SOC was very low. This was demonstrated by the result that model performance was improved with lower retention of digestate carbon in soil (CUE
BD). Regarding the mean annual amounts of applied digestate carbon of roughly 1000 and 1100 kg ha
−1 yr
−1 in CR6 and CR10, respectively, this result seems to be rather surprising. However, these inputs correspond to 2% of the total SOC stocks and a detectable and significant increase due to these inputs might need even more time than the 16 years of this experiment’s runtime. Digestate fertilization might also potentially increased the decomposition of SOC (positive priming effect), counteracting most of the additional carbon inputs [
43]. Very low carbon retention from organic fertilizers (liquid dairy manure) was reported in other studies in mono-cropped cereals [
44] and silage maize [
30]. The fertilizer-induced carbon retention was at least 20% in cropping systems including perennials in these studies. This value includes both the direct (fertilizer carbon) and indirect (plant carbon) effects. The observed slopes of annual SOC increments in the fertilized crop rotations of this study were 23 (CR10) and 28% (CR6) higher compared to the unfertilized treatments. Assuming that carbon sequestration at the study site will continue for several years until a new equilibrium is reached, the positive effect of biogas digestate application on SOC stocks might be more pronounced in the future.
3.2. Modelling SOC Turnover
Using a simple carbon balance model, we could demonstrate that, of the plant-related parameters, those affecting the partitioning between above- and belowground carbon had the biggest impact on model results, thereby confirming hypothesis iii. In general, our model of SOC pool turnover explained the measured δ13C values well, even when pre-defined parameter values were used. This was particularly the case for the topsoil in the unfertilized treatments. The effects of both plant carbon on subsoil SOC and digestate carbon on topsoil SOC were overestimated by the preliminary model settings.
Within each of the different plant-related parameter groups (CUE, pf, BCA
0–30, f
HR, ysi), the maize-specific parameter was most sensitive regarding its effect on model error. Interestingly, exactly the same reduction in model error was achieved regardless of whether CUE
SM, pf
SM, or ysi
SM were optimized. In contrast, optimizing the parameter for the allocation of maize carbon between topsoil and subsoil (BCA
0–30_SM) hardly reduced the model error despite its larger sensitivity compared with CUE
SM and pf
SM. The non-measurable yield-specific silage maize index (ysi
SM) was selected in the first step of the optimization procedure due to its extremely high sensitivity, in combination with the strong reduction in model error. In addition, the parameters CUE and pf were intentionally kept constant in this model. The partitioning factor for above- and belowground biomass (pf) was taken from literature [
45] and is based on the distinct difference between annual and perennial crops. Since all inputs of plant biomass to the soil originated from roots and stubbles, we assumed that the carbon use efficiency of this material did not significantly differ between crops.
The reason for the high sensitivities of maize-specific parameters was that the biggest errors in the preliminary model originated from the large overestimation of δ
13C in fertilized CR6. Thus, the observed large fertilization effect on maize yields resulted in the estimates of maize carbon inputs to the soil being too high when a constant partitioning to roots was assumed, although this was already much lower compared to the perennial crops. The reduced yield-specific silage maize index resulted in decreasing fractions of belowground carbon allocation with increasing maize yield, and thus a strongly improved model fit for the fertilized CR6. In contrast to Hirte et al. [
31], who proposed yield-independent functions for estimating belowground carbon inputs for both maize and wheat, our results suggest crop-specific functions are needed to describe the interactions between yield and belowground carbon allocation.
In the case of maize, our results highlight the minor effect of increased maize yields on soil carbon inputs, as also shown in other studies [
46], and confirm the limited contribution maize crops make to the reproduction of soil organic matter and thus carbon sequestration [
13]. However, it is very likely that maize, as a valuable crop for forage and biogas production, will remain one of the prevalent crops in Germany, particularly in dairy regions. Therefore, more sustainable maize cropping systems need to be developed. These should involve maize cultivation in rotations with grass-clover leys [
15], rather than maize mono-cropping. Our results show that significant carbon sequestration can also be achieved in maize cropping systems if maize is undersown with perennial crops and permanent soil cover is ensured. These systems represent a compromise since undersown clover-grass improves the capacity for carbon sequestration while maize improves biogas yields. To increase the attractiveness of these systems to farmers, research should focus on the effects of different living mulches and undersown crops on maize yields, with a focus on deep-rooting legumes (e.g., lucerne) and herbs (e.g., chicory), as it was shown that these can increase yields of subsequent crops via nutrient transfer from deeper soil layers to the topsoil [
47]. The positive effect of such systems will be even larger on sloping sites because undersown leys strongly reduce soil losses due to erosion, and in turn SOC losses. This not only reduces soil losses during the ley phase but also during the subsequent maize phase [
48,
49].