Global warming is expected to increase the frequency and intensity of droughts in the twenty-first century, but the relative contributions from changes in moisture supply (precipitation) versus evaporative demand (potential evapotranspiration; PET) have not been comprehensively assessed. Using output from a suite of general circulation model (GCM) simulations from phase 5 of the Coupled Model Intercomparison Project, projected twentyfirst century drying and wetting trends are investigated using two offline indices of surface moisture balance: the Palmer Drought Severity Index (PDSI) and the Standardized Precipitation Evapotranspiration Index (SPEI). PDSI and SPEI projections using precipitation and PenmanMonteith based PET changes from the GCMs generally agree, showing robust cross-model drying in western North America, Central America, the Mediterranean, southern Africa, and the Amazon and robust wetting occurring in the Northern Hemisphere high latitudes and east Africa (PDSI only). The SPEI is more sensitive to PET changes than the PDSI, especially in arid regions such as the Sahara and Middle East. Regional drying and wetting patterns largely mirror the spatially heterogeneous response of precipitation in the models, although drying in the PDSI and SPEI calculations extends beyond the regions of reduced precipitation. This expansion of drying areas is attributed to globally widespread increases in PET, caused by increases in surface net radiation and the vapor pressure deficit.Increased PET not only intensifies drying in areas where precipitation is already reduced, it also drives areas into drought that would otherwise experience little drying or even wetting from precipitation trends alone. This PET amplification effect is largest in the Northern Hemisphere mid-latitudes, and is especially pronounced in western North America, Europe, and southeast China. Compared to PDSI projections using precipitation changes only, the projections incorporating both precipitation and PET changes increase the percentage of global land area projected to experience at least moderate drying (PDSI standard deviation of B-1) by the end of the twenty-first century from 12 to 30 %. PET induced moderate drying is even more severe in the SPEI projections (SPEI standard deviation of B-1; 11 to 44 %), although this is likely less meaningful because much of the PET induced drying in the SPEI occurs in the aforementioned arid regions. Integrated accounting of both the supply and demand sides of the surface moisture balance is therefore critical for characterizing the full range of projected drought risks tied to increasing greenhouse gases and associated warming of the climate system.

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.

Historical trends in the tropical Pacific zonal sea surface temperature gradient (SST gradient) are analyzed herein using 41 climate models (83 simulations) and 5 observational data sets. A linear inverse model is trained on each simulation and observational data set to assess if trends in the SST gradient are significant relative to the stationary statistics of internal variability, as would suggest an important role for external forcings such as anthropogenic greenhouse gasses. None of the 83 simulations have a positive trend in the SST gradient, a strengthening of the climatological SST gradient with more warming in the western than eastern tropical Pacific, as large as the mean trend across the five observational data sets. If the observed trends are anthropogenically forced, this discrepancy suggests that state‐of‐the‐art climate models are not capturing the observed response of the tropical Pacific to anthropogenic forcing, with serious implications for confidence in future climate projections. There are caveats to this interpretation, however, as some climate models have a significant strengthening of the SST gradient between 1900 and 2013 Common Era, though smaller in magnitude than the observational data sets, and the strengthening in three out of five observational data sets is insignificant. When combined with observational uncertainties and the possibility of centennial time scale internal variability not sampled by the linear inverse model, this suggests that confident validation of anthropogenic SST gradient trends in climate models will require further emergence of anthropogenic trends. Regardless, the differences in SST gradient trends between climate models and observational data sets are concerning and motivate the need for process‐level validation of the atmosphere‐ocean dynamics relevant to climate change in the tropical Pacific.

