INTRODUCTION
Heat stress has been recognised as a serious problem around the world, but especially on dairy farms using breeds with a high potential of milk production (Arias et al 2008, Silanikove 2000). Usually, these animals fail to maintain an adequate thermal balance causing negative consequences for milk production and animal welfare (Collier et al 2006, Kadzere et al 2002). The poor performance and efficiency in milk yield have been associated mainly with a decrease in dry matter intake (West 2003), resulting in significant economic losses (St-Pierre et al 2003), estimated to be $897 million only for the dairy industry at the USA. A summary of the major behavioural responses includes changes on feeding, defecating and urinating frequency, water intake, lying time, standing time, shade seeking behaviour (Ratnakaran et al 2017). In addition, body temperature as well behaviour has been affected by heat stress (Jara et al 2016, Tucker et al 2008).
In Chile, there is limited information regarding the effect of heat stress on dairy cattle. However, there is a risk of heat stress for the central region of Chile, but with a moderate to low risk in the southern regions (Arias and Mader 2010). The aforementioned is relevant considering that dairy farms located in the central region of Chile have a greater degree of intensification (total mixed rations and high production genotypes) and higher milk yield per cow than those located in the southern regions, where the systems are based on pasture. On the other hand, this subject is acquiring major attention and concern among producers due the changes in animal welfare requirements as well as changes observed on climate, with an increasing proportion of heat waves and extended periods of drought. We hypothesised that cows of the central region of Chile experience physiological and behavioural changes as consequence of the summer weather conditions. Therefore, the objectives of the study were to assess in two dairy genotypes changes on patterns of tympanic temperature, panting scores, respiration rates and utilisation of artificial shade during summer conditions in the central region of Chile.
MATERIAL AND METHODS
The study was conducted during the summer of 2011 at the dairy research farm of the Pontificia Universidad Católica de Chile, located in the foothills of the Andes mountains of the central region of Chile (33°40’22.4”S 70°35’33.8”W, 654 meters above sea level). The study consisted of three periods of data collection: January 19 to 25th; February 5 to 8th; and March 20 to 24th. In each period, 6 Holstein (H) and 6 Holstein x Montbeliarde (HM) cows were selected to receive a data logger device to collect tympanic temperatures (TT) and to study its behaviour. Thus, a total of 36 mature multiparous and healthy cows (67 ± 6 d on milk), body condition score = 2.75 to 3.00 (1 to 5 scale) at the beginning of the experiment, were randomly selected from the commercial herd. All the cows were fed with the same diet three times a day (08:00 h, 11:30 h, and 18:00 h), which included (all expressed as DM/cow/ day) corn silage (4.0 kg), alfalfa soiling (4.9 kg), barley brewers (2.9 kg), and supplements (0.25 kg).
Animals were kept in open pens (80 m x 80 m) with ad libitum access to water and shade, and milked three times per day (03:00, 12:00 and 18:00 h). However, milk yield was not considered in this analysis. The shading structure consisted of treated wood poles of approximately 4.0 m high, with a black raschel mesh (80% shade, 2.1 m wide) located over the wood poles, providing 2.8 m2 of shade by cow.
TYMPANIC TEMPERATURE AND ANIMAL BEHAVIOR
Each cow within each period received an iButton data logger device (Maxim Integrated Products Inc., CA, USA) located manually in its tympanic canal. The devices were programmed to collect TT at 10 minute intervals and sub sequently compiled into hourly readings. The same cows used to collect TT were observed within each period to collect behavioural data. One trained observer recorded the proportion of cows using shade on the pen and also collected data of respiration rates (RR; breaths per minute; bpm). The RR were estimated by timing and counting 10 flank movements in each animal. Measurements of the panting score (PS) were also collected following the description of Mader and Davis (2002). First, the observer identified cows with the ibutton device within the herd, recording its location (shaded or unshaded), then recorded the count of flank movement and PS score independent of cow’s location in the pen.
ENVIRONMENTAL DATA COLLECTION AND THERMAL COMFORT INDICES
Ambient temperature (AT, °C), wind speed (WS, m/s), relative humidity (RH,%), solar radiation (SR, W/m2), and precipitation (P, mm/d)) were collected continuously at 10 minute intervals by using a weather station (U30 Hobo Onset, MA, USA) located at the dairy farm. Later, these data were compiled into hourly values to match TT dataset. Likewise, these climatic data were used to estimate three thermal comfort indices: temperature-humidity index (THI), adjusted THI (THIadj) and comprehensive climate index (CCI) by using the following equations:
Where:
FRH corresponds to the correction factor for AT due to relative humidity;
FWS corresponds to the correction factor for AT due to wind speed;
and FSR corresponds to the correction factor for AT due to solar radiation.
The risk of thermal stress is given by the following categories of CCI (Mader et al 2010): No stress (CCI ≤ 25); Mild (> 25 and ≤ 30); Moderate (> 30 and ≤ 35); Severe (> 35 and ≤ 40); Extreme (> 40 and ≤ 45); and Extreme danger (CCI > 45).
STATISTICAL ANALYSIS
During the study, eight data logger devices were lost (one per period for H genotype and 4 and 1 for HM genotype in periods 2 and 3, respectively). In addition, data from 2 animals were removed from data set due to incomplete records of TT (less than 24 h/day or less than 4 days/ period), resulting in an unbalanced number of animals per period (total n = 26). Meteorological and TT data were analysed by using a complete randomised experimental design arrangement, with each animal as experimental and observational unit. Additionally, a dummy variable (CCI condition) was created based on threshold of CCI defined by Mader et al (2010). Thus, hour of each day with CCI ≤ 25 °C were considered as “Normal”, otherwise it was considered as “Stressful”. Thus, genotype, moment of day (AM vs. PM), and CCI condition (Normal vs. Stressful) were considered like study factors in a factorial arrangement. In addition, TT was modeled using a repeated measurements analysis by using the MIXED procedure of SAS 9.4 (SAS Institute, Cary, NC) with TT as the de pendent variable with genotype, hour and its interaction as independent variables in the model. The hour was the repeated measurement and the animal (genotype) statement was used as random effect. The period was not included in the model because different animals were used across the periods. Finally, categorical data were analysed by using a Chi square test (Likelihood ratio test). The level of significance for all the statistical analyses was 0.05.
RESULTS
A summary of meteorological variables, thermal indices, and TT by period, type and moment of day are presented on table 1. Most of climatic variables showed a decrease on periods 2 and 3 regarding period 1, with the exception of HR. In addition, even though AT reached values >30 °C during all days in periods 1 and 2, there was a considerable drop during night-time (mean 8.7 °C), representing over 20 points of fluctuation (Figure 1).
TT = Tympanic temperature (genotypes pooled); AT= Ambient temperature; RH= Relative Humidity; WS= Wind speed; SR= Solar radiation; THI= Temperature humidity index; THIadj= Adjusted temperature humidity index; CCI = Comprehensive Climate Index.
Different letters within columns indicate significant differences between periods (P<0.05). In addition, all variables showed significant differences when compared Normal vs. Stressful and AM vs. PM (P<0.001).
Period 1 = January 19 to 25; Period 2 = February 05 to 08; and Period 3 = March 20 to 24 of 2011.
Figure 2 shows interaction period x CCI condition x genotype (P=0.003). Cows of both genotypes showed the highest TT during stressful hours within each period. There were also interactions for genotype x CCI condition (P<0.001), with HM cows showing greater TT than H cows in both CCI conditions; and for genotype x period (P<0.001), with higher TT on period 1 in both genotypes. Finally, there was also an interaction of genotype x hour on TT (P<0.001; Figure 3), with H cows showing lower TT during great part of the day (from 21:00 to 11:00; P<0.100).
Over 75% of cows were under shade during daytime. However, no differences were observed for morning vs. afternoon (P=0.543) for pooled data. Nevertheless, the proportion of cows using shade increased when the hourly CCI category moved from “Normal” to “Mild “ or “Moderate” stress (P=0.049), but without differences between the last two categories. Similarly, this proportion also was higher (P=0.007) during stressful condition when compared with normal condition (89.0 vs. 75.3%). As expected, the mean value of CCI for a stressful condition (27.96 ± 0.21 °C) was greater than the normal condition (21.40 ± 0.28 °C; P<0.001). Cows showed a trend to reduce the use of shade during the third period of study (P=0.094), coinciding with the end of the summer season and a lower SR.
The proportion of cows under shade across the periods was slightly higher during the morning vs. afternoon (86% AM vs. 79% PM). By the contrary, cows had higher RR during the afternoon when compared with the morning (67.51 ± 1.13 vs. 56.5 ± 0.99 bpm; P<0.001). Likewise, RR increased (P<0.001) from 59.98 ± 1.00 (Normal CCI condition) to 68.11 ± 1.35 bpm (Stressful CCI condition). When RR was compared among the hourly CCI categories showed an increase (P<0.001) from 60.16 ± 0.95 bpm (Normal) to 69.13 ± 2.73 and 67.76 ± 1.58 for “Mild” and “Moderate” CCI categories, but without differences between them (P>0.100). The PS changed across the day with a greater proportion of cows showing a PS=0 during the morning (57.4%), but decreasing during the afternoon (26.4%), whereas cows on PS=1 increased from 41.8% in the morning to 64.8% in the afternoon. The proportion of cows on PS=2 increased more than forty folds, from 0.2% to 8.8% when comparing morning vs. afternoon. A similar pattern was observed when PS of cows was compared across periods, but no effects of shade on PS were observed (P=0.765).
DISCUSSION
Heat stress has tremendous effect on production, health, and reproduction in lactating dairy cows, and has been largely ignored from a management standpoint (Tao and Dahl 2013). However, in many countries there is a lack of information about the impacts of heat stress on animal production and wellbeing. The present study is one of the firsts in Chile addressing this topic.
A possible explanation to the difference observed on TT between the genotypes (Figure 3) could be the better adaptation of the H cows to the environment when compa red with the HM cows, that represent a newer genotype in the area of study. The higher TT during nighttime of HM cows could be the resultant of a lower capacity to dissipate the heat accumulated during daytime, as demonstrated by Aharoni et al (2006) who compared H vs. HM dairy cows during the summer season in Jordan. These authors reported differences between H and HM in their diurnal patterns of heat production, suggesting that heat tolerance of HM was lower than H cows. In our study, we did not measure the heat production.
The difference of TT between the morning vs. after noon is in agreement with those reported by Vickers et al (2010), who measured vaginal temperatures of dairy cows. But differs from those reported by Aharoni et al (2005) who found no differences for rectal temperatures and RR. In the ruminants, rectal temperatures are considerably affected by rumen metabolism and do not represent a re liable index of the regulated temperature (Berman 1971). In fact, Burfeind et al (2010) concluded that measures of rectal temperature could be influenced by the procedure itself, type of thermometer, and the penetration depth into the rectum. Thus, some care is required when rectal temperatures are used as estimator of body temperature.
Although values of TT on “stressful” conditions were higher than under the “normal” condition, these are still within the normal range (Hillman 2009). We speculate that the cold nights observed during all the periods of data collection provide to the animals the opportunity to dissipate the heat accumulated during the daytime. In this context, some researchers have concluded that cattle that fail to or that do not cool down at night are prone to achieving greater body temperatures during hot days. Cattle that are prone to getting hot but can cool at night can keep peak body temperatures at or near those of cattle that tend to consistently maintain lower body temperatures (Mader and Johnson 2010, Mader et al 2010). The big drop in AT during night time, typical of the Andes foothills in that region (≥ 20 °C of fluctuation), apparently allowed the cows to maintain TT within a normal range and to cope in better way the impact of heat load.
In this study, THI and THIadj ranged from 57.0 to 67.3 and 62.9 to 73.8 respectively. These values that are under the threshold at which milk production (du Preez et al 1990, Johnson 1985) and dry matter intake (Johnson 1985) begin to decrease (THI≥72). Nevertheless, Zimbelman et al (2009) proposed a new THI threshold of 68 for high producing dairy cows (> 35 L/day). The previous is in agreement with those reported by Markovich (2012) in New Zealand, where milk production began to decrease at THI of 65 for H and 75 for Jerseys. Bernabucci et al (2014) also reported a reduction in milk production at THI from 65 to 76 depending of the parity of cows. There are other reports showing a range of multiple values of THI thresholds ranging from 69 to 78.2 (Bohmanova et al 2007, Bouraoui et al 2002, Dikmen and Hansen 2009, Johnson et al 1962, Ravagnolo et al 2000). However, it is important to mention that most of those reports are associated with the same genetic (Holstein Friesian). Nevertheless, other breeds are less sensitive to thermal stress than Holstein Friesian (Sharma et al 1983). Thus, differences in response to heat stress between genotypes can be attributed to varying levels of adaptability to hot environments (Kadzere et al 2002).
There is strong evidence that lactating cows increase RR during heat stress events or even they are able to reduce RR when they are properly cooled (Hillman et al 2001, Smith et al 2006). In our study, cows had RR slightly higher than those reported by Kendall et al (2007) for the shade treatment in a study conducted in New Zealand, with similar weather conditions, but it was lower than the uncooled treatment. The differences of RR between stressful vs. normal conditions herein reported are in agreement with those reported by Muller et al (1994), who found differences in RR for cows in hot days (AT ≥ 25.1 °C), but not for cool days (AT ≤ 25 °C). In our study, RR were also influenced by moment of day (AM vs. PM). The differences in RR herein presented as well as those reported by Bouraoui et al (2002), and by Muller et al (1994) suggest that cows were primarily subject to heat stress during those hours of the day when AT and SR reach maximum values, resulting in significant increases of RR. Thus, SR and AT are two important drivers and triggers of RR (Berman et al 1985, Harris et al 1960). The RR can be used as a valuable tool to assess the level of heat stress on cows. Finally, the high proportion of cows using shade, regardless of the moment of day, implies that cows uses this mechanism to cope with heat load of daytime. In addition, when data were analysed by period, there was a trend to seek shade during the afternoons (P=0.073) in period 1 (January) which coincides with the highest values of SR. Based on the results of this study we can indicate that cows of both genotypes showed a slight degree of heat stress in Central region of Chile, particularly during daytime. However, they were able to cope with it because of the strong drop on AT during nighttime as well as by the using of shade.