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: • This study combines crop, livestock and economic models.
• Models interaction is through use of relevant input and output variables.
• Future climate change will result in increased grass and wheat dry matter yields.
• Changes in grass, wheat and milk yields in future reduce farm emissions intensity.
• Changes in future dry matter yields and emissions lead to increased profitability.
There is a scientific consensus that the future climate change will affect grass and crop dry matter (DM) yields. Such yield changes may entail alterations to farm management practices to fulfill the feed requirements and reduce the farm greenhouse gas (GHG) emissions from dairy farms. While a large number of studies have focused on the impacts of projected climate change on a single farm output (e.g. GHG emissions or economic performance), several attempts have been made to combine bio-economic systems models with GHG accounting frameworks. In this study, we aimed to determine the physical impacts of future climate scenarios on grass and wheat DM yields, and demonstrate the effects such changes in future feed supply may have on farm GHG emissions and decision-making processes. For this purpose, we combined four models: BASGRA and CSM-CERES-Wheat models for simulating forage grass DM and wheat DM grain yields respectively; HolosNor for estimating the farm GHG emissions; and JORDMOD for calculating the impacts of changes in the climate and management on land use and farm economics. Four locations, with varying climate and soil conditions were included in the study: south-east Norway, south-west Norway, central Norway and northern Norway. Simulations were carried out for baseline (1961–1990) and future (2046–2065) climate conditions (projections based on two global climate models and the Special Report on Emissions Scenarios (SRES) A1B GHG emission scenario), and for production conditions with and without a milk quota. The GHG emissions intensities (kilogram carbon dioxide equivalent: kgCO2e emissions per kg fat and protein corrected milk: FPCM) varied between 0.8 kg and 1.23 kg CO2e (kg FPCM)− 1, with the lowest and highest emissions found in central Norway and south-east Norway, respectively. Emission intensities were generally lower under future compared to baseline conditions due mainly to higher future milk yields and to some extent to higher crop yields. The median seasonal above-ground timothy grass yield varied between 11,000 kg and 16,000 kg DM ha− 1 and was higher in all projected future climate conditions than in the baseline. The spring wheat grain DM yields simulated for the same weather conditions within each climate projection varied between 2200 kg and 6800 kg DM ha− 1. Similarly, the farm profitability as expressed by total national land rents varied between 1900 million Norwegian krone (NOK) for median yields under baseline climate conditions up to 3900 million NOK for median yield under future projected climate conditions.