An improved model to simulate rice yield

Rice is the staple food for about half of the world’s population. Although global production has more than doubled in the last 40 years, food security problems still persist and need to be managed based on early and reliable forecasting activities. This is especially true since the frequency of extreme weather events is forecasted to increase by the intergovernmental panel on climate change (IPCC). The most advanced crop yield forecasting systems are based on simulation models. However, examples of operational systems implementing models which are suitable for reproducing the peculiarities of paddy rice, especially on small scales, are missing. The rice model WARM is used within the crop yield forecasting system of the European Commission. In this article we evaluated the WARM model for the simulation of rice growth under flooded and unflooded conditions in China and Italy. The WARM model simulates crop growth and development, floodwater effect on the vertical thermal profile, blast disease, cold-shock induced spikelet sterility during the pre-flowering period and hydrological peculiarities of paddy soils. We identified the most relevant model parameters through sensitivity analyses carried out using the Sobol’ method and then calibrated using the simplex algorithm. Data from 11 published experiments, covering 13 locations and 10 years, were used. Two groups of rice varieties were identified for each country. Our results show that the model was able to reproduce rice growth in both countries. Specifically, the average relative root mean square error calculated on aboveground biomass curves was 21.9% for the calibration and 23.6% for validation. The parameters of the linear regression equation between measured and simulated values were always satisfactory. Indeed, intercept and slope were always close to their optima and R2 was always higher than 0.79. For some of the combinations of country and simulated variable, the indices of agreement calculated for the validation datasets were better then the corresponding ones computed at the end of the calibration, indirectly proving the robustness of the modeling approach. WARM’s robustness and accuracy, combined with the low requirements in terms of inputs and the implementation of modules for reproducing biophysical processes strongly influencing the year-to-year yield variation, make the model suitable for forecasting rice yields on regional, national and international scales.

[1]  C. Vignolles,et al.  A methodology for a combined use of normalised difference vegetation index and CORINE land cover data for crop yield monitoring and forecasting. A case study on Spain , 2001 .

[2]  John R. Williams,et al.  A method for estimating the direct and climatic effects of rising atmospheric carbon dioxide on growth and yield of crops: Part I--Modification of the EPIC model for climate change analysis , 1992 .

[3]  S. Sorooshian,et al.  Effective and efficient global optimization for conceptual rainfall‐runoff models , 1992 .

[4]  J. Monteith The quest for balance in crop modeling , 1996 .

[5]  Jeffrey G. Arnold,et al.  Rice Parameters Describing Crop Performance of Four U.S. Cultivars , 2001 .

[6]  B. Bouman,et al.  Effect of water-saving irrigation on rice yield and water use in typical lowland conditions in Asia , 2004 .

[7]  Neil M.J. Crout,et al.  A stochastic modelling approach for real-time forecasting of winter wheat yield , 1999 .

[8]  Carlos H. Díaz-Ambrona,et al.  Uncertainties in projected impacts of climate change on European agriculture and terrestrial ecosystems based on scenarios from regional climate models , 2007 .

[9]  John R. Williams,et al.  A modeling approach to determining the relationship between erosion and soil productivity [EPIC, Erosion-Productivity Impact Calculator, mathematical models] , 1984 .

[10]  Mirco Boschetti,et al.  ESTIMATION OF PARAMETERS DESCRIBING MORPHO-PHYSIOLOGICAL FEATURES OF MEDITERRANEAN RICE VARIETIES FOR MODELLING PURPOSES STIMA DI PARAMETRI MORFO-FISIOLOGICI DI INTERESSE MODELLISTICO PER VARIETÀ DI RISO MEDITERRANEE , 2006 .

[11]  Daniel Rodriguez,et al.  LINGRA-CC : a sink-source model to simulate the impact of climate change and management on grassland productivity , 1999 .

[12]  B. Bouman,et al.  Exploring options to combine high yields with high nitrogen use efficiencies in irrigated rice in China , 2007 .

