|
|
||||||||
Department of Dairy Science University of Wisconsin, Madison WI 53706
Corresponding author:
R. W. Palmer; e-mail:
rwpalmer{at}facstaff.wisc.edu.
| ABSTRACT |
|---|
|
|
|---|
Key Words: cow behavior automatic milking system feeding behavior milk yield
Abbreviation key: parlor = double-6 herringbone parlor, robot = automated milking system
| INTRODUCTION |
|---|
|
|
|---|
Initial reports with automatic milking systems indicate that cow feeding and milking activity is more evenly distributed across a 24-h period. However, visits to a robotic milker are lower during the night and early morning hours (Wendl et al., 2000; Olofsson et al., 2000). Wendl et al. (2000) found that the highest number of milking visits to the robotic milker were in the morning. Olofsson et al. (2000) reported that feed intakes were lower during the night for cows in an automatic milking system. OConnell et al. (1993) reviewed a study which found increased cow feeding activity occurred during the day and early evening hours (9 a.m. to 3 p.m. and 5 p.m. to 1 a.m.) vs. late night and early morning (1 a.m. to 7 a.m.), when cows were studied in confinement. The type of milking system used was not described but likely was not an automatic system.
Association with an opportunity to eat may motivate cows to visit robotic milking machines more frequently. A one-way gate system forces a cow to pass through a robotic milking unit to get to the feed bunk. Approximately 30 to 56% of total visits to a robotic milker were nonmilking visits, with the highest frequency occurring in the afternoon (Wendl et al., 2000; Morita et al., 2000). Therefore, cows do apparently visit the robot in order to reach the feed bunk.
An issue that has not been addressed in previous studies with automated milking systems is cow behavior and milk response during hot weather. Maintaining milk yield in dairy cattle is essential for a dairy farms profitability during periods of thermal stress. Monitoring changes in milk yield is important to ensure production levels are maintained. Cow characteristics, such as DIM and parity; management practices such as use of bST and frequency of milking; and weather conditions, such as temperature and relative humidity, influence milk yield.
Previous studies have examined effects of thermal stress on milk yield. Maximum, minimum, mean temperatures, relative humidity, and calculated temperature humidity indices have been evaluated as variables representing heat stress. Ravagnolo and Misztal (2000) and Maust et al. (1972) found maximum temperature to be a good measure of the weathers effect on milk yield. High temperatures and temperature humidity index values (>72) negatively affected production performance (Linvill and Pardue, 1992; Armstrong, 1994; West, 1994; Ravagnolo and Misztal, 2000). Previous days weather conditions have been considered as a possible factor affecting test day yields by Ravagnolo and Misztal (2000). Maust et al. (1972) found weather conditions 2, 3, and 4 d previous to be significantly correlated with a current milk yield.
Milking frequency can also directly affect milk yield and is a significant variable in predicting yields (Ravagnolo and Misztal, 2000). Research has focused on the yield response of two times per day milking vs. three times per day, as well as differences between parity groups. A review of published studies found a relatively fixed milk yield response of 7.7 lb/d (3.5 kg), by increasing milking frequency from two to three milkings per day for both primiparous and multiparous cows (Erdman and Varner, 1995). Effects of milking frequency and parity on yield response of selected (AI sires) vs. nonselected cows were reported by Barnes et al. (1990). Their results showed an increase in second parity yield with three-times-a-day milking over first parity yield in both genetic groups.
The objectives of the current study were to: 1) compare feeding behavior of cows milked in a conventional milking parlor (parlor) with those in an automatic milking (robot) system; 2) monitor cow milking behavior for a robotic milking system; 3) determine the effect on milk yield of cows milked with a conventional herringbone milking parlor 2x and cows milked in a robotic milking system; and 4) determine effects of ambient temperature during hot weather on daily milk yields within the two milking systems.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Video Monitoring Equipment
Four HOBO (Onset Computer Corporation, Pocasset, MA) H8 Pro Series sensors recorded temperature every 15 min. One sensor was located in each quadrant of the barn, northwest (NW), northeast (NE), southwest (SW), southeast (SE). Sensors were 38 to 40 ft (12.6 to 13.2 m) from each end wall, and 9 to 10 ft (2.9 to 3.3 m) above the alley floor surface over the front of the interior row of stalls. Data were downloaded from the sensors via a HOBO Shuttle. Videotape images were observed once per hour, as close to the beginning of the hour as the camera sequence would allow, on the 15th of each month from May 2001 to January 2002.
Animal Grouping
The south side of the barn housed cows milked in a conventional double-6 herringbone parlor, and the northside-housed cows milked with a Bou-matic (Bou-matic Dairy Equipment Company, Madison, WI) robot. Barn characteristics and videotape observation methods are described in more detail by Wagner-Storch et al. (2002). Briefly, the robot pen was equipped with a series of one-way gates, which only allowed cows to access the feed alley through the robot. The robot maintained a 10-d rolling average hourly milk level for each cow and used it to estimate the amount of milk expected from the cow based on the length of time since she was last milked. If expected milk was greater than the milk threshold level of 13.2 lb (6 kg), the cow was milked, if not she was released to the feed bunk area. The robot and adjacent holding area were visible 24 h per day because the area near the robot was lighted. The status of the robot, in-use (cow being milked) or not in use, and the number of cows in the robot holding area were recorded. A cow was assumed to be in the robot holding area if she was past the northeast waterer or if held in the gated area near the robot.
During the 9-mo study period, the number of cows milked in the parlor ranged from 43 to 58 cows, with a mean of 52 cows, and the number of cows milked by the robot ranged from 30 to 43 cows, with a mean of 36 cows.
Feed Bunk Monitoring
The number of cows at the feed bunk was recorded for both sides of the barn. A cow was counted as "at the feed bunk" if her head was through a self-locking manger stall (self-locks) opening. Wooden support posts blocked the view of approximately 2 ft (0.6 m) of the feed bunk on both sides of the barn. Both groups of cows were fed once daily (at approximately 11:00 a.m.), and cows were fed enough to support an 8% feed refusal.
Scheduled milking (6:00 a.m. and 5:00 p.m) and feeding (11:00 a.m.) times were set for cows milked in the parlor, but actual times varied. For cows milked by the robot, there were a scheduled wash cycle (2:00 p.m.) of the robot and a daily scheduled human-intervention time (7:00 a.m.) that influenced cow behavior. Human-intervention time refers to the scheduled time of day when an employee would identify cows not milked in the previous 15 h or cows that did not routinely visit the robot without human-intervention. These cows were manually guided to the holding area in front of the robot.
Long-day lighting in the barn prevented accurate collection of feed-bunk data during the scheduled dark period, 9:30 p.m. to 3:30 a.m., and feeding behavior data from those times are not included. Feed-bunk data were summarized as an average percentage of total cows at the feed bunk, by management group, for each hour, with calculated standard errors of those percentages. (Example calculation for average percentage of total cows at the south side (parlor) feed bunk at 2 p.m. for 9 d over 9 mo.: [(5/43) + (8/48) + (5/54) + (6/54) + (7/53) + (6/54) + (29/52) + (9/53) + (12/52)] / 9 = 0.188 x 100 = 18.8 ± 4.8%). The status of the robot and number of cows in the robot holding area were summarized as average hourly percentages ± standard errors. A paired t-test (SAS, 1999) was run to compare feed bunk percentages at each hour from 4 a.m. to 9 p.m. for each management group. A transformation, arcsine (square root [average proportion]) was done to account for unequal variances. Data from 12 a.m. to 3 a.m. and from 10 p.m. to 11 p.m. were not tested due to insufficient observations.
Milk Yield Model
Temperature, milk yield, and cow characteristic data were collected from July 9 to August 16, 2001. This period, during the hot summer conditions, was specifically selected to study the effect of parity, DIM, milking frequency (parlor = 2x, robot = 1x to 4x), and high temperatures on daily milk yields.
Temperature data used in the analysis were maximum value per day for each side of the barn. The maximum daily temperature for cows on the parlor side of the barn was determined from the higher of SW and SE sensors, and the maximum daily temperature for cows on the robot side was determined from the higher of NW and NE sensors.
Cow milk yield, parity, DIM, and milking frequency were collected from electronic records. Records were removed from the analysis if determined not to be plausible: milking frequency unequal to 2 on the parlor side, greater than 4 on the robot side or milk yield less than 10 lb (4.54 kg).
Data from 122 different cows were used in the analysis. Six cows had records for both sides of the barn, because they were moved from one side to the other during the study. Therefore, data from 84 cows in the parlor group and 44 cows in the robot group were available. Total number of cow-day observations in the analysis was 2982.
Daily milk yield data were analyzed with the MIXED Procedure in Statistical Analysis Software (SAS, 1999). The following model was tested: yhijk = ch, + mfi + pj + b1dimhijk + b2thijk + b3lag1thijk + b4lag2thijk + b5lag3thijk + ehijk where yhijk is daily milk yield for cow h with milking frequency i, parity j, on day of lactation k. Continuous variables were days in milk (dim), maximum temperature on observation day (t), maximum temperature from one day previous (lag1t), maximum temperature from 2 d previous (lag2t), maximum temperature from three days previous (lag3t), and residual. A repeated statement, type = AR(1), was used to account for autocorrelation errors due to repeated measures on each cow (Littell et al., 1998). Milking frequency, parity, and cow were categorical variables. The main objective of this research was to measure robot and parlor differences. However, it should be noted that due to physical constraints of the facility, the side of the barn was confounded with milking frequency and milking system, and therefore not included in the model. Replication of barns with robot and parlor milking systems on both sides was not possible. The average temperature on the south side (parlor) was almost 0.5°F (0.28°C) warmer in the summer and over 1.0°F (0.56°C) warmer in the winter than on the north side (robot). This was not accounted for in the milk yield model and may have biased the results shown.
During July, 50% (parlor) and 57% (robot) of the cows received bST treatments. In August, 48% (parlor) and 52% (robot) of the cows received bST treatments. However, bST use was similar in both groups and variation from bST use was included in the calculated residual value, making significance tests more conservative.
Models testing maximum and mean temperature effects separately, and maximum and mean temperatures together, were analyzed for the observation day, and 1, 2, and 3 d previous to that day.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
The paired-t test on the robot side showed feeding percentages early in the morning (4 to 6 a.m.) were lower (P < 0.05) than after the 7 a.m. scheduled human-intervention time. Few differences were found during the day and early evening (7 a.m. to 9 p.m.), except the highest numerical percentage (10 a.m.) was higher than late afternoon (5 to 6 p.m.), before the percentage of cows at the feed bunk increased again. The percentage of cows at the feed bunk at 9 a.m. was less than that at 10 a.m., which shows the cow numbers at the bunk continued to rise to a peak about 3 h after humanintervention and again from 7 to 9 p.m., after which feeding activity decreased. The 7 a.m. percentage was different from 10 a.m. and 12 p.m., when feed bunk activity was high. It appears that cows started going through the robot after the human-intervention time, causing the cow numbers to increase at the feed bunk until about 10 a.m., remained relatively constant until 5 p.m., and then rose again as cows came to the robot to be milked a second time. No difference was found in feeding behavior associated with the scheduled wash cycle of the robot or the scheduled feeding time for cows in the robot group.
The only significant differences (P< 0.05) between the parlor and robot groups in percentages of cows in the feeding area were at 10 a.m. and 9 p.m. At 10 a.m., the average percentage of total cows at the feed bunk in the robot group increased to approximately 20% after the cow build-up, which was initiated by the earlier human-intervention time. The average percentage of total cows at the parlor feed bunk decreased because cows had eaten after being milked and were lying down. At 9 p.m., cows in the robot group were returning to the feed bunk after being milked, whereas cows on the parlor side had been milked at 5 p.m., had eaten, and again were lying down.
Variability in feed bunk activity for the parlor group increased the hour after the scheduled milking, as indicated by the standard error bars. There was less variability on the robot side than on the parlor side because the parlor cows had a less structured movement and activity schedule.
Activity Within the Automatic Milking System
The percentage ± S.E. of total cows in the robot or robot holding area by hour are summarized in Figure 2
. The robot was used extensively from 7 a.m. to 7 p.m. with the average number of cows in the robot ranging from approximately 0.5 to 1.0 cows or about 2% of the group. A smaller average percentage of total cows were observed in the robot during the evening, night and early morning, ranging from approximately 0 to < 2% from 8 p.m. to 6 a.m. A paired t-test was run to compare average percentage of cows in the robot at each hour of the day for the robot management group. No cows were observed in the robot at the 5 a.m. time, so the robots 5 a.m. average usage was different from all other times. Activity at 12 a.m. was less than the two high activity times of 8 a.m. to 1 p.m. and 4 to 7 p.m. (P < 0.05) and less than 3 p.m. (P < 0.10). Activity increased at 1 a.m. with a greater average percentage than 3 and 5 a.m. and lower than 10 a.m. The 2 a.m. observation showed greater activity than at 5 a.m. and lower activity than at 8 a.m., 10 a.m., and 7 p.m. The 3 a.m. average was lower than 7 a.m. to 10 p.m. The 4 a.m. average was lower than 8 a.m. to 12 p.m. and 7 p.m. The decrease in activity at the scheduled 2 p.m. wash time was only different from the robot usage at 7 p.m. Based on these results it appears the highest robot activity was from 8 a.m. to 1 p.m. and 3 to 7 p.m. The first period was initiated by and followed human intervention and the second about 8 h later when many of the cows entered the robot a second time. Activity increased again, but to a lesser degree, at 1 to 2 a.m. and was assumed to be associated with cows being milked a third time.
|
Average percentage of cows in the robot holding area was greatest from 8 to 11 a.m. at about 12% and 2 to 7 p.m. with approximately 9 to 11% of the cows waiting to use the robot. The average percentage of cows in the robot holding area was smaller in the evening, night, and early morning, ranging from approximately 1 to 3% of the cows. A paired t-test was run to compare average percentage of cows in the robot holding area at each hour of the day for the robot management group. The 12 a.m. average was lower than the 7 a.m. to 7 p.m. averages (P < 0.05). The 1 a.m. average was lower than the 7 a.m. to 3 p.m and 6 to 7 p.m. averages, but not different from 4 and 5 p.m. because of their large variation. The 2 and 3 a.m. averages were lower than the 7 a.m. to 7 p.m. averages. The 4 a.m. average was lower than the 8 a.m., 11 a.m., and 2 to 6 p.m. averages. The 5 a.m. average was lower than the 6 a.m. to 7 p.m. and the 10 p.m. averages. The 6 a.m. observation was lower than the 8 a.m. and 6 p.m. averages. The 7 a.m. average was less than 8 a.m. and greater than the 11 p.m. average. Observations from 8 a.m. to 7 p.m. were greater than the 11 p.m. average. These results show that activity in the robot holding area increased after the 7 a.m. human-intervention. The increase observed at 7 a.m. was assumed to be caused by variability of actual vs. scheduled human-intervention times.
The average number of cows in the robot holding area was greatest at 8 a.m. to 11 a.m. (four to five cows), and again at 2 p.m. to 7 p.m. (three to four cows). The average number of cows waiting in the robot holding area was less than two from 8 p.m. to 6 a.m. Activity in the robot holding area increased after human-intervention time and remained high until these cows went through the robot and increased again later in the day as cows were waiting to be milked a second time. Had there been closer to the recommended 60 cows for the robot system, presence of cows in the holding area may have had a different or more pronounced distribution pattern.
Variability in the number of cows in the robot holding area increased during daytime hours and was highest from 8 to 11 a.m. and 3 to 5 p.m. This increased variation may have been caused by differences in human-intervention times or cows may have been waiting to use the robot, not for milking, but to gain access to the feed bunk, as suggested by earlier research (Wendl et al., 2000; Morita et al., 2000).
Milk Yield Comparisons Between Milking Systems
Daily average milk yields from July 9 to August 16, 2001, were analyzed to examine the effects of milking system (and frequency of milking), parity, DIM, and ambient temperature on milk yield. The parlor side had 84 different cows and the robot side had 44 different cows during this 39-d period. Mean Start-DIM, defined as the DIM for a cow at the beginning of the cows respective observation period, were similar between robot and parlor groups and averaged 221.9 ± 7.9 DIM, indicating that most cows were in mid to late lactation during the summer observations. However, cows in the robot system had lower parity (1.6 vs. 2.4 ± 0.1), higher milking frequency (2.4 ± 0.02 vs. 2.0 ± 0.0), and slightly higher milk production (58.1 vs. 56.7 ± 0.5 lb, or 26.4 vs. 25.8 ± 0.2 kg) than parlor cows (P < 0.05; Table 1
). Lower average parity of cows in the robot group was partially the result of the animal selection process, where cows with uneven udders were assigned to the parlor side to avoid robot unit attachment problems. The lower average parity of cows in the robot management group implies a higher percentage of cows in the robot management group were parity one animals. This would partially explain why milk yield for robot cows was only 1.4 lb (0.6 kg) higher than for the parlor group even though milking frequency was higher.
|
|
|
Because the average DIM of the parlor cows and robot cows were not significantly different, it appears that the increased milking frequency by the robot had the greatest influence on the milk yield differences seen. This supports the idea that an automatic milking system can be used to increase milking frequency and daily milk yields. These data also showed that although average milking frequency was increased by 0.4x (2.0x to 2.4x), the increase in production was not 40% of the expected difference in production from 2x to 3x. This can be partially attributed to the higher percentage of parity one cows, which can be expected to have lower daily milk yields.
The model showed the maximum temperature one day before the milk yield observation had more impact on milk yield than the temperature on the day of or 2 d before observation day. This result was in agreement with previous research (Maust et al., 1972; Ravagnolo and Misztal, 2000). Cows exposed to high temperatures during summer weather conditions experienced a decrease in milk production (Linvill and Pardue, 1992; Armstrong, 1994). The main reason for decreased production may be associated with a decrease in DMI caused by high temperatures. Feed intakes were not measured and could not be included in the model.
| CONCLUSIONS |
|---|
|
|
|---|
Although the robot management group contained fewer than the recommended 60 cows per robot, it achieved an average of only 2.4 milking per cow per day. This could have been caused by one of the following: 1) cows were conditioned to 2x milking before entering this pen; 2) cows were in the later part of their lactation and did not have sufficient udder pressure or desire to eat as early lactation cows would have; or 3) the low stocking rate could have reduced their natural herding behavior. Milking frequency best explained variation in milk yield and automatic milking systems may effectively increase daily milk yield by increasing milking frequency. Although milking frequency was increased by 0.4x (2.0x to 2.4x) for robot-milked cows, the increase in average milk yield was not 40% of the expected difference in production from 2x to 3x, likely because of a higher percentage of first parity cows in the robot group. Finally, maximum daily temperature from the observation day and from 1 and 2 d previous were the best temperature measures to predict the effects of hot weather on milk yield.
Future research should consider the effects of other management practices, DMI, and milking interval lengths on milk yield and implications of animal movement on animal handling systems for automatic milking systems.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication May 14, 2002. Accepted for publication September 26, 2002.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. M. Svennersten-Sjaunja and G. Pettersson Pros and cons of automatic milking in Europe J Anim Sci, March 1, 2008; 86(13_suppl): 37 - 46. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Bach, C. Iglesias, S. Calsamiglia, and M. Devant Effect of Amount of Concentrate Offered in Automatic Milking Systems on Milking Frequency, Feeding Behavior, and Milk Production of Dairy Cattle Consuming High Amounts of Corn Silage J Dairy Sci, November 1, 2007; 90(11): 5049 - 5055. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Speroni, G. Pirlo, and S. Lolli Effect of Automatic Milking Systems on Milk Yield in a Hot Environment J Dairy Sci, December 1, 2006; 89(12): 4687 - 4693. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. L. Mentink and N. B. Cook Short Communication: Feed Bunk Utilization in Dairy Cows Housed in Pens with Either Two or Three Rows of Free Stalls J Dairy Sci, January 1, 2006; 89(1): 134 - 138. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Abeni, L. Calamari, F. Calza, M. Speroni, G. Bertoni, and G. Pirlo Welfare Assessment Based on Metabolic and Endocrine Aspects in Primiparous Cows Milked in a Parlor or with an Automatic Milking System J Dairy Sci, October 1, 2005; 88(10): 3542 - 3552. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. DeVries and M. A. G. von Keyserlingk Time of Feed Delivery Affects the Feeding and Lying Patterns of Dairy Cows J Dairy Sci, February 1, 2005; 88(2): 625 - 631. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Muller and L. Schrader Individual Consistency of Dairy Cows' Activity in Their Home Pen J Dairy Sci, January 1, 2005; 88(1): 171 - 175. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Miron, M. Nikbachat, A. Zenou, D. Ben-Ghedalia, R. Solomon, E. Shoshani, I. Halachmi, N. Livshin, A. Antler, and E. Maltz Lactation Performance and Feeding Behavior of Dairy Cows Supplemented Via Automatic Feeders with Soy Hulls or Barley Based Pellets J Dairy Sci, November 1, 2004; 87(11): 3808 - 3815. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Miron, E. Yosef, M. Nikbachat, A. Zenou, E. Maltz, I. Halachmi, and D. Ben-Ghedalia Feeding Behavior and Performance of Dairy Cows Fed Pelleted Nonroughage Fiber Byproducts J Dairy Sci, May 1, 2004; 87(5): 1372 - 1379. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. DeVries, M. A. G. von Keyserlingk, and K. A. Beauchemin Short Communication: Diurnal Feeding Pattern of Lactating Dairy Cows J Dairy Sci, December 1, 2003; 86(12): 4079 - 4082. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |