Next Article in Journal
Understanding the Impacts of Financial Flows in the Landscape
Next Article in Special Issue
Can Cooperative Supports and Adoption of Improved Technologies Help Increase Agricultural Income? Evidence from a Recent Study
Previous Article in Journal / Special Issue
Making Agriculture Carbon Neutral Amid a Changing Climate: The Case of South-Western Australia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impacts of Climate Change on Livestock Location in the US: A Statistical Analysis

1
School of Public Finance and Taxation, Southwestern University of Finance and Economics, Chengdu 611130, China
2
Department of Agricultural Economics, Texas A&M University, 2124 TAMU, 600 John Kimbrough Blvd, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Submission received: 15 September 2021 / Revised: 13 November 2021 / Accepted: 17 November 2021 / Published: 18 November 2021
(This article belongs to the Special Issue Agricultural Land Use, Economics and Climate Change)

Abstract

:
Livestock production is a valuable part of US agriculture as it contributes 50% of total agricultural value. Climate change is likely a threat to livestock production, but research regarding the impact of climate change on livestock sectors is limited. This paper examines how climate change affects livestock mix and location. Specifically, we examine climate effects on grazing animals and, in particular, on beef cattle, dairy cattle, goats, and sheep. We examine this in the US based on county-level data by using fractional multinomial logit econometrics. Our results show that climate is an influential determinant of where livestock herds are located and species mix. The impacts of climate vary by species and region. We also find significant influences from geographic characteristics and animal product prices. Subsequently, we project how climate change would influence future livestock mix and location. It reveals a likely growth in beef cow land shares across most of the US with the largest gains in the northwest. We also find substitutions between species as climate change progresses with dairy cows exhibiting the largest reduction.

1. Introduction

For the past half-century, the climate has been changing and exhibits increased temperatures, altered precipitation patterns, and increased frequency of extreme events [1], and this is expected to continue. Agriculture is sensitive to climate change, and effects have already been felt across the landscape [2]. However, the effects are not uniform, with some cold regions likely benefiting and some hot regions losing [3].
Livestock production is a valuable part of US agriculture, contributing 50% of total value [4]. According to the 2017 Census of Agriculture, the US beef cattle industry had 32.5 million beef cows in inventory and USD 77.2 billion in sales, accounting for 20% of total US agricultural sales [5]. The poultry industries were the next largest livestock commodity in the US, with production valued at around USD 49.2 billion, followed by hogs and pigs at USD 26.3 billion. The sales of milk products from dairy cattle totaled USD 36.7 billion, accounting for 9.5% of total US agriculture sales. Goats and sheep contributed USD 1.03 billion, with goats totaling 2.64 million head and sheep totaling 5.2 million head.
Grazing animals are likely to be more affected by climate change than confined livestock since they are directly exposed to the climate. Some species are heat tolerant, and some are cold resistant. Some may have a wider range of tolerable temperatures while some have rather narrow ones. For example, dairy cows’ best range as measured by the thermo-neutral zone (TNZ) is estimated to range from 5 to 25 °C [6]. Conditions hotter than the TNZ have negative effects on milk production, animal growth, fertility, and health among other effects altering hormonal status, nutrient absorption, and immune functions [7]. The economic loss caused by heat stress was estimated to be as high as USD 897 million for the dairy industry [8]. Cold stress is also a concern, causing dairy cows to increase their metabolic rate to supply more body heat which may increase dietary requirements. Beef cows may also suffer from both heat and cold stresses, requiring more energy for maintenance. However, findings indicate that beef productivity loss caused by varying temperatures is smaller than those for other animals because beef cows have lower metabolic rates and a wider TNZ range [8]. However, excess winter precipitation may increase housing needs [9]. Goats and sheep are considered to be heat tolerant, but their body reserves are not sufficient for adapting to multiple stressors such as thermal, nutritional, and walking stresses [10]. A remarkable drop in body sizes has been observed for sheep and goats from the northern to southern areas of the Mediterranean area [11].
Thus, the question of how to reduce the economic loss associated with climate change is critical for livestock producers. One option is altering operations to accommodate the local climate, and today this does not only involve adjustment to an altered local climate but also involves an evolving climate. In the literature, the adjustment to a changing climate is called adaptation [12]. Livestock producers have been observed to adapt by altering management, land usage, livestock species mix, crop/grassland mix, feeding patterns, housing, and shade provision [13,14,15,16,17,18].
Conceptually, we can view livestock management as a multi-stage decision process [19]. In the first stage, farmers choose whether or not to raise livestock. In the second stage, they choose the species mix to utilize and where to locate them. In the third stage, they choose stocking rates, diets to use, grazing intensity, supplemental feeding, housing, etc. All of this, we presume, is conducted in pursuit of profit maximization. Adaptive activities involve climate change accommodating altercations in the decisions across these stages including species mix, the spatial distribution of animals, and animal density/stocking rate.
Here, we explore part of the adaptive decision-making space by examining how US farmers adapt grazing livestock species mix to alternative climate conditions. In conducting this study, we examine US county-level data by using a multinomial logit model for land use shares of four major species—dairy cows, beef cows, goats, and sheep. We also project how future livestock mix and location alter under different climate change scenarios. Our resultant hypotheses guiding this study are presented in Table 1.
This paper is organized as follows. The following section reviews prior literature researching farmers’ adaption associated with livestock management. Section 3 provides a detailed description of the multinomial logit model and data used in the estimation. Section 4 presents the estimation results and projection results under different climate change scenarios. Section 5 concludes this paper.

2. Literature Review

The effect of climate change on crop production locations has been frequently examined [20,21,22,23]. The prevalent finding is that climate change is causing crop production to shift to higher latitudes and/or elevations [18].
A few studies have examined effects on livestock land use and location. Some found that grassland would substitute for cropland when the climate is not beneficial for crop production [24,25,26]. Mu et al. [16] examined how climate affected land allocation to crops versus grazing at the US state level, finding U and inverted U shaped relationships with temperature and precipitation. Namely, as temperatures warm from low levels or rainfall increases from very dry conditions, land shifts from grassland to crops, but when temperatures or rainfall rise past thresholds, then land shifts out of croplands to grass. Seo and Mendelsohn [19] found that warming and drying would increase livestock ownership but that an increase in precipitation would decrease it. Mu et al. [27] projected that climate change would cause an increase in grazing land share in the Southern US and the inland pacific northwest but a decrease in the Northern US. However, US-wide studies of land allocation among livestock species have not been conducted.
Climate influences on species mix have been examined in a few regions [19,28,29]. Seo and Mendelsohn [19] found that in the probability of choosing beef and dairy cattle decreases as temperature increases in 10 African countries, with the probability of owning goats and sheep increasing. They also found that the probability of choosing chickens exhibits a hill-shaped response to temperature and that as precipitation increases, cattle and sheep decrease while goats and chickens increase. A related study conducted in South America generally found that a drying climate would increase sheep while decreasing beef cattle, dairy cattle, pigs, and chickens [29].
Breed choices have also been discussed as adaptation strategies [15,30,31]. Empirically, Zhang et al. [17] found substitutions between Bos indicus and more heat tolerant Bos Taurus breeds. Mader et al. [30] found positive impacts on beef cattle with increasing temperatures and negative impacts on milk production. They also suggested that animals can adapt to small changes in climate, but adjustments in management strategies are needed under larger changes.
Climate also influences decisions on herd size [32,33,34,35]. Farmers in Australia adapted to a hot arid climate by increasing sheep and beef cattle herd size [34]. In the US, the cattle stocking rate was found to decrease when summer temperature-humidity index increased [16]. However, no studies have used observed US herd data to examine the impact of climate change on farmers’ decisions on livestock species mix and location. Moreover, most previous studies used either US state-level data to investigate adaptation [16] or county-level data to investigate regional adaptations [17,36]. Few studies have utilized county-level data to examine the climate effects nationwide. This paper will fill that gap.

3. Methods and Data

3.1. Estimation Approach

In order to examine how climate influences livestock land shares, we estimate an econometric model across grazing species. Specifically, we follow Seo et al. [29] and assume that farmers maximize profit π * as follows:
π * = M a x j π j C , G , E
where π j denotes the profit arising from raising species j, which is a function of climate (C), geographic characteristics (G), and economic conditions (E).
In turn, the profit of species j is calculated as follows:
π j = π C , G , E | D j = 1 = P j Q j C , G , E c o s t j C ,   G ,   E
where D j is a dummy variable indicating the species j is chosen, P j is the price of the jth type of livestock product, Q j is a production function where production depends on C, G, and E, and c o s t j C ,   G ,   E is the production cost. In turn, the optimum production level ( Q j * ) can be written as follows:
Q j * = a r g m a x   π j C , G , E
where Q j * is the optimal livestock population that yields the highest profit.
In turn, the optimal proportional land share for livestock j in county i ( s i j * ) is as follows:
s i j * = Q j * A U E j G R i j Q j * A U E j G R i
where A U E j is the animal unit equivalent for livestock species j. The AUE parameter adjusts for the fact that different grazing livestock species differ in size and converts livestock numbers to a common measure of animal units. One animal unit (AUE = 1) is defined as a 454 kg cow with an unweaned calf. In turn, following Redfearn and Bidwell [37], the amount of forage required by seven goats is equal to that used by the cow-calf pair (AUE = 0.17 for a goat) as is the amount used by five sheep (AUE = 0.2 for a sheep). G R i is the grazing rate for county i = 1 , , N , representing the grassland size per animal unit. The grazing rates on a state basis are drawn from FASOMGHG [38]. The numerator of the land share equation calculates the amount of grazing land that livestock species j requires in county i. The denominator calculates the total amount of grazing land required by all livestock species in the county. Thus, s i j provides the share of grazing land used by each livestock species and, by construction, equals one when summed across all grazing livestock types.
The proportional land shares of livestock reflect the farmers’ selection preferences toward species including climate adaptation. Adaptation to climate can be measured by the change in land share in response to a change in climate. If climate change causes one species to become less profitable relative to another, the farmer could adjust land share between species. Product prices and geographic factors would also have an influence. Given this, we estimate s i j * as a function of C, G, and E.
To estimate this function, we use a fractional multinomial logistic regression that we estimate via the quasi-maximum likelihood (QML) following [16,18,39,40]. The resultant conditional means for the land shares equal the following:
E s i j t | x i t = G j x i t ; β = exp x i t β j k = 1 J exp x i t β k ,     j = 1 , , J
where s i j t   is the optimal land share of livestock type j in year t = 1 , , T , and county i = 1 , , N . x i t is a vector of the C, G, and E explanatory variables.
If we set livestock species J to be the baseline livestock species by normalizing β J = 0 , the conditional mean for the land share of the baseline livestock species J is calculated as follows.
E ( s i J t | x i t ) = G J x i t ; β = 1 1 + k = 1 J 1 exp x i t β k
Moreover, the conditional mean for the land share for other livestock is as follows.
E ( s i j t | x i t ) = G j x i t ; β = exp x i t β j 1 + k = 1 J 1 exp x i t β k ,   j = 1 , , J 1
The above equations cause the conditional expected land shares to add up to one ( j s i j t = 1 ) and to fall between zero and one ( s i j t 0 , 1 ) given that Pr s i j t = 0 | x i t 0 and Pr s i j t = 1 | x i t 0 for j = 1 , 2 , , J ).
Additionally, we use the fractional multinomial logit since it avoids issues with the independence of irrelevant alternatives [40] and extreme zero and one values [41].
In the estimation, the quasi-maximum log-likelihood function is utilized following Papke and Wooldridge [42].
l i t β = j J s i j t log G j x i t ; β
To obtain the QML estimator for β , we maximize l i t β resulting in the first order condition for species j land share in county i as follows.
l i β β j = i = 1 N x i t s i j t G j x i t ; β = 0
That expression provides the marginal effects that a one-unit change in x has on the log of the probability of choosing a livestock species over the base species. In order to develop direct interpretations, following Long and Freese [43], we calculate the average marginal effects (AME). AME is interpreted as the marginal impact of one unit change in the explanatory variables on the land shares. Note that since the land shares always add up to one, this means that an increase in the land share for one species is accompanied by decreases in shares for the other species. Thus, the sum of marginal effects across all species is always zero. For continuous explanatory variables, the AME is given by the following:
E [ s i j t | x i t ] x i t m = N 1 i = 1 N β j m G j x i t ; β G j x i t ; β k = 1 J 1 G k x i t ; β β k m
where x i t m is the value of mth continuous explanatory variables in county i for year t. β j m is the mth coefficient within the estimated equation for livestock species j. For discrete explanatory variables, AME is calculated as follows.
Δ E [ s i j t | x i t ] Δ x i t m = N 1 i = 1 N G x i t m β j m + β j m G x i t m β j m
Although the fractional multinomial logit model does not explicitly deal with the panel structure of the data, we include both regional and time-fixed effect dummy variables as conducted in the panel model.
E ( s i j t | x i t ; β ) = exp x i t β j + c i + z t k = 1 J exp x i t β k + c i + z t ,     j = 1 , , J
Spatial fixed effects variables are included at a more aggregate regional level as adding them for all counties would eliminate much of the variation in the data and compromise our ability to estimate the effects of climate [44]. The regions we used are given in Table 2 and mapped in Figure 1. The region classification coincides with those used in FASOMGHG [38] and is also consistent with the boundaries of USDA farm production regions [45]. Regional fixed effects are included to account for regionally specific agriculture-related characteristics such as soils, production resources, transportation access, and policy factors that may influence the relative share of each livestock species. We also include a time-fixed effect to capture region-invariant but time-specific confounding factors such as technological progress.

3.2. Data

