Improved descriptions of soil hydrology in crop models: The elephant in the room?

Soil-crop simulation models are widely used to assess the impacts of soil management and climate change on soil water balance, solute transport and crop production. In this context, it is important that hydrological processes in the soil-crop system are accurately modelled. We suggest here that empirical treatments of soil water flow, water uptake by plant mots and transpiration limit the applicability of crop models and increase prediction errors. We further argue that this empiricism is to a large extent unnecessary, as parsimonious physics-based descriptions of these water flow processes in the soil-crop system are now available. Recent reviews and opinion articles, whilst strongly advocating the need for improvements to crop models, fail to mention the significant role played by accurate treatments of soil hydrology. It seems to us that empirical models of soil water flow have become the elephant in the room

Impact of spatial soil and climate input data aggregation on regional yield simulations

We show the error in water-limited yields simulated by crop models which is associated with spatially aggregated soil and climate input data. Crop simulations at large scales (regional, national, continental) frequently use input data of low resolution. Therefore, climate and soil data are often generated via averaging and sampling by area majority. This may bias simulated yields at large scales, varying largely across models. Thus, we evaluated the error associated with spatially aggregated soil and climate data for 14 crop models. Yields of winter wheat and silage maize were simulated under water-limited production conditions. We calculated this error from crop yields simulated at spatial resolutions from 1 to 100 km for the state of North Rhine-Westphalia, Germany. Most models showed yields biased by <15% when aggregating only soil data. The relative mean absolute error (rMAE) of most models using aggregated soil data was in the range or larger than the inter-annual or inter-model variability in yields. This error increased further when both climate and soil data were aggregated. Distinct error patterns indicate that the rMAE may be estimated from few soil variables. Illustrating the range of these aggregation effects across models, this study is a first step towards an ex-ante assessment of aggregation errors in large-scale simulations

Aggregation of soil and climate input data can underestimate simulated biomass loss and nitrate leaching under climate change

Predicting areas of severe biomass loss and increased N leaching risk under climate change is critical for applying appropriate adaptation measures to support more sustainable agricultural systems. The frequency of annual severe biomass loss for winter wheat and its coincidence with an increase in N leaching in a temperate region in Germany was estimated including the error from using soil and climate input data at coarser spatial scales, using the soil-crop model CoupModel. We ran the model for a reference period (1980-2010) and used climate data predicted by four climate model(s) for the Representative Concentration Pathways (RCP) 2.6, 4.5 and 8.5. The annual median biomass estimations showed that for the period 2070-2100, under the RCP8.5 scenario, the entire region would suffer from severe biomass loss almost every year. Annual incidence of severe biomass loss and increased N leaching was predicted to increase from RCP4.5 to the 8.5 scenario. During 2070-2100 for RCP8.5, in more than half of the years an area of 95% of the region was projected to suffer from both severe biomass loss and increased N leaching. The SPEI3 predicted a range of 32 (P3 RCP4.5) to 55% (P3 RCP8.5) of the severe biomass loss episodes simulated in the climate change scenarios. The simulations predicted more severe biomass losses than by the SPEI index which indicates that soil water deficits are important in determining crop losses in future climate scenarios. There was a risk of overestimating the area where "no severe biomass loss + increased N leaching" occurred when using coarser aggregated input data. In contrast, underestimation of situations where "severe biomass loss + increased N leaching" occurred when using coarser aggregated input data. Larger annual differences in biomass estimations compared to the finest resolution of input data occurred when aggregating climate input data rather than soil data. The differences were even larger when aggregating both soil and climate input data. In half of the region, biomass could be erroneously estimated in a single year by more than 40% if using soil and climate coarser input data. The results suggest that a higher spatial resolution of especially climate input data would be needed to predict reliably annual estimates of severe biomass loss and N leaching under climate change scenarios

Effects of input data aggregation on simulated crop yields in temperate and Mediterranean climates

The modelling exercise for this study was highly supported by partner universities and research institutes in the framework of the MACSUR project and financially supported by the German Federal Ministry of Education and Research BMBF (FKZ 2815ERA01J) in the framework of the funding measure “Soil as a Sustainable Resource for the Bioeconomy – BonaRes”, project “BonaRes (Module B): BonaRes Centre for Soil Research (FKZ BOMA03037514, 031B0026A and 031A608A) and by the Ministry of Agriculture and Food (BMEL) in the framework of the MACSUR project (FKZ 2815ERA01J). In addition, the relevant co-authors from the partner institutes are separately financed by their respective projects. AV, EC, and EL were supported by The Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning (220-2007-1218) and by the strategic funding ‘Soil-Water-Landscape’ from the faculty of Natural Resources and Agricultural Sciences (Swedish University of Agricultural Sciences). JC thank the INRA ACCAF metaprogramm for funding. KCK, CN, XS and TS were supported by MACSUR2 (FKZ 031B0039C). MK thanks for the funding by the UK BBSRC (BB/N004922/1) and the MAXWELL HPC team of the University of Aberdeen for providing equipment and support for the DailyDayCent simulations. FE acknowledges support by the German Science Foundation (project EW 119/5-1). GRM, TG, and FE thank Andreas Enders and Gunther Krauss (INRES, University of Bonn) for support. The authors also would like to acknowledge the support provided by the BMBF and the valuable comments of the scientists of the Institut für Nutzpflanzenwissenschaften und Ressourcenschutz (INRES), University of Bonn, Germany.Peer reviewedPostprin

Impact of Spatial Soil and Climate Input Data Aggregation on Regional Yield Simulations

This work was financially supported by the German Federal Ministry of Food and Agriculture (BMEL) through the Federal Office for Agriculture and Food (BLE), (2851ERA01J). FT and RPR were supported by FACCE MACSUR (3200009600) through the Finnish Ministry of Agriculture and Forestry (MMM). EC, HE and EL were supported by The Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning (220-2007-1218) and by the strategic funding ‘Soil-Water-Landscape’ from the faculty of Natural Resources and Agricultural Sciences (Swedish University of Agricultural Sciences) and thank professor P-E Jansson (Royal Institute of Technology, Stockholm) for support. JC, HR and DW thank the INRA ACCAF metaprogramm for funding and Eric Casellas from UR MIAT INRA for support. CB was funded by the Helmholtz project “REKLIM—Regional Climate Change”. CK was funded by the HGF Alliance “Remote Sensing and Earth System Dynamics” (EDA). FH was funded by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG) under the Grant FOR1695. FE and SS acknowledge support by the German Science Foundation (project EW 119/5-1). HH, GZ, SS, TG and FE thank Andreas Enders and Gunther Krauss (INRES, University of Bonn) for support. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Peer reviewedPublisher PD

A framework for modelling soil structure dynamics induced by biological activity

Acknowledgments: This work was funded by the Swedish Research Council for Sustainable Development (FORMAS) in the project “Soil structure and soil degradation: improved model tools to meet sustainable development goals under climate and land use change” (grant no. 2018-02319). We would also like to thank Mikael Sasha Dooha for carrying out the measurements for the water retention curves shown in figure 4.Peer reviewedPublisher PD

Impact analysis of climate data aggregation at different spatial scales on simulated net primary productivity for croplands

For spatial crop and agro-systems modelling, there is often a discrepancy between the scale of measured driving data and the target resolution. Spatial data aggregation is often necessary, which can introduce additional uncertainty into the simulation results. Previous studies have shown that climate data aggregation has little effect on simulation of phenological stages, but effects on net primary production (NPP) might still be expected through changing the length of the growing season and the period of grain filling. This study investigates the impact of spatial climate data aggregation on NPP simulation results, applying eleven different models for the same study region (∼34,000 km2), situated in Western Germany. To isolate effects of climate, soil data and management were assumed to be constant over the entire study area and over the entire study period of 29 years. Two crops, winter wheat and silage maize, were tested as monocultures. Compared to the impact of climate data aggregation on yield, the effect on NPP is in a similar range, but is slightly lower, with only small impacts on averages over the entire simulation period and study region. Maximum differences between the five scales in the range of 1–100 km grid cells show changes of 0.4–7.8% and 0.0–4.8% for wheat and maize, respectively, whereas the simulated potential NPP averages of the models show a wide range (1.9–4.2 g C m−2 d−1 and 2.7–6.1 g C m−2 d−1 for wheat and maize, respectively). The impact of the spatial aggregation was also tested for shorter time periods, to see if impacts over shorter periods attenuate over longer periods. The results show larger impacts for single years (up to 9.4% for wheat and up to 13.6% for maize). An analysis of extreme weather conditions shows an aggregation effect in vulnerability up to 12.8% and 15.5% between the different resolutions for wheat and maize, respectively. Simulations of NPP averages over larger areas (e.g. regional scale) and longer time periods (several years) are relatively insensitive to climate data aggregation. However, the scale of climate data is more relevant for impacts on annual averages of NPP or if the period is strongly affected or dominated by drought stress. There should be an awareness of the greater uncertainty for the NPP values in these situations if data are not available at high resolution. On the other hand, the results suggest that there is no need to simulate at high resolution for long term regional NPP averages based on the simplified assumptions (soil and management constant in time and space) used in this study

Effects of climate input data aggregation on modelling regional crop yields

Crop models can be sensitive to climate input data aggregation and this response may differ among models. This should be considered when applying field-scale models for assessment of climate change impacts on larger spatial scales or when coupling models across scales. In order to evaluate these effects systematically, an ensemble of ten crop models was run with climate input data on different spatial aggregations ranging from 1, 10, 25, 50 and 100 km horizontal resolution for the state of North Rhine-Westphalia, Germany. Models were minimally calibrated to typical sowing and harvest dates, and crop yields observed in the region, subsequently simulating potential, water-limited and nitrogen-limited production of winter wheat and silage maize for 1982-2011. Outputs were analysed for 19 variables (yield, evapotranspiration, soil organic carbon, etc.). In this study the sensitivity of the individual models and the model ensemble in response to input data aggregation is assessed for crop yield. Results show that the mean yield of the region calculated from climate time series of 1 km horizontal resolution changes only little when using climate input data of higher aggregation levels for most models. However, yield frequency distributions change with aggregation, resembling observed data better with increasing resolution. With few exceptions, these results apply to the two crops and three production situations (potential, water-, nitrogen-limited) and across models including the model ensemble, regardless of differences among models in simulated yield levels and spatial yield patterns. Results of this study improve the confidence of using crop models at varying scales

The chaos in calibrating crop models

Calibration, the estimation of model parameters based on fitting the model to experimental data, is among the first steps in many applications of system models and has an important impact on simulated values. Here we propose and illustrate a novel method of developing guidelines for calibration of system models. Our example is calibration of the phenology component of crop models. The approach is based on a multi-model study, where all teams are provided with the same data and asked to return simulations for the same conditions. All teams are asked to document in detail their calibration approach, including choices with respect to criteria for best parameters, choice of parameters to estimate and software. Based on an analysis of the advantages and disadvantages of the various choices, we propose calibration recommendations that cover a comprehensive list of decisions and that are based on actual practices.HighlightsWe propose a new approach to deriving calibration recommendations for system modelsApproach is based on analyzing calibration in multi-model simulation exercisesResulting recommendations are holistic and anchored in actual practiceWe apply the approach to calibration of crop models used to simulate phenologyRecommendations concern: objective function, parameters to estimate, software usedCompeting Interest StatementThe authors have declared no competing interest