Next Article in Journal
Towards the Fulfillment of a Knowledge Gap: Wood Densities for Species of the Subtropical Atlantic Forest
Previous Article in Journal
Semantic Earth Observation Data Cubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Data Descriptor

Correlations between Environmental Factors and Milk Production of Holstein Cows

Department of Biotechnology and Animal Health, Dnipro State Agrarian and Economics University, Dnipro 49100, Ukraine
*
Author to whom correspondence should be addressed.
Submission received: 5 July 2019 / Revised: 14 July 2019 / Accepted: 16 July 2019 / Published: 19 July 2019

Abstract

:
Global climate change is a challenge for dairy farming. In this regard, identifying reliable correlations between environmental parameters and animals’ physiological responses is a starting point for the mathematical modeling of their effects on the future welfare and milk production of cows. The aim of the study was to examine the relationship between environmental parameters and the milk production of cows in hot period. Archival data from the Ukrainian Hydrometeorological Center were used to study the state of insolation conditions (IC), wind direction (WD), wind strength (WS), air temperature (AT), and relative humidity (RH). The temperature–humidity index (THI) (Kibler, 1964) and temperature–humidity index in the hangar-type cowshed (THICHT) (Mylostyvyi et al., 2019) served as integral indicators of the state of the cowshed’s microclimate. The daily milk yield (DMY), yield of milk fat (MF) and milk protein (MP), and percentage of milk fat (PMF) and protein (PMP) were taken into account by the DairyComp 305 herd management system (VAS, USA). Statistical data processing was performed using the mathematical functions of Microsoft Excel (Microsoft Inc.) and Statistica 10 (StatSoft Inc.). There was a weak correlation between IC and DMY at r = −0.2, between RH and DMY at r = +0.4, and between RH and MF at r = +0.2. Between DMY, MF, MP, and WS made up r = –0.2 to 0.4. Between DMY, MF, MP, and AT made up r = −0.2 to 0.5 (p < 0.05). The effects of weather factors on animal productivity will be the subject of further research.
Dataset: Available in Supplementary Materials.
Dataset License: CC-BY

1. Introduction

High air temperatures adversely affect the welfare and productivity of dairy cows during the summer, not only in tropical regions, but on most of the European continent. This is the cause of significant losses of livestock [1,2]. However, high temperature should not be considered the only factor leading to decreased cow milk yield in the hot period of the year. Humidity, air velocity, and insolation can significantly affect animal physiology, acting together on their thermoregulation [3,4]. Special indices are used for comprehensive assessment of the effects of these environmental factors on animals. Of these, the most common is the temperature–humidity index (THI), which takes into account the effects of air temperature and relative humidity on animals [5,6,7]. Indices that take into account the mobility of air, the intensity of insolation, and other factors besides temperature and humidity were also proposed not so long ago [8,9,10]. As a rule, they are all built taking into account the close connection (correlation) between environmental factors and individual indicators of the functional state of cows, which act as specific predictors [11,12]. Milk yield and milk components can also be indicators of animal comfort, since they are affected by the physiological state of the organism subject to high temperature [13]. Many researchers have taken into account data of meteorological stations located near farms when assessing the impact of weather on animal productivity [14,15]. Such an approach is acceptable when evaluating the effect of heat on animals, not only during grazing but also in uninsulated barns in which the microclimate is similar to the state of the environment [16]. Despite the general availability of weather data (for example, on national weather sites), they contain a large array of records and their systematization before statistical computer processing takes a lot of time. Private information (in our case, the productivity of cows in a separate dairy complex) can be secret and inaccessible in most cases, which limits its use (verification) by other researchers. In this regard, the aim of the work was to study the correlation between environmental parameters and milk production indicators of Holstein cows during the warm period of the year. The results and data presented in this paper will be used for mathematical modeling of the influence of weather conditions on the milk productivity of cows during periods of heat in our future studies.

2. Data Description

Data on weather conditions and milk productivity of cows were collected during the warm season. The analysis was carried out on the basis of weather data provided by the Ukrainian Hydrometeorological Center in the warm months of the year (May through August 2017); they are publicly available on the Internet at Meteo.ua (https://meteo.ua/) in the weather archive section. The data include 2912 records (in the Supplementary Materials Data_1.xls file), and log indicators such as insolation conditions (IC), wind direction (WD), wind strength (WS), air temperature (AT), and relative humidity (RH). Temperature–humidity indices (THI and THICHT) were calculated by taking into account AT and RH using the appropriate formulas (described in the Methods section). These data have been added to the Supplementary Materials Data_2.xls file. The average values (for each day) for all indicators were calculated using the built-in mathematical formulas in Microsoft Excel for the four warm months. Accordingly, there were 123 entries for each weather indicator (IC, WD, WS, AT, RH, THI, and THICHT). These data are given in the Supplementary Materials Data_3.xls file.
Data on the productivity of cows included daily milk yield (DMY), milk fat (MF) and milk protein (MP), and percentage of milk fat (PMF) and protein (PMP), located in the Supplementary Materials Data_4.xls file.
Averages were taken for mathematical processing of all indicators (in the Supplementary Materials Data_5.xls file). These data contain 123 complete records of the mean values that were taken for the correlation analysis in this study.

3. Methods

The dairy complex is located in an open area (48°28′44″ N, 35°36′46″ E) near the city of Pavlograd. Cows are kept without restraint on the dairy complex. This large dairy complex is designed for many milk cows. All the animals are divided into technological groups (early, middle, and late lactation). The study was conducted on cows in middle lactation (91 to 210 days). The average number of months ranged from 748.4 to 772.5. The uninsulated barns, in the form of hangars, are equipped with canvas curtains that remain open during warm periods. Animals rest in four rows of cubicles with an area of 2.24 m2 per cow. The area per cow in the barn is 3.4 m2 (without a feeding alley). Sand is used as bedding. The cows are fed the same type of nutrient-balanced feed mixture all year round. The feed rationing is made up of silage and haylage from corn, straw from winter and spring crops, and concentrates. The composition and energy value of food depends on the physiological state and productivity of the cows. The animals have free access to water (group drinkers for animals) and feed (table for feeding). The climate features in the barn we described previously [17] are similar to the state of the environment. Large-diameter fans work indoors around the clock during the warm period. The average air velocity is 0.5–0.9 m/s at the animal resting place. Since the feeding of animals is the same throughout the year, this factor can be leveled to a certain extent.
The distance between the cowsheds and the weather station located in Pavlograd city does not exceed a straight 25 km. Geographically, this area belongs to the steppe of Ukraine, which is characterized by relatively constant weather conditions over a large area. Therefore, we considered the data taken from a nearby meteorological station to be quite relevant. Systematization of various weather data was carried out using a key (codes). Insolation conditions (IC) were characterized as follows: clear weather: 4; overcast prevails: 3; cloudy and raining: 2; cloudy constant: 1. This rating from the brightest illumination (4 points) to less bright (1 point) was chosen due to the lack of a quantitative characteristic of the intensity of solar radiation (in W/m2 or derived units) at the weather station. The wind was considered according to cardinal direction: northern: 1; northeastern: 2; eastern: 3; southeastern: 4; southern: 5; southwestern: 6; western: 7; northwestern: 8. Such coding may not have been entirely correct from a logical point of view (since the average data obtained could not be quantified), but we used this for mathematical data processing anyway. Wind speed (strength) was quantified by points [18]: 0–0.5 m/s: 0; 0.6–1.7 m/s: 1; 1.8–3.3 m/s: 2; 3.4–5.2 m/s: 3; 5.3–7.4 m/s: 4; 7.5–9.8 m/s: 5; 12.5–18.2 m/s: 6; ≥18.3 m/s: 7. If a necessary criterion with which to characterize wind strength (for example, wind speed of 10 m/s) was absent, it was assigned to the nearest value (in this example, 5). The cowsheds are located from north to south in relation to the Earth’s directions. In the warm season, the curtains of the uninsulated cowsheds are left lowered, which could additionally affect the air exchange in the rooms despite the operating fans. Therefore, we took into account the wind speed and direction outside the building. Air temperature (°C) and relative humidity (%) were recorded in corresponding values. Initially, the weather data was supposed to be recorded every hour to calculate the average values per day (if they were in the weather archive). Calculations of the temperature–humidity index and temperature–humidity index in the cowshed hangar were made according to the following equations:
THI = 1.8 × T − (1 − RH/100) × (T − 14.3) + 32
THICHT = 46.00549 + 1.04460 × T
where THI is temperature–humidity index [19], THICHT is temperature–humidity index in the hangar-type cowshed, T is ambient air temperature (°C), and RH is relative humidity (%). THICHT is an indicator we suggested previously [17], on the basis of repeated temperature and humidity measurements in an animal room throughout the year. This indicator was used because the animals in the study were kept in exactly a hangar-type barn. THICHT allows determination of the temperature–humidity index for only one indicator (external temperature), taking into account the design features of the cowshed. Data on cow productivity (milk yield per herd, number of dairy cows, daily milk yield, milk fat and protein yield, as well as fat percentage and protein in milk) were obtained in the conditions of the dairy complex by the DairyComp 305 herd management system (VAS, USA). The relationship between the state of the environment and the productivity of dairy cows was assessed using Statistica 10 software (StatSoft Inc.). Correlation analysis was done using Spearman’s rank correlation coefficient. Significance was determined to be p < 0.05.

4. Results

4.1. Weather Conditions

The weather conditions were studied from 1 May to 31 August 2017. We estimated the weather for a duration of 2912 h, which is about 98.6% of the warm period of the year, lasting 123 days (2952 h). It was found (Table 1) that there were 2255 h with clear weather (77.4%), it was mostly cloudy for 527 h (18.1%), it was rainy for 95 h (3.3%), and continuous cloudiness lasted 35 h (1.2%). In May, it was cloudy more often; in July, it was rainier; and the number of clear days was greatest in August.
It was established (Figure 1) that the prevailing winds were from the north (30.9%), the west (16.8%), the east (15.2%), and the northeast (15.0%) in the warm season. The southeast wind direction was not fixed in one of the cases.
As shown by month in a separate table (Table 2), the prevailing winds were from the north in May, June, and July. In August there was more wind from the east (the north and northeast also had a fairly large share).
The total amount of time that the wind strength (Table 3) was equal to 0 points was 221 h (7.6%); 1 point, 97 h (3.3%); 2 points, 1022 h (35.1%); 3 points, 956 h (32.8%); 4 points, 494 h (17.0%); and 5 points, 122 h (4.2%). Wind strength of 6 and 7 points was not recorded.
The windiest time was in August, when wind force from 3 to 5 points was 507 h, or 69% of the total time. However, the duration of “calm” was also longest this month, at 5.3%. Average values of air temperature and relative humidity are given as background indicators in the warm period of the year (Table 4).
However, the average values of the indicators of air temperature and relative humidity, as well as indices calculated taking them into account, do not reflect the duration of the hot period. Accordingly, the distribution of their values in time is given more specifically (Table 5).
These data show that THI and THICHT values were high enough for dairy cows to be uncomfortable [20], especially from June to August when the values increased tenfold.

4.2. Milk Production of Cows

Analysis of the milk production of cows was conducted from 1 May to 31 August 2017. The average daily milk yield and the content of its main components were taken into account for the herd (Table 6).
We assumed that May was the most comfortable for cows in terms of weather conditions. In the summer months, there was a noticeable decrease in the components of milk compared with May. We observed a decrease in milk fat yield by 17–47 g, milk protein yield by 9–40 g, milk fat content by 0.1–0.15%, and milk protein by 0.06–0.12%. Differences between May and July–August were significant (p < 0.05–0.001). We attribute such changes in the milk composition to the influence of environmental factors on animals, because all the cows had the same lactation period.

4.3. Correlations between the State of the Environment and Milk Production Indicators

Considering the relationship between environmental factors and the productive qualities of cows (Table 7), we noted a reliable correlation. The correlation was positive between IC and DMY (r = –0.2), as well as between RH and DMY (r = +0.4) and between RH and MF (r = +0.2). A negative correlation of different densities was obtained between DMY, MF, MP, and WS (r = –0.2 to 0.4). The milk production indicators correlated negatively with AT (r = –0.2 to 0.5) in all cases, and, for THICHT, the correlation between the productivity of cows was significantly negative with a similar strength.
Thus, weather conditions such as the intensity of solar radiation, wind strength, air temperature, and humidity should be taken into account when assessing the influence of the environment on the milk productivity of cows.

5. Conclusions

This paper demonstrates the possibility of using meteorological data available to a wide range of people in order to identify the relationship between the state of the environment and the productive qualities of dairy cows, as a specific example. Some approaches to summarizing individual meteorological data before conducting a correlation analysis will be useful for researchers working in toward this end. The dataset presented here can be reproduced by other researchers, which will serve as a basis for new ideas and interpretations.

6. User Notes

The key (codes) listed in the sections above will allow researchers to reproduce the results obtained in the paper (using Excel files in the public domain).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2306-5729/4/3/103/s1, Text S1: Weather and cow productivity data.

Author Contributions

R.M. performed methodology, software validation, and original draft preparation; O.C. performed writing—review and editing.

Funding

This research received no external funding.

Acknowledgments

The authors are grateful to Halyna Novokshonova for the materials submitted for the study, as well as to students of the biotechnology faculty, Anastasia Sokolan and Alyna Popovych, for collecting primary data.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Das, R.; Sailo, L.; Verma, N.; Bharti, P.; Saikia, J.; Imtiwati; Kumar, R. Impact of heat stress on health and performance of dairy animals: A review. Vet. World 2016, 9, 260–268. [Google Scholar] [CrossRef] [PubMed]
  2. Sejian, V.; Bhatta, R.; Gaughan, J.B.; Dunshea, F.R.; Lacetera, N. Review: Adaptation of animals to heat stress. Animal 2018, 12, 431–444. [Google Scholar] [CrossRef] [PubMed]
  3. Abdela, N.; Jilo, K. Impact of climate change on livestock health: A review. Global Vet. 2016, 16, 419–424. [Google Scholar]
  4. Laporta, J.; Fabris, T.F.; Skibiel, A.L.; Powell, J.L.; Hayen, M.J.; Horvath, K.; Miller-Cushon, E.K.; Dahl, G.E. In utero exposure to heat stress during late gestation has prolonged effects on the activity patterns and growth of dairy calves. J. Dairy Sci. 2017, 100, 2976–2984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bohmanova, J.; Misztal, I.; Cole, J.B. Temperature-humidity indices as indicators of milk production losses due to heat stress. J. Dairy Sci. 2007, 90, 1947–1956. [Google Scholar] [CrossRef] [PubMed]
  6. Kumar, P.; Upadyay, R.C.; Kumar, R.; Chaudhary, P.K.; Maurya, P.K.; Kumar, R.; Yadav, V.; Kumar, M. Evaluation and comparison of heat stress indices for cattle and buffaloes. Int. J. Curr. Microbiol. App. Sci. 2018, 7, 292–298. [Google Scholar]
  7. Mylostyvyi, R.V.; Sejian, V. Welfare of dairy cattle in conditions of global climate change. Theor. Appl. Vet. Med. 2019, 7, 47–55. [Google Scholar] [CrossRef]
  8. Mader, T.L.; Davis, M.S.; Brown-Brandl, T. Environmental factors influencing heat stress in feedlot cattle. J. Anim. Sci. 2006, 84, 712–719. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, X.; Gao, H.; Gebremedhin, K.G.; Bjerg, B.S.; Van Os, J.; Tucker, C.B.; Zhang, G. A predictive model of equivalent temperature index for dairy cattle (ETIC). J. Thermal. Biol. 2018, 76, 165–170. [Google Scholar] [CrossRef] [PubMed]
  10. Herbut, P.; Angrecka, S.; Walczak, J. Environmental parameters to assessing of heat stress in dairy cattle—A review. Int. J. Biometeorol. 2018, 62, 2089–2097. [Google Scholar] [CrossRef] [PubMed]
  11. Kim, W.S.; Lee, J.-S.; Jeon, S.W.; Peng, D.Q.; Kim, Y.S.; Bae, M.H.; Jo, Y.H.; Lee, H.G. Correlation between blood, physiological and behavioral parameters in beef calves under heat stress. Asian Australas J. Anim. Sci. 2018, 31, 919–925. [Google Scholar] [CrossRef] [PubMed]
  12. Daltro, D.S.; Fischer, V.; Alfonzo, E.P.M.; Dalcin, V.C.; Stumpf, M.T.; Kolling, G.J.; Silva, M.V.G.B.; McManus, C. Infrared thermography as a method for evaluating the heat tolerance in dairy cows. Rev. Bras. Zootec. 2017, 46, 374–383. [Google Scholar] [CrossRef] [Green Version]
  13. Behera, R.; Chakravarty, A.K.; Sahu, A.; Kashyap, N.; Rai, S.; Mandal, A. Identification of best temperature humidity index model for assessing impact of heat stress on milk constituent traits in Murrah buffaloes under subtropical climatic conditions of Northern India. Indian J. Anim. Res. 2017, 52, 13–19. [Google Scholar] [CrossRef] [Green Version]
  14. Beux, S.; Cassandro, M.; Nogueira, A.; Waszczynskyj, N. Effect of THI on milk coagulation properties of Holstein-Friesian dairy cattle. Rev. Bras. Zootec. 2017, 46, 429–432. [Google Scholar] [CrossRef] [Green Version]
  15. Kekana, T.W.; Nherera-Chokuda, F.V.; Muya, M.C.; Manyama, K.M.; Lehloenya, K.C. Milk production and blood metabolites of dairy cattle as influenced by thermal-humidity index. Trop. Anim. Health Prod. 2018, 50, 921–924. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Lambertz, C.; Sanker, C.; Gauly, M. Climatic effects on milk production traits and somatic cell score in lactating Holstein-Friesian cows in different housing systems. J. Dairy Sci. 2014, 97, 319–329. [Google Scholar] [CrossRef] [PubMed]
  17. Mylostyvyi, R.; Chernenko, O.; Lisna, A. Prediction of comfort for dairy cows, depending on the state of the environment and the type of barn. In Development of Modern Science: The Experience of European Countries and Prospects for Ukraine: Monograph, 3rd ed.; Baltija Publishing: Riga, Latvia, 2019; pp. 394–410. [Google Scholar] [CrossRef]
  18. Antonenko, P.P.; Dorovskych, A.V.; Vysokos, M.P.; Mylostyvyi, R.V.; Kalinichenko, O.O.; Vasilenko, T.O. Methodological Bases and Methods of Scientific Research in Veterinary Hygiene, Sanitary and Expertise; Svidler, A.L.: Dnipro, Ukraine, 2018; p. 270. [Google Scholar]
  19. Kibler, H.H. Thermal effects of various temperature-humidity combinations on Holstein cattle as measured by eight physiological responses. Environmental physiology and shelter engineering. Res. Bull. Missouri. Agric. Exp. Stn. 1964, 862, 1–42. [Google Scholar]
  20. Hammami, H.; Bormann, J.; M’hamdi, N.; Montaldo, H.H.; Gengler, N. Evaluation of heat stress effects on production traits and somatic cell score of Holsteins in a temperate environment. J. Dairy Sci. 2013, 96, 1844–1855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Directions of prevailing winds during the study period.
Figure 1. Directions of prevailing winds during the study period.
Data 04 00103 g001
Table 1. Characteristics of weather conditions.
Table 1. Characteristics of weather conditions.
MonthClear WeatherOvercast PrevailsCloudy and RainingCloudy Constant
Hours%Hours%Hours%Hours%
May52971.816522.4273.6162.2
June56378.713418.7182.6
July53473.713819.0456.281.1
August62985.69012.250.7111.5
Table 2. Direction of the wind.
Table 2. Direction of the wind.
Wind DirectionMayJuneJulyAugust
Hours%Hours%Hours%Hours%
Northern23331.621029.427537.918224.8
Northeastern14820.1537.4729.916522.4
Eastern9212.5638.8709.721729.5
Southeastern
Southern415.6395.4354.8
Southwestern689.2659.1385.2121.6
Western10414.114720.613318.410414.2
Northwestern516.913819.310214.1557.5
Table 3. Wind strength.
Table 3. Wind strength.
PointsMayJuneJulyAugust
Hours%Hours%Hours%Hours%
0567.6456.38111.2395.3
1344.6294.0263.681.1
226936.527838.929440.618124.6
323331.623833.320928.827637.6
411115.110214.39813.518324.9
5344.6233.2172.3486.5
Table 4. Average temperature and relative humidity condition of air, mean ± standard deviation (SD). THI, temperature–humidity index; THICHT, temperature–humidity index in hangar-type cowshed.
Table 4. Average temperature and relative humidity condition of air, mean ± standard deviation (SD). THI, temperature–humidity index; THICHT, temperature–humidity index in hangar-type cowshed.
Month Air TemperatureRelative HumidityTHITHICHT
nMeanSDMeanSDMeanSDMeanSD
May3115.60.6459.92.4358.80.8762.30.67
June3020.90.5958.31.7966.30.8067.80.62
July3121.60.6263.91.8367.40.7968.50.65
August3124.20.8250.62.6669.91.0071.50.86
Table 5. Distribution of temperature and humidity index values (THI and THICHT) over time in the warm months of the year (h).
Table 5. Distribution of temperature and humidity index values (THI and THICHT) over time in the warm months of the year (h).
MonthTHITHICHT
<6868.0–71.972.0–79.980.0–89.9<6868.0–71.972.0–79.980.0–89.9
May64973156068744
June42915213314181231686
July386146188539113118419
August28814824554253156206120
Note: The gradation of temperature and humidity index values was made as follows: values below 68 correspond to comfortable conditions for cows, 68–71 corresponds to slight stress, 72–79 to moderate stress, and 80–89 to strong stress [17].
Table 6. Milk productivity of cows, mean ± SD.
Table 6. Milk productivity of cows, mean ± SD.
MonthDMY (kg)MF (kg)MP (kg)PMF (%)PMP (%)
MeanSDMeanSDMeanSDMeanSDMeanSD
May23.70.030.8510.0050.7650.0043.600.0203.230.018
June23.80.020.8340.0040.7560.0053.500.0143.170.020
July24.00.050.8330.0030.7460.0063.470.0103.110.023
August23.30.050.8040.0050.7250.0063.450.0183.110.022
DMY, daily milk yield; MF, yield of milk fat; MP yield of milk protein; PMF, percentage of milk fat; PMP percentage of milk protein.
Table 7. Correlations (r) between the state of the environment and milk production of cows.
Table 7. Correlations (r) between the state of the environment and milk production of cows.
State of the EnvironmentDMYMFMPPMFPMP
IC−0.195 *−0.168−0.110−0.095−0.046
WD+0.090+0.015−0.006−0.023−0.035
WS−0.409 *−0.206 *−0.189 *−0.046−0.059
AT−0.186 *−0.465 *−0.366 *−0.443 *−0.333 *
RH+0.399 *+0.225 *+0.126+0.066−0.009
THI−0.113−0.447 *−0.354 *−0.457 *−0.345 *
THICHT−0.187 *−0.466 *−0.367 *−0.444 *−0.333 *
* Reliable values are shown (p < 0.05).

Share and Cite

MDPI and ACS Style

Mylostyvyi, R.; Chernenko, O. Correlations between Environmental Factors and Milk Production of Holstein Cows. Data 2019, 4, 103. https://0-doi-org.brum.beds.ac.uk/10.3390/data4030103

AMA Style

Mylostyvyi R, Chernenko O. Correlations between Environmental Factors and Milk Production of Holstein Cows. Data. 2019; 4(3):103. https://0-doi-org.brum.beds.ac.uk/10.3390/data4030103

Chicago/Turabian Style

Mylostyvyi, Roman, and Olexandr Chernenko. 2019. "Correlations between Environmental Factors and Milk Production of Holstein Cows" Data 4, no. 3: 103. https://0-doi-org.brum.beds.ac.uk/10.3390/data4030103

Article Metrics

Back to TopTop