Advertisement
Research| Volume 100, ISSUE 8, P6516-6526, August 2017

Effects of changing freestall area on lameness, lying time, and leg injuries on dairy farms in Alberta, Canada

Open ArchivePublished:June 07, 2017DOI:https://doi.org/10.3168/jds.2016-12467

      ABSTRACT

      Cow comfort is of increasing importance in the dairy industry, due to an increased focus on animal welfare. However, whether producer changes to the cows' environment affect cow comfort has not been well characterized. Our objectives were to: (1) quantify the effect of freestall area changes on the prevalence of lameness, leg injuries, and average lying time; and (2) compare cow comfort outcomes on farms that had never had an assessment of cow comfort to farms that had had a previous assessment of cow comfort. A sample of 60 Holstein-Friesian cows were selected on each of 15 farms that made changes to the freestall area after an assessment of cow comfort (change, CHG); 15 farms that did not make changes to the freestall area after an assessment of cow comfort (no change, NC); and 14 farms that had yet to be evaluated (new farms, NF). Cows in NC and NF were lame 1.50 and 1.71 times more often, respectively, than cows on CHG farms. Additionally, daily lying time was 0.33 and 0.62 h/d lower in NC and NF, respectively, than on CHG farms. The prevalence of hock and knee injuries was not different among the 3 groups of farms. No differences were detected in the parameters of interest when comparing NF with NC farms; therefore, we concluded that the NC group was not biased by a previous assessment of cow comfort. Farms in the CHG group had a lower prevalence of lame cows and greater lying time than the NC and NF groups.

      Key words

      INTRODUCTION

      Worldwide, cow comfort and animal welfare are important topics for the dairy industry (
      • Barkema H.W.
      • von Keyserlingk M.A.G.
      • Kastelic J.P.
      • Lam T.J.G.M.
      • Luby C.
      • Roy J.-P.
      • LeBlanc S.J.
      • Keefe G.P.
      • Kelton D.F.
      Invited review: Changes in the dairy industry affecting dairy health and welfare.
      ). Quantifying animal-based measures (e.g., prevalence of lameness and injuries, lying time, and production information), evaluating environmental factors (e.g., barn design and stall dimensions), and determining management practices (e.g., record keeping and management training) are proven methods of assessing animal welfare (
      • National Farm Animal Care Council
      Code of practice for the care and handling of dairy cattle.
      ;
      • 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.
      ;
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). Furthermore, prevalence and risk factors for animal-based measures of cow comfort—including prevalence of lameness (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ), prevalence of injuries (
      • 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 lying behavior (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Nash C.G.R.
      • Orsel K.
      Associations between lying behavior and lameness in Canadian Holstein-Friesian cows housed in freestall barns.
      )—on freestall dairy farms have been well characterized, providing baseline information.
      Lameness reduces cow comfort and welfare, because it causes pain (
      • Rushen J.
      • Haley D.
      • de Passille A.M.
      Effect of softer flooring in tie stalls on resting behavior and leg injuries of lactating cows.
      ), negatively affects lying time (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Nash C.G.R.
      • Orsel K.
      Associations between lying behavior and lameness in Canadian Holstein-Friesian cows housed in freestall barns.
      ), and is associated with low BCS (
      • 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.
      ). Lameness prevalence varies by region, housing system, and management practice (
      • Cook N.B.
      • Nordlund K.V.
      The influence of the environment on dairy cow behavior, claw health and herd lameness dynamics.
      ), ranging from 0 to 55% in North America (
      • Cook N.B.
      • Bennett T.B.
      • Nordlund K.V.
      Effect of freestall surface on daily activity patterns in dairy cows with relevance to lameness prevalence.
      ;
      • von Keyserlingk M.A.G.
      • Barrientos A.
      • Ito K.
      • Galo E.
      • Weary D.M.
      Benchmarking cow comfort on North American freestall dairies: Lameness, leg injuries, lying time, facility design, and management for high-producing Holstein dairy cows.
      ;
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). Hock and knee injuries are important animal welfare concerns, because they can cause pain and discomfort (
      • Huxley J.N.
      • Whay H.R.
      Current attitudes of cattle practitioners to pain and the use of analgesics in cattle.
      ). In a recent Canadian study, 47 and 24% of cows had hock and knee injuries, respectively (
      • 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.
      ). Lying time is considered a good indicator of cow comfort, because it is optimal for cows to lie down 12 to 13 h/d (
      • Weary D.M.
      • Huzzey J.M.
      • von Keyserlingk M.A.G.
      Using behavior to predict and identify ill health in animals.
      ), and suboptimal stall design reduces lying time (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Nash C.G.R.
      • Orsel K.
      Associations between lying behavior and lameness in Canadian Holstein-Friesian cows housed in freestall barns.
      ). Lying time can also be linked with animal welfare: cows with restricted lying time have a higher risk of aggressive behavior, decreased productivity, poor hoof health, and compromised overall health (
      • Ito K.
      ;
      • Nechanitzky K.
      • Starke A.
      • Vidondo B.
      • Müller H.
      • Reckardt M.
      • Friedli K.
      • Steiner A.
      Analysis of behavioral changes in dairy cows associated with claw horn lesions.
      ;
      • Wang F.X.
      • Shao D.F.
      • Li S.L.
      • Wang Y.J.
      • Azarfar A.
      • Cao Z.J.
      Effects of stocking density on behavior, productivity, and comfort indices of lactating dairy cows.
      ).
      The Dairy Farmers of Canada and the National Farm Animal Care Council created the Canadian Dairy Code of Practice, containing recommended practices and requirements for Canadian dairy producers (
      • 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.
      ). To supplement this code, a cow comfort assessment was developed to help producers assess their compliance with the Canadian Dairy Code of Practice (
      • 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.
      ). These assessments were implemented in several Canada-wide studies (
      • 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.
      ;
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ), allowing participating producers to determine how their farms compared with regulations in the Canadian Dairy Code of Practice.
      Although information about cow comfort and welfare is available for Canadian dairy producers (e.g., the Canadian Dairy Code of Practice and conferences), it is not known whether they use this information to make changes that would improve cow comfort and welfare on their farm. Additionally, effects on measures of cow comfort after implementation of changes have not yet been reported. The objectives of this study were to: (1) quantify the effect of freestall area changes on lameness, leg injuries, and average lying time; and (2) compare cow comfort outcomes on farms that had never had an assessment of cow comfort to outcomes on farms that had had a previous assessment of cow comfort.

      MATERIALS AND METHODS

      The method of assessing cow comfort has been described in
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      and
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Nash C.G.R.
      • Orsel K.
      Associations between lying behavior and lameness in Canadian Holstein-Friesian cows housed in freestall barns.
      . Methods included evaluations of animal-based measures, environmental measures, and management factors (
      • 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.
      ). Standard operating procedures used were consistent with previous studies and are reported on the Canadian Dairy Research Portal (http://www.dairyresearch.ca/animal-comfort-tool.php).

      Farm Selection and Farm Visits

      The 91 freestall dairy farms in Alberta, Canada, that participated in a previous evaluation of comfort and lameness (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ) were invited to participate in the present study, with a target of 15 farms per study group. Thirty of these previously assessed farms were selected for participation because they met all the criteria. Remaining farms were not included because of logistic and scheduling conflicts. As well, 45 Alberta freestall dairy farms that did not participate in the
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      study were also invited to participate as a contemporary control group. Fifteen of those 45 farms were selected for the present study (Figure 1), because some of the farms that met the inclusion criteria did not fit into the farm visit schedule. Inclusion criteria are presented in Table 1 and were consistent with the previous study, as described by
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      .
      Figure thumbnail gr1
      Figure 1Selection process for participating farms in their respective study groups.
      Table 1Inclusion criteria for the 3 groups of farms
      Study groupInclusion criteria
      New farms (NF)Enrolled in DHI; ≥60 lactating Holstein-Friesian cows; freestall facility has been at least 1 yr in use; pasture access ≤2 h/d; no automatic milking systems; not enrolled in 2011 study
      Solano et al. (2015).
      Freestall change farms (CHG)Same criteria as for NF herds, plus changes made to one or more of the following after 2011 assessment
      Solano et al. (2015).
      : stall dimensions, stall base type, stall bedding type, quantity of stall bedding, number of stalls, barn flooring type
      No freestall change farms (NC)Same criteria as for NF herds, plus no changes made to the following after 2011 assessment
      Solano et al. (2015).
      : stall dimensions, stall base type, stall bedding type, quantity of stall bedding, number of stalls, barn flooring type
      1
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      .
      A total of 44 freestall dairy farms in Alberta were visited between April and December of 2015. Participating farms included 15 that made changes to their freestall area (change, CHG); 15 that did not make changes, or made changes that were not related to the freestall area (no change, NC); and 14 that were not enrolled in the
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      study and had had no previous cow comfort assessment (new farms, NF). Initially, the NF group consisted of 15 farms, but 1 was excluded from the analyses because of a disruption in data collection (Figure 1). Each farm was visited twice, with an interval of 5 to 7 d between visits. During the first visit, both environmental and animal-based measures were collected. During the second visit, a face-to-face questionnaire was conducted, along with any data collection (animal-based and environmental) that had not been completed during the first visit. Data were collected by trained research personnel from the University of Calgary (AB, Canada), including 1 graduate student (EM) and 2 research assistants. All protocols and questionnaires were approved by the University of Calgary Animal Care Committee (AC14–0216) and Research Ethics Board (REB14–2120).

      Cow-Based Measures

      Forty cows between 10 and 120 DIM were selected on each farm, because this period is known to be a critical one for lameness (
      • Green L.E.
      • Hedges V.J.
      • Schukken Y.H.
      • Blowey R.W.
      • Packington A.J.
      The impact of clinical lameness on the milk yield of dairy cattle.
      ). If the farm did not have 40 lactating cows between 10 and 120 DIM, cows >120 DIM were added, reflective of the parity distribution on farm. As well, another 20 cows between 121 and 310 DIM were selected to ensure that all stages of lactation were represented. Pen was not considered in the selection process, but pens were excluded if they did not contain freestalls. This occurred on 2 farms, and cows from these pens were excluded from the study. Additional cows from alternative pens with freestalls were selected according to the selection criteria.
      Selected cows were assessed for lameness, knee injuries, hock injuries, lying time, claw length, leg cleanliness, and BCS using standard operating procedures as described on the Canadian Dairy Research Portal (https://www.dairyresearch.ca/animal-comfort-tool.php). In short, cows were video-recorded for lameness as they exited the milking parlor. Recordings were subsequently viewed and cows were scored for lameness based on the presence of a head bob, asymmetrical steps, and a limp. Cows were considered lame if they had a limp (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ;
      • 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.
      ). Knee and hock injuries, claw length, leg cleanliness, and BCS were assessed in the milking parlor or in the freestall area (if data could not reliably be collected in the milking parlor). Hock injuries were scored on a scale of 0 to 3, separately for the left and right joints (
      • Gibbons J.
      • Vasseur E.
      • Rushen J.
      • de Passillé A.M.
      A training programme to ensure high repeatability of injury scoring of dairy cows.
      ). Lying time was measured using electronic data loggers (HOBO Pendant G Acceleration Data Loggers; Onset Computer Corp., Pocasset, MA). The loggers were placed on the cows' hind leg to measure lying time for 4 full days. Claw length was estimated by assessing the angle at which each of the hind claws contacted the ground. Claws were scored as overgrown if the angle of their hoof was <45°, or not overgrown if the angle of the hoof was ≥45° (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). For leg cleanliness, the lower hind right leg was scored on a scale of 0 (clean) to 3 (caked manure on >50% of the lower leg), a method adapted from . The BCS was scored in 0.25 increments (
      • 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.
      ). Data including DIM, parity, and milk production were provided by CanWest DHI (Guelph, ON, Canada).

      Environmental Measures

      Environmental measures in the freestall area have been described by
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      , and included measurements of pens, feed alleys, feeding area, watering area, stall dimensions, number of usable stalls, number of cows within the pen, cleanliness of the pens, and slipperiness of the pen floor. Additionally, the stall was assessed for type of base, type of bedding, depth of bedding, and cleanliness (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ).

      Management

      A questionnaire was administered during farm visits. This questionnaire captured cow comfort related to management practices on the farm (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). In short, the management questionnaire captured information on how stalls were maintained (e.g., freestall bedding frequency, freestall cleaning frequency, freestall floor cleaning frequency). The questionnaire is available on the Canadian Dairy Research Portal (https://www.dairyresearch.ca/animal-comfort-tool.php). Additionally, in 2015, questions were added to capture changes that had been made on the farms, including whether farmers had made changes after the assessment in 2011, and what kind of changes they made. The questions asked in the initial phone interview about changes aided in placing farms in their respective study groups.

      Statistical Analyses

      Data were entered into Microsoft Access (Microsoft Corp., Redmond, WA), and statistical analyses were performed using Stata 13.1 (2013; StataCorp, College Station, TX). For all analyses, a P-value <0.05 was considered significant. Descriptive statistics (mean, median, interquartile range) were used to describe both herd- and cow-level characteristics. Variables were tested for collinearity, and continuous variables were tested for normal distribution. If distributions were abnormal, data were log-transformed. If the log-transformation did not meet the preferred distribution, variables were categorized.
      Cow-level descriptive statistics for hock and knee injuries were completed by dichotomizing the variable to not injured (scores 0 and 1) and injured (scores 2 and 3). In agreement with the study by
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      , analyses of injury scores were completed using the higher score of the 2 limbs. Claw data were also dichotomous, analyzing the higher of the 2 scores. Two farms in the CHG group, 1 farm in the NC group, and 2 farms in the NF group had >1 pens with different environmental characteristics. On these farms, the pen with the largest number of cows in the study was used for analysis. Multiple different characteristics within each pen were classified as “other.” Additionally, if environmental variables were uncommon, such as recycled rubber for stall base type (which existed on only 1 farm), they were also classified as “other.” If predictor variables were correlated at r ≥0.70 (i.e., stall base and bedding type), only the variable with the strongest association or the higher number of observations was included in the analysis (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). These descriptive data were analyzed using Student's t-tests for continuous variables and Pearson's chi-squared analysis for categorical variables.
      A total of 4 final models were constructed for each of the outcome measures selected for cow comfort and welfare: lameness, lying time, hock injuries, and knee injuries. Logistic regression models were used for lameness and hock and knee injuries; a linear regression model was constructed for lying time. To construct each model, univariable analyses were performed for each of these outcomes to identify which explanatory variables had an association (P ≤ 0.25) with the outcome variable. Next, a multilevel regression model was constructed for each category of explanatory variables: cow-level variables, environmental variables, and management variables (GLMM in Stata) resulting in 3 models per outcome variable. Variables from each of the 3 models associated with the outcome at P < 0.10 were included in the final model. Backward elimination was used to build the final model, keeping any confounding variables (when the removal of a variable changed the coefficient of another variable by 30%) and any variables significant at P ≤ 0.05. Interactions were tested between remaining variables that made biological sense or when they were part of the lameness causal web (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). In all models, herd was a random effect and study groups were added as fixed effects.

      RESULTS

      Of the farms contacted, the 15 CHG farms (Figure 1) are described in Table 2. The most frequent changes were changing stall base to geomattress (n = 7), increasing bedding quantity (n = 7), and grooving crossover alleys (n = 6). The changes included in the study were reported by farmers, although 2 farms ended up switching groups due to discrepancies in what the farmer originally stated and what was measured by the researchers. The 44 farms included had an average herd size of 158 lactating cows. We observed no difference in herd size, average parity, DIM, or daily milk yield between the CHG, NC, and NF groups (Table 3). Claws were less often overgrown in CHG than in NC and NF. Cows in CHG herds had cleaner legs than cows in NC and NF herds (Table 3). The CHG herds had fewer cows with a BCS ≤2.5 than NC and NF herds, and more cows with a BCS of 2.75–3.25 than NF herds. The NF herds had more cows with a BCS of ≤3.5 than NC herds (Table 3).
      Table 2Summary of changes to the freestall area after initial assessment of cow comfort (15 change herds, CHG)
      FarmFreestall area changes
      1Replaced rubber mat stall base with geomattress, changed bedding type from straw to shavings, increased bedding quantity and frequency, grooved crossover alleys
      2Removed previous stall base, changed bedding type from shavings to composted manure, increased bedding quantity
      3Replaced rubber mat stall base with geomattress, decreased stall length, increased bedding frequency, grooved freestall area flooring, grooved crossover alleys
      4Replaced concrete stall base with rubber mat, changed bedding type from deep bedded straw to shavings, decreased bedding quantity
      5Replaced concrete stall base with geomattress, increased bedding quantity
      6Increased number of stalls, increased width of stalls, replaced old rubber mat stall base with new rubber mats, increased bedding frequency
      7Changed bedding type from straw to shavings, grooved crossover alleys
      8Increased bedding quantity, increased scraping frequency in the freestall area
      9Decreased number of stalls, increased stall width, decreased stall length, replaced rubber mat stall base with geomattress
      10Replaced waterbed stall base with geomattress, increased bedding quantity
      11Replaced concrete stall base with geomattress, changed bedding type from straw to shavings
      12Increased bedding quantity, grooved crossover alleys
      13Replaced rubber mat stall base with geomattress
      14Increased bedding quantity and frequency, grooved crossover alleys
      15Replaced stall base with geomattress, changed bedding type from straw to shavings, grooved crossover alleys
      Table 3Distribution of continuous (mean ± SD) and categorical [no. (%)] cow-level variables on farms that made changes to the freestall area (CHG), made no changes to the freestall area (NC), and new farms (NF)
      Herd or study group characteristicCHG (n = 15)NC (n = 15)NF (n = 14)
      Herd size (no. of lactating cows)
      Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
      160 ± 75
      Within a row, percentages without a common superscript differed (P < 0.05).
      162 ± 75
      Within a row, percentages without a common superscript differed (P < 0.05).
      151 ± 61
      Within a row, percentages without a common superscript differed (P < 0.05).
      Herd parity
      Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
      2.5 ± 1.4
      Within a row, percentages without a common superscript differed (P < 0.05).
      2.4 ± 1.1
      Within a row, percentages without a common superscript differed (P < 0.05).
      2.3 ± 1.1
      Within a row, percentages without a common superscript differed (P < 0.05).
      Herd DIM
      Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
      114 ± 34
      Within a row, percentages without a common superscript differed (P < 0.05).
      110 ± 38
      Within a row, percentages without a common superscript differed (P < 0.05).
      111 ± 27
      Within a row, percentages without a common superscript differed (P < 0.05).
      Parity
      Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
      2.7 ± 1.6
      Within a row, percentages without a common superscript differed (P < 0.05).
      2.3 ± 1.4
      Within a row, percentages without a common superscript differed (P < 0.05).
      2.2 ± 1.3
      Within a row, percentages without a common superscript differed (P < 0.05).
      DIM
      Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
      Median and interquartile range.
      95 ± 99
      Within a row, percentages without a common superscript differed (P < 0.05).
      93 ± 98
      Within a row, percentages without a common superscript differed (P < 0.05).
      94 ± 95
      Within a row, percentages without a common superscript differed (P < 0.05).
      Daily milk yield (kg)
      Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
      37 ± 10
      Within a row, percentages without a common superscript differed (P < 0.05).
      33 ± 7
      Within a row, percentages without a common superscript differed (P < 0.05).
      34 ± 7
      Within a row, percentages without a common superscript differed (P < 0.05).
      DIM
      Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
       10–60221 (26)
      Within a row, percentages without a common superscript differed (P < 0.05).
      207 (25)
      Within a row, percentages without a common superscript differed (P < 0.05).
      186 (26)
      Within a row, percentages without a common superscript differed (P < 0.05).
       61–120268 (32)
      Within a row, percentages without a common superscript differed (P < 0.05).
      231 (28)
      Within a row, percentages without a common superscript differed (P < 0.05).
      220 (30)
      Within a row, percentages without a common superscript differed (P < 0.05).
       121–180142 (17)
      Within a row, percentages without a common superscript differed (P < 0.05).
      161 (20)
      Within a row, percentages without a common superscript differed (P < 0.05).
      114 (16)
      Within a row, percentages without a common superscript differed (P < 0.05).
       181–240127 (15)
      Within a row, percentages without a common superscript differed (P < 0.05).
      117 (14)
      Within a row, percentages without a common superscript differed (P < 0.05).
      112 (15)
      Within a row, percentages without a common superscript differed (P < 0.05).
       241–30081 (10)
      Within a row, percentages without a common superscript differed (P < 0.05).
      102 (12)
      Within a row, percentages without a common superscript differed (P < 0.05).
      96 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Claw length
       Not overgrown714 (85)
      Within a row, percentages without a common superscript differed (P < 0.05).
      640 (78)
      Within a row, percentages without a common superscript differed (P < 0.05).
      544 (75)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Overgrown125 (15)
      Within a row, percentages without a common superscript differed (P < 0.05).
      178 (22)
      Within a row, percentages without a common superscript differed (P < 0.05).
      184 (25)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Leg cleanliness
       Clean747 (89)
      Within a row, percentages without a common superscript differed (P < 0.05).
      677 (83)
      Within a row, percentages without a common superscript differed (P < 0.05).
      566 (78)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Dirty92 (12)
      Within a row, percentages without a common superscript differed (P < 0.05).
      141 (17)
      Within a row, percentages without a common superscript differed (P < 0.05).
      162 (22)
      Within a row, percentages without a common superscript differed (P < 0.05).
      BCS
       ≤2.5260 (31)
      Within a row, percentages without a common superscript differed (P < 0.05).
      295 (36)
      Within a row, percentages without a common superscript differed (P < 0.05).
      270 (37)
      Within a row, percentages without a common superscript differed (P < 0.05).
       2.75–3.25487 (58)
      Within a row, percentages without a common superscript differed (P < 0.05).
      442 (54)
      Within a row, percentages without a common superscript differed (P < 0.05).
      364 (50)
      Within a row, percentages without a common superscript differed (P < 0.05).
       ≥3.592 (11)
      Within a row, percentages without a common superscript differed (P < 0.05).
      81 (9)
      Within a row, percentages without a common superscript differed (P < 0.05).
      94 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      a,b Within a row, percentages without a common superscript differed (P < 0.05).
      1 Data from 2015. Source: CanWest DHI (Guelph, ON, Canada).
      2 Median and interquartile range.

      Environmental and Management Variables

      Study groups did not differ in terms of stall dimensions, stall bedding type, or bedding cleanliness (Table 4). The CHG herds were less likely to have rubber mats and more likely to have geomattresses as a stall base than NC herds, and they also had stalls with bedding ≥2 cm more often than NC and NF herds. The NF herds had wetter stalls than CHG and NC herds. The CHG farms bedded stalls more frequently and had less slippery flooring than did NC and NF (Table 4).
      Table 4Distribution of continuous (mean ± SD) and categorical [no. (%)] herd-level environmental and management variables for freestall change farms (CHG), no freestall change farms (NC), and new farms (NF)
      Previous farms not included in the total, because they were not included in the final models.
      VariableCHG (n = 15)NC (n = 15)NF (n = 14)
      Average stocking density1.02 ± 0.11
      Within a row, percentages without a common superscript differed (P < 0.05).
      1.11 ± 0.20
      Within a row, percentages without a common superscript differed (P < 0.05).
      1.25 ± 0.17
      Within a row, percentages without a common superscript differed (P < 0.05).
      Width (cm)119 ± 7
      Within a row, percentages without a common superscript differed (P < 0.05).
      116 ± 7
      Within a row, percentages without a common superscript differed (P < 0.05).
      114 ± 5
      Within a row, percentages without a common superscript differed (P < 0.05).
      Bed length (cm)180 ± 10
      Within a row, percentages without a common superscript differed (P < 0.05).
      178 ± 11
      Within a row, percentages without a common superscript differed (P < 0.05).
      178 ± 16
      Within a row, percentages without a common superscript differed (P < 0.05).
      Brisket board (cm)11 ± 5
      Within a row, percentages without a common superscript differed (P < 0.05).
      12 ± 4
      Within a row, percentages without a common superscript differed (P < 0.05).
      12 ± 7
      Within a row, percentages without a common superscript differed (P < 0.05).
      Neck rail height (cm)117 ± 9
      Within a row, percentages without a common superscript differed (P < 0.05).
      117 ± 8
      Within a row, percentages without a common superscript differed (P < 0.05).
      115 ± 7
      Within a row, percentages without a common superscript differed (P < 0.05).
      Neck rail to rear curb (cm)172 ± 8
      Within a row, percentages without a common superscript differed (P < 0.05).
      170 ± 9
      Within a row, percentages without a common superscript differed (P < 0.05).
      167 ± 10
      Within a row, percentages without a common superscript differed (P < 0.05).
      Curb height (cm)20 ± 6
      Within a row, percentages without a common superscript differed (P < 0.05).
      22 ± 4
      Within a row, percentages without a common superscript differed (P < 0.05).
      22 ± 6
      Within a row, percentages without a common superscript differed (P < 0.05).
      Obstruction in lunge space
       Yes9 (60)
      Within a row, percentages without a common superscript differed (P < 0.05).
      11 (73)
      Within a row, percentages without a common superscript differed (P < 0.05).
      10 (71)
      Within a row, percentages without a common superscript differed (P < 0.05).
       No6 (40)
      Within a row, percentages without a common superscript differed (P < 0.05).
      4 (27)
      Within a row, percentages without a common superscript differed (P < 0.05).
      4 (29)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Stall base type
       Concrete1 (7)
      Within a row, percentages without a common superscript differed (P < 0.05).
      4 (27)
      Within a row, percentages without a common superscript differed (P < 0.05).
      4 (29)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Rubber mat1 (7)
      Within a row, percentages without a common superscript differed (P < 0.05).
      2 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      5 (36)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Geomattress9 (60)
      Within a row, percentages without a common superscript differed (P < 0.05).
      6 (40)
      Within a row, percentages without a common superscript differed (P < 0.05).
      3 (21)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Waterbed1 (7)
      Within a row, percentages without a common superscript differed (P < 0.05).
      1 (7)
      Within a row, percentages without a common superscript differed (P < 0.05).
      0 (0)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Other3 (20)
      Within a row, percentages without a common superscript differed (P < 0.05).
      2 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      2 (14)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Stall bedding type
       Straw2 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      2 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      3 (21)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Sawdust1 (7)
      Within a row, percentages without a common superscript differed (P < 0.05).
      2 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      3 (21)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Wood shavings9 (60)
      Within a row, percentages without a common superscript differed (P < 0.05).
      7 (47)
      Within a row, percentages without a common superscript differed (P < 0.05).
      5 (36)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Sand1 (7)
      Within a row, percentages without a common superscript differed (P < 0.05).
      0 (0)
      Within a row, percentages without a common superscript differed (P < 0.05).
      0 (0)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Other2 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      4 (27)
      Within a row, percentages without a common superscript differed (P < 0.05).
      3 (21)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Bedding depth (cm)
       >210 (67)
      Within a row, percentages without a common superscript differed (P < 0.05).
      6 (41)
      Within a row, percentages without a common superscript differed (P < 0.05).
      6 (43)
      Within a row, percentages without a common superscript differed (P < 0.05).
       ≤25 (33)
      Within a row, percentages without a common superscript differed (P < 0.05).
      9 (60)
      Within a row, percentages without a common superscript differed (P < 0.05).
      8 (57)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Bedding cleanliness
       Clean14 (93)
      Within a row, percentages without a common superscript differed (P < 0.05).
      12 (80)
      Within a row, percentages without a common superscript differed (P < 0.05).
      11 (79)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Dirty1 (7)
      Within a row, percentages without a common superscript differed (P < 0.05).
      3 (20)
      Within a row, percentages without a common superscript differed (P < 0.05).
      3 (21)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Bedding dryness
       Dry12 (80)
      Within a row, percentages without a common superscript differed (P < 0.05).
      13 (87)
      Within a row, percentages without a common superscript differed (P < 0.05).
      10 (71)
      Within a row, percentages without a common superscript differed (P < 0.05).
       Wet3 (20)
      Within a row, percentages without a common superscript differed (P < 0.05).
      2 (13)
      Within a row, percentages without a common superscript differed (P < 0.05).
      4 (29)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Bedding frequency
       >Once/wk11 (73)
      Within a row, percentages without a common superscript differed (P < 0.05).
      8 (53)
      Within a row, percentages without a common superscript differed (P < 0.05).
      6 (43)
      Within a row, percentages without a common superscript differed (P < 0.05).
       ≤Once/wk4 (27)
      Within a row, percentages without a common superscript differed (P < 0.05).
      7 (47)
      Within a row, percentages without a common superscript differed (P < 0.05).
      8 (57)
      Within a row, percentages without a common superscript differed (P < 0.05).
      Slipperiness
       Less slippery (≤3% slip)14 (90)
      Within a row, percentages without a common superscript differed (P < 0.05).
      12 (80)
      Within a row, percentages without a common superscript differed (P < 0.05).
      10 (71)
      Within a row, percentages without a common superscript differed (P < 0.05).
       More slippery (>3% slip)1 (10)
      Within a row, percentages without a common superscript differed (P < 0.05).
      3 (20)
      Within a row, percentages without a common superscript differed (P < 0.05).
      4 (29)
      Within a row, percentages without a common superscript differed (P < 0.05).
      a–c Within a row, percentages without a common superscript differed (P < 0.05).
      1 Previous farms not included in the total, because they were not included in the final models.

      Outcome Variables

      The CHG herds had a lower prevalence of lameness than NC herds and NF herds (Table 5). The prevalence of knee injuries was lower on CHG herds than on NF herds, and CHG herds had a higher average lying time than NC and NF herds (Table 5).
      Table 5Cow-level prevalence [no. (%)] of lameness, hock, and knee injuries and average lying time for freestall change farms (CHG), no freestall change farms (NC), and new farms (NF)
      VariableCHG (n = 15)NC (n = 15)NF (n = 14)
      Lameness
       Not lame747 (89)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      660 (81)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      560 (78)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
       Lame92 (11)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      158 (19)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      160 (22)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      Hock injuries
       Not injured (score 0–1)788 (94)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      775 (95)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      677 (93)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
       Injured (score 2–3)51 (6)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      43 (5)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      51 (7)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      Knee injuries
       Not injured (score 0–1)814 (97)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      769 (94)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      662 (91)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
       Injured (score 2–3)25 (3)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      49 (6)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      66 (9)
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      Average hours lying per day (mean ± SD)12.6 ± 2.3
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      10.7 ± 2.2
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      10.2 ± 3.1
      Within a row, percentages or estimates without a common superscript differed (P < 0.05).
      a,b Within a row, percentages or estimates without a common superscript differed (P < 0.05).

      Multivariable Analyses

      Lameness

      The final model for lameness included study group, parity, DIM, BCS, hock injury, claw overgrowth, and bedding quantity (Table 6). None of the tested interactions were significant at P < 0.05. Cows in NC and NF herds had higher odds of being lame than cows in CHG herds. Cows in parity 2, 3, 4, and higher had higher odds of lameness than first-parity heifers. Cows at 60–120 and 121–180 DIM had lower odds of being lame than cows in their first 60 DIM or cows later in lactation (181–300 DIM). Cows with a low BCS, injured hocks, and overgrown claws had higher odds of lameness. (Table 6).
      Table 6Final logistic regression models for lameness and hock and knee injuries, including cow and herd-level variables from 44 dairy farms in Alberta, Canada
      VariableLamenessHock injuriesKnee injuries
      Odds ratio (95% CI)P-valueOdds ratio (95% CI)P-valueOdds ratio (95% CI)P-value
      Study group
      CHG = farms that made changes to the freestall area after a previous cow comfort assessment; NC = farms that did not make changes to the freestall area after a previous cow comfort assessment; NF = farms that had not been exposed to the previous cow comfort assessment.
       CHGReferentReferentReferent
       NC1.50 (1.17–2.12)0.031.14 (0.92–1.75)0.531.04 (0.81–1.53)0.77
       NF1.71 (1.23–2.35)0.021.25 (0.99–2.01)0.561.17 (0.89–1.77)0.37
      Parity
       1ReferentReferentReferent
       21.51 (1.11–1.89)<0.0011.00 (0.77–1.52)0.450.98 (0.64–1.64)0.52
       32.20 (1.96–2.33)<0.0010.94 (0.51–1.67)0.691.26 (0.96–2.09)0.13
       4+3.10 (2.85–3.52)<0.0011.16 (0.92–1.90)0.101.67 (1.42–2.01)0.02
      DIM
       1–60ReferentReferent
       61–1200.40 (0.19–0.82)<0.0010.92 (0.84–1.58)0.36
       121–1800.64 (0.27–0.94)0.041.28 (0.98–1.68)0.09
       181–2400.99 (0.67–0.94)0.161.64 (1.24–2.11)0.01
       241–3001.02 (0.82–1.49)0.252.14 (1.72–2.45)<0.001
      BCS
       ≤2.5Referent
       2.75–3.250.57 (0.23–0.89)<0.001
       ≥3.50.56 (0.25–0.94)0.03
      Lameness
       Not lameReferent
       Lame1.57 (1.13–1.97)0.01
      Hock injury
       Not injuredReferent
       Injured1.78 (1.47–2.22)<0.001
      Claw length
       No overgrowthReferent
       Overgrowth1.30 (1.13–1.85)0.04
      Stall base
       ConcreteReferent
       Rubber mat1.14 (0.72–2.36)0.20
       Geomattress1.44 (1.01–1.72)0.05
       Waterbed1.07 (0.66–1.39)0.52
       Other0.60 (0.29–1.04)0.06
      Bedding depth (cm)
       >2ReferentReferent
       ≤21.50 (1.15–2.01)0.021.40 (1.18–1.61)0.04
      Slipperiness
       Less slippery (≤3% slip)Referent
       More slippery (>3% slip)2.11 (1.87–2.37)<0.001
      1 CHG = farms that made changes to the freestall area after a previous cow comfort assessment; NC = farms that did not make changes to the freestall area after a previous cow comfort assessment; NF = farms that had not been exposed to the previous cow comfort assessment.

      Hock and Knee Injuries

      Explanatory variables included in the final hock injury model included study group, DIM, lameness, and stall base (Table 6). The prevalence of hock injury was not different among study groups. Cows at 181–240 DIM and cows at >240 DIM had higher odds of a hock injury than cows at <181 DIM. Lame cows had higher odds of a hock injury than non-lame cows. Cows using stalls with a geomattress base had higher odds of a hock injury than cows using stalls with a concrete base.
      The final model for knee injury included study group, parity, stall base, bedding quantity, and slipperiness (Table 6). We observed no difference in the prevalence of knee injuries among the study groups. Cows in parity 4+ had higher odds of a knee injury than cows in parity 1. Additionally, cows housed in stalls with bedding ≤2 cm had higher odds of knee injury. Finally, cows housed in facilities with slippery floors had higher odds of knee injury (Table 6).

      Lying Time

      The final model for lying time included study group, parity, DIM, lameness, stall curb height, and stall width (Table 7). Cows in NC and NF herds had a lower average daily lying time than cows on CHG herds. Across all groups, cows in parity 3 and 4+ had a higher lying time than parity 1 and 2 cows. Cows that were 61–120, 121–180, 181–240, or 241–300 DIM had higher lying times than cows that were 0–60 DIM. The lying time of lame cows was higher than that of non-lame cows. Lying time was higher for every 1 cm over 22 cm of increased curb height. Last, wider stalls resulted in more prolonged lying time.
      Table 7Final generalized linear model for lying time including cow- and herd-level factors in 44 Alberta dairy herds with cow within-herd as the experimental unit
      VariableEstimate95% CIP-value
      Intercept10.118.55 to 11.03
      Study group
      CHG = farms that made changes to the freestall area after a previous cow comfort assessment; NC = farms that did not make changes to the freestall area after a previous cow comfort assessment; NF = farms that had not been exposed to the previous cow comfort assessment.
       CHGReferent
       NC−0.33−0.59 to −0.220.001
       NF−0.62−0.71 to −0.450.001
      Parity
       1Referent
       20.12−0.18 to 0.370.35
       30.550.36 to 1.11<0.001
       4+0.830.58 to 1.16<0.001
      DIM
       0–60Referent
       61–1200.410.18 to 0.80<0.001
       121–1801.110.81 to 1.34<0.001
       181–2401.170.86 to 1.41<0.001
       241–3001.621.44 to 1.88<0.001
      Lameness
       Not lameReferent
       Lame0.440.29 to 0.56<0.001
      Stall curb height >22 cm (per 1-cm increase)0.050.01 to 0.070.03
      Bedding depth (cm)
       ≥2Referent
       <2−0.11−0.16 to 0.070.004
      Stall width (cm)
       <114Referent
       ≥1140.410.11 to 0.750.048
      1 CHG = farms that made changes to the freestall area after a previous cow comfort assessment; NC = farms that did not make changes to the freestall area after a previous cow comfort assessment; NF = farms that had not been exposed to the previous cow comfort assessment.

      DISCUSSION

      The farms included in this study were either part of a study conducted in 2011 that assessed cow comfort in Alberta (n = 30), or had never been exposed to a cow comfort assessment (n = 14). In all statistical models, we observed no difference between the NC and NF herds. Therefore, NC groups were not biased by being exposed to a previous evaluation of cow comfort, and results were comparable to farms newly participating in 2015.
      Cows in CHG herds had a lower prevalence of lameness than cows in NC and NF herds. This was not unexpected, because many of the known risk factors for lameness have been identified as environmental factors in the freestall area (
      • Cook N.B.
      • Nordlund K.V.
      The influence of the environment on dairy cow behavior, claw health and herd lameness dynamics.
      ;
      • Chapinal N.
      • de Passillé A.M.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      • Rushen J.
      Using gait score, walking speed, and lying behavior to detect hoof lesions in dairy cows.
      ;
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). Additionally, the present study was not the first to demonstrate that monitoring and assessments can decrease lameness prevalence on farms (
      • Chapinal N.
      • Weary D.M.
      • Collings L.
      • von Keyserkingk M.A.G.
      Lameness and hock injuries improve on farms participating in an assessment program.
      ). In the present study, risk factors for lameness were similar to those reported by
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      , although the results for certain risk factors differed from this study and other literature (
      • Whay H.R.
      • Barker Z.E.
      • Leach K.A.
      • Main D.C.
      Promoting farmer engagement and activity in the control of dairy cattle lameness.
      ;
      • Chapinal N.
      • Liang Y.
      • Weary D.M.
      • Wang Y.
      • von Keyserlingk M.A.G.
      Risk factors for lameness and hock injuries in Holstein herds in China.
      ). For example, stall base was not associated with lameness prevalence. Due to a limited number of farms with specific stall bases, many stall base options were classified as “other,” possibly masking the effect of certain bases on lameness. In Canada, studies using the same standard operating procedures reported similar results in early lactation (
      • 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.
      ;
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ), although they were not able to capture the decreased prevalence of lameness mid-lactation and increased prevalence of lameness in the last 181–300 DIM. Although this was a novel result for Canadian studies, this trend has been reported elsewhere (
      • Warnick L.D.
      • Janssen D.
      • Guard C.L.
      • Gröhn Y.T.
      The effect of lameness on milk production in dairy cows.
      ).
      Cows in CHG herds had longer lying times than cows in NC and NF herds. This was consistent with our hypothesis, because factors affecting lying time (similar to lameness) include both cow factors and freestall environmental factors (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Zaffino Heyerhoff J.C.
      • Orsel K.
      Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
      ). An important factor to highlight from the lying time model was the broad range of DIM. In the present study, lying time increased with increasing DIM. This is particularly interesting when considering that lameness prevalence was higher in the last 181–300 DIM, because we know that lameness and lying time have a complicated relationship. In that regard, although prolonged lying time can demonstrate good cow comfort, it is also highly prevalent among severely lame animals (
      • Solano L.
      • Barkema H.W.
      • Pajor E.A.
      • Mason S.
      • LeBlanc S.J.
      • Nash C.G.R.
      • Orsel K.
      Associations between lying behavior and lameness in Canadian Holstein-Friesian cows housed in freestall barns.
      ). However, we observed an increase in lameness prevalence at the end of lactation, and the increased prevalence of lameness may have caused the increased lying time at this stage. However, our study was cross-sectional, and this finding should be substantiated with longitudinal studies.
      The prevalence of both hock and knee injuries was not different on CHG herds compared to NC and NF herds. This finding was in contrast to a study by
      • Chapinal N.
      • Weary D.M.
      • Collings L.
      • von Keyserkingk M.A.G.
      Lameness and hock injuries improve on farms participating in an assessment program.
      , which found a decrease in hock injuries after participation in an assessment program. We inferred that changes made to the freestall area might be less successful at improving leg injuries, even though the causal diagram specifies that stall dimensions, bedding type, base, and management affect these injuries (
      • 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.
      ). However, this study did not address how specific changes affected various outcome measures. Producers made many changes on their farm at different times, making it difficult to determine which changes, or combination of changes, were having the greatest effect. However, the most popular changes that occurred on selected farms in this study included changing the stall base to geomattress, increasing bedding quantity, and grooving crossover alleys. According to the literature, stall base is an important risk factor for leg injuries, because they can be hard and abrasive (
      • Huxley J.N.
      • Whay H.R.
      Current attitudes of cattle practitioners to pain and the use of analgesics in cattle.
      ;
      • Kielland C.
      • Ruud L.E.
      • Zanella A.J.
      • Osteras O.
      Prevalence and risk factors for skin lesions on legs of dairy cattle housed in freestalls in Norway.
      ). However, several farms changed stall base (switched to mattresses). The surface of mattresses are typically abrasive to reduce slip, and, along with rubber mats and concrete, have been known to increase the prevalence of leg injuries (
      • Huxley J.N.
      • Whay H.R.
      Current attitudes of cattle practitioners to pain and the use of analgesics in cattle.
      ). Although many farmers increased bedding depth, we did not analyze how much deeper the bedding was; we only recorded whether the bedding was above or below 2 cm. As well, we did not see many farmers change to ideal bedding types to prevent hock and knee injuries. Deep bedded straw or, more preferably, sand, decrease leg injuries (
      • Huxley J.N.
      • Whay H.R.
      Current attitudes of cattle practitioners to pain and the use of analgesics in cattle.
      ). Additionally, stall dimensions have an effect, with longer resting areas being optimal for decreasing leg injuries (
      • Kielland C.
      • Ruud L.E.
      • Zanella A.J.
      • Osteras O.
      Prevalence and risk factors for skin lesions on legs of dairy cattle housed in freestalls in Norway.
      ;
      • 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.
      ). This change was not frequently observed in study farms.
      A longitudinal study following specific changes would be beneficial for understanding the effects of specific changes on various animal-based measures of cow comfort. Such a study would require more cows to understand how each change affected the herd and would allow for a greater variety in parity and DIM. Additionally, it would be beneficial for a future study to consider how cows are housed during their dry period, because a cow may have been injured or lame before joining the milking herd. Factors associated with hoof trimming should also be considered in the future, because this study did not include changes made to hoof trimming protocols. Depending on location, temperature should also be assessed. Temperature was not assessed in this paper, because dairy cows in Alberta are not exposed to heat stress. Heat stress can have a negative effect on cow comfort by causing many physiological problems (
      • Kadzere C.T.
      • Murphy M.R.
      • Silanikove N.
      • Maltz E.
      Heat stress in lactating dairy cows: A review.
      ). Finally, it would be beneficial to use qualitative assessments to help identify unknown barriers and motivators for improvements to risk factors associated with cow comfort.
      • Leach K.A.
      • Whay H.R.
      • Maggs C.M.
      • Barker Z.E.
      • Paul E.
      • Main D.C.J.
      Working towards a reduction in cattle lameness: 1. Understanding barriers to lameness control on dairy farms.
      ,
      • Leach K.A.
      • Whay H.R.
      • Maggs C.M.
      • Barker Z.E.
      • Paul E.
      • Main D.C.J.
      Working towards a reduction in cattle lameness: 2. Understanding dairy farmers' motivations.
      identified barriers that prevent farmers from improving lameness, such as lack of labor, and identified motivators, including taking pride in having a healthy herd. However, farmers often underestimate the prevalence of animal-based measures of cow comfort, such as lameness, in their herds (
      • Whay H.R.
      • Main D.C.J.
      • Green L.E.
      • Webster A.J.F.
      Farmer perceptions of lameness prevalence.
      ). If farmers participate in cow comfort assessments and are exposed to the results, they may be motivated to improve overall herd health on their farms.

      CONCLUSIONS

      Overall, CHG herds had a lower prevalence of lameness and higher average daily lying time compared with NC and NF herds. No variables differed between the NF and NC herds, demonstrating that they were comparable groups. The CHG herds had better cow comfort than the NC and NF herds.

      ACKNOWLEDGMENTS

      The authors thank participating farmers, collaborators, technicians, students, and CanWest DHI for their support with this study. We are especially grateful for the contributions of Casey Jacobs, Colin Nichols (University of Calgary, Calgary, AB, Canada), and Raisa Engelbert (Utrecht University, Utrecht, the Netherlands) in data collection. We thank Steve Mason for establishing the database. We also thank John Kastelic for editing this manuscript. This study was funded by the Dairy Farmers of Canada (Ottawa, ON, Canada).

      REFERENCES

        • Barkema H.W.
        • von Keyserlingk M.A.G.
        • Kastelic J.P.
        • Lam T.J.G.M.
        • Luby C.
        • Roy J.-P.
        • LeBlanc S.J.
        • Keefe G.P.
        • Kelton D.F.
        Invited review: Changes in the dairy industry affecting dairy health and welfare.
        J. Dairy Sci. 2015; 98 (26342982): 7426-7445
        • Chapinal N.
        • de Passillé A.M.
        • Weary D.M.
        • von Keyserlingk M.A.G.
        • Rushen J.
        Using gait score, walking speed, and lying behavior to detect hoof lesions in dairy cows.
        J. Dairy Sci. 2009; 92: 4365-4374
        • Chapinal N.
        • Liang Y.
        • Weary D.M.
        • Wang Y.
        • von Keyserlingk M.A.G.
        Risk factors for lameness and hock injuries in Holstein herds in China.
        J. Dairy Sci. 2014; 97 (24835971): 4309-4316
        • Chapinal N.
        • Weary D.M.
        • Collings L.
        • von Keyserkingk M.A.G.
        Lameness and hock injuries improve on farms participating in an assessment program.
        Vet. J. 2014; 202 (25447801): 646-648
        • Cook N.B.
        Footbath alternatives.
        (Accessed Jan. 26, 2017.)
        • Cook N.B.
        • Bennett T.B.
        • Nordlund K.V.
        Effect of freestall surface on daily activity patterns in dairy cows with relevance to lameness prevalence.
        J. Dairy Sci. 2004; 87 (15375052): 2912-2922
        • 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 (17983785): 360-369
        • 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: 379-388
        • Green L.E.
        • Hedges V.J.
        • Schukken Y.H.
        • Blowey R.W.
        • Packington A.J.
        The impact of clinical lameness on the milk yield of dairy cattle.
        J. Dairy Sci. 2002; 85 (12362457): 2250-2256
        • Huxley J.N.
        • Whay H.R.
        Current attitudes of cattle practitioners to pain and the use of analgesics in cattle.
        Vet. Rec. 2006; 159 (17099174): 662-668
        • Ito K.
        Assessing cow comfort using lying behaviour and lameness. MSc Thesis. University of British Columbia, Vancouver, Canada2009
        • Kadzere C.T.
        • Murphy M.R.
        • Silanikove N.
        • Maltz E.
        Heat stress in lactating dairy cows: A review.
        Livest. Prod. Sci. 2002; 77: 59-91
        • Kielland C.
        • Ruud L.E.
        • Zanella A.J.
        • Osteras O.
        Prevalence and risk factors for skin lesions on legs of dairy cattle housed in freestalls in Norway.
        J. Dairy Sci. 2009; 92 (19841211): 5487-5496
        • Leach K.A.
        • Whay H.R.
        • Maggs C.M.
        • Barker Z.E.
        • Paul E.
        • Main D.C.J.
        Working towards a reduction in cattle lameness: 1. Understanding barriers to lameness control on dairy farms.
        Res. Vet. Sci. 2010; 89 (20363487): 311-317
        • Leach K.A.
        • Whay H.R.
        • Maggs C.M.
        • Barker Z.E.
        • Paul E.
        • Main D.C.J.
        Working towards a reduction in cattle lameness: 2. Understanding dairy farmers' motivations.
        Res. Vet. Sci. 2010; 89 (20413137): 318-323
        • National Farm Animal Care Council
        Code of practice for the care and handling of dairy cattle.
        (Accessed Feb. 28, 2017.)
        • Nechanitzky K.
        • Starke A.
        • Vidondo B.
        • Müller H.
        • Reckardt M.
        • Friedli K.
        • Steiner A.
        Analysis of behavioral changes in dairy cows associated with claw horn lesions.
        J. Dairy Sci. 2016; 99 (26874422): 2904-2914
        • Rushen J.
        • Haley D.
        • de Passille A.M.
        Effect of softer flooring in tie stalls on resting behavior and leg injuries of lactating cows.
        J. Dairy Sci. 2007; 90 (17638975): 3647-3651
        • Solano L.
        • Barkema H.W.
        • Pajor E.A.
        • Mason S.
        • LeBlanc S.J.
        • Nash C.G.R.
        • Orsel K.
        Associations between lying behavior and lameness in Canadian Holstein-Friesian cows housed in freestall barns.
        J. Dairy Sci. 2016; 99 (26805982): 2086-2101
        • Solano L.
        • Barkema H.W.
        • Pajor E.A.
        • Mason S.
        • LeBlanc S.J.
        • Zaffino Heyerhoff J.C.
        • Orsel K.
        Prevalence of lameness and associated risk factors in Canadian Holstein-Friesian cows housed in freestall barns.
        J. Dairy Sci. 2015; 98 (26254526): 6978-6991
        • 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 (23660141): 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 (25465625): 698-708
        • von Keyserlingk M.A.G.
        • Barrientos A.
        • Ito K.
        • Galo E.
        • Weary D.M.
        Benchmarking cow comfort on North American freestall dairies: Lameness, leg injuries, lying time, facility design, and management for high-producing Holstein dairy cows.
        J. Dairy Sci. 2012; 95 (23063152): 7399-7408
        • Wang F.X.
        • Shao D.F.
        • Li S.L.
        • Wang Y.J.
        • Azarfar A.
        • Cao Z.J.
        Effects of stocking density on behavior, productivity, and comfort indices of lactating dairy cows.
        J. Dairy Sci. 2016; 99 (26971155): 3709-3717
        • Warnick L.D.
        • Janssen D.
        • Guard C.L.
        • Gröhn Y.T.
        The effect of lameness on milk production in dairy cows.
        J. Dairy Sci. 2001; 84 (11573778): 1988-1997
        • Weary D.M.
        • Huzzey J.M.
        • von Keyserlingk M.A.G.
        Using behavior to predict and identify ill health in animals.
        J. Anim. Sci. 2009; 87 (18952731): 770-777
        • Whay H.R.
        • Barker Z.E.
        • Leach K.A.
        • Main D.C.
        Promoting farmer engagement and activity in the control of dairy cattle lameness.
        Vet. J. 2012; 193 (22892183): 617-621
        • Whay H.R.
        • Main D.C.J.
        • Green L.E.
        • Webster A.J.F.
        Farmer perceptions of lameness prevalence.
        in: Shearer J.K. Proc. 12th Int. Symp. Lameness in Ruminants, Orlando, FL. 12th Int. Symp. Lameness in Ruminants, Organizing Committee, Orlando, FL2002: 355-358
        • 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 (24359942): 173-184