We estimate our land share equation by using USDA county-level data on the inventory numbers of beef cows, dairy cows, goats, and sheep. In performing this, we realize that shifts in livestock species are longer-term decisions and occur in response to longer-term developments (as argued in Mu et al. [16]). Thus, we use longer time intervals between decisions, allowing time for adjustment. As a consequence, we use census data that are collected every five years for herd, geographic, and economic data [46] and three-year moving average climate data [47]. We believe this allows us to observe long-term species land use shifts in response to long-term climate movements.
The US county-level livestock inventory data were drawn from the Census of Agriculture that is released every five years and available in the USDA National Agricultural Statistics Service (NASS) [46]. The census years we used span from 1997 to 2017 (5 census years fell in that interval). For goats, 1997 data were not available. Land shares for dairy cows are only considered in the counties east of the Mississippi River since dairy cows are generally raised in confinement in the Western US. As a consequence, we also performed separate estimations for the Eastern and Western US. Table 3 summarizes the average land share by species and region. Beef cows are the most common species in most regions, followed by dairy cows, and then sheep, with goats being the smallest groups.
For climate data, we used county-level weather station data. In order to prepare the data, we proceeded as follows. First, we obtained weather station data from the United States Historical Climatology Network (USHCN) [47] and extracted monthly maximum and minimum temperatures (in °C) and monthly precipitation (in mm) for the years from 1995 to 2017. Then, we computed county-level measures of the following: (a) for counties with multiple stations, we took the average; (b) for counties with only one station, we used data for that station; and (c) for counties without a station, we linearly interpolated values from neighboring counties using ArcGIS. Moreover, in order to use data more representative of long-term climate rather than short-term weather, we computed 3-year moving averages by using data for the three years preceding each census year. Note that in our estimating equation we also included squared terms for the climate variables, allowing the equation to incorporate nonlinear climate effects on livestock mix and allowing identification of potential critical thresholds.
For economic factors, we used county-level measures of species product sale possibilities. For the cow-calf operations, we used the average calf sale value which we derived by dividing total sales revenue by the number of calves sold in order to obtain the average calf sale value. For dairy cows, we multiplied the milk price by reported milk yield per cow in order to obtain the average milk revenue per dairy cow. For goats, we divided sales revenue by the number of animals sold. For sheep, we calculated the wool sales per animal by multiplying wool price by wool production per animal and then added in the average lamb sales revenue per ewe. The data for this were drawn from the Census of Agriculture [46] and the Animal Products database of the USDA Economic Research Service (ERS) [48].
Variables descriptive of regional geographic conditions that we used included the amount of land classified as grassland in a county and a regional fixed effect dummy for the market regions shown above. Grassland available is used as a proxy of the resource base available for raising livestock. Regional dummies for 10 market regions are used to control for systematic regional fixed effects. Year dummies for 5 census years are used to control for year-specific effects.
Table 4 summarizes the descriptive statistics for the explanatory variables. Table 5 presents core explanatory variables values by region. There, we observe that Western US had less precipitation and a wider temperature range. We also observe that the sales value of livestock does not vary much across regions. However, the western region has much larger grassland endowments.

4. Results

4.1. Livestock Species Mix

Since we assume that the dairy cows in the east predominantly use grazing land while those in the west are generally kept in confined conditions, we conducted estimations for the east and west separately. The results where beef cows are the baseline species are presented in the Appendix A. The coefficients reported there provide the effects on the log of the land share for each species relative to that for beef cows. Computed average marginal effects are provided in Table 6 and Table 7.
The coefficients here provide the effects on land share for each species. The results show that most explanatory variables exert significant effects and that all four species exhibit nonlinear relationships with at least one climate variable, which confirms Hypothesis 1.
Figure 2 and Figure 3 show graphs of the estimated relationship for each species that relate the grazing land share with climate variables, holding all other variables at their mean level. Since beef cows have larger land shares, a secondary y-axis for beef cows is displayed on the right.
We first examine average marginal land share effects in Eastern US. There, we find maximum temperature exerts positive effects on beef cow land share and negative effects on dairy cow and goat land shares, indicating a potential substitution between beef and dairy cows/goats. The beef cow share approaches one when the annual maximum temperature reaches high levels, and dairy cow and goat shares increase when the annual maximum temperature falls to low levels. Sheep display inverse U shape curves with their land share increasing from low temperature levels and then peaking at 22 °C and decreasing from thereon. This shows differential responses by species verifying Hypothesis 2.
For annual minimum temperature, we observe smaller beef cow land shares as the minimum becomes higher. Unlike beef cows, the sheep land share increases as the minimum temperature increases. An inverse U-shape land share curve is found for dairy cows and goats peaking at 9 °C and 6 °C, respectively. By examining the maximum and minimum temperature figures together, we find that the optimal temperature ranges for dairy cows and sheep are rather narrow. The land shares of dairy cows and sheep are larger under simultaneous high minimum and relatively low maximum temperature, but at extreme temperatures, almost all lands switch to beef cows. Goats exhibit flatter curves and maintain stable proportions as the minimum temperature changes, showing that goats are less sensitive to temperature. Again, Hypothesis 2 is confirmed here.
With respect to precipitation, we find a U-shape beef cow land share curve with an inflection point of 1000 mm with the share increasing when the precipitation rises or falls from that. All other species display inverse U shape curves with dairy cow land shares peaking at 900 mm, sheep land shares peaking at 1100 mm, and goat shares peaking at 1200 mm. We find that in the east, dairy cows substitute for beef cows as they have almost opposite curves. Moreover, we observed that dairy cow shares increase in warm regions with relatively high minimum temperature but decrease in cold regions with low minimum temperature or hot regions with high maximum temperature. The land share of beef cows approaches one when the minimum temperature drops and the maximum temperature rises, reflecting the higher range of tolerance of beef animals.
Figure 3 shows land shares for beef cows, goats, and sheep in the western US (as explained above grazing dairy cows are not included there). Here, we see different effects on land shares compared to the east. Namely, Figure 3a shows that beef cow shares demonstrate an inverse U shape relationship with the maximum temperature peaking at 25 °C and beef cow shares decreasing afterward. The opposite is found for goats and sheep with the threshold of around 25 °C and 27 °C, respectively, with the curve for goats being flatter. With respect to annual minimum temperature, we find a U shape curve for beef cow shares and inverse U shape curves for goat and sheep shares with the inflection point for annual minimum temperature being 2 °C The effect of precipitation shows a similar pattern. An inverse U shape curve is found for beef cow shares peaking at 1500 mm while U shape curves are observed for goats and sheep shares with inflection points of 1600 mm and 1400 mm. Thus, beef cow shares display opposite curves to other species, signifying substitution effects between beef cows and the other species.
In terms of sales value, in the east, the land share of beef cows is positively impacted by the calf sales price but is negatively impacted by the sales market values for other species. For dairy cows, we observe a positive response to milk, goat, and sheep sales value but a negative response to that for beef cows again, showing a land share substitution relationship between beef and dairy cows. Both goats and sheep shares react positively to beef sales value but negatively to milk sales value. Goat shares are not significantly affected by their own market value while sheep are positively correlated with theirs. In the west, beef cow shares are negatively affected by calf sales value but positively by goat sales value. We also observed that sheep shares are sensitive to market values but goat shares are not. To be specific, sheep react positively to sheep and beef cow sales value but negatively to that of goats. Goats are not significantly affected by the market value.
The amount of grassland in a county has a positive impact on the land share of beef and dairy cows in the east but no significant effect in the west. The possible reason for this is that grassland in the east is more valuable [49], and farmers tend to increase cow land shares when grassland is abundant. According to the report of Land Values 2021 Summary released by USDA NASS [49], grassland in the east is valued at USD 2993 per acre while that in the west is valued at USD 1209 per acre. In the east, both sheep and goat shares react negatively to grassland amount, showing that sheep and goats are more likely to be raised in counties with smaller grassland availability. In the west, goats are not significantly affected by grassland availability while sheep are positively impacted.
Most regional dummies are highly significant, suggesting the existence of region-specific effects. In terms of year-specific effects, we observe an increasing share of beef cows since the time dummy coefficients are positive and increasing as we move to more recent years. The coefficients for dairy cow shares decrease with a similar magnitude, demonstrating shrinking land shares over time perhaps due to industry trends toward larger animals with greater milk yields and an assumption of constant AUM requirements. Sheep and goat shares exhibit less unique year-specific effects due to the insignificance of coefficients.

4.2. Climate Scenarios

We projected the effects of climate change on grazing livestock land shares by decades using projections from the Coupled Model Intercomparison Project Phase 5 (CMIP5) [50]. We use four Representative Concentration Pathways (RCPs) to reflect different degrees of climate change mitigation: RCP 2.6, RCP 4.5, RCP 6.0, and RCP 8.5. As explained in IPCC [1], RCP 2.6 scenario implies the least future warming, and RCP 8.5 implied the most.
We obtained the projected temperature and precipitation by averaging gridded data from six different climate models following Cho and McCarl [18]: CanESM2, CCSM4, CESM1-CAM5, GFDL-CM3, HadGEM2-ES, and MPI-ESM-MR. The specific projections were obtained from the Technical Service Center, US Bureau of Reclamation [50]. We calculated county-level climate data by averaging grid-point values. All climate scenarios exhibit increasing temperature and precipitation (Figure 4). Temperature does not differ very much among RCPs between now and 2050, whereas there is about a 4 °C discrepancy between RCP 8.5 and RCP 2.6 by 2090. The rainfall predictions are more variable. In general, the greater the radiative forcing, the higher the temperature, and the more the precipitation.
Since we focus on climate effects, we keep economic and fixed effect variables fixed at their 2017 level so that the projections only reflect the marginal effects of climate change and grassland change. Grassland change is projected from FASOMGHG [38]. Figure 5 shows the impacts of these climate scenarios on species land use shares. The baseline is 2020. Here, we observed that dairy and beef cow shares experience substantial changes in opposite directions under all four RCP scenarios with the size of the shifts increasing as the radiative forcing is larger. Under RCP 8.5, the 2090 share for beef cows increases by 9%, while that for dairy cows decreases by 8%. The goat land share shows a downward trend, and the sheep land share increases under all scenarios with RCP 8.5 displaying the greatest changes.
Figure 6 shows regional changes in land shares under RCP 8.5 for 2090 and illustrates results relevant to Hypothesis 3. Note that in the projection dairy cows share effects are only present east of the Mississippi River. We observed significant regional similarities in the effects for all species. For example, the most significant increase in beef cow shares primarily occurs in Northwestern US with a few counties exhibiting shares proceeding up to more than 50%. Beef cow land shares in Central and Northeastern US are stable, varying within a 5% range. Dairy cows exhibit a declining land proportion almost everywhere except for counties in the lake states and the northeast. The land share changes of goats are minimal, with most counties adjusting within a 5% range. However, goat shares in the lake states and northeast region decrease significantly due to substitution with dairy cows. Mixed effects occur for sheep land shares, with decreases in most western counties and increases in Eastern US. Sheep shares decrease significantly in the northwest with beef cows increasing.

5. Conclusions

This study explores how farmers adapt grazing livestock species land use shares to climate. We estimate a fractional multinomial logit econometric model to examine this for dairy cows, beef cows, goats, and sheep. Dairy cows were only considered in the east because they are mostly raised in confinement in the west. The data show that beef cows are dominant in most counties with many land shares close to one. Meanwhile, goats and sheep typically have relatively small shares. A multinomial fractional logit is utilized for the estimation because it is able to deal with a data set with a nontrivial frequency of boundary values (0 and 1) and is free of the independence of irrelevant alternatives (IIA) problem.
We find in eastern US regions that dairy cow and sheep land shares increase under mild weather conditions with relatively high minimum temperatures, low maximum temperatures, and moderate precipitation. At extreme conditions we find beef cow shares increase often with land shares approaching one. In Western US, the beef cow shares display a nonlinear relationship with climate variables with land shares for goats and sheep substituting for land share with beef cows. Goats exhibit flatter curves in both the east and the west. Precipitation also has substantial effects on the east and west land shares of all species. We find beef cows demonstrate a U-shape curve east of the Mississippi River and an inverse U-shape curve in the west with all other species reacting oppositely.
We simulate future livestock share shifts under RCP 2.6, RCP 4.5, RCP 6.0, and RCP 8.5. The results under all four RCP scenarios exhibit similar patterns with the shifts being largest under the highest-emission scenario (RCP 8.5). The most significant change occurs with cow land use shares. The results generally show increasing shares for beef cows with that in Northwestern US increasing the most. We also observe decreasing dairy cow shares except in the lake states and northeast. Goat land use shares exhibit small shifts in most of the US except for the lake states and northeast where they are substituted with dairy cows. Sheep land use shares show spatial clustering effects with decreasing land share primarily observed in the west and increasing land share in the east.
As previous US studies focused on crops or the transition from cropland to grassland, there are no comparable US studies. However, our results can be compared with studies in South America [29], Africa [19], and Australia [33,34]. Our results are consistent with Australian studies where we both find increases in beef cows and sheep along with decreases in dairy cows. Studies in South America and Africa also found fewer dairy cows and a greater percentage of sheep and goats but a smaller share of beef cows in contrast with our findings.
There are some likely reasons for the difference in findings. First, livestock management differs. Our study considers only grazing animals while other studies include confined animals that usually are less affected by climate change because confinement allows more control over climate exposure with confined pigs and chickens increasing. Second, geography and local climate patterns play a role, with substantial portions of Africa and South America in the tropic exhibiting higher temperatures and more varied precipitation. Cattle there may suffer from greater heat stress, forcing farmers to switch to more heat-tolerant species such as sheep and goats. On the other hand, Australia and US exhibit more land in temperate zones, which may favor the cattle industry. Third, we believe that the composition of livestock species and local food preferences causes some differences as well. For example, goats account for less than 1% of livestock product values in the US, but Africa contributes 35.7% of world production, suggesting preference differences for goats and goat products. On the other hand, the beef cattle industry is the largest livestock sector in the US, with high expected profit and various breeds. Fourth, adaptation through changing breeds may be a force but we do not have data to examine that [17].
Taking these considerations together, we do find that land shares of livestock species have a nonlinear relationship with climate. Extreme climate conditions reduce livestock performance; therefore, farmers adapt by changing land use shares. Note that dairy cows are relatively vulnerable to climate change among the species we considered. The vulnerability of dairy cows should be noted by policymakers as the dairy industry contributed around 20% of livestock values.
There are several limitations to this study. First, we only explored the spatial variation of data by including year-fixed effects because only five census years were available (four for goats). Therefore, an extension would be to use more available data over a longer period of time. Second, we have performed some weather and inventory data interpolation by using values from nearby counties or other years and also some aggregation over livestock sales values. More complete data would be desirable. Third, we tried to include control variables on product sales values and geographic characteristics to estimate other systematic factors behind livestock land use share shifts but there might be important omitted variables such as carbon dioxide concentration, livestock production costs, and region-specific agricultural policies that could improve our ability to explain the shifts. For example, increased concentration of carbon dioxide may boost pasture yield or alter grass quality and partially offset the adverse effects of climate change.

Author Contributions

