Severe and persistent 21st-century drought in southwestern North America (SWNA) motivates comparisons to medieval megadroughts and questions about the role of anthropogenic climate change. We use hydrological modeling and new 1200-year tree-ring reconstructions of summer soil moisture to demonstrate that the 2000–2018 SWNA drought was the second driest 19-year period since 800 CE, exceeded only by a late-1500s megadrought. The megadrought-like trajectory of 2000–2018 soil moisture was driven by natural variability superimposed on drying due to anthropogenic warming. Anthropogenic trends in temperature, relative humidity, and precipitation estimated from 31 climate models account for 47% (model interquartiles of 35 to 105%) of the 2000–2018 drought severity, pushing an otherwise moderate drought onto a trajectory comparable to the worst SWNA megadroughts since 800 CE.

The tree-ring-based North American Drought Atlas (NADA), Monsoon Asia Drought Atlas (MADA), and Old World Drought Atlas (OWDA) collectively yield a near-hemispheric gridded reconstruction of hydroclimate variability over the last millennium. To test the robustness of the large-scale representation of hydroclimate variability across the drought atlases, the joint expression of seasonal climate variability and teleconnections in the NADA, MADA, and OWDA are compared against two global, observation-based PDSI products. Predominantly positive (negative) correlations are determined between seasonal precipitation (surface air temperature) and collocated tree-ring-based PDSI, with average Pearson’s correlation coefficients increasing in magnitude from boreal winter to summer. For precipitation, these correlations tend to be stronger in the boreal winter and summer when calculated for the observed PDSI record, while remaining similar for temperature. Notwithstanding these differences, the drought atlases robustly express teleconnection patterns associated with El Niño–Southern Oscillation (ENSO), the North Atlantic Oscillation (NAO), the Pacific decadal oscillation (PDO), and the Atlantic multidecadal oscillation (AMO). These expressions exist in the drought atlas estimates of boreal summer PDSI despite the fact that these modes of climate variability are dominant in boreal winter, with the exception of the AMO. ENSO and NAO teleconnection patterns in the drought atlases are particularly consistent with their well-known dominant expressions in boreal winter and over the OWDA domain, respectively. Collectively, the findings herein confirm that the joint Northern Hemisphere drought atlases robustly reflect large-scale patterns of hydroclimate variability on seasonal to multidecadal time scales over the twentieth century and are likely to provide similarly robust estimates of hydroclimate variability prior to the existence of widespread instrumental data.

Abstract. Water availability is fundamental to societies and ecosystems, but our understanding of variations in hydroclimate (including extreme events, flooding, and decadal periods of drought) is limited because of a paucity of modern instrumental observations that are distributed unevenly across the globe and only span parts of the 20th and 21st centuries. Such data coverage is insufficient for characterizing hydroclimate and its associated dynamics because of its multidecadal to centennial variability and highly regionalized spatial signature. High-resolution (seasonal to decadal) hydroclimatic proxies that span all or parts of the Common Era (CE) and paleoclimate simulations from climate models are therefore important tools for augmenting our understanding of hydroclimate variability. In particular, the comparison of the two sources of information is critical for addressing the uncertainties and limitations of both while enriching each of their interpretations. We review the principal proxy data available for hydroclimatic reconstructions over the CE and highlight the contemporary understanding of how these proxies are interpreted as hydroclimate indicators. We also review the available last-millennium simulations from fully coupled climate models and discuss several outstanding challenges associated with simulating hydroclimate variability and change over the CE. A specific review of simulated hydroclimatic changes forced by volcanic events is provided, as is a discussion of expected improvements in estimated radiative forcings, models, and their implementation in the future. Our review of hydroclimatic proxies and last-millennium model simulations is used as the basis for articulating a variety of considerations and best practices for how to perform proxy-model comparisons of CE hydroclimate. This discussion provides a framework for how best to evaluate hydroclimate variability and its associated dynamics using these comparisons and how they can better inform interpretations of both proxy data and model simulations. We subsequently explore means of using proxy-model comparisons to better constrain and characterize future hydroclimate risks. This is explored specifically in the context of several examples that demonstrate how proxy-model comparisons can be used to quantitatively constrain future hydroclimatic risks as estimated from climate model projections.