A major issue in forecasting the economic effects of climate change is the extent to which people can adapt to different weather patterns in the long term. If weather affects production or investment, adaptation can take several forms: adjusting the mix of inputs, choosing alternative technologies, or even developing a new technology that is more suited to the new climate. In addition, long-term changes in weather may lead to ecological changes, which might be amplified by human action, especially if property rights are not well-established (Diamond Reference Diamond1997; Libecap Reference Libecap2007). In any case, understanding the scope of long-term adjustment is key for measuring the impact of climate change.
While early approaches to this question relied on simulations of agronomic models, recent econometric work has used hedonic/cross-sectional or fixed-effect/panel methods. Hedonic models account for the price of land across areas with various factors, including climate variables. The estimated coefficients are interpreted as the discounted marginal value of these factors to land productivity. Important studies were conducted by Robert Mendelsohn, William Nordhaus, and Daigee Shaw (Reference Mendelsohn, Nordhaus and Shaw1994) and Wolfram Schlenker, Michael Hanemann, and Anthony Fisher (Reference Schlenker, Michael Hanemann and Fisher2005, Reference Schlenker, Michael Hanemann and Fisher2006). If areas have had enough time to adapt to permanent features of their climate, then hedonic estimates can measure the long-term effects of climate. One limitation of such cross-sectional comparisons, however, is that they measure the effect of climate on land value plus the value of any omitted variables that happen to be correlated with the measures of climate employed. Moreover, the hedonic approach cannot provide direct evidence of adaptation: hedonic methods do not actually observe adaptation, but rather assume it. A complementary approach to this question uses panel and fixed-effect models that consider the impact of cross-year weather variations within an area (for example, Deschênes and Greenstone Reference Deschênes and Greenstone2007; Schlenker and Roberts Reference Schlenker and Roberts2009; Deschênes, Greenstone, and Guryan Reference Deschênes, Greenstone and Guryan2009; Guitieras 2009; Dell, Jones, and Olken Reference Dell, Jones and Olken2012; Burke, Hsiang, and Miguel Reference Burke, Hsiang and Miguel2015). Both methods have advantages but are limited by the range of technology during their sample window and, especially for panel methods, by the scope for adjustment over the time range implicit in their research design. Meaningful adaptation occurs over long periods of time; thus, direct evidence of adaptation requires observation over longer periods.
This study adopts aspects of both methodologies but from a longer-term perspective; we examine U.S. historical data spanning about 140 years. We investigate geospatial data on average weather conditions as well as historical weather information from measurement stations on temperature and rainfall. These data are matched to economic variables at the county and state levels. The United States is well-suited for this analysis for two reasons: first, the country has sufficient weather variability across areas to estimate these models; second, conformable data are available on both weather and economic variables over a long stretch of time.
We first consider the effect of weather on farm productivity. In the late 1800s, the cross-section of farm value per acre by average temperature had an inverse U-shape. Below the median temperature, higher average temperature indicated more valuable farmland. Above the median, however, a warmer climate indicated lower farmland values. On the other hand, this negative effect for heat has attenuated, if not reversed, in recent decades. Specifically, in the second half of the twentieth century, farm value per acre increased monotonically with average temperature. A similar pattern is seen for rainfall: above the median, more average rain predicted lower farm values in the nineteenth century but higher values by the middle of the twentieth century.
Next, we use panel/fixed-effect methods to estimate the response of farm productivity to variations in temperature and rainfall. We match a long-term weather variable (measured at the decade level) with farm value, and short-term weather variables (measured by annual averages in each census year) with farm output value. The response of both farm-productivity measures to high temperature or rainfall is estimated to be more negative in the early parts of the sample than in recent years. These results suggest a substantial economic effect of weather at various horizons but also point to a change in both the short- and long-term relationships between weather and farm productivity, particularly in warm and wet areas.
In the next section, we demonstrate instability over time in the relationship between weather in early life and later-life human capital, measured by occupation-based adult income, using a census sample of cohorts born from 1860 to 1960. For those born prior to c.1920, being born in warmer or wetter areas was associated with substantially lower occupational income score in adulthood, but this correlation disappears among those born more recently. Similarly, using a fixed-effect estimator, we find effects on adult occupational income of interannual weather fluctuations in early childhood. However, this effect was strongest among those born earlier in the sample and is sharply attenuated among those born in the past 75 years. Once again, we find relationships between weather and occupational income that were strong in the past but are now considerably weaker.
The purpose of this study is to assess instability in the weather–productivity relationship over time, which measures adaptation to climate. We therefore do not characterize specific determinants of productivity change such as government policies and investments, market access, local infrastructure, and so on. Such analyses are quite limited due to a lack of historical variables and endogeneity problems. More importantly, those determinants could be part of adaptation to climate, and controlling for them would obscure the measurement of historical adaptation.
What is the relevance of these results to more contemporary issues? The timescale of anthropogenic climate change is best measured in centuries, with the greatest impacts occurring in the far future. This study suggests a simple thought experiment. Suppose we were in 1910 but were nevertheless charged with forecasting the impacts of climate change. We show that estimating standard models with circa 1900 data yields results quite different from those obtained from estimating such models with circa 2000 data. Therefore, estimates from the turn of one century would be poor guides to forecasting climate-change impacts in the century that follows. However, the recent literature is doing something similar by using estimates from the late twentieth century to forecast the systematic response of economic variables to weather or climate over the course of the twenty-first and twenty-second centuries. The results that follow suggest that greater caution is called for when conducting such exercises.
BACKGROUND AND RELATED LITERATURE
The effect of climate change on agriculture was extensively studied in the 1980s amid the increasing worldwide concern about global warming (Houghton, Jenkins, and Ephraums Reference Houghton, Jenkins and Ephraums1990). The early approach was based on agronomic models of production functions (Adams Reference Adams1989; Adams et al. Reference Adams, Rosenzweig and Pearl1990; inter alia). This traditional approach models a function of crop production from empirical or experimental results. Climatic components such as temperature and precipitation are used as key input variables. Researchers simulate the change in crop yield by varying the input variables. Using this production-function approach, various studies have predicted that U.S. agriculture would experience significant hardship because of global warming. This method is useful and accurate if the response of crop yield to weather is constant across space and time. However, it can be misleading if the response changes over time. For example, new technologies (including plant varieties) and economic conditions might reduce the negative effects of drought and extreme temperature on crop production—that is, adaptation.
Mendelsohn, Nordhaus, and Shaw (Reference Mendelsohn, Nordhaus and Shaw1994) emphasized the scope of adaptations and adjustments farmers ordinarily make in response to different economic and environmental conditions. For example, adaptive and profit-maximizing farmers will switch their main crop from wheat to corn as local temperature increases. Alan Olmstead and Paul Rhode (Reference Olmstead and Rhode2008, Reference Olmstead and Rhode2011) observed the converse in U.S. economic history, where farmers have adopted innovative varieties that dramatically increased the range over which certain crops can be grown. Thus, without considering adaptation or substitution, the production-function approach can overestimate the agricultural damage from climate change.
As an alternative approach, Mendelsohn, Nordhaus, and Shaw (Reference Mendelsohn, Nordhaus and Shaw1994) developed a hedonic price model. Based on David Ricardo's value theory, the hedonic approach assumes that the economic value of adaptive activities for climate change is reflected in land rent. In other words, if farms had enough time to adapt to climate change, the land rent will be equal to the net yield of the highest and best use of the land.Footnote 1 They measured land rent with farm value, which should be the present value of future rent. By running a regression of farm value on climatic, environmental, economic, and agronomic variables, they estimated the discounted marginal value or price sensitivity of these factors to land productivity. They found a significantly lower impact of global warming on U.S. agriculture than did the traditional production-function approach. Although the newly developed method and findings have influenced subsequent studies on climate change, the study does not actually observe the adaptation to climate or dynamic adjustment because it conducted cross-sectional comparisons in c.1980. Instead, it assumes adaptation.
Other limitations of the hedonic approach have been reported by several studies. William Cline (Reference Cline1996) questioned the validity of applying cross-sectional analyses of current land values to predict future global-warming impacts because the relative price of grain will change due to the global reduction in grain yield. Roy Darwin (Reference Darwin1999) argued that the model produced estimates that violate basic agricultural principles, as it insufficiently captured the effect of irrigation on farmland value. Similarly, Schlenker, Hanemann, and Fisher (Reference Schlenker, Michael Hanemann and Fisher2005) showed that the economic effects of climate change on agriculture must be assessed differently in dry land and irrigated areas, emphasizing the role of irrigation and water supply in the hedonic approach.
The validity of the hedonic regression was carefully examined by Olivier Deschênes and Michael Greenstone (Reference Deschênes and Greenstone2007), who recognized that unobserved characteristics are spatially correlated, such that the standard ordinary least squares (OLS) formulas used in the hedonic regression can be incorrect. In this case, the hedonic approach can confound climate with other regional factors that happen to be correlated with the measures of climate employed. As a possible solution to the omitted variable bias issue, Deschênes and Greenstone (Reference Deschênes and Greenstone2007) proposed a fixed-effects model that uses the across-year variation in temperature and precipitation within counties to estimate the effect of climate change on agricultural profits and yields. More recently, Schlenker and Michael Roberts (Reference Schlenker and Roberts2009) employed a large panel of crop yields and daily weather variables spanning most U.S. counties from 1950 to 2005. Using a fixed-effects model, they estimated the temperature thresholds above which crop yields sharply decline. The findings suggest that global warming might severely damage U.S. crop yields within the century, contrary to the prediction of Deschênes and Greenstone (Reference Deschênes and Greenstone2007).
Each of those previous studies has an advantage for understanding the effect of climate change on agriculture. One advantage of the hedonic approach is that it can account for adaptive behavior (for example, crop choice and developing new agricultural technologies), whereas the fixed-effects model can resolve the omitted variable bias problem. Nevertheless, the significant unknown in existing studies is the scope of adaptation. As noted earlier, this is not actually observed but is assumed in the hedonic approach. The studies based on the fixed-effects model are limited by the time window employed in their research design.
Meaningful adaptation takes place over a long time span; therefore, direct evidence of adaptation requires observation over longer periods of time. One strength of this study is that it harnesses the advantages of hedonic and fixed-effects methods. However, its greatest strength is that it extends the time window to span 140 years, throughout which humans have created various technological adaptations to the weather.
Relatedly, Olmstead and Rhode (Reference Olmstead and Rhode2011) provided evidence of specific adaptations by examining how North American wheat farmers overcame significant climatic challenges over the past two centuries. Richard Hornbeck (Reference Hornbeck2012) examined how agricultural adjustment occurred in the counties severely eroded by the 1930s American Dust Bowl throughout the Depression to the 1950s. This study investigates general patterns of adaptation to climate instead; therefore, their work is complementary to ours.
Finally, a few studies have examined the effect of weather in terms of non-agricultural outcomes. Most have focused on identifying the causal effect of climatic variables on outcomes such as mortality, health, educational attainment, and income (Deschênes, Greenstone, and Guryan Reference Deschênes, Greenstone and Guryan2009; Maccini and Yang Reference Maccini and Yang2009), which we will discuss at greater length in a later section. Recently, Alan Barreca et al. (Reference Barreca, Clay and Deschênes2015) estimated that the negative impact of extreme heat on mortality significantly attenuated from 1900 to 2004 by employing cross-sectional and panel methods. They argued that existing studies have ignored long-term adaptation to changes in climatic variables that evolve over decades or more and emphasize that more adaptation becomes available in the future.
DATA
This study merges historical weather variables with measures of farm productivity at the county level and of income (proxied by occupational scores) at the individual level. This section briefly explains how each variable is constructed or obtained. Detailed information on data sources and related issues is provided in Appendix 1.
First, the construction of historical county weather variables is based on the Nineteenth-Century U.S. Climate Data Set Project developed by the National Climate Data Center and the Long-Term Daily and Monthly Climate Records from the Stations across the Contiguous United States provided by the U.S. Historical Climatology Network. Both datasets report monthly mean temperature and monthly accumulated precipitation according to thousands of weather stations that have existed over the past two centuries. However, weather information is not available for all counties. Weather variables for counties without historical weather stations in certain periods need to be estimated. We use a geostatistical technique called “Kriging,” an algorithm that imputes the climatic value at a given (target) location as a weighted sum of the data values at surrounding locations. Kriging assigns weights with inverse distance squared among locations; it assigns a lower weight to surrounding weather stations located at a greater distance from the target location (Stein Reference Stein1999). Using this method, we estimate monthly mean temperatures and monthly accumulated precipitation by county and for every month from 1860 to 2000.Footnote 2 Then, we utilize the monthly estimation results to calculate annual, decadal, and centurial mean temperature and precipitation at the county- or state-level according to the research designs in later sections.
Second, this study uses two measures of county farm productivity found in Historical, Demographic, Economic, and Social Data: The United States, 1790–2002 (Haines and ICPSR Reference Haines2010): farm value per farmland acre and farm output value per farmland acre. Farm value is defined as the value of all farmland, housing, and outbuildings at the time of census enumeration. Although it is available from the 1850 agricultural census, we use the variable only from 1870 to 2000, considering the periods when the quality of estimated weather variables was acceptable. Farm output value is the total value of all farm products within the year prior to the enumeration day.Footnote 3 However, the variable is not found in the earlier source for the census years of 1910, 1920, and 1930.
Third, the changes of county boundaries matter when county fixed effects are considered. To partially fix the problem, we adjust these county variables from 1880 to 2000 on the 1870 county boundary, using the area-weighted average method.Footnote 4 We use every 10 years from 1870 on, for convenience, to refer to each census year. For some years, the actual year when each census was conducted is different from the referred year in this article. For example, the 2000 agricultural census was surveyed in 1997.
Fourth, as a measure of individual income, we use the occupational income score, which is available for a large number of censuses. The variable represents the median total income by disaggregated occupational categories that were calibrated using data from the 1950 census. Occupation has been recorded by the census for more than a century; thus, the income proxy is available for a substantial stretch of cohorts. We use micro samples from the 1880 and 1900–1990 IPUMS-Integrated Public Use Microdata Series (Ruggles et al. Reference Ruggles, Genadek and Goeken2015) and then select adult white males aged 20 to 65 who were born between 1860 and 1960. We choose 1860 as the earliest year of birth in view of the low quality of the weather data in earlier years. The units of observation of the unbalanced panel are cohorts, defined by year of birth and state of birth.
WEATHER AND FARM PRODUCTIVITY
Long-Term Weather Conditions and County Farm Value
We first investigate the historical pattern of the relationships between the county average of farm value and its weather conditions using scatter plots. This cross-sectional analysis is intended to show that county farm value has adapted to the permanent features of local weather over time in the United States. Although the range of temperature and precipitation is substantial in the United States, the results that follow do not specifically address the over-time instability of such relationships outside the United States.
For the six selected decades, the upper panel of Figure 1 plots the logarithm of each county's average farm value per acre against its decadal average of annual mean temperature (Fahrenheit), which is calculated for the 10 years prior to the year of census.Footnote 5 We also put a linear trend line and a lowess fit curve on each decade's scatter plot. The bivariate diagrams suggest that the relevance of temperature to farm value gradually changed across decades. High temperatures depressed farm value in the second half of the nineteenth century, but this negative association weakened throughout the early twentieth century. In recent decades, farm value has been (marginally) increasing with temperature. In addition, a quadratic relationship is evident in the late nineteenth century, with a peak of around 50oF. However, this inverse U-shape relationship has attenuated over time and appears non-significant in recent decades.
A similar historical pattern is observed between county farm value and decadal averages of annual accumulated precipitation, which is presented in the lower panel of Figure 1. An inverse U-shape relationship is observed until the early twentieth century. For the counties whose decadal average of annual accumulated precipitation was below about 35 inches, farm value increased with precipitation. However, a high level of precipitation above this threshold depressed farm value. The nonlinear pattern indicates a negative correlation between county farm value and precipitation in the late nineteenth century, and the negative slope became flat and non-significant throughout the early twentieth century. The attenuation of the negative relationship and inverse U-shape relationship accelerated toward the late twentieth century; county farm value has been rapidly increasing with precipitation levels in recent decades.
This bivariate cross-sectional analysis shows that (a) adaptation to hot and wet weather can be clearly observed if we look at adaptation across longer time periods and that (b) the pattern of adaptation has occurred very gradually. However, this hedonic approach does not determine whether this adaptation to permanent weather conditions resulted from a change in the climate effect or a change in the effect of other local agricultural, socioeconomic, or demographic characteristics that happen to be correlated with weather variables.
Fixed-Effects Model
As a complementary approach, we use a fixed-effects model and compare it with the cross-sectional results described earlier. Equation (1) estimates the historical pattern of the effect of long-term weather on farm value, controlling for state- and year-specific improvements in agricultural technologies and time-invariant county factors:Footnote 6
where Yijt is the logarithm of county i (in state j)'s average farm value per acre measured at census year t, Wijt is the decadal average of each weather variable (temperature or precipitation) in county i, Xijt are vectors of county-level agricultural controls and constant terms, Dt are dummies indicating year t, δt are year fixed effects, δjt are state-by-year fixed effects, and δi are county fixed effects. The latter two sets of fixed effects appear only in some specifications, as noted later. As county-level controls, we use county population density, the ratio of white population to county population, the ratio of farmland to total available county area, and the estimated number of farmers per farmland acre. Sources and the construction of each control variable are discussed in Appendix 2. Equation (1) considers the impact of cross-decade variation in weather within a county and focuses only on the linear relationship between county farm value and weather variables. In the estimation model, the reference year is 1870. The coefficients βt estimate the extent to which the response of farm value at census year t to local temperature or precipitation differed from its level in 1870. Finally, the equation is estimated using weighted ordinary least squares regressions. We calculated standard errors of estimated weather values using the Kriging process and used the inverse of the standard errors as regression weights.
Each panel in Figure 2 graphically presents the coefficient βt for temperature and precipitation. Dotted lines represent the results of estimating equation (1) including only year fixed effects from three types of fixed effects; dashed lines represent the results with additional control for state-by-year fixed effects; solid lines represent the estimation results with county fixed effects.Footnote 7
Figure 2 shows that the effect of variation in long-term temperature and rainfall on farm value rapidly changed until the mid-twentieth century. The increasing trend of βt over the years implies that high temperatures and rainfall became more beneficial to farm value as time passed. Then, βt remains stable after the period between 1960 and 1970, suggesting that adaptation to hot and rainy weather slowed in recent decades. The estimated coefficients depend on the use of fixed effects. County fixed-effects models estimate smaller effects of cross-decade variation in weather on farm value relative to 1870 levels than do the year or state-by-year fixed-effects models. Accordingly, a slower adaption to hot and rainy weather is estimated when county fixed effects are used. This means that the historical pattern observed in Figure 1 can be explained largely by non-weather unobservable county characteristics that are correlated with county weather conditions.Footnote 8
Measuring Adaptation to the Weather across Centuries
The key findings noted earlier are as follows: the effect of hot and rainy weather on farm value substantially changed over the past two centuries, and the turnover of the historical pattern occurred chiefly throughout the first half of the twentieth century. Thus, if we compare the effect between the beginning (late nineteenth century) and end (late twentieth century) periods of the sample time window, the magnitude of the adaptation can be estimated more accurately.
For the comparison, we select two periods: 1870 to 1900 and 1970 to 2000. In columns (1) and (2) of Table 1, we run weighted regressions of county farm value on decadal weather variables, the standard controls used earlier, and state dummies (in other words, state fixed effects) for each century. Column (3) estimates the significance of the difference in marginal effects between the two centuries by examining interactions between all control variables and a dummy that indicates the census year in the nineteenth century (D 19). The estimation model also controls for century-by-state fixed effects. Panels A and B estimate the marginal effect of temperature (denoted by T in Table 1) and precipitation (P) on farm value, respectively. We report only the coefficients and standard errors for weather variables in Table 1. The results of columns (1) and (2) reflect the long-term relationship between decadal weather conditions and farm value shown in Figure 1. The coefficient of the weather variables interacted with D 19 in column (3) shows that the cross-century change in the effect is statistically significant.
∗ = Significant at the 90 percent level.
∗∗ = Significant at the 95 percent level.
∗∗∗ = Significant at the 99 percent level.
Notes: We estimate the century difference in the response of county farm value to long-term weather conditions. For the regressions, we select the representative periods of two centuries: the decades ending with 1870, 1880, 1890, and 1900 for the nineteenth century, and the decades ending with 1970, 1980, 1990, and 2000 for the twentieth century. The long-term weather conditions are measured by the 10-year average of annual mean temperature (T) and annual accumulated precipitation (P) from each census year, respectively. Tm and Pm denote the sample mean of temperature and precipitation variables, respectively. D 19 denotes the dummy variable that indicates the census years in the nineteenth century (in other words, 1870, 1880, 1890, or 1900). In all the regression models, we control for each county's population density, ratio of white populations, ratio of farmland to total county area, and number of farmers per farmland acre in each census year. Columns (1) and (2) estimate the marginal effect of weather on farm value in each century. Columns (3) and (4) estimate whether the marginal effect is significantly different across the centuries by pooling both centuries and adding the interaction terms between weather variables and D 19, per equation (2). However, we use different fixed-effect models in columns (3) and (4). Weighted regressions are used in all the estimations, where the regression weight is the inverse value of the weather estimates' standard errors. Panels A and B employ temperature and precipitation, respectively; Panel C employs both variables together. The table reports only the coefficients of key variables. The standard errors of regression coefficients, reported in parentheses, are clustered on county.
Sources: Authors' calculations.
Column (4) uses state-by-year fixed effects rather than state-by-century fixed effects and also controls for county fixed effects, as follows:
where i, j, and t denote county, state, and census year, respectively. This slightly changes the marginal effect of weather in each century and the magnitude of its difference across centuries (β in equation [2]). However, the difference remains significant and similar to that seen in Figure 2.
The purpose of the analysis throughout is not to characterize the agricultural and policy determinants of farm productivity but rather to assess instability in the weather–productivity relationship over time, which we argue measures adaptation (defined broadly) over the centuries. Many endogenous factors affect farm productivity, such as government policies and investments, local infrastructure like railroads, and so on. Of course, fixed-effects models absorb all the fixed county-level and state-by-year characteristics, be they endogenous or exogenous (Figure 2 and Table 2 show similar coefficients with or without such fixed effects). However, our use of additional controls is limited by several considerations. First, few potential control variables are consistently measured across the censuses. Second, because of a lack of plausible instruments, we would not be able to correct the endogeneity problems for the policy and infrastructure variables, which are likely biased by reverse causality and omitted variables. Third and most importantly, such factors are investments in adapting to circumstances, including the climate. We would not want to remove such adaptation from the estimates, because they are part of the story.Footnote 9
∗ = Significant at the 90 percent level.
∗∗ = Significant at the 95 percent level.
∗∗∗ = Significant at the 99 percent level.
Notes: We estimate the century differences in the response of county farm output value to short-term weather conditions. The short-term weather conditions are measured by the annual mean temperature (T) and accumulated precipitation (P) in the census year. The specification in each regression is the same as in Table 2. The standard errors of regression coefficients, reported in parentheses, are clustered on county.
Sources: Authors' calculations.
The two weather variables are frequently interrelated. For example, southern counties' climates are characterized not only by high temperatures but also by wet and humid weather. The impact of one weather component on farm value can be compounded by another. We consider this issue in Panel C by running a regression for county farm value on decadal temperature and precipitation and the demeaned interaction of the two weather variables. Other specifications are the same as those used in Panels A and B across the columns. In terms of the coefficients in column (4), which uses stricter controls, the interaction effect and its change over the century appear small or less significant.
To gain a better sense of the implications of the estimated coefficients, we imagine two counties with different levels of annual mean temperature: 55.4oF for the average county and 62.5oF for a hot county, which is one standard deviation higher than the mean temperature.Footnote 10 We assume that the levels of precipitation and other agricultural conditions are identical between counties and across centuries. Then, using the coefficient of T×D 19 in column (4) of Panel C, the cross-century difference in farm value as a ratio (in other words, the level of farm value in the late twentieth century relative to that in the late nineteenth century) is estimated to be 13.7 for the average county and 19.1 for the hot county.Footnote 11 Because there is little difference in county farm value by local temperature today, this simulation suggests that adaptation occurred more rapidly in hotter areas over time. Similarly, using the coefficient of P×D 19 in column (4) of Panel C, the cross-century difference in farm value as a ratio is estimated to be 1.6 for a county with average annual precipitation (41.4 inches) and 1.9 for a wet county with 51.9 inches of annual precipitation, which is one standard deviation higher than the average precipitation. The rate of adaptation seems to be smaller than is that of temperature. However, it still suggests that adaptation to precipitation was evident in wetter areas.
Short-Term Weather Fluctuations and Farm Output Value
In this subsection, we employ an alternative measure of farm productivity: the value of county farm output. This measures how the growth of crops is affected by local weather conditions.Footnote 12 Unfortunately, farm output value is not available in the census years of 1910, 1920, and 1930. Thus, we use the panel regression described by equation (2), which compares 1870–1900 to 1970–2000, and to which we add state-by-year and county fixed effects. In addition, agricultural production depends more on weather conditions in the year of production (in other words, short-term weather conditions) than decadal weather conditions. Thus, we estimated the short-term weather variables as average temperature and accumulated precipitation in one year before each census year.Footnote 13
In Table 2, we examine the marginal effect of short-term weather conditions on county farm output value and its change over the centuries. The regression model and specification in each cell are the same as those used in Table 1. The main findings are that hot and rainy short-term weather was less favorable for producing agricultural output in the late nineteenth century than it was in the late twentieth century and that the adaptation of farm output value to weather occurred more-or-less steadily across decades.Footnote 14 Moreover, the magnitude of the adaptation to the weather (in other words, the difference in marginal effect across the centuries) is similar to what we estimated using farm value and long-term weather variables in the previous subsection. Again, these results suggest that a forecast of the impact of climate change during one century made with data from the prior century would be problematic.
Robustness
In Table 3, we test the significance of the agricultural adaptation with alternative weather variables and by scale of agriculture. First, we employ alternative weather variables in columns (2) and (3). Crops are cultivated even in winter in some regions, and most of the precipitation comes from the rainy winter season in others. However, weather conditions in winter are less significant for agriculture in many regions. To reflect this fact, we control for annual and decadal average weather variables calculated only for the months of March to October in column (2). On the other hand, we employ average weather variables calculated over the whole time horizon (1860–2000) in column (3) to eliminate the bias caused by the noisy nature of temperature and precipitation. The estimation results in both models are similar to those of the baseline estimation in column (1).
∗ = Significant at the 90 percent level.
∗∗ = Significant at the 95 percent level.
∗∗∗ = Significant at the 99 percent level.
Notes: We estimate the level of agricultural adaptation to weather, which we examined in Tables 1 and 2, using alternative weather variables, specification, or county groups, as explained in the head of columns. Panels A and B estimate the relationship between county farm value and decadal weather variables, and Panels C and D regard county farm output value and annual weather variables. The specification of each regression is the same as for column (4) of Panel A or B in Tables 1 and 2, where we adopted the baseline results in this table. We report only the coefficient of weather variables interacting with the dummy variable indicating the census years in the nineteenth century. Because column (3) uses the average weather variables calculated for the whole time horizon, the estimation of short-term weather variables in Panels C and D is unfeasible. The regression weights, except in column (4), are the inverse value of the decadal, annual or centurial weather estimates' standard errors, which is calculated from the Kriging estimation. The standard errors of regression coefficients, reported in parentheses, are clustered on county.
Sources: Authors' calculations.
Second, the extent of agriculture varies across counties and decades. More agricultural counties can be more active in adapting to weather than small-scale agricultural counties. Thus, the relationship between agricultural outcomes and weather could be influenced by counties with large-scale agriculture. To test this possibility, we conduct three additional analyses. In column (4), we run weighted regressions with the ratio of farmland to total county area in each census year as regression weights. This up-weights counties that are more agricultural. The result is similar to that of the baseline estimation. Alternatively, we restrict the regressions to counties with large-scale agriculture both in 1870 and 2000 in column (5) and those with small-scale agriculture in both census years in column (6). In both columns, we utilize the average ratio of county farmland in each census year to divide large- from small-scale agriculture. In columns (5) and (6), the adaptation is substantially estimated not only for large-scale agricultural counties but also small-scale ones. In particular, the adaptation in the relationship between annual temperature and farm output value is estimated more substantially for small-scale agricultural counties, suggesting that the adaptation to weather widely occurred across the broad scale of agriculture.
Third, some may argue that a substantial response to extreme weather conditions would be to not settle the regions with extreme conditions. In this case, many unsettled areas, particularly in earlier sample years, could skew the analysis by resulting in high farm value per farmland acre even though county land was not valuable on average. Then, the value per farmland acre could fall over time as more places become tolerable due to technological adaptation. We may not be able to completely test and reflect this possibility. Considering data limitations, one way of testing is to employ farm value or output value per county area rather than those per farmland acre.
The results of using those alternative dependent variables are reported in column (2) of Table 4. The coefficients of precipitation interacted with the nineteenth century dummy are estimated to be somewhat smaller than those of the baseline estimation in Panels B and D. This may suggest that much farmland in the late nineteenth century was unsettled or not valuable. If the pattern was more frequent in dry counties (more land-improved) than in wet counties (less land-improved), the nineteenth-century slope of agricultural outcomes per county acre against average precipitation would be flatter than as seen in Figure 1; however, this pattern seems to be less likely in terms of temperature. The estimated coefficients of temperature in column (2) are similar to those of the baseline.
∗ = Significant at the 90 percent level.
∗∗ = Significant at the 95 percent level.
∗∗∗ = Significant at the 99 percent level.
Notes: In this table, we employ alternative measures of county farm value and output value. Column (2) uses county farm value or output value per county acre rather than per county farmland acre. In column (3), to reflect the potential double counting of farm output value, we re-calculated county farm output value by excluding 50 percent of corn production and all hay production for nineteenth-century censuses and 15 percent of corn production and all hay production for twentieth-century censuses. We discounted crop values by 10 percent in column (4). Columns (3) and (4) use only the data for 1870, 1900, 1970, and 2000 due to data limitation. Other specification of estimation is the same as that of Table 3. The regression weights are the inverse value of the decadal, annual or centurial weather estimates' standard errors, which is calculated from the Kriging estimation. The standard errors of regression coefficients, reported in parentheses, are clustered on county.
Sources: Authors' calculations.
Finally, there could be a measurement error problem due to double counting in farm output value. Regarding this issue, the 1880 census volume states the following (U.S. Bureau of Census 1883, p. 26):
A large part of the corn, and a still greater portion of the hay, returned in the census are consumed for the purpose of the annual product of animal food. If the values of both the vegetable and the animal products are counted, there will be duplication to this extent…An investigation of the distribution and consumption of the corn crop of 1882, undertaken by the statistician of the department of agriculture, made the consumption for free of cattle and swine, for flesh-making or fattening purposes, 46.6 percent of the total crop.Footnote 15
The double-counting problem can be partially addressed by constructing separate crop and livestock output measures. However, it seems to be difficult to reflect the double counting problem in our estimation systematically. Calculating a county's total value of farm output was based on farmers' reports in the census. Because farmers generally report only the value of their ultimate product, we are not able to calculate the value of duplications. Handling the double-counting problem becomes more complex at the county level because feed crops produced in a county are not always consumed in the same county.
In addition, as farms are less specialized, the double-counting problem is more serious; thus, the problem would be greater for data in the nineteenth and early twentieth centuries (Gardner Reference Gardner2009). As discussed earlier, about 46.6 percent of corn production was estimated as the proportion of double counting; the ratio declined to 14.2 percent in the 1982 agricultural census and to 13.9 percent in the 1992 census.Footnote 16
We can partially test how double counting may affect our estimation. First, we recalculated county farm output values reflecting the potential double-counting problem.Footnote 17 We assume that 50 percent of corn production was used for feeding livestock in the late nineteenth century, and 15 percent in the late twentieth century, following the discussion earlier. We also assume that hay would be used only for feeding livestock. Second, we discount county crop value by 10 percent in calculating total farm output value.Footnote 18 Columns (3) and (4) in Table 4 employ recalculated farm output values as dependent variables. Their results are comparable to the baseline results, even though only data for 1870, 1900, 1970, and 2000 were used due to data limitations.
EARLY-LIFE WEATHER AND ADULT INCOME
Significance of Weather Conditions in Early Life
Many studies have revealed that climate and climate change affect activities and the lifetime health of individuals by causing malnutrition from crop failure (Ó Gráda Reference Ó Gráda2007), spreading infectious diseases (Lafferty Reference Lafferty2009; Bleakley Reference Bleakley2007; Hong Reference Hong2007, Reference Hong2011, Reference Hong2013), distorting health conditions over the course of life (WHO, WMO, and UNEP 2003; Patz et al. Reference Patz, Campbell-Lendrum and Holloway2005), changing ecological systems (Thomas et al. Reference Thomas, Cameron and Green2003), and so on. They generally predict that the overall impact of climate change like global warming would be detrimental to human life. From an economic perspective, the channels noted earlier are closely related to the reduction of individuals' economic productivity.Footnote 19
Recently, Deschênes, Greenstone, and Jonathan Guryan (Reference Deschênes, Greenstone and Guryan2009) demonstrated that extreme temperature in utero can increase the risk of having babies with low birth weights (see also Murray et al. Reference Murray, O'Reilly and Betts2000; Lawlor, Leon, and Smith Reference Lawlor, Leon and Smith2005). Although further evidence is required, they suggest that this can be caused by the association between weather in utero and fetal nutrient intake and stress. The relationship between weather in utero and birth outcome is significant because adverse birth outcomes can affect life-time human capital accumulation (for example, schooling, disability, adult income), as suggested by the fetal origins hypothesis, which has been supported by many studies (Barker Reference Barker1994; Almond Reference Almond2006; Black, Devereux, and Salvanes Reference Black, Devereux and Salvanes2007). In addition, it has been reported consistently that maternal exposure to infectious diseases promoted by weather—most notably malaria—increases the probability of having babies with low birth weights (Desai et al. Reference Desai, Ter Kuile and Nosten2007) and that early-life exposure to malaria can reduce cognitive ability and thus the acquisition of human capital (Holding and Snow Reference Holding and Snow2001; Sachs and Malaney Reference Sachs and Malaney2002; Bleakley Reference Bleakley2010; Barreca Reference Barreca2010; Hong Reference Hong2011).
Cross-Sectional Evidence
In consideration of this mechanism, our study focuses on the significance of early-life weather for later economic outcomes and relevant changes throughout U.S. history. We first look for cross-sectional evidence demonstrating the relationship between weather and adult occupational income from a longer-term perspective using simple bivariate analysis. The adult income variable used is the average occupational income score by state-of-birth cohorts aggregated from the 1880 and 1900 to 1990 IPUMS. We limited individuals to white males born from 1860 to 1960 and observed during working ages (20 to 65) in the Census IPUMS of 1880 to 1990. Using the county weather variables imputed by the Kriging interpolation technique, we construct the 1860–2000 average of annual mean temperatures and annual accumulated precipitation by state. The weather variables represent the overall weather conditions that each cohort might have experienced in early life rather than at a specific point in time.
In Figure 3, we plot average adult occupational income score against long-term average weather variables in the state of birth by two cohort groups: 1860 to 1899 and 1931 to 1960. The figure shows that adult income among the cohort born before 1900 decreases with increasing temperature and precipitation in the state of birth.Footnote 20 In other words, people who were born or had spent their early lives in hotter and wetter states attained lower paying occupations, on average, later in life. However, this relationship almost completely disappears for the cohort born after 1930; there is little difference in adult occupational income by early-life average weather conditions. (The y-axis scale is the same across cohorts for the same outcome; note the compression of the points for the later cohorts.) This suggests the possibility that technological adaptation took place throughout the early twentieth century, mitigating the adverse effects of high temperature and rainfall on individuals' economic productivity.
The Diminishing Effect of Short-Term Weather Fluctuations
In addition to the historical change in the relationship between long-term weather in early life and adult income, we find instability in the effect of short-term weather fluctuations around the year of birth on adult income using fixed-effects models. The use of long-term average weather not only captures the effect of weather around the time of birth but also measures the accumulated or lifetime effect of weather if a large proportion of each state-of-birth cohort has lived in the same state for a considerable amount of time. The use of short-term weather variables can more clearly identify the effect of early-life weather on adult income.
In the following estimation model, we search for the specific point in time around birth at which weather conditions have significant effects on adult income by running a regression predicting the cohort average occupational income score from annual mean temperature and accumulated precipitation at the year of birth, at the three years after birth, and at two years before birth, per equation (3):
where the subscripts j, k, and t denote state of birth, year of birth, and census year, respectively.
The regression is conducted for cohort samples defined by state of birth and year of birth from 1860 to 1960. The model also contains dummies for each cell of state of birth times census year (δjt ), and dummies for each cell of year of birth times census year (δkt ). These fixed effects will capture omitted and unmeasured characteristics in state of birth, year of birth, and census year that might commonly affect those in a given cohort.
We summarize the key results of the regression for equation (3) in Figure 4. The y axis in each plot displays the estimated regression coefficient on the interaction of the indicated weather variable and a dummy for the calendar year minus year of birth—thus, βl or γl ; the x axis in each plot denotes the calendar year minus the year of birth. The error bars reflect 95 percent confidence intervals for each coefficient.
The main finding is that adult occupational income scores have been highly associated with cross-year weather fluctuations for specific years in early life. Weather conditions before birth do not have meaningful effects. The year in which temperature has the most substantial, significant effect on adult occupational income is one year after birth; the effect declines over the next two years. The effect of precipitation is most substantial at the year of birth and then becomes non-significant.Footnote 21
Because the state-of-birth fixed effects capture state-specific climatic characteristics, the earlier estimation measures the effect of within-state across-year variation in weather conditions around birth. Although many studies suggest a negative effect of hot and wet weather in utero on birth outcomes, some studies have found a positive association between hot/wet weather in early life and conditions at birth through improvement in maternal nutritional status (Costa Reference Costa2004). From the perspective of the fetal origins hypothesis, the positive effect of temperature on birth outcomes can lead to a higher level of human capital accumulation and thus adult income (Maccini and Yang Reference Maccini and Yang2009; Barreca Reference Barreca2010).
Following this line of research, it is interesting that meaningful associations among weather, nutrition in utero, birth outcome, and later outcome are frequently found in studies of historical experiences and developing countries. This association might be attenuated in modern, developed countries by a sufficient year-round food supply, public programs, and technological adaptation. This suggests that the coefficients presented in Figure 4 may vary across cohorts. To investigate this possibility, we look at how the effect of short-term weather fluctuations on adult income has changed across cohorts.
We measure this long-term pattern by estimating equation (3) with a 40-year-wide moving window. The sample for each regression includes cohorts born 20 years before and after the target year; the target years are thus between 1880 and 1940. Figure 5 presents the estimated regression coefficient for the interaction between the indicated weather variable and a dummy for the year chosen according to the results in Figure 4 (in other words, one year after the year of birth for temperature and the year of birth for precipitation). We also add the 95 percent confidence intervals for each coefficient.
The result shows that the effect of short-term weather fluctuations diminished across cohorts. The effect is generally positive among cohorts born in the second half of the nineteenth century, and its across-cohort variation is substantial. However, as the target year approaches 1910 to 1920, the magnitude, significance, and across-cohort variation of the effect attenuate. Weather in early life became less relevant to adult income.
In the upper panel of Figure 5, the estimated coefficient of temperature for 1890 suggests that a one-degree increase in annual mean temperature at age 1 (from its state average) resulted in a 2.4 percent increase in adult income (measured by occupational income score) during the 1870 to 1910 period. The coefficient for 1930 implies that a one-degree increase in annual mean temperature led to a 0.47 percent increase in adult income during the 1910 to 1950 period, but the result is statistically insignificant. Similarly, in the lower panel of Figure 5, the coefficient of precipitation for 1890 suggests that a one-inch increase in annual accumulated precipitation at the year of birth (from its state average) resulted in a 0.7 percent increase in adult income during the 1870 to 1910 period. That of 1930 implies that a one-inch increase in precipitation led to a 0.12 percent increase in adult income during the 1910 to 1950 period, but the result is statistically insignificant.
The pattern of attenuation over time in Figure 5 is similar to what we observe in the cross-sectional analysis using long-term weather variables, although the role of weather in early life is estimated differently between the two analyses. Moreover, it is worthwhile to remind that the diminishing effect of weather was similarly found in the analysis of farm productivity in the previous section. The findings thus far suggest that adaptation to weather significantly occurred throughout the early twentieth century, not only in the agricultural sector but also at the individual level. Once again, we see different coefficients estimated over the course of U.S. history. In this case, the more dated coefficients drastically overstate the impact of climate and human capital.
CONCLUDING REMARKS
This study's findings suggest difficulties in forecasting the effects of climate change with either hedonic or fixed-effect models, particularly if they are estimated within a relatively narrow time window. An advantage of the cross-sectional (in other words, hedonic) comparison method is that the observed differences across areas partially reflect adaptation. However, a weakness in this approach is that such adaptation is observable only within the scope of currently available technology; this is clearly an incomplete account of adaptation. Stated otherwise, imagine being in 1910 and tasked with forecasting the effects of a warmer and wetter climate. A projection based on U.S. historical experience up to that point would vastly overstate the (negative) response by failing to account for the subsequent evolution of technology, which improved productivity in warm and wet areas. On the other hand, an advertised benefit of panel/fixed-effect methods is that they provide an upper bound to the magnitude of the response to climate change. This argument rests on the so-called “Le Chatelier principle”: more adaptation is possible in the long run; thus, sensitivity to a change in weather should lessen at longer time periods. Nevertheless, the scope of adaptation has not been sufficiently examined in existing studies, most of which examine the second half of the twentieth century.
The main contribution and novelty of this study is its extension of the time window to the past century-and-a-half, throughout which humankind has experienced substantial adaptation to the weather. This allows us to quantify the level of adaptation from various aspects (in other words, cross-sectional and panel/fixed approach, agricultural and individual productivity, and short-term and long-term adaptation). The main finding is that in the U.S. farm value and farm output value substantially declined as the region became hotter and wetter in the late nineteenth century, but the adverse effects have been almost completely attenuated in recent decades. The lesson from U.S. history is that human efforts to enhance agricultural technologies and improve ecological environments have been effective in overcoming the adverse effects of hotter and wetter weather. Although it is difficult to forecast the effects of climate change, our study suggests that the consideration of possible technological adaptations and a longer-term perspective are crucial in the forecasting.Footnote 22 At the very least, we note that today's relationships between weather or climate and economic outcomes may be poor guides to the future forms of such relationships, inasmuch as these relationships have been unstable over historical time scales.
Appendix 1: Weather Variable
Raw Data from Weather Stations
Monthly mean temperature and accumulated precipitation are obtained from two historical sources: the records of the Nineteenth-Century U.S. Climate Data Set Project developed by the National Climate Data Center and the Long-Term Daily and Monthly Climate Records from the Stations across the Contiguous United States provided by the U.S. Historical Climatology Network. The former dataset covers from 1822 to 1900 and includes 4,056 individual weather stations with different names. However, not all the stations existed throughout the nineteenth century. The number of weather stations was under 100 prior to 1854 and increased to about 400 in the 1860s, and then to about 1,000 after the 1880s. The later dataset has been utilized to obtain monthly weather information covering 1901 to 2000. The number of weather stations in the dataset is 1,221. Each weather station existed over most years of the twentieth century, so that each year's number of stations is over 1,000. Two datasets include information on weather stations' latitude and longitude, which is used to locate them on the map.
Spatial Averaging (“Kriging”)
A Kriging technique is used to interpolate monthly weather variables for the counties whose weather records are not found in the two sources noted earlier. The method is based on a linear least squares estimation algorithm, which minimizes the variance of the prediction error. Thus, a Kriging estimator is a linear combination of the values at nearby locations; the distance between neighboring points is generally employed as weights. In this study, the maximum distance between points is set as 200 miles. A sufficient number of weather stations within a 200-mile radius of the target location are required to obtain a reliable estimation result. The sizes of weather stations in the years prior to 1860 were considered insufficient; those years were thus dropped from the analysis.
We tested whether the Kriging estimation is valid by comparing estimated weather variables with actual ones. For the test, we selected 5 percent random counties out of counties with weather stations between 1890 and 1899 (60 out of 1,200 counties). We conducted the Kriging method to estimate monthly temperature and precipitation variables across counties, dropping the variables of 60 randomly selected counties. In Appendix Figure 1, we compare estimated values with actual values in terms of annual (1899) and decadal (1890–1899) weather variables for the 60 counties. Scatter plots are located around 45-degree lines. We also tested a null hypothesis that actual and estimated variables are statistically the same. The p-value for the test is reported in the figure, and the hypothesis is accepted in all cases.
Estimation of Weather Variables by County
To estimate the county-level monthly mean temperature and accumulated precipitation from 1860 to 2000, the geographical center of each county was used as the target location in the Kriging estimation, and the actual weather information of weather stations within a 200-mile buffer around the target location was utilized. The county-level decadal (long-term) or annual (short-term) average temperature is calculated by averaging monthly estimated results over the specified periods. The decadal precipitation is the average value of annual accumulated precipitation for the 10 years prior to each census year; the annual precipitation is calculated by summing up monthly precipitation over the year. Appendix Table 1 presents the mean and standard error of the estimated long- and short-term weather variables throughout the sample periods.
Appendix 2: County-Level Agricultural, Economic, Demographic, and Environmental Data
Farm Productivity
Two variables are used to measure U.S. farm productivity: farm value per farmland acre and farm output value per farmland acre. Farm value is the value of all land, housing, and outbuildings within the farm at the time of census enumeration; this can be interpreted as land price because the value of land on the farm accounts for a large portion of farm value, following hedonic studies. Farm output value is the total value of all farm products, such as crop and livestock products, within the year prior to the enumeration day. Both variables are generally obtained from Historical, Demographic, Economic, and Social Data: The United States, 1790–2002 (Haines and ICPSR Reference Haines2010), which digitized the key variables in the U.S. federal population and agricultural censuses from 1850 to 1930 and U.S. county data books from 1940 to 2000. We tried to construct the variables for the 14 decades from 1870 to 2000, but the variables of county-level farm output value are missing from the ICPSR dataset for the three decades from 1910 to 1930. All values are converted in constant 2000 U.S. dollars using the BLS GDP deflator.
Controls
The long-period analysis substantially limits the use of various controls that might account for historical changes in the relationship between climate and farm productivity. We constructed a panel of four variables at the county level from 1870 to 2000: population density, the ratio of the white population to the total county population, the ratio of farmland to the total available county area, and the estimated number of farmers per farmland acre. The first three variables were obtained from the dataset of Haines and ICPSR (Reference Haines2010) mentioned earlier. To estimate the number of farmers per farmland acre, we first calculated the county ratio of farmers from the IPUMS dataset, a census subsample, and applied the ratio to the actual population per farmland acre. Sample means of each farm-productivity measure and controls by census year are reported in Appendix Table 2.
Boundaries of Consistent Counties
Changes in county boundaries over time can make fixed-effects analyses by county problematic. We thus adjusted all the county variables from 1880 to 2000 on the 1870 county boundary by tracking the overlap of county boundaries across decades and using the (overlapped) area-weighted average method. Appendix 3 shows that the use of county boundary-adjusted data does not cause critical problems in estimating fixed-effects models.
Appendix 3: Supplementary Analyses
Fixed-Effects Analysis Using Data with the County Boundary Adjustment
Throughout this study, we use county boundary-adjusted data in estimating county fixed-effects regressions, as discussed in Appendix 2. Using boundary-adjusted data, columns (1) and (3) in Appendix Table 3 present the estimated coefficients of the interaction terms between each weather variable and each year dummy, and their standard errors, clustered on county, after regressing equation (1) with county fixed effects. These coefficients are also graphically depicted by solid lines in Figure 2. In columns (2) and (4), we conduct the same regressions but use data without the county boundary adjustment. No substantial difference is found between the estimation results with and without the adjustment.
Effect of County Boundary Change
The county boundary change adjustment may miss the fact that the sample area and data can change much over the period of this study, with different types of climate in new areas. To see whether this can affect the key findings of this study, we constructed a balanced panel including counties whose boundaries have not changed from 1870 to 2000. By analyzing the change of county boundaries on GIS, we found 654 counties out of 2,323 in 1870 that satisfy the condition and constructed a balanced panel. We report the estimation result of using the balance panel in column (2) of Appendix Table 4. The result shows that the adaptation has occurred more substantially but that the magnitude of adaptation is estimated as somewhat smaller than that of the baseline estimation in column (1), particularly in the effect on farm output value. This may be because a small number of southern counties are contained in the balanced panel (only 24 out of 654 counties).
In fact, the earlier result is very similar to that in column (3) of Appendix Table 4, where we run the same regression only with non-southern counties from the baseline sample (in other words, boundary-adjusted sample). This implies that the change of county boundaries does not critically affect the level of agricultural adaptation or the implications of this article.