Ensemble yield simulations: crop and climate uncertainties, sensitivity to temperature and genotypic adaptation to climate change

Estimates of the response of crops to climate change rarely quantify the uncertainty inherent in the simulation of both climate and crops. We present a crop simulation ensemble for a location in India, perturbing the response of both crop and climate under both baseline (12 720 simulations) and doubled-CO2 (171720 simulations) climates. Some simulations used parameter values representing genotypic adaptation to mean temperature change. Firstly, observed and simulated yields in the baseline climate were compared. Secondly, the response of yield to changes in mean temperature was examined and compared to that found in the literature. No consistent response to temperature change was found across studies. Thirdly, the relative contribution of uncertainty in crop and climate simulation to the total uncertainty in projected yield changes was examined. In simulations without genotypic adaptation, most of the uncertainty came from the climate model parameters. Comparison with the simulations with genotypic adaptation and with a previous study suggested that the relatively low crop parameter uncertainty derives from the observational constraints on the crop parameters used in this study. Fourthly, the simulations were used, together with an observed dataset and a simple analysis of crop cardinal temperatures and thermal time, to estimate the potential for adaptation using existing cultivars. The results suggest that the germplasm for complete adaptation of groundnut cultivation in western India to a doubled-CO2 environment may not exist. In conjunction with analyses of germplasm and local management

Quantifying uncertainty in the simulation of crop-climate processes

The impacts of climate change on crop productivity are often assessed using simulations from a numerical climate model as an input to a crop simulation model. The precision of these predictions reflects the uncertainty in both models. We examined how uncertainty in a climate (HadAM3) and crop General Large-Area Model (GLAM) for annual crops model affects the mean and standard deviation of crop yield simulations in present and doubled carbon dioxide (CO2) climates by perturbation of parameters in each model. The climate sensitivity parameter (λ, the equilibrium response of global mean surface temperature to doubled CO2) was used to define the control climate. Observed 1966–1989 mean yields of groundnut (Arachis hypogaea L.) in India were simulated well by the crop model using the control climate and climates with values of λ near the control value. The simulations were used to measure the contribution to uncertainty of key crop and climate model parameters. The standard deviation of yield was more affected by perturbation of climate parameters than crop model parameters in both the present-day and doubled CO2 climates. Climate uncertainty was higher in the doubled CO2 climate than in the present-day climate. Crop transpiration efficiency was key to crop model uncertainty in both present-day and doubled CO2 climates. The response of crop development to mean temperature contributed little uncertainty in the present-day simulations but was among the largest contributors under doubled CO2. The ensemble methods used here to quantify physical and biological uncertainty offer a method to improve model estimates of the impacts of climate change.

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.

Adaptation of crops to climate change through genotypic responses to mean and extreme temperatures

The importance of temperature in the determination of the yield of an annual crop (groundnut; Arachis hypogaea L. in India) was assessed. Simulations from a regional climate model (PRECIS) were used with a crop model (GLAM) to examine crop growth under simulated current (1961-1990) and future (2071-2100) climates. Two processes were examined: the response of crop duration to mean temperature and the response of seed-set to extremes of temperature. The relative importance of, and interaction between, these two processes was examined for a number of genotypic characteristics, which were represented by using different values of crop model parameters derived from experiments. The impact of mean and extreme temperatures varied geographically, and depended upon the simulated genotypic properties. High temperature stress was not a major determinant of simulated yields in the current climate, but affected the mean and variability of yield under climate change in two regions which had contrasting statistics of daily maximum temperature. Changes in mean temperature had a similar impact on mean yield to that of high temperature stress in some locations and its effects were more widespread. Where the optimal temperature for development was exceeded, the resulting increase in duration in some simulations fully mitigated the negative impacts of extreme temperatures when sufficient water was available for the extended growing period. For some simulations the reduction in mean yield between the current and future climates was as large as 70%, indicating the importance of genotypic adaptation to changes in both means and extremes of temperature under climate change. (c) 2006 Elsevier B.V. All rights reserved.

Simulation of the impact of high temperature stress on annual crop yields

