Calibration and bias correction of climate projections for crop modelling: an idealised case study over Europe

Producing projections of future crop yields requires careful thought about the appropriate use of atmosphere-ocean global climate model (AOGCM) simulations. Here we describe and demonstrate multiple methods for ‘calibrating’ climate projections using an ensemble of AOGCM simulations in a ‘perfect sibling’ framework. Crucially, this type of analysis assesses the ability of each calibration methodology to produce reliable estimates of future climate, which is not possible just using historical observations. This type of approach could be more widely adopted for assessing calibration methodologies for crop modelling. The calibration methods assessed include the commonly used ‘delta’ (change factor) and ‘nudging’ (bias correction) approaches. We focus on daily maximum temperature in summer over Europe for this idealised case study, but the methods can be generalised to other variables and other regions. The calibration methods, which are relatively easy to implement given appropriate observations, produce more robust projections of future daily maximum temperatures and heat stress than using raw model output. The choice over which calibration method to use will likely depend on the situation, but change factor approaches tend to perform best in our examples. Finally, we demonstrate that the uncertainty due to the choice of calibration methodology is a significant contributor to the total uncertainty in future climate projections for impact studies. We conclude that utilising a variety of calibration methods on output from a wide range of AOGCMs is essential to produce climate data that will ensure robust and reliable crop yield projections.

Trait physiology and crop modelling as a framework to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding for complex traits because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. In this study, we explored whether physiological dissection and integrative modelling of complex traits could link phenotype complexity to underlying genetic systems in a way that enhanced the power of molecular breeding strategies. A crop and breeding system simulation study on sorghum, which involved variation in 4 key adaptive traits-phenology, osmotic adjustment, transpiration efficiency, stay-green-and a broad range of production environments in north-eastern Australia, was used. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages assuming gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies in the data. Based on the analyses of gene effects, a range of marker-assisted selection breeding strategies was simulated. It was shown that the inclusion of knowledge resulting from trait physiology and modelling generated an enhanced rate of yield advance over cycles of selection. This occurred because the knowledge associated with component trait physiology and extrapolation to the target population of environments by modelling removed confounding effects associated with environment and gene context dependencies for the markers used. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate genetic regions.

Trait physiology and crop modelling to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. We explore whether a crop growth and development modelling framework can link phenotype complexity to underlying genetic systems in a way that strengthens molecular breeding strategies. We use gene-to-phenotype simulation studies on sorghum to consider the value to marker-assisted selection of intrinsically stable QTLs that might be generated by physiological dissection of complex traits. The consequences on grain yield of genetic variation in four key adaptive traits – phenology, osmotic adjustment, transpiration efficiency, and staygreen – were simulated for a diverse set of environments by placing the known extent of genetic variation in the context of the physiological determinants framework of a crop growth and development model. It was assumed that the three to five genes associated with each trait, had two alleles per locus acting in an additive manner. The effects on average simulated yield, generated by differing combinations of positive alleles for the traits incorporated, varied with environment type. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages with gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies. We simulated a marker-assisted selection (MAS) breeding strategy based on the analyses of gene effects. When marker scores were allocated based on the contribution of gene effects to yield in a single environment, there was a wide divergence in rate of yield gain over all environments with breeding cycle depending on the environment chosen for the QTL analysis. It was suggested that knowledge resulting from trait physiology and modelling would overcome this dependency by identifying stable QTLs. The improved predictive power would increase the utility of the QTLs in MAS. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate QTLs.

Adapting wheat sowing dates to projected climate change in the Australian subtropics: analysis of crop water use and yield

