Negative association between high temperature-humidity index and milk performance and quality in Korean dairy system: big data analysis

The aim of this study was to investigate the effects of heat stress on milk traits in South Korea using comprehensive data (dairy production and climate). The dataset for this study comprised 1,498,232 test-day records for milk yield, fat- and protein-corrected milk, fat yield, protein yield, milk urea nitrogen (MUN), and somatic cell score (SCS) from 215,276 Holstein cows (primiparous: n = 122,087; multiparous: n = 93,189) in 2,419 South Korean dairy herds. Data were collected from July 2017 to April 2020 through the Dairy Cattle Improvement Program, and merged with meteorological data from 600 automatic weather stations through the Korea Meteorological Administration. The segmented regression model was used to estimate the effects of the temperature-humidity index (THI) on milk traits and elucidate the break point (BP) of the THI. To acquire the least-squares mean of milk traits, the generalized linear model was applied using fixed effects (region, calving year, calving month, parity, days in milk, and THI). For all parameters, the BP of THI was observed; in particular, milk production parameters dramatically decreased after a specific BP of THI (p < 0.05). In contrast, MUN and SCS drastically increased when THI exceeded BP in all cows (p < 0.05) and primiparous cows (p < 0.05), respectively. Dairy cows in South Korea exhibited negative effects on milk traits (decrease in milk performance, increase in MUN, and SCS) when the THI exceeded 70; therefore, detailed feeding management is required to prevent heat stress in dairy cows.