Changes in individual climate variables have been widely documented over the past century. However, assessments that consider changes in the collective interaction amongst multiple climate variables are relevant for understanding climate impacts on ecological and human systems yet are less well documented than univariate changes. We calculate annual multivariate climate departures during 1958–2017 relative to a baseline 1958–1987 period that account for covariance among four variables important to Earth’s biota and associated systems: annual climatic water deficit, annual evapotranspiration, average minimum temperature of the coldest month, and average maximum temperature of the warmest month. Results show positive trends in multivariate climate departures that were nearly three times that of univariate climate departures across global lands. Annual multivariate climate departures exceeded two standard deviations over the past decade for approximately 30% of global lands. Positive trends in climate departures over the last six decades were found to be primarily the result of changes in mean climate conditions consistent with the modeled effects of anthropogenic climate change rather than changes in variability. These results highlight the increasing novelty of annual climatic conditions viewed through a multivariate lens and suggest that changes in multivariate climate departures have generally outpaced univariate departures in recent decades.
Climate change represents a key challenge to Earth’s biota. These changes will require understanding climate change impacts across multiple timescales and variables1. Place-based climate impacts to human-environment systems are often assessed with respect to local climate variability for a variable of interest (e.g., mean annual temperature), as systems are often adapted to such variability2,3,4. Individual climate attributes have already departed from historical ranges of variability in some regions, thereby contributing to climate impacts5,6,7,8,9. However, climate is intrinsically multivariate; many organisms and socio-ecological systems are adapted to and impacted by multiple climate variables and their interactions10,11,12,13,14,15. While climate change is typically portrayed by changes in individual climate variables, often considered independently of one another, changes in multivariate climatic conditions may be more appropriate for anticipating impacts to many systems.
Climate variables exhibit covariance on a variety of timescales through dynamic processes and land-surface interactions16,17. Covariance is particularly acute for water balance variables that are coupled through surface energy-water fluxes10,18 which are widely recognized as important factors controlling the distribution and productivity of agriculture and ecosystems, ecosystem disturbance, and water scarcity19,20,21,22. The covariance structure among climate variables results in combinations (e.g., warm and dry summers) that occur more frequently than if the variables were independent23. In contrast, variable combinations that are orthogonal to the canonical axis of variability (e.g., warm and wet summers) are more unusual in a climatic context and represent larger climate departures (Supplemental Methods). Covariance between variables is apparent on short time scales and contributes to meteorological extremes such as heat waves, and at longer timescales and contributes to climatological extremes such as droughts16. Herein, climate departures refer to absolute departures from a reference climate state in a univariate or multivariate sense. We argue that climate departures may be a valuable metric for tracking potential climate impacts that arise with anomalies of either sign from reference conditions24,25.
Here, we examine trends in multivariate annual climate departures for global land surfaces during 1958–2017 relative to 1958–1987 baseline conditions. We focus on four variables that encapsulate the climatic basis for Earth’s biota and have distinct impacts on human and natural systems: the coldest monthly average minimum temperature (Tn,min), warmest monthly average maximum temperature (Tx,max), annual actual evapotranspiration (AET), and annual climatic water deficit (D). These four variables collectively span thermal and moisture climatic axes and define the climatic niche of many species as well as climate impacts thereof. For instance, Tn,min influences the life history and distributions of species due to the physiological, ecological, and evolutionary impacts26 and is the basis of hardiness zones in the agricultural sector, while Tx,max has been shown to have impacts on electricity demand27 and crop yields28. AET and D are a reduced set of biologically relevant climate variables that represent the supply and unmet atmospheric demand components of the water balance29. AET is a proxy for productivity in natural systems30, while D is widely-used metric for ecosystem disturbance and the terrestrial carbon cycle19,31. In addition, various combinations of these variables are predictor variables for modeling species distributions32,33 and climate impacts34,35.
We use Mahalanobis distance and its standardization based on the Chi distribution3,36 to compute annual multivariate climate departures (σd). The Mahalanobis distance compactly quantifies multiple variables, accounts for their covariance structure, and measures the distance in multivariate space away from a centroid through principal components analysis of standardized anomalies. σd can be interpreted as a multivariate z-score which scales climate departure with respect to local interannual climate variability. We compare σd to its univariate counterpart for individual variables using standardized Euclidean distance (e.g., σTx,max). To quantify the relative contribution of anthropogenic climate forcing to trends in climate departures, we develop a counterfactual simulation that removes the modeled climate change signal from observations. Two sensitivity experiments are used to elucidate the relative contribution of changes in means and changes in intra-annual to inter-annual variability on observed changes in climate departures. Both experiments remove the linear trends for variables during 1958–2017; the first experiment excludes linear trends in annual mean temperature, precipitation, and reference evapotranspiration, while the second only excludes linear trends in annual precipitation. Lastly, we compare observed trends in climate departures with those from internal climate variability using a 500 year control climate simulation.
Positive trends in climate departures are evident in recent decades, particularly for σd (Fig. 1a,b). The four individual calendar years with the greatest global median σd coincided with the four warmest years globally over the period of analysis (2010, 2015–2017). The 60-year median trend in σd across global terrestrial surfaces was +0.78σ. Positive trends in climate departures over the 60-year period were observed for Tx,max, Tn,min, AET and D, but the magnitude of global median trends in σd was approximately three times greater than trends in climate departures for individual variables (Fig. 1e). Similarly, the geographic extent of land with σd exceeding 2 standard deviations increased markedly from a baseline of ~5% of land during the reference period (1958–1987) to ~30% during the most recent decade (2008–2017); increases in the geographic coverage of land >2 standard deviations for σd outpaced increases for individual variables (Fig. 1c,d,f). Similar results were seen using a truncated non-parametric approach for transforming individual variables, although the magnitude of trends was reduced due to the conservative approach (Supplemental Methods). In addition to observed trends, we found a weak correlation (r = 0.21, p = 0.1) between global median σd and the Multivariate ENSO Index37 averaged over Jan-Jun, highlighting the tendency for larger climate departures during El Niño years.
Widespread positive trends in σd were observed over the past six decades across most terrestrial surfaces with statistically significant increases present for 58% of lands and the largest increases in southern Europe, southeast Asia, Africa, and the Amazon (Fig. 2a; Supplementary Fig. 1). Significant positive trends in σTx,max and σTn,min were seen across 16–18% of global lands (Fig. 2b,c). Significant positive trends in σAET were seen for 16% of global lands, primarily across high latitudes of the Northern Hemisphere and equatorial regions. Significant positive trends in σD were found for 11% of global lands. In contrast, significant negative trends in climate departures were confined to small geographic areas of 1–5% of land area depending on variable.
Geographic hotspots of large positive trends often occurred in areas with low interannual variability. This is further demonstrated by examining trends in climate departures and the magnitude of interannual variability by latitude and biome. Latitudinal patterns show large positive trends in σd and σTn,min near the equator co-located with regions with low variability during the reference period (Fig. 3a). Additionally, larger positive σd trends were seen in the Northern Hemisphere than the Southern Hemisphere (Fig. 3a). Trends for temperature-based departures generally exceeded moisture-based departures from 40°S to 50°N. At higher-latitudes, trends in σAET outpaced trends in temperature departures consistent with higher temperature variability at higher latitudes (Fig. 3b). Trends in climate departures by biome were generally largest where interannual variability among the four climate variables was low (Fig. 3c,d), with the largest positive σd trend in the Tropical Forest biome and the smallest positive trend in the Temperate Grassland biome.
Larger positive σd trends compared to univariate climate departure trends can be understood by considering the covariance structure and trends in the underlying climate variables. AET and D exhibited strong negative correlations on interannual timescales across most terrestrial surfaces, with Tx,max exhibiting positive correlation to D (Supplementary Fig. 3). By contrast, Tn,min was weakly correlated to the other climate variables, consistent with Tn,min being seasonally out of phase with the primary climatic factors that influence surface water availability. The principal components (PC) used in the Mahalanobis distance calculation largely reflect these relationships with the first loading projecting onto patterns of water limitations, loadings 2 and 3 projecting onto orthogonal dimensions related to Tn,min and Tx,max, respectively, and loading 4 projecting to same signed anomalies of AET and D (Supplementary Fig. 4). Positive trends in Tx,max and Tn,min were seen across most terrestrial surfaces with AET and D trends showing more spatial variability but generally positive trends during 1958–2017 (Supplementary Fig. 5). Trends in σd, by definition, arise from trends in Euclidean distances of PC scores. While no single loading predominantly accounted for σd trends, trends in PC4 and PC3 accounted for a majority of positive σd trends for 22% and 16% of the globe, respectively (Supplementary Table 1). Positive trends in AET and D were jointly observed for 46% of terrestrial surfaces that map onto the PC4 loading. Likewise, joint significant increases in Tx,max and Tn,min were observed for 24% of the globe and contributed to the positive trend in σd as Tx,max and Tn,min were poorly correlated on interannual timescales during the reference period outside of the tropics.
Approximately half of the observed increase in σd can be accounted for by anthropogenic climate change (Fig. 1e, Supplementary Fig. 2a). The counterfactual simulation produced: 1) σd trends with a global median of +0.36σ and were spatially similar to that observed, 2) mostly non-significant trends in temperature based departures, and 3) trends in moisture based departures that were approximately a third the magnitude of observed trends. Similarly, excluding the effect of anthropogenic climate change led to a 45% reduction in the extent of global land area with σd > 2σ (Fig. 1f).
Several factors account for the σd increase after excluding the modeled influence of anthropogenic climate change, including changes in climate variability and divergence between observed and modeled changes in climate. The sensitivity experiment that excluded trends in annual temperature, precipitation, and reference evapotranspiration had σd trends that were approximately a third of those observed and trends in climate departures for individual variables that were less than 30% of those observed (Fig. 1e, Supplementary Fig. 1c). By contrast, the sensitivity experiment that excluded trends in annual precipitation had nominal influence in attenuating observed trends in climate departures thus highlighting the importance of changes in temperature and evaporative demand (Supplementary Fig. 1b). Our results further suggest that most of the observed increase in climate departures have materialized through changes in mean climate conditions rather than changes in climate variability. Changes in the interannual variability of individual variables showed confounding signals spatially and across variables (Supplemental Methods; Supplementary Table 2). The most notable and widespread change was a global median 15% reduction in the standard deviation of Tn,min during 1958–2017. The larger fraction of trends accounted for by explicitly detrending observations versus those using the counterfactual model is expected given that observed σd trends include changes arising from internal variability.
Lastly, the magnitude of observed changes in climate departures exceeded those estimated from natural variability in a stationary climate. A null model based on 500-years of climate model output run under pre-industrial control simulations38 yielded a global median σd trend comparable to that of the detrended sensitivity experiment (average trend +0.24σ over 60-years, maximum of +0.29σ; Supplementary Fig. 6). Internal decadal variability along with the use of a reference period from which climate departures are calculated and associated sampling biases generally results in positive σd trends even in a stationary climate. However, the substantially larger magnitude of observed changes versus those obtained using the null model suggest that the observed increases in climate departures over our study period are unlikely due to internal climate variability.
We demonstrate that annual climate conditions over the past three decades have departed substantially from previous decades, that these climate departures exceed what would be expected from natural climate variability, and that approximately half of the magnitude of these changes can be attributed to anthropogenic forcings. While previous studies have highlighted the importance of multivariate extremes and changes thereof11,13, we find that multivariate climate departures assessed at annual time scales have rapidly outpaced univariate departures, highlighting the potential latent impacts imposed by observed change to systems adapted to multiple interacting climate variables. The largest positive σd trends occurred in regions of historically low variance (e.g., equatorial regions)8,39,40, in regions such as southern Europe that have seen large changes in climate trends in the variables considered41, and in boreal regions that saw joint increases in AET and D that are orthogonal to the historical covariance of these variables. The growing extent of lands with annual multivariate climate departures>2σ in recent decades complements other studies that have considered univariate departures42. While our results were specific to the four variables than span moisture and thermal constraints in a given year, trends in multivariate climate departures may differ as a function of the variables chosen.
Trends in multivariate climate departures were largely accounted for by changes in climate means rather than changes in climate variability. While there is a perception that observed climate change has resulted in heightened variability, observational and modeling studies generally do not support widespread detectable changes to date43,44. Nonetheless, while observed changes in variability have generally been small, some studies suggest increased hydroclimatic variability or intensification over the 21st century that would further exacerbate increased climate departures45,46,47.
Our framing of univariate and multivariate departures builds on previous studies3,39,48 while emphasizing more nuanced ways in which climate has differed from reference conditions on annual timescales. Only considering changes in mean conditions or changes for multi-decadal timescales can obscure departures that occur annually. This type of interannual variability is critical for assessing climate change impacts. For example, years with concurrent high D and Tx,max enable fire activity in flammability-limited forests19, while years with high moisture availability (high AET, low D) and low Tx,max promote recruitment pulses in water-limited forests49. Warm summers accompanied by surface water availability (high Tx,max, high AET) may impact human health directly through heat stress and indirectly through promoting environmental conditions that favor certain vector-borne diseases50,51. Annual climate departures are also relevant to impacts on agricultural systems, land-use planning, and built infrastructure52,53.
As the climate becomes increasingly unfamiliar to the flora and fauna at a specific location, organisms, including humans, must adapt to changing local conditions or move to maintain suitable climates54. Indeed, Earth’s biota is already adapting to climate change. For example, migration and thermophilization of biotic communities has been observed55,56. Additionally, humans are increasingly adopting crop varieties, seed provenances, and horticultural practices from other locations with climates that are anticipated to resemble future conditions3,57. Unfortunately, the adaptive capacity of systems is not necessarily commensurate with their climate change exposure. For example, thermal specialists in the tropics may be uniquely vulnerable to the large climate departures seen in these regions58,59. Collectively, our results suggest that annual multivariate climate departures have changed more dramatically in the last 60 years than univariate estimates of departures for individual climate variables.
Monthly data from TerraClimate during 1958–2017 were used for our primary analysis60. TerraClimate combines high-resolution multi-decadal climatologies from WorldClim61 with coarser-scale temporally varying data from reanalyses and CRU Ts4.062 to create 4-km (1/24°) spatial resolution surfaces of monthly first-order climate variables (e.g., temperature, precipitation), as well as surface water balance products such as AET and D. The water balance model is an accounting-based approach that simulates moisture fluxes using a simple snow model and soil water balance approach for a reference vegetation type. This modeling framework has been widely used in hydrological and ecological studies10,63. Other gridded datasets provide similar outputs, including those using more sophisticated surface hydrologic models. However, they do not cover the 60-year period of observational record or are of much coarser spatial resolution. Nonetheless, to interrogate the sensitivity of our results with respect to structural uncertainty of climate data sources, we replicate all analyses for the period 1958–2016 using gridded data from the Princeton Global Meteorological Forcing dataset (version 3) at a 0.25° spatial resolution64. Results from these data are presented in Supplemental Methods.
Four climate metrics were chosen given their established links to factors that influence Earth’s biota including species occurrence and impacts, and their frequent usage in species distribution models10,29,32,65. These metrics included: (1) average maximum temperature of the warmest month for the calendar year (Tx,max), (2) average minimum temperature of the coldest month for the calendar year (Tn,min), (3) calendar year cumulative AET, and (4) calendar year cumulative D. These four variables provide complementary information pertinent to ecological and agricultural systems. For example, Tn,min can be a limiting factor for some species and agriculture and modify overwinter morality rates of some organisms66,67, while D and AET provide a reduced set of biologically relevant and physically based variables that account for the concurrent availability of both water and energy important for both ecological and agricultural impacts29,68.
We calculated σd and its univariate equivalent standardized Euclidean distances for each variable during 1958–2017. The first 30-years (1958–1987) were used to define a baseline period that determine both the centroid (i.e. the mean) and covariance structure using principal components analysis to calculate σd (Supplemental Methods). A given 30-year period comprises a sample of a population of climate statistics that may differ from other baseline periods in some regions due to decadal variability69. Likewise, anthropogenic climate forcing is evident throughout the observational record limiting the ability to procure a true ‘natural’ reference state. However, for the basis of this study, the earliest 30-years covers a time period when anthropogenic forcing was substantially less than in recent decades.
We used Mahalanobis distances and its standardization based on the Chi distribution to calculate σd. This approach accounts for covariance among variables, the dimensionality of the data (number of variables), and local inter-annual climate variability thus facilitating comparisons across space and variable combinations. Mahalanobis distances were calculated on standardized data (i.e., normal distributions based on means and standard deviation calculated during 1958–1987). This approach assumes variables adhere to a multivariate Gaussian distribution, which may be reasonable for some variables, but could be problematic for zero-bounded data (e.g., D) and for temperature extremes in some portions of the globe. We conducted a sensitivity analysis that used clamped non-parametric distributions to assess whether results were substantially altered by these assumptions (Supplemental Methods). The choice of data transformations did not substantially alter the general results of the study, although the truncated nature of the non-parametric distributions reduced overall magnitudes of departures. More advanced multivariate approaches such as copulas may provide additional nuance beyond the linear approaches used herein70. Trends for σd and climate departure for individual variables during 1958–2017 were calculated using Sen-Theil slope estimator and were considered significant using the Mann-Kendall trend test at p < 0.05. Similarly, standard trends using the same approach were calculated on the raw climate variables.
A counterfactual simulation of terrestrial climate was developed that excludes the modeled influence of anthropogenic climate forcing during 1958–2017. The counterfactual simulation removed first-order modeled trends in monthly climate from observed data similar to that of previous studies71,72,73. We used a pattern scaling approach that applies the median response among 23 different climate models to proximate modeled anthropogenic changes (Supplemental Methods)73. We also considered two sensitivity tests that use detrended observations to decompose changes in climate departures into those associated with trends versus those associated with changes in intraannual-to-interannual variability. One approach removed only annual precipitation trends, while the other approach removed annual trends in temperature, precipitation, and reference evapotranspiration. Detrending was facilitated by calculating a linear least square fit on annual data. This linear fit was subtracted from the observed time series leaving data for 1958 unaltered. By only removing annual trends, we allow monthly trends to differ in a relative sense from annual counterparts. Data for the counterfactual simulation and sensitivity experiments were run through the water balance model to calculate AET and D.
We consider two additional filters in an effort to avoid misinterpreting results in locations of reduced data quality or where the baseline observations were heavily positively-skewed. First, we omit calculations in locations where annual D or AET was 0 for a majority of the years in the baseline period. Second, we used the data quality flags in TerraClimate inherited from CRU TS4.0 to select pixels where at least four stations contributed to monthly precipitation, temperature, and vapor pressure fields for at least 75% of the period of record. Fig. S4 shows the data from Fig. 2 with data poor locations masked out. Summarized timeseries and reported trends for global land surfaces exclude pixels not meeting either of these criteria.
Lastly, we use 500-years (model years 400–899) of pre-industrial control climate simulated from the LENS experiment which uses NCAR’s Community Earth System Model version 1 (CESM1) with CAM5.2 as its atmospheric model38 to develop a null model for changes in σd that result purely from internal climate variability. Monthly output from LENS was subsequently used in the water balance model, albeit at the native spatial resolution of LENS output. We subsequently calculated σd using non-overlapping moving 30-year blocks over the simulation (e.g., 430–459) and calculated both forward (e.g., 430–489) and backward (e.g., 400–459) linear trends in global terrestrial median σd and the fraction of land surfaces with σd > 2σ. The slope of trends calculated from backward samples was inverted. As with observations, trends were calculated using Sen-Theil slope estimator.
The primary datasets analyzed in the current study are available through the Northwest Knowledge Network data repository at https://climate.northwestknowledge.net/TERRACLIMATE-DATA/.
Harris, R. M. B. et al. Biological responses to the press and pulse of climate trends and extreme events. Nat. Clim. Chang. 8, 579 (2018).
Kawecki, T. J. & Ebert, D. Conceptual issues in local adaptation. Ecol. Lett. 7, 1225–1241 (2004).
Mahony, C. R., Cannon, A. J., Wang, T. & Aitken, S. N. A closer look at novel climates: new methods and insights at continental to landscape scales. Glob. Chang. Biol. (2017).
Flannigan, M. D. & Harrington, J. B. A study of the relation of meteorological variables to monthly provincial area burned by wildfire in Canada (1953-80). J. Appl. Meteorol. 27, 441–452 (1988).
Piao, S. et al. The impacts of climate change on water resources and agriculture in China. Nature 467, 43 (2010).
Duffy, P. B. et al. Strengthened scientific support for the Endangerment Finding for atmospheric greenhouse gases. Science 363, eaat5982 (2019).
Abatzoglou, J. T., Williams, A. P. & Barbero, R. Global emergence of anthropogenic climate change in fire weather indices. Geophys. Res. Lett. 46, 326–336 (2019).
Mahlstein, I., Knutti, R., Solomon, S. & Portmann, R. W. Early onset of significant local warming in low latitude countries. Environ. Res. Lett. 6, 34009 (2011).
Frame, D., Joshi, M., Hawkins, E., Harrington, L. J. & de Roiste, M. Population-based emergence of unfamiliar climates. Nat. Clim. Chang. 7, 407 (2017).
Dobrowski, S. Z. et al. The climate velocity of the contiguous United States during the 20th century. Glob. Chang. Biol. 19, 241–251 (2013).
Hao, Z., AghaKouchak, A. & Phillips, T. J. Changes in concurrent monthly precipitation and temperature extremes. Environ. Res. Lett. 8, 34014 (2013).
Mazdiyasni, O. & AghaKouchak, A. Substantial increase in concurrent droughts and heatwaves in the United States. Proc. Natl. Acad. Sci. 112, 11484–11489 (2015).
Sarhadi, A., Ausín, M. C., Wiper, M. P., Touma, D. & Diffenbaugh, N. S. Multidimensional risk in a nonstationary climate: Joint probability of increasingly severe warm and dry conditions. Sci. Adv. 4, eaau3487 (2018).
Zhou, S., Zhang, Y., Williams, A. P. & Gentine, P. Projected increases in intensity, frequency, and terrestrial carbon costs of compound drought and aridity events. Sci. Adv. 5, eaau5740 (2019).
Mora, C. et al. Broad threat to humanity from cumulative climate hazards intensified by greenhouse gas emissions. Nat. Clim. Chang. 8, 1062–1071 (2018).
Trenberth, K. E. & Shea, D. J. Relationships between precipitation and surface temperature. Geophys. Res. Lett. 32 (2005).
Seneviratne, S. I. et al. Investigating soil moisture–climate interactions in a changing climate: A review. Earth-Science Rev. 99, 125–161 (2010).
Willmott, C. J., Rowe, C. M. & Mintz, Y. Climatology of the terrestrial seasonal water cycle. J. Climatol. 5, 589–606 (1985).
Abatzoglou, J. T., Williams, A. P., Boschetti, L., Zubkova, M. & Kolden, C. A. Global patterns of interannual climate-fire relationships. Glob. Chang. Biol. 24, 5164–5175 (2018).
Allen, C. D. et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manage. 259, 660–684 (2010).
van Mantgem, P. J. et al. Climatic stress increases forest fire severity across the western U nited S tates. Ecol. Lett. 16, 1151–1156 (2013).
Wada, Y. et al. Multimodel projections and uncertainties of irrigation water demand under climate change. Geophys. Res. Lett. 40, 4626–4632 (2013).
Zscheischler, J. & Seneviratne, S. I. Dependence of drivers affects risks associated with compound events. Sci. Adv. 3, e1700263 (2017).
Burke, M., Hsiang, S. M. & Miguel, E. Global non-linear effect of temperature on economic production. Nature 527, 235 (2015).
Schlenker, W. & Roberts, M. J. Nonlinear temperature effects indicate severe damages to US crop yields under climate change. Proc. Natl. Acad. Sci. 106, 15594–15598 (2009).
Inouye, D. W. The ecological and evolutionary significance of frost in the context of climate change. Ecol. Lett. 3, 457–463 (2000).
Apadula, F., Bassini, A., Elli, A. & Scapin, S. Relationships between meteorological variables and monthly electricity demand. Appl. Energy 98, 346–356 (2012).
Deryng, D., Conway, D., Ramankutty, N., Price, J. & Warren, R. Global crop yield response to extreme heat stress under multiple climate change futures. Environ. Res. Lett. 9, 34011 (2014).
Stephenson, N. L. Climatic control of vegetation distribution: the role of the water balance. Am. Nat. 649–670 (1990).
Taylor, P. G. et al. Temperature and rainfall interact to control carbon cycling in tropical forests. Ecol. Lett. 20, 779–788 (2017).
Anderegg, W. R. L. et al. Pervasive drought legacies in forest ecosystems and their implications for carbon cycle models. Science 349, 528–532 (2015).
Lutz, J. A., van Wagtendonk, J. W. & Franklin, J. F. Climatic water deficit, tree species ranges, and climate change in Yosemite National Park. J. Biogeogr. 37, 936–950 (2010).
Parks, S. A., Parisien, M.-A., Miller, C. & Dobrowski, S. Z. Fire activity and severity in the western US vary along proxy gradients representing fuel amount and fuel moisture. Plos One 9, e99699 (2014).
Westerling, A. L., Turner, M. G., Smithwick, E. A. H., Romme, W. H. & Ryan, M. G. Continued warming could transform Greater Yellowstone fire regimes by mid-21st century. Proc. Natl. Acad. Sci. USA 108, 13165–13170 (2011).
Jackson, S. T., Betancourt, J. L., Booth, R. K. & Gray, S. T. Ecology and the ratchet of events: Climate variability, niche dimensions, and species distributions. Proc. Natl. Acad. Sci. 106, 19685 LP–19692 (2009).
Fitzpatrick, M. C. & Dunn, R. R. Contemporary climatic analogs for 540 North American urban areas in the late 21st century. Nat. Commun. 10, 614 (2019).
Wolter, K. & Timlin, M. S. El Niño/Southern Oscillation behaviour since 1871 as diagnosed in an extended multivariate ENSO index (MEI. ext). Int. J. Climatol. 31, 1074–1087 (2011).
Kay, J. E. et al. The Community Earth System Model (CESM) large ensemble project: A community resource for studying climate change in the presence of internal climate variability. Bull. Am. Meteorol. Soc. 96, 1333–1349 (2015).
Mora, C. et al. The projected timing of climate departure from recent variability. Nature 502, 183 (2013).
Hawkins, E. & Sutton, R. Time of emergence of climate signals. Geophys. Res. Lett. 39, (2012).
Cook, B. I., Anchukaitis, K. J., Touchan, R., Meko, D. M. & Cook, E. R. Spatiotemporal drought variability in the Mediterranean over the last 900 years. J. Geophys. Res. Atmos. 121, 2060–2074 (2016).
Coumou, D. & Robinson, A. Historic and future increase in the global land area affected by monthly heat extremes. Environ. Res. Lett. 8, 34018 (2013).
Huntingford, C., Jones, P. D., Livina, V. N., Lenton, T. M. & Cox, P. M. No increase in global temperature variability despite changing regional patterns. Nature 500, 327 (2013).
Screen, J. A. Arctic amplification decreases temperature variance in northern mid-to high-latitudes. Nat. Clim. Chang. 4, 577–582 (2014).
Swain, D. L., Langenbrunner, B., Neelin, J. D. & Hall, A. Increasing precipitation volatility in twenty-first-century California. Nat. Clim. Chang. 8, 427 (2018).
Diffenbaugh, N. S. & Ashfaq, M. Intensification of hot extremes in the United States. Geophys. Res. Lett. 37, (2010).
Seager, R., Naik, N. & Vogel, L. Does global warming cause intensified interannual hydroclimate variability? J. Clim. 25, 3355–3372 (2012).
Mahony, C. R. & Cannon, A. J. Wetter summers can intensify departures from natural variability in a warming climate. Nat. Commun. 9, 783 (2018).
Davis, K. T. et al. Wildfires and climate change push low-elevation forests across a critical climate threshold for tree regeneration. Proc. Natl. Acad. Sci. 116, 6193–6198 (2019).
Sherwood, S. C. & Huber, M. An adaptability limit to climate change due to heat stress. Proc. Natl. Acad. Sci. 107, 9552–9555 (2010).
Ebi, K. L. & Nealon, J. Dengue in a changing climate. Environ. Res. 151, 115–123 (2016).
Hatfield, J. et al. Climate change impacts in the United States: The third national climate assessment. Washington, DC 150–174 (2014).
Sivakumar, M. V. K., Das, H. P. & Brunini, O. Impacts of present and future climate variability and change on agriculture and forestry in the arid and semi-arid tropics. in increasing climate variability and change 31–72 (Springer, 2005).
Parmesan, C. Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. Syst. 37, 637–669 (2006).
Chen, I.-C., Hill, J. K., Ohlemüller, R., Roy, D. B. & Thomas, C. D. Rapid range shifts of species associated with high levels of climate warming. Science 333, 1024–1026 (2011).
Bertrand, R. et al. Changes in plant community composition lag behind climate warming in lowland forests. Nature 479, 517 (2011).
Thomas, C. D. Translocation of species, climate change, and the end of trying to recreate past ecological communities. Trends Ecol. Evol. 26, 216–221 (2011).
Tewksbury, J. J., Huey, R. B. & Deutsch, C. A. Putting the heat on tropical animals. Science 320, 1296–1297 (2008).
Huey, R. B. et al. Predicting organismal vulnerability to climate warming: roles of behaviour, physiology and adaptation. Philos. Trans. R. Soc. B Biol. Sci. 367, 1665–1679 (2012).
Abatzoglou, J. T., Dobrowski, S. Z., Parks, S. A. & Hegewisch, K. C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 5, 170191 (2018).
Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol., https://doi.org/10.1002/joc.5086 (2017).
Harris, I., Jones, P. D., Osborn, T. J. & Lister, D. H. Updated high-resolution grids of monthly climatic observations – the CRU TS3.10 Dataset. Int. J. Climatol. 34, 623–642 (2014).
Gleick, P. H. The development and testing of a water balance model for climate impact assessment: modeling the Sacramento basin. Water Resour. Res. 23, 1049–1061 (1987).
Sheffield, J., Goteti, G. & Wood, E. F. Development of a 50-year high-resolution global dataset of meteorological forcings for land surface modeling. J. Clim. 19, 3088–3111 (2006).
Parks, S. A. et al. Wildland fire deficit and surplus in the western United States, 1984–2012. Ecosphere 6, 1–13 (2015).
Parker, L. E. & Abatzoglou, J. T. Projected changes in cold hardiness zones and suitable overwinter ranges of perennial crops over the United States. Environ. Res. Lett. 11, 34001 (2016).
Williams, C. M., Henry, H. A. L. & Sinclair, B. J. Cold truths: how winter drives responses of terrestrial organisms to climate change. Biol. Rev. 90, 214–235 (2015).
Wriedt, G., V der Velde, M., Aloe, A. & Bouraoui, F. Estimating irrigation water requirements in Europe. J. Hydrol. 373, 527–544 (2009).
Hulme, M. & New, M. Dependence of large-scale precipitation climatologies on temporal and spatial sampling. J. Clim. 10, 1099–1113 (1997).
AghaKouchak, A., Cheng, L., Mazdiyasni, O. & Farahmand, A. Global warming and changes in risk of concurrent climate extremes: Insights from the 2014 California drought. Geophys. Res. Lett. 41, 8847–8852 (2014).
Abatzoglou, J. T. & Williams, A. P. Impact of anthropogenic climate change on wildfire across western US forests. Proc. Natl. Acad. Sci. 113, 11770–11775 (2016).
Williams, A. P. et al. Contribution of anthropogenic warming to California drought during 2012–2014. Geophys. Res. Lett. 42, 6819–6828 (2015).
Mitchell, T. D. Pattern scaling: an examination of the accuracy of the technique for describing future climates. Clim. Change 60, 217–242 (2003).
J.T.A. was partially supported by the National Science Foundation under award DMS-1520873. S.Z.D. was partly funded by National Science Foundation Grant BCS-1461576 and the USDA National Institute of Food and Agriculture, McIntire Stennis program, project 1012438.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Abatzoglou, J.T., Dobrowski, S.Z. & Parks, S.A. Multivariate climate departures have outpaced univariate changes across global lands. Sci Rep 10, 3891 (2020). https://doi.org/10.1038/s41598-020-60270-5
Environmental Research Letters (2021)
Climate and biotic evolution during the Permian-Triassic transition in the temperate Northern Hemisphere, Kuznetsk Basin, Siberia, Russia
Palaeogeography, Palaeoclimatology, Palaeoecology (2021)
Nature Ecology & Evolution (2021)
Scientific Reports (2020)