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 degrees 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 developmen

Adaptation of crops to climate change has to be addressed locally due to the variability of soil, climate and the specific socioeconomic settings influencing farm management decisions. Adaptation of rainfed cropping systems in the Mediterranean is especially challenging due to the projected decline in precipitation in the coming decades, which will increase the risk of droughts. Methods that can help explore uncertainties in climate projections and crop modelling, such as impact response surfaces (IRSs) and ensemble modelling, can then be valuable for identifying effective adaptations. Here, an ensemble of 17 crop models was used to simulate a total of 54 adaptation options for rainfed winter wheat (Triticum aestivum) at Lleida (NE Spain). To support the ensemble building, an ex post quality check of model simulations based on several criteria was performed. Those criteria were based on the "According to Our Current Knowledge" (AOCK) concept, which has been formalized here. Adaptations were based on changes in cultivars and management regarding phenology, vernalization, sowing date and irrigation. The effects of adaptation options under changed precipitation (P), temperature (T), [CO 2 ] and soil type were analysed by constructing response surfaces, which we termed, in accordance with their specific purpose, adaptation response surfaces (ARSs). These were created to assess the effect of adaptations through a range of plausible P, T and [CO 2 ] perturbations. The results indicated that impacts of altered climate were predominantly negative. No single adaptation was capable of overcoming the detrimental effect of the complex interactions imposed by the P, T and [CO 2 ] perturbations except for supplementary irrigation (sI), which reduced the potential impacts under most of the perturbations. Yet, a combination of adaptations for dealing with climate change demonstrated that effective adaptation is possible at Lleida. Combinations based on a cultivar without vernalization requirements showed good and wide adaptation potential. Few combined adaptation options

Agro-meteorological quantities are often in the form of time series, and knowledge about their temporal scaling properties is fundamental for transferring locally measured fluctuations to larger scales and vice versa. However, the scaling analysis of these quantities is complicated due to the presence of localized trends and nonstationarities. The objective of this study was to characterise scaling properties (i.e. statistical self-similarity) of the chosen agro-meteorological quantities through multifractal detrended fluctuation analysis (MFDFA). For this purpose, MFDFA was performed with 11 322 measured time series (31 yr) of daily air temperature, wind velocity, relative air humidity, global radiation and precipitation from stations located in Finland, Germany, Poland and Spain. The empirical singularity spectra indicated their multifractal structure. The richness of the studied multifractals was evaluated by the width of their spectrum, indicating considerable differences in dynamics and development. In log-log plots of the cumulative distributions of all meteorological parameters the linear functions prevailed for high values of the response, indicating that these distributions were consistent with power-law asymptotic behaviour. Additionally, we investigated the type of multifractality that underlies the q-dependence of the generalized Hurst exponent by analysing the corresponding shuffled and surrogate time series. For most of the studied meteorological parameters, the multifractality is due to different long-range correlations for small and large fluctuations. Only for precipitation does the multifractality result mainly from broad probability function. This feature may be especially valuable for assessing the effect of change in climate dynamics.

The daily air temperature and precipitation time series recorded between January 1, 1980 and December 31, 2010 in four European sites (Jokioinen, Dikopshof, Lleida and Lublin) from different climatic zones were modeled and forecasted. In our forecasting we used the methods of the Box-Jenkins and Holt- Winters seasonal auto regressive integrated moving-average, the autoregressive integrated moving-average with external regressors in the form of Fourier terms and the time series regression, including trend and seasonality components methodology with R software. It was demonstrated that obtained models are able to capture the dynamics of the time series data and to produce sensible forecasts.