Abstract: Climate change impact assessments are plagued with uncertainties from many sources, such as climate projections or the inadequacies in structure and parameters of the impact model. Previous studies tried to account for the uncertainty from one or two of these. Here, we developed a triple-ensemble probabilistic assessment using seven crop models, multiple sets of model parameters and eight contrasting climate projections together to comprehensively account for uncertainties from these three important sources. We demonstrated the approach in assessing climate change impact on barley growth and yield at Jokioinen, Finland in the Boreal climatic zone and Lleida, Spain in the Mediterranean climatic zone, for the 2050s. We further quantified and compared the contribution of crop model structure, crop model parameters and climate projections to the total variance of ensemble output using Analysis of Variance (ANOVA). Based on the triple-ensemble probabilistic assessment, the median of simulated yield change was -4% and +16%, and the probability of decreasing yield was 63% and 31% in the 2050s, at Jokioinen and Lleida, respectively, relative to 1981-2010. The contribution of crop model structure to the total variance of ensemble output was larger than that from downscaled climate projections and model parameters. The relative contribution of crop model parameters and downscaled climate projections to the total variance of ensemble output varied greatly among the seven crop models and between the two sites. The contribution of downscaled climate projections was on average larger than that of crop model parameters. This information on the uncertainty from different sources can be quite useful for model users to decide where to put the most effort when preparing or choosing models or parameters for impact analyses. We concluded that the triple-ensemble probabilistic approach that accounts for the uncertainties from multiple important sources provide more comprehensive information for quantifying uncertainties in climate change impact assessments as compared to the conventional approaches that are deterministic or only account for the uncertainties from one or two of the uncertainty sources.

Abstract: The data set for this study comprised 1,488,474 test-day records for milk, fat, and protein yields and fat and protein percentages from 191,012 first-, second-, and third-parity Holstein cows from 484 farms. Data were collected from 2001 through 2007 and merged with meteorological data from 35 weather stations. A linear model (M1) was used to estimate the effects of the temperature-humidity index (THI) on production traits. Least squares means from M1 were used to detect the THI thresholds for milk production in all parities by using a 2-phase linear regression procedure (M2). A multiple-trait repeatability test-model (M3) was used to estimate variance components for all traits and a dummy regression variable (t) was defined to estimate the production decline caused by heat stress. Additionally, the estimated variance components and M3 were used to estimate traditional and heat-tolerance breeding values (estimated breeding values, EBV) for milk yield and protein percentages at parity 1. An analysis of data (M2) indicated that the daily THI at which milk production started to decline for the 3 parities and traits ranged from 65 to 76. These THI values can be achieved with different temperature/humidity combinations with a range of temperatures from 21 to 36°C and relative humidity values from 5 to 95%. The highest negative effect of THI was observed 4 d before test day over the 3 parities for all traits. The negative effect of THI on production traits indicates that first-parity cows are less sensitive to heat stress than multiparous cows. Over the parities, the general additive genetic variance decreased for protein content and increased for milk yield and fat and protein yield. Additive genetic variance for heat tolerance showed an increase from the first to third parity for milk, protein, and fat yield, and for protein percentage. Genetic correlations between general and heat stress effects were all unfavorable (from -0.24 to -0.56). Three EBV per trait were calculated for each cow and bull (traditional EBV, traditional EBV estimated with the inclusion of THI covariate effect, and heat tolerance EBV) and the rankings of EBV for 283 bulls born after 1985 with at least 50 daughters were compared. When THI was included in the model, the ranking for 17 and 32 bulls changed for milk yield and protein percentage, respectively. The heat tolerance genetic component is not negligible, suggesting that heat tolerance selection should be included in the selection objectives.

Abstract: This study investigated the mortality of dairy cows during heat waves. Mortality data (46,610 cases) referred to dairy cows older than 24 mo that died on a farm from all causes from May 1 to September 30 during a 6-yr period (2002-2007). Weather data were obtained from 12 weather stations located in different areas of Italy. Heat waves were defined for each weather station as a period of at least 3 consecutive days, from May 1 to September 30 (2002-2007), when the daily maximum temperature exceeded the 90th percentile of the reference distribution (1971-2000). Summer days were classified as days in heat wave (HW) or not in heat wave (nHW). Days in HW were numbered to evaluate the relationship between mortality and length of the wave. Finally, the first 3 nHW days after the end of a heat wave were also considered to account for potential prolonged effects. The mortality risk was evaluated using a case-crossover design. A conditional logistic regression model was used to calculate odds ratio and 95% confidence interval for mortality recorded in HW compared with that recorded in nHW days pooled and stratified by duration of exposure, age of cows, and month of occurrence. Dairy cows mortality was greater during HW compared with nHW days. Furthermore, compared with nHW days, the risk of mortality continued to be higher during the 3 d after the end of HW. Mortality increased with the length of the HW. Considering deaths stratified by age, cows up to 28 mo were not affected by HW, whereas all the other age categories of older cows (29-60, 61-96, and >96 mo) showed a greater mortality when exposed to HW. The risk of death during HW was higher in early summer months. In particular, the highest risk of mortality was observed during June HW. Present results strongly support the implementation of adaptation strategies which may limit heat stress-related impairment of animal welfare and economic losses in dairy cow farm during HW.