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: The Mediterranean region has always shown a marked inter-annual variability in seasonal weather, creating uncertainty in decisional processes of cultivation and livestock breeding that should not be neglected when modeling farmers’ adaptive responses. This is especially relevant when assessing the impact of climate change (CC), which modifies the atmospheric variability and generates new uncertainty conditions, and the possibility of adaptation of agriculture. Our analysis examines this aspect reconstructing the effects of inter-annual climate variability in a diversified farming district that well represents a wide range of rainfed and irrigated agricultural systems in the Mediterranean area. We used a Regional Atmospheric Modelling System and a weather generator to generate 150 stochastic years of the present and near future climate. Then, we implemented calibrated crop and livestock models to estimate the corresponding productive responses in the form of probability distribution functions (PDFs) under the two climatic conditions. We assumed these PDFs able to represent the expectations of farmers in a discrete stochastic programming (DSP) model that reproduced their economic behaviour under uncertainty conditions. The comparison of the results in the two scenarios provided an assessment of the impact of CC, also taking into account the possibility of adjustment allowed by present technologies and price regimes. The DSP model is built in blocks that represent the farm typologies operating in the study area, each one with its own resource endowment, decisional constraints and economic response. Under this latter aspect, major differences emerged among farm typologies and sub-zones of the study area. A crucial element of differentiation was water availability, since only irrigated C3 crops took full advantage from the fertilization effect of increasing atmospheric CO2 concentration. Rainfed crop production was depressed by the expected reduction of spring rainfall associated to the higher temperatures. So, a dualism emerges between the smaller impact on crop production in the irrigated plain sub-zone, equipped with collective water networks and abundant irrigation resources, and the major negative impact in the hilly area, where these facilities and resources are absent. However intensive dairy farming was also negatively affected in terms of milk production and quality, and cattle mortality because of the increasing summer temperatures. This provides explicit guidance for addressing strategic adaptation policies and for framing farmers’ perception of CC, in order to help them to develop an awareness of the phenomena that are already in progress, which is a prerequisite for effective adaptation responses.