Abstract: THIS PAPER AIMS: (i) to identify at national scale areas where crop yield formation is currently most prone to climate-induced stresses, (ii) to evaluate how the severity of these stresses is likely to develop in time and space, and (iii) to appraise and quantify the performance of two strategies for adapting crop cultivation to a wide range of (uncertain) climate change projections. To this end we made use of extensive climate, crop, and soil data, and of two modelling tools: N-AgriCLIM and the WOFOST crop simulation model. N-AgriCLIM was developed for the automatic generation of indicators describing basic agroclimatic conditions and was applied over the whole of Finland. WOFOST was used to simulate detailed crop responses at four representative locations. N-AgriCLIM calculations have been performed nationally for 3829 grid boxes at a 10 × 10 km resolution and for 32 climate scenarios. Ranges of projected shifts in indicator values for heat, drought and other crop-relevant stresses across the scenarios vary widely – so do the spatial patterns of change. Overall, under reference climate the most risk-prone areas for spring cereals are found in south-west Finland, shifting to south-east Finland towards the end of this century. Conditions for grass are likely to improve. WOFOST simulation results suggest that CO2 fertilization and adjusted sowing combined can lead to small yield increases of current barley cultivars under most climate scenarios on favourable soils, but not under extreme climate scenarios and poor soils. This information can be valuable for appraising alternative adaptation strategies. It facilitates the identification of regions in which climatic changes might be rapid or otherwise notable for crop production, requiring a more detailed evaluation of adaptation measures. The results also suggest that utilizing the diversity of cultivar responses seems beneficial given the high uncertainty in climate change projections.

Abstract: Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model (DVM) with 10 unknown parameters, using the DREAM((ZS)) (DiffeRential Evolution Adaptive Metropolis) Markov chain Monte Carlo (MCMC) sampler. We focus on comparing model inversions, considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a priori or jointly inferred together with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root mean square errors (RMSEs) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19, 1.04 to 1.56 g C m(-2) day(-1) and 0.50 to 1.28 mm day(-1), respectively. For the calibration period, using a homoscedastic eddy covariance residual error model resulted in a better agreement between measured and modelled data than using a heteroscedastic residual error model. However, a model validation experiment showed that CARAIB models calibrated considering heteroscedastic residual errors perform better. Posterior parameter distributions derived from using a heteroscedastic model of the residuals thus appear to be more robust. This is the case even though the classical linear heteroscedastic error model assumed herein did not fully remove heteroscedasticity of the GPP residuals. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling. Besides the residual error treatment, differences between model parameter posterior distributions among the four grassland sites are also investigated. It is shown that the marginal distributions of the specific leaf area and characteristic mortality time parameters can be explained by site-specific ecophysiological characteristics.

Abstract: The upland rice (UR) cropped area in Brazil has decreased in the last decade. Importantly, a portion of this decrease can be attributed to the current UR breeding programme strategy, according to which direct grain yield selection is targeted primarily to the most favourable areas. New strategies for more-efficient crop breeding under non-optimal conditions are needed for Brazil’s UR regions. Such strategies should include a classification of spatio-temporal yield variations in environmental groups, as well as a determination of prevalent drought types and their characteristics (duration, intensity, phenological timing, and physiological effects) within those environmental groups. This study used a process-based crop model to support the Brazilian UR breeding programme in their efforts to adopt a new strategy that accounts for the varying range of environments where UR is currently cultivated. Crop simulations based on a commonly grown cultivar (BRS Primavera) and statistical analyses of simulated yield suggested that the target population of environments can be divided into three groups of environments: a highly favorable environment (HFE, 19% of area), a favorable environment (FE, 44%), and least favourable environment (LFE, 37%). Stress-free conditions dominated the HFE group (69% likelihood) and reproductive stress dominated the LFE group (68% likelihood), whereas reproductive and terminal drought stress were found to be almost equally likely to occur in the FE group. For the best and worst environments, we propose specific adaptation focused on the representative stress, while for the FE, wide adaptation to drought is suggested. ‘Weighted selection’ is also a possible strategy for the FE and LFE environment groups.