Projected increases in atmospheric carbon dioxide concentration ([CO2]) and air temperature associated with future climate change are expected to affect crop development, crop yield, and, consequently, global food supplies. They are also likely to change agricultural production practices, especially those related to agricultural water management and sowing date. The magnitude of these changes and their implications to local production systems are mostly unknown. The objectives of this study were to: (i) simulate the effect of projected climate change on spring wheat (Triticum aestivum L. cv. Lang) yield and water use for the subtropical environment of the Darling Downs, Queensland, Australia; and (ii) investigate the impact of changing sowing date, as an adaptation strategy to future climate change scenarios, on wheat yield and water use. The multimodel climate projections from the IPCC Coupled Model Intercomparison Project (CMIP3) for the period 2030–2070 were used in this study. Climate scenarios included combinations of four changes in air temperature (08C, 18C, 28C, and 38C), three [CO2] levels (380 ppm, 500 ppm, and 600 ppm), and three changes in rainfall (–30%, 0%, and +20%), which were superimposed on observed station data. Crop management scenarios included a combination of six sowing dates (1 May, 10 May, 20 May, 1 June, 10 June, and 20 June) and three irrigation regimes (no irrigation (NI), deficit irrigation (DI), and full irrigation (FI)). Simulations were performed with the model DSSAT4.5, using 50 years of daily weather data.Wefound that: (1) grain yield and water-use efficiency (yield/evapotranspiration) increased linearly with [CO2]; (2) increases in [CO2] had minimal impact on evapotranspiration; (3) yield increased with increasing temperature for the irrigated scenarios (DI and FI), but decreased for the NI scenario; (4) yield increased with earlier sowing dates; and (5) changes in rainfall had a small impact on yield for DI and FI, but a high impact for the NI scenario.

Integrating crop improvement technologies for rapid genetic advance

This project developed a novel approach to integrating enhanced gene mapping technologies with crop modelling to enhance the rate of improvement in sorghum yield.

Modelling temperature, photoperiod and vernalization responses of Brunonia australis (Goodeniaceae) and Calandrinia sp (Portulacaceae) to predict flowering time

Crop models for herbaceous ornamental species typically include functions for temperature and photoperiod responses, but very few incorporate vernalization, which is a requirement of many traditional crops. This study investigated the development of floriculture crop models, which describe temperature responses, plus photoperiod or vernalization requirements, using Australian native ephemerals Brunonia australis and Calandrinia sp. A novel approach involved the use of a field crop modelling tool, DEVEL2. This optimization program estimates the parameters of selected functions within the development rate models using an iterative process that minimizes sum of squares residual between estimated and observed days for the phenological event. Parameter profiling and jack-knifing are included in DEVEL2 to remove bias from parameter estimates and introduce rigour into the parameter selection process. Development rate of B. australis from planting to first visible floral bud (VFB) was predicted using a multiplicative approach with a curvilinear function to describe temperature responses and a broken linear function to explain photoperiod responses. A similar model was used to describe the development rate of Calandrinia sp., except the photoperiod function was replaced with an exponential vernalization function, which explained a facultative cold requirement and included a coefficient for determining the vernalization ceiling temperature. Temperature was the main environmental factor influencing development rate for VFB to anthesis of both species and was predicted using a linear model. The phenology models for B. australis and Calandrinia sp. described development rate from planting to VFB and from VFB to anthesis in response to temperature and photoperiod or vernalization and may assist modelling efforts of other herbaceous ornamental plants. In addition to crop management, the vernalization function could be used to identify plant communities most at risk from predicted increases in temperature due to global warming.

Crop yield reduction in the tropics under climate change: Processes and uncertainties

Many modelling studies examine the impacts of climate change on crop yield, but few explore either the underlying bio-physical processes, or the uncertainty inherent in the parameterisation of crop growth and development. We used a perturbed-parameter crop modelling method together with a regional climate model (PRECIS) driven by the 2071-2100 SRES A2 emissions scenario in order to examine processes and uncertainties in yield simulation. Crop simulations used the groundnut (i.e. peanut; Arachis hypogaea L.) version of the General Large-Area Model for annual crops (GLAM). Two sets of GLAM simulations were carried out: control simulations and fixed-duration simulations, where the impact of mean temperature on crop development rate was removed. Model results were compared to sensitivity tests using two other crop models of differing levels of complexity: CROPGRO, and the groundnut model of Hammer et al. [Hammer, G.L., Sinclair, T.R., Boote, K.J., Wright, G.C., Meinke, H., and Bell, M.J., 1995, A peanut simulation model: I. Model development and testing. Agron. J. 87, 1085-1093]. GLAM simulations were particularly sensitive to two processes. First, elevated vapour pressure deficit (VPD) consistently reduced yield. The same result was seen in some simulations using both other crop models. Second, GLAM crop duration was longer, and yield greater, when the optimal temperature for the rate of development was exceeded. Yield increases were also seen in one other crop model. Overall, the models differed in their response to super-optimal temperatures, and that difference increased with mean temperature; percentage changes in yield between current and future climates were as diverse as -50% and over +30% for the same input data. The first process has been observed in many crop experiments, whilst the second has not. Thus, we conclude that there is a need for: (i) more process-based modelling studies of the impact of VPD on assimilation, and (ii) more experimental studies at super-optimal temperatures. Using the GLAM results, central values and uncertainty ranges were projected for mean 2071-2100 crop yields in India. In the fixed-duration simulations, ensemble mean yields mostly rose by 10-30%. The full ensemble range was greater than this mean change (20-60% over most of India). In the control simulations, yield stimulation by elevated CO2 was more than offset by other processes-principally accelerated crop development rates at elevated, but sub-optimal, mean temperatures. Hence, the quantification of uncertainty can facilitate relatively robust indications of the likely sign of crop yield changes in future climates. (C) 2007 Elsevier B.V. All rights reserved.

Use of a crop model ensemble to quantify CO2 stimulation of water-stressed and well-watered crops

Increased atmospheric concentrations of carbon dioxide (CO2) will benefit the yield of most crops. Two free air CO2 enrichment (FACE) meta-analyses have shown increases in yield of between 0 and 73% for C3 crops. Despite this large range, few crop modelling studies quantify the uncertainty inherent in the parameterisation of crop growth and development. We present a novel perturbed-parameter method of crop model simulation, which uses some constraints from observations, that does this. The model used is the groundnut (i.e. peanut; Arachis hypogaea L.) version of the general large-area model for annual crops (GLAM). The conclusions are of relevance to C3 crops in general. The increases in yield simulated by GLAM for doubled CO2 were between 16 and 62%. The difference in mean percentage increase between well-watered and water-stressed simulations was 6.8. These results were compared to FACE and controlled environment studies, and to sensitivity tests on two other crop models of differing levels of complexity: CROPGRO, and the groundnut model of Hammer et al. [Hammer, G.L., Sinclair, T.R., Boote, K.J., Wright, G.C., Meinke, H., Bell, M.J., 1995. A peanut simulation model. I. Model development and testing. Agron. J. 87, 1085-1093]. The relationship between CO2 and water stress in the experiments and in the models was examined. From a physiological perspective, water-stressed crops are expected to show greater CO2 stimulation than well-watered crops. This expectation has been cited in literature. However, this result is not seen consistently in either the FACE studies or in the crop models. In contrast, leaf-level models of assimilation do consistently show this result. An analysis of the evidence from these models and from the data suggests that scale (canopy versus leaf), model calibration, and model complexity are factors in determining the sign and magnitude of the interaction between CO2 and water stress. We conclude from our study that the statement that 'water-stressed crops show greater CO2 stimulation than well-watered crops' cannot be held to be universally true. We also conclude, preliminarily, that the relationship between water stress and assimilation varies with scale. Accordingly, we provide some suggestions on how studies of a similar nature, using crop models of a range of complexity, could contribute further to understanding the roles of model calibration, model complexity and scale. (C) 2008 Elsevier B.V. All rights reserved.

Measuring and predicting canopy nitrogen nutrition in wheat using a spectral index—The canopy chlorophyll content index (CCCI).

Varying the spatial distribution of applied nitrogen (N) fertilizer to match demand in crops has been shown to increase profits in Australia. Better matching the timing of N inputs to plant requirements has been shown to improve nitrogen use efficiency and crop yields and could reduce nitrous oxide emissions from broad acre grains. Farmers in the wheat production area of south eastern Australia are increasingly splitting N application with the second timing applied at stem elongation (Zadoks 30). Spectral indices have shown the ability to detect crop canopy N status but a robust method using a consistent calibration that functions across seasons has been lacking. One spectral index, the canopy chlorophyll content index (CCCI) designed to detect canopy N using three wavebands along the "red edge" of the spectrum was combined with the canopy nitrogen index (CNI), which was developed to normalize for crop biomass and correct for the N dilution effect of crop canopies. The CCCI-CNI index approach was applied to a 3-year study to develop a single calibration derived from a wheat crop sown in research plots near Horsham, Victoria, Australia. The index was able to predict canopy N (g m-2) from Zadoks 14-37 with an r2 of 0.97 and RMSE of 0.65 g N m-2 when dry weight biomass by area was also considered. We suggest that measures of N estimated from remote methods use N per unit area as the metric and that reference directly to canopy %N is not an appropriate method for estimating plant concentration without first accounting for the N dilution effect. This approach provides a link to crop development rather than creating a purely numerical relationship. The sole biophysical input, biomass, is challenging to quantify robustly via spectral methods. Combining remote sensing with crop modelling could provide a robust method for estimating biomass and therefore a method to estimate canopy N remotely. Future research will explore this and the use of active and passive sensor technologies for use in precision farming for targeted N management.

