Abstract: Rice models have been widely used in simulating and predicting rice phenology in contrasting climate zones, however the uncertainties from model structure (different equations or models) and/or model parameters were rarely investigated. Here, five rice phenological models/modules (Le., CERES-Rice, ORYZA2000, RCM, Beta Model and SIMRIW) were applied to simulate rice phenology at 23 experimental stations from 1992 to 2009 in two major rice cultivation regions of China: the northeastern China and the southwestern China. To investigate the uncertainties from model biophysical parameters, each model was run with randomly perturbed 50 sets of parameters. The results showed that the median of ensemble simulations were better than the simulation by most models. Models couldn’t simulate well in some specific years despite of parameters optimization, suggesting model structure limit model performance in some cases. The models adopting accumulative thermal time function (e.g., CERES-Rice and ORYZA2000) had better performance in the southwestern China, in contrast, those adopting exponential function (e.g., Beta model and RCM model) had better performance in the northeastern China. In northeastern China, the contribution of model structure and model parameters to model total variance was, respectively, about 55.90% and 44.10% in simulating heading date, and about 75.43% and 24.57% in simulating maturity date. In the southwestern China, the contribution of model structure and model parameters to model total variance was, respectively, about 79.97% and 27.03% in simulating heading date, about 92.15% and 7.85% in simulating maturity date. Uncertainty from model structure was the most relevant source. The results highlight that the temperature response functions of rice development rate under extreme climate conditions should be improved based on environment-controlled experimental data.

Abstract: The available evidence suggests that the current increasing trend in global surface temperatures will continue during this century, which will be accompanied by a greater frequency of extreme events. The IPCC has projected that higher temperatures may outscore the known optimal and maximum temperatures for maize. The purpose of this study was to improve the ability of the maize model CSM-IXIM to simulate crop development, growth, and yield under hot conditions, especially with regards to the impact of above-optimal temperatures around anthesis. Field and greenhouse experiments that were performed over three years (2014-2016) using the same short-season hybrid, PR37N01 (FAO 300), provided the data for this work. Maize was sown at a target population density of 5 plants M-2 on two sowing dates in 2014 and 2015 and on one in 2016 at three locations in Spain (northern, central, and southern Spain) with a well-defined thermal gradient. The same hybrid was also sown in two greenhouse chambers with daytime target temperatures of approximately 25 and above 35 degrees C. During the nighttime, the temperature in both chambers was allowed to equilibrate with the outside temperature. The greenhouse treatments consisted of moving 18 plants at selected phenological stages (V4, V9, anthesis, lag phase, early grain filling) from the cool chamber to the hot chamber over a week and then returning the plants back to the cool chamber. An additional control treatment remained in the cool chamber all season, and in 2015 and 2016, one treatment remained permanently in the hot chamber. Two maize models in the Decision Support System for Agrotechnology Transfer (DSSAT) V4.6 were compared, namely CERES and IXIM. The HUM version included additional components that were previously developed to improve the crop N simulation and to incorporate the anthesis-silking interval (ASI). A new thermal time calculation, a heat stress index, the impact of pollen-sterilizing temperatures, and the explicit simulation of male and female flowering as affected by the daily heat conditions were added to IXIM. The phenology simulation in field experiments by IXIM improved substantially. The RMSE for silking and maturity in CERES were 7.9 and 13.7 days, decreasing in DCIM to 2.8 and 7.3 days, respectively. Similarly, the estimated kernel numbers, kernel weight, grain yield and final biomass were always closer to the measurements in HUM than in CERES. The worst simulations were for kernel weight, and for that reason, the differences in grain yield between the models were small (the RMSE in CERES was 1219 kg ha(-1) vs. 1082 kg ha(-1) in IXIM). The greenhouse results also supported the improved estimations of crop development by IXIM (RMSE of 2.6 days) relative to CERES (7.4 days). The impact of the heat treatments on grain yield was consistently overestimated by CERES, while HUM captured the general trend. The new HUM model improved the CERES simulations when elevated temperatures were included in the evaluation data. Additional model testing with measurements from a wider latitudinal range and relevant heat conditions are required.

Abstract: Despite widespread application in studying climate change impacts, most crop models ignore complex interactions among air temperature, crop and soil water status, CO2 concentration and atmospheric conditions that influence crop canopy temperature. The current study extended previous studies by evaluating Tc simulations from nine crop models at six locations across environmental and production conditions. Each crop model implemented one of an empirical (EMP), an energy balance assuming neutral stability (EBN) or an energy balance correcting for atmospheric stability conditions (EBSC) approach to simulate Tc. Model performance in predicting Tc was evaluated for two experiments in continental North America with various water, nitrogen and CO2 treatments. An empirical model fit to one dataset had the best performance, followed by the EBSC models. Stability conditions explained much of the differences between modeling approaches. More accurate simulation of heat stress will likely require use of energy balance approaches that consider atmospheric stability conditions.