[13]  Luca Bechini,et al.  Parameterization of a crop growth and development simulation model at sub-model components level. An example for winter wheat (Triticum aestivum L.) , 2006, Environ. Model. Softw..

[14]  M. Dingkuhn,et al.  RICE CROP DURATION AND LEAF APPEARANCE RATE IN A VARIABLE THERMAL ENVIRONMENT. II. COMPARISON OF GENOTYPES , 1998 .

[15]  Luigi Mariani,et al.  Analysis and modelling of water and near water temperatures in flooded rice (Oryza sativa L.) , 2005 .

[16]  Confalonieri Roberto,et al.  Exploratory Sensitivity Analysis of CropSyst, WARM and WOFOST: a Case-Study with Rice Biomass Simulations , 2006 .

[17]  M. Dingkuhn,et al.  Rice crop duration and leaf appearance rate in a variable thermal environment.: I. Development of an empirically based model , 1998 .

[18]  Stefano Bocchi,et al.  Evaluation of CropSyst for simulating the yield of flooded rice in northern Italy , 2005 .

[19]  H. H. Laar,et al.  ORYZA 1 : an ecophysiological model for irrigated rice production , 1994 .

[20]  J. Goudriaan,et al.  ON APPROACHES AND APPLICATIONS OF THE WAGENINGEN CROP MODELS , 2003 .

[21]  Confalonieri Roberto,et al.  Optimization Algorithms for Calibrating Cropping Systems Simulation Models - A Case Study with Simplex-derived Methods Integrated in the WARM Simulation Environment , 2006 .

[22]  P. Aggarwal,et al.  Climate Change and Rice Yields in Diverse Agro-Environments of India. I. Evaluation of Impact Assessment Models , 2002 .

[23]  F. E. Satterthwaite An approximate distribution of estimates of variance components. , 1946, Biometrics.

[24]  Welch Bl THE GENERALIZATION OF ‘STUDENT'S’ PROBLEM WHEN SEVERAL DIFFERENT POPULATION VARLANCES ARE INVOLVED , 1947 .

[25]  Andrej Ceglar,et al.  Crop growth simulation model WOFOST. , 2010 .

[26]  M. Dingkuhn,et al.  Growth and yield potential of Oryza sativa and O. glaberrima upland rice cultivars and their interspecific progenies , 1998 .

[27]  C. A. van Diepen,et al.  Spatial resolution of precipitation and radiation: The effect on regional crop yield forecasts , 2005 .

[28]  C. N. Bezuidenhout,et al.  Operational forecasting of South African sugarcane production: Part 1 – System description , 2007 .

[29]  B. Bouman,et al.  Exploring options to grow rice using less water in northern China using a modelling approach: I. Field experiments and model evaluation , 2007 .

[30]  Ilya M. Sobol,et al.  Sensitivity Estimates for Nonlinear Mathematical Models , 1993 .

[31]  J. Heilman,et al.  Seasonal variation in radiation use efficiency of irrigated rice , 2001 .

[32]  Satterthwaite Fe An approximate distribution of estimates of variance components. , 1946 .

[33]  Michaël Dingkuhn,et al.  Relationships between upland rice canopy characteristics and weed competitiveness , 1999 .

[34]  J. Goudriaan,et al.  Monitoring rice reflectance at field level for estimating biomass and LAI , 1998 .

[35]  Stefano Bocchi,et al.  The CropSyst Model to Simulate the N Balance of Rice for Alternative Management , 2006 .

[36]  B. Bouman,et al.  Performance of aerobic rice varieties under irrigated conditions in North China , 2006 .

[37]  R. Pachauri Climate change 2007. Synthesis report. Contribution of Working Groups I, II and III to the fourth assessment report , 2008 .

[38]  H. van Keulen,et al.  Modelling of agricultural production: Weather, soils, and crops , 1986 .