Next Article in Journal
A Nonlinear Double Model for Multisensor-Integrated Navigation Using the Federated EKF Algorithm for Small UAVs
Next Article in Special Issue
Classification of Smoke Contaminated Cabernet Sauvignon Berries and Leaves Based on Chemical Fingerprinting and Machine Learning Algorithms
Previous Article in Journal
Exploring the Effects of Various Polymeric Backbones on the Performance of a Hydroxyaromatic 1,2,3-Triazole Anion Sensor
Previous Article in Special Issue
A Machine Learning-based Algorithm for Water Network Contamination Source Localization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Artificial Intelligence Applied to a Robotic Dairy Farm to Model Milk Productivity and Quality based on Cow Data and Daily Environmental Parameters

1
Digital Agriculture, Food, and Wine Group, Faculty of Veterinary and Agricultural Sciences, The University of Melbourne, Parkville, VIC 3010, Australia
2
Agricultural Production System Modelling Group, Faculty of Veterinary and Agricultural Sciences, The University of Melbourne, Parkville, VIC 3010, Australia
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(10), 2975; https://doi.org/10.3390/s20102975
Submission received: 27 April 2020 / Revised: 8 May 2020 / Accepted: 22 May 2020 / Published: 24 May 2020

Abstract

:
Increased global temperatures and climatic anomalies, such as heatwaves, as a product of climate change, are impacting the heat stress levels of farm animals. These impacts could have detrimental effects on the milk quality and productivity of dairy cows. This research used four years of data from a robotic dairy farm from 36 cows with similar heat tolerance (Model 1), and all 312 cows from the farm (Model 2). These data consisted of programmed concentrate feed and weight combined with weather parameters to develop supervised machine learning fitting models to predict milk yield, fat and protein content, and actual cow concentrate feed intake. Results showed highly accurate models, which were developed for cows with a similar genetic heat tolerance (Model 1: n = 116, 456; R = 0.87; slope = 0.76) and for all cows (Model 2: n = 665, 836; R = 0.86; slope = 0.74). Furthermore, an artificial intelligence (AI) system was proposed to increase or maintain a targeted level of milk quality by reducing heat stress that could be applied to a conventional dairy farm with minimal technology addition.

1. Introduction

Robotic dairy farms or Automated Milking Systems (AMS) are the result of the implementation of state of the art technology related to robotics to increase milk yield through increased efficiency and automation [1,2]. These technologies are developed in response to the increasing market opportunities for the dairy industry globally, which is projected to grow by 35% by 2030 [3]. However, global demands will also be accompanied by 14 million traditional dairy farms shutting down production due to increased competitiveness and requirements for guaranteed milk quality and animal welfare [4]. The latter is considered a growing concern for consumers, which is achieved by AMS since it is based on the “milking when they like” system increasing wellbeing and welfare of cows [5]. Further potential advances to AMS technologies have been researched in recent years through the implementation of biometrics monitoring of animals to assess physiological changes in production systems [6]. Some of these technologies are noninvasive using visible (RGB) imagery/video, and infrared thermal imagery for heart rate, respiration rate, and body temperature assessments. These technologies could result in improvements in the monitoring of heat stress in farm animals.
Modeling heat stress in AMS has concentrated recently on the rumination and milking performance [7], identifying specific thresholds with production factors [8] and thermal comfort indices [9], mainly through the calculation of the temperature-humidity index (THI) using several models [10]. According to a study by Nascimento et al. [11], who compared nine different models to calculate THI, the equation from Berman et al. [12] (Equation (10) below) was significantly correlated with physiological data of cows such as respiration rate, heart rate, rectal, and skin temperatures. However, all previous methods use deterministic mathematical equations with minimal animal information in the analysis, and the noncontact biometric analysis could be cost-prohibitive for the near-future application to conventional dairy farms.
Artificial intelligence (AI) applied to Digital Agriculture deals with the implementation and integration of digital data, sensors, and tools on agricultural applications from the farm to consumers [13]. These technologies can include big data, sensor technology, sensor networks, remote sensing, robotics, and unmanned aerial vehicles (UAV). Data processing is performed using new and emerging technologies, such as computer vision, machine learning, and AI, among others. The implementation of AI not only should benefit high technological systems, such as AMS, but also conventional dairy farms to increase their competitiveness in the future.
This research was based on machine learning modeling using ubiquitous environmental data obtained from automatic meteorological stations and cow information available by all dairy farms as inputs. Target information related to important parameters related to milk productivity, milk quality, and actual feed of dairy cows was obtained from an AMS belonging to The University of Melbourne, Australia. High accurate machine learning (ML) models that can be applied to any dairy farm from AMS to conventional were obtained. Furthermore, this paper proposes an AI system model to be implemented in any dairy farm to automatically assess and ameliorate heat stress by implementing ML models developed with an automated sorting and gate system.