Conceptualization, B.A.M.; Data curation, M.W.; Formal analysis, M.W. and B.A.M.; Methodology, M.W. and B.A.M.; Software, M.W.; Validation, B.A.M.; Writing—original draft, M.W.; Writing—review & editing, M.W. and B.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities (Grant no. JBK2101049).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in USDA Quick Stats at https://quickstats.nass.usda.gov [46] (accessed on 17 November 2021), USHCN at https://www.ncei.noaa.gov/products/land-based-station/us-historical-climatology-network [47] (accessed on 17 November 2021), USDA ERS at https://www.ers.usda.gov/data-products/milk-cost-of-production-estimates/ [48] (accessed on 17 November 2021), and Downscaled CMIP3 and CMIP5 climate projections at ftp://gdo-dcp.ucllnl.org/pub/dcp/archive/cmip5/bcsd [50] (accessed on 17 November 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Estimation results for livestock land use share in the east.
Table A1. Estimation results for livestock land use share in the east.
VARIABLESDairy CowsGoatsSheep
annual_tmax0.6183 ***−1.3994 ***1.8614 ***
(0.2014)(0.3438)(0.2810)
annual_tmaxsq−0.0238 ***0.0194 ***−0.0493 ***
(0.0040)(0.0069)(0.0058)
annual_tmin0.2782 ***0.2702 ***0.3547 ***
(0.0147)(0.0371)(0.0270)
annual_tminsq0.0010−0.0077 ***0.0086 ***
(0.0013)(0.0026)(0.0019)
annual_prcp0.4886 ***0.9598 ***0.3736 ***
(0.0947)(0.1903)(0.1224)
annual_prcpsq−0.0261 ***−0.0396 ***−0.0168 ***
(0.0041)(0.0077)(0.0052)
Value_cattle−0.1377 ***0.0560 ***0.1060 ***
(0.0128)(0.0217)(0.0162)
Value_milkcow0.1070 ***−0.0051−0.0070
(0.0061)(0.0126)(0.0094)
Value_goat0.1862 ***0.2120 *0.0966
(0.0657)(0.1136)(0.0775)
Value_sheep0.2469 ***0.15280.4406 ***
(0.0547)(0.0958)(0.0677)
GrassLand0.0500−1.5575 ***−0.9170 ***
(0.0395)(0.2975)(0.2736)
Region = LS0.5977 ***0.2032 *0.3438 ***
(0.0684)(0.1220)(0.0824)
Region = NE0.7395 ***0.7070 ***0.8266 ***
(0.0504)(0.0878)(0.0603)
Region = SC−0.3843 ***0.0447−0.2985 **
(0.0608)(0.1236)(0.1518)
Region = SE−0.6774 ***0.6278 ***0.4474 ***
(0.0608)(0.1231)(0.1445)
Year = 2002−0.2696 ***−0.2121 *−0.0867
(0.0520)(0.1165)(0.0844)
Year = 2007−0.6154 ***−0.1672−0.3636 ***
(0.0719)(0.1442)(0.1171)
Year = 2012−1.0237 ***−0.3843 *−0.5037 ***
(0.1073)(0.2090)(0.1709)
Year = 2017−1.9417 ***−0.3564−0.6432 ***
(0.1320)(0.2485)(0.1991)
Constant−5.2240 *14.2166 ***−20.5704 ***
(2.7325)(4.7980)(3.6909)
Observations994199419941
Note: (1) The base case is beef cows. (2) Robust standard errors in parentheses. (3) *** p < 0.01, ** p < 0.05, * p < 0.1.
Table A2. Estimation results for on livestock land use share in the west.
Table A2. Estimation results for on livestock land use share in the west.
VARIABLESGoatsSheep
annual_tmax−0.7557 ***−1.5047 ***
(0.2080)(0.1584)
annual_tmaxsq0.0153 ***0.0277 ***
(0.0037)(0.0032)
annual_tmin0.0699 **0.0253
(0.0314)(0.0195)
annual_tminsq−0.0199 ***−0.0088 ***
(0.0019)(0.0015)
annual_prcp−0.2730 ***−0.2359 ***
(0.0590)(0.0381)
annual_prcpsq0.0086 ***0.0082 ***
(0.0027)(0.0017)
Value_cattle0.00890.0616 **
(0.0362)(0.0243)
Value_goat−0.2097−0.3525 **
(0.1843)(0.1435)
Value_sheep−0.13200.1141 **
(0.2168)(0.0577)
GrassLand−0.01480.0180 ***
(0.0097)(0.0044)
Region = PNW0.14900.1975
(0.2740)(0.1775)
Region = PSW−0.02870.8701 ***
(0.2985)(0.1895)
Region = RM−0.28950.4385 ***
(0.2384)(0.1278)
Region = SW−0.1887−0.6635 ***
(0.2492)(0.1495)
Year = 2002−0.2117−0.2974 ***
(0.1659)(0.1072)
Year = 2007−0.0803−0.4100 ***
(0.1803)(0.1203)
Year = 2012−0.3573−0.7970 ***
(0.2418)(0.1619)
Year = 2017−0.0330−0.6097 ***
(0.2218)(0.1590)
Constant7.9269 **18.6982 ***
(3.1065)(2.0474)
Observations52685268
Note: (1) The base case is beef cows. (2) Robust standard errors in parentheses. (3) *** p < 0.01, ** p < 0.05, * p < 0.1.

References

  1. Masson-Delmotte, V.; Zhai, P.; Pirani, P.; Connors, S.L.; Péan, C.; Berger, S.; Caud, N.; Chen, Y.; Goldfarb, L.; Gomis, M.I.; et al. IPCC, 2021: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2021; in press. [Google Scholar]
  2. Shukla, P.R.; Skea, J.; Calvo Buendia, E.; Masson-Delmotte, V.; Pörtner, H.O.; Roberts, D.C.; Zhai, P.; Slade, R.; Connors, S.; Van Diemen, R. IPCC, 2019: Climate Change and Land: An IPCC Special Report on Climate Change, Desertification, Land Degradation, Sustainable Land Management, Food Security, and Greenhouse Gas Fluxes in Terrestrial Ecosystems; Intergovernmental Panel on Climate Change (IPCC): Geneva, Switzerland, 2019. [Google Scholar]
  3. Mendelsohn, R.O.; Nordhaus, W.D.; Shaw, D. The Impact of Global Warming on Agriculture: A Ricardian Analysis. Am. Econ. Rev. 1994, 84, 753–771. [Google Scholar]
  4. USDA. 2017 Census of Agriculture Highlights; US Department of Agriculture: Washington, DC, USA, 2017.
  5. USDA. 2017 Census of Agriculture; US Department of Agriculture: Washington, DC, USA, 2019.
  6. Roenfeldt, S. You Can’t Afford to Ignore Heat Stress. Dairy Herd Manag. 1998, 35, 6–12. [Google Scholar]
  7. Hristov, A.N.; DeGaetano, A.T.; Rotz, C.A.; Hoberg, E.; Skinner, R.H.; Felix, T.; Li, H.; Patterson, P.H.; Roth, G.; Hall, M.; et al. Climate change effects on livestock in the Northeast US and strategies for adaptation. Clim. Chang. 2017, 146, 33–45. [Google Scholar] [CrossRef] [Green Version]
  8. St-Pierre, N.; Cobanov, B.; Schnitkey, G. Economic Losses from Heat Stress by US Livestock Industries. J. Dairy Sci. 2003, 86, E52–E77. [Google Scholar] [CrossRef] [Green Version]
  9. Nienaber, J.; Hahn, G. Livestock production system management responses to thermal challenges. Int. J. Biometeorol. 2007, 52, 149–157. [Google Scholar] [CrossRef] [Green Version]
  10. Sejian, V.; Maurya, V.P.; Kumar, K.; Naqvi, S.M.K. Effect of multiple stresses on growth and adaptive capability of Malpura ewes under semi-arid tropical environment. Trop. Anim. Health Prod. 2012, 45, 107–116. [Google Scholar] [CrossRef]
  11. Nardone, A.; Ronchi, B.; Lacetera, N.; Ranieri, M.; Bernabucci, U. Effects of climate changes on animal production and sustainability of livestock systems. Livest. Sci. 2010, 130, 57–69. [Google Scholar] [CrossRef]
  12. Intergovernmental Panel on Climate Change. Climate Change 2014: Impacts, Adaptation, and Vulnerability-Part B: Regional Aspects-Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Barros, V.R., Field, C.B., Dokken, D.J., Mastrandrea, M.D., Mach, K.J., Bilir, T.E., Chatterjee, M., Ebi, K.L., Estrada, Y.O., Genova, R.C., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014; ISBN 978-1-107-68386-0. [Google Scholar]
  13. Intergovernmental Panel on Climate Change. Climate Change 2001: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Third Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  14. Herrero, M.; Thornton, P.; Kruska, R.; Reid, R. Systems dynamics and the spatial distribution of methane emissions from African domestic ruminants to 2030. Agric. Ecosyst. Environ. 2008, 126, 122–137. [Google Scholar] [CrossRef]
  15. Hoffmann, I. Climate change and the characterization, breeding and conservation of animal genetic resources. Anim. Genet. 2010, 41, 32–46. [Google Scholar] [CrossRef]
  16. Mu, J.E.; McCarl, B.A.; Wein, A.M. Adaptation to climate change: Changes in farmland use and stocking rate in the U.S. Mitig. Adapt. Strat. Glob. Chang. 2013, 18, 713–730. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, Y.W.; Hagerman, A.D.; McCarl, B.A. Influence of climate factors on spatial distribution of Texas cattle breeds. Clim. Chang. 2013, 118, 183–195. [Google Scholar] [CrossRef]
  18. Cho, S.J.; McCarl, B.A. Climate change influences on crop mix shifts in the United States. Sci. Rep. 2017, 7, 40845. [Google Scholar] [CrossRef] [Green Version]
  19. Seo, S.N.; Mendelsohn, R.O. Measuring Impacts and Adaptations to Climate Change: A Structural Ricardian Model of African Livestock Management. Agric. Econ. 2008, 38, 151–165. [Google Scholar] [CrossRef]
  20. Tubiello, F.; Rosenzweig, C.; Goldberg, R.; Jagtap, S.; Jones, J. Effects of climate change on US crop production: Simulation results using two different GCM scenarios. Part I: Wheat, potato, maize, and citrus. Clim. Res. 2002, 20, 259–270. [Google Scholar] [CrossRef] [Green Version]
  21. Ovalle-Rivera, O.; Läderach, P.; Bunn, C.; Obersteiner, M.; Schroth, G. Projected Shifts in Coffea arabica Suitability among Major Global Producing Regions Due to Climate Change. PLoS ONE 2015, 10, e0124155. [Google Scholar] [CrossRef] [Green Version]
  22. Hijmans, R.J. The effect of climate change on global potato production. Am. J. Potato Res. 2003, 80, 271–279. [Google Scholar] [CrossRef]
  23. Machovina, B.; Feeley, K.J. Climate change driven shifts in the extent and location of areas suitable for export banana production. Ecol. Econ. 2013, 95, 83–95. [Google Scholar] [CrossRef]
  24. Havlík, P.; Valin, H.; Mosnier, A.; Obersteiner, M.; Baker, J.S.; Herrero, M.; Rufino, M.C.; Schmid, E. Crop Productivity and the Global Livestock Sector: Implications for Land Use Change and Greenhouse Gas Emissions. Am. J. Agric. Econ. 2012, 95, 442–448. [Google Scholar] [CrossRef]
  25. Mu, J.E.; Sleeter, B.M.; Abatzoglou, J.T.; Antle, J.M. Climate impacts on agricultural land use in the USA: The role of socio-economic scenarios. Clim. Chang. 2017, 144, 329–345. [Google Scholar] [CrossRef] [Green Version]
  26. Ghahramani, A.; Kingwell, R.; Maraseni, T.N. Land use change in Australian mixed crop-livestock systems as a transformative climate change adaptation. Agric. Syst. 2020, 180, 102791. [Google Scholar] [CrossRef]
  27. Mu, J.E.; McCarl, B.A.; Sleeter, B.; Abatzoglou, J.T.; Zhang, H. Adaptation with climate uncertainty: An examination of agricultural land use in the United States. Land Use Policy 2018, 77, 392–401. [Google Scholar] [CrossRef]
  28. Kabubo-Mariara, J. Climate change adaptation and livestock activity choices in Kenya: An economic analysis. Nat. Resour. Forum 2008, 32, 131–141. [Google Scholar] [CrossRef]
  29. Seo, S.N.; McCarl, B.A.; Mendelsohn, R. From beef cattle to sheep under global warming? An analysis of adaptation by livestock species choice in South America. Ecol. Econ. 2010, 69, 2486–2494. [Google Scholar] [CrossRef]
  30. Mader, T.L.; Frank, K.L.; Harrington, J.A., Jr.; Hahn, L.G.; Nienaber, J.A. Potential climate change effects on warm-season livestock production in the Great Plains. Clim. Chang. 2009, 97, 529–541. [Google Scholar] [CrossRef] [Green Version]
  31. Craine, J.M.; Elmore, A.; Olson, K.C.; Tolleson, D. Climate change and cattle nutritional stress. Glob. Chang. Biol. 2010, 16, 2901–2911. [Google Scholar] [CrossRef]
  32. Hahn, G.L. Dynamic responses of cattle to thermal heat loads. J. Anim. Sci. 1997, 77, 10–20. [Google Scholar] [CrossRef]
  33. Seo, S.N.; McCarl, B. Managing Livestock Species under Climate Change in Australia. Animals 2011, 1, 343–365. [Google Scholar] [CrossRef] [Green Version]
  34. Seo, S.N. Adapting to extreme climates: Raising animals in hot and arid ecosystems in Australia. Int. J. Biometeorol. 2014, 59, 541–550. [Google Scholar] [CrossRef] [PubMed]
  35. Godde, C.; Dizyee, K.; Ash, A.; Thornton, P.; Sloat, L.; Roura, E.; Henderson, B.; Herrero, M. Climate change and variability impacts on grazing herds: Insights from a system dynamics approach for semi-arid Australian rangelands. Glob. Chang. Biol. 2019, 25, 3091–3109. [Google Scholar] [CrossRef] [Green Version]
  36. Schlenker, W.; Hanemann, W.M.; Fisher, A.C. The Impact of Global Warming on U.S. Agriculture: An Econometric Analysis of Optimal Growing Conditions. Rev. Econ. Stat. 2006, 88, 113–125. [Google Scholar] [CrossRef]
  37. Redfearn, D.D.; Bidwell, T.G. Stocking Rate: The Key to Successful Livestock Production. 2006. Available online: https://extension.okstate.edu/fact-sheets/print-publications/pss/stocking-rate-the-key-to-successful-livestock-production-pss-2871.pdf (accessed on 17 November 2021).
  38. Adams, D.M.; Alig, R.J.; McCarl, B.A.; Murray, B.C. FASOMGHG Conceptual Structure, and Specification: Documentation 2005. Available online: https://docplayer.net/62907502-Fasomghg-conceptual-structure-and-specification-documentation.html (accessed on 17 November 2021).
  39. Mullahy, J. Multivariate Fractional Regression Estimation of Econometric Share Models. J. Econ. Methods 2015, 4, 71–100. [Google Scholar] [CrossRef]
  40. Murteira, J.M.R.; Ramalho, J. Regression Analysis of Multivariate Fractional Data. Econ. Rev. 2014, 35, 515–552. [Google Scholar] [CrossRef]
  41. Ramalho, E.A.; Ramalho, J.J.S.; Murteira, J.M.R. Alternative Estimating and Testing Empirical Strategies for Fractional Regression Models. J. Econ. Surv. 2011, 25, 19–68. [Google Scholar] [CrossRef]
  42. Papke, L.E.; Wooldridge, J.M. Econometric Methods for Fractional Response Variables with an Application to 401(K) Plan Participation Rates. J. Appl. Econom. 1996, 11, 619–632. [Google Scholar] [CrossRef] [Green Version]
  43. Long, J.S.; Freese, J. Regression Models for Categorical Dependent Variables Using Stata; Stata Press: College Station, TX, USA, 2006; Volume 7. [Google Scholar]
  44. Fisher, A.C.; Hanemann, W.M.; Roberts, M.J.; Schlenker, W. The Economic Impacts of Climate Change: Evidence from Agricultural Output and Random Fluctuations in Weather: Comment. Am. Econ. Rev. 2012, 102, 3749–3760. [Google Scholar] [CrossRef]
  45. USDA’s Economic Research Service USDA ERS Farm Resource Regions. Available online: https://www.ers.usda.gov/webdocs/publications/42298/32489_aib-760_002.pdf?v=42487 (accessed on 23 October 2021).
  46. US Department of Agriculture. Quick Stats 2.0; National Agricultural Statistics Service, US Department of Agriculture: Washington, DC, USA, 2018.
  47. Menne, M.J.; Durre, I.; Vose, R.S.; Gleason, B.E.; Houston, T.G. An Overview of the Global Historical Climatology Network-Daily Database. J. Atmos. Ocean. Technol. 2012, 29, 897–910. [Google Scholar] [CrossRef]
  48. USDA’s Economic Research Service Milk Cost of Production Estimates. Available online: https://www.ers.usda.gov/data-products/milk-cost-of-production-estimates/ (accessed on 21 August 2021).
  49. USDA. Land Values 2021 Summary; US Department of Agriculture: Washington, DC, USA, 2021; p. 22.
  50. US Department of the Interior, Bureau of Reclamation, Technical Services Center. Downscaled CMIP3 and CMIP5 Climate and Hydrology Projections: Release of Downscaled CMIP5 Climate Projections, Comparison with Preceding Information, and Summary of User Needs; US Department of the Interior Denver: Lakewood, CO, USA, 2013; p. 47.
