Probabilistic simulations of crop yield over western India using the DEMETER seasonal hindcast ensembles

Process-based integrated modelling of weather and crop yield over large areas is becoming an important research topic. The production of the DEMETER ensemble hindcasts of weather allows this work to be carried out in a probabilistic framework. In this study, ensembles of crop yield (groundnut, Arachis hypogaea L.) were produced for 10 2.5 degrees x 2.5 degrees grid cells in western India using the DEMETER ensembles and the general large-area model (GLAM) for annual crops. Four key issues are addressed by this study. First, crop model calibration methods for use with weather ensemble data are assessed. Calibration using yield ensembles was more successful than calibration using reanalysis data (the European Centre for Medium-Range Weather Forecasts 40-yr reanalysis, ERA40). Secondly, the potential for probabilistic forecasting of crop failure is examined. The hindcasts show skill in the prediction of crop failure, with more severe failures being more predictable. Thirdly, the use of yield ensemble means to predict interannual variability in crop yield is examined and their skill assessed relative to baseline simulations using ERA40. The accuracy of multi-model yield ensemble means is equal to or greater than the accuracy using ERA40. Fourthly, the impact of two key uncertainties, sowing window and spatial scale, is briefly examined. The impact of uncertainty in the sowing window is greater with ERA40 than with the multi-model yield ensemble mean. Subgrid heterogeneity affects model accuracy: where correlations are low on the grid scale, they may be significantly positive on the subgrid scale. The implications of the results of this study for yield forecasting on seasonal time-scales are as follows. (i) There is the potential for probabilistic forecasting of crop failure (defined by a threshold yield value); forecasting of yield terciles shows less potential. (ii) Any improvement in the skill of climate models has the potential to translate into improved deterministic yield prediction. (iii) Whilst model input uncertainties are important, uncertainty in the sowing window may not require specific modelling. The implications of the results of this study for yield forecasting on multidecadal (climate change) time-scales are as follows. (i) The skill in the ensemble mean suggests that the perturbation, within uncertainty bounds, of crop and climate parameters, could potentially average out some of the errors associated with mean yield prediction. (ii) For a given technology trend, decadal fluctuations in the yield-gap parameter used by GLAM may be relatively small, implying some predictability on those time-scales.

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

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.

Drought, pod yield, pre-harvest Aspergillus infection and aflatoxin contamination on peanut in Niger

Soil moisture and soil temperature affect pre-harvest infection with Aspergillus flavus and production of aflatoxin. The objectives of our field research in Niger, West Africa, were to: (i) examine the effects of sowing date and irrigation treatments on pod yield, infection with A. flavus and aflatoxin concentration; and (ii) to quantify relations between infection, aflatoxin concentration and soil moisture stress. Seed of an aflatoxin susceptible peanut cv. JL24 was sown at two to four different sowing dates under four irrigation treatments (rainfed and irrigation at 7, 14 and 21 days intervals) between 1991 and 1994, giving 40 different 'environments'. Average air and soil temperatures of 28-34 degrees C were favourable for aflatoxin contamination. CROPGRO-peanut model was used to simulate the occurrence of moisture stress. The model was able to simulate yields of peanut well over the 40 environments (r(2) = 0.67). In general, early sowing produced greater pod yields, as well as less infection and lower aflatoxin concentration. There were negative linear relations between infection (r(2) = 0.62) and the average simulated fraction of extractable soil water (FESW) between flowering and harvest, and between aflatoxin concentration (r(2) = 0.54) and FESW in the last 25 days of pod-filling. This field study confirms that infection and aflatoxin concentration in peanut can be related to the occurrence of soil moisture stress during pod-filling when soil temperatures are near optimal for A. flavus. These relations could form the basis of a decision-support system to predict the risk of aflatoxin contamination in peanuts in similar environments. (c) 2005 Elsevier B.V. All rights reserved.

Effect of plant density and initial crown size on growth, development and yield in strawberry cultivars Elsanta and Bolero

The effects of density (plant spacing) and initial plant size on vegetative growth, flowering and fruiting were studied in the strawberry cultivars Elsanta and Bolero in their first and second years of cropping. The influence of these factors on light use and dry-matter partitioning was investigated. The size of planting material in 'Elsanta' and 'Bolero' slightly affected plant growth and yield, but this effect was not consistent and radiation use efficiency (RUE) and harvest index were unaltered. Plant spacing did not significantly affect the early stages of crop growth, but was important in determining growth and yield later in the season, this effect being more significant in the second year of cropping. Plant growth and yield per plant increased as plant spacing increased from 20 to 30 cm in both 'Elsanta' and 'Bolero', but the highest harvest index and yield per square metre were obtained at the closest spacing. Increased plant spacing also resulted in a greater leaf area and leaf area index. However, light was used less efficiently resulting in a lower RUE and lower harvest index (HI).

Influence of foliar diseases and their control by fungicides on grain yield and quality in wheat

Twenty-eight field experiments on sandy-loam soils in the UK (1982-2003) are reviewed by relating the extension of the green area duration of the flag leaf (GLADF) by fungicides to effects on yield and quality of winter wheat. Over all experiments mean grain yield = 8.85t ha(-1) at 85% DM. With regards quality, mean values were: thousand grain weight (TGW) = 44.5 g; specific weight (SWT) = 76.9 kg hl(-1); crude protein concentration (CP (N x 5.7)) = 12.5 % DM; Hagberg falling number (HFN) = 285 s; and sodium dodecyl sulphate (SDS)-sedimentation volume = 69ml. For each day (d) that fungicides increased GLADF there were associated average increases in yield (0.144 1 ha(-1) d(-1), se 0.0049, df = 333), TGW (0.56 gd(-1), se = 0.017) and SWT (0.22 kg hl(-1) d(-1), se 0.011). Some curvature was evident in all these relationships. When GLADF was delayed beyond 700 degrees Cd after anthesis, as was possible in cool wet seasons, responses were curtailed, or less reliable. Despite this apparent terminal sink limitation, fungicide effects on sink size, eg endosperm cell numbers or maximum water mass per grain, were not prerequisites for large effects on grain yield, TGW or SWT. Fungicide effects on CP were variable. Although the average response of CP was negative (-0.029%DM/d; se = 0.00338), this depended on cultivar and disease controlled. Controlling biotrophs such as rusts, (Puccinia spp.) tended to increase CP, whereas controlling a more necrotrophic pathogen (Septoria tritici) usually reducedCP. Irrespective of pathogen controlled, delaying senescence of the flag leaf was associated with increased nitrogen yields in the grain (averaging 2.24 kg N ha-1 d(-1), se = 0.0848) due to both increased N uptake into the above ground crop, and also more efficient remobilisation of N from leaf laminas. When sulphur availability appeared to be adequate, fungicide x cultivar interactions were similar on S as for CP, although N:S ratios tended to decline (i.e. improve for bread making) when S. tritici was controlled. On average, SDS-sedimentation volume declined (-0. 18 ml/d, se = 0.027) with increased GLADF, broadly commensurate with the average effect on CP. Hagberg falling number decreased as fungicide increased GLADF (-2.73 s/d, se = 0.178), indicating an increase in alpha-amylase activity.