2. Materials and Methods

2.1. Site, Robotic Dairy Farm, and Data Acquisition

The study was conducted in a dairy farm located at The University of Melbourne Dookie College, Victoria, Australia (36°22’48” S, 145°42’36” E). This region had an average annual rainfall of 537 mm (monthly extremes: 30.5–57.6 mm) and mean daily solar exposure of 17 MJ m2 −1 (extremes: 7.3–27.3 m2 −1) from 1991–2019; data obtained from the Bureau of Meteorology (BoM) Dookie Agricultural College station 081013. The farm consists of 43 ha of irrigated pastures based on perennial ryegrass (Lolium perenne) and annual ryegrass (Lolium multiflorum). The herd in this site consists of Holstein-Friesian cows. The farm contains three Lely Astronaut robotic milking machines (Lely Holding S.à.r.l., Maassluis, The Netherlands), with a capacity of 60 cows per machine (maximum capacity of 180 cows) that move voluntarily for milking. As described by Dunshea et al. [14], cows wear an identification transponder neck collar (Lely Holding S.à.r.l., Maassluis, The Netherlands), which records the cows’ activity. The robotic milking system can automatically record parameters such as lactation days counted from day 0 at calving up to the time of next calving including the dry cow period, lactation number, milking frequency per day, milk yield (kg day−1), milk protein (%), milk fat (%) and somatic cells, programmed concentrate feed (kg day−1), concentrate feed intake (kg day−1), and liveweight (kg). Records of these data from June 2016 to March 2019 were used for this study.
Weather data were obtained from the meteorological station (Adcon Telemetry GmbH, Klosterneuburg, Austria), located at the Dookie Agricultural College, which provides data every 15 min for each day of the year. Parameters obtained were (i) temperature (T; °C), (ii) relative humidity (RH; %), (iii) rainfall as daily running total (mm), (iv) wind speed (km h−1), and (v) wind direction (°). Based on these data, other variables such as dewpoint temperature (Tdp; °C; Equation (1); [15]), wet bulb temperature (Twet; °C), and THI were calculated. The latter was calculated using the following nine different equations (Equations (2)–(10); [11]):
T d p = 243.5   ( 17.67 × T 243.5 + T + l n R H 100 ) 17.67 ( 17.67 × T 243.5 + T + l n R H 100 )
T H I 1 = 0.4 × ( T + T w e t ) × 1.8 + 32 + 15
T H I 2 = ( 0.15 × T + 0.85 × T w e t ) × 1.8 + 32
T H I 3 = ( T × 0.35 + T w e t × 0.65 ) × 1.8 + 32
T H I 4 = 0.72 × ( T + T w e t ) + 40.6
T H I 5 = ( 1.8 × T + 32 ) [ ( 0.55 0.0055 × R H ) × ( 1.8 + T 26 ) ]
T H I 6 = ( 0.55 × T + 0.2 × T d p ) × 1.8 + 32 + 17.5
T H I 7 = T + ( 0.36 × T d p ) + 41.2
T H I 8 = ( 0.8 × T ) + ( R H 100 ) × ( T 14.4 ) + 46.4
T H I 9 = 3.43 + 1.058 × T 0.293 × R H + 0.0164 × T × R H + 35.7
where Twet was calculated in batch using a customized code written in MATLAB® R2020a (Mathworks Inc., Natick, MA, USA; [16]), calculations were based on T, Tdp, and surface pressure and the bisection search method.

2.2. Statistical Data and Machine Learning Modeling

Mean values of THI calculated with Equation (10) along with milk yield, milk protein, and fat content, and concentrate feed intake were obtained and plotted to visualize the effects of the different seasons on each parameter. Statistical data obtained from the inputs and targets consisted of minimum, maximum, and mean values of each parameter.
Two ML models were developed based on artificial neural networks (ANN) using the Bayesian Regularization training algorithm. The latter was chosen as it showed the best accuracy and performance as well as no over or underfitting [17] after testing 17 different algorithms using a customized code written in MATLAB® R2020a. The inputs for the models (Figure 1) were based on the maximum values per day of the weather data (i) T, (ii) RH, (iii) rainfall, (iv) wind speed, (v) wind direction, (vi) Tdp, (vii) Twet, (viii–xvi) THI calculated with the nine equations, and some data obtained from the robotic milking system, (xvii) programmed concentrate feed, (xviii) lactation days, (xix) lactation number, (xx) milking frequency, and (xxi) liveweight. The targets were also obtained from the robotic milking system. They consisted of (i) milk yield, (ii) milk protein, (iii) milk fat, and (iv) concentrate feed intake (i.e., cereal grain-based pellets fed to cows during milking, making up approximately 40% of cows diet). All data were normalized from −1 to 1. Model 1 was constructed using the data of cows with a similar heat tolerance (N = 36; heat tolerance range: 93–112) determined by estimation of Australian genomic breeding values for heat tolerance [18] following genotyping of each cow using hair follicle samples as per the commercial procedure (CLARIFIDE for dairy, Zoetis Australia Pty Ltd, Banyo, QLD, Australia). The genotyping experiment was approved by the University of Melbourne Faculty of Veterinary and Agricultural Science (FVAS) Animal Ethics Committee (AEC ID 1814645.1). In general, for heat stress, cows with Australian breeding values < 100 are less tolerant to hot, humid conditions than the average, while the cows with values > 100 are more tolerant than the average. Specifically, cows with breeding values of 93 will be 7% less heat tolerant than an average cow, and a cow with heat tolerance breeding values of 110 would be 10% more heat tolerant as compared to an average cow. In contrast, Model 2 was developed using data from all cows (N = 312) independent of their heat tolerance to create a general model. Samples were divided randomly as 70% for training and 30% for testing using a default derivative function. Ten neurons were chosen as the best number giving the highest accuracy and best performance based on the means squared error (MSE).

3. Results

Figure 2 shows the mean values per season of each year for THI9 and the four parameters used as targets in the ML models to represent the effect of different weather patterns on those variables. As expected, the highest THI were obtained in the summer seasons of all years (77.5–79.7) and the lowest in winter of all years (47.6–49.1). The highest average milk yield per cow was observed in winter (33.4 kg day−1) and spring 2017 (33.5 kg day−1) with the lowest yield in summer 2018–2019 (23.4 kg day−1). The latter season also presented the lowest protein content in milk (3.1%) and concentrate feed intake (4.3 kg day−1). Spring 2018 and autumn 2018 had the lowest (3.9%) and the highest milk fat content (4.6%), respectively.
Table 1 shows the minimum, maximum, and mean values per year of each parameter used as inputs to construct the ML models. The lowest mean temperature (19.3 °C) was observed during 2016, which, at the same time, presented the lowest mean THI1THI9 (58.1–72.0), highest mean RH (95.6%), and daily rainfall (3.9 mm). On the contrary, 2019 had the highest maximum temperature (44.9 °C) and, until March, the lowest mean RH (69.2%), and daily rainfall (0.3 mm), as well as the highest mean THI1THI9 (68.6–82.8). Data for lactation days = 0 are the day the calf was born, and milk production commenced. Due to the voluntary milking system on the farm, there are some days when cows are not milked (i.e., milking frequency = 0). Furthermore, there were cows on the farm with extended lactations (>600 days). These were ‘carryover’ cows that were in an extended lactation because they failed to get pregnant in a timely manner.
Table 2 shows the minimum, maximum, and mean values of the parameters used as targets for the ML models. It can be observed that 2017 presented the highest milk yield per cow on average (30.7 kg day−1), although 2016 had the highest maximum milk yield per cow (65.4 kg day−1). Likewise, for milk protein, 2017 had the highest maximum value (6.1%), while 2018 presented the highest mean value (3.4%). Regarding milk fat content, 2019 had the highest maximum and mean values (10.9% and 4.3%, respectively). In 2019, the lowest average concentrate feed intake (4.0 kg day−1) was observed, while 2017 presented the highest mean (7.4 kg day−1) and the highest maximum value (24.3 kg day−1).
Table 3 shows the statistical results of both models to predict milk yield, milk fat, and protein content, and concentrate feed intake. It can be observed that both models presented similar results with high overall correlation coefficients (Model 1: R = 0.87; Model 2: 0.86; Figure 3). None of the models showed any signs of overfitting as the correlation coefficient of all stages was the same, and the performance of training (Model 1: MSE = 0.0186; Model 2: MSE = 0.0154) was lower than the testing stage (Model 1: MSE = 0.0189; Model 2: MSE = 0.0157). According to the 95% confidence bounds, Model 1 presented 3.88% outliers (4513 out of 116,456) and Model 2 presented 3.60% (23,998 out of 665,836).

4. Discussion

4.1. Seasonality and Milk Yield

During the four years included in this study (2016–2019), there was a clear variation within seasons reflected by environmental parameters (THI) and milk productivity parameters (Figure 2). Higher heat stress risks for cows were observed in the summer of 2018–2019. Even though the THI parameter had a higher tendency, it was not significantly greater compared to the THI of summers belonging to 2017 and 2016 (THI = 79.7 compared to 78.7 and 77.5, respectively). However, milk yield and quality parameters were lower for 2018 compared with previous years. The high variability among all parameters shown through the years considered for this study can be considered as an advantage for ML modeling. These differences can be further supported by the data presented in Table 1 and Table 2 with more specific data per year. Prolonged periods of high temperature and relative humidity have shown to be detrimental to dairy cows performance due to heat stress [19]. This makes more critical the development of cost-effective methodologies to measure and alleviate heat stress during these periods of high THI [20].

4.2. Machine Learning Models

By investigating thermotolerance in cows from a genetic point of view, it could help to decrease economic losses associated with lower milk productivity, quality, and animal welfare [21,22]. Other methods have been based on the physical modification of the environment, such as shade and shelters, and dietary interventions to reduce heat stress effects, such as grape residue [23], açai [24], betaine [14,25], slowly fermentable grains [26], and other types of feed [27,28].
The ML models developed in this research (Model 1 and Model 2) do not differ much when considering 36 genetically similar cows for heat tolerance compared to a total of 320 cows. There is a slight difference in the slope for the general model considering all cows (Model 2; slope = 0.74) compared to Model 1 (slope = 0.76). Considering highly heat stress-tolerant cows helps to decrease underestimations made by Model 1 compared to Model 2. However, it can be considered that these differences are minimal when considering the number of cows deemed for Model 1 (n = 36) compared to Model 2 (n = 312). Furthermore, Model 1 presented a slightly higher percentage of outliers, considering them as outside the 95% confidence bounds, with 3.88% compared to 3.60% for Model 2 (Figure 3), this difference is minimal and small for both models considering the number of observations in each model (Table 3).

4.3. Artificial Intelligence to Manage Heat Stress and Milk Productivity

Physical modification of the environment to reduce ambient temperature or increase heat loss from the animal body, such as shading and fans, have been previously applied for lactating buffaloes with positive results [29], and in dairy cows using mixed-flow fans [30]. However, one of the most effective methods found is spraying water over animals using sprinkler systems [31,32,33,34]. This paper proposed the implementation of Model 2 with an automated system based on an individual cow assessment combined with environmental factors obtained from an automatic meteorological station (AME) (Figure 4). The AME can be easily connected to a processing unit (microprocessor or smartphone App) that can read the RFID from cows that are going to be milked to obtain cow information required by the model (Figure 1). The model outputs can be automatically set to specific thresholds for volume and milk quality that is desired by the dairy farm. The system can then automatically control gates to direct individual cows either to a cooling system with water sprinklers, the cows to reduce heat stress or to normal milking sections. The heat-stressed cows will be assessed the next day again, if they continue to be heat stressed, they will go to the sprinkler system and get milked to avoid mastitis.
The technical advantages of the proposed system (Figure 4) are: (i) ML modeling is based on readily available environmental data by most of the dairy farms and from government services with meteorological stations close to the farms; (ii) the environmental data can be automatically extracted from government services, such as the Bureau of Meteorology (BoM, Australia) [35] or by direct connectivity of a nearby automatic meteorological station to the RFID & ML Processing Unit (Figure 4); (iii) the digital database per cow can be implemented as part of the system to incorporate data such as programmed concentrate feed, lactation days and number, milking frequency, and liveweight. This information will need to be updated by the dairy farm personnel; (iv) cows can be identified by the system with normal RFID systems to extract cow data automatically from databases, and (v) the system requires an automated gate system to draft cows to the heat stress sprinkler system or the normal milking facilities.
The managerial advantages that could be obtained by implementing the system proposed are: (i) milk volume and quality information available in real-time, per cow, and according to daily environmental conditions; (ii) prediction of actual concentrate feed intake per cow for feed monitoring management compared to programmed concentrate feed; (iii) real-time information to manage heat stress in a per cow basis to increase efficiency and maintain milk volumes and quality set as objectives, and (iv) data recorded from specific dairy farms can be incorporated in the model to increase the accuracy of target predictions.
With these considerations, an AI system for dairy farms can be implemented with reasonable investment affordable to small and medium dairy farmers. An alternative or complementary approach to an engineering solution may be to introduce dietary interventions such as betaine or antioxidants to cows likely to experience heat stress [14,28]. However, the time lag before the tissue concentrations of these nutrients are optimized could reduce the immediacy of this approach.
It should be noted that individual pasture intake could not be included in the model as the cows grazed as a single herd, so it was not measured. While this could no doubt add precision to the model, individual pasture intake cannot be measured under commercial grazing systems, and inclusion in the model would reduce its commercial utility.

5. Conclusions

The machine learning models developed in this research may be applied to assess automatically animal welfare, milk productivity, and quality. Based on the inputs of the models, this machine learning modeling technique can be applied to any dairy farm. Implementation of Artificial Intelligence in dairy farms and the ML models developed here will require minimal technological additions, automated gate, and cooling systems. This paper has shown a practical application of AI using detailed information from a robotic dairy farm for the benefit of small and medium dairy farms to increase competitiveness in an increasingly demanding international market.

Author Contributions

Conceptualization, S.F. and C.G.V.; data curation, B.C.; formal analysis, S.F. and C.G.V.; investigation, S.F., C.G.V., B.C., S.S.C., and F.R.D.; methodology, S.F. and C.G.V.; project administration, B.C. and F.R.D.; resources, B.C. and F.R.D.; software, S.F.; supervision, B.C. and F.R.D.; validation, S.F., C.G.V., B.C., S.S.C., and F.R.D.; visualization, S.F., C.G.V., and E.T.; writing—original draft, S.F. and C.G.V.; writing—review and editing, B.C., S.S.C., and F.R.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Britt, J.; Cushman, R.; Dechow, C.D.; Dobson, H.; Humblot, P.; Hutjens, M.; Jones, G.; Ruegg, P.; Sheldon, I.; Stevenson, J. Invited review: Learning from the future—A vision for dairy farms and cows in 2067. J. Dairy Sci. 2018, 101, 3722–3741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hostiou, N.; Fagon, J.; Chauvat, S.; Turlot, A.; Kling-Eveillard, F.; Boivin, X.; Allain, C. Impact of precision livestock farming on work and human-animal interactions on dairy farms. A review. Biotechnol. Agron. Soc. Environ. 2017, 21, 268–275. [Google Scholar]
  3. McCullough, C. Global dairy outlook by 2030: Industry. The Dairy Mail 2019, 26, 12–15. [Google Scholar]
  4. Bear, C.; Holloway, L. Beyond resistance: Geographies of divergent more-than-human conduct in robotic milking. Geoforum 2019, 104, 212–221. [Google Scholar] [CrossRef]
  5. Heyden, T. The Cows That Queue up to Milk Themselves. Available online: https://robohub.org/the-cows-that-queue-up-to-milk-themselves-bbc-news/ (accessed on 3 May 2020).
  6. Jorquera-Chavez, M.; Fuentes, S.; Dunshea, F.R.; Warner, R.D.; Poblete, T.; Jongman, E.C. Modelling and Validation of Computer Vision Techniques to Assess Heart Rate, Eye Temperature, Ear-Base Temperature and Respiration Rate in Cattle. Animals 2019, 9, 1089. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ji, B.; Banhazi, T.; Ghahramani, A.; Bowtell, L.; Wang, C.; Li, B. Modelling of heat stress in a robotic dairy farm. Part 3: Rumination and milking performance. Biosyst. Eng. 2020. [Google Scholar] [CrossRef]
  8. Ji, B.; Banhazi, T.; Ghahramani, A.; Bowtell, L.; Wang, C.; Li, B. Modelling of heat stress in a robotic dairy farm. Part 2: Identifying the specific thresholds with production factors. Biosyst. Eng. 2019. [Google Scholar] [CrossRef]
  9. Ji, B.; Banhazi, T.; Ghahramani, A.; Bowtell, L.; Wang, C.; Li, B. Modelling of heat stress in a robotic dairy farm. Part 1: Thermal comfort indices as the indicators of production loss. Biosyst. Eng. 2019. [Google Scholar] [CrossRef]
  10. Ekine-Dzivenu, C.; Mrode, R.A.; Ojango, J.M.; Okeyo Mwai, A. Evaluating the impact of heat stress as measured by temperature-humidity index (THI) on test-day milk yield of dairy cattle in Tanzania. In Proceedings of the Seventh All Africa conference on Animal Agriculture, Accra, Ghana, 29 July–2 August 2019. [Google Scholar]
  11. Nascimento, F.G.d.O.; Aguiar, H.C.P.; Rodrigues, G.M.; Guimarães, E.C.; Nascimento, M.R.B.d.M. What is the best temperature-humidity index equation to indicate heat stress in crossbred dairy calves in a tropical environment? Ciência Rural 2019, 49, e20180132. [Google Scholar] [CrossRef]
  12. Berman, A.; Horovitz, T.; Kaim, M.; Gacitua, H. A comparison of THI indices leads to a sensible heat-based heat stress index for shaded cattle that aligns temperature and humidity stress. Int. J. Biometeorol. 2016, 60, 1453–1462. [Google Scholar] [CrossRef]
  13. Carolan, M. Automated agrifood futures: Robotics, labor and the distributive politics of digital agriculture. J. Peasant Stud. 2020, 47, 184–207. [Google Scholar] [CrossRef]
  14. Dunshea, F.R.; Oluboyede, K.; DiGiacomo, K.; Leury, B.J.; Cottrell, J.J. Betaine improves milk yield in grazing dairy cows supplemented with concentrates at high temperatures. Animals 2019, 9, 57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Bolton, D. The computation of equivalent potential temperature. Mon. Weather Rev. 1980, 108, 1046–1053. [Google Scholar] [CrossRef] [Green Version]
  16. Goodwin, R. Wet Bulb, Matlab Central File Exchange. Available online: https://www.mathworks.com/matlabcentral/fileexchange/50785-wet-bulb (accessed on 20 April 2020).
  17. Gonzalez Viejo, C.; Torrico, D.; Dunshea, F.; Fuentes, S. Development of Artificial Neural Network Models to Assess Beer Acceptability Based on Sensory Properties Using a Robotic Pourer: A Comparative Model Approach to Achieve an Artificial Intelligence System. Beverages 2019, 5, 33. [Google Scholar] [CrossRef] [Green Version]
  18. Nguyen, T.T.; Bowman, P.J.; Haile-Mariam, M.; Nieuwhof, G.J.; Hayes, B.J.; Pryce, J.E. Implementation of a breeding value for heat tolerance in Australian dairy cattle. J. Dairy Sci. 2017, 100, 7362–7367. [Google Scholar] [CrossRef] [Green Version]
  19. Ouellet, V.; Cabrera, V.; Fadul-Pacheco, L.; Charbonneau, É. The relationship between the number of consecutive days with heat stress and milk production of Holstein dairy cows raised in a humid continental climate. J. Dairy Sci. 2019, 102, 8537–8545. [Google Scholar] [CrossRef]
  20. Gunn, K.M.; Holly, M.A.; Veith, T.L.; Buda, A.R.; Prasad, R.; Rotz, C.A.; Soder, K.J.; Stoner, A.M. Projected heat stress challenges and abatement opportunities for US milk production. PLoS ONE 2019, 14. [Google Scholar] [CrossRef] [Green Version]
  21. Summer, A.; Lora, I.; Formaggioni, P.; Gottardo, F. Impact of heat stress on milk and meat production. Anim. Front. 2019, 9, 39–46. [Google Scholar] [CrossRef]
  22. Sigdel, A.; Abdollahi-Arpanahi, R.; Aguilar, I.; Peñagaricano, F. Whole Genome Mapping Reveals Novel Genes and Pathways Involved in Milk Production Under Heat Stress in US Holstein Cows. Front. Genet. 2019, 10, 928. [Google Scholar] [CrossRef]
  23. Alba, D.F.; Campigotto, G.; Cazarotto, C.J.; dos Santos, D.S.; Gebert, R.R.; Reis, J.H.; Souza, C.F.; Baldissera, M.D.; Gindri, A.L.; Kempka, A.P. Use of grape residue flour in lactating dairy sheep in heat stress: Effects on health, milk production and quality. J. Therm. Biol. 2019, 82, 197–205. [Google Scholar] [CrossRef]
  24. dos Santos, D.d.S.; Klauck, V.; Campigotto, G.; Alba, D.F.; Dos Reis, J.H.; Gebert, R.R.; Souza, C.F.; Baldissera, M.D.; Schogor, A.L.B.; Santos, I.D. Benefits of the inclusion of açai oil in the diet of dairy sheep in heat stress on health and milk production and quality. J. Therm. Biol. 2019, 84, 250–258. [Google Scholar] [CrossRef] [PubMed]
  25. Hall, L.; Dunshea, F.; Allen, J.; Rungruang, S.; Collier, J.; Long, N.; Collier, R.J. Evaluation of dietary betaine in lactating Holstein cows subjected to heat stress. J. Dairy Sci. 2016, 99, 9745–9753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Gonzalez-Rivas, P.A.; Sullivan, M.; Cottrell, J.J.; Leury, B.J.; Gaughan, J.B.; Dunshea, F.R. Effect of feeding slowly fermentable grains on productive variables and amelioration of heat stress in lactating dairy cows in a sub-tropical summer. Trop. Anim. Health Prod. 2018, 50, 1763–1769. [Google Scholar] [CrossRef] [PubMed]
  27. Coppock, C.E. Reduce Heat Stress in High-Producing Dairy Cows with Feeding and Nutrition Management. In Dairy Science Handbook; CRC Press: Boca Raton, FL, USA, 2019; pp. 315–324. [Google Scholar]
  28. Dunshea, F.R.; Leury, B.J.; Fahri, F.; DiGiacomo, K.; Hung, A.; Chauhan, S.; Clarke, I.J.; Collier, R.; Little, S.; Baumgard, L. Amelioration of thermal stress impacts in dairy cows. Anim. Prod. Sci. 2013, 53, 965–975. [Google Scholar] [CrossRef]
  29. Ahmad, M.; Bhatti, J.A.; Abdullah, M.; Ullah, R.; ul Ain, Q.; Hasni, M.S.; Ali, M.; Rashid, A.; Qaisar, I.; Rashid, G. Different ambient management intervention techniques and their effect on milk production and physiological parameters of lactating NiliRavi buffaloes during hot dry summer of subtropical region. Trop. Anim. Health Prod. 2019, 51, 911–918. [Google Scholar] [CrossRef]
  30. Yao, C.; Shi, Z.; Zhao, Y.; Ding, T. Effect of Mixed-Flow Fans with a Newly Shaped Diffuser on Heat Stress of Dairy Cows Based on CFD. Energies 2019, 12, 4315. [Google Scholar] [CrossRef] [Green Version]
  31. Tresoldi, G.; Schütz, K.E.; Tucker, C.B. Cooling cows with sprinklers: Effects of soaker flow rate and timing on behavioral and physiological responses to heat load and production. J. Dairy Sci. 2019, 102, 528–538. [Google Scholar] [CrossRef] [Green Version]
  32. Sruthi, S.; Sasidharan, M.; Anil, K.; Harikumar, S.; Simon, S. Effect of automated intermittent wetting and forced ventilation on the physiological parameters and milk production of Murrah buffaloes in humid tropics. Pharma Innov. J. 2019, 8, 315–319. [Google Scholar]
  33. Imbabi, T.; Hassan, T.; Radwan, A.; Soliman, A. Production, haematological and biochemical metabolites of Egyptian buffaloes (Bubalus bubalis) during the hot summer months in Egypt. Slovak J. Anim. Sci. 2019, 52, 152–159. [Google Scholar]
  34. Pinto, S.; Hoffmann, G.; Ammon, C.; Heuwieser, W.; Levit, H.; Halachmi, I.; Amon, T. Effect of two cooling frequencies on respiration rate in lactating dairy cows under hot and humid climate conditions. Annal. Anim. Sci. 2019, 19, 821–834. [Google Scholar] [CrossRef] [Green Version]
  35. Meteorology, B.O. BoM Automated Data Services. Available online: http://www.bom.gov.au/catalogue/data-feeds.shtml (accessed on 10 February 2020).
