Advertisement
Research Article| Volume 98, ISSUE 11, P7865-7885, November 2015

Association of bedding types with management practices and indicators of milk quality on larger Wisconsin dairy farms

Open AccessPublished:August 19, 2015DOI:https://doi.org/10.3168/jds.2015-9866

      Abstract

      The objective of this study was to identify associations of bedding type and selected management practices with bulk milk quality and productivity of larger Wisconsin dairy farms. Dairy herds (n = 325) producing ≥11,340 kg of milk daily were surveyed during a single farm visit. Monthly bulk milk SCC and total bacteria counts were obtained from milk buyers for 255 farms for a 2-yr period. Of farms with the same type of bedding in all pens during the study period, most used inorganic bedding (IB), followed by organic nonmanure bedding (OB) and manure products (MB). Almost all bulk milk total bacterial counts were <10,000 cfu/mL and total bacterial count was not associated with bedding type. Bulk milk somatic cell score (BMSCS) was least for farms using IB, varied seasonally, and was greatest in the summer. The BMSCS was reduced when new bedding was added to stalls at intervals greater than 1 wk and when teats were dried before attaching the milking unit. The BMSCS for farms using OB was reduced when bedding in the backs of stalls was removed and replaced regularly and when fewer cows with nonfunctioning mammary quarters were present. The BMSCS for farms using MB was reduced when the proportion of cows with milk discarded was less. The rolling herd average (RHA) of herds using IB was 761 and 1,153 kg greater than the RHA of herds using OB and MB, respectively. The RHA was 353 kg greater on farms where farmers understood subclinical mastitis and 965 kg greater on farms milking 3 times daily. Each 1% increase of cows with nonfunctioning mammary quarters was associated with a decrease of 57 kg of RHA. The BMSCS, proportions of cows with milk discarded and proportion of cows with nonfunctioning mammary quarters were least for herds using IB and were associated with increased productivity. Large Wisconsin dairy farms that used inorganic bedding had greater productivity and better milk quality compared with herds using other bedding types.

      Key words

      Introduction

      Demand for dairy products is driven by consumer desire for inexpensive food that is safe, high quality, and produced using socially responsible management practices (
      • Rollin B.E.
      Annual meeting keynote address: Animal agriculture and emerging social ethics for animals.
      ;
      • Fulwider W.K.
      • Grandin T.
      • Rollin B.E.
      • Engle T.E.
      • Dalsted N.L.
      • Lamm W.D.
      Survey of dairy management practices on one hundred thirteen north central and northeastern United States dairies.
      ;
      • Sapp S.G.
      • Arnot C.
      • Fallon J.
      • Fleck T.
      • Soorholtz D.
      • Sutton-Vermeulen M.
      • Wilson J.J.H.
      Consumer trust in the U.S. food system: An examination of the recreancy theorem.
      ). Bulk milk somatic cell count (BMSCC) and total bacterial count (BMTBC) are international reference indicators of raw milk quality (
      • Costello M.
      • Rhee M.
      • Bates M.P.
      • Clark S.
      • Luedecke L.O.
      • Kang D.
      Eleven-year trends of microbiological quality in bulk tank milk.
      ). Production of high quality milk from healthy animals is a requirement for market access (

      Europa. 2009. Regulation (EC) No. 853/2004 of the European Parliament and of the Council of 29 April 2004 laying down specific hygiene rules for food of animal origin (OJ L 226, 25.6.2004, p. 22) [with successive amendments and corrections consolidated]. Accessed Mar. 21, 2015. http://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=CONSLEG:2004R0853:20071114:EN:PDF

      ;

      Wisconsin Administrative Code. 2013. Chapter ATCP 60 Dairy Farms. Accessed May 13, 2014. http://www.legis.state.wi.us/rsb/code/atcp/atcp060.pdf

      ), and processors are increasingly demanding milk that meets international standards for these indicators.
      Bulk milk SCC is associated with prevalence of IMI (
      • Schukken Y.H.
      • Wilson D.J.
      • Welcome F.
      • Garrison-Tikofsky L.
      • Gonzalez R.N.
      Monitoring udder health and milk quality using somatic cell counts.
      ) and is often used as a determinant for payment of premiums by milk purchasers (
      • Rowbotham R.F.
      The value of quality milk to the processor.
      ;
      • Jayarao B.M.
      • Pillai S.R.
      • Sawant A.A.
      • Wolfgang D.R.
      • Hegde N.V.
      Guidelines for monitoring bulk tank milk somatic cell and bacterial counts.
      ). Increased premiums for milk with low BMSCC directly increases milk price and gross farm revenue. Reduced prevalence of IMI (often measured as BMSCC) is associated with increased productivity and profitability and improved animal welfare (
      • Ott S.L.
      • Novak P.R.
      Association of herd productivity and bulk-tank somatic cell counts in US dairy herds in 1996.
      ;
      • Schukken Y.H.
      • Wilson D.J.
      • Welcome F.
      • Garrison-Tikofsky L.
      • Gonzalez R.N.
      Monitoring udder health and milk quality using somatic cell counts.
      ). Somatic cells migrate to the mammary gland in response to infection, so reducing the prevalence of IMI is crucial to reducing BMSCC (
      • Deluyker H.A.
      • Gay J.M.
      • Weaver L.D.
      Interrelationships of somatic cell count, mastitis, and milk yield in a low somatic cell count herd.
      ;
      • Schukken Y.H.
      • Wilson D.J.
      • Welcome F.
      • Garrison-Tikofsky L.
      • Gonzalez R.N.
      Monitoring udder health and milk quality using somatic cell counts.
      ). Through adoption of modern milking practices, many farms have controlled transmission of contagious mastitis pathogens and have reduced the prevalence of intramammary infections caused by these organisms (
      • Makovec J.A.
      • Ruegg P.L.
      Results of milk samples submitted for microbiological examination in Wisconsin from 1994 to 2001.
      ). Currently in the United States, environmental mastitis pathogens account for the greatest proportion of IMI, whereas only 2 decades ago most IMI were caused by contagious pathogens (
      • Lago A.
      • Godden S.M.
      • Bey R.
      • Ruegg P.L.
      • Leslie K.
      The selective treatment of clinical mastitis based on on-farm culture results: I. Effects on antibiotic use, milk withholding time, and short-term clinical and bacteriological outcomes.
      ;
      • Pinzón-Sánchez C.
      • Ruegg P.L.
      Risk factors associated with short-term post- treatment outcomes of clinical mastitis.
      ;
      • Oliveira L.
      • Hulland C.
      • Ruegg P.L.
      Characterization of clinical mastitis occurring in cows on 50 large dairy herds in Wisconsin.
      ).
      Bacterial exposure at the teat end is a primary source of exposure to potential mastitis pathogens, and reducing this exposure is an important aspect of controlling environmental mastitis. Teats become contaminated with environmental bacteria through contact with bedding and other environmental risk factors (
      • Rendos J.J.
      • Eberhart R.J.
      • Kesler E.M.
      Microbial populations of teat ends of dairy cows, and bedding materials.
      ;
      • Schreiner D.A.
      • Ruegg P.L.
      Relationship between udder and leg hygiene scores and subclinical mastitis.
      ;
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      ;
      • DeVries T.J.
      • Aarnoudse M.G.
      • Barkema H.W.
      • Leslie K.E.
      • von Keyserlingk M.A.G.
      Associations of dairy cow behavior, barn hygiene, cow hygiene, and risk of elevated somatic cell count.
      ). Teats may be in direct contact with bedding materials for 12 to 14 h per day, making bedding a primary reservoir for environmental pathogens (
      • Tucker C.B.
      • Weary D.M.
      Bedding on geotextile mattresses: How much is needed to improve cow comfort?.
      ;
      • Cook N.B.
      • Bennett T.B.
      • Nordlund K.V.
      Monitoring indices of cow comfort in free-stall-housed dairy herds.
      ;
      • Hogan J.
      • Smith K.L.
      Managing environmental mastitis.
      ). The number of bacteria on teat skin has been shown to represent the level of exposure to some mastitis pathogens (
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      ), and increased exposure has been associated with increased rates of clinical mastitis (
      • Hogan J.S.
      • Smith K.L.
      • Hoblet K.H.
      • Todhunter D.A.
      • Schoenberger P.S.
      • Hueston W.D.
      • Pritchard D.E.
      • Bowman G.L.
      • Heider L.E.
      • Brockett B.L.
      • Conrad H.R.
      Bacterial counts in bedding materials used on nine commercial dairies.
      ).
      Bulk milk total bacterial count is a measure of the microbial content of milk, is associated with animal hygiene and milk quality, and is frequently used as a basis for determining milk quality premiums (
      • Pantoja J.C.F.
      • Reinemann D.J.
      • Ruegg P.L.
      Associations among milk quality indicators in raw bulk milk.
      ;
      • Velthuis A.G.J.
      • van Asseldonk M.A.P.M.
      Process audits versus product quality monitoring of bulk milk.
      ). Bacteria in raw milk can originate from IMI, contamination during milking, inadequate cleaning of milking equipment, or failure to properly cool and store milk (
      • Reinemann D.J.
      • Mein G.A.
      • Bray D.R.
      • Reid D
      • Britt J.S.
      Troubleshooting high bacteria counts in farm milk.
      ).
      For several decades the proportion of milk produced by large dairy farms in the United States has been increasing. Herds containing ≥200 milk cows currently produce 75% of all milk and those containing ≥500 cows produce 63% of US milk (

      USDA-NASS. 2014b. Quick Stats Tools. Accessed Apr. 29, 2014. http://www.nass.usda.gov/Quick_Stats/

      ). This segment of herds will continue to produce the majority of milk sold in the United States, so it is imperative that we understand the effect of housing and management decisions used on these larger farms. The objective of this study was to determine associations of bedding type and selected management practices with bulk milk quality and productivity of larger Wisconsin dairy farms.

      Materials and Methods

      Herd Recruitment and Inclusion Criteria

      A list of licensed dairy producers was purchased from the

      Wisconsin Department of Administration. 2012. Bureau of Document Services Document Sales & Distribution, Madison, WI.

      . Farms on this list had been previously categorized by size from processor records or farm visits. Inclusion criteria required farms to be producing at least 11,340 kg of milk (1/2 of a milk tanker truck) daily and to be actively selling milk from May 1, 2010, to April 30, 2012. Between June 2012 and October 2013, researchers mailed postcards to potentially eligible farms (n = 428), advising them that they would be soon be visited by a researcher. These farms were visited between June 26, 2012, and December 1, 2013, resulting in the identification of 387 eligible farms, after excluding farms not meeting inclusion criteria (n = 41). Survey administrators also excluded 27 farms after failing to make contact with farm owners after 3 or more visits. Thirty-five farmers refused to participate, resulting in participation of 325 (84.0%) of eligible farms in the survey (Table 1).
      Table 1Description of Wisconsin dairy herds enrolled and eligible for model entry
      CharacteristicFarms (n)
      Wisconsin licensed dairy herds in 201012,770
      Excluded due to producing <11,340 kg milk per day
       Initial exclusions using processor records12,342
       Exclusions based on personal interviews41
      Potential study participants387
      Excluded
       Excluded, unable to contact farm owner27
       Excluded, farmer declined participation35
      Farms that were surveyed325
       Used different bedding in different pens or over time39
      Farms with the same bedding in all pens during the study period286
       Inorganic bedding
      Inorganic bedding includes fresh sand (n=155), recycled sand (n=28), mixture of sands (n=6), washed sand (n=3), fresh sand and lime (n=2), and field-grade lime (n=1).
      195
       Organic nonmanure bedding
      Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n=45); straw (n=7); mixtures of organic beddings (n=5); straw or sawdust over lime or sand (n=3); and recycled paper products (n=2).
      62
       Manure products for bedding
      Manure based bedding includes digester solids (n=13), biosolids (n=4), composted manure (n=4), manure solids (n=4), composted digester solids (n=3), and separated solids (n=1).
      29
      Farms eligible for analysis of BMSCS by bedding type
       Farms with the same bedding in all pens during the study286
        Excluded, milk quality test results not received from processors56
       Farms included in multivariate models230
       Bedding used on enrolled farms included in analysis
        Inorganic bedding
      Inorganic bedding includes fresh sand (n=121), recycled sand (n=23), mixture of sands (n=6), washed sand (n=3), fresh sand and lime (n=2), and field-grade lime (n=1).
      156
        Organic nonmanure bedding
      Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n=34); straw (n=7); mixtures of organic beddings (n=4); straw or sawdust over lime or sand (n=2); and recycled paper products (n=2).
      49
        Manure products for bedding
      Manure based bedding includes digester solids (n=11), composted manure (n=4), biosolids (n=3), composted digester solids (n=3), manure solids (n=3), and separated solids (n=1).
      25
      Farms eligible for analysis of rolling herd average (RHA) by bedding type
       Farms with the same bedding in all pens during the study286
        Excluded, not enrolled in DHI52
        Excluded, did not report RHA7
       Farms included in analysis227
      1 Inorganic bedding includes fresh sand (n = 155), recycled sand (n = 28), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 2), and field-grade lime (n = 1).
      2 Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n = 45); straw (n = 7); mixtures of organic beddings (n = 5); straw or sawdust over lime or sand (n = 3); and recycled paper products (n = 2).
      3 Manure based bedding includes digester solids (n = 13), biosolids (n = 4), composted manure (n = 4), manure solids (n = 4), composted digester solids (n = 3), and separated solids (n = 1).
      4 Inorganic bedding includes fresh sand (n = 121), recycled sand (n = 23), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 2), and field-grade lime (n = 1).
      5 Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n = 34); straw (n = 7); mixtures of organic beddings (n = 4); straw or sawdust over lime or sand (n = 2); and recycled paper products (n = 2).
      6 Manure based bedding includes digester solids (n = 11), composted manure (n = 4), biosolids (n = 3), composted digester solids (n = 3), manure solids (n = 3), and separated solids (n = 1).

      Data Collection

      Survey of Management Practices

      Survey administrators (n = 5) were trained as a group and practiced administering questionnaires during a 1-d training session. Following training, each survey administrator observed the administration of the questionnaires by the trainer (R.F.R.) on at least 2 farms followed by being observed on at least 2 farms by the trainer before administering questionnaires on their own. During the 18 mo when the survey was conducted, survey administrators met as a group monthly to discuss progress and to avoid procedural drift.
      Five trained survey administrators (including the first author) administered a 60-question survey instrument (http://milkquality.wisc.edu/association-of-milk-quality-to-bedding/) including herd profile, inventory, and production (6 questions); bulk milk cooling and cleaning (8 questions); use of DHI, California Mastitis Test (CMT), and bacteriologic culturing of milk samples (6 questions); milking management and routine (17 questions); level of management expertise (5 questions); types and management of bedding (8 questions); and parlor design and management (10 questions). Most questions were close-ended, but some questions (such as those collecting information on brands and trade names of teat dips) were open-ended with responses categorized postsurvey.

      Collection of Milk Quality Test Results

      For each farm, survey administrators requested a release of bulk milk test results from milk buyers for the period from May 1, 2010, to April 30, 2012. Information obtained from processors included date of shipment, milk weight, BMSCC, and BMTBC. Bulk milk test results were delivered to the senior author (P.L.R.) who removed all potentially identifying information and merged the survey data with the bulk milk results. This blinding was used to avoid any potential conflict of interest due to the first author’s employment within the dairy processing industry.
      Bulk milk test results received for farms using the same bedding in all pens during the entire study (n = 230) were used to analyze associations among indicators of milk quality and bedding. Some processors reported individual BMSCC and milk weights for each delivery of milk, whereas others reported monthly milk weighted average BMSCC for each farm. When individual BMSCC results were reported, researchers calculated a milk weighted BMSCC for each month. Seasons were defined as spring (March to May), summer (June to August), fall (September to November), and winter (December to February). For each farm-season combination, a simple average of the 3 monthly BMSCC results was calculated. Seasonal average BMSCC was transformed into bulk milk somatic cell score (BMSCS) using the equation: BMSCS = base 2 logarithm (BMSCC in cells per mL/100,000) + 3 (
      • Shook G.E.
      Genetic improvement of mastitis through selection on somatic cell count.
      ), and used for analysis.
      Some processors reported individual BMTBC results for each milk pickup, whereas others reported a single BMTBC for each farm each month. For farms reporting >1 monthly BMTBC, one BMTBC was randomly selected from all BMTBC for each farm-month combination. Based on the 85th percentile of total bacterial count, monthly BMTBC were categorized as good (≤10,000 cfu/mL) or bad (>10,000 cfu/mL).

      Statistical Analysis

      Participating farms were categorized by the type of bedding used: inorganic bedding (IB), manure products (MB), other organic bedding (OB), or different bedding in different pens or at different times during the 2-yr study period (Table 1).
      Statistical analyses were performed using SAS 9.4 (
      SAS Institute
      SAS/STAT User’s Guide 9.4: Statistics..
      ). Farm was the unit of analysis. Descriptive statistics using PROC MEANS, PROC UNIVARIATE, PROC FREQ, and PROC SGPLOT were used to verify data integrity, observe frequency distributions, and identify missing data. Biologically relevant categories of independent variables with <5 responses were combined for analysis. Continuous variables were stratified by type of bedding and analyzed using PROG GLM or PROC NPAR1WAY (when assumptions of linear models were not met). Multiple comparisons were adjusted using Tukey-Kramer for linear models and the Dwass, Steel, Critchlow-Fligner method for nonparametric models. Categorized variables were stratified by type of bedding and analyzed using PROC FREQ using a chi-squared test or Fisher’s exact test (when the expected value of any cell in the contingency table was <5). All means are presented as mean ± SE.

      Development of Models

      All multivariate models were built using a manual process that incorporated the biological and statistical relevance of each variable. Initially, all biologically relevant variables (survey responses) were assessed by screening for unconditional associations between each survey response and the relevant outcome variable. For each model, risk factors that were unconditionally associated with the outcome variable at P ≤ 0.25 were tested for associations among predictor variables using chi-squared (for categorical variables) or correlation analysis (for continuous variables). If variables were highly associated (R2 ≥ 0.60), the most biologically relevant variable was selected for further analysis. For each model, risk factors that were unconditionally associated with the outcome variable at P ≤ 0.25 were offered for further multivariate modeling. Both forward and backward variable selection procedures were used to select the variables that remained in the final models. Confounding was assessed by examining the effect of each variable on the estimates of other explanatory variables. No variables included in any final model resulted in substantial changes among estimates of other explanatory variables (≥20%;
      • Dohoo I.
      • Martin W.
      • Stryhn H.
      ), indicating that confounding was not a problem. Biologically relevant first-order interactions among variables were offered for backward and forward variable selection to construct the final multivariate regression models. The Akaike information criterion (AIC;
      • Akaike H.
      Fitting autoregressive models for prediction.
      ) was chosen for comparing models based on goodness of fit. All P-values for multiple comparisons within models were Tukey adjusted.

      Models for Bulk Milk Somatic Cell Score

      Initial modeling demonstrated that bedding type was associated with BMSCS (Table 2) and bedding management practices differed by bedding type, so 3 separate models were used to determine associations of bedding type (IB, OB, MB) with possible explanatory variables. Associations between selected survey responses (explanatory variables; Table 3) and BMSCS (response variable) were tested using repeated measures mixed models using PROC MIXED with an autoregressive covariance structure that included the intercept, year as a random effect, season, and the explanatory variable of interest. Explanatory variables unconditionally associated with BMSCS, and biologically relevant first-order interactions were used to develop one multivariate repeated measures mixed model for each of the 3 bedding types.
      Table 2Adjusted estimates of mean bulk milk SCS
      Bulk milk somatic cell score=log2(bulk milk SCC/100,000) + 3.
      for all explanatory variables that remained in the final multivariate mixed repeated measures regression model for Wisconsin dairy farms (n = 222) producing ≥11,340 kg of milk and using the same bedding in all pens during the study
      Explanatory variablenEstimateSEType III95% CISCC
      Back transformation to SCC, not used for analysis.
      P-valueLowerUpper(1,000 cells/mL)
      Intercept3.700.107162
      Cows with nonfunctioning mammary quarters, %2220.01420.00660.0300.00130.0271
      Categorical distributionLevels
       Year (random)12220.050.050−0.050.14
      2222−0.050.050−0.140.05
       BeddingIB
      Inorganic bedding.
      1523.90
      Values within a row with different superscripts differ (P<0.05).
      0.080<0.0013.754.06187
      MB
      Manure bedding.
      234.19
      Values within a row with different superscripts differ (P<0.05).
      0.0813.964.42228
      OB
      Organic, nonmanure bedding.
      474.20
      Values within a row with different superscripts differ (P<0.05).
      0.0804.024.39230
       SeasonSpring2224.02
      Values within a row with different superscripts differ (P<0.05).
      0.081<0.0013.864.18203
      Summer2224.25
      Values within a row with different superscripts differ (P<0.05).
      0.0804.094.41238
      Fall2224.12
      Values within a row with different superscripts differ (P<0.05).
      0.0803.964.28217
      Winter2224.01
      Values within a row with different superscripts differ (P<0.05).
      0.0803.854.17201
       Type of postdipIodine1584.18
      Values within a row with different superscripts differ (P<0.05).
      0.0810.0154.024.33227
      Noniodine644.02
      Values within a row with different superscripts differ (P<0.05).
      0.0913.844.20203
       Written milking protocolPresent1644.24
      Values within a row with different superscripts differ (P<0.05).
      0.087<0.0013.784.03236
      Absent583.96
      Values within a row with different superscripts differ (P<0.05).
      0.0934.064.53195
       Wipe off predip before milkingYes2053.90
      Values within a row with different superscripts differ (P<0.05).
      0.066<0.0013.784.03187
      No174.24
      Values within a row with different superscripts differ (P<0.05).
      0.0874.064.53236
      Interaction
       WMP × beddingPresent − IB1204.09
      Values within a row with different superscripts differ (P<0.05).
      0.078<0.0013.944.24213
      Present − MB154.50
      Values within a row with different superscripts differ (P<0.05).
      0.1444.224.78283
      Present − OB294.14
      Values within a row with different superscripts differ (P<0.05).
      0.1043.934.34220
      Absent − IB323.72
      Values within a row with different superscripts differ (P<0.05).
      0.1013.523.91165
      Absent − MB83.89
      Values within a row with different superscripts differ (P<0.05).
      0.1603.584.20185
      Absent – OB184.27
      Values within a row with different superscripts differ (P<0.05).
      0.1254.034.52241
      a–c Values within a row with different superscripts differ (P < 0.05).
      1 Bulk milk somatic cell score = log2(bulk milk SCC/100,000) + 3.
      2 Back transformation to SCC, not used for analysis.
      3 Inorganic bedding.
      4 Manure bedding.
      5 Organic, nonmanure bedding.
      Table 3Variables assessed for unconditional association with bulk milk SCS
      VariableTypeCategories, nResponses for categorical variables
      Rolling herd average, kgContinuous
      Cows with milk discarded, %Continuous
      Cows milking on fewer than 4 quarters, %Continuous
      SeasonCategorical4Spring, summer, fall, winter
      Farm size, cowsCategorical4200 to 499, 500 to 999, 1,000 to 1,999, ≥2,000
      Parlor typeCategorical3Parallel, rotary, herringbone or other
      Frequency of DHI testingCategorical2Monthly, less than once monthly
      CMT
      California Mastitis Test.
      used on farm
      Categorical2Yes, no
      Frequency of CMT useCategorical4Daily, >1× weekly, occasionally, never
      CMT used on fresh cowsCategorical2Most or all, few or none
      CMT used on cows with clinical mastitisCategorical2Most or all, few or none
      CMT used on cows with elevated SCCCategorical2Most or all, few or none
      Manager presence in parlor during milkingCategorical2Always or portions of all milkings, only when problems are suspected or never
      Written milking protocolCategorical2Present, absent
      Type of predipCategorical2Iodine, other
      Type of postdipCategorical2Iodine, other
      Frequency of forestrippingCategorical2Always, sometimes or never
      Dry teats before attaching milking unitCategorical2Yes, no
      Use mechanical teat scrubberCategorical2Yes, no
      Understanding of subclinical mastitisCategorical2Understand, limited knowledge
      Segregate cows with elevated SCCCategorical2Yes, no
      Culture bulk milkCategorical2≥Once monthly, less than monthly
      Culture fresh cowsCategorical2Most or all, few or none
      Culture cows with clinical mastitisCategorical2Most or all, few or none
      Culture cows with elevated SCCCategorical2Most or all, few or none
      Bedding depthCategorical2Deep, shallow
      Frequency of adding bedding, hContinuous
      Frequency of adding beddingCategorical2More than once weekly, every 7 d or more
      Regularly remove and replace some bedding
      All bedding in the back one-third of stalls removed and replaced at least once weekly.
      Categorical2Yes, no
      Regularly remove and replace all bedding
      Farmers indicated that all bedding was removed from stalls and replaced at regular intervals.
      Categorical2Yes, no
      Type of bedding groomingCategorical2Rake, mechanical
      Use mechanical groomerCategorical2Yes, no
      Frequency of grooming bedding in backs of stalls, hContinuous
      Type of sand (inorganic bedding model only)Categorical2Fresh, recycled
      1 California Mastitis Test.
      2 All bedding in the back one-third of stalls removed and replaced at least once weekly.
      3 Farmers indicated that all bedding was removed from stalls and replaced at regular intervals.

      Model for Rolling Herd Average

      Associations between selected survey responses (explanatory variables; Table 4) and rolling herd average (RHA; response variable) were tested using a univariate regression model using PROC MIXED, which included the intercept and the explanatory variable of interest. Explanatory variables unconditionally associated with RHA and biologically relevant first-order interactions were used to develop the final multivariate model for RHA.
      Table 4Variables assessed for unconditional association with rolling herd average
      VariableTypeCategories, nResponses for categorical variables
      Cows with milk discarded, %Continuous
      Cows milking on fewer than 4 quarters, %Continuous
      Frequency of adding bedding, hContinuous
      Frequency of grooming back 1/3 of stalls, hContinuous
      Type of beddingCategorical3Inorganic, manure-based, organic nonmanure
      Milkings per dayCategorical22×, 3×
      Farm size, cowsCategorical4200 to 499, 500 to 999, 1,000 to 1,999, ≥2,000
      Parlor typeCategorical3Parallel, rotary, herringbone or other
      CMT
      California Mastitis Test.
      used on farm
      Categorical2Yes, no
      Manager presence in parlor during milkingCategorical2Always or portions of all milkings, only when problems are suspected or never
      Written milking protocolCategorical2Present, absent
      Type of predipCategorical2Iodine, other
      Type of postdipCategorical2Iodine, other
      Frequency of forestrippingCategorical2Always, sometimes or never
      Complete milking routineCategorical2Yes, no
      Parlor cleanliness after milkingCategorical3Clean, neither clean nor dirty, dirty
      Understanding of subclinical mastitisCategorical2Understand, limited knowledge
      Segregate cows with elevated SCCCategorical2Yes, no
      Culture bulk milkCategorical2≥Once monthly, less than monthly
      Culture fresh cowsCategorical2Most or all, few, or none
      Culture cows with clinical mastitisCategorical2Most or all, few, or none
      Culture cows with elevated SCCCategorical2Most or all, few, or none
      Bedding depthCategorical2Deep, shallow
      1 California Mastitis Test.

      Model for Bulk Milk Total Bacteria Count

      Association between bedding type (explanatory variable) and categorized BMTBC (response variable) was tested using repeated measures mixed logistic model using PROC GLIMMIX with an autoregressive covariance structure that included the intercept, year as a random effect, month, and bedding type (IB, OB, MB).

      Results

      Characteristics of Surveyed Farms

      Milk Production

      Herds (n = 325) milked between 255 and 8,100 lactating cows (868 ± 44) and contained a total of 282,235 lactating cows. Farms with ≥2,000 milking cows (n = 20) represented 6.2% of farms, 23.5% of cows, and 24.0% of milk. Farms with 1,000 to 1,999 milking cows (n = 61) represented 18.8, 28.7, and 28.4% of the farms, cows, and milk, respectively. Farms with 500 to 999 milking cows (n = 140) represented 43.1, 33.1, and 32.5% of the farms, cows, and milk, respectively. The remaining 104 farms had <500 milking cows and represented 32.0, 14.8, and 15.1% of the farms, cows, and milk, respectively.
      Daily milk production ranged from 11,340 to 317,151 kg (32,149 ± 1.712 kg) per farm with 37.0 ± 0.3 kg of milk sold per cow per day (Table 5). Milk was sold from 98.2 ± 0.07% of lactating cows (Table 5). Individual cow production was tested through DHI on 81.5% (n = 265) of herds. Testing frequency varied from 1- to 12-mo intervals (1.8 ± 0.16). Rolling herd average on farms enrolled in DHI (n = 257) ranged from 8,165 to 15,422 kg (12,584 ± 75 kg) with 8 farmers enrolled in DHI testing not reporting RHA.
      Table 5Characteristics of Wisconsin dairy farms (n = 325) producing ≥11,340 kg of milk daily and participating in survey
      OutcomeFarmsMinimumMeanMaximumSE
      Mean herd size, n3252558688,10044.28
      Total daily production, kg32511,34032,149317,5151,712
      DHI testing frequency,
      60 herds did not participate in DHI.
      mo
      26511.8120.16
      Rolling herd average,
      8 of 265 farmers participating in DHI did not report rolling herd average.
      kg
      2578,16512,58415,42275.33
      Milk sold per cow per day, kg32521.537.049.90.26
      Proportion of cows with milk not sold, %3250.01.86.70.07
      Cows milking on fewer than 4 quarters,
      13 farmers did not report percent of cows with nonfunctioning mammary quarters.
      %
      3120.04.730.00.23
      Times bulk milk cultured annually,
      62 farmers did not culture bulk milk.
      n
      263120.33653.09
      Milk pickup temperature, °C32503.670.05
      Milk pickup frequency, h325219.5480.50
      Bulk tank or trailer wash frequency, h325221.6480.52
      1 60 herds did not participate in DHI.
      2 8 of 265 farmers participating in DHI did not report rolling herd average.
      3 13 farmers did not report percent of cows with nonfunctioning mammary quarters.
      4 62 farmers did not culture bulk milk.

      Milking Management

      The majority of farms had a written milking routine that was available in the first language of the milking technicians (Table 6). Most farmers had adopted a complete milking routine consisting of premilking teat sanitation, removal of foremilk, drying teats with individual towels, and application of germicidal postdip (Table 6). Liquid pre- and postteat dips were used by 75 and 97% of farmers, respectively, with the remaining farms using foam teat dips. Few (10%) farmers used teat scrubbers, and of those, only 18% dried teats before machine attachment (compared with 99.7% of other farms). Most managers were present in the milking parlor during portions of milking shifts, but some managers did not supervise milking during any shifts (Table 6).
      Table 6Milking routine and mastitis management practices on Wisconsin dairy farms (n = 325) producing ≥11,340 kg of milk daily and participating in survey
      PracticeTotalWithin category
      n%n%
      Predip use, type, and practices
       Always predip32299.1
       Never predip30.9
      Drying of teats by farmers using predip
       Do not dry teats before machine attachment268.1
       Dry teats before machine attachment29691.9
      Use a teat scrubber
       Yes3310.2
       No29289.8
      Forestrip as part of milking routing
       Always28186.5
       Sometimes298.9
       Never154.6
      Postdip use, type, and practices
       Always postdip32499.7
       Sometimes postdip10.3
      Complete milking routine
      Defined as use of a milking routine that always includes forestrip, predip, dry teats, and postdip.
      26380.9
      Written milking protocol
       Yes24073.8
       No8526.2
      Manager in the parlor during milking
       Always3611.1
       Portions of every milking shift9328.6
       Only during daytime milkings11936.6
       Only when problems are suspected206.2
       Occasionally or never5717.5
      DHIA enrollment
       Used DHI testing26581.5
       Did not use DHI testing6018.5
      1 Defined as use of a milking routine that always includes forestrip, predip, dry teats, and postdip.

      Management of Mastitis

      Most farms had a written protocol for mastitis treatments (62.5%). When asked an open-ended question about how they identify subclinical mastitis (SM), 73.5% (n = 239) of all farmers demonstrated an understanding of SM without further prompting; of these, 31 did not check cows for SM. Of the 86 farmers who did not demonstrate knowledge of SM without prompting, 71 demonstrated understanding when prompted with common examples of methods used to detect SM; of these farmers, 48 indicated that they use individual cow SCC to detect SM. Two farmers believed that they did not have cows with SM present on their farms. Almost all farms enrolled in DHI (97.7%) received individual cow SCC results on each test date. Cows with high SCC were segregated during milking on 32% of farms (n = 104) with 75% of those farms using separate pens, 17.3% using hospital barns, and 7.7% using separate milking units. On farms with separate pens for cows with high SCC, 94.9% milked these pens last.
      Many farmers used bacteriological culturing to identify mastitis pathogens in both bulk milk and milk samples from individual cows. Approximately one-half (53%) of farmers received bulk milk cultures once or more monthly, whereas 66% regularly received culture results from individual cow samples. Very few farmers used bacteriological culture for all fresh cows (7%) or cows with elevated SCC (9%), whereas about one-third of farmers received bacteriological culture results from all cows detected to have clinical mastitis.
      The percent of cows with nonfunctioning mammary quarters was reported by 96.0% of surveyed farmers (n = 312). The proportion of cows with nonfunctioning mammary quarters ranged from 0 to 30%. The 25th, 50th, 75th, and 90th percentiles of nonfunctioning mammary quarters were 2, 4, 6, and 10%, respectively.

      Bedding Type and Management

      Inorganic bedding was used most commonly (60% of farms) followed by OB (19%) and MB (9%). Twelve percent of surveyed farms used different bedding in different pens or at different times during the study. Bedding on farms using the same type of bedding in all pens during the entire study period included IB: fresh sand (n = 121), recycled sand (n = 23), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 2), and field-grade lime (n = 1); OB: sawdust, shavings, and grain hulls (n = 34), straw (n = 7), mixtures of organic beddings (n = 4), straw or sawdust over lime or sand (n = 2), and recycled paper products (n = 2); and MB: digester solids (n = 11), composted manure (n = 4), biosolids (n = 3), composted digester solids (n = 3), manure solids (n = 3), and separated solids (n = 1; Table 1). Farms on which MB was used contained more lactating cows (1,502 ± 187) than farms using IB (849 ± 62) or OB (706 ± 53; P < 0.001) and had a greater proportion of cows with milk not sold (2.4 ± 0.30%) than farms bedding with IB (1.7 ± 0.08%; P = 0.005; Table 7). A tendency was observed for farms using MB to have a greater proportion of cows with nonfunctioning mammary quarters (P = 0.055; Table 7). As compared with farms with MB or OB, farms with IB had greater milk sold per cow per day (P = < 0.001), greater RHA (P < 0.001), and a lesser proportion of cows with milk not sold (Table 7).
      Table 7Characteristics of Wisconsin dairy farms (n = 286) producing ≥11,340 kg of milk daily and using the same bedding type during the study period stratified by bedding type
      OutcomenBedding during study periodP-value
      IB inorganic
      Inorganic bedding includes fresh sand (n=155), recycled sand (n=28), mixture of sands (n=6), washed sand (n=3), fresh sand and lime (n=2), and field-grade lime (n=1).
      OB organic nonmanure
      Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n=45); straw (n=7); mixtures of organic beddings (n=5); straw or sawdust over lime or sand (n=3); and recycled paper products (n=2).
      MB manure products
      Manure based bedding includes digester solids (n=13), biosolids (n=4), composted manure (n=4), manure solids (n=4), composted digester solids (n=3), and separated solids (n=1).
      Farms,
      39 herds used different bedding in different pens or changed bedding over time.
      n
      2861956229
      Mean herd size, n286849
      Values within a row with different superscripts differ (P<0.05).
      706
      Values within a row with different superscripts differ (P<0.05).
      1502
      Values within a row with different superscripts differ (P<0.05).
      <0.001
      Total daily production, kg28632,331
      Values within a row with different superscripts differ (P<0.05).
      24,492
      Values within a row with different superscripts differ (P<0.05).
      53,175
      Values within a row with different superscripts differ (P<0.05).
      <0.001
      DHI testing frequency,
      Includes IB (n=161), OB (n=47), MB (n=26).
      mo
      2341.72.42.00.275
      Rolling herd average,
      Includes IB (n=159), OB (n=44), MB (n=24).
      kg
      22712,870
      Values within a row with different superscripts differ (P<0.05).
      12,025
      Values within a row with different superscripts differ (P<0.05).
      11,779
      Values within a row with different superscripts differ (P<0.05).
      <0.001
      Milk sold per cow per day, kg28638.0
      Values within a row with different superscripts differ (P<0.05).
      34.9
      Values within a row with different superscripts differ (P<0.05).
      35.4
      Values within a row with different superscripts differ (P<0.05).
      <0.001
      Proportion of cows with milk not sold, %2861.6
      Values within a row with different superscripts differ (P<0.05).
      1.9
      Values within a row with different superscripts differ (P<0.05).
      2.4
      Values within a row with different superscripts differ (P<0.05).
      0.005
      Cows with nonfunctioning mammary quarters,
      Includes IB (n=190), OB (n=61), MB (n=24).
      %
      2754.54.86.30.055
      Bulk milk somatic cell count (BMSCC),
      Includes IB (n=156), OB (n=49), MB (n=25).
      cells/mL
      230198
      Values within a row with different superscripts differ (P<0.05).
      220
      Values within a row with different superscripts differ (P<0.05).
      248
      Values within a row with different superscripts differ (P<0.05).
      0.004
      Bulk milk total bacteria count (BMTBC),
      Includes IB (n=156), OB (n=49), MB (n=25); percent of BMTBC ≤10,000 cfu/mL.
      %
      23078.275.176.10.088
      Times bulk milk cultured annually,
      Includes IB (n=129), OB (n=36), MB (n=22).
      n
      18722.212.638.50.565
      Milk pickup temperature, °C2863.63.73.30.150
      Milk pickup frequency, h28619.6
      Values within a row with different superscripts differ (P<0.05).
      21.6
      Values within a row with different superscripts differ (P<0.05).
      13.3
      Values within a row with different superscripts differ (P<0.05).
      <0.001
      Bulk tank or trailer wash frequency, h28621.8
      Values within a row with different superscripts differ (P<0.05).
      22.8
      Values within a row with different superscripts differ (P<0.05).
      14.7
      Values within a row with different superscripts differ (P<0.05).
      <0.001
      a,b Values within a row with different superscripts differ (P < 0.05).
      1 Inorganic bedding includes fresh sand (n = 155), recycled sand (n = 28), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 2), and field-grade lime (n = 1).
      2 Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n = 45); straw (n = 7); mixtures of organic beddings (n = 5); straw or sawdust over lime or sand (n = 3); and recycled paper products (n = 2).
      3 Manure based bedding includes digester solids (n = 13), biosolids (n = 4), composted manure (n = 4), manure solids (n = 4), composted digester solids (n = 3), and separated solids (n = 1).
      4 39 herds used different bedding in different pens or changed bedding over time.
      5 Includes IB (n = 161), OB (n = 47), MB (n = 26).
      6 Includes IB (n = 159), OB (n = 44), MB (n = 24).
      7 Includes IB (n = 190), OB (n = 61), MB (n = 24).
      8 Includes IB (n = 156), OB (n = 49), MB (n = 25).
      9 Includes IB (n = 156), OB (n = 49), MB (n = 25); percent of BMTBC ≤10,000 cfu/mL.
      10 Includes IB (n = 129), OB (n = 36), MB (n = 22).
      The frequency of adding new bedding to stalls varied by bedding type and depth of bedding in the stalls. Farms with IB were 15.7 and 1.8 times more likely to have deep bedding as compared with farms with OB and MB, respectively (P < 0.001). For stalls containing mattresses as a base, new MB or OB was added to stalls approximately twice as frequently as compared with IB (Table 8; P < 0.001). Farms with IB added bedding to stalls less frequently (median of 168 h) than farms with OB or MB (medians of 78 and 72 h, respectively). Farms which used OB, removed and replaced bedding in the backs of stalls, or completely removed and replaced all bedding in stalls more often than those using IB (P < 0.001; Table 8). For deep-bedded stalls, farmers added MB approximately twice as frequently as IB (P < 0.001). For stalls with tires or sand traps, farmers added IB less frequently than in stalls without a base bedding or mattress (P = 0.023; Table 8).
      Table 8Bedding management practices on Wisconsin dairy farms (n = 286) producing ≥ 11,340 kg of milk daily and using the same type of bedding for the study period stratified by bedding type, h
      ItemInorganic,
      Inorganic bedding includes fresh sand (n=155), recycled sand (n=28), mixture of sands (n=6), washed sand (n=3), fresh sand and lime (n=2), and field-grade lime (n=1).
      mean ± SE (n)
      Manure products,
      Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n=45); straw (n=7); mixtures of organic beddings (n=5); straw or sawdust over lime or sand (n=3); and recycled paper products (n=2).
      mean ± SE (n)
      Organic nonmanure,
      Manure based bedding includes digester solids (n=13), biosolids (n=4), composted manure (n=4), manure solids (n=4), composted digester solids (n=3), and separated solids (n=1).
      mean ± SE (n)
      P-value for row
      Grooming back of stalls,
      9 did not groom backs of stalls.
      h
      10 ± 0.9 (188)12 ± 2.1 (28)10 ± 1.5 (61)0.754
      Mechanical grooming,
      188 did not use mechanical bedding groomers.
      h
      94 ± 12.5 (79)49 ± 24.4 (18)59 ± ‡ (1)
      Excluded from analysis due to sparse data.
      0.109
      Remove and replace back one-third of stall,
      100 did not remove and replace bedding in back one-third of stalls.
      h
      416
      Values within a row with different superscripts differ (P<0.05).
      ± 60.5 (55)
      168
      Values within a row with different superscripts differ (P<0.05).
      ± 139.1 (23)
      39
      Values within a row with different superscripts differ (P<0.05).
      ± 85.7 (23)
      <0.001
      Remove and replace all,
      259 did not remove and replace all bedding in stalls.
      h
      450
      Values within a row with different superscripts differ (P<0.05).
      ± 65.6 (11)
      205
      Values within a row with different superscripts differ (P<0.05).
      ± 80.3 (10)
      41
      Values within a row with different superscripts differ (P<0.05).
      ± 60.7 (10)
      <0.001
      Bedding addition, h161
      Values within a row with different superscripts differ (P<0.05).
      ± 5.2 (195)
      71
      Values within a row with different superscripts differ (P<0.05).
      ± 12.3 (62)
      79
      Values within a row with different superscripts differ (P<0.05).
      ± 8.9 (29)
      <0.001
      Frequency of adding new bedding stratified by base bedding material and bedding management practices
       Base bedding
        Mattress133
      Values within a row with different superscripts differ (P<0.05).
      ,
      Values within a column within a category with different superscripts differ (P<0.05).
      ± 11.6 (21)
      67
      Values within a row with different superscripts differ (P<0.05).
      ± 13.1 (15)
      76
      Values within a row with different superscripts differ (P<0.05).
      ± 8.3 (58)
      <0.001
        Tires or sand traps201
      Values within a column within a category with different superscripts differ (P<0.05).
      ± 32.5 (18)
      A dash indicates no farms within category.
      192 ± 168.0 (2)
      Excluded from analysis due to sparse data.
        None160
      Values within a row with different superscripts differ (P<0.05).
      ,
      Values within a column within a category with different superscripts differ (P<0.05).
      ± 5.5 (156)
      75
      Values within a row with different superscripts differ (P<0.05).
      ± 11.5 (14)
      40 ± 16.0 (2)
      Excluded from analysis due to sparse data.
      <0.001
        P-value for column0.0230.374
       Method of replacing bedding
        Never replace all bedding163
      Values within a row with different superscripts differ (P<0.05).
      ± 7.2 (139)
      80
      Values within a row with different superscripts differ (P<0.05).
      ± 10.7 (21)
      88
      Values within a row with different superscripts differ (P<0.05).
      ± 11.1 (39)
      <0.001
        Replace in backs of stalls162
      Values within a row with different superscripts differ (P<0.05).
      ± 8.6 (44)
      54 ± 30.0 (2)
      Excluded from analysis due to sparse data.
      47
      Values within a row with different superscripts differ (P<0.05).
      ± 12.8 (13)
      <0.001
        Replace all bedding129
      Values within a row with different superscripts differ (P<0.05).
      ± 15.1 (11)
      42
      Values within a row with different superscripts differ (P<0.05).
      ± 9.7 (6)
      86
      Values within a row with different superscripts differ (P<0.05).
      ± 31.5 (10)
      0.007
        Differs by pen or over time120‡ (1)
      Excluded from analysis due to sparse data.
        P-value for column0.3840.7920.087
       Type of grooming
        Rake173
      Values within a row with different superscripts differ (P<0.05).
      ,
      Values within a column within a category with different superscripts differ (P<0.05).
      ± 8.2 (114)
      69
      Values within a row with different superscripts differ (P<0.05).
      ± 16.7 (11)
      80
      Values within a row with different superscripts differ (P<0.05).
      ± 9.2 (61)
      <0.001
        Mechanical groomer143
      Values within a row with different superscripts differ (P<0.05).
      ,
      Values within a column within a category with different superscripts differ (P<0.05).
      ± 6.5 (78)
      7
      Values within a row with different superscripts differ (P<0.05).
      ± 9.9 (18)
      8
      Standard error not reported for single observations.
      (1)
      Excluded from analysis due to sparse data.
      <0.001
        Do not groom bedding168 ± 0.0 (3)
      Excluded from analysis due to sparse data.
        P-value for column0.0080.647
      a,b Values within a row with different superscripts differ (P < 0.05).
      y,z Values within a column within a category with different superscripts differ (P < 0.05).
      1 Inorganic bedding includes fresh sand (n = 155), recycled sand (n = 28), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 2), and field-grade lime (n = 1).
      2 Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n = 45); straw (n = 7); mixtures of organic beddings (n = 5); straw or sawdust over lime or sand (n = 3); and recycled paper products (n = 2).
      3 Manure based bedding includes digester solids (n = 13), biosolids (n = 4), composted manure (n = 4), manure solids (n = 4), composted digester solids (n = 3), and separated solids (n = 1).
      4 9 did not groom backs of stalls.
      5 188 did not use mechanical bedding groomers.
      6 100 did not remove and replace bedding in back one-third of stalls.
      7 259 did not remove and replace all bedding in stalls.
      8 A dash indicates no farms within category.
      Standard error not reported for single observations.
      * Excluded from analysis due to sparse data.
      Frequency of grooming bedding in the back of stalls was not associated with bedding type (P = 0.754) presence or type of base bedding materials (P ≥ 0.667), or type of stall grooming (P ≥ 0.062). For all types of bedding, the median frequency of stall grooming was 8 h (corresponding to milking frequency). Mechanical groomers were used on 41% of farms with IB in contrast to 29 and 3% of farms using MB and OB, respectively (Table 6). No difference was found in the frequency of mechanical grooming by type of bedding (P ≥ 0.087), deep-bedded versus shallow-bedded stalls (P ≥ 0.579), or the frequency of removing and replacing bedding in stalls (P ≥ 0.087). For farms using the same bedding in all pens during the entire study period, no association was observed of bedding type with frequency of DHI testing, milk pickup temperature, or the frequency bacteriological cultures of bulk milk (P ≥ 0.15; Table 7).

      Analyses of BMSCS by Bedding Type

      For the 2-yr study period, processors reported a median of 656 BMSCC results per farm. During the study period, all farms were in compliance with US SCC regulations. For 212 of 230 farms (92.5%) using the same bedding in all pens during the study, all 3-mo geometric mean BMSCC were under the EU limit of 400,000 cells/mL.
      Based on relationships between explanatory variables and BMSCS, variables eligible for inclusion in the 3 multivariate models (IB, OB, MB) for BMSCS included season (eligible for IB, OB, MB), proportion of cows with milk discarded or not sold (eligible for IB, MB), proportion of cows with nonfunctioning mammary quarters (eligible for IB, OB, MB), manager presence during milking (eligible for IB, OB), written milking protocol (WMP; eligible for IB, OB, MB), type of predip (eligible for IB), type of postdip (eligible for IB, OB), dry teats (eligible for IB), segregate cows with elevated SCC (eligible for IB, MB), culture bulk milk (eligible for IB), culture cows with clinical mastitis (eligible for IB), bedding depth (eligible for IB), frequency of adding bedding (eligible for IB), removal and replacement of bedding in back one-third of stalls (eligible for OB), type of bedding grooming (eligible for IB), and type of sand (eligible for IB).

      Inorganic Bedding

      Data from 151 farms using IB were used to determine associations of eligible management factors and BMSCS. Five of 156 farms originally considered eligible for this model were excluded because 2 farms did not groom bedding and 3 farms did not have complete BMSCS results. Of explanatory variables eligible for entry into the BMSCS model, season, frequency of adding bedding, WMP, drying teats before attaching the milking unit, type of grooming of bedding, and depth of bedding remained in the final model. No interactions remained in the final model (Table 9). Bulk milk SCS was least in winter and spring, intermediate in fall, and greatest in summer (P < 0.001). The BMSCC for farms using IB was 0.15 log2 units greater when bedding was added more than once weekly, 0.35 log2 units greater when teats were not dried before machine attachment, and 0.35 log2 units lesser when a WMP was not present on the farm (Table 9). All 13 farms with IB which did not dry teats before milking unit attachment used teat scrubbers. Of farms using IB and teat scrubbers (n = 14), only 1 farm also dried teats before milking unit attachment. A tendency was found for greater BMSCS when stalls were deep bedded (P = 0.065). Bulk milk SCS was not associated with type of grooming of bedding (P = 0.74), but model fit was improved with the inclusion of this variable (AIC = 16.4 vs. 24.9). All model assumptions were satisfied.
      Table 9Adjusted estimates of mean bulk milk somatic cell score for all explanatory variables that remained in the final multivariate mixed repeated measures regression model for Wisconsin dairy farms (n = 151) producing ≥11,340 kg of milk daily and using inorganic
      Inorganic bedding includes fresh sand (n=120), recycled sand (n=20), mixture of sands (n=6), washed sand (n=3), fresh sand and lime (n=1), and field-grade lime (n=1).
      bedding in all pens during the study period (linear score)
      Explanatory variableLevelsFarms, nEstimateSEType III P-value95% CISCC
      Back transformation to SCC not used for analysis.
      LowerUpper(1,000 cells/mL)
      Intercept3.640.106156
      Year (random)11510.060.062−0.060.18
      2151−0.060.062−0.180.06
      Season
      Spring (March – May), summer (June – August), fall (September – November), and winter (December – February).
      Spring1513.82
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.095<0.0013.734.11177
      Summer1514.05
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0953.644.01207
      Fall1513.91
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0953.874.24188
      Winter1513.81
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0953.634.00175
      Frequency of adding beddingEvery 7 or more d1073.83
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0940.0393.644.01178
      More than 1× weekly443.98
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.1083.774.19197
      Written milking protocolPresent1194.08
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.092<0.0013.894.26211
      Absent323.73
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.1123.513.95166
      Dry teats before attaching milking unitYes1383.73
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0810.0023.573.89166
      No134.08
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.1333.824.34211
      Type of grooming of beddingRaked943.890.0950.7403.704.08185
      Mechanical573.910.1063.714.12188
      Depth of beddingDeep bedded1253.990.0900.0653.814.16199
      Shallow bedded263.820.1173.594.05177
      a–c Estimates within a predictor with different superscripts differ. Tukey-adjusted P < 0.05.
      1 Inorganic bedding includes fresh sand (n = 120), recycled sand (n = 20), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 1), and field-grade lime (n = 1).
      2 Back transformation to SCC not used for analysis.
      3 Spring (March – May), summer (June – August), fall (September – November), and winter (December – February).

      Nonmanure Organic Bedding

      Data from 48 farms using OB was used to determine associations of eligible management factors with BMSCS. One of 49 farms originally eligible was excluded from analysis because complete BMSCS results were not available. Of explanatory variables eligible for entry into the BMSCS model, percent of cows with nonfunctioning mammary quarters, season, presence of a manger during milking, and regular removal and replacement of bedding in the back one-third of stalls remained in the final multivariate model. No interactions remained in the final model. Bulk milk SCS was least in the winter and spring, intermediate in the fall, and greatest in the summer (P < 0.001). The BMSCC for farms using OB was 0.37 log2 units greater when managers were present during all milking shifts and 0.33 log2 units greater when the back one-third of stalls were not cleaned out at least once weekly (Table 10). Increasing the percentage of cows with nonfunctioning mammary quarters by 1% tended to be associated with increased BMSCC (P = 0.069). All model assumptions were satisfied and model fit was satisfactory (AIC = 7.0).
      Table 10Adjusted estimates of mean bulk milk somatic cell score for all explanatory variables that remained in the final multivariate mixed repeated measures regression model for Wisconsin dairy farms (n = 48) producing ≥11,340 kg of milk daily and using nonmanure organic
      Nonmanure organic bedding includes sawdust, shavings, and grain hulls (n=33); straw (n=7); mixtures of organic beddings (n=4); straw or sawdust over lime or sand (n=2); and recycled paper products (n=2).
      bedding in all pens during the study period (linear score)
      Explanatory variableLevelsFarms, nEstimateSEType III P-value95% CISCC
      Back transformation to SCC not used for analysis.
      LowerUpper(1,000 cells/mL)
      Intercept3.500.117141
      Cows with nonfunctioning mammary quarters, %480.0200.0110.069−0.00150.0406
      Year (random)1480.00
      2480.00
      Season
      Spring (March – May), summer (June – August), fall (September – November), and winter (December – February).
      Spring483.95
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.064<0.0013.824.07193
      Summer484.20
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0654.074.33230
      Fall484.09
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0653.964.21213
      Winter483.95
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0643.824.07193
      Manager presence during milkingNever or when problems suspected323.86
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0720.0033.724.00182
      Always or portions of all shifts164.23
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.1014.034.43235
      Replace bedding in backs of stallsRegularly
      Farmers replaced bedding in the back one-third of stalls at least once weekly.
      173.88
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.1000.0083.694.08184
      Not regularly314.21
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0734.074.35231
      a–c Estimates within a predictor with different superscripts differ. Tukey-adjusted P < 0.05.
      1 Nonmanure organic bedding includes sawdust, shavings, and grain hulls (n = 33); straw (n = 7); mixtures of organic beddings (n = 4); straw or sawdust over lime or sand (n = 2); and recycled paper products (n = 2).
      2 Back transformation to SCC not used for analysis.
      3 Spring (March – May), summer (June – August), fall (September – November), and winter (December – February).
      4 Farmers replaced bedding in the back one-third of stalls at least once weekly.

      Manure-Based Organic Bedding

      Data from 23 farms using MB was used to determine associations of eligible management factors with BMSCS. Two of 25 farms eligible for screening models were excluded from analysis because BMSCS results were not available in each period. Of explanatory variables eligible for entry into the BMSCS model, the percentage of cows with milk discarded, season, and presence of a WMP remained in the final multivariate model. No interactions remained in the final model (Table 11). Bulk milk SCS was greater in the summer than other seasons (P = 0.014) and BMSCS was 0.42 log2 units greater with the presence of a WMP (P = 0.021). Each 1% increase in the percent of cows with milk discarded was associated with a 0.14 increase in BMSCS (P = 0.014). All model assumptions were satisfied, and model fit was satisfactory (AIC = 30.1).
      Table 11Adjusted estimates of mean bulk milk somatic cell score for all explanatory variables that remained in the final multivariate mixed repeated measures regression model for Wisconsin dairy farms (n = 23) producing ≥13,400 kg of milk daily and using manure-based
      Manure based bedding includes digester solids (n=10), composted manure (n=4), biosolids (n=3), composted digester solids (n=2), manure solids (n=3), and separated solids (n=1).
      bedding in all pens during the study period (linear score)
      Explanatory variableLevelsFarms, nEstimateSEType III P-value95% CISCC
      Back transformation to SCC not used for analysis.
      LowerUpper(1,000 cells/mL)
      Intercept3.970.188196
      Cows with milk discarded, %230.140.0550.0140.0270.244
      Year (random)1230.0
      2230.0
      Season
      Spring (March – May), summer (June – August), fall (September – November), and winter (December – February).
      Spring234.09
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.089<0.0013.914.26213
      Summer234.25
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0914.074.43238
      Fall234.15
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0903.984.33222
      Winter234.11
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.0893.944.29217
      Written milking protocol (WMP)Present114.36
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.1000.0213.654.23257
      Absent123.94
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      0.1463.654.23192
      a,b Estimates within a predictor with different superscripts differ. Tukey-adjusted P < 0.05.
      1 Manure based bedding includes digester solids (n = 10), composted manure (n = 4), biosolids (n = 3), composted digester solids (n = 2), manure solids (n = 3), and separated solids (n = 1).
      2 Back transformation to SCC not used for analysis.
      3 Spring (March – May), summer (June – August), fall (September – November), and winter (December – February).

      Association of Bedding Type with Bulk Milk Total Bacteria Count

      Processors reported BMTBC results for 230 farms and 72,978 individual BMTBC test results were received. No farms were out of compliance with PMO or EU standards for BMTBC during the 2-yr period. The distribution of total bacterial count values was ≤100,000 cfu/mL (98%), ≤ 25,000 cfu/mL (93%), and ≤10,000 cfu/mL (85%). Bulk tank total bacterial count category was not associated with bedding type (P = 0.787), bedding management practices (P ≥ 0.356), or herd management practices (P ≥ 0.113) included in the survey.

      Association of Bedding Type with Rolling Herd Average

      Data from 227 farms was used to determine associations of bedding type with RHA. Of explanatory variables eligible for entry into the RHA model (Table 12), the proportion of cows with nonfunctioning mammary quarters, bedding type, understanding of subclinical mastitis, and number of times the cows were milked each day remained in the final multivariate model. No interactions remained in the final model (Table 13). As compared with herds using IB, RHA were 771 and 1,153 kg less for herds using OB or MB, respectively (P < 0.001). The RHA was 353 kg greater for farms where the farmer understood subclinical mastitis as compared with farms where farmers had a limited understanding of subclinical mastitis (P = 0.017). Each 1% increase in cows with nonfunctioning mammary quarters was associated with a 57-kg decrease in RHA (P = 0.001).
      Table 12Unadjusted estimates and mean rolling herd average
      Rolling herd average from DHI records.
      for all explanatory variables that were eligible for inclusion in multivariate models for Wisconsin dairy farms (n = 227) producing ≥11,340 kg of milk daily (kg)
      Explanatory variableLevelsnEstimateSEP-value
      Cows with nonfunctioning mammary quarters, %227−76.0617.99<0.001
      BeddingIB
      Inorganic bedding includes fresh sand (n=122), recycled sand (n=26), mixture of sands (n=6), washed sand (n=3), fresh sand and lime (n=1), and field-grade lime (n=1).
      15912,87086.68<0.001
      OB
      Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n=30); straw (n=6); mixtures of organic beddings (n=5); straw or sawdust over lime or sand (n=2); and recycled paper products (n=1).
      4412,0251,644.770.031
      MB
      Manure based bedding includes digester solids (n=14), biosolids (n=2), composted manure (n=3), manure solids (n=2), composted digester solids (n=2), and separated solids (n=1).
      2411,779223.1
      California Mastitis Test used on the farmYes17112,298155.190.022
      No5612,68688.81
      Frequency of forestrippingAlways20112,65481.81
      Sometimes or never2612,096227.46
      Type of predip
      One farm did not use predip.
      Iodine14212,52395.870.149
      Noniodine8412,750124.65
      Type of postdip
      One farm did not use postdip.
      Iodine16012,46790.990.008
      Noniodine6512,921142.75
      Understanding of subclinical mastitisUnderstood15312,73993.270.006
      Limited knowledge7412,283134.11
      Milkings per day22211,907245.550.004
      320512,66480.44
      1 Rolling herd average from DHI records.
      2 Inorganic bedding includes fresh sand (n = 122), recycled sand (n = 26), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 1), and field-grade lime (n = 1).
      3 Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n = 30); straw (n = 6); mixtures of organic beddings (n = 5); straw or sawdust over lime or sand (n = 2); and recycled paper products (n = 1).
      4 Manure based bedding includes digester solids (n = 14), biosolids (n = 2), composted manure (n = 3), manure solids (n = 2), composted digester solids (n = 2), and separated solids (n = 1).
      5 One farm did not use predip.
      6 One farm did not use postdip.
      Table 13Adjusted estimates and mean rolling herd average
      Rolling herd average from DHI records.
      for all explanatory variables that remained in the final multivariate regression model for Wisconsin dairy farms (n = 220) producing ≥11,340 kg of milk daily and enrolled in DHI (kg)
      Explanatory variableLevelsnEstimateSEP-value95% CI
      LowerUpper
      Intercept22012,590186.5
      Cows with nonfunctioning mammary quarters, %220−56.716.30.001−88.9−24.6
      BeddingIB
      Inorganic bedding includes fresh sand (n=120), recycled sand (n=25), mixture of sands (n=6), washed sand (n=3), fresh sand and lime (n=1), and field-grade lime (n=1).
      15612,411
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      126.5<0.00112,16112,660
      OB
      Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n=29); straw (n=6); mixtures of organic beddings (n=5); straw or sawdust over lime or sand (n=2); and recycled paper products (n=1).
      4311,650
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      180.311,29512,006
      MB
      Manure-based bedding includes digester solids (n=11), biosolids (n=2), composted manure (n=1), manure solids (n=4), composted digester solids (n=2), and separated solids (n=1).
      2111,258
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      251.610,76311,754
      Understanding of subclinical mastitisUnderstood14911,950
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      147.70.01711,65912,241
      Limited knowledge7111,597
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      165.411,27111,923
      Milkings per day321912,256
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      99.3<0.00112,06012,452
      22111,291
      Estimates within a predictor with different superscripts differ. Tukey-adjusted P<0.05.
      237.310,82311,795
      a,b Estimates within a predictor with different superscripts differ. Tukey-adjusted P < 0.05.
      1 Rolling herd average from DHI records.
      2 Inorganic bedding includes fresh sand (n = 120), recycled sand (n = 25), mixture of sands (n = 6), washed sand (n = 3), fresh sand and lime (n = 1), and field-grade lime (n = 1).
      3 Organic nonmanure bedding includes sawdust, shavings, and grain hulls (n = 29); straw (n = 6); mixtures of organic beddings (n = 5); straw or sawdust over lime or sand (n = 2); and recycled paper products (n = 1).
      4 Manure-based bedding includes digester solids (n = 11), biosolids (n = 2), composted manure (n = 1), manure solids (n = 4), composted digester solids (n = 2), and separated solids (n = 1).

      Discussion

      At the time that the study was conducted, the herds that participated in this study contained about 22.2% of the Wisconsin dairy herd (1,267,000 cows;

      USDA-NASS. 2014a. Milk Production. ISSN: 1949–1557 Released Feb. 20, 2014. Accessed Mar. 17, 2015. https://usda.mannlib.cornell.edu/MannUsda/viewDocumentInfo.do?documentID=1103

      ) and represented a considerable proportion of milk produced in Wisconsin. In contrast to previous research using large Wisconsin farms (
      • Pantoja J.C.F.
      • Reinemann D.J.
      • Ruegg P.L.
      Associations among milk quality indicators in raw bulk milk.
      ;
      • Oliveira L.
      • Hulland C.
      • Ruegg P.L.
      Characterization of clinical mastitis occurring in cows on 50 large dairy herds in Wisconsin.
      ), eligibility criteria for this study was based simply on the amount of milk produced per farm per day. The single eligibility criteria required farms to sell at least 11,340 kg of milk per day. This criterion was selected as it represented approximately one-half of the capacity of a milk tanker truck. Farms of this size are producing an increasing proportion of the total US milk supply. Dairy herds in the United States containing ≥200 cows produced 76% of US milk in 2012 (

      USDA-NASS. 2014b. Quick Stats Tools. Accessed Apr. 29, 2014. http://www.nass.usda.gov/Quick_Stats/

      ). These farms use more intensive management strategies than smaller herds, and a better understanding of how bedding management on these farms influences productivity and milk quality is needed.
      The high rate of participation (84%) for farms meeting the eligibility criteria reduced selection bias and likely extends the reference population to comparably sized farms that use similar management practices in other regions. Several characteristics of study farms indicate that this group of larger herds used more intensive and technologically advanced management practices as compared with the general US dairy herd population. In contrast to the study herds, only 22 and 4% of US dairy farms are estimated to use IB and MB, respectively (

      USDA. 2008a. Dairy 2007, Facility Characteristics and Cow Comfort on U.S. Dairy Operations, 2007 USDA–APHIS–VS, CEAH. Fort Collins, CO #524.1210.

      ). Of study farms, only 8% milked cows twice daily in contrast to approximately 59% of US dairy herds containing ≥500 milking cows (

      USDA. 2008b. Dairy 2007, Part III: Reference of Dairy Cattle Health and Management Practices in the United States, 2007 USDA–APHIS–VS, CEAH. Fort Collins, CO #N482.0908.

      ). Standardized milking practices are highly adopted by Wisconsin dairy farmers of all herd sizes, indicating wide acceptance of these practices (
      • Rodrigues A.C.O.
      • Caraviello D.Z.
      • Ruegg P.L.
      Management of Wisconsin dairy herds enrolled in milk quality teams.
      ). Recommended best management practices for milking were highly adopted in study herds and were similar or greater than adoption by US dairy farms containing >500 milking cows (

      USDA. 2008b. Dairy 2007, Part III: Reference of Dairy Cattle Health and Management Practices in the United States, 2007 USDA–APHIS–VS, CEAH. Fort Collins, CO #N482.0908.

      ). Similar proportions of herds used forestripping (86.5% of study herds versus 83.5% of US herds milking ≥500 cows). However, a greater proportion of study herds adopted both pre- (99.1%) and postdipping (99.7%) as compared with pre- (80.2%) and postdipping (about 95%) performed on US herds milking ≥500 cows (

      USDA. 2008b. Dairy 2007, Part III: Reference of Dairy Cattle Health and Management Practices in the United States, 2007 USDA–APHIS–VS, CEAH. Fort Collins, CO #N482.0908.

      ).
      Dairy farms need to produce high quality milk to maintain market access (
      • Rowbotham R.F.
      The value of quality milk to the processor.
      ;
      • Norman H.D.
      • Lombard J.E.
      • Wright J.R.
      • Kopral C.A.
      • Rodriguez J.M.
      • Miller R.H.
      Consequence of alternative standards for bulk tank somatic cell count of dairy herds in the United States.
      ;
      • McCarty K.K.
      What processors really need for milk quality—And why.
      ). In the United States, raw milk quality is regulated by the Pasteurized Milk Ordinance (PMO;

      US-FDA. 2011. Grade “A” Pasteurized Milk Ordinance, 2011 revision. Accessed Mar. 13, 2014. http://www.fda.gov/downloads/Food/GuidanceRegulation/UCM291757.pdf

      ). The PMO establishes limits of 750,000 cells/mL and 100,000 cfu/mL for BMSCC and BMTBC, respectively.
      Milk quality standards of many dairy processors in Wisconsin are generally much more stringent than regulatory limits set in the PMO because minimizing BMSCC and BMTBC dramatically affects the final quality of dairy products (
      • Boor K.J.
      • Murphy S.C.
      Microbiology of market milks.
      ;
      • Barbano D.M.
      • Ma Y.
      • Santos M.V.
      Influence of raw milk quality on fluid milk shelf life.
      ) and allows US processors to gain access to export markets, which accounted for 15.4% of all sales of domestic production in 2013–2014 (

      US Dairy Export Council. 2015. NDM/SMP and whey products pace May exports. Accessed Apr. 14, 2015. www.usdec.org/Documents/Supplier%20Site/ResearchandData/USExportData/SummaryofU.S.ExportsYTD.pdf

      ). Although some Wisconsin milk manufacturers do not export their core products such as cheese and butter, most of the products made from the byproduct streams of manufacturing processes (such as whey) are exported or comingled with product which may be exported. To maintain the ability to export byproducts, individual farms must maintain a USDA export certificate that mandates compliance with EU milk quality regulations. European Union milk quality guidelines require that farms maintain a geometric mean SCC <400,000 cells/mL for 3 consecutive months and a geometric mean total bacterial count <100,000 cfu/mL for 2 consecutive months (

      Europa. 2009. Regulation (EC) No. 853/2004 of the European Parliament and of the Council of 29 April 2004 laying down specific hygiene rules for food of animal origin (OJ L 226, 25.6.2004, p. 22) [with successive amendments and corrections consolidated]. Accessed Mar. 21, 2015. http://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=CONSLEG:2004R0853:20071114:EN:PDF

      ). Thus, the majority of milk produced in Wisconsin is required to comply with EU standards, effectively making the FDA BMSCC and BMTBC requirements in the PMO superfluous for many farms.
      In this study, 100% of farms complied with US milk quality regulations for the 2-yr study period. Whereas 100% of farms complied with EU BMTBC regulations, 3 farms (2 using MB and 1 using OB) were potentially out of compliance with EU BMSCC standards. The 3 farms that were potentially out of compliance could have requested derogations from USDA for up to 1 yr while they worked to improve BMSCC. The 1.2% of study herds that were potentially out of compliance with EU standards for BMSCC is far less than previous estimates (7.8 to 16.1%;
      • Norman H.D.
      • Lombard J.E.
      • Wright J.R.
      • Kopral C.A.
      • Rodriguez J.M.
      • Miller R.H.
      Consequence of alternative standards for bulk tank somatic cell count of dairy herds in the United States.
      ). Better milk quality and a greater level of compliance with EU milk quality standards by study herds is likely due partly to greater adoption of best management practices by larger herds (
      • Wells S.J.
      • Ott S.L.
      What is the quality of the U.S. milk supply?.
      ;
      • Jayarao B.M.
      • Pillai S.R.
      • Sawant A.A.
      • Wolfgang D.R.
      • Hegde N.V.
      Guidelines for monitoring bulk tank milk somatic cell and bacterial counts.
      ;
      • Norman H.D.
      • Lombard J.E.
      • Wright J.R.
      • Kopral C.A.
      • Rodriguez J.M.
      • Miller R.H.
      Consequence of alternative standards for bulk tank somatic cell count of dairy herds in the United States.
      ) and partly in response to pressure from milk processors as they prepared to ensure compliance with EU standards that began on January 1, 2012.
      The BMSCC is well known to be greatest in summer and least during winter and spring primarily due to environmental conditions that encourage growth of mastitis pathogens during warmers seasons (
      • Ruegg P.L.
      • Tabone T.J.
      The relationship between antibiotic residue violations and somatic cell counts in Wisconsin dairy herds.
      ;
      • Olde Riekerink R.G.M.
      • Barkema H.W.
      • Stryhn H.
      The effect of season on somatic cell count and the incidence of clinical mastitis.
      ;
      • Archer S.C.
      • Mc Coy F.
      • Wapenaar W.
      • Green M.J.
      Association of season and herd size with somatic cell count for cows in Irish, English, and Welsh dairy herds.
      ). As expected, seasonality was significant in all 3 BMSCS models, indicating that seasonality in BMSCC is not dependent upon the type of bedding used within the study population.
      Drying of teats before attaching milking units is a recommended best management practice and was highly adopted by herds that participated in this study (
      • Moxley J.E.
      • Kennedy B.W.
      • Downey B.R.
      • Bowman J.S.T.
      Survey of milking hygiene practices and their relationships to somatic cell counts and milk production.
      ;
      • Galton D.M.
      • Petersson L.G.
      • Merrill W.G.
      Effects of premilking udder preparation practices on bacterial counts in milk and on teats.
      ;
      • Johnson A.P.
      A proper milking routine: The key to quality milk.
      ). For study herds that used IB, drying teats before attaching milking units was associated with a 0.35 log2 units lower BMSCC. Although it is likely that a similar association existed on herds using other bedding types, the high adoption of this practice in the smaller number of herds using other types of bedding did not allow for testing this association. Only 1 farm using OB and 2 farms using MB reported that they did not dry teats before attaching milking units. Use of teat scrubbers and not drying teats were correlated (r = 0.96) on farms with IB where all farms (n = 13) not drying teats also used teat scrubbers and only 1 farm using teat scrubbers dried teats. The current study was not able to model both of these correlated explanatory variables. Further research focusing on the effectiveness of teat scrubbers with or without teat drying is warranted.
      For herds using both IB and MB, BMSCS was greater for farms that had a WMP. Greater BMSCS on farms with written protocols, which outline work procedures and define best practices, is counterintuitive. Other research has also found contradictory associations with this variable. A previous study of 180 Wisconsin herds enrolled in milk quality management teams (
      • Rodrigues A.C.O.
      • Caraviello D.Z.
      • Ruegg P.L.
      Management of Wisconsin dairy herds enrolled in milk quality teams.
      ) reported that WMP were present in a greater proportion of herds with low (<250,000 cells/mL) or high (>400,000 cells/mL) BMSCC as compared with herds with medium BMSCC (250,000 to 400,000 cells/mL). However, no herds enrolled in the current study had BMSCC >400,000 cells/mL. It is likely that the presence of a WMP is irrelevant for herds that do not adequately implement the protocol and assessment of the actual milking procedures would be a more useful variable to include in future models. It is also possible that this is a spurious finding as cross-sectional observational studies can identify the significance and strength of associations between variables but cannot define the direction of association nor assign causality (
      • Dohoo I.
      • Martin W.
      • Stryhn H.
      ). It is possible that farmers who experienced milk quality problems worked with consultants to develop a WMP as a mechanism to reduce mastitis and improve milk quality. The presence of the WMP, without proper implementation would not be expected to have had positive results and it is plausible that farms with increased BMSCC in the past were more likely to have a WMP.
      On farms using OB, increased BMSCS was unexpectedly associated with the presence of managers during milking. It is possible that on farms with well-trained milking technicians, management presence in the parlors on a frequent basis has a lesser effect on BMSCC as compared with farms that have poorer results. Our survey instrument did not contain questions that allowed us to further evaluate this relationship and additional research about the role of training, use of written protocols, and labor management is warranted.
      Bedding is a primary reservoir for growth of most environmental mastitis pathogens and serves as a point of exposure of teats, so it is important to understand the relationship between pathogen populations in bedding and on teat skin (
      • Hogan J.S.
      • Smith K.L.
      • Todhunter D.A.
      • Schoenberger P.S.
      Bacterial counts associated with recycled newspaper bedding.
      ;
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      ;
      • Proietto R.L.
      • Hinckley L.S.
      • Fox L.K.
      • Andrew S.M.
      Evaluation of a clay-based acidic bedding conditioner for dairy cattle bedding.
      ). In vitro studies of the ability of sterilized bedding materials to support pathogen growth have shown that recycled manure can support the log growth phase for 0 to 24 h, 6 to 24 h, or no log growth for coliforms, Streptococcus uberis and Enterococcus faecium, respectively (
      • Zehner M.M.
      • Farnsworth R.J.
      • Appleman R.D.
      • Larntz K.
      • Springer J.A.
      Growth of environmental mastitis pathogens in various bedding materials.
      ;
      • Godden S.
      • Bey R.
      • Lorch K.
      • Farnsworth R.
      • Rapnicki P.
      Ability of organic and inorganic bedding materials to promote growth of environmental bacteria.
      ). These studies reported log growth of coliforms for up to 48 h (followed by a stationary phase for up to 102 h) in clean sand and wood products, but no log or stationary growth phases in these same bedding materials for streptococci (which began a death phase at the time of inoculation). These growth patterns demonstrate the ability of both organic and inorganic bedding to supply the necessary nutrients for pathogen growth for short periods of time.
      • Kristula M.A.
      • Dou Z.
      • Toth J.D.
      • Smith B.I.
      • Harvey N.
      • Sabo M.
      Evaluation of free-stall mattress bedding treatments to reduce mastitis bacterial growth.
      reported that the number of total coliforms, Klebsiella spp., streptococci, and Escherichia coli were least at the time of placement of wood shavings into stalls (with mattresses) and continued to increase (when measured at 12 h intervals) up to 48 h after placement. In a study of mastitis pathogens in sand and sawdust bedding (
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      ), growth in sawdust was similar to
      • Kristula M.A.
      • Dou Z.
      • Toth J.D.
      • Smith B.I.
      • Harvey N.
      • Sabo M.
      Evaluation of free-stall mattress bedding treatments to reduce mastitis bacterial growth.
      , with exponential growth of bacteria after placement through d 1. After the first day, bacterial growth continued at a slower rate through the rest of the trial periods. Our study found an association between decreased BMSCS and removal of OB from the back one-third of stalls at least once weekly. This relationship may have resulted from reduced exposure to mastitis pathogens in clean bedding. For farms using IB and OB, associations of regular removal and replacement of bedding and BMSCS could not be evaluated due to small sample sizes. Only 5 farms using IB replaced bedding at least once monthly and only 1 farm using MB replaced bedding more than 4 times per year.
      • Kristula M.A.
      • Rogers W.
      • Hogan J.S.
      • Sabo M.
      Comparison of bacteria populations in clean and recycled sand used for bedding in dairy facilities.
      measured potential mastitis pathogens in clean and recycled sand bedding and reported that quantities were least at time of placement of bedding into stalls and increased logarithmically until 1 d after placement.
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      also reported that the number of potential mastitis pathogens cultured from sand bedding peaked 1 d after placement of bedding in stalls. Interestingly, the number of Streptococcus spp. remained constant through the end of 1-wk trial periods (n = 3) whereas the number of coliforms and Klebsiella declined by the end of the trial periods. Thus maximum exposure to pathogens in IB occurs approximately 1 to 2 d after new bedding is placed in stalls. In our study, the addition of new IB to stalls more than once weekly was associated with increased BMSCS. It is possible that frequent addition of new bedding supplied nutrients that maintained organisms in a logarithmic growth phase, rather than allowing the pathogens to consume available nutrients and enter a stationary or death phase, as suggested by
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      . On farms using IB, increased frequency of adding new bedding was associated with increased BMSCS; however, on farms using OB, increased frequency of replacement was associated with decreased BMSCS. These seemingly contradictory associations are representative of the complexity of bedding management on large dairy farms. Of farms using OB, 35% removed and replaced the bedding in the backs of stalls at least weekly. In contrast, of farms using IB, only 3% removed and replaced bedding at least weekly. Most farms using IB simply added new bedding to bedding already present in stalls. Differences in bedding management practices of farms using different types of bedding (Table 8) and differences in BMSCS related to bedding management practices suggest that further research about these complex relationships is necessary.
      Several studies have shown strong positive correlations between the number and type of bacteria in bedding and on teat skin (
      • Hogan J.S.
      • Smith K.L.
      • Todhunter D.A.
      • Schoenberger P.S.
      Bacterial counts associated with recycled newspaper bedding.
      ,
      • Hogan J.S.
      • Bogacz V.L.
      • Thompson L.M.
      • Romig S.
      • Schoenberger P.S.
      • Weiss W.P.
      • Smith K.L.
      Bacterial counts associated with sawdust and recycled manure bedding treated with commercial conditioners.
      ;
      • Hogan J.S.
      • Smith K.L.
      Bacteria counts in sawdust bedding.
      ). Correlations between pathogen populations in bedding and teat skin have ranged from 0.36 (coliforms in sand bedding;
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      ) to 0.97 (streptococci in sawdust plus an acidic conditioner;
      • Proietto R.L.
      • Hinckley L.S.
      • Fox L.K.
      • Andrew S.M.
      Evaluation of a clay-based acidic bedding conditioner for dairy cattle bedding.
      ). The number and distribution of organisms found in bedding profoundly affects colonization of teat skin (
      • Rendos J.J.
      • Eberhart R.J.
      • Kesler E.M.
      Microbial populations of teat ends of dairy cows, and bedding materials.
      ;
      • Hogan J.S.
      • Smith K.L.
      Bacteria counts in sawdust bedding.
      ). In general, bacterial counts on teat skin have been reported to be greater for cows housed on organic bedding as compared with cows housed on inorganic bedding (
      • Fairchild T.P.
      • McArthur B.J.
      • Moore J.H.
      • Hylton W.E.
      Coliform counts in various bedding materials.
      ;
      • Janzen J.J.
      • Bishop R.R.
      • Bodine A.B.
      • Caldwell C.A.
      • Johnson D.W.
      Composted dairy waste solids and crushed limestone as bedding in free stalls.
      ;
      • Zdanowicz M.
      • Shelford J.A.
      • Tucker C.B.
      • Weary D.M.
      • von Keyserlingk M.A.G.
      Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
      ). Studies of the association between bacterial populations in bedding and on teat skin have usually been of short duration, with cows exposed to each bedding for 1 to 3 wk (
      • Hogan J.S.
      • Smith K.L.
      • Todhunter D.A.
      • Schoenberger P.S.
      Bacterial counts associated with recycled newspaper bedding.
      ,
      • Hogan J.S.
      • Bogacz V.L.
      • Thompson L.M.
      • Romig S.
      • Schoenberger P.S.
      • Weiss W.P.
      • Smith K.L.
      Bacterial counts associated with sawdust and recycled manure bedding treated with commercial conditioners.
      ;
      • Hogan J.S.
      • Smith K.L.
      Bacteria counts in sawdust bedding.
      ). Although these short duration studies are sufficient to determine associations between bacterial populations of bedding and on teat skin, the exposure period is not long enough to be able to identify associations with the risk of mastitis. Results of our study indicate that bedding type is associated with indicators of mastitis including BMSCS and the percentage of cows with nonfunctioning mammary quarters (Table 2). Further research is needed to improve bedding management and reduce the risk of mastitis for herds that are using organic (OB or MB) bedding.
      As compared with herds included in previous studies (
      • Borneman D.L.
      • Stiegert K.
      • Ingham S.
      Limitations of on-site dairy farm regulatory debits as milk quality predictors.
      ;
      • O’Connell A.
      • McParland S.
      • Ruegg P.L.
      • O’Brien B.
      • Gleeson D.
      Seasonal trends in milk quality in Ireland between 2007 and 2011.
      ), herds included in this study had exceptionally good BMTBC and were similar to an earlier study of large Wisconsin dairies (
      • Pantoja J.C.F.
      • Reinemann D.J.
      • Ruegg P.L.
      Associations among milk quality indicators in raw bulk milk.
      ).
      • Borneman D.L.
      • Stiegert K.
      • Ingham S.
      Limitations of on-site dairy farm regulatory debits as milk quality predictors.
      reported that the 75th percentile of BMTBC for 1,301 Wisconsin grade A farms subjected to regulatory survey in 2012 was 25,000 cells/mL. However, farms included in that study were a subset of all Wisconsin grade A dairy farms and were not selected based on farm size. Only 7% of BMTBC in our study were >25,000 cfu/mL and the technological advantages of management practices used by larger dairy herds on BMTBC is apparent. In contrast, in a study that included almost all Irish dairy herds (n = 10,819),
      • O’Connell A.
      • McParland S.
      • Ruegg P.L.
      • O’Brien B.
      • Gleeson D.
      Seasonal trends in milk quality in Ireland between 2007 and 2011.
      reported that the 75th percentile for BMTBC herds ranged from 25,000 to 30,000 cfu/mL. However, the average Irish dairy herd contained only 55 milking cows, often stored milk in bulk tanks on farms for several days and used less rigorous methods of washing milking equipment. The quality of milk produced on Wisconsin and US farms may continue to improve as processor standards become more rigorous to meet consumer and export demands and milk production continues to shift from smaller to larger farms which tend to adopt newer technologies sooner.
      The type of bedding used on farms was not associated with BMTBC. Management practices that result in production of high quality milk with minimal bacterial contamination are well documented and include proper udder hygiene and pre-milking sanitization, postmilking sanitization, effective milking machine cleaning, adequate cooling, and storage of milk between 0 and 4.4°C (
      • Reinemann D.J.
      • Mein G.A.
      • Bray D.R.
      • Reid D
      • Britt J.S.
      Troubleshooting high bacteria counts in farm milk.
      ;
      • Reinemann D.J.
      Troubleshooting high bacteria counts in bulk milk: What needs cleaning, the machine or the cows?.
      ). Maintaining optimal BMTBC is generally less complicated than maintaining optimal BMSCC. Unlike increased BMSCC, resolution of problems with increased BMTBC is typically accomplished by correction of equipment failures (such as malfunctioning cleaning processes) or procedural lapses (such as poor premilking sanitization of teats;
      • Reinemann D.J.
      • Mein G.A.
      • Bray D.R.
      • Reid D
      • Britt J.S.
      Troubleshooting high bacteria counts in farm milk.
      ). Improvements in milking procedures, cleaning of equipment, and proper milk cooling can result in immediate and long-lasting improvements in BMTBC. In contrast to the long-term effect of mastitis on BMSCC (
      • Keefe G.
      Update on control of Staphylococcus aureus and Streptococcus agalactiae for management of mastitis.
      ), resolution of the underlying problem causing increased BMTBC has little or no residual effect. Improvements in equipment, cooling technologies, or hygiene often result in BMTBC values that meet or exceed processor thresholds for maximizing quality premiums.
      • Borneman D.L.
      • Stiegert K.
      • Ingham S.
      Limitations of on-site dairy farm regulatory debits as milk quality predictors.
      were unable to show associations of official BMTBC counts and results of regulatory inspection visits that are used to verify compliance with the PMO, mostly likely because most milk processors in Wisconsin enforce more rigorous standards for BMTBC than required by regulators and herds that exceed processor threshold for BMTBC are often promptly visited by processor representatives. Bulk milk total bacterial count values that almost universally meet or exceed market and regulatory expectations indicate that most of the large dairy farms included in this study have implemented management practices necessary to control BMTBC to the point that it is no longer a milk quality issue for most of the surveyed farms.
      The use of IB was associated with considerable increases in RHA, which would result in increased gross revenue from milk sales. The RHA of farms using IB were 761 and 1,153 kg greater than the RHA of farms using OB and MB, respectively. The average farm included in this study contained 853 milking cows and the average milk price from May 2010 to April 2012 for Wisconsin dairy farms was $18.52 per 46.36 kg of milk (

      USDA-AMS. 2015. Federal Order Mailbox Prices, Upper Midwest Milk Marketing Area. Accessed May 13, 2015. http://www.fmma30.com/Homepage/FO30-MailboxPrices.htm

      ). Using these assumptions, the difference in RHA for herds using IB equates to additional annual gross revenue per cow of $304 and $461 as compared with herds using OB or MB, respectively. Using the same assumptions, average gross revenue per cow for the 220 herds used in the RHA model was $5,074, and this potential difference represents an additional potential of 6 to 9% in average annual gross milk sales per cow or $392,894 per year for an average study farm. The large potential financial opportunity associated with use of IB warrants further research of the possible economic implications. Potential benefits associated with increased gross revenue need to be modeled against cost differences associated with management of different bedding materials.
      Results of this study demonstrate the potential effect of mastitis on herd productivity and illustrate potential negative consequences on animal health based on bedding type. Rolling herd average of cows on farms using IB was greater than on farms using organic bedding (OB or MB). The data included in this study suggest that at least part of the differences in productivity associated with different bedding types may be attributed to better udder health on farms using IB as contrasted with farms using OB or MB. Cows free of mastitis produce more milk than infected cows (
      • Seegers H.
      • Fourichon C.
      • Beaudeau F.
      Production effects related to mastitis and mastitis economics in dairy cattle herds.
      ;
      • Bar D.
      • Tauer L.W.
      • Bennett G.
      • González R.N.
      • Hertl J.A.
      • Schukken Y.H.
      • Schulte H.F.
      • Welcome F.L.
      • Gröhn Y.T.
      The cost of generic clinical mastitis in dairy cows as estimated by using dynamic programming.
      ). Study herds bedded with IB had a lesser BMSCS, indicating that the prevalence of subclinical mastitis was likely greater in herds using organic bedding (
      • Schukken Y.H.
      • Wilson D.J.
      • Welcome F.
      • Garrison-Tikofsky L.
      • Gonzalez R.N.
      Monitoring udder health and milk quality using somatic cell counts.
      ). Farms using MB had a greater percentage of cows with milk discarded, likely indicating a greater incidence of clinical mastitis, which necessitates discarding milk. Greater incidence of clinical mastitis on farms using MB contributes to less efficient milk production and may ultimately lead to farmers drying off individual mammary gland quarters of chronically affected cows (as evidenced by the tendency for a greater percentage of cows with nonfunctioning mammary quarters on farms using MB as compared with IB).
      Possible causes of nonfunctioning mammary quarters include prepartum suckling, termination of production due to chronic mastitis, trauma, or intentional cessation of lactation due to chronic mastitis (
      • Duraes M.C.
      • Wilcox C.J.
      • Head H.H.
      • Van Horn H.H.
      Frequency and effects on production of blind quarters in first lactation dairy cows.
      ). Except for a single study that focused on quarters which were nonfunctional at first parturition (
      • Duraes M.C.
      • Wilcox C.J.
      • Head H.H.
      • Van Horn H.H.
      Frequency and effects on production of blind quarters in first lactation dairy cows.
      ), studies of the effects of nonfunctional quarters on milk production are lacking. In the present study, at the herd level, decreasing the prevalence of cows with nonfunctional mammary quarters by 1% was associated with a 57-kg increase in RHA. The authors hypothesize that when nutritional management is optimized, the total number of functional alveoli in the udder may be a limiting factor for milk production, thus resulting in reduced production from animals with at least one nonfunctioning mammary quarter. Using the results from the multivariate RHA model, a reduction in the proportion of cows milking with fewer than 4 functioning mammary quarters by 1% would result in approximately $23.00 in increased milk sales per cow annually. Approximately 5% of study herds reported >10% of milking cows with at least one nonfunctioning mammary quarter. For these herds, reducing this proportion to the study median of 4% would result in at least $136 in additional annual milk sales per cow ($115,926 per herd) for an average farm. Reducing the proportion of cows with at least one nonfunctioning mammary quarter could be accomplished through improved mastitis prevention strategies.
      Due to the large scale of the farms that enrolled in this study, even modest improvements in productivity result in large changes in gross farm income. These improvements may justify investment in capital or labor to implement new processes or procedures (such as changing the type of bedding that is used). Although a complete financial analysis of the implications of changing the type of bedding is beyond the scope of this study, potential returns from implementing some of the practices associated with increased RHA or decreased BMSCS could result in considerably increased gross milk value and warrant consideration by farm owners.

      Conclusions

      Large Wisconsin dairy farms produce high-quality milk with minimal BMTBC and low BMSCC in part due to almost universal adoption of best management practices. Farms using IB (as compared with use of organic bedding) had greater RHA and better milk quality. Several management practices were associated with improved BMSCS. Farms using IB had lesser BMSCS when bedding was added at >7-d intervals and lesser BMSCS when teats were dried before attaching the milking unit. Production of higher quality milk on farms using IB had the potential to increase the value of each unit of milk produced. Further research is warranted into the relative economic value of using different bedding materials and the effect of bedding choice on animal health and well-being.

      Acknowledgments

      The authors thank dairy owners and milk buyers for their participation in this project and the Grande Cheese Company (Brownsville, WI) for project funding.

      References

        • Akaike H.
        Fitting autoregressive models for prediction.
        Ann. Inst. Stat. Math. 1969; 21: 243-247
        • Archer S.C.
        • Mc Coy F.
        • Wapenaar W.
        • Green M.J.
        Association of season and herd size with somatic cell count for cows in Irish, English, and Welsh dairy herds.
        Vet. J. 2013; 196: 515-521
        • Bar D.
        • Tauer L.W.
        • Bennett G.
        • González R.N.
        • Hertl J.A.
        • Schukken Y.H.
        • Schulte H.F.
        • Welcome F.L.
        • Gröhn Y.T.
        The cost of generic clinical mastitis in dairy cows as estimated by using dynamic programming.
        J. Dairy Sci. 2008; 91: 2205-2214
        • Barbano D.M.
        • Ma Y.
        • Santos M.V.
        Influence of raw milk quality on fluid milk shelf life.
        J. Dairy Sci. 2006; 89: E15-E19
        • Boor K.J.
        • Murphy S.C.
        Microbiology of market milks.
        in: Robinson R.K. Dairy Microbiology Handbook. 3rd. John Wiley & Sons Inc., Hoboken, NJ2002: 91-122
        • Borneman D.L.
        • Stiegert K.
        • Ingham S.
        Limitations of on-site dairy farm regulatory debits as milk quality predictors.
        J. Dairy Sci. 2015; 98: 1526-1538
        • Cook N.B.
        • Bennett T.B.
        • Nordlund K.V.
        Monitoring indices of cow comfort in free-stall-housed dairy herds.
        J. Dairy Sci. 2005; 88: 3876-3885
        • Costello M.
        • Rhee M.
        • Bates M.P.
        • Clark S.
        • Luedecke L.O.
        • Kang D.
        Eleven-year trends of microbiological quality in bulk tank milk.
        Food Prot. Trends. 2003; 23: 393-400
        • Deluyker H.A.
        • Gay J.M.
        • Weaver L.D.
        Interrelationships of somatic cell count, mastitis, and milk yield in a low somatic cell count herd.
        J. Dairy Sci. 1993; 76: 3445-3452
        • DeVries T.J.
        • Aarnoudse M.G.
        • Barkema H.W.
        • Leslie K.E.
        • von Keyserlingk M.A.G.
        Associations of dairy cow behavior, barn hygiene, cow hygiene, and risk of elevated somatic cell count.
        J. Dairy Sci. 2012; 95: 5730-5739
        • Dohoo I.
        • Martin W.
        • Stryhn H.
        Methods in Epidemiologic Research. Ver. Inc., Charlottetown, Prince Edward Island, Canada2012
        • Duraes M.C.
        • Wilcox C.J.
        • Head H.H.
        • Van Horn H.H.
        Frequency and effects on production of blind quarters in first lactation dairy cows.
        J. Dairy Sci. 1982; 65: 1804-1807
      1. Europa. 2009. Regulation (EC) No. 853/2004 of the European Parliament and of the Council of 29 April 2004 laying down specific hygiene rules for food of animal origin (OJ L 226, 25.6.2004, p. 22) [with successive amendments and corrections consolidated]. Accessed Mar. 21, 2015. http://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=CONSLEG:2004R0853:20071114:EN:PDF

        • Fairchild T.P.
        • McArthur B.J.
        • Moore J.H.
        • Hylton W.E.
        Coliform counts in various bedding materials.
        J. Dairy Sci. 1982; 65: 1029-1035
        • Fulwider W.K.
        • Grandin T.
        • Rollin B.E.
        • Engle T.E.
        • Dalsted N.L.
        • Lamm W.D.
        Survey of dairy management practices on one hundred thirteen north central and northeastern United States dairies.
        J. Dairy Sci. 2008; 91: 1686-1692
        • Galton D.M.
        • Petersson L.G.
        • Merrill W.G.
        Effects of premilking udder preparation practices on bacterial counts in milk and on teats.
        J. Dairy Sci. 1986; 69: 260-266
        • Godden S.
        • Bey R.
        • Lorch K.
        • Farnsworth R.
        • Rapnicki P.
        Ability of organic and inorganic bedding materials to promote growth of environmental bacteria.
        J. Dairy Sci. 2008; 91: 151-159
        • Hogan J.
        • Smith K.L.
        Managing environmental mastitis.
        Vet. Clin. North Am. Food Anim. Pract. 2012; 28: 217-224
        • Hogan J.S.
        • Bogacz V.L.
        • Thompson L.M.
        • Romig S.
        • Schoenberger P.S.
        • Weiss W.P.
        • Smith K.L.
        Bacterial counts associated with sawdust and recycled manure bedding treated with commercial conditioners.
        J. Dairy Sci. 1999; 82: 1690-1695
        • Hogan J.S.
        • Smith K.L.
        Bacteria counts in sawdust bedding.
        J. Dairy Sci. 1997; 80: 1600-1605
        • Hogan J.S.
        • Smith K.L.
        • Hoblet K.H.
        • Todhunter D.A.
        • Schoenberger P.S.
        • Hueston W.D.
        • Pritchard D.E.
        • Bowman G.L.
        • Heider L.E.
        • Brockett B.L.
        • Conrad H.R.
        Bacterial counts in bedding materials used on nine commercial dairies.
        J. Dairy Sci. 1989; 72: 250-258
        • Hogan J.S.
        • Smith K.L.
        • Todhunter D.A.
        • Schoenberger P.S.
        Bacterial counts associated with recycled newspaper bedding.
        J. Dairy Sci. 1990; 73: 1756-1761
        • Janzen J.J.
        • Bishop R.R.
        • Bodine A.B.
        • Caldwell C.A.
        • Johnson D.W.
        Composted dairy waste solids and crushed limestone as bedding in free stalls.
        J. Dairy Sci. 1982; 65: 1025-1028
        • Jayarao B.M.
        • Pillai S.R.
        • Sawant A.A.
        • Wolfgang D.R.
        • Hegde N.V.
        Guidelines for monitoring bulk tank milk somatic cell and bacterial counts.
        J. Dairy Sci. 2004; 87: 3561-3573
        • Johnson A.P.
        A proper milking routine: The key to quality milk.
        in: Proc. Natl. Mastitis Council Mtg., Reno, NV Natl. Mastitis Counc. Inc., Verona, WI2000: 123
        • Keefe G.
        Update on control of Staphylococcus aureus and Streptococcus agalactiae for management of mastitis.
        Vet. Clin. Food Anim. 2012; 28: 203-216
        • Kristula M.A.
        • Dou Z.
        • Toth J.D.
        • Smith B.I.
        • Harvey N.
        • Sabo M.
        Evaluation of free-stall mattress bedding treatments to reduce mastitis bacterial growth.
        J. Dairy Sci. 2008; 91: 1885-1892
        • Kristula M.A.
        • Rogers W.
        • Hogan J.S.
        • Sabo M.
        Comparison of bacteria populations in clean and recycled sand used for bedding in dairy facilities.
        J. Dairy Sci. 2005; 88: 4317-4325
        • Lago A.
        • Godden S.M.
        • Bey R.
        • Ruegg P.L.
        • Leslie K.
        The selective treatment of clinical mastitis based on on-farm culture results: I. Effects on antibiotic use, milk withholding time, and short-term clinical and bacteriological outcomes.
        J. Dairy Sci. 2011; 94: 4441-4456
        • Makovec J.A.
        • Ruegg P.L.
        Results of milk samples submitted for microbiological examination in Wisconsin from 1994 to 2001.
        J. Dairy Sci. 2003; 86: 3466-3472
        • McCarty K.K.
        What processors really need for milk quality—And why.
        in: Natl. Mastitis Counc. Mtg. Proc., Memphis, TN Natl. Mastitis Counc. Inc., Verona, WI2015: 92-94
        • Moxley J.E.
        • Kennedy B.W.
        • Downey B.R.
        • Bowman J.S.T.
        Survey of milking hygiene practices and their relationships to somatic cell counts and milk production.
        J. Dairy Sci. 1978; 61: 1637-1644
        • Norman H.D.
        • Lombard J.E.
        • Wright J.R.
        • Kopral C.A.
        • Rodriguez J.M.
        • Miller R.H.
        Consequence of alternative standards for bulk tank somatic cell count of dairy herds in the United States.
        J. Dairy Sci. 2011; 94: 6243-6256
        • O’Connell A.
        • McParland S.
        • Ruegg P.L.
        • O’Brien B.
        • Gleeson D.
        Seasonal trends in milk quality in Ireland between 2007 and 2011.
        J. Dairy Sci. 2015; 98: 3778-3790
        • Olde Riekerink R.G.M.
        • Barkema H.W.
        • Stryhn H.
        The effect of season on somatic cell count and the incidence of clinical mastitis.
        J. Dairy Sci. 2007; 90: 1704-1715
        • Oliveira L.
        • Hulland C.
        • Ruegg P.L.
        Characterization of clinical mastitis occurring in cows on 50 large dairy herds in Wisconsin.
        J. Dairy Sci. 2013; 96: 7538-7549
        • Ott S.L.
        • Novak P.R.
        Association of herd productivity and bulk-tank somatic cell counts in US dairy herds in 1996.
        J. Am. Vet. Med. Assoc. 2001; 218: 1325-1330
        • Pantoja J.C.F.
        • Reinemann D.J.
        • Ruegg P.L.
        Associations among milk quality indicators in raw bulk milk.
        J. Dairy Sci. 2009; 92: 4978-4987
        • Pinzón-Sánchez C.
        • Ruegg P.L.
        Risk factors associated with short-term post- treatment outcomes of clinical mastitis.
        J. Dairy Sci. 2011; 94: 3397-3410
        • Proietto R.L.
        • Hinckley L.S.
        • Fox L.K.
        • Andrew S.M.
        Evaluation of a clay-based acidic bedding conditioner for dairy cattle bedding.
        J. Dairy Sci. 2013; 96: 1044-1053
        • Reinemann D.J.
        Troubleshooting high bacteria counts in bulk milk: What needs cleaning, the machine or the cows?.
        Countdown Meet, Melbourne, Australia2011
        • Reinemann D.J.
        • Mein G.A.
        • Bray D.R.
        • Reid D
        • Britt J.S.
        Troubleshooting high bacteria counts in farm milk.
        in: Natl. Mastitis Counc. Mtg. Proc., Albuquerque, NM Natl. Mastitis Counc. Inc., Verona, WI1997: 16-19
        • Rendos J.J.
        • Eberhart R.J.
        • Kesler E.M.
        Microbial populations of teat ends of dairy cows, and bedding materials.
        J. Dairy Sci. 1975; 58: 1492-1500
        • Rodrigues A.C.O.
        • Caraviello D.Z.
        • Ruegg P.L.
        Management of Wisconsin dairy herds enrolled in milk quality teams.
        J. Dairy Sci. 2005; 88: 2660-2671
        • Rollin B.E.
        Annual meeting keynote address: Animal agriculture and emerging social ethics for animals.
        J. Anim. Sci. 2004; 82: 955-964
        • Rowbotham R.F.
        The value of quality milk to the processor.
        in: Natl. Mastitis Counc. Mtg. Proc., Cleveland, OH Natl. Mastitis Counc. Inc., Verona, WI2000: 1-6
        • Ruegg P.L.
        • Tabone T.J.
        The relationship between antibiotic residue violations and somatic cell counts in Wisconsin dairy herds.
        J. Dairy Sci. 2000; 83: 2805-2809
        • Sapp S.G.
        • Arnot C.
        • Fallon J.
        • Fleck T.
        • Soorholtz D.
        • Sutton-Vermeulen M.
        • Wilson J.J.H.
        Consumer trust in the U.S. food system: An examination of the recreancy theorem.
        Rural Sociol. 2009; 74: 525-545
        • SAS Institute
        SAS/STAT User’s Guide 9.4: Statistics..
        SAS Institute Inc., Cary, NC2014
        • Schreiner D.A.
        • Ruegg P.L.
        Relationship between udder and leg hygiene scores and subclinical mastitis.
        J. Dairy Sci. 2003; 86: 3460-3465
        • Schukken Y.H.
        • Wilson D.J.
        • Welcome F.
        • Garrison-Tikofsky L.
        • Gonzalez R.N.
        Monitoring udder health and milk quality using somatic cell counts.
        Vet. Res. 2003; 34: 579-596
        • Seegers H.
        • Fourichon C.
        • Beaudeau F.
        Production effects related to mastitis and mastitis economics in dairy cattle herds.
        Vet. Res. 2003; 34: 475-491
        • Shook G.E.
        Genetic improvement of mastitis through selection on somatic cell count.
        Vet. Clin. North Am. Food Anim. Pract. 1993; 9: 563-581
        • Tucker C.B.
        • Weary D.M.
        Bedding on geotextile mattresses: How much is needed to improve cow comfort?.
        J. Dairy Sci. 2004; 87: 2889-2895
      2. US Dairy Export Council. 2015. NDM/SMP and whey products pace May exports. Accessed Apr. 14, 2015. www.usdec.org/Documents/Supplier%20Site/ResearchandData/USExportData/SummaryofU.S.ExportsYTD.pdf

      3. USDA. 2008a. Dairy 2007, Facility Characteristics and Cow Comfort on U.S. Dairy Operations, 2007 USDA–APHIS–VS, CEAH. Fort Collins, CO #524.1210.

      4. USDA. 2008b. Dairy 2007, Part III: Reference of Dairy Cattle Health and Management Practices in the United States, 2007 USDA–APHIS–VS, CEAH. Fort Collins, CO #N482.0908.

      5. USDA-AMS. 2015. Federal Order Mailbox Prices, Upper Midwest Milk Marketing Area. Accessed May 13, 2015. http://www.fmma30.com/Homepage/FO30-MailboxPrices.htm

      6. USDA-NASS. 2014a. Milk Production. ISSN: 1949–1557 Released Feb. 20, 2014. Accessed Mar. 17, 2015. https://usda.mannlib.cornell.edu/MannUsda/viewDocumentInfo.do?documentID=1103

      7. USDA-NASS. 2014b. Quick Stats Tools. Accessed Apr. 29, 2014. http://www.nass.usda.gov/Quick_Stats/

      8. US-FDA. 2011. Grade “A” Pasteurized Milk Ordinance, 2011 revision. Accessed Mar. 13, 2014. http://www.fda.gov/downloads/Food/GuidanceRegulation/UCM291757.pdf

        • Velthuis A.G.J.
        • van Asseldonk M.A.P.M.
        Process audits versus product quality monitoring of bulk milk.
        J. Dairy Sci. 2011; 94: 235-249
        • Wells S.J.
        • Ott S.L.
        What is the quality of the U.S. milk supply?.
        in: Proc. 37th Annu. Mtg. Natl. Mastitis Counc., St. Louis, MO Natl. Mastitis Counc., Verona, WI1998: 10-19
      9. Wisconsin Administrative Code. 2013. Chapter ATCP 60 Dairy Farms. Accessed May 13, 2014. http://www.legis.state.wi.us/rsb/code/atcp/atcp060.pdf

      10. Wisconsin Department of Administration. 2012. Bureau of Document Services Document Sales & Distribution, Madison, WI.

        • Zdanowicz M.
        • Shelford J.A.
        • Tucker C.B.
        • Weary D.M.
        • von Keyserlingk M.A.G.
        Bacterial populations on teat ends of dairy cows housed in free stalls and bedded with either sand or sawdust.
        J. Dairy Sci. 2004; 87: 1694-1701
        • Zehner M.M.
        • Farnsworth R.J.
        • Appleman R.D.
        • Larntz K.
        • Springer J.A.
        Growth of environmental mastitis pathogens in various bedding materials.
        J. Dairy Sci. 1986; 69: 1932-1941