Abstract: Renewed interest in heat stress effects on livestock productivity derives from climate change, which is expected to increase temperatures and the frequency of extreme weather events. This study aimed at evaluating the effect of temperature and humidity on milk production in highly selected dairy cattle populations across 3 European regions differing in climate and production systems to detect differences and similarities that can be used to optimize heat stress (HS) effect modeling. Milk, fat, and protein test day data from official milk recording for 1999 to 2010 in 4 Holstein populations located in the Walloon Region of Belgium (BEL), Luxembourg (LUX), Slovenia (SLO), and southern Spain (SPA) were merged with temperature and humidity data provided by the state meteorological agencies. After merging, the number of test day records/cows per trait ranged from 686,726/49,655 in SLO to 1,982,047/136,746 in BEL. Values for the daily average and maximum temperature-humidity index (THIavg and THImax) ranges for THIavg/THImax were largest in SLO (22-74/28-84) and shortest in SPA (39-76/46-83). Change point techniques were used to determine comfort thresholds, which differed across traits and climatic regions. Milk yield showed an inverted U-shaped pattern of response across the THI scale with a HS threshold around 73 THImax units. For fat and protein, thresholds were lower than for milk yield and were shifted around 6 THI units toward larger values in SPA compared with the other countries. Fat showed lower HS thresholds than protein traits in all countries. The traditional broken line model was compared with quadratic and cubic fits of the pattern of response in production to increasing heat loads. A cubic polynomial model allowing for individual variation in patterns of response and THIavg as heat load measure showed the best statistical features. Higher/lower producing animals showed less/more persistent production (quantity and quality) across the THI scale. The estimated correlations between comfort and THIavg values of 70 (which represents the upper end of the THIavg scale in BEL-LUX) were lower for BEL-LUX (0.70-0.80) than for SPA (0.83-0.85). Overall, animals producing in the more temperate climates and semi-extensive grazing systems of BEL and LUX showed HS at lower heat loads and more re-ranking across the THI scale than animals producing in the warmer climate and intensive indoor system of SPA.

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: Highlights • Regularization techniques for regression outperformed the classical regression techniques in predicting crop yields. • Different regression techniques with similar prediction accuracy showed different responses of major climatic variables to crop yield. • The regression models showed some responses of crop yield to climatic conditions that is mostly absent in process based crop models. Abstract Statistical regression models represent alternatives to process-based dynamic models for predicting the response of crop yields to variation in climatic conditions. Regression models can be used to quantify the effect of change in temperature and precipitation on yields. However, it is difficult to identify the most relevant input variables that should be included in regression models due to the high number of candidate variables and to their correlations. This paper compares several regression techniques for modeling response of winter oilseed rape yield to a high number of correlated input variables. Several statistical regression methods were fitted to a dataset including 689 observations of winter oilseed rape yield from replicated field experiments conducted in 239 sites in Denmark, covering nearly all regions of the country from 1992 to 2013. Regression methods were compared by cross-validation. The regression methods leading to the most accurate yield predictions were Lasso and Elastic Net, and the least accurate methods were ordinary least squares and stepwise regression. Partial least squares and ridge regression methods gave intermediate results. The estimated relative yield change for a +1°C temperature increase during flowering was estimated to range between 0 and +6 %, depending on choice of regression method. Precipitation was found to have an adverse effect on yield during autumn and winter. It was estimated that an increase in precipitation of +1 mm/day would result in a relative yield change ranging from 0 to −4 %. Soil type was also important for crop yields with lower yields on sandy soils compared to loamy soils. Later sowing was found to result in increased crop yield. The estimated effect of climate on yield was highly sensitive to the chosen regression method. Regression models showing similar performance led in some cases to different conclusions with respect to effect of temperature and precipitation. Hence, it is recommended to apply an ensemble of regression models, in order to account for the sensitivity of the data driven models for projecting crop yield under climate change.