Predicting spatiotemporal soil organic carbon responses to management using EPIC-IIASA meta-models.

The management of Soil Organic Carbon (SOC) is a critical component of both nature-based solutions for climate change mitigation and global food security. Agriculture has contributed substantially to a reduction in global SOC through cultivation, thus there has been renewed focus on management practices which minimize SOC losses and increase SOC gain as pathways towards maintaining healthy soils and reducing net greenhouse gas emissions. Mechanistic models are frequently used to aid in identifying these pathways due to their scalability and cost-effectiveness. Yet, they are often computationally costly and rely on input data that are often only available at coarse spatial resolutions. Herein, we build statistical meta-models of a multifactorial crop model in order to both (a) obtain a simplified model response and (b) explore the biophysical determinants of SOC responses to management and the geospatial heterogeneity of SOC dynamics across Europe. Using 5600 unique simulations of crop growth from the gridded Environmental Policy Integrated Climate-based Gridded Agricultural Model (EPIC-IIASA GAM) covering 86,000 simulation units across Europe, we build multiple polynomial regression ensemble meta-models for unique combinations of climate and soil across Europe in order to predict SOC responses to varying management intensities. We find that our biophysically-explicit meta models are highly accurate (R2 = 0.97) representations of the full mechanistic model and can be used in lieu of the full EPIC-IIASA GAM model for the estimation of SOC responses to cropland management. Model stratification by means of climate and soil clustering improved the performance of the meta-models compared to the full EU-scale model. In regional and local validations of the meta-model predictions, we find that the meta-models largely capture broad SOC dynamics such as the linear nature of SOC responses to residue application, yet they often underestimate the magnitude of SOC responses to management. Furthermore, we find notable differences between the results from the biophysically-specific models throughout Europe, which point to spatially-distinct SOC responses to management choices such as nitrogen fertilizer application rates and residue retention that illustrate the potential for these models to be used for future management applications. While more accurate input data, calibration, and validation will be needed to accurately predict SOC change, we demonstrate the use of our meta-models for biophysical cluster and field study scale analyses of broad SOC dynamics with basically zero fine-tuning of the models needed. This work provides a framework for simplifying large-scale agricultural models and identifies the opportunities for using these meta-models for assessing SOC responses to management at a variety of scales.

Dynamic soil functions assessment employing land use and climate scenarios at regional scale.

Soils as key component of terrestrial ecosystems are under increasing pressures. As an advance to current static assessments, we present a dynamic soil functions assessment (SFA) to evaluate the current and future state of soils regarding their nutrient storage, water regulation, productivity, habitat and carbon sequestration functions for the case-study region in the Lower Austrian Mostviertel. Carbon response functions simulating the development of regional soil organic carbon (SOC) stocks until 2100 are used to couple established indicator-based SFA methodology with two climate and three land use scenarios, i.e. land sparing (LSP), land sharing (LSH), and balanced land use (LBA). Results reveal a dominant impact of land use scenarios on soil functions compared to the impact from climate scenarios and highlight the close link between SOC development and the quality of investigated soil functions, i.e. soil functionality. The soil habitat and soil carbon sequestration functions on investigated agricultural land are positively affected by maintenance of grassland under LSH (20% of the case-study region), where SOC stocks show a steady and continuous increase. By 2100 however, total regional SOC stocks are higher under LSP compared to LSH or LBA, due to extensive afforestation. The presented approach may improve integrative decision-making in land use planning processes. It bridges superordinate goals of sustainable development, such as climate change mitigation, with land use actions taken at local or regional scales. The dynamic SFA broadens the debate on LSH and LSP and can reduce trade-offs between soil functions through land use planning processes.

Verifiable soil organic carbon modelling to facilitate regional reporting of cropland carbon change: A test case in the Czech Republic.

Regional monitoring, reporting and verification of soil organic carbon change occurring in managed cropland are indispensable to support carbon-related policies. Rapidly evolving gridded agronomic models can facilitate these efforts throughout Europe. However, their performance in modelling soil carbon dynamics at regional scale is yet unexplored. Importantly, as such models are often driven by large-scale inputs, they need to be benchmarked against field experiments. We elucidate the level of detail that needs to be incorporated in gridded models to robustly estimate regional soil carbon dynamics in managed cropland, testing the approach for regions in the Czech Republic. We first calibrated the biogeochemical Environmental Policy Integrated Climate (EPIC) model against long-term experiments. Subsequently, we examined the EPIC model within a top-down gridded modelling framework constructed for European agricultural soils from Europe-wide datasets and regional land-use statistics. We explored the top-down, as opposed to a bottom-up, modelling approach for reporting agronomically relevant and verifiable soil carbon dynamics. In comparison with a no-input baseline, the regional EPIC model suggested soil carbon changes (~0.1-0.5 Mg C ha-1 y-1) consistent with empirical-based studies for all studied agricultural practices. However, inaccurate soil information, crop management inputs, or inappropriate model calibration may undermine regional modelling of cropland management effect on carbon since each of the three components carry uncertainty (~0.5-1.5 Mg C ha-1 y-1) that is substantially larger than the actual effect of agricultural practices relative to the no-input baseline. Besides, inaccurate soil data obtained from the background datasets biased the simulated carbon trends compared to observations, thus hampering the model's verifiability at the locations of field experiments. Encouragingly, the top-down agricultural management derived from regional land-use statistics proved suitable for the estimation of soil carbon dynamics consistently with actual field practices. Despite sensitivity to biophysical parameters, we found a robust scalability of the soil organic carbon routine for various climatic regions and soil types represented in the Czech experiments. The model performed better than the tier 1 methodology of the Intergovernmental Panel on Climate Change, which indicates a great potential for improved carbon change modelling over larger political regions.

African crop yield reductions due to increasingly unbalanced Nitrogen and Phosphorus consumption.

The impact of soil nutrient depletion on crop production has been known for decades, but robust assessments of the impact of increasingly unbalanced nitrogen (N) and phosphorus (P) application rates on crop production are lacking. Here, we use crop response functions based on 741 FAO maize crop trials and EPIC crop modeling across Africa to examine maize yield deficits resulting from unbalanced N : P applications under low, medium, and high input scenarios, for past (1975), current, and future N : P mass ratios of respectively, 1 : 0.29, 1 : 0.15, and 1 : 0.05. At low N inputs (10 kg ha(-1)), current yield deficits amount to 10% but will increase up to 27% under the assumed future N : P ratio, while at medium N inputs (50 kg N ha(-1)), future yield losses could amount to over 40%. The EPIC crop model was then used to simulate maize yields across Africa. The model results showed relative median future yield reductions at low N inputs of 40%, and 50% at medium and high inputs, albeit with large spatial variability. Dominant low-quality soils such as Ferralsols, which are strongly adsorbing P, and Arenosols with a low nutrient retention capacity, are associated with a strong yield decline, although Arenosols show very variable crop yield losses at low inputs. Optimal N : P ratios, i.e. those where the lowest amount of applied P produces the highest yield (given N input) where calculated with EPIC to be as low as 1 : 0.5. Finally, we estimated the additional P required given current N inputs, and given N inputs that would allow Africa to close yield gaps (ca. 70%). At current N inputs, P consumption would have to increase 2.3-fold to be optimal, and to increase 11.7-fold to close yield gaps. The P demand to overcome these yield deficits would provide a significant additional pressure on current global extraction of P resources.

Spatially-explicit land use change emissions and carbon payback times of biofuels under the Carbon Offsetting and Reduction Scheme for International Aviation (CORSIA).

The Carbon Offsetting and Reduction Scheme for International Aviation (CORSIA) requires airlines to offset their greenhouse gas (GHG) emissions above 2019 levels by either buying carbon offsets or using Sustainable Aviation Fuels (SAFs). These are drop-in jet fuels made from biomass or other renewable resources that reduce GHG emissions by at least 10 % compared to kerosene and meet certain sustainability criteria. This study assesses the direct land use change (DLUC) emissions of SAF, i.e., GHG emissions from on-site land conversion from previous uses (excluding primary forests, peatlands, wetlands, and protected and biodiversity-rich areas) into alternative feedstocks, considering spatial variability in global yields and land carbon stocks. The results provide DLUC values and carbon payback times at 0.5-degree resolution for six SAF pathways, with and without irrigation and a medium-input intensity, according to CORSIA sustainability criteria. When excluding CORSIA non-compliant areas, soybean SAF shows the highest mean DLUC factor (31.9 ± 20.7 gCO2/MJ), followed by reed canary grass and maize. Jatropha SAF shows the lowest mean DLUC factor (3.6 ± 31.4 gCO2/MJ), followed by miscanthus and switchgrass. The latter feedstocks show potential for reducing GHG emissions over large areas but with relatively greater variability. Country-average DLUC values are higher than accepted ILUC ones for all pathways except for maize. To ensure the GHG benefits of CORSIA, feedstocks must be produced in areas where not only carbon stocks are relatively low but also where attainable yields are sufficiently high. The results help identify locations where the combination of these two factors may be favourable for low-DLUC SAF production. Irrigated miscanthus offers the highest SAF production potential (2.75 EJ globally) if grown on CORSIA-compliant cropland and grassland areas, accounting for ∼1/5 of the total kerosene used in 2019. Quantifying other environmental impacts of SAFs is desirable to understand sustainability trade-offs and financial constraints that may further limit production potentials.

Climate impacts on global agriculture emerge earlier in new generation of climate and crop models.

Potential climate-related impacts on future crop yield are a major societal concern. Previous projections of the Agricultural Model Intercomparison and Improvement Project's Global Gridded Crop Model Intercomparison based on the Coupled Model Intercomparison Project Phase 5 identified substantial climate impacts on all major crops, but associated uncertainties were substantial. Here we report new twenty-first-century projections using ensembles of latest-generation crop and climate models. Results suggest markedly more pessimistic yield responses for maize, soybean and rice compared to the original ensemble. Mean end-of-century maize productivity is shifted from +5% to -6% (SSP126) and from +1% to -24% (SSP585)-explained by warmer climate projections and improved crop model sensitivities. In contrast, wheat shows stronger gains (+9% shifted to +18%, SSP585), linked to higher CO2 concentrations and expanded high-latitude gains. The 'emergence' of climate impacts consistently occurs earlier in the new projections-before 2040 for several main producing regions. While future yield estimates remain uncertain, these results suggest that major breadbasket regions will face distinct anthropogenic climatic risks sooner than previously anticipated.

The Global Gridded Crop Model Intercomparison phase 1 simulation dataset.

The Global Gridded Crop Model Intercomparison (GGCMI) phase 1 dataset of the Agricultural Model Intercomparison and Improvement Project (AgMIP) provides an unprecedentedly large dataset of crop model simulations covering the global ice-free land surface. The dataset consists of annual data fields at a spatial resolution of 0.5 arc-degree longitude and latitude. Fourteen crop modeling groups provided output for up to 11 historical input datasets spanning 1901 to 2012, and for up to three different management harmonization levels. Each group submitted data for up to 15 different crops and for up to 14 output variables. All simulations were conducted for purely rainfed and near-perfectly irrigated conditions on all land areas irrespective of whether the crop or irrigation system is currently used there. With the publication of the GGCMI phase 1 dataset we aim to promote further analyses and understanding of crop model performance, potential relationships between productivity and environmental impacts, and insights on how to further improve global gridded crop model frameworks. We describe dataset characteristics and individual model setup narratives.

Parameterization-induced uncertainties and impacts of crop management harmonization in a global gridded crop model ensemble.

Global gridded crop models (GGCMs) combine agronomic or plant growth models with gridded spatial input data to estimate spatially explicit crop yields and agricultural externalities at the global scale. Differences in GGCM outputs arise from the use of different biophysical models, setups, and input data. GGCM ensembles are frequently employed to bracket uncertainties in impact studies without investigating the causes of divergence in outputs. This study explores differences in maize yield estimates from five GGCMs based on the public domain field-scale model Environmental Policy Integrated Climate (EPIC) that participate in the AgMIP Global Gridded Crop Model Intercomparison initiative. Albeit using the same crop model, the GGCMs differ in model version, input data, management assumptions, parameterization, and selection of subroutines affecting crop yield estimates via cultivar distributions, soil attributes, and hydrology among others. The analyses reveal inter-annual yield variability and absolute yield levels in the EPIC-based GGCMs to be highly sensitive to soil parameterization and crop management. All GGCMs show an intermediate performance in reproducing reported yields with a higher skill if a static soil profile is assumed or sufficient plant nutrients are supplied. An in-depth comparison of setup domains for two EPIC-based GGCMs shows that GGCM performance and plant stress responses depend substantially on soil parameters and soil process parameterization, i.e. hydrology and nutrient turnover, indicating that these often neglected domains deserve more scrutiny. For agricultural impact assessments, employing a GGCM ensemble with its widely varying assumptions in setups appears the best solution for coping with uncertainties from lack of comprehensive global data on crop management, cultivar distributions and coefficients for agro-environmental processes. However, the underlying assumptions require systematic specifications to cover representative agricultural systems and environmental conditions. Furthermore, the interlinkage of parameter sensitivity from various domains such as soil parameters, nutrient turnover coefficients, and cultivar specifications highlights that global sensitivity analyses and calibration need to be performed in an integrated manner to avoid bias resulting from disregarded core model domains. Finally, relating evaluations of the EPIC-based GGCMs to a wider ensemble based on individual core models shows that structural differences outweigh in general differences in configurations of GGCMs based on the same model, and that the ensemble mean gains higher skill from the inclusion of structurally different GGCMs. Although the members of the wider ensemble herein do not consider crop-soil-management interactions, their sensitivity to nutrient supply indicates that findings for the EPIC-based sub-ensemble will likely become relevant for other GGCMs with the progressing inclusion of such processes.

Impacts and Uncertainties of +2°C of Climate Change and Soil Degradation on European Crop Calorie Supply.

Even if global warming is kept below +2°C, European agriculture will be significantly impacted. Soil degradation may amplify these impacts substantially and thus hamper crop production further. We quantify biophysical consequences and bracket uncertainty of +2°C warming on calories supply from 10 major crops and vulnerability to soil degradation in Europe using crop modeling. The Environmental Policy Integrated Climate (EPIC) model together with regional climate projections from the European branch of the Coordinated Regional Downscaling Experiment (EURO-CORDEX) was used for this purpose. A robustly positive calorie yield change was estimated for the EU Member States except for some regions in Southern and South-Eastern Europe. The mean impacts range from +30 Gcal ha-1 in the north, through +25 and +20 Gcal ha-1 in Western and Eastern Europe, respectively, to +10 Gcal ha-1 in the south if soil degradation and heat impacts are not accounted for. Elevated CO2 and increased temperature are the dominant drivers of the simulated yield changes in high-input agricultural systems. The growth stimulus due to elevated CO2 may offset potentially negative yield impacts of temperature increase by +2°C in most of Europe. Soil degradation causes a calorie vulnerability ranging from 0 to 50 Gcal ha-1 due to insufficient compensation for nutrient depletion and this might undermine climate benefits in many regions, if not prevented by adaptation measures, especially in Eastern and North-Eastern Europe. Uncertainties due to future potentials for crop intensification are about 2-50 times higher than climate change impacts.

Analyzing and modelling the effect of long-term fertilizer management on crop yield and soil organic carbon in China.

This study analyzes the influence of various fertilizer management practices on crop yield and soil organic carbon (SOC) based on the long-term field observations and modelling. Data covering 11 years from 8 long-term field trials were included, representing a range of typical soil, climate, and agro-ecosystems in China. The process-based model EPIC (Environmental Policy Integrated Climate model) was used to simulate the response of crop yield and SOC to various fertilization regimes. The results showed that the yield and SOC under additional manure application treatment were the highest while the yield under control treatment was the lowest (30%-50% of NPK yield) at all sites. The SOC in northern sites appeared more dynamic than that in southern sites. The variance partitioning analysis (VPA) showed more variance of crop yield could be explained by the fertilization factor (42%), including synthetic nitrogen (N), phosphorus (P), potassium (K) fertilizers, and fertilizer NPK combined with manure. The interactive influence of soil (total N, P, K, and available N, P, K) and climate factors (mean annual temperature and precipitation) determine the largest part of the SOC variance (32%). EPIC performs well in simulating both the dynamics of crop yield (NRMSE = 32% and 31% for yield calibration and validation) and SOC (NRMSE = 13% and 19% for SOC calibration and validation) under diverse fertilization practices in China. EPIC can assist in predicting the impacts of different fertilization regimes on crop growth and soil carbon dynamics, and contribute to the optimization of fertilizer management for different areas in China.

Uncertainty in soil data can outweigh climate impact signals in global crop yield simulations.

Global gridded crop models (GGCMs) are increasingly used for agro-environmental assessments and estimates of climate change impacts on food production. Recently, the influence of climate data and weather variability on GGCM outcomes has come under detailed scrutiny, unlike the influence of soil data. Here we compare yield variability caused by the soil type selected for GGCM simulations to weather-induced yield variability. Without fertilizer application, soil-type-related yield variability generally outweighs the simulated inter-annual variability in yield due to weather. Increasing applications of fertilizer and irrigation reduce this variability until it is practically negligible. Importantly, estimated climate change effects on yield can be either negative or positive depending on the chosen soil type. Soils thus have the capacity to either buffer or amplify these impacts. Our findings call for improvements in soil data available for crop modelling and more explicit accounting for soil variability in GGCM simulations.