Advertisement
Research Article| Volume 99, ISSUE 5, P3732-3743, May 2016

Cow- and farm-level risk factors for lameness on dairy farms with automated milking systems

Open ArchivePublished:February 25, 2016DOI:https://doi.org/10.3168/jds.2015-10414

      Abstract

      Lameness is a major concern to animal health and welfare within the dairy industry. Our objectives were to describe the prevalence of lameness in high-producing cows on farms with automated milking systems (AMS) and to identify the main risk factors for lameness at the animal and farm level. We visited 36 AMS farms across Canada and Michigan. Farm-level factors related to stall design, bedding use, flooring, and stocking rates were recorded by trained observers. Cows were scored for lameness, leg injuries, body condition (BCS), and body size (hip width and rump height; n = 1,378; 25–40 cows/farm). Mean herd prevalence of clinical lameness was 15% (range = 2.5–46%). Stall width relative to cow size and parity was found to be the most important factor associated with lameness. Not fitting the average stall width increased the odds of being lame 3.7 times in primiparous cows. A narrow feed alley [<430 cm; odds ratio (OR) = 1.9], obstructed lunge space (OR = 1.7), a low BCS (OR = 2.1 for BCS ≤2.25 compared with BCS 2.75–3.0), and presence of hock lesions (OR = 1.6) were also identified as important risk factors for lameness. Only 1 of 36 farms had stalls of adequate width and length for the cows on their farm. For lameness prevention, it can be concluded that more emphasis needs be placed on either building stalls of appropriate width or selecting for smaller-framed cows that fit the existing stalls.

      Key words

      Introduction

      Lameness remains one of the most important concerns to animal health and welfare within the dairy industry. Lameness is considered to be painful and is associated with reduced milk yields, negative effects on reproductive performance, and higher culling rates, resulting in substantial production losses (
      • Huxley J.N.
      Impact of lameness and claw lesions in cows on health and production.
      ). One key to successful milking on farms with automatic milking systems (AMS) is voluntary attendance at the milking robot. Prevention of lameness is therefore important in herds with AMS also from an economic point of view, as lame cows voluntarily approach the AMS less frequently (
      • Borderas T.F.
      • Fournier A.
      • Rushen J.
      • de Passillé A.M.B.
      Effect of lameness on dairy cows’ visits to automatic milking systems.
      ;
      • Miguel-Pacheco G.G.
      • Kaler J.
      • Remnant J.
      • Cheyne L.
      • Abbott C.
      • French A.P.
      • Pridmore T.P.
      • Huxley J.N.
      Behavioural changes in dairy cows with lameness in an automatic milking system.
      ) and additional labor is required to fetch lame cows for milking (
      • Bach A.
      • Dinarés M.
      • Devant M.
      • Carré X.
      Associations between lameness and production, feeding and milking attendance of Holstein cows milked with an automatic milking system.
      ).
      An important risk factor for lameness in freestall systems with parlor milking is a lack of cow comfort, resulting in increased time spent standing by the cows, which is especially problematic if flooring is hard (
      • Bell N.J.
      • Bell M.J.
      • Knowles T.G.
      • Whay H.R.
      • Main D.J.
      • Webster A.J.F.
      The development, implementation and testing of a lameness control programme based on HACCP principles and designed for heifers on dairy farms.
      ). This lack of cow comfort could be due to inadequate stall design or bedding management. These risk factors are likely to be present also in AMS herds. However, some potential lameness risk factors are specific to AMS herds, such as overstocking of cows at the AMS or management of cow traffic, as these also affect cows’ time budget (
      • Lexer D.
      • Hagen K.
      • Palme R.
      • Troxler J.
      • Waiblinger S.
      Time budgets and adrenocortical activity of cows milked in a robot or a milking parlour: interrelationships and influence of social rank.
      ;
      • Helmreich S.
      • Hauser R.
      • Jungbluth T.
      • Wechsler B.
      • Gygax L.
      Time-budget constraints for cows with high milking frequency on farms with automatic milking systems.
      ). Studies reporting on prevalence of lameness and risk factors for lameness in AMS herds are, however, limited. The objective of the current study was to describe the prevalence of lameness in high-producing cows on AMS farms within Canada and Michigan, and to identify potential risk factors for lameness related to animal-based measures and to facility design and management practices on AMS farms.

      Materials and Methods

      This study was part of a larger study aimed at evaluating the comfort of cows on AMS farms across Canada and Michigan (
      • Westin R.
      • Vaughan A.
      • Passillé A.M.d.
      • DeVries T.J.
      • Pajor E.A.
      • Pellerin D.
      • Siegford J.M.
      • Vasseur E.
      • Rushen J.
      Lying times of lactating cows on dairy farms with automatic milking systems and the relation to lameness, leg lesions and body condition score.
      ). The study was approved by the Institutional Animal Care Committees and Research Ethics Boards at Université Laval, the University of Guelph, the University of Calgary, and Michigan State University.

      Herd Selection

      Between June 2010 and November 2012, we visited 36 farms with AMS in Quebec (QC; n = 10), Ontario (ON; n = 10), British Columbia (BC; n = 4), and Alberta (AB; n = 5), Canada, as well as Michigan (MI; n = 7). Farms had to have operated the AMS for at least 6 mo. Farms were invited by mail to participate in the study. When letters were returned indicating the willingness of the producer to participate, they were interviewed by telephone to determine if they met the additional study inclusion criteria, which included having cows housed in freestalls in their present barn for at least 1 yr and no access to outdoor exercise area or pasture. The mean (±SD) number of lactating cows in the participating farms was 125 ± 108 (maximum 495 cows) and the mean annual milk production was 9,346 ± 773 kg [retrieved from Valacta (Sainte-Anne-de-Bellevue, QC, Canada) and CanWest DHI Herd Recording (Guelph, ON, Canada) data, available for 17 farms]. The number of AMS ranged from 1 to 8 per farm, with the majority of farms having 1 (47%) or 2 (30%) AMS units. Twenty-eight farms had a free cow traffic system to access the AMS, whereas 8 farms had a directed traffic system.

      Cow Selection

      In each herd, we selected 40 lactating Holstein focal cows for detailed cow-based measurements. Where possible, we maximized the number of focal cows that were ≤120 DIM. If the milking herd had less than 40 cows ≤120 DIM, the selection criterion was expanded above 120 DIM until a sample of 40 cows was obtained. The sample of focal cows was also chosen to reflect the ratio of primiparous to multiparous cows in the herd. In 5 herds, fewer than 40 lactating cows were available; thus, fewer cows were selected as focal animals (range = 25–39 cows for these 5 herds). Data were successfully obtained from 1,378 cows in total.

      Animal-Based Measures

      All focal cows were assessed for lameness, for the presence of lesions on the hock (tarsus joint) and knees (carpal joints), leg cleanliness, BCS, hip width, and rump height. In addition, the parity and DIM of each focal cow were retrieved from farm records. Standard operating procedures were used to collect each animal based measure as described by
      • Vasseur E.
      • Gibbons J.
      • Rushen J.
      • de Passillé A.M.
      Development and implementation of a training program to ensure high repeatability of body condition scoring of dairy cows.
      ;
      • Vasseur E.
      • Gibbons J.
      • Rushen J.
      • Pellerin D.
      • Pajor E.
      • Lefebvre D.
      • de Passillé A.M.
      An assessment tool to help producers improve cow comfort on their farms.
      ), and these are available on the Canadian Dairy Research Portal (https://www.dairyresearch.ca/animal-comfort-tool.php). Measures were taken while cows were locked in self-locking head gates at the feed alley. Cows were scored for hock and knee injuries according to the criteria described by
      • Gibbons J.
      • Vasseur E.
      • Rushen J.
      • de Passillé A.M.
      A training programme to ensure high repeatability of injury scoring of dairy cows.
      and
      • Zaffino Heyerhoff J.C.
      • LeBlanc S.J.
      • DeVries T.J.
      • Nash C.G.R.
      • Gibbons J.
      • Orsel K.
      • Barkema H.W.
      • Solano L.
      • Rushen J.
      • de Passillé A.M.
      • Haley D.B.
      Prevalence of and factors associated with hock, knee, and neck injuries on dairy cows in freestall housing in Canada.
      . Both the left and right limbs were scored for hock and knee injuries on a 4-point scale (Table 1). Cleanliness on the lateral lower hind right leg (from the top of the claw to the middle of the hock) was assessed using a 4-point scoring system adapted from . Scores were based on the level of manure contamination where 0 = fresh manure for <50% of the area; 1 = fresh manure for >50% of the area; 2 = dried caked and fresh manure for <50% of the area; and 3 = entire area with dried caked manure. Body condition scores were recorded on a 5-point scale in 0.25 increments using the Elanco Animal Health body condition scoring chart for dairy cattle [
      Elanco Animal Health
      ; based on
      • Wildman E.E.
      • Jones G.M.
      • Wagner P.E.
      • Boman R.L.
      • Troutt Jr., H.F.
      • Lesch T.N.
      A dairy cow body condition scoring system and its relationship to selected production characteristics.
      and
      • Ferguson J.D.
      • Galligan D.T.
      • Thomsen N.
      Principal descriptors of body condition score in Holstein cows.
      , adapted by
      • Vasseur E.
      • Gibbons J.
      • Rushen J.
      • de Passillé A.M.
      Development and implementation of a training program to ensure high repeatability of body condition scoring of dairy cows.
      ]. Claw length was assessed by estimating the angle of the dorsal surface of the left and right lateral claws in relation to the ground. We defined claw length as no overgrowth (angle ≥45°) or overgrown (angle <45°). Hip width was measured between the points of the 2 hook bones with a flexible measuring tape when the cow was in a standing position. Rump height was measured from the ground to the spine parallel to the hook bone using a height stick. To identify lame cows, all cows were individually video recorded while walking in the feeding alley after they were released, one at a time, from the head locks. Cows were later assessed using the video to determine whether they had a head bob or had an obvious limp, defined as uneven weight bearing of one or more limbs, adapted from
      • Flower F.C.
      • Weary D.M.
      Effect of hoof pathologies on subjective assessments of dairy cow gait.
      .
      Table 1Scoring scale for carpal and tarsal joint injury in dairy cattle (
      • Gibbons J.
      • Vasseur E.
      • Rushen J.
      • de Passillé A.M.
      A training programme to ensure high repeatability of injury scoring of dairy cows.
      )
      ItemScore 0Score 1Score 2Score 3
      Tarsal jointNo swelling. No hair is missing. Thinning of hair or broken hair.No swelling or minor swelling (<1 cm). Bald area on hock.Medium swelling (1–2.5 cm) or lesion on bald area.Major swelling (>2.5 cm). May have bald area/lesion.
      Carpal jointNo skin change.Hairless patch.Lesion or scab with or without medium swelling (<2.5 cm). May have hairless patch.Major swelling (>2.5 cm) with or without lesion or hairless patch.
      Data were collected during a single visit by 2 assessors. Each assessor took the same measures on all cows. In total, 11 assessors collected data on the different farms. To ensure high repeatability, all assessors had taken part in the same training program before the study started, which is described in detail by
      • Gibbons J.
      • Vasseur E.
      • Rushen J.
      • de Passillé A.M.
      A training programme to ensure high repeatability of injury scoring of dairy cows.
      and
      • Vasseur E.
      • Gibbons J.
      • Rushen J.
      • de Passillé A.M.
      Development and implementation of a training program to ensure high repeatability of body condition scoring of dairy cows.
      . Briefly, the assessors underwent an intensive training program and the repeatability for each assessor was assessed against 2 standard trainers, with only assessors that reached the target level of a weighted kappa (Kw) ≥0.6 being used (
      • Gibbons J.
      • Vasseur E.
      • Rushen J.
      • de Passillé A.M.
      A training programme to ensure high repeatability of injury scoring of dairy cows.
      ). Once field data collection was in progress, the trainees were reassessed twice (via a refresher course 3 to 4 wk after initial training and a midway assessment) to ensure that they remained objective and repeatable in their scoring.

      Farm-Level Measurements

      Pen Design and Stocking Density

      Farm-level measures were taken in all pens containing focal cows. In 34 farms, all focal cows were located in a single pen. On the remaining 2 farms, focal cows were spread across 2 and 3 pens. In each pen, overall pen area (m2/cow) was defined as the total length by the width of the pen divided by the number of cows in the pen. Stall stocking density was calculated as the number of cows per useable lying stall. The length of the feed bunk was measured, and feeding space per cow was calculated by dividing the feed bunk space by the number of cows in the pen. The type of feed barrier and flooring type in the cow alley adjacent to the feed bunk were recorded. Cleanliness of the feed alley floor was assessed by walking the entire length of the feed bunk alley 20 min before and after scraping and measuring the height of manure that collected on the heel of a clean rubber boot. Floor cleanliness was scored as clean (<1 cm of manure) or dirty (≥1 cm of manure). If the scraping system was manual, feed alley cleanliness was scored at the beginning and at the end of the visit. The width of the alley where cows stood to feed was measured as the distance from the feed bunk to the rear curb or the front of the lying stalls.

      Stall Design and Bedding Maintenance

      Data on 8 dimensions per stall, including stall length, bed length (distance from rear curb to brisket board), brisket board height, neck rail height, curb height, distance of neck rail from rear curb, height of upper edge of bottom divide rail, and lunge space, were estimated for an average of 6 ± 1 stalls per pen (range = 4 to 10 stalls). The measures were obtained from the first and the last usable stalls in 3 representative rows in each pen containing focal cows. If <3 rows were present, the middle stall in each row was also assessed. In the same stalls, bedding quantity, cleanliness, and dryness were also assessed. Stall width was measured as the average width of 3 consecutive stalls on either side of the middle stall of each row that was measured (range 9 to 24 stalls per pen).
      To take these measures, we again used the standard operating procedures described by
      • Zaffino Heyerhoff J.C.
      • LeBlanc S.J.
      • DeVries T.J.
      • Nash C.G.R.
      • Gibbons J.
      • Orsel K.
      • Barkema H.W.
      • Solano L.
      • Rushen J.
      • de Passillé A.M.
      • Haley D.B.
      Prevalence of and factors associated with hock, knee, and neck injuries on dairy cows in freestall housing in Canada.
      and
      • Vasseur E.
      • Gibbons J.
      • Rushen J.
      • Pellerin D.
      • Pajor E.
      • Lefebvre D.
      • de Passillé A.M.
      An assessment tool to help producers improve cow comfort on their farms.
      , which are available on the Canadian Dairy Research Portal (www.dairyresearch.ca/animal-comfort-tool.php). In short, lunge space was assessed as adequate if no obstruction was present in the bob-zone, ≤76 cm forward from the top of the brisket board as proposed by
      • Nordlund K.
      • Cook N.B.
      A flowchart for evaluating dairy cow freestalls.
      . If no brisket board was present, lunge space was measured horizontally below the neck rail and 76 cm forward, 10 cm above the stall surface. Bedding quantity was evaluated as none (unable to rake bedding), ≤2 cm (when raked flat for organic bedding, or below the curb for sand bedding), or >2 cm (when raked flat for organic bedding, or even with or above the curb for sand bedding). Bedding cleanliness was assessed on a 5-point scale, where 0 = none; 1 = little manure or visible wet areas; 2 = manure-free area larger than contaminated area; 3 = contaminated area larger than manure-free area; and 4 = entire area contaminated. Dryness was scored as 0 = dry, 1 = wet, and 2 = very wet. The type of bedding and stall base (surface under bedding) were also recorded. If different types of stall bases were used in the same pen, the number of stalls with each base type was counted.

      General Management

      During the visit, all producers were interviewed about their general management routines including, for example, stall and pen management and herd health issues such as lameness monitoring and hoof trimming practices. They were also asked how many times per day cows were fetched to be milked. The complete management questionnaire is available on the Canadian Dairy Research Portal (https://www.dairyresearch.ca/animal-comfort-tool.php).

      Data Handling and Statistical Analysis

      Data Handling

      Data were entered into a relational database (Access 2010; Microsoft Corp., Redmond, WA), and then exported to Excel files (Excel 2010; Microsoft Corp.) and into Stata/IC 12 (StataCorp LP, College Station, TX) for statistical analyses. Due to missing data, the number of cows on each farm varied from measure to measure (Table 2). Data on parity and knee lesions were only obtained from 34 and 33 herds, respectively. If a categorical predictor had too few observations in any category, categories were combined based on biological cut-points or to equalize the number of observations across categories. All farm-level categorical predictors considered for analyses had at least 5 pens or farms per category (Table 3). Days in milk were categorized in quartiles. Parity was categorized as 1 and ≥2. Feed alley width was dichotomized as <430 and ≥430 cm, as a width of ≥430 cm is recommend in the Canadian Dairy Code of Practice (

      Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

      ). In alignment with these recommendations, feeding space was dichotomized as <60 and ≥60 cm. A cow was classified to not fit the stall width if the average stall width was narrower than 2× the cow’s hip width, as proposed by
      • Ceballos A.
      • Sanderson D.
      • Rushen J.
      • Weary D.M.
      Improving stall design: Use of 3-D kinematics to measure space use by dairy cows when lying down.
      and
      • Anderson N.
      . A cow was assessed to not fit the bed length if the average bed length was narrower than 1.2× the cow’s rump height as proposed by
      • Anderson N.
      . Cows were categorized as lame if they had an obvious limp, and severe lameness was categorized as having both an obvious limp and a head bob. Cows were classified as having a hock or a knee lesion if they had a score of ≥2 on at least one leg.
      Table 2Distribution of all cow-level explanatory variables hypothesized to be associated with lameness as measured on 1,378 cows from 36 dairy farms with automated milking systems
      VariableClassificationCows, n
      Does not always equal 1,387 cows per variable because of missing observations.
      (%)
      Parity
      Data only obtained from 34 farms.
      1425 (34)
      ≥2838 (66)
      DIM≤57319 (25)
      58–98316 (25)
      99–163323 (25)
      ≥165317 (25)
      BCS≤2.25185 (13)
      2.5375 (27)
      2.75–3.0451 (33)
      ≥3.25364 (27)
      Overgrown clawsNo718 (66)
      Yes373 (34)
      Leg cleanlinessClean1,150 (84)
      Dirty (score ≥2)225 (16)
      Hock lesionNo929 (70)
      Yes (score ≥2)398 (30)
      Knee lesion
      Data only obtained from 33 farms.
      No896 (73)
      Yes (score ≥2)329 (27)
      Fits stall widthCow fits503 (38)
      Cow does not fit808 (62)
      Fits bed lengthCow fits759 (58)
      Cow does not fit554 (42)
      1 Does not always equal 1,387 cows per variable because of missing observations.
      2 Data only obtained from 34 farms.
      3 Data only obtained from 33 farms.
      Table 3Distribution of all pen and farm-level explanatory variables hypothesized to be associated with lameness as measured on 1,378 cows in 39 pens from 36 dairy farms with automated milking systems
      VariableClassificationPens, n
      Does not always equal 39 pens per variable because of missing observations.
      (%)
      Herd, n (%)MeanSD
      Stall baseRubber or geotextile mat19 (49)
      Sand10 (26)
      Other10 (26)
      Bedding typeSand10 (26)
      Straw9 (23)
      Sawdust8 (29)
      Wood shavings6 (15)
      None5 (13)
      Bedding quantity≥2 cm213 (33)
      None or <2 cm326 (67)
      Bedding cleanliness
      Not tested in statistical analysis because of too few observations in one category.
      Clean (score 0–1)
      In ≥75% of measured stalls.
      38 (97)
      Dirty (score ≥2)1 (3)
      Bedding dryness
      Not tested in statistical analysis because of too few observations in one category.
      Dry
      In ≥75% of measured stalls.
      35 (97)
      Wet (score ≥1)1 (3)
      Lunge spaceAdequate
      In ≥75% of measured stalls.
      29 (74)
      Obstructed
      In >25% of measured stalls.
      10 (26)
      Raking stall frequency
      Not tested in statistical analysis because of too few observations in one category.
      ≥Once a day35 (88)
      <Once a day4 (10)
      Type of feed barrierPost-and-rail13 (33)
      Headlock or diagonal bars22 (56)
      Mix4 (10)
      Feeding space<60 cm/cow21 (54)
      ≥60 cm/cow18 (46)
      Floor at feed bunkConcrete25 (64)
      Rubber14 (36)
      Width of feed alley<430 cm27 (69)
      ≥430 cm12 (31)
      Feed alley cleanliness 1 (before scraping)Clean (score 0–1)31 (86)
      Dirty (score 2–3)8 (21)
      Feed alley cleanliness 2
      Not tested in statistical analysis because of too few observations in one category.
      (after scraping)
      Clean (score 0–1)37 (100)
      Dirty (score 2–3)0 (0)
      AMS systemFree traffic28 (78)
      Forced traffic8 (22)
      Herd size<100 cows20 (56)
      ≥100 cows16 (44)
      Stall dimensions
       Widthcm391204
       Lengthcm3924314
       Bed lengthcm3917810
       Brisket board heightcm24144
       Height of neck railcm391187
       Distance of neck rail from rear curbcm3917111
       Curb heightcm39204
      No. of cows per robotNo. of cows39549
      Stocking densityCows per stall390.930.16
      Overall pen aream2/cow3993
      1 Does not always equal 39 pens per variable because of missing observations.
      2 Not tested in statistical analysis because of too few observations in one category.
      3 In ≥75% of measured stalls.
      4 In >25% of measured stalls.
      Lunge space was classified as adequate at the pen level if no obstruction was present in ≥75% of the measured stalls. Similarly, bedding quantity, cleanliness, and dryness was classified as ≥2 cm, clean (score 0–1), and dry (score 0) if ≥75% of the measured stalls in the pen received the corresponding scores.

      Statistical Analysis

      Model building was performed in 2 steps. The outcome of interest was the presence of lameness (i.e., having an obvious limp) at the cow level. First, simple logistic regression analysis was performed to assess the association between lameness and each predictor variable, applying the LOGIT procedure in Stata (StataCorp LP). Hypothesized and tested variables are listed in Tables 2 and 3. Predictors with an association with lameness at P ≤ 0.25 were considered for further modeling (listed in Table 4).
      Table 4Univariable associations expressed as odds ratios of the presence of lameness with individual and herd-level factors in 1,378 cows on 36 dairy farms with automated milking system
      Results of univariable logistic regression models. Only associations at P≤0.25 are shown.
      Variablen
      Does not always equal 1,378 cows per variable because of missing observations.
      Odds ratioSEP-value
      Individual measures
      BCS1,375
       2.75–3.01.0
      Reference category.
       2.0–2.252.420.59<0.001
       2.51.270.280.26
       3.25–4.50.870.200.55
      Does not fit stall width
      Cow does not fit if average stall width <2× cow hip width.
      1,3111.790.31<0.001
      Leg dirty (score ≥2)1,3751.340.260.14
      Parity (≥2)
      Data from 34 farms only.
      1,2631.720.320.002
      Presence of hock lesion1,3271.480.240.02
      Presence of knee lesion
      Data from 33 farms only.
      1,2251.840.31<0.001
      Housing and management
       Bedding quantity
      In ≥75% of measured stalls.
      (≥2 cm)
      1,3780.720.120.06
       Brisket board height (per cm increase)8880.960.020.12
       Herd size (≥100 cows)1,3781.200.180.23
       Feeding space (<60 cm/cow)1,3781.240.190.16
       Narrow feed alley (<430 cm)1,3781.520.270.01
       No cows per robot (per 10 cows increase)1,3781.150.100.09
       Obstructed lunge space
      In >25% of measured stalls.
      1,3781.540.250.01
       Overall pen area (per 1 m2/cow increase)1,3780.930.030.005
       Raking stall frequency (<once a day)1,3781.580.380.06
      Stall base1,378
       Rubber or geotextile mat1.0
      Reference category.
       Sand0.670.140.06
       Other1.180.200.34
      Stall stocking density (per 10% increase)1,3781.120.050.02
      1 Results of univariable logistic regression models. Only associations at P ≤ 0.25 are shown.
      2 Does not always equal 1,378 cows per variable because of missing observations.
      3 Reference category.
      4 Cow does not fit if average stall width <2× cow hip width.
      5 Data from 34 farms only.
      6 Data from 33 farms only.
      7 In ≥75% of measured stalls.
      8 In >25% of measured stalls.
      Next, we built 2 separate multivariable mixed-effects logistic regression models with the selected predictor variables (Table 4). One model considered cow-level measures and the other examined management and explanatory variables at the farm and pen level. As the focal cows were nested within pen and within farm, the random effects of pen within farm were forced into the models using the XTMELOGIT procedure in Stata (StataCorp LP). Both models were constructed using manual backward stepwise elimination, retaining variables with P ≤ 0.05. Continuous variables were centered on the overall mean and tested with their quadratic terms. The variable stall base was tested using only 2 categories (sand vs. all others) as no difference was noted between the non-sand categories in the univariable analyses. In addition, the fixed effect of geographic location (AB and BC, ON, QC, and MI) was tested in both models but not retained (P > 0.05). Finally, previously excluded variables were forced into the model again and retained if P ≤ 0.05 or if they were judged as confounders (changed the estimate of any other variable by >30%). Biologically plausible 2-way interactions between the variables retained in the models were tested. An interaction between the variables parity and fit stall width was found to be significant (P ≤ 0.05) and was therefore kept in the final model concerning cow-level measures. Coefficient estimates were transformed and are presented as odds ratios (OR). The final models were validated by examination of influential covariate patterns and by running the model with and without these. We found no reason to exclude any observations, and the model fits were considered satisfactory.

      Results

      Lameness Prevalence

      The mean herd prevalence of lameness was 15%, but ranged from 2.5 to 46% among the participating farms. Severe lameness was less common, averaging 4%; however, in 5 herds, 10 to 12% of the focal cows were severely lame (i.e., obviously limping and having a head bob). All producers with >20% lameness in their herd (8 farms) stated that they had a major or moderate problem with lameness when answering the questionnaire. However, 7 out of 8 producers with a herd prevalence between 15 and 20% answered that they had a minor problem with lameness. Of the 36 participating farms, 10 kept records of all identified cases of lameness, whereas on 13 farms treated cases only were recorded. Four farms used hoof trimmer records as the single documentation of lameness. The remaining 9 farms did not keep any record of lameness at all. Almost all farms trimmed their cows routinely (≥1 time per year). Only 2 farmers stated that they called the trimmer out or trimmed themselves only when a cow needed it. Producers stated that they fetched cows to the AMS unit on average 3 times/d (range = 0–12). A positive correlation was noted between the herd prevalence of lameness and the number of times per day producers stated that cows were fetched to be milked (r = 0.45, P = 0.007).

      Cow-Level Factors

      Results of single variable analyses (P ≤ 0.25) are presented in Table 4 and results of the multivariable analysis in Table 5. Body condition score, presence of hock lesion, and the interaction between parity and not fitting the average stall width were associated with lameness in the final model concerning cow-level factors (Table 5). A BCS of ≤2.25 was associated with twice the odds of being lame (OR = 2.11, P = 0.005) compared with a BCS of 2.75 to 3.0. In total, 13% of the cows had a BCS ≤2.25 but the proportion ranged from 0 to 40% across the herds. Farms situated in QC had, on average, 29% of their focal cows within BCS category ≤2.25, compared with 10% in ON and MI and 3% in AB and BC. Only 19 cows (1%) had a BCS of 2.0. The proportion of lame cows within different BCS categories is presented in Figure 1.
      Table 5Results of a multivariable mixed-effects logistic regression model
      Farm and pen included as random effects.
      used to investigate the effect of animal-based measures on the probability of lameness in 1,173 cows in 34 dairy farms with automatic milking systems
      VariableOdds ratio95% CIP-value
      Body condition score
       2.75–3.01.0
      Reference category.
       2.0–2.252.111.25–3.560.005
       2.51.140.71–1.830.60
       3.25–4.50.660.38–1.160.15
      Presence of hock lesion1.631.10–2.410.014
      Parity × fit stall width
      Cow does not fit stall width if average stall width <2× cow hip width.
      interaction
       Primiparous – Fits stall width1.0
      Reference category.
       Primiparous – Does not fit stall width3.731.54–9.070.004
       Multiparous – Fits stall width3.491.41–8.630.007
       Multiparous – Does not fit stall width4.471.96–10.23<0.001
      1 Farm and pen included as random effects.
      2 Reference category.
      3 Cow does not fit stall width if average stall width <2× cow hip width.
      Figure thumbnail gr1
      Figure 1Proportion of lame cows within different BCS categories (n = 1,375 cows). A BCS of ≤2.25 was associated with twice the odds of being lame (odds ratio = 2.1, P = 0.005) compared with a BCS of 2.75 to 3.0.
      Cows had a mean (±SD) hip width of 62 ± 7 cm (n = 1,311) and a mean rump height of 148 ± 5 cm (n = 1,313). The width of the stalls was 120 ± 4 cm, on average, and mean bed length was 178 ± 10 cm (Table 3). As a result, 62% of the cows were too big to fit the average stall width (2× the cow’s hip width needed) and 42% did not fit the average bed length (1.2× the cow’s rump height needed). Only 30% of the multiparous cows fit the stall width. Due to a significant interaction with parity, not fitting the average stall width increased the odds of being lame 3.7 times in primiparous (P = 0.004) and 1.3 times in multiparous cows (P < 0.001; Table 5; Figure 2). We found no association between lameness and the cow not fitting the stall bed length (P = 0.83). Out of 36 farms, 10 farms succeeded in fitting ≥75% of the focal animals to stall width and 19 farms in fitting cows to bed length. Only 1 farm had stalls wide and long enough to fit all focal cows. It was cow size (hip width), rather than stall width, that varied across farms with a different proportion of cows fitting the stall width (Table 6).
      Figure thumbnail gr2
      Figure 2Proportion of lame cows with different parity presented in relation to whether the cow fits the average stall width or not (n = 1,215 cows). Not fitting the average stall width increased the odds of being lame by 3.7 times in primiparous (P = 0.004) and 1.3 times in multiparous cows (P < 0.001).
      Table 6Description of stall width, bed length and cow size in 36 farms with different proportions of cows fitting the stalls
      ItemNo. of

      farms
      Mean stall

      width ± SD
      Mean hip

      width ± SD
      Mean difference

      (stall width – 2× cow’s hip width)
      Proportion of cows fitting mean stall width
       <25%17119 ± 466 ± 2−13 ± 4
       25–49%7122 ± 264 ± 1−5 ± 1
       50–74%2116 ± 259 ± 2−1 ± 2
       75–100%10121 ± 354 ± 113 ± 3
      Proportion of cows fitting mean bed length
       <25%11167 ± 4149 ± 2−12 ± 5
       25–49%2174 ± 2148 ± 1−4 ± 1
       50–74%4182 ± 1150 ± 12 ± 1
       75–100%19186 ± 8148 ± 38 ± 7

      Farm-Level Factors

      Of the housing and management predictors tested, a narrow feed alley (<430 cm; OR = 1.85; P = 0.006) and an obstructed lunge space (OR = 1.71, P = 0.02) were associated with higher odds of lameness in the final multivariable model (Table 7). The width of the feed alley averaged 406 cm (range = 301–640 cm across pens), and the proportion of lame cows within pens with different feed alley widths is presented in Figure 3. Twenty-seven farms had adequate lunge space in ≥75% their measured stalls. The mean prevalence of lameness in these farms was 13% compared with 19% in farms that had more stalls with obstructed lunge spaces.
      Table 7Results of a multivariable mixed-effects logistic regression model
      Farm and pen included as random effects.
      used to investigate the effect of farm-level factors on the probability of lameness in 1,378 cows in 36 dairy farms with automatic milking systems
      VariableOdds ratio95% CIP-value
      Narrow feed alley (<430 cm)1.851.19–2.870.006
      Obstructed lunge space
      In >25% of measured stalls.
      1.671.08–2.600.021
      Stall base, sand vs. all others0.630.38–1.030.065
      1 Farm and pen included as random effects.
      2 In >25% of measured stalls.
      Figure thumbnail gr3
      Figure 3Proportion of lame cows in 39 pens from 36 dairy farms with automatic milking systems, presented in relation to width of the feed alley. Dotted line shows minimum feed alley width recommended in the Canadian Dairy Code of Practice (

      Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

      ).
      Stall base was associated with lameness in the screening of single variables (P = 0.04), and a tendency was noted for sand bedding to be associated with lower odds of lameness in the final multivariable model (P = 0.06; Table 7). More pen space per cow was associated with decreased odds, and increasing stall stocking density was associated with higher odds for lameness in the simple analysis only (Table 4). On 10 farms, stocking density was higher than 1 cow/lying stall, which is the recommended best practice according to the Canadian Dairy Code of Practice (

      Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

      ). Only 1 farm exceeded the Canadian Dairy Code requirement of a maximum stocking density of 1.2 cows/stall.
      Neither the AMS traffic system used nor the number of cows per AMS unit was related to lameness. Half of the farms provided the recommended 60 cm of feeding space for each cow, but higher feed bunk density was not associated with lameness.

      Discussion

      Prevalence of Lameness in AMS Farms

      The mean prevalence of lameness in our study was 15%. This is similar to the 14% reported for AMS farms in Denmark (
      • Klaas I.C.
      • Rousing T.
      • Foussing C.
      • Hindhede J.
      • Sørensen J.T.
      Is lameness a welfare problem in dairy farms with automatic milking systems?.
      ), but lower than lameness prevalence on 141 Canadian freestall farms with a milking parlor (21%;
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Nash C.G.R.
      • Haley D.B.
      • Vasseur E.
      • Pellerin D.
      • Rushen J.
      • de Passillé A.M.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ) and 100 Canadian tiestall farms (24%;
      • Charlton G.L.
      • Bouffard V.
      • Gibbons J.
      • Vasseur E.
      • Haley D.B.
      • Pellerin D.
      • Rushen J.
      • de Passillé A.M.
      Can automated measures of lying time help assess lameness and leg lesions on tie-stall dairy farms?.
      ). The latter 2 studies followed the same data collection protocol as described in the current paper. Our results, thus, indicate that prevalence of lameness seems to be lower in AMS herds in comparison to Canadian tiestall and freestall herds without AMS. This may be due to factors associated with the AMS, such as desynchronized milking times, or it may reflect the management styles of producers who adopt AMS, including a greater awareness of the importance of controlling lameness. On the other hand, visited AMS farms were few and it is possible that these herds had a low prevalence of lameness before transition to AMS. Also, cows on visited AMS farms were often housed in newer facilities compared with the majority of participating freestall farms without AMS. In other studies, differences in locomotion scores associated with converting to AMS have not been found when holding management and facility design constant (
      • Hillerton J.E.
      • Dearing J.
      • Dale J.
      • Poelarends J.J.
      • Neijenhuis F.
      • Sampimon O.C.
      • Miltenburg J.D.H.M.
      • Fossing C.
      Impact of automatic milking on animal health.
      ;
      • Vosika B.
      • Lexer D.
      • Stanek C.
      • Troxler J.
      • Waiblinger S.
      The influence of an automatic milking system on claw health and lameness of dairy cows.
      ).
      In the present study we found a positive correlation between the herd prevalence of lameness and the number of times the farmer stated that cows had to be fetched for milking. The number of times cows were fetched was only stated by the farmers, but not measured directly. However, we have no reason to believe that farmers did not answer this question truthfully with estimates close to the actual number; therefore, we consider these data valid for concluding there is a correlation between lameness and fetching. Fetching cows is very time consuming, and the increased workload possibly makes herd managers more motivated to reduce lameness in AMS compared with non-AMS farms. Some farms in the current study were indeed very successful at controlling lameness, having ≤5% lame cows. This demonstrates that a lameness prevalence of 5% or less is an achievable goal.

      Cow-Level Factors Associated with Lameness

      Not fitting the average stall width was the strongest animal-based risk factor for lameness, increasing the odds of lameness 3.7 times in primiparous cows. Difficulties in performing lying down movements may make the cow more hesitant to lie down, possibly resulting in greater standing time. It is also possible that rest is disrupted by neighboring cows more often if stalls are too narrow. A 3-dimensional biomechanical study by
      • Ceballos A.
      • Sanderson D.
      • Rushen J.
      • Weary D.M.
      Improving stall design: Use of 3-D kinematics to measure space use by dairy cows when lying down.
      revealed that the lateral displacement of a cow during lying down movements reaches up to 180% of the hip width. A stall width that is 200% of the hip width is therefore proposed as adequate to accommodate most lying movements (
      • Ceballos A.
      • Sanderson D.
      • Rushen J.
      • Weary D.M.
      Improving stall design: Use of 3-D kinematics to measure space use by dairy cows when lying down.
      ;
      • Anderson N.
      ).
      • Tucker C.B.
      • Weary D.M.
      • Fraser D.
      Free-stall dimensions: Effects on preference and stall usage.
      found that cows reduced their lying time by 1.2 h/d in narrow stalls (112 vs. 132 cm wide). Reduced lying time acts as an exacerbating factor in the development of claw lesions (
      • Leonard F.C.
      • O'Connell J.M.
      • O'Farrell K.J.
      Effect of overcrowding on claw health in first-calved friesian heifers.
      ), and poor cow comfort resulting in reduced lying time is believed to have a critical role in both the development of lameness and the speed of recovery (
      • Cook N.B.
      • Nordlund K.V.
      The influence of the environment on dairy cow behavior, claw health and herd lameness dynamics.
      ).
      Stalls in our study were, on average, 120 cm wide and the majority of both primiparous and multiparous cows were too large to fit these properly. In the Canadian Dairy Code of Practice, 117 cm is the recommended minimum stall width for cows weighing 545 kg, 122 cm for cows weighing 636 kg, and 127 cm for 727 kg cows (

      Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

      ). In 2014, Quebec herd-average cow weight was 664 kg (including first lactation cows at 619 kg;

      Valacta. 2015. L'évolution de la production laitière Québécoise 2014. Quebec. Accessed Aug. 21, 2015. www.valacta.com/FR/Nos-publications/Documents/%C3%89VOLUTION%20LAITI%C3%88RE/%C3%89volution%202014_FINAL.pdf

      ). The fact that cows too large for their stalls in our study were at greater risk of being lame suggests that the guidelines within the Canadian Dairy Code of Practice (

      Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

      ) are reasonably accurate. Experimental studies are needed to explore the association between lameness and stall width in detail. As most farms failed to provide stalls of adequate width to accommodate the size of their cows, correct body size does not appear to be taken into account when building stalls. This could be addressed by constructing wider stalls or, alternatively, by selecting for smaller-framed cows that would fit into the pre-existing stalls.
      Lower BCS (≤2.25) were associated with increased odds for lameness. A correlation between low BCS and lameness has been demonstrated in numerous studies (reviewed by
      • Huxley J.N.
      Impact of lameness and claw lesions in cows on health and production.
      ). Historically, it was assumed that lame cows lost weight as a consequence of reduced feed intake. A growing number of studies, however, suggest that a low BCS may predispose cattle to lameness (
      • Green L.E.
      • Huxley J.N.
      • Banks C.
      • Green M.J.
      Temporal associations between low body condition, lameness and milk yield in a UK dairy herd.
      ;
      • Lim P.Y.
      • Huxley J.N.
      • Willshire J.A.
      • Green M.J.
      • Othman A.R.
      • Kaler J.
      Unravelling the temporal association between lameness and body condition score in dairy cattle using a multistate modelling approach.
      ;
      • Randall L.V.
      • Green M.J.
      • Chagunda M.G.G.
      • Mason C.
      • Archer S.C.
      • Green L.E.
      • Huxley J.N.
      Low body condition predisposes cattle to lameness: An 8-year study of one dairy herd.
      ). Body condition score is positively correlated with the thickness of the claw’s digital cushion, which in turn is a strong predictor of lameness (
      • Bicalho R.C.
      • Machado V.S.
      • Caixeta L.S.
      Lameness in dairy cattle: A debilitating disease or a disease of debilitated cattle? A cross-sectional study of lameness prevalence and thickness of the digital cushion.
      ). Also, evidence exists that thinner cows may have a lower social rank and, therefore, are less able to compete for lying space and access to feed (
      • Hohenbrink S.
      • Meinecke-Tillmann S.
      Influence of social dominance on the secondary sex ratio and factors affecting hierarchy in Holstein dairy cows.
      ). Cows in AMS herds also have to compete for access to the milking robot and low-ranking cows in general spend more time queuing in the waiting area (
      • Melin M.
      • Hermans G.G.N.
      • Pettersson G.
      • Wiktorsson H.
      Cow traffic in relation to social rank and motivation of cows in an automatic milking system with control gates and an open waiting area.
      ;
      • Lexer D.
      • Hagen K.
      • Palme R.
      • Troxler J.
      • Waiblinger S.
      Time budgets and adrenocortical activity of cows milked in a robot or a milking parlour: interrelationships and influence of social rank.
      ). Simulations of waiting time show that low-ranking cows may wait to milk about 1 h on a typical day, compared with 10 min for a middle-rank animal and 3.5 min for the dominant ones (
      • Halachmi I.
      Simulating the hierarchical order and cow queue length in an automatic milking system.
      ). In crowded situations, simulations show that a middle-ranking cow queues for approximately 1.5 h, whereas a low-ranking cow may wait about 7 h (
      • Halachmi I.
      Simulating the hierarchical order and cow queue length in an automatic milking system.
      ). The daily lying time of the cows participating in this study was up to 1 h shorter in individuals with BCS ≤2.25 compared with cows with BCS ≥3.25 (
      • Westin R.
      • Vaughan A.
      • Passillé A.M.d.
      • DeVries T.J.
      • Pajor E.A.
      • Pellerin D.
      • Siegford J.M.
      • Vasseur E.
      • Rushen J.
      Lying times of lactating cows on dairy farms with automatic milking systems and the relation to lameness, leg lesions and body condition score.
      ). Currently, dairy producers in Canada are required to take corrective action for animals at a BCS of 2 or lower (

      Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

      ). Only 1% of examined cows fell within this category; thus, our results show that at a BCS of 2.25, animals are at higher odds of being lame. To take corrective action sooner is therefore warranted.
      An association between lameness and presence of hock lesions has been found in several studies of freestall systems with milking parlors (
      • Brenninkmeyer C.
      • Dippel S.
      • Brinkmann J.
      • March S.
      • Winckler C.
      • Knierim U.
      Hock lesion epidemiology in cubicle housed dairy cows across two breeds, farming systems and countries.
      ;
      • Kester E.
      • Holzhauer M.
      • Frankena K.
      A descriptive review of the prevalence and risk factors of hock lesions in dairy cows.
      ;
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Nash C.G.R.
      • Haley D.B.
      • Vasseur E.
      • Pellerin D.
      • Rushen J.
      • de Passillé A.M.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). As risk factors and the outcome of interest were measured at the same time in our study, we cannot draw conclusions about causality. The presence of hock lesions is, however, strongly related to time spent lying on abrasive surfaces and prolonged high local pressure or friction of the hock on hard surfaces, and it is suggested that lameness aggravates the risk of hock lesions rather than the other way around (
      • Kester E.
      • Holzhauer M.
      • Frankena K.
      A descriptive review of the prevalence and risk factors of hock lesions in dairy cows.
      ).

      Housing and Management Factors Associated with Lameness

      Of the other housing and management factors tested, the width of the feed alley, lunge space, and stall base (sand vs. all others) were the most important in relation to probability of being lame. To allow easy movement of cattle in the feeding area, a feed alley width of 430 cm is recommended in the Dairy Code of Practice (

      Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

      ). Farms fulfilling this criterion in our study had fewer lame cows, suggesting that this recommendation is relevant. If walking alleys are too narrow, cows cannot move around without interference from other cows. Reducing the general space per cow in cubicle systems has been shown to increase the number of agonistic interactions among cows (
      • Menke C.
      • Waiblinger S.
      • Fölsch D.W.
      • Wiepkema P.R.
      Social behaviour and injuries of horned cows in loose housing systems.
      ;
      • Fregonesi J.A.
      • Leaver J.D.
      Influence of space allowance and milk yield level on behaviour, performance and health of dairy cows housed in strawyard and cubicle systems.
      ). Agonistic behavior is common around the feed bunk (
      • Miller K.
      • Wood-Gush D.G.M.
      Some effects of housing on the social behaviour of dairy cows.
      ); therefore, it seems likely that aggressive behaviors will occur more frequently if space in the feeding alley is limited. According to
      • Sarjokari K.
      • Kaustell K.O.
      • Hurme T.
      • Kivinen T.
      • Peltoniemi O.A.T.
      • Saloniemi H.
      • Rajala-Schultz P.J.
      Prevalence and risk factors for lameness in insulated free stall barns in Finland.
      , agonistic interactions possibly induce claw trauma as a result of rapidly moving to avoid or give way to dominant cows. In alignment with our findings,
      • Sarjokari K.
      • Kaustell K.O.
      • Hurme T.
      • Kivinen T.
      • Peltoniemi O.A.T.
      • Saloniemi H.
      • Rajala-Schultz P.J.
      Prevalence and risk factors for lameness in insulated free stall barns in Finland.
      also found an association between a narrow feed alley and an increased prevalence of lameness (<320 cm compared with >340 cm). Other recent, large studies of herd-level risk factors for lameness seem to have overlooked the width of the feed alley (
      • Barker Z.E.
      • Leach K.A.
      • Whay H.R.
      • Bell N.J.
      • Main D.C.J.
      Assessment of lameness prevalence and associated risk factors in dairy herds in England and Wales.
      ;
      • Chapinal N.
      • Barrientos A.K.
      • von Keyserlingk M.A.G.
      • Galo E.
      • Weary D.M.
      Herd-level risk factors for lameness in freestall farms in the northeastern United States and California.
      ;
      • Pérez-Cabal M.A.
      • Alenda R.
      Clinical lameness and risk factors in a Spanish Holstein population.
      ). We therefore, propose that width of the feed alley be included as a potential risk factor for lameness in future studies.
      Stall designs that fail to provide for the movements associated with lying and rising and the need of a soft, cushioned surface reduce the use of the stall and thus increase the risk of lameness (
      • Cook N.B.
      • Nordlund K.V.
      The influence of the environment on dairy cow behavior, claw health and herd lameness dynamics.
      ;
      • Bicalho R.C.
      • Oikonomou G.
      Control and prevention of lameness associated with claw lesions in dairy cows.
      ). Therefore, it is no surprise that in our study an obstructed lunge space was associated with higher odds for lameness and sand bedding in stalls was associated with lower odds for lameness. Deep-bedded sand, because of its ability to supply cushion and traction, allows cows to perform the processes of rising and lying down more easily, without fear of slipping (
      • Cook N.B.
      • Nordlund K.V.
      The influence of the environment on dairy cow behavior, claw health and herd lameness dynamics.
      ). Sand-bedded stalls also increase the daily lying time, especially in lame cows (
      • Gomez A.
      • Cook N.B.
      Time budgets of lactating dairy cattle in commercial freestall herds.
      ;
      • van Gastelen S.
      • Westerlaan B.
      • Houwers D.J.
      • van Eerdenburg F.J.C.M.
      A study on cow comfort and risk for lameness and mastitis in relation to different types of bedding materials.
      ).
      It is widely accepted that having 60 cows/AMS unit is the optimal stocking density (
      • Deming J.A.
      • Bergeron R.
      • Leslie K.E.
      • DeVries T.J.
      Associations of housing, management, milking activity, and standing and lying behavior of dairy cows milked in automatic systems.
      ). The majority of farms in our study had fewer than 60 cows/unit and the number of cows per AMS unit was not associated with lameness. The type of traffic management surrounding the AMS (free vs. directed traffic) was also not associated with lameness in our study. Only 8 of 36 farms had a directed traffic system in the current study; further, the design of the directed traffic systems was highly variable. More detailed studies are therefore needed to fully understand the effect of traffic system on lameness in AMS herds.
      Stall stocking density was only positively correlated with lameness in the univariable analysis. It has previously been demonstrated that overstocking reduces lying time in a linear fashion (
      • Fregonesi J.A.
      • Tucker C.B.
      • Weary D.M.
      Overstocking reduces lying time in dairy cows.
      ), which could explain a possible link to lameness. However, in the present study, all but one farm followed the Canadian guidelines of ≤1.2 cows/stall. Thus, data from overcrowded herds are needed to more fully evaluate the risk level for lameness due to overstocking.

      Conclusions

      The observed prevalence of lameness in our study population of AMS herds was lower than that observed previously in Canadian tiestall and freestall herds without AMS. The main risk factors for lameness in AMS farms involve stall design and, thus, appear to be similar to those on freestall farms with conventional milking parlors, rather than being related to management of the AMS. In particular, stall width relative to cow size was found to be the most important factor associated with lameness. Only 1 of 36 farms had stalls of both adequate width and length for the cows on their farm. Narrow feed alleys were also identified as a risk factor. The strong association between stall width and prevalence of lameness indicates that more emphasis needs be placed on either building stalls of appropriate width or selecting for smaller framed cows which fit the existing stalls. Experimental studies are needed to specifically explore the association between lameness and stall width.

      Acknowledgments

      We thank the participating farmers from the United States and Canada and the collaborators, technicians, students, and co-op students from Michigan State University, the University of British Columbia, University of Calgary, University of Guelph, Laval University, and Valacta Inc. (Sainte-Anne-de-Bellevue, QC, Canada). We thank Jenny Gibbons (AHDB Dairy, Warwickshire, UK) and Gemma Charlton (Harper Adams College, Newport, UK) for invaluable help with establishing the database. This study was funded by the Dairy Farmers of Canada (Ottawa, ON, Canada), the Canadian Dairy Commission (Ottawa, ON, Canada), Agriculture and Agri-Food Canada (Ottawa, ON, Canada), as part of the Dairy Science Cluster initiative, Agriculture and Agri-Food Canada and the BC Ministry of Agriculture (Victoria, BC, Canada) through the Canada-BC Agri-Innovation Program under Growing Forward 2, a federal-provincial-territorial initiative (delivered by the Investment Agriculture Foundation of BC), Alberta Milk (Edmonton, AB, Canada), the Alberta Livestock and Meat Agency (Edmonton, AB, Canada), and the Natural Sciences and Engineering Research Council of Canada (Ottawa, ON, Canada). We also thank the Michigan State University Department of Animal Science for supporting this research through an award from the Elwood Kirkpatrick Dairy Research Science Endowment to Janice Siegford (East Lansing, MI).

      References

        • Anderson N.
        Free Stall Dimensions. Infosheet - Livestock Technology. Ministry of Agriculture, Food and Rural Affairs, Fergus, Ontario, Canada2007
        • Bach A.
        • Dinarés M.
        • Devant M.
        • Carré X.
        Associations between lameness and production, feeding and milking attendance of Holstein cows milked with an automatic milking system.
        J. Dairy Res. 2007; 74 (http://dx.doi.org/10.1017/S0022029906002184): 40-46
        • Barker Z.E.
        • Leach K.A.
        • Whay H.R.
        • Bell N.J.
        • Main D.C.J.
        Assessment of lameness prevalence and associated risk factors in dairy herds in England and Wales.
        J. Dairy Sci. 2010; 93 (http://dx.doi.org/10.3168/jds.2009-2309): 932-941
        • Bell N.J.
        • Bell M.J.
        • Knowles T.G.
        • Whay H.R.
        • Main D.J.
        • Webster A.J.F.
        The development, implementation and testing of a lameness control programme based on HACCP principles and designed for heifers on dairy farms.
        Vet. J. 2009; 180 (http://dx.doi.org/10.1016/j.tvjl.2008.05.020): 178-188
        • Bicalho R.C.
        • Machado V.S.
        • Caixeta L.S.
        Lameness in dairy cattle: A debilitating disease or a disease of debilitated cattle? A cross-sectional study of lameness prevalence and thickness of the digital cushion.
        J. Dairy Sci. 2009; 92 (http://dx.doi.org/10.3168/jds.2008-1827): 3175-3184
        • Bicalho R.C.
        • Oikonomou G.
        Control and prevention of lameness associated with claw lesions in dairy cows.
        Livest. Sci. 2013; 156 (http://dx.doi.org/10.1016/j.livsci.2013.06.007): 96-105
        • Borderas T.F.
        • Fournier A.
        • Rushen J.
        • de Passillé A.M.B.
        Effect of lameness on dairy cows’ visits to automatic milking systems.
        Can. J. Anim. Sci. 2008; 88 (http://dx.doi.org/10.4141/CJAS07014): 1-8
        • Brenninkmeyer C.
        • Dippel S.
        • Brinkmann J.
        • March S.
        • Winckler C.
        • Knierim U.
        Hock lesion epidemiology in cubicle housed dairy cows across two breeds, farming systems and countries.
        Prev. Vet. Med. 2013; 109 (http://dx.doi.org/10.1016/j.prevetmed.2012.10.014): 236-245
        • Ceballos A.
        • Sanderson D.
        • Rushen J.
        • Weary D.M.
        Improving stall design: Use of 3-D kinematics to measure space use by dairy cows when lying down.
        J. Dairy Sci. 2004; 87 (http://dx.doi.org/10.3168/jds.S0022-0302(04)70022-3): 2042-2050
        • Chapinal N.
        • Barrientos A.K.
        • von Keyserlingk M.A.G.
        • Galo E.
        • Weary D.M.
        Herd-level risk factors for lameness in freestall farms in the northeastern United States and California.
        J. Dairy Sci. 2013; 96 (http://dx.doi.org/10.3168/jds.2012-5940): 318-328
        • Charlton G.L.
        • Bouffard V.
        • Gibbons J.
        • Vasseur E.
        • Haley D.B.
        • Pellerin D.
        • Rushen J.
        • de Passillé A.M.
        Can automated measures of lying time help assess lameness and leg lesions on tie-stall dairy farms?.
        Appl. Anim. Behav. Sci. 2015; (http://dx.doi.org/10.1016/j.applanim.2015.02.011, In press)
      1. Cook, N. B. 2006. Footbath alternatives. Accessed May 19, 2015. http://www.vetmed.wisc.edu/dms/fapm/fapmtools/6lame/Footbath_Alternatives.pdf

        • Cook N.B.
        • Nordlund K.V.
        The influence of the environment on dairy cow behavior, claw health and herd lameness dynamics.
        Vet. J. 2009; 179 (http://dx.doi.org/10.1016/j.tvjl.2007.09.016): 360-369
        • Deming J.A.
        • Bergeron R.
        • Leslie K.E.
        • DeVries T.J.
        Associations of housing, management, milking activity, and standing and lying behavior of dairy cows milked in automatic systems.
        J. Dairy Sci. 2013; 96 (http://dx.doi.org/10.3168/jds.2012-5985): 344-351
      2. Dairy Farmers of Canada and the National Farm Animald Care Council (DFC-NFACC). 2009. Code of practices for the care and handling of dairy cattle. Dairy Farmers of Canada, Ottawa, Ontarioa, Canada. Accessed May 10, 2015. http://www.nfacc.ca/codes-of-practice

        • Elanco Animal Health
        Body Condition Scoring. Bulletin AI 8478, Rev. 9/96. Elanco Animal Health, Indianapolis, IN1996
        • Ferguson J.D.
        • Galligan D.T.
        • Thomsen N.
        Principal descriptors of body condition score in Holstein cows.
        J. Dairy Sci. 1994; 77 (http://dx.doi.org/10.3168/jds.S0022-0302(94)77212-X): 2695-2703
        • Flower F.C.
        • Weary D.M.
        Effect of hoof pathologies on subjective assessments of dairy cow gait.
        J. Dairy Sci. 2006; 89 (http://dx.doi.org/10.3168/jds.S0022-0302(06)72077-X): 139-146
        • Fregonesi J.A.
        • Leaver J.D.
        Influence of space allowance and milk yield level on behaviour, performance and health of dairy cows housed in strawyard and cubicle systems.
        Livest. Prod. Sci. 2002; 78 (http://dx.doi.org/10.1016/S0301-6226(02)00097-0): 245-257
        • Fregonesi J.A.
        • Tucker C.B.
        • Weary D.M.
        Overstocking reduces lying time in dairy cows.
        J. Dairy Sci. 2007; 90 (http://dx.doi.org/10.3168/jds.2006-794): 3349-3354
        • Gibbons J.
        • Vasseur E.
        • Rushen J.
        • de Passillé A.M.
        A training programme to ensure high repeatability of injury scoring of dairy cows.
        Anim. Welf. 2012; 21 (http://dx.doi.org/10.7120/09627286.21.3.379): 379-388
        • Gomez A.
        • Cook N.B.
        Time budgets of lactating dairy cattle in commercial freestall herds.
        J. Dairy Sci. 2010; 93 (http://dx.doi.org/10.3168/jds.2010-3436): 5772-5781
        • Green L.E.
        • Huxley J.N.
        • Banks C.
        • Green M.J.
        Temporal associations between low body condition, lameness and milk yield in a UK dairy herd.
        Prev. Vet. Med. 2014; 113 (http://dx.doi.org/10.1016/j.prevetmed.2013.10.009): 63-71
        • Halachmi I.
        Simulating the hierarchical order and cow queue length in an automatic milking system.
        Biosyst. Eng. 2009; 102 (http://dx.doi.org/10.1016/j.biosystemseng.2009.01.010): 453-460
        • Helmreich S.
        • Hauser R.
        • Jungbluth T.
        • Wechsler B.
        • Gygax L.
        Time-budget constraints for cows with high milking frequency on farms with automatic milking systems.
        Livest. Sci. 2014; 167 (http://dx.doi.org/10.1016/j.livsci.2014.06.014): 315-322
        • Hillerton J.E.
        • Dearing J.
        • Dale J.
        • Poelarends J.J.
        • Neijenhuis F.
        • Sampimon O.C.
        • Miltenburg J.D.H.M.
        • Fossing C.
        Impact of automatic milking on animal health.
        in: Meijering A. Hogeveen H. Koning C.J.A.M.d. Automatic Milking-A Better Understanding. Wageningen Academic Publishers, Wageningen, the Netherlands2004: 125-134 ()
        • Hohenbrink S.
        • Meinecke-Tillmann S.
        Influence of social dominance on the secondary sex ratio and factors affecting hierarchy in Holstein dairy cows.
        J. Dairy Sci. 2012; 95 (http://dx.doi.org/10.3168/jds.2011-5281): 5694-5701
        • Huxley J.N.
        Impact of lameness and claw lesions in cows on health and production.
        Livest. Sci. 2013; 156 (http://dx.doi.org/10.1016/j.livsci.2013.06.012): 64-70
        • Kester E.
        • Holzhauer M.
        • Frankena K.
        A descriptive review of the prevalence and risk factors of hock lesions in dairy cows.
        Vet. J. 2014; 202 (http://dx.doi.org/10.1016/j.tvjl.2014.07.004): 222-228
        • Klaas I.C.
        • Rousing T.
        • Foussing C.
        • Hindhede J.
        • Sørensen J.T.
        Is lameness a welfare problem in dairy farms with automatic milking systems?.
        Anim. Welf. 2003; 12: 599-603
        • Leonard F.C.
        • O'Connell J.M.
        • O'Farrell K.J.
        Effect of overcrowding on claw health in first-calved friesian heifers.
        Br. Vet. J. 1996; 152 (http://dx.doi.org/10.1016/S0007-1935(96)80040-6): 459-472
        • Lexer D.
        • Hagen K.
        • Palme R.
        • Troxler J.
        • Waiblinger S.
        Time budgets and adrenocortical activity of cows milked in a robot or a milking parlour: interrelationships and influence of social rank.
        Anim. Welf. 2009; 18: 73-80
        • Lim P.Y.
        • Huxley J.N.
        • Willshire J.A.
        • Green M.J.
        • Othman A.R.
        • Kaler J.
        Unravelling the temporal association between lameness and body condition score in dairy cattle using a multistate modelling approach.
        Prev. Vet. Med. 2015; 118 (http://dx.doi.org/10.1016/j.prevetmed.2014.12.015): 370-377
        • Melin M.
        • Hermans G.G.N.
        • Pettersson G.
        • Wiktorsson H.
        Cow traffic in relation to social rank and motivation of cows in an automatic milking system with control gates and an open waiting area.
        Appl. Anim. Behav. Sci. 2006; 96 (http://dx.doi.org/10.1016/j.applanim.2005.06.013): 201-214
        • Menke C.
        • Waiblinger S.
        • Fölsch D.W.
        • Wiepkema P.R.
        Social behaviour and injuries of horned cows in loose housing systems.
        Anim. Welf. 1999; 8: 243-258
        • Miguel-Pacheco G.G.
        • Kaler J.
        • Remnant J.
        • Cheyne L.
        • Abbott C.
        • French A.P.
        • Pridmore T.P.
        • Huxley J.N.
        Behavioural changes in dairy cows with lameness in an automatic milking system.
        Appl. Anim. Behav. Sci. 2014; 150 (http://dx.doi.org/10.1016/j.applanim.2013.11.003): 1-8
        • Miller K.
        • Wood-Gush D.G.M.
        Some effects of housing on the social behaviour of dairy cows.
        Anim. Sci. 1991; 53 (http://dx.doi.org/10.1017/S0003356100020262): 271-278
        • Nordlund K.
        • Cook N.B.
        A flowchart for evaluating dairy cow freestalls.
        Bov. Pract. 2003; 37: 89-96
        • Pérez-Cabal M.A.
        • Alenda R.
        Clinical lameness and risk factors in a Spanish Holstein population.
        Livest. Sci. 2014; 164 (http://dx.doi.org/10.1016/j.livsci.2014.03.012): 168-174
        • Randall L.V.
        • Green M.J.
        • Chagunda M.G.G.
        • Mason C.
        • Archer S.C.
        • Green L.E.
        • Huxley J.N.
        Low body condition predisposes cattle to lameness: An 8-year study of one dairy herd.
        J. Dairy Sci. 2015; 98 (http://dx.doi.org/10.3168/jds.2014-8863): 3766-3777
        • Sarjokari K.
        • Kaustell K.O.
        • Hurme T.
        • Kivinen T.
        • Peltoniemi O.A.T.
        • Saloniemi H.
        • Rajala-Schultz P.J.
        Prevalence and risk factors for lameness in insulated free stall barns in Finland.
        Livest. Sci. 2013; 156 (http://dx.doi.org/10.1016/j.livsci.2013.06.010): 44-52
        • Solano L.
        • Barkema H.W.
        • Pajor E.A.
        • Mason S.
        • LeBlanc S.J.
        • Zaffino Heyerhoff J.C.
        • Nash C.G.R.
        • Haley D.B.
        • Vasseur E.
        • Pellerin D.
        • Rushen J.
        • de Passillé A.M.
        • Orsel K.
        Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
        J. Dairy Sci. 2015; 98 (http://dx.doi.org/10.3168/jds.2015-9652): 6978-6991
        • Tucker C.B.
        • Weary D.M.
        • Fraser D.
        Free-stall dimensions: Effects on preference and stall usage.
        J. Dairy Sci. 2004; 87 (http://dx.doi.org/10.3168/jds.S0022-0302(04)73271-3): 1208-1216
      3. Valacta. 2015. L'évolution de la production laitière Québécoise 2014. Quebec. Accessed Aug. 21, 2015. www.valacta.com/FR/Nos-publications/Documents/%C3%89VOLUTION%20LAITI%C3%88RE/%C3%89volution%202014_FINAL.pdf

        • van Gastelen S.
        • Westerlaan B.
        • Houwers D.J.
        • van Eerdenburg F.J.C.M.
        A study on cow comfort and risk for lameness and mastitis in relation to different types of bedding materials.
        J. Dairy Sci. 2011; 94 (http://dx.doi.org/10.3168/jds.2010-4019): 4878-4888
        • Vasseur E.
        • Gibbons J.
        • Rushen J.
        • de Passillé A.M.
        Development and implementation of a training program to ensure high repeatability of body condition scoring of dairy cows.
        J. Dairy Sci. 2013; 96 (http://dx.doi.org/10.3168/jds.2012-6359): 4725-4737
        • Vasseur E.
        • Gibbons J.
        • Rushen J.
        • Pellerin D.
        • Pajor E.
        • Lefebvre D.
        • de Passillé A.M.
        An assessment tool to help producers improve cow comfort on their farms.
        J. Dairy Sci. 2015; 98 (http://dx.doi.org/10.3168/jds.2014-8224): 698-708
        • Vosika B.
        • Lexer D.
        • Stanek C.
        • Troxler J.
        • Waiblinger S.
        The influence of an automatic milking system on claw health and lameness of dairy cows.
        in: Meijering A. Hogeveen H. Koning C.J.A.M.d. Automatic Milking—A Better Understanding. Wageningen Academic Publishers, Wageningen, the Netherlands2004: 155-160
        • Westin R.
        • Vaughan A.
        • Passillé A.M.d.
        • DeVries T.J.
        • Pajor E.A.
        • Pellerin D.
        • Siegford J.M.
        • Vasseur E.
        • Rushen J.
        Lying times of lactating cows on dairy farms with automatic milking systems and the relation to lameness, leg lesions and body condition score.
        J. Dairy Sci. 2016; 99 (http://dx.doi.org/10.3168/jds.2015-9737): 551-561
        • Wildman E.E.
        • Jones G.M.
        • Wagner P.E.
        • Boman R.L.
        • Troutt Jr., H.F.
        • Lesch T.N.
        A dairy cow body condition scoring system and its relationship to selected production characteristics.
        J. Dairy Sci. 1982; 65 (http://dx.doi.org/10.3168/jds.S0022-0302(82)82223-6): 495-501
        • Zaffino Heyerhoff J.C.
        • LeBlanc S.J.
        • DeVries T.J.
        • Nash C.G.R.
        • Gibbons J.
        • Orsel K.
        • Barkema H.W.
        • Solano L.
        • Rushen J.
        • de Passillé A.M.
        • Haley D.B.
        Prevalence of and factors associated with hock, knee, and neck injuries on dairy cows in freestall housing in Canada.
        J. Dairy Sci. 2014; 97 (http://dx.doi.org/10.3168/jds.2012-6367): 173-184