Economic case study of ABCD cane management practices in the Burdekin River Irrigation Area (BRIA) region

A case study was undertaken to determine the economic impact of a change in management class as detailed in the A, B, C and D management class framework. This document focuses on the implications of changing from D to C, C to B and B to A class management in the Burdekin River irrigation area (BRIA) and if the change is worthwhile from an economic perspective. This report provides a guide to the economic impact that may be expected when undertaking a particular change in farming practices and will ultimately lead to more informed decisions being made by key industry stakeholders. It is recognised that these management classes have certain limitations and in many cases the grouping of practices may not be reflective of the real situation. The economic case study is based on the A, B, C and D management class framework for water quality improvement developed in 2007/2008 for the Burdekin natural resource management region. The framework for the Burdekin is currently being updated to clarify some issues and incorporate new knowledge since the earlier version of the framework. However, this updated version is not yet complete and so the Paddock to Reef project has used the most current available version of the framework for the modelling and economics. As part of the project specification, sugarcane crop production data for the BRIA was provided by the APSIM model. The information obtained from the APSIM crop modelling programme included sugarcane yields and legume grain yield (legume grain yield only applies to A class management practice). Because of the complexity involved in the economic calculations, a combination of the FEAT, PiRisk and a custom made spreadsheet was used for the economic analysis. Figures calculated in the FEAT program were transferred to the custom made spreadsheet to develop a discounted cash flow analysis. The marginal cash flow differences for each farming system were simulated over a 5-year and 10-year planning horizon to determine the net present value of changing across different management practices. PiRisk was used to test uncertain parameters in the economic analysis and the potential risk associated with a change in value.

Economic case study of ABCD cane management practices in the Burdekin Delta region.

A case study was undertaken to determine the economic impact of a change in management class as detailed in the A, B, C and D management class framework. This document focuses on the implications of changing from D to C, C to B and B to A class management in the Burdekin Delta region and if the change is worthwhile from an economic perspective. This report provides a guide to the economic impact that may be expected when undertaking a particular change in farming practices and will ultimately lead to more informed decisions being made by key industry stakeholders. It is recognised that these management classes have certain limitations and in many cases the grouping of practices may not be reflective of the real situation. The economic case study is based on the A, B, C and D management class framework for water quality improvement developed in 2007/2008 for the Burdekin natural resource management region. The framework for the Burdekin is currently being updated to clarify some issues and incorporate new knowledge since the earlier version of the framework. However, this updated version is not yet complete and so the Paddock to Reef project has used the most current available version of the framework for the modelling and economics. As part of the project specification, sugarcane crop production data for the Burdekin Delta region was provided by the APSIM model. The information obtained from the APSIM crop modelling programme included sugarcane yields and legume grain yield (legume grain yield only applies to A class management practice). Because of the complexity involved in the economic calculations, a combination of the FEAT, PiRisk and a custom made spreadsheet was used for the economic analysis. Figures calculated in the FEAT program were transferred to the custom made spreadsheet to develop a discounted cash flow analysis. The marginal cash flow differences for each farming system were simulated over a 5-year and 10-year planning horizon to determine the Net Present Value of changing across different management practices. PiRisk was used to test uncertain parameters in the economic analysis and the potential risk associated with a change in value.

Processing Tomato Production in the Burdekin: Opportunities and Risks for Growers