During the Medieval Climate Anomaly (MCA), Western North America experienced episodes of intense aridity that persisted for multiple decades or longer. These megadroughts are well documented in many proxy records, but the causal mechanisms are poorly understood. General circulation models (GCMs) simulate megadroughts, but do not reproduce the temporal clustering of events during the MCA, suggesting they are not caused by the time history of volcanic or solar forcing. Instead, GCMs generate megadroughts through (1) internal atmospheric variability, (2) sea-surface temperatures, and (3) land surface and dust aerosol feedbacks. While no hypothesis has been definitively rejected, and no GCM has accurately reproduced all features (e.g., timing, duration, and extent) of any specific megadrought, their persistence suggests a role for processes that impart memory to the climate system (land surface and ocean dynamics). Over the 21st century, GCMs project an increase in the risk of megadrought occurrence through greenhouse gas forced reductions in precipitation and increases in evaporative demand. This drying is robust across models and multiple drought indicators, but major uncertainties still need to be resolved. These include the potential moderation of vegetation evaporative losses at higher atmospheric [CO 2 ], variations in land surface model complexity, and decadal to multidecadal modes of natural climate variability that could delay or advance onset of aridification over the the next several decades. Because future droughts will arise from both natural variability and greenhouse gas forced trends in hydroclimate, improving our understanding of the natural drivers of persistent multidecadal megadroughts should be a major research priority. © 2016 Wiley Periodicals, Inc. How to cite this article:WIREs Clim Change 2016Change , 7:411-432. doi: 10.1002 INTRODUCTION R ecurrent droughts are a normal part of climate variability in Western North America, and recent events (e.g., California, the Southwest) have highlighted the vulnerability of people and ecosystems to the capricious nature of water availability in this region. [1][2][3][4] Despite these challenges, however, there is robust evidence from the paleoclimate record that Western North America experienced even worse droughts (megadroughts) over the last two 411millennia. 5,6 These events, many lasting multiple decades, had profound impacts on the contemporary indigenous societies, [7][8][9][10] vegetation, [11][12][13] and landscape. 12,14,15 Megadroughts are not precisely defined in the literature, but typically refer to persistent drought events in the preindustrial period with durations longer than a decade. These droughts stand in sharp contrast to the shorter length of more recent 20th-century events such as the Dust Bowl (1932Bowl ( -1939 and the 1950s drought (1948)(1949)(1950)(1951)(1952)(1953)(1954)(1955)(1956)(1957). Megadroughts are thus primarily differentiated from more recent instrumental-era droughts in terms of their duration and, occasionally, ...

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: 258-42.58N, 1258-1058W). 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.

Multidecadal drought periods in the North American Southwest (258-42.58N, 1258-1058W), 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 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 western United States was affected by several megadroughts during the last 1200 years, most prominently during the Medieval Climate Anomaly (MCA; 800 to 1300 CE). A null hypothesis is developed to test the possibility that, given a sufficiently long period of time, these events are inevitable and occur purely as a consequence of internal climate variability. The null distribution of this hypothesis is populated by a linear inverse model (LIM) constructed from global sea surface temperature anomalies and self-calibrated Palmer drought severity index data for North America. Despite being trained only on seasonal data from the late twentieth century, the LIM produces megadroughts that are comparable in their duration, spatial scale, and magnitude to the most severe events of the last 12 centuries. The null hypothesis therefore cannot be rejected with much confidence when considering these features of megadrought, meaning that similar events are possible today, even without any changes to boundary conditions. In contrast, the observed clustering of megadroughts in the MCA, as well as the change in mean hydroclimate between the MCA and the 1500–2000 period, are more likely to have been caused by either external forcing or by internal climate variability not well sampled during the latter half of the twentieth century. Finally, the results demonstrate that the LIM is a viable tool for determining whether paleoclimate reconstructions events should be ascribed to external forcings or to “out of sample” climate mechanisms, or if they are consistent with the variability observed during the recent period.

The temporal stationarity of the teleconnection between the tropical Pacific Ocean and North America (NA) is analyzed in atmosphere‐only, and coupled last‐millennium, historical, and control runs from the Coupled Model Intercomparison Project Phase 5 data archive. The teleconnection, defined as the correlation between December‐January‐February (DJF) tropical Pacific sea surface temperatures (SSTs) and DJF 200 mb geopotential height, is found to be nonstationary on multidecadal timescales. There are significant changes in the spatial features of the teleconnection over NA in continuous 56‐year segments of the last millennium and control simulations. Analysis of atmosphere‐only simulations forced with observed SSTs indicates that atmospheric noise cannot account for the temporal variability of the teleconnection, which instead is likely explained by the strength of, and multidecadal changes in, tropical Pacific Ocean variability. These results have implications for teleconnection‐based analyses of model fidelity in simulating precipitation, as well as any reconstruction and forecasting efforts that assume stationarity of the observed teleconnection.