Figure 1. Diagram of the two-layer feedforward regression models with a tan-sigmoid function in the hidden layer and linear transfer function in the output layer. Abbreviations: THI: Temperature-humidity index; W: Weights; b: Bias.
Figure 1. Diagram of the two-layer feedforward regression models with a tan-sigmoid function in the hidden layer and linear transfer function in the output layer. Abbreviations: THI: Temperature-humidity index; W: Weights; b: Bias.
Sensors 20 02975 g001
Figure 2. Mean values per season of each year for temperature-humidity index (THI9) and the four parameters used as targets in the machine learning (ML) models to represent the effect of different weather patterns on potential heat stress, milk productivity, and quality.
Figure 2. Mean values per season of each year for temperature-humidity index (THI9) and the four parameters used as targets in the machine learning (ML) models to represent the effect of different weather patterns on potential heat stress, milk productivity, and quality.
Sensors 20 02975 g002
Figure 3. Overall regression graphs of (a) Model 1: Using the 36 cows with similar heat tolerance (93–112), and (b) Model 2: Using data from 312 cows.
Figure 3. Overall regression graphs of (a) Model 1: Using the 36 cows with similar heat tolerance (93–112), and (b) Model 2: Using data from 312 cows.
Sensors 20 02975 g003
Figure 4. Proposed artificial intelligence (AI) application based on the automated processing of meteorological station and radio frequency identification system (RFID) for specific cow data input and machine learning (ML) processing. This system activates the gate system to draft cows to a cooling system or normal milking.
Figure 4. Proposed artificial intelligence (AI) application based on the automated processing of meteorological station and radio frequency identification system (RFID) for specific cow data input and machine learning (ML) processing. This system activates the gate system to draft cows to a cooling system or normal milking.
Sensors 20 02975 g004
Table 1. Minimum, maximum, and mean values of the parameters used as inputs to develop the machine learning models.
Table 1. Minimum, maximum, and mean values of the parameters used as inputs to develop the machine learning models.
Parameter/Year2016 *201720182019 *
MinMaxMeanSDMinMaxMeanSDMinMaxMeanSDMinMaxMeanSD
T (°C)7.937.819.37.158.342.022.28.278.943.322.67.9016.344.931.76.07
RH (%)66.010095.66.2056.210092.39.2744.210087.911.6739.292.669.212.13
Tdp (°C)3.824.911.63.402.122.811.53.971.622.110.23.754.321.213.53.79
Twet (°C)6.625.313.83.566.825.314.74.205.924.514.13.9711.324.118.63.11
Rainfall (mm day−1)0.034.03.96.760.031.61.864.550.037.81.44.070.05.00.30.99
Wind speed (km h−1)5.238.315.15.225.834.315.35.195.138.016.25.829.439.819.56.14
Wind direction (°)127.7360.0344.026.54247.2360.0345.123.76112.2360.0341.832.81241.7360.0338.930.62
THI157.289.670.67.3758.494.673.38.7258.792.973.38.3467.396.382.86.25
THI244.578.658.17.0845.481.260.28.4244.579.259.68.0154.980.568.66.06
THI344.781.160.08.2446.086.862.99.7946.383.862.69.3456.387.873.26.99
THI450.883.264.27.3752.088.266.98.7252.386.566.98.3460.989.976.46.25
THI547.182.263.48.1547.586.566.59.1349.084.366.98.4860.487.576.25.55
THI659.191.572.07.3360.297.074.78.7560.594.474.68.3968.798.784.36.56
THI750.883.663.97.4051.989.166.68.8452.286.566.58.4760.690.976.36.62
THI847.182.063.48.0947.586.366.59.0749.084.066.88.4160.487.276.05.49
THI933.486.658.812.4533.492.663.713.7336.589.464.412.6155.493.778.17.83
Programmed concentrate feed (kg day−1)0.015.08.93.030.023.08.53.110.015.77.83.150.08.05.12.30
Lactation number1.06.02.70.971.07.03.01.241.07.02.31.611.08.03.01.75
Lactation days0.0736.0225158.170.0668.0198.1139.280.0705.0228.3142.260.0755.0227.5144.01
Milking frequency (per day)0.05.02.40.710.06.02.50.750.06.02.40.840.05.01.90.81
Liveweight (kg)373.0938.0677.782.85428.0951.0668.278.25335.0959.0655.484.57410.0896.0629.571.86
* Values from 2016 cover from June to December and 2019 cover from January to March. Abbreviations: Min: Minimum; Max: Maximum; T: Temperature; RH: Relative humidity; Tdp: Dewpoint temperature; Twet: Wet-bulb temperature; THI: Temperature-humidity index; SD: Standard deviation.
Table 2. Minimum, maximum, and mean values of the parameters used as targets to develop the machine learning models.
Table 2. Minimum, maximum, and mean values of the parameters used as targets to develop the machine learning models.
Parameter/Year2016 *201720182019 *
MinMaxMeanMinMaxMeanMinMaxMeanMinMaxMean
Milk yield (kg day−1)0.065.428.10.060.230.70.061.228.80.052.121.2
Milk protein (%)1.85.83.31.86.13.22.25.83.40.94.93.1
Milk fat (%)1.010.74.20.810.24.00.710.34.20.710.94.3
Concentrate feed intake (kg day−1)0.019.57.30.024.37.40.018.86.70.010.64.0
* Values from 2016 cover from June to December and 2019 cover from January to March. Abbreviations: Min: Minimum; Max: Maximum.
Table 3. Statistical results of each stage of the machine learning models.
Table 3. Statistical results of each stage of the machine learning models.
StageSamples (Cows x Days)Observations (Samples x Targets)RbPerformance (MSE)
Model 1
Training20,38081,5200.870.760.0186
Testing873434,9360.860.760.0189
Overall29,114116,4560.870.76-
Model 2
Training116,521466,0840.860.740.0154
Testing49,938199,7520.860.740.0157
Overall166,459665,8360.860.74-
Abbreviations: R: Correlation coefficient; b: Slope; MSE: Means squared error.

Share and Cite

MDPI and ACS Style

Fuentes, S.; Gonzalez Viejo, C.; Cullen, B.; Tongson, E.; Chauhan, S.S.; Dunshea, F.R. Artificial Intelligence Applied to a Robotic Dairy Farm to Model Milk Productivity and Quality based on Cow Data and Daily Environmental Parameters. Sensors 2020, 20, 2975. https://doi.org/10.3390/s20102975

AMA Style

Fuentes S, Gonzalez Viejo C, Cullen B, Tongson E, Chauhan SS, Dunshea FR. Artificial Intelligence Applied to a Robotic Dairy Farm to Model Milk Productivity and Quality based on Cow Data and Daily Environmental Parameters. Sensors. 2020; 20(10):2975. https://doi.org/10.3390/s20102975

Chicago/Turabian Style

Fuentes, Sigfredo, Claudia Gonzalez Viejo, Brendan Cullen, Eden Tongson, Surinder S. Chauhan, and Frank R. Dunshea. 2020. "Artificial Intelligence Applied to a Robotic Dairy Farm to Model Milk Productivity and Quality based on Cow Data and Daily Environmental Parameters" Sensors 20, no. 10: 2975. https://doi.org/10.3390/s20102975

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop