Paleoclimate Model-Data Comparisons of Hydroclimate over North America with a Focus on Megadroughts

For the first time in the history of the Intergovernmental Panel on Climate Change, the Paleoclimate and Coupled Model Intercomparison Projects Phases 3 and 5 produced an ensemble of forced transient simulations of the last millennium. This wealth of model output, when combined with a growing collection of high spatial and temporal resolution pa- leoclimate estimates of past climate variability, represents an important and unprecedented source of information on climate variability over decades to centuries. This dissertation thus combines paleoclimate evidence with climate modeling to define a physical and statistical paradigm through which to analyze these combined sources of information and subsequently to characterize the features of the North American climate system that cannot be sufficiently understood using instrumental data alone. This includes features that have long timescales of variability or that are rare, and by consequence have few degrees of freedom over the short instrumental interval (1850 C.E. to Present), as well as interannual dynamical relationships that, while potentially well characterized by observations, are non-stationary. An integrative approach to analyzing these features or relationships serves two fundamental purposes: 1) It provides a more comprehensive characterization of past climate variability, albeit with the caveat of model bias, to clarify understanding of the dynamics that produce these features or relationships in the real world; and 2) it assesses whether coupled general circulation models (CGCMs) are able to simulate these features or relationships, which is necessary to determine that state-of-the-art CGCMs can accurately constrain the risk of future climate change. The focus herein will be on multidecadal hydroclimate change, or megadroughts, in the paleoclimate record of the American Southwest to better inform our understanding of the risk of future hydroclimate change over the region. Two fundamental understandings are derived from this work. Firstly, CGCMs are successful at simulating discrete periods of mul- tidecadal hydroclimate change that are characteristic in length, magnitude, and frequency of occurrence of megadroughts in the paleoclimate record. The simulated megadroughts are not tied in any coherent way to exogenous forcing, however, suggesting that CGCMs simulate large-magnitude internal variability on multidecadal timescales. Secondly, the dynamical characteristics of CGCMs are important in determining the atmosphere-ocean variability that drives multidecadal hydroclimate change. The dynamical characteristics of relevance include teleconnection realism and stationarity, the magnitude of ocean variability, and the relative magnitudes of different modes of atmosphere-ocean and purely atmospheric vari- ability. Additionally, a new understanding of real-world megadrought dynamics is derived herein, with the characteristics of some CGGMs providing a better representation of these dynamics

Model-dependent spatial skill in pseudoproxy experiments testing climate field reconstruction methods for the Common Era

The spatial skill of four climate field reconstruction (CFR) methods is investigated using pseudoproxy experiments (PPEs) based on five last millennium and historical simulations from the Coupled and Paleo Model Intercomparison Projects Phases 5 and 3 (CMIP5/PMIP3) data archives. These simulations are used for the first time in a PPE context, the frameworks of which are constructed to test a recently assembled multiproxy network and multiple CFR techniques. The experiments confirm earlier findings demonstrating consistent methodological performance across the employed methods and spatially dependent reconstruction errors in all of the derived CFRs. Spectral biases in the reconstructed fields demonstrate that CFR methods can alone alter the ratio of spectral power at all locations in the field, independent of whether there are any spectral biases inherent in the underlying pseudoproxy series. The patterns of spectral biases are model dependent and indicate the potential for regions in the derived CFRs to be biased by changes in either low or high-frequency spectral power. CFR methods are also shown to alter the pattern of mean differences in the tropical Pacific during the Medieval Climate Anomaly and the Little Ice Age, with some model experiments indicating that CFR methodologies enhance the statistical likelihood of achieving larger mean differences between independent 300-year periods in the region. All of the characteristics of CFR performance are model dependent, indicating that CFR methods must be evaluated across multiple models and that conclusions from PPEs should be carefully connected to the spatial statistics of real-world climatic fields

North American Pancontinental Droughts in Model Simulations of the Last Millennium

Pancontinental droughts in North America, or droughts that simultaneously affect a large percentage of the geographically and climatically distinct regions of the continent, present significant on-the-ground management challenges and, as such, are an important target for scientific research. The methodology of paleoclimate-model data comparisons is used herein to provide a more comprehensive understanding of pancontinental drought dynamics. Models are found to simulate pancontinental drought with the frequency and spatial patterns exhibited by the paleoclimate record. They do not, however, agree on the modes of atmosphere–ocean variability that produce pancontinental droughts because simulated El Niño–Southern Oscillation (ENSO), Pacific decadal oscillation (PDO), and Atlantic multidecadal oscillation (AMO) dynamics, and their teleconnections to North America, are different between models and observations. Despite these dynamical differences, models are able to reproduce large-magnitude centennial-scale variability in the frequency of pancontinental drought occurrence—an important feature of the paleoclimate record. These changes do not appear to be tied to exogenous forcing, suggesting that simulated internal hydroclimate variability on these time scales is large in magnitude. Results clarify our understanding of the dynamics that produce real-world pancontinental droughts while assessing the ability of models to accurately characterize future drought risks

Are Simulated Megadroughts in the North American Southwest Forced?

Multidecadal drought periods in the North American Southwest (25°–42.5°N, 125°–105°W), so-called megadroughts, are a prominent feature of the paleoclimate record over the last millennium (LM). Six forced transient simulations of the LM along with corresponding historical (1850–2005) and 500-yr preindustrial control runs from phase 5 of the Coupled Model Intercomparison Project (CMIP5) are analyzed to determine if atmosphere–ocean general circulation models (AOGCMs) are able to simulate droughts that are similar in persistence and severity to the megadroughts in the proxy-derived North American Drought Atlas. Megadroughts are found in each of the AOGCM simulations of the LM, although there are intermodel differences in the number, persistence, and severity of these features. Despite these differences, a common feature of the simulated megadroughts is that they are not forced by changes in the exogenous forcing conditions. Furthermore, only the Community Climate System Model (CCSM), version 4, simulation contains megadroughts that are consistently forced by cooler conditions in the tropical Pacific Ocean. These La Niña–like mean states are not accompanied by changes to the interannual variability of the El Niño–Southern Oscillation system and result from internal multidecadal variability of the tropical Pacific mean state, of which the CCSM has the largest magnitude of the analyzed simulations. Critically, the CCSM is also found to have a realistic teleconnection between the tropical Pacific and North America that is stationary on multidecadal time scales. Generally, models with some combination of a realistic and stationary teleconnection and large multidecadal variability in the tropical Pacific are found to have the highest incidence of megadroughts driven by the tropical Pacific boundary conditions

The improbable but unexceptional occurrence of megadrought clustering in the American West during the Medieval Climate Anomaly

The five most severe and persistent droughts in the American West (AW) during the Common Era occurred during a 450 year period known as the Medieval Climate Anomaly (MCA—850–1299 C.E.). Herein we use timeseries modeling to estimate the probability of such a period of hydroclimate change occurring. Clustering of severe and persistent drought during an MCA-length period occurs in approximately 10% of surrogate timeseries that were constructed to have the same characteristics as a tree-ring derived estimate of AW hydroclimate variability between 850 and 2005 C.E. Periods of hydroclimate change like the MCA are thus expected to occur in the AW, although not frequently, with a recurrence interval of approximately 11 000 years. Importantly, a shift in mean hydroclimate conditions during the MCA is found to be necessary for drought to reach the severity and persistence of the actual MCA megadroughts. This result has consequences for our understanding of the atmosphere-ocean dynamics underlying the MCA and a persistently warm Atlantic Multidecadal Oscillation is suggested to have played an important role in causing megadrought clustering during this period

Paleoclimate constraints on the spatiotemporal character of past and future droughts

Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 33(22), (2020): 9883-9903, https://doi.org/10.1175/JCLI-D-20-0004.1.Machine-learning-based methods that identify drought in three-dimensional space–time are applied to climate model simulations and tree-ring-based reconstructions of hydroclimate over the Northern Hemisphere extratropics for the past 1000 years, as well as twenty-first-century projections. Analyzing reconstructed and simulated drought in this context provides a paleoclimate constraint on the spatiotemporal characteristics of simulated droughts. Climate models project that there will be large increases in the persistence and severity of droughts over the coming century, but with little change in their spatial extent. Nevertheless, climate models exhibit biases in the spatiotemporal characteristics of persistent and severe droughts over parts of the Northern Hemisphere. We use the paleoclimate record and results from a linear inverse modeling-based framework to conclude that climate models underestimate the range of potential future hydroclimate states. Complicating this picture, however, are divergent changes in the characteristics of persistent and severe droughts when quantified using different hydroclimate metrics. Collectively our results imply that these divergent responses and the aforementioned biases must be better understood if we are to increase confidence in future hydroclimate projections. Importantly, the novel framework presented herein can be applied to other climate features to robustly describe their spatiotemporal characteristics and provide constraints on future changes to those characteristics.This material is based upon work supported by the National Center for Atmospheric Research, which is a major facility sponsored by the National Science Foundation under Cooperative Agreement 1852977. JAF was also supported by the Regional and Global Model Analysis (RGMA) component of the Earth and Environmental System Modeling Program of the U.S. Department of Energy's Office of Biological & Environmental Research (BER) via National Science Foundation IA 1844590. JS was supported in part by the U.S. National Science Foundation through Grants AGS-1602920 and AGS-1805490, and by the National Oceanic and Atmospheric Administration by Grant NA20OAR4310425. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1) for producing and making available their model output. For CMIP, the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portal. We thank the editor and two reviewers for comments that greatly improved the quality of this manuscript. This is SOEST Publication No. 11116 and LDEO Publication No. 8450.2021-04-1

Internal ocean-atmosphere variability drives megadroughts in Western North America

Multidecadal droughts that occurred during the Medieval Climate Anomaly represent an important target for validating the ability of climate models to adequately characterize drought risk over the near-term future. A prominent hypothesis is that these megadroughts were driven by a centuries-long radiatively forced shift in the mean state of the tropical Pacific Ocean. Here we use a novel combination of spatiotemporal tree ring reconstructions of Northern Hemisphere hydroclimate to infer the atmosphere-ocean dynamics that coincide with megadroughts over the American West and find that these features are consistently associated with 10–30 year periods of frequent cold El Niño–Southern Oscillation conditions and not a centuries-long shift in the mean of the tropical Pacific Ocean. These results suggest an important role for internal variability in driving past megadroughts. State-of-the-art climate models from the Coupled Model Intercomparison Project Phase 5, however, do not simulate a consistent association between megadroughts and internal variability of the tropical Pacific Ocean, with implications for our confidence in megadrought risk projections

Projected drought risk in 1.5°C and 2°C warmer climates

The large socioeconomic costs of droughts make them a crucial target for impact assessments of climate change scenarios. Using multiple drought metrics and a set of simulations with the Community Earth System Model targeting 1.5°C and 2°C above preindustrial global mean temperatures, we investigate changes in aridity and the risk of consecutive drought years. If warming is limited to 2°C, these simulations suggest little change in drought risk for the U.S. Southwest and Central Plains compared to present day. In the Mediterranean and central Europe, however, drought risk increases significantly for both 1.5°C and 2°C warming targets, and the additional 0.5°C of the 2°C climate leads to significantly higher drought risk. Our study suggests that limiting anthropogenic warming to 1.5°C rather than 2°C, as aspired to by the Paris Climate Agreement, may have benefits for future drought risk but that such benefits may be regional and in some cases highly uncertain

Megadroughts in Southwestern North America in ECHO-G Millennial Simulations and Their Comparison to Proxy Drought Reconstructions

Simulated hydroclimate variability in millennium-length forced transient and control simulations from the ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G) coupled atmosphere–ocean general circulation model (AOGCM) is analyzed and compared to 1000 years of reconstructed Palmer drought severity index (PDSI) variability from the North American Drought Atlas (NADA). The ability of the model to simulate megadroughts in the North American southwest is evaluated. (NASW: 25°–42.5°N, 125°–105°W). Megadroughts in the ECHO-G AOGCM are found to be similar in duration and magnitude to those estimated from the NADA. The droughts in the forced simulation are not, however, temporally synchronous with those in the paleoclimate record, nor are there significant differences between the drought features simulated in the forced and control runs. These results indicate that model-simulated megadroughts can result from internal variability of the modeled climate system rather than as a response to changes in exogenous forcings. Although the ECHO-G AOGCM is capable of simulating megadroughts through persistent La Niña–like conditions in the tropical Pacific, other mechanisms can produce similarly extreme NASW moisture anomalies in the model. In particular, the lack of low-frequency coherence between NASW soil moisture and simulated modes of climate variability like the El Niño–Southern Oscillation, Pacific decadal oscillation, and Atlantic multidecadal oscillation during identified drought periods suggests that stochastic atmospheric variability can contribute significantly to the occurrence of simulated megadroughts in the NASW. These findings indicate that either an expanded paradigm is needed to understand multidecadal hydroclimate variability in the NASW or AOGCMs may incorrectly simulate the strength and/or dynamics of the connection between NASW hydroclimate variability and the tropical Pacific