Figure 1. Map of US market regions used in fixed effects. Source: Adams et al. [38].
Figure 1. Map of US market regions used in fixed effects. Source: Adams et al. [38].
Land 10 01260 g001
Figure 2. Average marginal effect of climate variables on land shares in the east: (a) Maximum Temperature; (b) Minimum Temperature; (c) Precipitation.
Figure 2. Average marginal effect of climate variables on land shares in the east: (a) Maximum Temperature; (b) Minimum Temperature; (c) Precipitation.
Land 10 01260 g002
Figure 3. Average marginal effect of climate variables on land shares in the west: (a) Maximum Temperature; (b) Minimum Temperature; (c) Precipitation.
Figure 3. Average marginal effect of climate variables on land shares in the west: (a) Maximum Temperature; (b) Minimum Temperature; (c) Precipitation.
Land 10 01260 g003
Figure 4. Projected temperature and precipitation by 2090 under climate scenarios.
Figure 4. Projected temperature and precipitation by 2090 under climate scenarios.
Land 10 01260 g004
Figure 5. Projected land share change of livestock by 2090 relative to 2020 under climate scenarios.
Figure 5. Projected land share change of livestock by 2090 relative to 2020 under climate scenarios.
Land 10 01260 g005
Figure 6. Regional change of land share of livestock in 2090 relative to 2020 under RCP 8.5: (a) Beef cows; (b) Dairy cows; (c) Goats; (d) Sheep.
Figure 6. Regional change of land share of livestock in 2090 relative to 2020 under RCP 8.5: (a) Beef cows; (b) Dairy cows; (c) Goats; (d) Sheep.
Land 10 01260 g006
Table 1. Research hypotheses.
Table 1. Research hypotheses.
HypothesesReasoning
1. The share of each livestock species has a nonlinear relationship with climate factors.Animals have optimal temperature ranges for growth and comfort known as the thermoneutral zone (TNZ). Livestock yield the highest productivity within the TNZ but may lose profitability if the climate falls above or below that zone. Thus, the share of each species is affected by climate factors.
2. Magnitudes of sensitivity to climate vary by species.Different species have different TNZs, and their land use shifts differentially under altered environmental factors. For example, dairy cows’ milk yields may drop significantly under hot and wet conditions while goats are more heat tolerant.
3. Livestock share locations will shift under climate change.With climate change, the climate of some regions may become favorable for a specific species while that of others may become likely to cause greater heat stress or cold stress. Local farmers will adapt to climate change by changing species land use shares.
Table 2. Definition of US market regions used in fixed effects.
Table 2. Definition of US market regions used in fixed effects.
RegionMarket Region States
EastCorn BeltCBIL, IN, IA, MO, OH
Lake StatesLSMI, MN, WI
NortheastNECT, DE, ME, MD, MA, NH, NJ, NY, PA, RI, VT, WV
South CentralSCAL, AR, KY, LA, MS, TN
SoutheastSEFL, GA, NC, SC, VA
WestGreat PlainsGPKS, NE, ND, SD
Pacific NorthwestPNWOR, WA
Pacific SouthwestPSWCA
Rocky MountainsRMAZ, CO, ID, MT, NV, NM, UT, WY
SouthwestSWOK, TX
Table 3. Average land share of livestock by region.
Table 3. Average land share of livestock by region.
RegionMarket RegionBeef CowDairy CowGoatSheep
EastCB75.36%19.05%1.86%3.73%
LS39.80%53.76%1.94%4.50%
NE42.29%46.03%4.07%7.61%
SC92.64%4.08%2.22%1.05%
SE83.25%8.09%5.60%3.05%
Total72.53%20.87%3.13%3.47%
WestGP94.56%0.00%1.54%3.90%
PNW89.81%0.00%2.50%7.68%
PSW84.18%0.00%3.64%12.18%
RM88.62%0.00%1.98%9.40%
SW92.55%0.00%4.04%3.41%
Total91.52%0.00%2.56%5.91%
Table 4. Descriptive statistics.
Table 4. Descriptive statistics.
VariableExplanationUnitNMeanSt. Dev.MinMax
annual_prcp3-year moving average of annual temperature and precipitation and their squares100 mm15,2409.813.480.8024.68
annual_tmax°C15,24026.703.4217.5736.52
annual_tmin°C15,240−0.954.93−14.3815.26
annual_prcpsq°C215,240108.7665.520.64609.16
annual_tmaxsq°C215,240724.34180.05308.831334.03
annual_tminsq10,000 mm215,24025.1831.670.00233.00
Value_beefcowThe market value per unit of livestockUSD15,240678.97246.35152.122080.00
Value_dairycowUSD15,2402927.30811.091602.224875.65
Value_goatUSD15,240132.0425.626.91891.35
Value_sheepUSD15,240161.6030.8356.291816.87
GrasslandGrassland available1000 acres15,240151.69335.250.026151.50
RegionRegional dummiesNA15,240NANA0.001.00
YearYear dummiesNA15,240NANA0.001.00
Table 5. Explanatory variables values by region.
Table 5. Explanatory variables values by region.
EastWest
VariableCBLSNESCSETotalGPPNWPSWRMSWTotal
annual_prcp10.388.3511.4813.4812.3711.546.419.816.153.718.396.55
annual_tmax25.3921.324.0928.5828.8826.3326.0322.5928.7725.0931.4927.4
annual_tmin−2.86−6.93−3.112.383.73−0.53−5.82−1.123.89−5.363.86−1.81
Value_cattle7.77.326.045.685.526.48.8676.977.136.827.53
Value_milkcow28.7731.3729.5823.6831.3128.5229.133.6731.4132.0330.3330.7
Value_goat1.381.461.431.191.231.311.331.441.531.341.281.34
Value_sheep1.641.661.681.471.571.591.691.761.891.741.521.66
Grassland3.582.121.724.312.923.1623.6023.3926.4763.6636.9438.37
Table 6. Average marginal effects on land share in Eastern US.
Table 6. Average marginal effects on land share in Eastern US.
VARIABLESBeef CowsDairy CowsGoatsSheep
annual_tmax−0.0729 ***0.0639 ***−0.0471 ***0.0562 ***
(0.0265)(0.0238)(0.0102)(0.0090)
annual_tmaxsq0.0031 ***−0.0025 ***0.0008 ***−0.0014 ***
(0.0005)(0.0005)(0.0002)(0.0002)
annual_tmin−0.0428 ***0.0286 ***0.0059 ***0.0083 ***
(0.0020)(0.0018)(0.0011)(0.0009)
annual_tminsq−0.00010.0001−0.0002 ***0.0003 ***
(0.0002)(0.0002)(0.0001)(0.0001)
annual_prcp−0.0807 ***0.0500 ***0.0249 ***0.0059
(0.0119)(0.0114)(0.0056)(0.0040)
annual_prcpsq0.0040 ***−0.0028 ***−0.0010 ***−0.0002
(0.0005)(0.0005)(0.0002)(0.0002)
Value_beefcow0.0110 ***−0.0182 ***0.0023 ***0.0048 ***
(0.0016)(0.0015)(0.0006)(0.0005)
Value_milkcow−0.0110 ***0.0131 ***−0.0008 **−0.0013 ***
(0.0008)(0.0007)(0.0004)(0.0003)
Value_goat−0.0264 ***0.0204 ***0.00500.0009
(0.0091)(0.0075)(0.0032)(0.0022)
Value_sheep−0.0387 ***0.0245 ***0.00250.0116 ***
(0.0073)(0.0064)(0.0027)(0.0020)
Grassland0.0487 ***0.0248 ***−0.0449 ***−0.0285 ***
(0.0117)(0.0061)(0.0087)(0.0087)
Region = LS−0.0912 ***0.0859 ***0.00100.0044 *
(0.0110)(0.0112)(0.0028)(0.0026)
Region = NE−0.1327 ***0.0989 ***0.0136 ***0.0201 ***
(0.0081)(0.0082)(0.0026)(0.0024)
Region = SC0.0476 ***−0.0455 ***0.0030−0.0052
(0.0082)(0.0074)(0.0029)(0.0036)
Region = SE0.0395 ***−0.0827 ***0.0227 ***0.0205 ***
(0.0093)(0.0066)(0.0043)(0.0060)
Year = 20020.0388 ***−0.0358 ***−0.00400.0009
(0.0076)(0.0075)(0.0034)(0.0028)
Year = 20070.0825 ***−0.0783 ***−0.0002−0.0040
(0.0100)(0.0101)(0.0042)(0.0037)
Year = 20120.1315 ***−0.1235 ***−0.0039−0.0041
(0.0138)(0.0137)(0.0058)(0.0053)
Year = 20170.2059 ***−0.2051 ***0.0010−0.0017
(0.0148)(0.0136)(0.0073)(0.0061)
Observations9941994199419941
Note: (1) Average marginal effects are calculated from coefficients of multinomial logit model. (2) Heteroskedasticity-consistent robust standard errors [42] in parentheses. (3) *** p < 0.01, ** p < 0.05, * p < 0.1.
Table 7. Average marginal effect on land share in Western US.
Table 7. Average marginal effect on land share in Western US.
VARIABLESBeef CowsGoatsSheep
annual_tmax0.0952 ***−0.0159 ***−0.0794 ***
(0.0113)(0.0049)(0.0086)
annual_tmaxsq−0.0018 ***0.0003 ***0.0015 ***
(0.0002)(0.0001)(0.0002)
annual_tmin−0.0029 **0.0017 **0.0012
(0.0015)(0.0007)(0.0010)
annual_tminsq0.0009 ***−0.0005 ***−0.0004 ***
(0.0001)(0.0001)(0.0001)
annual_prcp0.0185 ***−0.0063 ***−0.0122 ***
(0.0029)(0.0015)(0.0020)
annual_prcpsq−0.0006 ***0.0002 ***0.0004 ***
(0.0001)(0.0001)(0.0001)
Value_beefcow−0.0034 *0.00010.0033 **
(0.0018)(0.0009)(0.0013)
Value_goat0.0231 **−0.0045−0.0185 **
(0.0101)(0.0044)(0.0076)
Value_sheep−0.0029−0.00340.0064 **
(0.0064)(0.0053)(0.0030)
Grassland−0.0006−0.0004 *0.0010 ***
(0.0004)(0.0002)(0.0002)
Region = PNW−0.01460.00390.0107
(0.0142)(0.0080)(0.0100)
Region = PSW−0.0619 ***−0.00340.0653 ***
(0.0195)(0.0077)(0.0159)
Region = RM−0.0201 *−0.00800.0281 ***
(0.0106)(0.0060)(0.0077)
Region = SW0.0296 ***−0.0039−0.0257 ***
(0.0100)(0.0066)(0.0061)
Year = 20020.0242 **−0.0046−0.0196 ***
(0.0095)(0.0041)(0.0074)
Year = 20070.0275 ***−0.0012−0.0263 ***
(0.0105)(0.0047)(0.0082)
Year = 20120.0510 ***−0.0070−0.0440 ***
(0.0125)(0.0057)(0.0096)
Year = 20170.0360 ***0.0005−0.0365 ***
(0.0129)(0.0059)(0.0098)
Observations526852685268
Note: (1) Average marginal effects are calculated from coefficients of multinomial logit model. (2) Heteroskedasticity-consistent robust standard errors [42] in parentheses. (3) *** p < 0.01, ** p < 0.05, * p < 0.1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, M.; McCarl, B.A. Impacts of Climate Change on Livestock Location in the US: A Statistical Analysis. Land 2021, 10, 1260. https://0-doi-org.brum.beds.ac.uk/10.3390/land10111260

AMA Style

Wang M, McCarl BA. Impacts of Climate Change on Livestock Location in the US: A Statistical Analysis. Land. 2021; 10(11):1260. https://0-doi-org.brum.beds.ac.uk/10.3390/land10111260

Chicago/Turabian Style

Wang, Minglu, and Bruce A. McCarl. 2021. "Impacts of Climate Change on Livestock Location in the US: A Statistical Analysis" Land 10, no. 11: 1260. https://0-doi-org.brum.beds.ac.uk/10.3390/land10111260

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop