Abstract: This study explored the utility of the impact response surface (IRS) approach for investigating model ensemble crop yield responses under a large range of changes in climate. IRSs of spring and winter wheat Triticum aestivum yields were constructed from a 26-member ensemble of process-based crop simulation models for sites in Finland, Germany and Spain across a latitudinal transect. The sensitivity of modelled yield to systematic increments of changes in temperature (-2 to +9°C) and precipitation (-50 to +50%) was tested by modifying values of baseline (1981 to 2010) daily weather, with CO2 concentration fixed at 360 ppm. The IRS approach offers an effective method of portraying model behaviour under changing climate as well as advantages for analysing, comparing and presenting results from multi-model ensemble simulations. Though individual model behaviour occasionally departed markedly from the average, ensemble median responses across sites and crop varieties indicated that yields decline with higher temperatures and decreased precipitation and increase with higher precipitation. Across the uncertainty ranges defined for the IRSs, yields were more sensitive to temperature than precipitation changes at the Finnish site while sensitivities were mixed at the German and Spanish sites. Precipitation effects diminished under higher temperature changes. While the bivariate and multi-model characteristics of the analysis impose some limits to interpretation, the IRS approach nonetheless provides additional insights into sensitivities to inter-model and inter-annual variability. Taken together, these sensitivities may help to pinpoint processes such as heat stress, vernalisation or drought effects requiring refinement in future model development.

Abstract: ARMOSA is a dynamic simulation model which was developed to simulate crop growth and development, water and nitrogen dynamics under different pedoclimatic conditions and cropping systems in the arable land. The model is meant to be a tool for the evaluation of the impact of different crop management practices on soil nitrogen and carbon cycles and groundwater nitrate pollution. A large data set collected over three to six years from six monitoring sites in Lombardia plain was used to calibrate and validate the model parameters. Measured meteorological data, soil chemical and physical characterizations, crop-related data of different cropping systems allowed for a proper parameterization. Fit indexes showed the reliability of the model in adequately predicting crop-related variables, such as above ground biomass (RRMSE=11.18, EF=0.94, r=0.97), Leaf Area Index maximum value (RRMSE=8.24, EF=0.37, r=0.72), harvest index (RRMSE=19.4, EF=0.32, r=0.74), and crop N uptake (RRMSE=20.25, EF=0.69, r=0.85). Using two different one-year data set from each monitoring site, the model was calibrated and validated, getting to encouraging results: RRMSE=6.28, EF=0.52, r=0.68 for soil water content at different depths, and RRMSE=34.89, EF=0.59, r=0.75 for soil NO3-N content along soil profile. The simulated N leaching was in full agreement with measured data (RRMSE=26.62, EF=0.88, r=0.98).

Abstract: Soil water percolation is a key process in the life cycle of water in fields, watersheds, and river basins. The Environmental Policy Integrated Climate (EPIC) and the Agricultural Policy/Environmental eXtender (APEX) are continuous models developed for evaluating the environmental effects of agricultural management. Traditionally, these models have simulated soil water percolation processes using a tipping-bucket approach, with the rate of flow limited by the saturated hydraulic conductivity. This simple approach often leads to inaccuracy in simulating elevated soil water conditions where soil water content (SWC) levels may remain above field capacity under prolonged wet weather periods or limited drainage. To overcome this deficiency, a new sub-model, the variable saturation hydraulic conductivity (VSHC) method, was developed for simulating soil water percolation processes using a nonlinear equation to estimate the effective hydraulic conductivity as a function of the SWC and soil properties. The VSHC method was evaluated at three sites in the United States and two sites in Europe. In addition, a numerical solution of the Richards equation was used as a benchmark for SWC comparison. Results show that the VSHC method substantially improves the accuracy of the SWC simulation in long-term simulations, particularly during wet periods. At the watershed scale, results on the Riesel Y2 watershed indicate that the VSHC method enhances model performance in the high-flow regime of channel peak flows because of the improved estimation of SWC, which implies that the improved SWC simulation at the field scale is beneficial to hydrologic modeling at the watershed scale.