Global Warming, Home Runs, and the Future of America’s Pastime

Christopher W. Callahan Program in Ecology, Evolution, Environment and Society, and Department of Geography, Dartmouth College, Hanover, New Hampshire;

Search for other papers by Christopher W. Callahan in
Current site
Google Scholar
PubMed
Close
,
Nathaniel J. Dominy Program in Ecology, Evolution, Environment and Society, and Department of Anthropology, Dartmouth College, Hanover, New Hampshire;

Search for other papers by Nathaniel J. Dominy in
Current site
Google Scholar
PubMed
Close
,
Jeremy M. DeSilva Program in Ecology, Evolution, Environment and Society, and Department of Anthropology, Dartmouth College, Hanover, New Hampshire;

Search for other papers by Jeremy M. DeSilva in
Current site
Google Scholar
PubMed
Close
, and
Justin S. Mankin Program in Ecology, Evolution, Environment and Society, and Department of Geography, and Department of Earth Sciences, Dartmouth College, Hanover, New Hampshire, and Ocean and Climate Physics, Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York

Search for other papers by Justin S. Mankin in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

Home runs in baseball—fair balls hit out of the field of play—have risen since 1980, driving strategic shifts in gameplay. Myriad factors likely account for these trends, with some speculating that global warming has contributed via a reduction in ballpark air density. Here we use observations from 100,000 Major League Baseball games and 240,000 individual batted balls to show that higher temperatures substantially increase home runs. We isolate human-caused warming with climate models, finding that >500 home runs since 2010 are attributable to historical warming. Several hundred additional home runs per season are projected due to future warming. Adaptations such as building domes on stadiums or shifting day games to night games reduce temperature’s effects on America’s pastime. Our results highlight the myriad ways that a warmer planet will restructure our lives, livelihoods, and recreation, some quantifiable and easily adapted to, as shown here, many others, not.

© 2023 American Meteorological Society. This published article is licensed under the terms of the default AMS reuse license. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Christopher W. Callahan, christopher.w.callahan.gr@dartmouth.edu

Abstract

Home runs in baseball—fair balls hit out of the field of play—have risen since 1980, driving strategic shifts in gameplay. Myriad factors likely account for these trends, with some speculating that global warming has contributed via a reduction in ballpark air density. Here we use observations from 100,000 Major League Baseball games and 240,000 individual batted balls to show that higher temperatures substantially increase home runs. We isolate human-caused warming with climate models, finding that >500 home runs since 2010 are attributable to historical warming. Several hundred additional home runs per season are projected due to future warming. Adaptations such as building domes on stadiums or shifting day games to night games reduce temperature’s effects on America’s pastime. Our results highlight the myriad ways that a warmer planet will restructure our lives, livelihoods, and recreation, some quantifiable and easily adapted to, as shown here, many others, not.

© 2023 American Meteorological Society. This published article is licensed under the terms of the default AMS reuse license. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Christopher W. Callahan, christopher.w.callahan.gr@dartmouth.edu

Home runs are an exhilarating component of baseball, occurring, typically, when a batter hits the ball entirely out of the field of play, scoring one or more runs for their team. A single home run can alter the outcome of a game, causing elation or despair for hundreds of thousands of fans, and the powerful hitters who produce scores of home runs annually are among the most valued athletes in professional sports. Major League Baseball (MLB), the highest professional baseball organization of the United States, has witnessed a long-term increase in home runs per game since the 1980s (Fig. 1a). The recent surge from 2015 to 2019 raised concern from baseball officials, fans, and players that the dominance of home runs may undermine other compelling parts of the ballgame (Albert et al. 2018, 2019; Schreck and Sickles 2018). Increases in home runs can alter baseball strategy and player acquisition and compensation, as well as the popular culture that sustains the sport. As such, understanding the source of home run trends is an important area of focus for baseball analysts.

Fig. 1.
Fig. 1.

Observed trends in home runs, temperature, and air density during baseball games. (a) Average number of home runs in Major League Baseball games from 1962 to 2019. (b) Average daily high temperature (°C: y axis; °F: internal y axis) at baseball stadiums on game days from 1962 to 2019. (c) Average air density (kg m−3) at baseball stadiums on game days from 1962 to 2019, calculated using game-day temperature, vapor pressure, and sea level pressure. In all plots, the thin line is the actual time series and the thick line is a 10-yr running mean.

Citation: Bulletin of the American Meteorological Society 104, 5; 10.1175/BAMS-D-22-0235.1

In addition to factors such as the construction of the baseball, performance-enhancing drugs, advanced technology, analytics, and player training, climate change has been raised as a potential contributor to home run trends (Dykstra 2012; Samenow 2012). The ideal gas law tells us that air density is inversely proportional to temperature (Clapeyron 1835). Ballistics tells us that the trajectory of a batted ball is influenced by temperature via its effect on density. All else being equal, warmer air is less dense, and a batted ball will carry farther (Adair 2002; Albert et al. 2018; Bahill et al. 2009). The well-documented rise in home runs (Fig. 1a) has coincided with a long-term increase in gametime temperatures at baseball stadiums (Fig. 1b) and a resulting decrease in air density during games (Fig. 1c). These trends have fueled spirited debate among commentators and sportswriters, including controversial arguments concerning the role of global warming (Dykstra 2012; Samenow 2012), but a formal analysis linking human-caused climate change and home run totals has not been performed to date.

Observers have leveraged the trove of baseball data to highlight correlations between game-day temperatures and batted ball distances (Kraft and Skeeter 1995) or home run totals (Koch and Panorska 2013), but many confounding factors complicate a causal interpretation of these relationships. Home runs have likely risen due to advanced analytics and data on pitcher tendencies, training, and changes to the construction of baseballs (Albert et al. 2018, 2019; Schreck and Sickles 2018; Vincent 2007), all of which happen to coincide with the warming trend. The variable dimensions and elevations of ballparks may further confound associations between temperature and home runs by making home runs easier or more difficult within and across parks. Additionally, the baseball season coincides with the annual cycle of temperature, and home runs may increase with spring warming as players become more comfortable hitting after the offseason. As such, estimating a causal effect of temperature on home runs requires an identification strategy that separates temperature from these other, potentially unobservable, drivers.

Here we empirically estimate the relationship between game-day temperatures and home runs across >100,000 Major League Baseball games between 1962 and 2019 and >240,000 individual batted balls between 2015 and 2019. Our model is a Poisson regression, which estimates the number of home runs in a game as a function of the temperature during that game [appendix, Eq. (A1)]. We use fixed effects to control for the spatially varying, long-term trending, and seasonally varying factors that could confound the temperature–home run relationship (Angrist and Pischke 2008; Hsiang 2016). The result from this framework is a plausibly causal estimate of the effect of an increase in temperature on the number of home runs hit in each game (appendix). We then couple our empirical estimates with experiments from climate models that separate human-caused warming from natural climate variations. In conjunction with our empirical model relating temperature to home runs, these model experiments allow us to quantify the influence of historical climate change on home run totals. They also allow us to project how home runs may change in the future with warming.

A 1°C increase in the daily high temperature on the day of a baseball game played in a stadium without a dome increases the number of home runs in that game by 1.96% [95% confidence interval (CI): 1.5–2.4] (Fig. 2a). The effect is larger (2.4% °C−1) for games played in the early afternoon when temperatures are highest (“day games”), and smaller (1.7% °C−1) for games played in the evening when temperatures are milder (“night games”). These coefficients are outside the range of estimates produced when game-day temperatures are randomized within each park (Fig. 2b, appendix), indicating that they are unlikely to be spurious (Hsiang and Jina 2014).

Fig. 2.
Fig. 2.

Empirical relationship between game-day temperature and home runs. (a) Relationship between temperature and home runs per game, derived from a Poisson model with park, year, and day-of-year fixed effects. Red line is the estimate and gray shading shows 95% confidence intervals, centered on the sample average temperature and home run number. Model is fitted to data from parks without domes. Lower histogram shows the density of temperature observations in the sample. (b) Poisson regression coefficients showing the percent change in home runs per 1°C for day and night games in parks with and without domes. Red estimates are from the actual data; black estimates are from a model where temperature values are randomized within park groups (same park, same home runs, different temperature). (c) Comparison of results from the main data (red) and Statcast data (blue). (top) Main estimate of percent change in home runs per 1°C from all games in parks without domes. (Middle) Percentage-point change in the probability of a batted ball being a home run per 1°C using a linear probability model fitted to the Statcast data (appendix). (bottom) Change in total number of home runs per game per 1°C for both data sources, converted to comparable estimates. In (b) and (c), dots correspond to central estimates while lines span the 95% confidence intervals.

Citation: Bulletin of the American Meteorological Society 104, 5; 10.1175/BAMS-D-22-0235.1

The effects described above include all games played in open-air stadiums and in retractable-roof stadiums when the roof is open (appendix), comprising 89% of our observations. In the remaining 11% of games, which are played under closed domes, we find only small and insignificant effects of temperature, as covered games are less exposed to ambient weather (Fig. 2b, Table ES1 in the online supplement).

Home runs respond similarly to temperature for both home and visiting teams independently, for subsets of the time period of analysis, and for alternative temperature datasets (Table ES2, Fig. ES1), again indicating the relationship is robust. Excluding years after 2000, in which the most recent home run and warming surges have occurred, does not alter our results, indicating that this recent period is not driving our identification. A negative binomial regression model, which relaxes the requirement for equal mean and variance in the Poisson model, yields similar results (Table ES3). When we control for precipitation, relative humidity, or wind speed, the temperature effect remains (Table ES4). Precipitation is unlikely to have an effect since games are usually canceled in inclement weather and relative humidity has only a minor effect on air density (Fig. ES2). Consistent with expectation, we do not find these additional variables to have significant effects. We also find no significant effect of wind speed (Table ES4), potentially because we do not have high-quality within-park measurements of wind speed and direction. All standard errors are clustered at the park and year levels, which accounts for spatial and temporal autocorrelation in the processes generating home runs.

Higher temperatures may affect home runs through multiple complex pathways beyond air density, such as heat stress on pitchers (Howe and Boden 2007). To clarify the mechanism, we perform two additional analyses. First, we directly allow air density to enter the regression alongside temperature (Table ES5). Independently, both air density and temperature have significant effects with comparable magnitudes. When both are included in the model, however, the effect of temperature becomes small and not significant, while the effect of air density remains (Table ES5). These results imply that temperature’s effect on home run occurs primarily by modifying air density.

To further confirm this mechanism, we analyze data from the Statcast system of high-speed cameras installed in MLB stadiums, which provides high-resolution tracking data on every fly ball from 2015 to 2019. These data allow us to control for the launch angle and launch speed of each batted ball (Albert et al. 2018). Such controls effectively hold the skill of the pitcher and hitter constant, comparing a ball leaving the bat at the same angle and speed on a warm day versus a cool day. As such, only the air density mechanism is left to account for changes in home run probability. Our game-level model [appendix, Eq. (A1)] yields a 1.83% increase in home runs per game per 1°C over 2015–19 (Fig. 2c, top). We fit a linear probability model to the Statcast data to estimate the effect of temperature on the probability that a fly ball is a home run, controlling for launch angle and speed [appendix, Eq. (A2)]. We find a 0.16-percentage-point increase in home run probability per 1°C (CI: 0.14–0.19) (Fig. 2c, middle). We can compare the game-level and batted-ball-level estimates with dimensional analysis: multiplying the game-level estimate of 1.83% by the average number of home runs in a game (2.38) yields 0.044 additional home runs for a 1°C temperature increase (CI: 0.031–0.056). Multiplying the Statcast probability change of 0.16 percentage points per fly ball by the average number of fly balls in a game (24.98) yields a very similar estimate of 0.041 home runs (CI: 0.034–0.048) (Fig. 2c, bottom). These two estimates are similar in magnitude and statistically indistinguishable, supporting our conclusion that home runs increase with temperature due to reduced air density.

Comparing temperature differences in climate model experiments that include greenhouse gas forcing versus those that exclude it allows us to quantify changes in home runs due to human-caused warming (appendix, Fig. 3). We find that human-caused climate change decreased home runs between 1962 and 1995 and increased them thereafter (Fig. 3a). Between 2010 and 2019, global warming led to an additional 58 home runs per year on average [ensemble standard deviation (SD): 20–96] and 577 cumulatively (SD: 195–959). The high atmospheric loading of anthropogenic aerosols cooled regional climate from the 1960s through the 1980s (Eyring et al. 2021). Policy changes aiming to improve air quality succeeded, so greenhouse gas forcing has dominated recent climate changes, accelerating home runs. This acceleration, while statistically distinguishable, is small—an increase of some 1% relative to the total number of home runs in 2019. Other factors such as changes in the height of the stitches on the baseball appear to have been more important in driving recent home run trends (Albert et al. 2018, 2019).

Fig. 3.
Fig. 3.

Home runs due to anthropogenic climate change. (a) Change in total home runs per season due to anthropogenic climate change over the historical period (black) and across four Shared Socioeconomic Pathway (SSP) emissions scenarios: SSP1–2.6 (teal), SSP2–4.5 (blue), SSP3–7.0 (orange), and SSP5–8.5 (red). Solid lines show the multimodel ensemble mean; shading spans the mean plus or minus the ensemble standard deviation. (b) Change in total home runs per year referenced to global mean surface temperature (GMST) change when the fraction of day and night games is held at historical averages for each park (black; Table ES6) and when all games are night games (blue). Solid line shows the ensemble mean and shading spans the mean plus or minus the ensemble standard deviation for the historical day-game fraction. Dots denote decadal averages of home runs and GMST change, following Pachauri and Meyer (2014). Inset shows cumulative additional home runs at each global warming level in thousands. Bar height shows ensemble mean and error bars span the mean plus or minus the ensemble standard deviation. (c) Additional home runs per year for each individual park at global warming levels of 1.5°, 2°, 3°, and 4°C. Parks are ranked from left to right based on the number of additional home runs at 4°C. Data are pooled across SSP scenarios in (b) and (c).

Citation: Bulletin of the American Meteorological Society 104, 5; 10.1175/BAMS-D-22-0235.1

Despite the modest historical change, home runs will—absent changes in gameplay—increase with future warming (Fig. 3a). Each degree of global warming is associated with an additional ∼95 home runs per baseball season (Fig. 3b). The magnitude and spatial pattern of warming, structural features like domes, and the fraction of day games determine ballpark-level responses (Fig. 3c, Table ES6).

Warming associated with a high-emissions pathway (SSP5–8.5) results in an additional 192 home runs per year by 2050 (SD: 131–253) and an additional 467 by 2100 (SD: 337–598), almost a 10% increase relative to 2000–19. Strong climate mitigation (SSP1–2.6) holds the increase to ∼130 by 2100. In context, these increases are large: the home run “explosion” between 2015 and 2019 resulted in ∼340 additional home runs per year (Schreck and Sickles 2018), and warming can generate a home run surge that exceeds previous such surges. Mitigation reduces these effects: limiting global warming to 1.5°C instead of 2°C would avoid an additional 1,865 home runs cumulatively (SD: 961–2,770). Limiting warming to 2°C instead of 3°C would avoid an additional 3,891 (SD: 2,434–5,347) (Fig. 3b, inset).

In the absence of mitigation, MLB could adapt to limit warming’s effects on home runs. Converting all day games into night games, for example, reduces exposure to daily high temperatures. While our main analysis holds the future frequency of day and night games fixed at the historical level for each park (appendix), a season entirely comprised of night games would cut the effect of warming by 15% on average (Fig. 3b). Other structural adaptations such as adding domes to stadiums would also dampen the influence of climate change: on average, domed parks experience only 45% of the home run increase that non-domed parks do (Fig. 3c). Other unobserved adaptations such as altering pitching strategy, in-game management, or changes in bat composition may further blunt the effects of warming.

While changes in technology and player skill will undoubtedly shape the projections we show here, our results highlight that MLB will need to contend with climate change’s influence on baseball performance. Steadily rising home runs may alter the incentives for player acquisition, offensive and defensive strategy, and public perception and engagement with the game, with consequences for the business of baseball and its on-field play.

Identifying, quantifying, and adapting to climate change requires data, and professional baseball is one of the most well-documented activities on the planet. While such data position us to assess warming impacts on home runs, they belie the far more difficult (and data-poor) task of assessing the acute human risks that accompany such gameplay changes, such as increased heat stress for ballpark staff, patrons, and players alike, or wider risks to the unequal cities in which these ballparks reside. More broadly, our findings are emblematic of the widespread influence anthropogenic global warming has already had on all aspects of life. Warming will continue to burden the poorest and most vulnerable among us, altering the risks of wildfires, heatwaves, droughts, and tropical cyclones (IPCC 2022). Our results point to the reality that even the elite billion-dollar sports industry is vulnerable to unexpected impacts. Greenhouse gas mitigation and climate adaptation are a priority not only to reduce the large-scale loss and damage associated with extreme climate events, but also to avoid pervasive (and sometimes subtle) changes to recreation and leisure activities enjoyed by people.

Acknowledgments.

This work was supported by National Science Foundation Graduate Research Fellowship 1840344 to C.W.C. and grants from Dartmouth’s Neukom Computational Institute and the Nelson A. Rockefeller Center to J.S.M. Author contributions: C.W.C. conceived the study and performed the analysis; C.W.C. and J.S.M. designed the analysis and interpreted the results; N.J.D. and J.M.D. provided insights to aid the analysis and interpretation of results; C.W.C. and J.S.M. wrote the paper with all authors providing feedback and edits. The authors declare no competing interests. The authors received no funding or endorsement from Major League Baseball or any affiliated organization.

Data availability statement.

All data and code that support the findings of this study are available at github.com/ccallahan45/Callahan-et-al_ClimateBaseball_2023.

Appendix: Materials and methods

Data.

The data used for the empirical analysis are a panel dataset of the number of home runs in each MLB game matched with daily high temperatures on the day of the game, made up of 114,417 observations spanning 1962–2019.

Our primary source of baseball data is the Retrosheet game logs database (www.retrosheet.org/gamelogs/index.html), which provides summary data for all MLB games played, reaching back to 1871. We use game log data starting in 1962, the first year in which the current 162-game season was adopted, and continuing through the end of the 2019 season. Data from 2020 or 2021 are excluded due to the unusual conditions induced by the COVID-19 pandemic. For each game, we extract the total number of home runs as well as identifying information including date, park, and the names of the home and visiting teams. For some games, manually recorded temperature within each ballpark at game time is also available from the Retrosheet event files (www.retrosheet.org/game.htm).

The primary temperature data are drawn from the Hadley Integrated Surface Dataset (HadISD) database v3.1.2 (Dunn 2019; Dunn et al. 2012, 2016; Smith et al. 2011), which provides hourly data for weather stations across the United States. Using each baseball stadium’s latitude and longitude, we find the five closest HadISD stations and average them, weighting by the inverse distance between each station and the stadium. We extract daily maximum and mean temperature, wind speed, vapor pressure, and sea level pressure from each station. Mean temperature, vapor pressure, and sea level pressure are used to calculate air density (Fig. 1c). Alternative temperature data come from ERA5 (C3S 2017), where we select daily maximum temperatures from the grid cell closest to each park, and the Retrosheet game-time temperature data. Daily total precipitation data over 1997–2015 come from the Global Precipitation Climatology Project (Adler et al. 2017), again matched by the closest grid cell for each stadium.

Most baseball stadiums are open air, but several have retractable or (in the case of Tropicana Field) fixed domes. In many of the games played in these stadiums, the Retrosheet data list whether the dome was open or closed in each game. We estimate separate regressions for games played in the open air and games played under domes (Fig. 2b) given the different exposure to ambient temperature in each situation. We drop observations where the Retrosheet data list the dome status as “unknown,” which removes ∼17% of observations in the retractable-roof parks but does not affect the data in open parks.

We download granular data on individual batted balls from 2015 to 2019 from the Statcast web portal (https://baseballsavant.mlb.com/statcast_search). We filter the data to all batted balls classified as “line drives” and “fly balls,” with the “outcome” as one of the following: home run, single, double, triple, field out, error, and sacrifice fly. We then match these batted ball data to our original game-day temperature data from HadISD, for a total of 245,175 observations.

Climate model data are drawn from the sixth phase of the Coupled Model Intercomparison Project (CMIP6) (Eyring et al. 2016; Gillett et al. 2016; O’Neill et al. 2016). We use daily maximum temperature data from the historical, historical-natural, and several Shared Socioeconomic Pathway (SSP) experiments: SSP1–2.6, SSP2–4.5, SSP3–7.0, and SSP5–8.5. Comparison of temperatures in the historical and historical-natural experiments enable an analysis of human contributions to warming thus far, and therefore how home runs have changed due to historical human-caused warming (Gillett et al. 2016). The SSP scenarios enable an analysis of future changes in home runs due to projected warming conditioned on different assumptions about emissions mitigation (O’Neill et al. 2016). The SSP scenarios used here are the “Tier 1” experiments from the Scenario Model Intercomparison Project and span a wide range of future climate outcomes, from aggressive mitigation (SSP1–2.6) to very high emissions (SSP5–8.5) (Tebaldi et al. 2021). We analyze data from 14 models across the historical and SSP experiments, each with between 1 and 50 realizations, for a total of 403 simulations (Tables ES7–ES11). The combination of multiple models and multiple realizations from individual models allows us to sample uncertainty from both model structure and model representations of internal climate variability, while using multiple scenarios allows us to sample uncertainty from emissions trajectories.

Park-level temperatures are calculated from each climate model by finding the grid cell closest to each park’s latitude and longitude and calculating the spatial average across that grid cell and the eight grid cells surrounding it.

Regression analysis.

The central goal of our analysis is to estimate a plausibly causal effect of temperature on home runs per game. We estimate the following Poisson fixed effects model that relates total home runs in each game (HR) to the daily high temperature on the day of the game (Tx), with games indexed by park p, day of year d, and year y:
HRpdy=exp(βTxpdy+μp+γd+δy).

The coefficient of interest is β, which identifies the percent change in home runs for a 1°C change in the daily high temperature on the day of the game. Parameter μ is a park fixed effect, which accounts for time-invariant average differences between parks such as their physical dimensions, latitude, or altitude. Parameter γ is a day-of-year fixed effect, which accounts for the seasonality of home runs and temperature, capturing changes in the comfort level of batters after the offseason as temperatures warm. Parameter δ is a year fixed effect, which accounts for league-wide changes across years, such as alterations to the structure of the baseball itself and the long-term global warming trend. These fixed effect terms are equivalent to demeaning each term before identifying; a park fixed effect, for example, means that the home run and temperature terms are estimated using deviations from each park’s long-term average number of home runs and temperature, respectively. We use a Poisson model because the home run data are integer values with an approximately exponential distribution.

We use two-way clustering for the standard errors on the estimates from Eq. (A1), clustering simultaneously at the park and year level to account for autocorrelation in the errors in both time and space. When uncertainty in the coefficients is sampled in the climate model analysis (see “Historical home run attribution” and “Future home run projections” sections), we generate a set of coefficients from a Gaussian distribution using the estimate of β as the mean and the clustered standard error on β as the standard deviation. We also randomize temperatures within parks in Fig. 2b to nonparametrically estimate a null distribution of coefficients, accounting for autocorrelation in home runs within parks across time.

Because higher temperatures could affect pitcher and hitter skill through mechanisms such as heat stress, air density may not be the only pathway by which global warming might alter home runs. To isolate the role of air density, we perform two analyses. First, we estimate models with and without air density as a separate term (Table ES5), which shows that when both temperature and air density are included, only air density has a statistically significant effect. Second, to control for batted-ball characteristics, we analyze the Statcast data at the level of individual batted balls. These data include both the initial speed (sometimes called “exit velocity”) and angle of line drives and fly balls as they come off the bat. Controlling for launch speed and launch angle holds the characteristics of the batted ball constant, so that environmental variables like temperature can only affect home runs after the ball is hit.

We estimate the following linear probability model for the binary outcome that a given fly ball is a home run (P), with balls indexed by b, A representing launch angle, and S representing launch speed:
Pbpdy=β1Txpdy+β2Abpdy+β3Abpdy2+β4Sbpdy+μp+γd+δy+ϵbpdy.

The coefficient of interest is β1, which is the change in probability that a ball is a home run for a 1°C change in temperature, after controlling for 1) the fixed characteristics of the ballpark, time of year, and year [as in Eq. (A1)] and 2) the initial characteristics of each batted ball before it must fly through the air and therefore be affected by air resistance. We use a quadratic in launch angle because batted balls that are both too low (line drives) and too high (pop-ups) are unlikely to be home runs, though using only a linear term for launch angle yields similar results for β1. We cluster standard errors at the game level because the treatment (temperature) varies by game, not by individual batted ball (Abadie et al. 2017).

Historical home run attribution.

To determine the contribution of historical global warming to observed home run trends, we compare the historical and historical-natural (“natural”) simulations from CMIP6 over 1962–2019 (Table ES7). Note that because the historical simulations extend only to 2014, whereas the natural simulations extend through 2019, we splice the historical simulations with years 2015–19 from each model’s corresponding SSP2–4.5 simulation, following the CMIP6 experimental protocol (Gillett et al. 2016).

We use a “delta method” to bias-correct the simulations (Diffenbaugh et al. 2021). For each day in each simulation, we take the difference between the historical and natural simulations. We then subtract this difference from the observed temperature on that day to calculate the “counterfactual” temperature. We then apply the regression coefficients to the observed and counterfactual temperatures on game days to calculate the difference in home runs due to warming.

We use distinct coefficients for day and night games and games played with and without domes (Table ES1). For games where the status of the retractable-roof dome is “unknown” (∼17% of games played in retractable-roof stadiums, ∼3.1% of all games), we make the conservative assumption that the dome was closed on those days. As such, our estimates of historical home runs due to warming are likely a lower bound on the true number.

Future home run projections.

Future home runs are projected by comparing daily maximum temperatures from the SSP scenarios over 2020–2100 to the 2000–19 average in each park. We preserve the seasonal cycle in the climatology by calculating the average 2000–19 temperature for each day of year for each park from the HadISD data. This comparison quantifies how future home runs evolve relative to a counterfactual scenario in which temperatures remain at their early-twenty-first-century average.

We use a similar delta method as in the historical analysis. We take the difference between the model-projected temperature in each game and subtract the model-simulated climatological 2000–19 temperature for the corresponding day of year. We then add this difference to the observed 2000–19 climatology. Additional home runs due to warming are calculated by comparing the bias-corrected model-projected temperature for each game day to the corresponding observed climatological temperature and applying the regression coefficients to calculate the effect of this difference on home runs.

Projections are referenced to global mean surface temperature (GMST) change where appropriate. GMST change is calculated as the average of global surface temperature relative to 1850–1900, spatially weighted by the cosine of latitude. Where we show global warming levels (GWLs), we define them as the first year that a 30-yr running mean of GMST change crosses each level, similar to the method used by the Intergovernmental Panel on Climate Change (Lee et al. 2021). As an analytical strategy, benchmarking projections to GWLs neglects the time scale and pace of future emissions. However, the quasi-linear relationship between temperature change and future home runs (Fig. 3b) supports our strategy of calculating cumulative home runs at GWLs regardless of the specific emissions scenario.

Uncertainty in the climate model analysis.

In both the historical attribution and future projections, uncertainty is sampled with a Monte Carlo analysis. In the historical attribution analysis, we run 1,000 simulations, where each simulation randomly selects one climate model realization and one set of regression coefficients from their respective distributions.

In the case of the future projections, it is unknown how future baseball games will be distributed throughout the baseball season (approximately April–September). Eighty-one games are played in each park in each regular season, but the effect of warming may depend on the times during the year that those games are played. As such, we simulate a series of future baseball seasons, randomly distributing 81 games on dates from April through September in each of the stadiums used by the 30 MLB teams (Table ES6). These schedule simulations are performed 250 times, and in each iteration, a random set of the games are chosen as day games and the rest are chosen as night games. While the exact dates chosen as day and night games is random, the fraction of day and night games is set to equal the average 2000–19 fraction in each park (Table ES6). Similarly, for retractable-roof stadiums, the fraction of games where the dome is closed is set to equal the average 2000–19 fraction (Table ES6). These simulated schedules do not perfectly represent how MLB schedules will be set in the future, which is unknown, but they allow us to quantify the magnitude of the change in home runs driven by warming assuming similar gameplay conditions as today.

Then, for each SSP scenario, we run 250 simulations where each simulation randomly selects one set of regression coefficients from their distributions, one climate model realization, and one simulated schedule. With four scenarios, this yields 1,000 simulations in total.

In both the historical and future cases, we adjust the sampling probability for models with multiple realizations so that each model is equally likely to be sampled and models with multiple realizations do not dominate the analysis. The coefficient samples are from a Gaussian distribution (see “Regression analysis”) and the schedule samples are from a uniform distribution.

Home run definition.

As defined by Major League Baseball, a home run is any fair ball where the batter makes a complete circuit of the bases. This definition includes both out-of-play home runs (balls hit over the outfield fence, the vast majority of home runs) and inside-the-park home runs (fair balls inside the field of play where a player is able to run around all the bases before being put out). Increased temperatures and their effects on air density are likely to affect out-of-play home runs substantially more than inside-the-park home runs, since inside-the-park home runs do not necessarily need to fly several hundred feet in the air.

The Retrosheet game summary data we use does not distinguish between these two types of home runs, making it imperfect for isolating the effect of temperature on home runs through the mechanism of reduced air density. However, because inside-the-park home runs are extremely rare—less than 1% of all home runs—conflating these two types of home runs in the data should not substantially affect the analysis. Additionally, because the more granular Statcast data allow us to limit our analysis to only fly ball home runs, the estimates using those data are not subject to this limitation.

References

  • Abadie, A., S. Athey, G. W. Imbens, and J. Wooldridge, 2017: When should you adjust standard errors for clustering? NBER Working Paper Series, Working Paper 24003, 30 pp., www.nber.org/system/files/working_papers/w24003/w24003.pdf.

  • Adair, R. K., 2002: The Physics of Baseball. 3rd ed. HarperCollins, 192 pp.

  • Adler, R., J.-J. Wang, M. Sapiano, G. Huffman, D. Bolvin, E. Nelkin, and NOAA CDR Program, 2017: Global Precipitation Climatology Project (GPCP) Climate Data Record (CDR), version 1.3 (daily). NOAA National Centers for Environmental Information, accessed 18 October 2018, www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C00999.

  • Albert, J., and Coauthors, 2018: Report of the Committee Studying Home Run Rates in Major League Baseball. Office of the Commissioner of Baseball, 84 pp., http://baseball.physics.illinois.edu/HRReport2018.pdf.

  • Albert, J., A. Hosoi, A. Nathan, and L. Smith, 2019: Preliminary report of the Committee Studying Home Run Rates in MLB. Office of the Commissioner of Baseball, 27 pp., http://baseball.physics.illinois.edu/HRReport2019.pdf.

  • Angrist, J. D., and J.-S. Pischke, 2008: Mostly Harmless Econometrics: An Empiricist’s Companion. Princeton University Press, 392 pp.

  • Bahill, A. T., D. G. Baldwin, and J. S. Ramberg, 2009: Effects of altitude and atmospheric conditions on the flight of a baseball. Int. J. Sports Sci. Eng., 3, 109128.

    • Search Google Scholar
    • Export Citation
  • Clapeyron, E., 1835: Mémoire sur la puissance motrice de la chaleur. J. Ec. Polytech., XIV, 153190.

  • C3S, 2017: ERA5: Fifth generation of ECMWF atmospheric reanalysis of the global climate. Copernicus Climate Change Service Climate Data Store, accessed 26 October 2019, https://cds.climate.copernicus.eu/cdsapp.

  • Diffenbaugh, N. S., F. V. Davenport, and M. Burke, 2021: Historical warming has increased U.S. crop insurance losses. Environ. Res. Lett., 16, 084025, https://doi.org/10.1088/1748-9326/ac1223.

    • Search Google Scholar
    • Export Citation
  • Dunn, R. J. H., 2019: HadISD version 3: Monthly updates. Met Office Hadley Centre Tech. Note 103, 8 pp., https://digital.nmla.metoffice.gov.uk/digitalFile_13890750-fb6f-42c7-92df-1c4504621fae/.

  • Dunn, R. J. H., K. M. Willett, P. W. Thorne, E. V. Woolley, I. Durre, A. Dai, D. E. Parker, and R. S. Vose, 2012: HadISD: A quality-controlled global synoptic report database for selected variables at long-term stations from 1973–2011. Climate Past, 8, 16491679, https://doi.org/10.5194/cp-8-1649-2012.

    • Search Google Scholar
    • Export Citation
  • Dunn, R. J. H., K. M. Willett, D. E. Parker, and L. Mitchell, 2016: Expanding HadISD: Quality-controlled, sub-daily station data from 1931. Geosci. Instrum., Methods Data Syst., 5, 473491, https://doi.org/10.5194/gi-5-473-2016.

    • Search Google Scholar
    • Export Citation
  • Dykstra, P., 2012: Is global warming causing more home runs in baseball? Scientific American, 30 April, www.scientificamerican.com/article/global-warming-more-baseball-home-runs/.

  • Eyring, V., S. Bony, G. A. Meehl, C. A. Senior, B. Stevens, R. J. Stouffer, and K. E. Taylor, 2016: Overview of the Coupled Model Intercomparison Project phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev., 9, 19371958. https://doi.org/10.5194/gmd-9-1937-2016.

    • Search Google Scholar
    • Export Citation
  • Eyring, V., and Coauthors, 2021: Human influence on the climate system. Climate Change 2021: The Physical Science Basis, V. Masson-Delmotte et al. Eds., Cambridge University Press, 423552, www.ipcc.ch/report/ar6/wg1/chapter/chapter-3/.

    • Search Google Scholar
    • Export Citation
  • Gillett, N., and Coauthors, 2016: The Detection and Attribution Model Intercomparison Project (DAMIP v1.0) contribution to CMIP6. Geosci. Model Dev., 9, 36853697, https://doi.org/10.5194/gmd-9-3685-2016.

    • Search Google Scholar
    • Export Citation
  • Howe, A. S., and B. P. Boden, 2007: Heat-related illness in athletes. Amer. J. Sports Med., 35, 13841395, https://doi.org/10.1177/0363546507305013.

    • Search Google Scholar
    • Export Citation
  • Hsiang, S., 2016: Climate econometrics. Annu. Rev. Resour. Econ., 8, 4375, https://doi.org/10.1146/annurev-resource-100815-095343.

  • Hsiang, S. M., and A. S. Jina, 2014: The causal effect of environmental catastrophe on long-run economic growth: Evidence from 6,700 cyclones. NBER Working Paper 20352, 69 pp., https://gspp.berkeley.edu/assets/uploads/research/pdf/The_Causal_Effect_of_Environmental_Catastrophe_on_Long-Run_Economic_Growth.pdf.

  • IPCC, 2022: Climate Change 2022: Impacts, Adaptation, and Vulnerability. Cambridge University Press, 3056 pp.

  • Koch, B. L. D., and A. K. Panorska, 2013: The impact of temperature on Major League Baseball. Wea. Climate Soc., 5, 359366, https://doi.org/10.1175/WCAS-D-13-00002.1.

    • Search Google Scholar
    • Export Citation
  • Kraft, M. D., and B. R. Skeeter, 1995: The effect of meteorological conditions on fly ball distances in North American Major League Baseball games. Geogr. Bull., 37, 4048.

    • Search Google Scholar
    • Export Citation
  • Lee, J.-Y., and Coauthors, 2021: Future global climate: Scenario-based projections and near-term information. Climate Change 2021: The Physical Science Basis, V. Masson-Delmotte et al., Eds., Cambridge University Press, 553672, www.ipcc.ch/report/ar6/wg1/chapter/chapter-4/.

    • Search Google Scholar
    • Export Citation
  • O’Neill, B. C., and Coauthors, 2016: The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6. Geosci. Model Dev., 9, 34613482, https://doi.org/10.5194/gmd-9-3461-2016.

    • Search Google Scholar
    • Export Citation
  • Pachauri, R. K., and L. A. Meyer, 2014: Climatic Change 2014: Synthesis Report. IPCC, 151 pp., www.ipcc.ch/report/ar5/syr/.

  • Samenow, J., 2012: Tim McCarver isn’t crazy: The home run and global warming connection. Washington Post, 30 April, www.washingtonpost.com/blogs/capital-weather-gang/post/tim-mccarver-may-not-be-crazy-the-home-run-and-global-warming-connection/2012/04/30/gIQA1hI1rT_blog.html.

  • Schreck, L., and E. Sickles, 2018: Offensive explosion: Trends in hitting production in Major League Baseball. Baseball Res. J., 47, 1928, https://sabr.org/journal/article/offensive-explosion-trends-in-hitting-production-in-major-league-baseball/.

    • Search Google Scholar
    • Export Citation
  • Smith, A., N. Lott, and R. Vose, 2011: The integrated surface database: Recent developments and partnerships. Bull. Amer. Meteor. Soc., 92, 704708, https://doi.org/10.1175/2011BAMS3015.1.

    • Search Google Scholar
    • Export Citation
  • Tebaldi, C., and Coauthors, 2021: Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6. Earth Syst. Dyn., 12, 253293, https://doi.org/10.5194/esd-12-253-2021.

    • Search Google Scholar
    • Export Citation
  • Vincent, D., 2007: Home Run: The Definitive History of Baseball’s Ultimate Weapon. Potomac Books, 304 pp.

Supplementary Materials

Save
  • Abadie, A., S. Athey, G. W. Imbens, and J. Wooldridge, 2017: When should you adjust standard errors for clustering? NBER Working Paper Series, Working Paper 24003, 30 pp., www.nber.org/system/files/working_papers/w24003/w24003.pdf.

  • Adair, R. K., 2002: The Physics of Baseball. 3rd ed. HarperCollins, 192 pp.

  • Adler, R., J.-J. Wang, M. Sapiano, G. Huffman, D. Bolvin, E. Nelkin, and NOAA CDR Program, 2017: Global Precipitation Climatology Project (GPCP) Climate Data Record (CDR), version 1.3 (daily). NOAA National Centers for Environmental Information, accessed 18 October 2018, www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C00999.

  • Albert, J., and Coauthors, 2018: Report of the Committee Studying Home Run Rates in Major League Baseball. Office of the Commissioner of Baseball, 84 pp., http://baseball.physics.illinois.edu/HRReport2018.pdf.

  • Albert, J., A. Hosoi, A. Nathan, and L. Smith, 2019: Preliminary report of the Committee Studying Home Run Rates in MLB. Office of the Commissioner of Baseball, 27 pp., http://baseball.physics.illinois.edu/HRReport2019.pdf.

  • Angrist, J. D., and J.-S. Pischke, 2008: Mostly Harmless Econometrics: An Empiricist’s Companion. Princeton University Press, 392 pp.

  • Bahill, A. T., D. G. Baldwin, and J. S. Ramberg, 2009: Effects of altitude and atmospheric conditions on the flight of a baseball. Int. J. Sports Sci. Eng., 3, 109128.

    • Search Google Scholar
    • Export Citation
  • Clapeyron, E., 1835: Mémoire sur la puissance motrice de la chaleur. J. Ec. Polytech., XIV, 153190.

  • C3S, 2017: ERA5: Fifth generation of ECMWF atmospheric reanalysis of the global climate. Copernicus Climate Change Service Climate Data Store, accessed 26 October 2019, https://cds.climate.copernicus.eu/cdsapp.

  • Diffenbaugh, N. S., F. V. Davenport, and M. Burke, 2021: Historical warming has increased U.S. crop insurance losses. Environ. Res. Lett., 16, 084025, https://doi.org/10.1088/1748-9326/ac1223.

    • Search Google Scholar
    • Export Citation
  • Dunn, R. J. H., 2019: HadISD version 3: Monthly updates. Met Office Hadley Centre Tech. Note 103, 8 pp., https://digital.nmla.metoffice.gov.uk/digitalFile_13890750-fb6f-42c7-92df-1c4504621fae/.

  • Dunn, R. J. H., K. M. Willett, P. W. Thorne, E. V. Woolley, I. Durre, A. Dai, D. E. Parker, and R. S. Vose, 2012: HadISD: A quality-controlled global synoptic report database for selected variables at long-term stations from 1973–2011. Climate Past, 8, 16491679, https://doi.org/10.5194/cp-8-1649-2012.

    • Search Google Scholar
    • Export Citation
  • Dunn, R. J. H., K. M. Willett, D. E. Parker, and L. Mitchell, 2016: Expanding HadISD: Quality-controlled, sub-daily station data from 1931. Geosci. Instrum., Methods Data Syst., 5, 473491, https://doi.org/10.5194/gi-5-473-2016.

    • Search Google Scholar
    • Export Citation
  • Dykstra, P., 2012: Is global warming causing more home runs in baseball? Scientific American, 30 April, www.scientificamerican.com/article/global-warming-more-baseball-home-runs/.

  • Eyring, V., S. Bony, G. A. Meehl, C. A. Senior, B. Stevens, R. J. Stouffer, and K. E. Taylor, 2016: Overview of the Coupled Model Intercomparison Project phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev., 9, 19371958. https://doi.org/10.5194/gmd-9-1937-2016.

    • Search Google Scholar
    • Export Citation
  • Eyring, V., and Coauthors, 2021: Human influence on the climate system. Climate Change 2021: The Physical Science Basis, V. Masson-Delmotte et al. Eds., Cambridge University Press, 423552, www.ipcc.ch/report/ar6/wg1/chapter/chapter-3/.

    • Search Google Scholar
    • Export Citation
  • Gillett, N., and Coauthors, 2016: The Detection and Attribution Model Intercomparison Project (DAMIP v1.0) contribution to CMIP6. Geosci. Model Dev., 9, 36853697, https://doi.org/10.5194/gmd-9-3685-2016.

    • Search Google Scholar
    • Export Citation
  • Howe, A. S., and B. P. Boden, 2007: Heat-related illness in athletes. Amer. J. Sports Med., 35, 13841395, https://doi.org/10.1177/0363546507305013.

    • Search Google Scholar
    • Export Citation
  • Hsiang, S., 2016: Climate econometrics. Annu. Rev. Resour. Econ., 8, 4375, https://doi.org/10.1146/annurev-resource-100815-095343.

  • Hsiang, S. M., and A. S. Jina, 2014: The causal effect of environmental catastrophe on long-run economic growth: Evidence from 6,700 cyclones. NBER Working Paper 20352, 69 pp., https://gspp.berkeley.edu/assets/uploads/research/pdf/The_Causal_Effect_of_Environmental_Catastrophe_on_Long-Run_Economic_Growth.pdf.

  • IPCC, 2022: Climate Change 2022: Impacts, Adaptation, and Vulnerability. Cambridge University Press, 3056 pp.

  • Koch, B. L. D., and A. K. Panorska, 2013: The impact of temperature on Major League Baseball. Wea. Climate Soc., 5, 359366, https://doi.org/10.1175/WCAS-D-13-00002.1.

    • Search Google Scholar
    • Export Citation
  • Kraft, M. D., and B. R. Skeeter, 1995: The effect of meteorological conditions on fly ball distances in North American Major League Baseball games. Geogr. Bull., 37, 4048.

    • Search Google Scholar
    • Export Citation
  • Lee, J.-Y., and Coauthors, 2021: Future global climate: Scenario-based projections and near-term information. Climate Change 2021: The Physical Science Basis, V. Masson-Delmotte et al., Eds., Cambridge University Press, 553672, www.ipcc.ch/report/ar6/wg1/chapter/chapter-4/.

    • Search Google Scholar
    • Export Citation
  • O’Neill, B. C., and Coauthors, 2016: The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6. Geosci. Model Dev., 9, 34613482, https://doi.org/10.5194/gmd-9-3461-2016.

    • Search Google Scholar
    • Export Citation
  • Pachauri, R. K., and L. A. Meyer, 2014: Climatic Change 2014: Synthesis Report. IPCC, 151 pp., www.ipcc.ch/report/ar5/syr/.

  • Samenow, J., 2012: Tim McCarver isn’t crazy: The home run and global warming connection. Washington Post, 30 April, www.washingtonpost.com/blogs/capital-weather-gang/post/tim-mccarver-may-not-be-crazy-the-home-run-and-global-warming-connection/2012/04/30/gIQA1hI1rT_blog.html.

  • Schreck, L., and E. Sickles, 2018: Offensive explosion: Trends in hitting production in Major League Baseball. Baseball Res. J., 47, 1928, https://sabr.org/journal/article/offensive-explosion-trends-in-hitting-production-in-major-league-baseball/.

    • Search Google Scholar
    • Export Citation
  • Smith, A., N. Lott, and R. Vose, 2011: The integrated surface database: Recent developments and partnerships. Bull. Amer. Meteor. Soc., 92, 704708, https://doi.org/10.1175/2011BAMS3015.1.

    • Search Google Scholar
    • Export Citation
  • Tebaldi, C., and Coauthors, 2021: Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6. Earth Syst. Dyn., 12, 253293, https://doi.org/10.5194/esd-12-253-2021.

    • Search Google Scholar
    • Export Citation
  • Vincent, D., 2007: Home Run: The Definitive History of Baseball’s Ultimate Weapon. Potomac Books, 304 pp.

  • Fig. 1.

    Observed trends in home runs, temperature, and air density during baseball games. (a) Average number of home runs in Major League Baseball games from 1962 to 2019. (b) Average daily high temperature (°C: y axis; °F: internal y axis) at baseball stadiums on game days from 1962 to 2019. (c) Average air density (kg m−3) at baseball stadiums on game days from 1962 to 2019, calculated using game-day temperature, vapor pressure, and sea level pressure. In all plots, the thin line is the actual time series and the thick line is a 10-yr running mean.

  • Fig. 2.

    Empirical relationship between game-day temperature and home runs. (a) Relationship between temperature and home runs per game, derived from a Poisson model with park, year, and day-of-year fixed effects. Red line is the estimate and gray shading shows 95% confidence intervals, centered on the sample average temperature and home run number. Model is fitted to data from parks without domes. Lower histogram shows the density of temperature observations in the sample. (b) Poisson regression coefficients showing the percent change in home runs per 1°C for day and night games in parks with and without domes. Red estimates are from the actual data; black estimates are from a model where temperature values are randomized within park groups (same park, same home runs, different temperature). (c) Comparison of results from the main data (red) and Statcast data (blue). (top) Main estimate of percent change in home runs per 1°C from all games in parks without domes. (Middle) Percentage-point change in the probability of a batted ball being a home run per 1°C using a linear probability model fitted to the Statcast data (appendix). (bottom) Change in total number of home runs per game per 1°C for both data sources, converted to comparable estimates. In (b) and (c), dots correspond to central estimates while lines span the 95% confidence intervals.

  • Fig. 3.

    Home runs due to anthropogenic climate change. (a) Change in total home runs per season due to anthropogenic climate change over the historical period (black) and across four Shared Socioeconomic Pathway (SSP) emissions scenarios: SSP1–2.6 (teal), SSP2–4.5 (blue), SSP3–7.0 (orange), and SSP5–8.5 (red). Solid lines show the multimodel ensemble mean; shading spans the mean plus or minus the ensemble standard deviation. (b) Change in total home runs per year referenced to global mean surface temperature (GMST) change when the fraction of day and night games is held at historical averages for each park (black; Table ES6) and when all games are night games (blue). Solid line shows the ensemble mean and shading spans the mean plus or minus the ensemble standard deviation for the historical day-game fraction. Dots denote decadal averages of home runs and GMST change, following Pachauri and Meyer (2014). Inset shows cumulative additional home runs at each global warming level in thousands. Bar height shows ensemble mean and error bars span the mean plus or minus the ensemble standard deviation. (c) Additional home runs per year for each individual park at global warming levels of 1.5°, 2°, 3°, and 4°C. Parks are ranked from left to right based on the number of additional home runs at 4°C. Data are pooled across SSP scenarios in (b) and (c).

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 34382 6258 936
PDF Downloads 21397 2570 117