An Econometric Analysis of Climatic Effects on Total Factor Productivity Across U.S. Dairy Counties
Kamil Bora Bolat, Merve Bolat, Boris E. Bravo-Ureta

TL;DR
This study examines how heat stress from climate change affects dairy productivity in U.S. counties, finding that technological progress largely offsets these effects.
Contribution
The study introduces a detailed statistical model to separate climatic effects from technological progress in dairy productivity.
Findings
Technological advancements drive 2.52% annual total factor productivity growth in U.S. dairy counties.
Heat stress reduces productivity growth by an average of 0.008% annually, with the largest reduction at 0.08%.
Adaptation strategies like cooling systems have effectively mitigated heat stress impacts on productivity.
Abstract
High-yielding dairy cows face increasing heat stress due to climate change. Our study analyzed over four decades of data from U.S. dairy counties to understand the economic implications of this challenge. Using a detailed statistical model, we separated the effects of technological progress from climatic pressures on total factor productivity. Our findings show that while technology is the primary driver of productivity growth, heat stress imposes a modest constraint on this growth. Although nationwide results indicate a low reduction in the rate of productivity growth from heat, our analysis reveals that this impact varies by region. More importantly, the slight overall effect suggests that the U.S. dairy industry has been highly effective in managing these climatic challenges through innovation. This finding implies the sector’s resilience and the success of locally implemented…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsEffects of Environmental Stressors on Livestock · Climate change impacts on agriculture · Agriculture Sustainability and Environmental Impact
1. Introduction
The agricultural, food, and related industries constitute the economic bedrock of many states in the U.S., where the dairy sector is a vital component. Contributing about 12% of the world’s milk output and generating annual revenues nearing USD 51 billion, the U.S. dairy sector is a cornerstone of many rural communities, supporting local services such as feed suppliers, veterinarians, and processors [1,2]. However, reliance on specialized, high-yielding breeds (Holstein cattle) exposes it to climatic risks [3,4,5,6]. Indeed, recurrent heat stress has already caused 14–35% declines in milk yields, while diminishing surface water availability—especially in the Southwest—is driving up cooling and irrigation costs [7,8]. Moreover, extreme climatic events threaten feed crops [9,10]. The convergence of these factors suggests a challenging financial future, with potential losses for the U.S. dairy industry projected to surpass USD 1.7 billion by 2050 [11]. Underscoring these climate-driven threats are ongoing market volatility and policy shifts, necessitating more localized research aimed at understanding their effects on productivity [12,13]. Robust county-level data would enable both producers and policymakers to invest in adaptive measures that enhance the resilience and sustainability of U.S. dairy production.
Genetic selection in modern dairy farming, while significantly increasing milk yields, has also exacerbated the physiological vulnerability of dairy cows to heat stress due to their narrow thermoneutral zone [4,14,15,16,17,18]. When ambient temperatures and humidity exceed an optimal range, cows must use considerable energy to maintain homeostasis, diverting resources that would otherwise support milk production [7,19].
Researchers often quantify these physiological challenges using the Temperature-Humidity Index (THI), widely recognized as the key metric for assessing heat stress in dairy cattle [20,21]. When the THI exceeds a threshold of 65, physiological stress symptoms begin to manifest and above a THI of 72 severe productivity losses are observed [22,23,24]. Prolonged exposure to high THI levels further reduces feed intake, milk yield, reproductive performance, and milk quality [25,26,27,28]. Ultimately, climate-induced heat stress poses substantial challenges to dairy cow welfare and to the overall industry productivity [5,29,30,31,32,33,34,35].
A growing body of agro-economic studies focuses on the effects of climatic factors on total factor productivity (TFP); however, most consider the agricultural sector as a whole [36]. Few studies go beyond technical efficiency (TE) at the farm level and decompose TFP into constituent components—such as technological progress (TP), scale efficiency change (SEC), TE change, and climatic or environmental effect indices (CEI)—to understand the underlying drivers of TFPG [37,38,39,40,41,42,43,44]. Research spanning various regions, including Australia, several Latin American countries, and the United States, suggests that TFPG is generally positive, with technological progress consistently identified as the primary driver [37,38,39,44,45]. However, there are notable exceptions, including Bravo-Ureta et al. [42], who reported negative TFPG in El Salvador, indicating that technological change combined with climatic conditions can reduce productivity.
Despite its potential importance, the impact of climatic factors on TFP in the dairy sector has received limited attention [36]. A consistent finding across this literature is that models omitting weather variables are rejected in comparative robustness tests, underscoring their explanatory power [46]. Omitting such variables can lead to biased results, as weather-induced output fluctuations may be attributed incorrectly to technical inefficiency, resulting in lower and more volatile TE scores [47,48]. Another common theme is that TFPG is positive, but climatic effects tend to slow down this growth [43,47,48,49]. Researchers have employed a range of variables to quantify climatic effects, yet the literature consistently concludes that climatic stressors reduce productivity growth, a trend projected to intensify as climate change continues [39,50,51,52]. Finally, several studies identify technological progress as the primary engine of TFP change, even as climatic factors work to offset its gains [40,46,53].
Although substantial evidence shows that climatic shocks reduce agricultural output and disrupt dairy cow physiology, the relationship between weather variability and TFP is often limited to broad regional analyses, which obscure local variation, leaving two primary gaps in the dairy productivity literature. First, most studies assess productivity at a regional level, thereby overlooking significant heterogeneity among U.S. counties [47,48,49,54]. Second, while the THI is commonly used as a heat-stress indicator, few studies utilize it to isolate the direct physiological effects of heat. Even fewer distinguish these effects from indirect climatic impacts on feed availability, which are addressed here using growing degree days (GDD) and precipitation.
This study is based on a 45-year (1978–2022) panel of publicly available, county-level data from the USDA National Agricultural Statistics Service (NASS). This extensive period and open accessibility provide a key advantage over typical farm-level analyses, which often utilize shorter periods and data with restricted access. Addressing these gaps, the primary objective of this study is to analyze TFPG across major U.S. dairy counties, with a special focus on the role of climatic factors. Specifically, we aim to assess the effects of alternative definitions of climatic factors on dairy farm productivity over time and space and identify the drivers of TFP by decomposing it into six elements: technological progress; scale efficiency change; technical efficiency change; an environmental index for direct heat stress (THI); an environmental index for indirect climatic impacts on feed ((GDD) and precipitation); and statistical noise. To achieve this, our analysis covers 179 major dairy counties distributed across the Northeast, Midwest, Southeast, and West, allowing us to capture significant geographic and environmental heterogeneity. Our study makes two primary contributions to the literature. First, we conduct a county-level analysis that decomposes dairy TFP into separate elements, explicitly distinguishing the direct effect of heat stress, measured by THI and indirect climatic impacts on feed captured by GDD and precipitation. Second, we report regional heterogeneity in these components across major U.S. dairy regions, insinuating the importance of localized adaptation policies.
The last four decades have seen extensive on-farm adaptation of new technologies and management practices in the U.S. dairy sector. On the biological side, technological progress has involved sustained genetic improvement in Holstein dairy cattle and, more recently, targeted breeding for greater resilience and thermotolerance [4,16,55,56,57], alongside the use of selected biotechnologies such as recombinant bovine somatotropin (rBST) [58]. On the physical and management side, adaptation has relied on investments in cooling and environmental control [7,8,52,59], and on advances in barn design and cow comfort that support both productivity and animal welfare [60,61,62]. Another dimension has been the diffusion of on-farm microclimate and herd-level monitoring tools that enable more precise heat-stress management [63,64]. Data limitations do not allow us to explicitly model individual technologies; their combined effects are reflected in the technological progress component and the climatic indices. Our analysis indicates that how individual biological, physical, and informational innovations interact with climatic stressors at the farm level is an important area for future work.
The remainder of this paper is organized as follows: Section 2 details the dataset and empirical methodology used. Section 3 presents the results and discussion, and Section 4 offers a summary of our findings and concluding remarks. Given the vital role of productivity in shaping agricultural policy, our findings provide insights for policymakers and stakeholders in the United States and beyond.
2. Materials and Methods
2.1. Materials
This study utilized county-level data from the USDA Census of Agriculture, published every five years by NASS. The dataset spans the years 1978 through 2022 and includes comprehensive county-level information for all U.S. counties with at least USD 1000 in annual farm sales [65]. These data provide detailed information and broad geographical coverage offering a valuable alternative to farm-level analysis that relies on data that is often difficult to access [58,66,67,68,69,70].
The county selection process followed the methodology of Njuki and Bravo-Ureta [69], which prioritizes the top 100 counties in terms of dairy cow inventories as reported in the “State and County Rankings” volume of each agricultural census. Following their approach, we also included counties where dairy sales constituted at least 30% of total agricultural sales for each census year included in the analysis. This procedure yielded 179 dairy counties (Figure 1), and after excluding observations with missing data, our final unbalanced panel comprised 1783 county-level observations. Those counties are distributed across four regions: Midwest, Northeast, South, and West, as shown in Figure 1, following the definitions used by the USDA NASS to collect agricultural census data [65].
Our analysis employed carefully defined output and input variables to ensure robust estimation. We consistently measured milk equivalent, derived by dividing the total value of dairy product sales by the relevant state-based unit milk price for the period. Counties in our sample produced an average of 264,424 tons of milk equivalent annually (Table 1).
We disaggregated conventional inputs as follows: dairy cows represent the average number of dairy cows per county, averaging 29,856 head. However, this national average mask significant variation across regions and census periods, as detailed in Table A1 of Appendix A. Notably, the West leads in dairy cow population, which has seen substantial growth over the past 45 years, while the South consistently exhibits the lowest numbers. Farm labor accounts for both hired and contract labor hours, averaging 1.1 million hours annually per county. Farm expenses encompass outlays for livestock inputs, fertilizers, and chemicals utilized explicitly in dairy farming, averaging USD 14.6 million per county (Table 1). We measured concentrate feed (16% protein) at an average of 120.8 tons per county. Finally, we represented machinery by the value of equipment associated with milk production, averaging USD 67.4 million per county per year (Table 1). We deflated all monetary figures, including farm expenses and machinery values, to constant 2022 dollars using the Price Indices and Implicit Quantities of Farm Output and Inputs dataset published by the USDA Economic Research Service (ERS) [71].
We utilized county-level weather data for the period 1978–2022 from the Visual Crossing Weather Database (VCWD) [72]. From this database, which provides consistent information for two-kilometer grid cells, we compiled daily minimum, maximum, and average temperatures, relative humidity, and monthly precipitation. We employed three alternative variables to represent climatic effects: growing degree days (GDD); annual precipitation; and the Temperature-Humidity Index (THI) load (Table 1).
GDD measures cumulative heat exposure for a crop within a specific temperature range, calculated by summing daily heat over the growing season (1 April to 30 September). Specifically, GDDi = − T_lb_}di, where T_max_ and T_min_ are daily maximum and minimum temperatures (°C), T_ub_ (30 °C) and T_lb_ (10 °C) are the upper and lower temperature bounds, and D is the number of days in the growing season [70]. Maize was selected as the reference crop for GDD calculations because corn silage is the primary forage source in U.S. dairy rations and is highly sensitive to thermal conditions [58]. The inclusion of GDD is critical because it directly reflects the thermal requirements for feed crops, is well-established in the scientific literature, and is essential for capturing the environmental heterogeneity that drives regional differences in dairy TFP [51,73]. Consequently, it provides a pertinent measure of the interplay between climatic effects and productivity. The annual county-level GDD average is 1578. We also included annual precipitation in our model as another key climatic variable. Sufficient precipitation ensures a reliable water supply for both direct cow consumption and the cultivation of feed crops [54,58]. Conversely, drought conditions severely constrain milk production by reducing forage yields, exacerbating heat stress, and increasing operational costs [58]. Therefore, this variable serves as a key indicator for feed availability and environmental stress, which are fundamental drivers of dairy productivity [74,75]. The average annual precipitation is 717 mm (Table 1).
The THI load is the cumulative sum of THI when it exceeds 72 annually and can be defined as Equation (1):
where T_dbit_ is the dry-bulb temperature (°C) and RH_it_ is the relative humidity (%). The THI load is defined as the cumulative sum of daily THI values above a threshold of 72. Formally,
where THI_max_ is the maximum daily THI on day d, and THI72 is a threshold value of 72 [16]. Because the physiological impact of heat stress is cumulative, simple annual averages can be misleading by masking the effect of intense, prolonged heat waves [22,55]. Therefore, following Key and Sneeringer [30], we employed the THI load, which quantifies both the duration and magnitude of heat exposure above a critical threshold, offering a more biologically relevant measure of the total heat burden experienced by dairy herds. The annual average THI load is 101 (Table 1). Table A2 of Appendix A presents the mean THI load values by year and region. Although the South exhibits a higher average THI load, the West has experienced a more pronounced increase, particularly over the last decade.
The empirical framework accounts for unobserved heterogeneity by grouping counties into the four regions and incorporating regional fixed effects (Table 1). To account for temporal shifts arising from technological progress, we also included time fixed effects for each census year (Table 1) [41,49].
2.2. Methods
2.2.1. Cobb–Douglas Stochastic Production Frontier
This study employs the stochastic production-frontier (SPF) approach to examine how different variables influence milk output and productivity across U.S. dairy counties. Since Aigner et al. [76] and Meeusen and van den Broeck [77] independently introduced the SPF framework, it has become a staple in agricultural economics, including dairy research [78,79,80,81,82]. Within this SPF framework, we adopt a Cobb–Douglas (CD) functional form rather than the more flexible translog (TL) specification. In principle, the TL can capture richer interaction effects among inputs, but the CD specification, under suitable parameter restrictions, guarantees global non-negativity and monotonicity, and its coefficients can be interpreted directly as output elasticities, which facilitates economic interpretation [38,53]. Moreover, TFP measures derived from a CD function satisfy the transitivity property, an essential requirement from index-number theory for consistent productivity comparisons across regions and over time [38,53]. So, we adopt the C-D considering these conceptual issues and its good performance in related stochastic frontier applications.
We make a distinction between production technology and the characteristics of the production environment. O’Donnell [53] defines the “production technology “as a way, method, or system that transforms inputs into outputs” (p. 2). In contrast, environmental characteristics are variables that are outside the firm’s control but physically involved in the production process (e.g., weather, geographic location, topography). Following O’Donnell [53], the production technology can be expressed as (Equation (3)):
To address the inflexibilities of the traditional CD model, which imposes fixed output elasticities across all observations, we apply a Random Parameters Cobb–Douglas specification [41]. We then employ the estimates of the random parameters (RP) to examine the effects of climatic variables on the productivity of dairy counties. Unlike traditional SPF models, the random parameters model also accounts for technological and environmental heterogeneity [83,84]. With this approach, we measure the heterogeneous responses of each dairy county to both conventional inputs and changes in climatic effects [41]. Following Julien et al. [85] and Bravo-Ureta et al. [42], a CD representation of the RP model can be expressed as:
where y_it_ is milk output for the ith county in year t, is a region level fixed effect, and X_nit_ are conventional inputs, namely number of cows, feed, labor, machinery-equipment, and farm expenses. Z_jit_ are climatic effects, namely GDD, annual precipitation, and THI load. represents time dummy variables corresponding to the specific census years: 1978, 1982, 1987, 1992, 1997, 2002, 2007, 2012, 2017, and 2022, which are included to capture technological change. The Greek characters β, ρ, and γ are parameters to be estimated. The random parameter α_i_ ≡ α(z_i_^^) in the model shows the effect of unobserved heterogeneity that does not change over time, while β_it_ ≡ β( ) and ρit ≡ ρ ( ) capture time-varying heterogeneity, allowing for county-specific responses to inputs and climatic effects. Their respective distributions are jit ∼ N ( j, ), ∼ N (β_n_, ) [42]. Finally, v_it_ is the statistical noise with a normal distribution, v_it_ ∼ N (0, σ^2^_v), and *u_it is an inefficiency term with a half-normal distribution, u_it_ ∼ N + (0, σ^2^u) [49,76]. We calculate the technical efficiency component for each observation given by TE_it_ = exp (−û_it_), following Jondrow et al. [86].
2.2.2. Total Factor Productivity Index (TFPI) and Components
Productivity is equal to output quantity relative to the total amount of inputs [49]. A TFPI is given by the ratio of the productivity of firm i in period t relative to firm k in period s and can be represented as:
In this analysis, we use the general index developed by O’Donnell [53] to calculate sources of productivity growth in the U.S. dairy counties. Following Bravo-Ureta et al. [42], the complete representation of TFP under the random parameter model is given by:
where the first component on the right-hand side of the equation is the output-oriented technological progress index (TP), which captures the change in productivity associated with shifts in the production frontier due to the discovery and application of new technologies over time. The second component is the output-oriented scale efficiency index (OSEI), reflecting changes due to economies of scale, where is an estimator of the mean of the distribution of β_nit_. The third component is the environmental index (EI), which accounts for weather variables and time-invariant region-specific heterogeneity. We express this index as two components: EI-THI, which captures changes in growth due to heat stress (THI > 72); and EI-FEED, indicating changes in productivity that come from GDD and precipitation. The fourth component is the output-oriented technical efficiency index (OTEI), which measures management performance and reflects movements toward the production frontier. The last component is the statistical noise index (SNI), which identifies changes in productivity that are not accounted for by the model.
To measure the change in TFP over time for all components, we take one county, Maricopa, AZ, in 1978, as the reference point. We selected this county systematically, as it is the first observation when our sample is sorted by the Federal Information Processing System code, thus ensuring an objective point of reference. It is important to reiterate that under O’Donnell’s [53] index-number framework, transitivity ensures that the choice of a reference observation does not affect the relative rates of change in TFP across counties or over time. The rate of change in TFP is revealed by comparing the other years and all remaining counties with the reference county in 1978.
3. Results and Discussion
Before presenting the detailed analysis of TFP measures and their components, this section outlines the model specifications we compared, which support selecting the preferred Random Parameters Stochastic Production Frontier (RP-SPF) model through hypothesis testing.
3.1. Definition of Models and Hypothesis Tests
Having estimated various models, we report the results for the two most relevant specifications.
Model 1: Pooled-SPF that includes the conventional inputs, three climatic indicators, regional fixed effects, and time dummy variables. Model 1 assumes that all dairy counties share the same production technology.
Model 2: RP-SPF that incorporates the conventional inputs, climatic effects, regional fixed effects, and time dummy variables. The parameters for cows, labor, farm expenses, feed, and machinery inputs, and GDD, precipitation, and THI load are assumed to be random, while the parameters of the regional fixed effects and time dummies are treated as non-random (Equation (4)). Model 2 extends Model 1 by incorporating random parameters, which account for unobserved county-level heterogeneity.
To test the hypothesis that the estimated coefficients of the conventional inputs and climatic effects are random, we compared Model 1 (Pooled-SPF) and Model 2 (RP-SPF). Model 1 is nested within Model 2, serving as the restricted (R) specification, while Model 2 is the unrestricted (UR) specification. Because the R model is nested in the UR model [87], we use a likelihood ratio (LR) test (Table 2). The LR statistic is calculated as LR = −2[lnL_R_ − lnL_UR_], which follows a chi-square (χ^2^) distribution with J degrees of freedom (the number of restrictions) [87]. We compare this statistic with the critical values from Kodde and Palm [88] and reject the null hypothesis (H_0_) of the restricted model if the LR value exceeds the critical value [89]. The likelihood values for Models 1 and 2 are 1826 and 2055, respectively (Table 3), with 8 degrees of freedom corresponding to the five conventional inputs and three climatic variables for which random parameters are estimated. The LR statistic is 458 (p < 0.001), which exceeds the 5% critical value; therefore, we reject H_0_ and select the unrestricted Model 2 for the analysis.
We clustered the standard errors at the county level and estimated all RP models using simulated maximum likelihood methods with 100 Halton draws in NLogit version 6 and then generated the means of the random parameter estimates for the analysis by averaging the individual parameters across all counties and time periods. Table 3 reports the estimated means of the random parameters, which incorporate the conventional inputs and the climatic effects. The table also provides model fit statistics such as the Akaike Information Criterion (AIC) and Log-Likelihood values as well as the parameters for the error components.
3.2. Parameter Estimates of the RP Models
In our preferred specification (Model 2), the parameters of the five conventional inputs are positive, as expected, and are significant at the 1% level. We used the Wald test for constant returns to scale (i.e., Σβ_m_ = 1) [87,88]. This hypothesis is rejected (χ^2^ = 49.90, p < 0.001), indicating mildly increasing returns to scale for Model 2 (Table 2). The parameters for all climatic variables are also significant. The time dummies are significant at the 1% level, as are those for the regional fixed effects in Models 1 and 2 (Table 3).
We also estimated alternative model specifications. One alternative model replaced regional fixed effects with county-specific latitude and longitude coordinates, which performed poorly and was discarded.
By contrast, Model 2 incorporates regional fixed effects instead of continuous coordinates [90,91]. This approach helps address collinearity and allows the model to account for time-invariant, region-specific factors that the geographical coordinates may not capture. We therefore treat Model 2 as the most robust specification and focus on it below.
The output elasticities from Model 2 confirm that the number of cows is the most important conventional input (0.720), consistent with the literature [49,92,93]. A 1% increase in herd size raises milk output by 0.72%. The elasticities for labor (0.105), farm expenses (0.066), feed (0.056), and machinery (0.084) are also positive. The elasticity of machinery exceeds that of concentrate feed, reflecting the capital-intensive nature of modern U.S. dairy systems [58]. Additional investments in machinery are therefore associated with stronger output responses than comparable increases in feed use.
The parameters for the climatic variables are also statistically significant. Favorable conditions for feed cultivation, such as a 1% increase in GDD and precipitation, are associated with output increases of 0.05% and 0.01%, respectively, which aligns with existing studies [49,58]. These factors indirectly affect milk output by influencing feed quantity and quality [58,94]. The estimated elasticity of output with respect to the THI load is −0.011, which appears small at the mean. Given that the average THI load is about 101 (SD = 183; Table 1), a one-standard-deviation increase in THI load is associated with an output reduction of about 2.0% (183 × 0.011 ≈ 2.0%). This result is consistent with earlier research showing that higher THI levels reduce milk output [30,47,55].
Finally, both regional and time-fixed effects (the South region and the 1978 census year as the excluded categories) are statistically significant. The likelihood ratio test confirms the importance of regional fixed effects (p = 0.000), with results showing that the West has the most favorable production environment (relative to the South baseline). This pattern is consistent with the long-run geographic transformation of the U.S. dairy sector, whereby cow inventories have shifted from traditional production regions in the Midwest and Northeast towards Western states [67]. Western dairy production has moved toward a highly industrialized model, characterized by larger herd sizes and substantial investment in large-scale, technologically advanced facilities [67].
Similarly, the positive time-fixed effects for each census period after 1978 point to robust technological progress with estimated productivity growth peaking at 74% in the 2017–2022 period relative to the 1978 baseline.
Regarding the model’s error structure, the parameter lambda (σ_u_/σ_v_) is equal to 2.99 (Table 3), indicating that the variance of the one-sided error component is slightly larger than that of the two-sided component [95]. This highlights the statistical significance of technical inefficiency on milk output changes [38]. The annual average TE score estimated from Model 2 is 92.4% (Table 4). This score is consistent with the high TE scores of previous dairy studies reported in the meta-analysis by Bravo-Ureta et al. [96] and Moreira Lopez and Bravo-Ureta [80].
The analysis indicates that U.S. dairy counties are, on average, highly efficient, achieving 92.4% of their maximum potential output given their resources and technology. This high average is driven by a large concentration of well-managed operations clustered at the top end of the performance scale, suggesting best practices in resource use [97,98].
However, as illustrated in Figure 2, the high average TE masks a critical efficiency gap. While the top-performing counties achieve high TE scores (>0.98), Figure 2 identifies a “tail” of underperforming operations, in several counties (e.g., Lafayette, FL; Orange, VT) operating at TE scores below 0.63, far below the national average. From an industry perspective, this “tail” of underperforming counties represents an untapped potential for productivity and output growth.
3.3. Total Factor Productivity (TFP) Measures and Its Components
The analysis of the TFP index (TFPI) and its components (Table 5) reveals a strong increasing trend in U.S. dairy productivity after 1978, peaking in 2022. Following O’Donnell [53], the TFPI is calculated as TFPI = TP × OSEI × EI-FEED × EI-THI × OTEI × SNI. For example, the U.S. TFPI for the base year of 1978 was calculated using the 1978 components from Table 5 as follows: TFPI = 1.000 × 0.966 × 0.931 × 1.018 × 0.999 × 1.037 = 0.949). The highest five-year growth rate of 15.5% was achieved between 2007 and 2012. This growth was primarily driven by TP, which displayed a consistent upward trend and reached its highest level (2.99%) in the final period. The largest surge in TP occurred between 1978 and 1982 at 60.1%.
This progress reflects an evolution of productivity–enhancing innovations promoting adaptation–oriented investments [56,58]. Earlier decades were characterized mainly by gains in genetics and mechanization [4,16], whereas from the 1980s onward the rapid expansion of cooling infrastructure, improved barn environments, and microclimate monitoring systems became critical for sustaining TFPG under rising heat loads [7,8,62,63,99]. The pivotal role of TP is underscored in the 2002–2007 period, where a temporary decline in TFPG was directly attributable to a 4.6% decrease in TP. The decline in the 2002–2007 period coincides with rapid structural change in the U.S. dairy sector, as production shifted toward larger capital-intensive farms in the West [100]. While these operations ultimately achieve lower unit costs, the transition process entails short-run adjustment—such as farm exits and capital reallocation—that can temporarily depress productivity. Our results are therefore consistent with an evolving structural reorganization rather than a permanent reversal of technological progress.
In contrast to the steady gains from technological progress, average OSEI fluctuated before peaking in 1982 and then exhibited a slight decreasing trend. Similarly, the average OTEI showed a slight decline over the 45 years, suggesting a modest decrease in managerial performance, with the largest drop (2.7%) occurring between 2017 and 2022.
The environmental indices reveal mixed effects on productivity. The average EI-FEED showed an overall increasing trend despite modest fluctuations, including a 0.8% decline in the 1987–1992 subperiod that corresponded with declines in GDD averages (Table A3 in Appendix A). On the other hand, EI-THI declines overall but exhibits sharp fluctuations over time (e.g., −1% in 1992–1997). The most pronounced decline occurs between 1992 and 2002. These swings—for instance, the drop between 1992 and 1997—are consistent with observed variation in THI load across regions (Table A2 in Appendix A), reflecting the influence of heat waves and interannual climate variability on the environmental component of TFP. Based on this overview, we now turn to the year-to-year dynamics of TFPI and its components and examine the main drivers of these changes.
Table 6 presents the average annual TFPG for the U.S. and its regions over the 45-year period. The average annual TFPG for the U.S. was 1.59%, driven primarily by TP at a 2.52% per year (Table 6). These results are consistent with other dairy studies. For example, Njuki et al. [49] found TFPG of 2.16% in Wisconsin, while a nationwide study by Njuki [58] identified annual TFPG of 2.51%.
While technology was the primary driver of TFPG, other components made smaller contributions. Climatic effects played a dual role; favorable feed conditions (EI-FEED) provided a modest annual increase of 0.02% to TFPG. This positive effect is directionally consistent with Njuki [58], who found a more pronounced contribution of 2.2% from feed-related climatic factors. The difference in magnitude likely stems from our use of county-level census data over a longer timeframe versus Njuki’s farm-level survey over a shorter period, underscoring the value of examining climatic impacts at different aggregation scales to capture both granular farm dynamics and broader long-run regional trends.
In contrast, the impact of heat stress (EI-THI) on the rate of productivity growth was minimal, reducing TFPG by 0.01% annually. This low magnitude is consistent with the hypothesis that widespread on-farm adaptation strategies have mitigated part of the negative effects of heat stress [101,102]. However, other studies show that these strategies have not fully neutralized the impact of heat, which still results in substantial economic losses for the U.S. dairy industry, with annual estimates ranging from USD 897 million to USD 1.5 billion [3,57,62,103]. Meanwhile, efficiency change has been slightly negative. OSEI decreased at an average rate of 0.01% per year, consistent with other recent studies [49,58].
The slightly negative trend in OTEI (−0.091%) suggests a modest weakening in average managerial performance rather than a sharp decline. Although our framework cannot identify specific mechanisms, recent evidence on U.S. dairy consolidation [104] points to plausible sector-level influences, including increased reliance on hired labor, uneven adoption of productivity-enhancing technologies, and the continued operation of high-cost farms. These structural features may contribute to the negative OTEI observed in our analysis. Finally, the SNI component averaged −0.83% (Table 6).
A detailed breakdown by region (Table 6) and by state (Table A4 in Appendix A) reveals significant performance disparities. The Midwest (1.74%) was the fastest-growing region, led by South Dakota (1.93%), in contrast with Iowa (1.46%), the slowest-growing state. This pattern is consistent with documented outreach efforts in the Northern Great Plains, where producers have been actively managed heat stress by optimizing environmental control systems to enhance cow comfort and productivity [105]. These initiatives focus on fan maintenance, airflow pattern detection, and animal monitoring, reflecting a proactive approach to mitigating heat stress impacts [105].
The Northeast (1.53%) shows wide variation, from Massachusetts (1.75%) to Vermont (1.20%), highlighting that farm-level strategies can be more decisive than broad climatic trends [106]. The South (1.33%) was the slowest-growing region [107], with Texas (1.80%) and Florida (0.58%), the nation’s slowest-growing state. The West (1.42%) showed the benefits of OSEI and adaptation, with New Mexico (1.88%) and Colorado (1.86%) achieving high growth despite climatic pressures, supported by sophisticated management practices [108] and advanced technologies [109].
Figure A1 of Appendix A demonstrates the counties with the highest and lowest TFPG along with a sample of selected counties. The county-level analysis highlights a wide range of outcomes, underscoring that local adaptation and resilience can play a significant role. A noteworthy pattern is that some of the fastest-growing counties also face significant climatic constraints. For example, Brookings, SD (Midwest), and Owyhee, ID (West), achieved some of the nation’s highest annual TFPG rates, at 2.23% and 2.19%, respectively (Figure A1a in Appendix A). This occurred even as they experienced notable TFP reductions from heat stress (−0.02% for Brookings and −0.06% for Owyhee) (Figure A1b in Appendix A), suggesting that good management and technology adoption can more than compensate for adverse environmental pressures. This contrasts with counties like Lafayette, FL (South), and Franklin, VT (Northeast), which saw negative overall TFPG, at −1.28% and −0.42%, respectively (Figure A1a in Appendix A). The most severe climatic challenge was observed in Weld, CO (West), where heat stress reduced annual TFPG by 0.08% (Figure A1b in Appendix A). However, the county still achieved a strong overall growth rate of 1.86% (Figure A1a in Appendix A), further demonstrating the decisive role of adaptation in mitigating climatic impacts.
Our findings show that the U.S. dairy sector has maintained positive TFP growth despite rising heat loads. The combination of robust technological progress (2.52% annually), high average technical efficiency (92.4%), and the relatively small constraints from climatic variables—particularly heat stress—suggests that farmers have implemented effective adaptation strategies [47]. Studies in the Southeast and Wisconsin demonstrate that direct mitigation measures, such as cooling systems, significantly reduce production losses and generate benefits that exceed their costs [47,110].
Our results also align with international evidence from other regions facing similar climatic pressures, strengthening the external validity of our conclusions. Research from Southern Europe, Australia, and Latin America consistently shows that, while heat stress and adverse weather slow growth, underlying productivity trends remain positive in well-adapted systems [19,28,37,44,111,112].
This resilience requires ongoing investment [59,110,113]. The limited overall effect of climate on TFP likely reflects both these direct adaptations and broader, often costly, managerial changes, including precision data systems and improved feeding strategies [59,108,109,110,113,114]. This pattern underscores the importance of management in balancing avoided losses with additional operating costs. Future research should move beyond aggregate indices to quantify the farm-level economic returns of specific adaptation strategies and identify which combinations of technologies and practices offer the greatest resilience per unit cost.
4. Conclusions
This study provides a comprehensive analysis of total factor productivity growth (TFPG) across 179 U.S. dairy counties from 1978 to 2022, with a specific focus on climatic factors. We utilize publicly available county-level panel data from USDA NASS, which offer a distinct advantage over comparable farm-level panel data with restricted access. We employ the RP-SPF framework to decompose TFPG. The analysis yielded two principal findings. First, technological progress has been the primary driver of productivity growth, averaging 2.52% annually. Second, while climatic factors pose challenges, their net impact on TFPG has been low. Our model disaggregated this into a modest positive contribution from feed-related conditions (EI-FEED) and, most notably, a low national-level drag from direct heat stress (EI-THI).
Our framework explains changes in total factor productivity; however, it does not isolate the specific contribution of individual strategies—such as cooling systems, precision feeding, or breeding for thermotolerance—from overall technological progress. Additionally, structural changes within the sector and potential implications for animal welfare are only indirectly reflected in our aggregate indices. Consequently, our conclusions regarding adaptation rely on indirect evidence.
In conclusion, this research indicates that the U.S. dairy sector has maintained a positive productivity trend despite increasing climatic pressures. The sector has successfully leveraged technological change and adaptive management to buffer the adverse effects of a warming climate. The long-term sustainability of these gains will depend on two key factors: first, future research using farm-level data to directly quantify the economic returns of specific adaptation schemes, and second, continued innovation and policies that support the diffusion of the most effective strategies across diverse U.S. dairy systems.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1UNFAO United Nations Food and Agriculture Organization Database 2025 Available online: https://www.fao.org/faostat/en/#data(accessed on 3 May 2025)
- 2USDA ERS Farm Income and Wealth Statistics—Cash Receipts by State. U.S. Department of Agriculture, Economic Research Service Database 2025 Available online: https://data.ers.usda.gov/reports.aspx?ID=4052#P 9ba 3fb 89275541 c 6a 571c 9d 216b 7745 c_2_17i T 0R 0x 0(accessed on 25 July 2025)
- 3Ouellet V. Cabrera V.E. Fadul-Pacheco L. CharbonneauÉ. The relationship between the number of consecutive days with heat stress milk production of Holstein dairy cows raised in a humid continental climate J. Dairy Sci.20191028537854510.3168/jds.2018-1606031255266 · doi ↗ · pubmed ↗
- 4Contreras-Correa Z.E. Sánchez-Rodríguez H.L. Arick M.A.II Muñiz-Colón G. Lemley C.O. Thermotolerance capabilities, blood metabolomics, and mammary gland hemodynamics and transcriptomic profiles of slick-haired Holstein cattle during mid lactation in Puerto Rico J. Dairy Sci.20241074017403210.3168/jds.2023-2387838246540 · doi ↗ · pubmed ↗
- 5Souza V.C. Moraes L.E. Baumgard L.H. Santos J.E.P. Mueller N.D. Rhoads R.P. Kebreab E. Modeling the effects of heat stress in animal performance and enteric methane emissions in lactating dairy cows J. Dairy Sci.20231064725473710.3168/jds.2022-2265837225587 · doi ↗ · pubmed ↗
- 6Kanca H. Milk Production and Composition in Ruminants Under Heat Stress Nutrients in Dairy and Their Implications on Health and Disease Watson R.R. Collier R.J. Preedy V.R. Academic Press San Diego, CA, USA 201797109978-012-809-762-5
- 7Fournel S. Ouellet V. CharbonneauÉ. Practices for alleviating heat stress of dairy cows in humid continental climates: A literature review Animals 201773710.3390/ani 705003728468329 PMC 5447919 · doi ↗ · pubmed ↗
- 8Smith T.R. Chapa A. Willard S. Herndon C. Williams R.J. Crouch J. Riley T. Pogue D. Evaporative Tunnel Cooling of Dairy Cows in the Southeast. II: Impact on Lactation Performance J. Dairy Sci.2006893915392310.3168/jds.S 0022-0302(06)72434-116960067 · doi ↗ · pubmed ↗
