Browsing by Author "Gaiser, Thomas"

Adaptation of crops to climate change has to be addressed locally due to the variability of soil, climate and the specific socio-economic settings influencing farm management decisions. Adaptation of rainfed cropping systems in the Mediterranean is especially challenging due to the projected decline in precipitation in the coming decades, which will increase the risk of droughts. Methods that can help explore uncertainties in climate projections and crop modelling, such as impact response surfaces (IRSs) and ensemble modelling, can then be valuable for identifying effective adaptations. Here, an ensemble of 17 crop models was used to simulate a total of 54 adaptation options for rainfed winter wheat (Triticum aestivum) at Lleida (NE Spain). To support the ensemble building, an ex post quality check of model simulations based on several criteria was performed. Those criteria were based on the “According to Our Current Knowledge” (AOCK) concept, which has been formalized here. Adaptations were based on changes in cultivars and management regarding phenology, vernalization, sowing date and irrigation. The effects of adaptation options under changed precipitation (P), temperature (T), [CO2] and soil type were analysed by constructing response surfaces, which we termed, in accordance with their specific purpose, adaptation response surfaces (ARSs). These were created to assess the effect of adaptations through a range of plausible P, T and [CO2] perturbations. The results indicated that impacts of altered climate were predominantly negative. No single adaptation was capable of overcoming the detrimental effect of the complex interactions imposed by the P, T and [CO2] perturbations except for supplementary irrigation (sI), which reduced the potential impacts under most of the perturbations. Yet, a combination of adaptations for dealing with climate change demonstrated that effective adaptation is possible at Lleida. Combinations based on a cultivar without vernalization requirements showed good and wide adaptation potential. Few combined adaptation options performed well under rainfed conditions. However, a single sI was sufficient to develop a high adaptation potential, including options mainly based on spring wheat, current cycle duration and early sowing date. Depending on local environment (e.g. soil type), many of these adaptations can maintain current yield levels under moderate changes in T and P, and some also under strong changes. We conclude that ARSs can offer a useful tool for supporting planning of field level adaptation under conditions of high uncertainty.

Modelling European Agriculture with Climate Change for Food Security (MACSUR) is a knowledge hub exploiting and improving data, methods and modelling tools for a detailed climate change risk assessment. The hub comprises 73 interacting agricultural (crop, livestock, trade) scientific and modelling research groups from 16 European countries and Israel. The crop modelling (CropM) component of MACSUR concentrates on overcoming weaknesses in crop modelling approaches and tools with specific attention to exploiting data on important European field crops, crop rotations and farming systems, and the modelling of diverse (mitigative) adaptation options. CropM outputs are scaled up to farm, regional and (supra-) national level as required for concerted integrated studies on the European agri-food sector and its contribution to global food security under climate change. The specific objectives of CropM are: (i) to conduct crop model intercomparisons to detect deficiencies, (ii) compile data in support of model improvements, (iii) advance scaling methods and model linkages, (iv) improve climate scenario data and impact uncertainty analysis, (v) build research capacity in these areas, and (vi) combine new knowledge and tools with those from livestock and trade modellers to allow interdisciplinary studies and interaction with a diverse range of stakeholders for climate change impact assessments. We identify requirements for improving model simulations, e.g. concerning impacts of heat and drought stress as well as intense rainfall and warm winters on crop yield. We show possibilities for enhancing methods of linking models and data of different resolutions. Finally, we give examples of how to improve quantification and reporting of crop impact uncertainties including the contribution from various sources (i.e. emission scenarios, climate modelling, downscaling of climate model data and crop impact modelling itself).

Crop growth simulation models are increasingly used for regionally assessing the effects of climate change and variability on crop yields. These models require spatially and temporally detailed, location-specific, environmental (weather and soil) and management data as inputs, which are often difficult to obtain consistently for larger regions. Aggregating the resolution of input data for crop model applications may increase the uncertainty of simulations to an extent that is not well understood. The present study aims to systematically analyse the effect of changes in the spatial resolution of weather input data on yields simulated by four crop models (LINTUL-SLIM, DSSAT-CSM, EPIC and WOFOST) which were utilized to test possible interactions between weather input data resolution and specific modelling approaches representing different degrees of complexity. The models were applied to simulate grain yield of spring barley in Finland for 12 years between 1994 and 2005 considering five spatial resolutions of daily weather data: weather station (point) and grid-based interpolated data at resolutions of 10 km × 10 km; 20 km × 20 km; 50 km × 50 km and 100 km × 100 km. Our results show that the differences between models were larger than the effect of the chosen spatial resolution of weather data for the considered years and region. When displaying model results graphically, each model exhibits a characteristic ‘fingerprint’ of simulated yield frequency distributions. These characteristic distributions in response to the inter-annual weather variability were independent of the spatial resolution of weather input data. Using one model (LINTUL-SLIM), we analysed how the aggregation strategy, i.e. aggregating model input versus model output data, influences the simulated yield frequency distribution. Results show that aggregating weather data has a smaller effect on the yield distribution than aggregating simulated yields which causes a deformation of the model fingerprint. We conclude that changes in the spatial resolution of weather input data introduce less uncertainty to the simulations than the use of different crop models but that more evaluation is required for other regions with a higher spatial heterogeneity in weather conditions, and for other input data related to soil and crop management to substantiate our findings. Our results provide further evidence to support other studies stressing the importance of using not just one, but different crop models in climate assessment studies.

Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (−2 to +9°C) and precipitation (−50 to +50%). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses. The model ensemble was used to simulate yields of winter and spring wheat at four sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern. The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description. Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index. Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.

Climate change impact assessments are plagued with uncertainties from many sources, such as climate projections or the inadequacies in structure and parameters of the impact model. Previous studies tried to account for the uncertainty from one or two of these. Here, we developed a triple‐ensemble probabilistic assessment using seven crop models, multiple sets of model parameters and eight contrasting climate projections together to comprehensively account for uncertainties from these three important sources. We demonstrated the approach in assessing climate change impact on barley growth and yield at Jokioinen, Finland in the Boreal climatic zone and Lleida, Spain in the Mediterranean climatic zone, for the 2050s. We further quantified and compared the contribution of crop model structure, crop model parameters and climate projections to the total variance of ensemble output using Analysis of Variance (ANOVA). Based on the triple‐ensemble probabilistic assessment, the median of simulated yield change was −4% and +16%, and the probability of decreasing yield was 63% and 31% in the 2050s, at Jokioinen and Lleida, respectively, relative to 1981–2010. The contribution of crop model structure to the total variance of ensemble output was larger than that from downscaled climate projections and model parameters. The relative contribution of crop model parameters and downscaled climate projections to the total variance of ensemble output varied greatly among the seven crop models and between the two sites. The contribution of downscaled climate projections was on average larger than that of crop model parameters. This information on the uncertainty from different sources can be quite useful for model users to decide where to put the most effort when preparing or choosing models or parameters for impact analyses. We concluded that the triple‐ensemble probabilistic approach that accounts for the uncertainties from multiple important sources provide more comprehensive information for quantifying uncertainties in climate change impact assessments as compared to the conventional approaches that are deterministic or only account for the uncertainties from one or two of the uncertainty sources.

The complexity of risks posed by climate change and possible adaptations for crop production has called for integrated assessment and modelling (IAM) approaches linking biophysical and economic models. This paper attempts to provide an overview of the present state of crop modelling to assess climate change risks to food production and to which extent crop models comply with IAM demands. Considerable progress has been made in modelling effects of climate variables, where crop models best satisfy IAM demands. Demands are partly satisfied for simulating commonly required assessment variables. However, progress on the number of simulated crops, uncertainty propagation related to model parameters and structure, adaptations and scaling are less advanced and lagging behind IAM demands. The limitations are considered substantial and apply to a different extent to all crop models. Overcoming these limitations will require joint efforts, and consideration of novel modelling approaches.

Diversification of crop rotations is considered an option to increase the resilience of European crop production under climate change. So far, however, many crop simulation studies have focused on predicting single crops in separate one-year simulations. Here, we compared the capability of fifteen crop growth simulation models to predict yields in crop rotations at five sites across Europe under minimal calibration. Crop rotations encompassed 301 seasons of ten crop types common to European agriculture and a diverse set of treatments (irrigation, fertilisation, CO2 concentration, soil types, tillage, residues, intermediate or catch crops). We found that the continuous simulation of multi-year crop rotations yielded results of slightly higher quality compared to the simulation of single years and single crops. Intermediate crops (oilseed radish and grass vegetation) were simulated less accurately than main crops (cereals). The majority of models performed better for the treatments of increased CO2 and nitrogen fertilisation than for irrigation and soil-related treatments. The yield simulation of the multi-model ensemble reduced the error compared to single-model simulations. The low degree of superiority of continuous simulations over single year simulation was caused by (a) insufficiently parameterised crops, which affect the performance of the following crop, and (b) the lack of growth-limiting water and/or nitrogen in the crop rotations under investigation. In order to achieve a sound representation of crop rotations, further research is required to synthesise existing knowledge of the physiology of intermediate crops and of carry-over effects from the preceding to the following crop, and to implement/improve the modelling of processes that condition these effects.

Climate change and its associated higher frequency and severity of adverse weather events require genotypic adaptation. Process-based ecophysiological modelling offers a powerful means to better target and accelerate development of new crop cultivars. Barley (Hordeum vulgare L.) is an important crop throughout the world, and a good model for study of the genetics of stress adaptation because many quantitative trait loci and candidate genes for biotic and abiotic stress tolerance have been identified in it. Here, we developed a new approach to design future crop ideotypes using an ensemble of eight barley simulation models (i.e. APSIM, CropSyst, HERMES, MCWLA, MONICA, SIMPLACE, SiriusQuality, and WOFOST), and applied it to design climate-resilient barley ideotypes for Boreal and Mediterranean climatic zones in Europe. The results showed that specific barley genotypes, represented by sets of cultivar parameters in the crop models, could be promising under future climate change conditions, resulting in increased yields and low inter-annual yield variability. In contrast, other genotypes could result in substantial yield declines. The most favorable climate-zone-specific barley ideotypes were further proposed, having combinations of several key genetic traits in terms of phenology, leaf growth, photosynthesis, drought tolerance, and grain formation. For both Boreal and Mediterranean climatic zones, barley ideotypes under future climatic conditions should have a longer reproductive growing period, lower leaf senescence rate, larger radiation use efficiency or maximum assimilation rate, and higher drought tolerance. Such characteristics can produce substantial positive impacts on yields under contrasting conditions. Moreover, barley ideotypes should have a low photoperiod and high vernalization sensitivity for the Boreal climatic zone; for the Mediterranean, in contrast, it should have a low photoperiod and low vernalization sensitivity. The drought-tolerance trait is more beneficial for the Mediterranean than for the Boreal climatic zone. Our study demonstrates a sound approach to design future barley ideotypes based on an ensemble of well-tested, diverse crop models and on integration of knowledge from multiple disciplines. The robustness of model-aided ideotypes design can be further enhanced by continuously improving crop models and enhancing information exchange between modellers, agro-meteorologists, geneticists, and plant breeders.

Understanding the drivers of yield levels under climate change is required to support adaptation planning and respond to changing production risks. This study uses an ensemble of crop models applied on a spatial grid to quantify the contributions of various climatic drivers to past yield variability in grain maize and winter wheat of European cropping systems (1984–2009) and drivers of climate change impacts to 2050. Results reveal that for the current genotypes and mix of irrigated and rainfed production, climate change would lead to yield losses for grain maize and gains for winter wheat. Across Europe, on average heat stress does not increase for either crop in rainfed systems, while drought stress intensifies for maize only. In low-yielding years, drought stress persists as the main driver of losses for both crops, with elevated CO2 offering no yield benefit in these years.

Early vigour in wheat is a trait that has received attention for its benefits reducing evaporation from the soil surface early in the season. However, with the growth enhancement common to crops grown under elevated atmospheric CO2 concentrations (e[CO2]), there is a risk that too much early growth might deplete soil water and lead to more severe terminal drought stress in environments where production relies on stored soil water content. If this is the case, the incorporation of such a trait in wheat breeding programmes might have unintended negative consequences in the future, especially in dry years. We used selected data from cultivars with proven expression of high and low early vigour from the Australian Grains Free Air CO2 Enrichment (AGFACE) facility, and complemented this analysis with simulation results from two crop growth models which differ in the modelling of leaf area development and crop water use. Grain yield responses to e[CO2] were lower in the high early vigour group compared to the low early vigour group, and although these differences were not significant, they were corroborated by simulation model results. However, the simulated lower response with high early vigour lines was not caused by an earlier or greater depletion of soil water under e[CO2] and the mechanisms responsible appear to be related to an earlier saturation of the radiation intercepted. Whether this is the case in the field needs to be further investigated. In addition, there was some evidence that the timing of the drought stress during crop growth influenced the effect of e[CO2] regardless of the early vigour trait. There is a need for FACE investigations of the value of traits for drought adaptation to be conducted under more severe drought conditions and variable timing of drought stress, a risky but necessary endeavour.

We assessed the weather data aggregation effect (DAE) on the simulation of cropping systems for different crops, response variables, and production conditions. Using 13 process-based crop models and the ensemble mean, we simulated 30 yr continuous cropping systems for 2 crops (winter wheat and silage maize) under 3 production conditions for the state of North Rhine-Westphalia, Germany. The DAE was evaluated for 5 weather data resolutions (i.e. 1, 10, 25, 50, and 100 km) for 3 response variables including yield, growing season evapotranspiration, and water use efficiency. Five metrics, viz. the spatial bias (Δ), average absolute deviation (AAD), relative AAD, root mean squared error (RMSE), and relative RMSE, were used to evaluate the DAE on both the input weather data and simulated results. For weather data, we found that data aggregation narrowed the spatial variability but widened the Δ, especially across mountainous areas. The DAE on loss of spatial heterogeneity and hotspots was stronger than on the average changes over the region. The DAE increased when coarsening the spatial resolution of the input weather data. The DAE varied considerably across different models, but changed only slightly for different production conditions and crops. We conclude that if spatially detailed information is essential for local management decision, higher resolution is desirable to adequately capture the spatial variability for heterogeneous regions. The required resolution depends on the choice of the model as well as the environmental condition of the study area.

We compared the precision of simple random sampling (SimRS) and seven types of stratified random sampling (StrRS) schemes in estimating regional mean of water-limited yields for two crops (winter wheat and silage maize) that were simulated by fourteen crop models. We found that the precision gains of StrRS varied considerably across stratification methods and crop models. Precision gains for compact geographical stratification were positive, stable and consistent across crop models. Stratification with soil water holding capacity had very high precision gains for twelve models, but resulted in negative gains for two models. Increasing the sample size monotonously decreased the sampling errors for all the sampling schemes. We conclude that compact geographical stratification can modestly but consistently improve the precision in estimating regional mean yields. Using the most influential environmental variable for stratification can notably improve the sampling precision, especially when the sensitivity behavior of a crop model is known.

AbstractWheat is the most widely grown food crop, with 761 Mt produced globally in 2020. To meet the expected grain demand by mid-century, wheat breeding strategies must continue to improve upon yield-advancing physiological traits, regardless of climate change impacts. Here, the best performing doubled haploid (DH) crosses with an increased canopy photosynthesis from wheat field experiments in the literature were extrapolated to the global scale with a multi-model ensemble of process-based wheat crop models to estimate global wheat production. The DH field experiments were also used to determine a quantitative relationship between wheat production and solar radiation to estimate genetic yield potential. The multi-model ensemble projected a global annual wheat production of 1050 ± 145 Mt due to the improved canopy photosynthesis, a 37% increase, without expanding cropping area. Achieving this genetic yield potential would meet the lower estimate of the projected grain demand in 2050, albeit with considerable challenges.

The spatial variability of soil properties is an important driver of yield variability at both field and regional scale. Thus, when using crop growth simulation models, the choice of spatial resolution of soil input data might be key in order to accurately reproduce observed yield variability. In this study we used four crop models (SIMPLACE, DSSAT-CSM, EPIC and DAISY) differing in the detail of modeling above-ground biomass and yield as well as of modeling soil water dynamics, water uptake and drought effects on plants to simulate winter wheat in two (agro-climatologically and geo-morphologically) contrasting regions of the federal state of North-Rhine-Westphalia (Germany) for the period from 1995 to 2008. Three spatial resolutions of soil input data were taken into consideration, corresponding to the following map scales: 1:50 000, 1:300 000 and 1:1 000 000. The four crop models were run for water-limited production conditions and model results were evaluated in the form of frequency distributions, depicted by bean-plots. In both regions, soil data aggregation had very small influence on the shape and range of frequency distributions of simulated yield and simulated total growing season evapotranspiration for all models. Further analysis revealed that the small influence of spatial resolution of soil input data might be related to: (a) the high precipitation amount in the region which partly masked differences in soil characteristics for water holding capacity, (b) the loss of variability in hydraulic soil properties due to the methods applied to calculate water retention properties of the used soil profiles, and (c) the method of soil data aggregation. No characteristic “fingerprint” between sites, years and resolutions could be found for any of the models. Our results support earlier recommendation to evaluate model results on the basis of frequency distributions since these offer quick and better insight into the distribution of simulation results as compared to summary statistics only. Finally, our results support conclusions from other studies about the usefulness of considering a multi-model approach to quantify the uncertainty in simulated yields introduced by the crop growth simulation approach when exploring the effects of scaling for regional yield impact assessments.