The research undertaken here was in response to a decision by a major food producer in about 2009 to consider establishing processing tomato production in northern Australia. This was in response to a lack of water availability in the Goulburn Valley region following the extensive drought that continued until 2011. The high price of water and the uncertainty that went with it was important in making the decision to look at sites within Queensland. This presented an opportunity to develop a tomato production model for the varieties used in the processing industry and to use this as a case study along with rice and cotton production. Following some unsuccessful early trials and difficulties associated with the Global Financial Crisis, large scale studies by the food producer were abandoned. This report uses the data that was collected prior to this decision and contrasts the use of crop modelling with simpler climatic analyses that can be undertaken to investigate the impact of climate change on production systems. Crop modelling can make a significant contribution to our understanding of the impacts of climate variability and climate change because it harnesses the detailed understanding of physiology of the crop in a way that statistical or other analytical approaches cannot do. There is a high overhead, but given that trials are being conducted for a wide range of crops for a variety of purposes, breeding, fertiliser trials etc., it would appear to be profitable to link researchers with modelling expertise with those undertaking field trials. There are few more cost-effective approaches than modelling that can provide a pathway to understanding future climates and their impact on food production.

Maize yield and rainfall on different spatial and temporal scales in Southern Brazil

This study aimed to establish relationships between maize yield and rainfall on different temporal and spatial scales, in order to provide a basis for crop monitoring and modelling. A 16-year series of maize yield and daily rainfall from 11 municipalities and micro-regions of Rio Grande do Sul State was used. Correlation and regression analyses were used to determine associations between crop yield and rainfall for the entire crop cycle, from tasseling to 30 days after, and from 5 days before tasseling to 40 days after. Close relationships between maize yield and rainfall were found, particularly during the reproductive period (45-day period comprising the flowering and grain filling). Relationships were closer on a regional scale than at smaller scales. Implications of the crop-rainfall relationships for crop modelling are discussed.

Predicting the effects of nitrogen and planting density on maize water use in semi-arid Kenya

Models are important tools to assess the scope of management effects on crop productivity under different climatic and soil regimes. Accordingly, this study developed and used a simple model to assess the effects of nitrogen fertiliser and planting density on the water use efficiency (q) of maize in semi-arid Kenya. Field experiments were undertaken at Sonning, Berkshire, UK, in 1996 (one sowing) and 1997 (two sowings). The results from the field experiments plus soil and weather data for Machakos, Kenya (1 degree 33'S, 37 degree 14'E and 1560 m above sea level), were then used to predict the effects that N application and planting density may have on water use by a maize crop grown in semi-arid Kenya. The increase in q due to N application was greater under irrigated (15%-19%) than rainfed (7%-8%) conditions. Also, high planting density increased q (by 13%) under irrigation but decreased q (by 17%) under rainfed conditions. The current study has shown the significance of crop modelling techniques in assessing the influence of N and planting density on maize production in one region of semi-arid Kenya where there is high variability of rainfall.

Decimal growth stages for precision wheat production in changing environments?

The utility of the decimal growth stage (DGS) scoring system for cereals is reviewed. The DGS is the most widely used scale in academic and commercial applications because of its comprehensive coverage of cereal developmental stages, the ease of use and definition provided and adoption by official agencies. The DGS has demonstrable and established value in helping to optimise the timing of agronomic inputs, particularly with regard to plant growth regulators, herbicides, fungicides and soluble nitrogen fertilisers. In addition, the DGS is used to help parameterise crop models, and also in understanding the response and adaptation of crops to the environment. The value of the DGS for increasing precision relies on it indicating, to some degree, the various stages in the development of the stem apex and spike. Coincidence of specific growth stage scores with the transition of the apical meristem from a vegetative to a reproductive state, and also with the period of meiosis, is unreliable. Nonetheless, in pot experiments it is shown that the broad period of booting (DGS 41–49) appears adequate for covering the duration when the vulnerability of meiosis to drought and heat stress is exposed. Similarly, the duration of anthesis (61–69) is particularly susceptible to abiotic stresses: initially from a fertility perspective, but increasingly from a mean grain weight perspective as flowering progresses to DGS 69 and then milk development. These associations with DGS can have value at the crop level of organisation: for interpreting environmental effects, and in crop modelling. However, genetic, biochemical and physiological analysis to develop greater understanding of stress acclimation during the vegetative state, and tolerance at meiosis, does require more precision than DGS can provide. Similarly, individual floret analysis is needed to further understand the genetic basis of stress tolerance during anthesis.