Brief periods of high temperature which occur near flowering can severely reduce the yield of annual crops such as wheat and groundnut. A parameterisation of this well-documented effect is presented for groundnut (i.e. peanut; Arachis hypogaeaL.). This parameterisation was combined with an existing crop model, allowing the impact of season-mean temperature, and of brief high-temperature episodes at various times near flowering, to be both independently and jointly examined. The extended crop model was tested with independent data from controlled environment experiments and field experiments. The impact of total crop duration was captured, with simulated duration being within 5% of observations for the range of season-mean temperatures used (20-28 degrees C). In simulations across nine differently timed high temperature events, eight of the absolute differences between observed and simulated yield were less than 10% of the control (no-stress) yield. The parameterisation of high temperature stress also allows the simulation of heat tolerance across different genotypes. Three parameter sets, representing tolerant, moderately sensitive and sensitive genotypes were developed and assessed. The new parameterisation can be used in climate change studies to estimate the impact of heat stress on yield. It can also be used to assess the potential for adaptation of cropping systems to increased temperature threshold exceedance via the choice of genotype characteristics. (c) 2005 Elsevier B.V. All rights reserved.

Design and optimisation of a large-area process-based model for annual crops

The formulation of a new process-based crop model, the general large-area model (GLAM) for annual crops is presented. The model has been designed to operate on spatial scales commensurate with those of global and regional climate models. It aims to simulate the impact of climate on crop yield. Procedures for model parameter determination and optimisation are described, and demonstrated for the prediction of groundnut (i.e. peanut; Arachis hypogaea L.) yields across India for the period 1966-1989. Optimal parameters (e.g. extinction coefficient, transpiration efficiency, rate of change of harvest index) were stable over space and time, provided the estimate of the yield technology trend was based on the full 24-year period. The model has two location-specific parameters, the planting date, and the yield gap parameter. The latter varies spatially and is determined by calibration. The optimal value varies slightly when different input data are used. The model was tested using a historical data set on a 2.5degrees x 2.5degrees grid to simulate yields. Three sites are examined in detail-grid cells from Gujarat in the west, Andhra Pradesh towards the south, and Uttar Pradesh in the north. Agreement between observed and modelled yield was variable, with correlation coefficients of 0.74, 0.42 and 0, respectively. Skill was highest where the climate signal was greatest, and correlations were comparable to or greater than correlations with seasonal mean rainfall. Yields from all 35 cells were aggregated to simulate all-India yield. The correlation coefficient between observed and simulated yields was 0.76, and the root mean square error was 8.4% of the mean yield. The model can be easily extended to any annual crop for the investigation of the impacts of climate variability (or change) on crop yield over large areas. (C) 2004 Elsevier B.V. All rights reserved.

Simulation of crop yields using ERA-40: Limits to skill and nonstationarity in weather-yield relationships

Reanalysis data provide an excellent test bed for impacts prediction systems. because they represent an upper limit on the skill of climate models. Indian groundnut (Arachis hypogaea L.) yields have been simulated using the General Large-Area Model (GLAM) for annual crops and the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-yr reanalysis (ERA-40). The ability of ERA-40 to represent the Indian summer monsoon has been examined. The ability of GLAM. when driven with daily ERA-40 data, to model both observed yields and observed relationships between subseasonal weather and yield has been assessed. Mean yields "were simulated well across much of India. Correlations between observed and modeled yields, where these are significant. are comparable to correlations between observed yields and ERA-40 rainfall. Uncertainties due to the input planting window, crop duration, and weather data have been examined. A reduction in the root-mean-square error of simulated yields was achieved by applying bias correction techniques to the precipitation. The stability of the relationship between weather and yield over time has been examined. Weather-yield correlations vary on decadal time scales. and this has direct implications for the accuracy of yield simulations. Analysis of the skewness of both detrended yields and precipitation suggest that nonclimatic factors are partly responsible for this nonstationarity. Evidence from other studies, including data on cereal and pulse yields, indicates that this result is not particular to groundnut yield. The detection and modeling of nonstationary weather-yield relationships emerges from this study as an important part of the process of understanding and predicting the impacts of climate variability and change on crop yields.

Quantification of physical and biological uncertainty in the simulation of the yield of a tropical crop using present-day and doubled CO2 climates

The impacts of climate change on crop productivity are often assessed using simulations from a numerical climate model as an input to a crop simulation model. The precision of these predictions reflects the uncertainty in both models. We examined how uncertainty in a climate (HadAM3) and crop General Large-Area Model (GLAM) for annual crops model affects the mean and standard deviation of crop yield simulations in present and doubled carbon dioxide (CO2) climates by perturbation of parameters in each model. The climate sensitivity parameter (lambda, the equilibrium response of global mean surface temperature to doubled CO2) was used to define the control climate. Observed 1966-1989 mean yields of groundnut (Arachis hypogaea L.) in India were simulated well by the crop model using the control climate and climates with values of lambda near the control value. The simulations were used to measure the contribution to uncertainty of key crop and climate model parameters. The standard deviation of yield was more affected by perturbation of climate parameters than crop model parameters in both the present-day and doubled CO2 climates. Climate uncertainty was higher in the doubled CO2 climate than in the present-day climate. Crop transpiration efficiency was key to crop model uncertainty in both present-day and doubled CO2 climates. The response of crop development to mean temperature contributed little uncertainty in the present-day simulations but was among the largest contributors under doubled CO2. The ensemble methods used here to quantify physical and biological uncertainty offer a method to improve model estimates of the impacts of climate change.

In vitro evaluation of fibrolytic enzymes as additives for maize (Zea mays L.) silage - I. Effects of ensiling temperature, enzyme source and addition level

A series of experiments was completed to investigate the impact of addition of enzymes at ensiling on in vitro rumen degradation of maize silage. Two commercial products, Depot 40 (D, Biocatalysts Ltd., Pontypridd, UK) and Liquicell 2500 (L, Specialty Enzymes and Biochemicals, Fresno, CA, USA), were used. In experiment 1, the pH optima over a pH range 4.0-6.8 and the stability of D and L under changing pH (4.0, 5.6, 6.8) and temperature (15 and 39 degreesC) conditions were determined. In experiment 2, D and L were applied at three levels to whole crop maize at ensiling, using triplicate 0.5 kg capacity laboratory minisilos. A completely randomized design with a factorial arrangement of treatments was used. One set of treatments was stored at room temperature, whereas another set was stored at 40 degreesC during the first 3 weeks of fermentation, and then stored at room temperature. Silages were opened after 120 days. Results from experiment I indicated that the xylanase activity of both products showed an optimal pH of about 5.6, but the response differed according to the enzyme, whereas the endoglucanase activity was inversely related to pH. Both products retained at least 70% of their xylanase activity after 48 h incubation at 15 or 39 degreesC. In experiment 2, enzymes reduced (P < 0.05) silage pH, regardless of storage temperature and enzyme level. Depol 40 reduced (P < 0.05) the starch contents of the silages, due to its high alpha-amylase activity. This effect was more noticeable in the silages stored at room temperature. Addition of L reduced (P < 0.05) neutral detergent fiber (NDF) and acid detergent fiber (ADF) contents. In vitro rumen degradation, assessed using the Reading Pressure Technique (RPT), showed that L increased (P < 0.05) the initial 6 h gas production (GP) and organic matter degradability (OMD), but did not affect (P > 0.05) the final extent of OMD, indicating that this preparation acted on the rumen degradable material. In contrast, silages treated with D had reduced (P < 0.05) rates of gas production and OMD. These enzymes, regardless of ensiling temperature, can be effective in improving the nutritive quality of maize silage when applied at ensiling. However, the biochemical properties of enzymes (i.e., enzymic activities, optimum pH) may have a crucial role in dictating the nature of the responses. (C) 2003 Elsevier B.V. All rights reserved.

Heat tolerance in groundnut

Tolerance to high soil and air temperature during the reproductive phase is an important component of adaptation to and and semi-arid cropping environments in groundnut. Between 10 and 22 genotypes were screened for tolerance to high air and soil temperature in controlled environments. To assess tolerance to high soil temperature, 10 genotypes were grown from start of podding to harvest at ambient (28 degrees) and high (38 degreesC) soil temperatures, and crop growth rate (CGR), pod growth rate (PGR) and partitioning (ratio PGR:CGR) measured. To assess tolerance to high air temperature during two key stages-microsporogenesis (3-6 days before flowering, DBF) and flowering, fruit-set was measured in two experiments. In the first experiment, 12 genotypes were exposed to short (3-6 days) episodes of high (38 degreesC) day air temperature at 6 DBF and at flowering. In the second experiment, 22 genotypes were exposed to 40 degreesC day air temperature for I day at 6 DBF, 3 DBF or at flowering. Cellular membrane thermostability (relative injury, RI) was also measured in these 22 genotypes. There was considerable variation among genotypes in response to high temperature, whether assessed by growth rates, fruit-set or RI. Pod weight at high soil temperature was associated with variation in CGR rather than partitioning. Flowering was more sensitive to high air temperature than microsporogenesis. Genotypes tolerant to high air temperature at microsporogenesis were not necessarily tolerant at flowering, and nor was tolerance correlated with RI. Six genotypes (796, 55-437, ICG 1236, ICGV 86021, lCGV 87281 and ICGV 92121) were identified as heat tolerant based on their performance in all tests. These experiments have shown that groundnut genotypes can be easily screened for reproductive tolerance to high air and soil temperature and that several sources of heat tolerance are available in groundnut germplasm. (C) 2003 Elsevier Science B.V. All rights reserved.

Climate change and the flowering time of annual crops

Crop production is inherently sensitive to variability in climate. Temperature is a major determinant of the rate of plant development and, under climate change, warmer temperatures that shorten development stages of determinate crops will most probably reduce the yield of a given variety. Earlier crop flowering and maturity have been observed and documented in recent decades, and these are often associated with warmer (spring) temperatures. However, farm management practices have also changed and the attribution of observed changes in phenology to climate change per se is difficult. Increases in atmospheric [CO2] often advance the time of flowering by a few days, but measurements in FACE (free air CO2 enrichment) field-based experiments suggest that elevated [CO2] has little or no effect on the rate of development other than small advances in development associated with a warmer canopy temperature. The rate of development (inverse of the duration from sowing to flowering) is largely determined by responses to temperature and photoperiod, and the effects of temperature and of photoperiod at optimum and suboptimum temperatures can be quantified and predicted. However, responses to temperature, and more particularly photoperiod, at supraoptimal temperature are not well understood. Analysis of a comprehensive data set of time to tassel initiation in maize (Zea mays) with a wide range of photoperiods above and below the optimum suggests that photoperiod modulates the negative effects of temperature above the optimum. A simulation analysis of the effects of prescribed increases in temperature (0-6 degrees C in + 1 degrees C steps) and temperature variability (0% and + 50%) on days to tassel initiation showed that tassel initiation occurs later, and variability was increased, as the temperature exceeds the optimum in models both with and without photoperiod sensitivity. However, the inclusion of photoperiod sensitivity above the optimum temperature resulted in a higher apparent optimum temperature and less variability in the time of tassel initiation. Given the importance of changes in plant development for crop yield under climate change, the effects of photoperiod and temperature on development rates above the optimum temperature clearly merit further research, and some of the knowledge gaps are identified herein.

Bias from farmer self-selection in genetically modified crop productivity estimates: Evidence from Indian data

In the continuing debate over the impact of genetically modified (GM) crops on farmers of developing countries, it is important to accurately measure magnitudes such as farm-level yield gains from GM crop adoption. Yet most farm-level studies in the literature do not control for farmer self-selection, a potentially important source of bias in such estimates. We use farm-level panel data from Indian cotton farmers to investigate the yield effect of GM insect-resistant cotton. We explicitly take into account the fact that the choice of crop variety is an endogenous variable which might lead to bias from self-selection. A production function is estimated using a fixed-effects model to control for selection bias. Our results show that efficient farmers adopt Bacillus thuringiensis (Bt) cotton at a higher rate than their less efficient peers. This suggests that cross-sectional estimates of the yield effect of Bt cotton, which do not control for self-selection effects, are likely to be biased upwards. However, after controlling for selection bias, we still find that there is a significant positive yield effect from adoption of Bt cotton that more than offsets the additional cost of Bt seed.

Whole crop cereals 2. Prediction of apparent digestibility and energy value from in vitro digestion techniques and near infrared reflectance spectroscopy and of chemical composition by near infrared reflectance spectroscopy

Samples of whole crop wheat (WCW, n = 134) and whole crop barley (WCB, n = 16) were collected from commercial farms in the UK over a 2-year period (2003/2004 and 2004/2005). Near infrared reflectance spectroscopy (NIRS) was compared with laboratory and in vitro digestibility measures to predict digestible organic matter in the dry matter (DOMD) and metabolisable energy (ME) contents measured in vivo using sheep. Spectral models using the mean spectra of two scans were compared with those using individual spectra (duplicate spectra). Overall NIRS accurately predicted the concentration of chemical components in whole crop cereals apart from crude protein. ammonia-nitrogen, water-soluble carbohydrates, fermentation acids and solubility values. In addition. the spectral models had higher prediction power for in vivo DOMD and ME than chemical components or in vitro digestion methods. Overall there Was a benefit from the use of duplicate spectra rather than mean spectra and this was especially so for predicting in vivo DOMD and ME where the sample population size was smaller. The spectral models derived deal equally well with WCW and WCB and Would he of considerable practical value allowing rapid determination of nutritive value of these forages before their use in diets of productive animals. (C) 2008 Elsevier B.V. All rights reserved.