Pan-Continental Droughts in North America over the Last Millennium

Regional droughts are common in North America, but pan-continental droughts extending across multiple regions, including the 2012 event, are rare relative to single-region events. Here, the tree-ring-derived North American Drought Atlas is used to investigate drought variability in four regions over the last millennium, focusing on pan-continental droughts. During the Medieval Climate Anomaly (MCA), the central plains (CP), Southwest (SW), and Southeast (SE) regions experienced drier conditions and increased occurrence of droughts and the Northwest (NW) experienced several extended pluvials. Enhanced MCA aridity in the SW and CP manifested as multidecadal megadroughts. Notably, megadroughts in these regions differed in their timing and persistence, suggesting that they represent regional events influenced by local dynamics rather than a unified, continental-scale phenomena. There is no trend in pan-continental drought occurrence, defined as synchronous droughts in three or more regions. SW, CP, and SE (SW+CP+SE) droughts are the most common, occurring in 12 percent of all years and peaking in prevalence during the twelfth and thirteenth centuries; patterns involving three other regions occur in about 8 percent of years. Positive values of the Southern Oscillation index (La Nina conditions) are linked to SW, CP, and SE (SW+CP+SE) droughts and SW, CP, and NW (SW+CP+NW) droughts, whereas CP, NW, and SE (CP+NW+SE) droughts are associated with positive values of the Pacific decadal oscillation and Atlantic multidecadal oscillation. While relatively rare, pan-continental droughts are present in the paleo record and are linked to defined modes of climate variability, implying the potential for seasonal predictability. Assuming stable drought teleconnections, these events will remain an important feature of future North American hydroclimate, possibly increasing in their severity in step with other expected hydroclimate responses to increased greenhouse gas forcing

Bridging Past and Future Climate across Paleoclimatic Reconstructions, Observations, and Models: A Hydroclimate Case Study

Potential biases in tree-ring reconstructed Palmer drought severity index (PDSI) are evaluated using Thornthwaite (TH), Penman–Monteith (PM), and self-calibrating Penman–Monteith (SC) PDSI in three diverse regions of the United States and tree-ring chronologies from the North American drought atlas (NADA). Minimal differences are found between the three PDSI reconstructions and all compare favorably to independently reconstructed Thornthwaite-based PDSI from the NADA. Reconstructions are bridged with model-derived PDSI_TH and PDSI_PM, which both closely track modeled soil moisture (near surface and full column) during the twentieth century. Differences between modeled moisture-balance metrics only emerge in twenty-first-century projections. These differences confirm the tendency of PDSI_TH to overestimate drying when temperatures exceed the range of the normalization interval; the more physical accounting of PDSI_PM compares well with modeled soil moisture in the projection interval. Remaining regional differences in the secular behavior of projected soil moisture and PDSI_PM are interpreted in terms of underlying physical processes and temporal sampling. Results demonstrate the continued utility of PDSI as a metric of surface moisture balance while additionally providing two recommendations for future work: 1) PDSI_PM (or similar moisture-balance metrics) compare well to modeled soil moisture and are an appropriate means of representing soil-moisture balance in model simulations and 2) although PDSI_PM is more physically appropriate than PDSI_TH, the latter metric does not bias tree-ring reconstructions of past hydroclimate variability and, as such, reconstructions targeting PDSI_TH can be used with confidence in data–model comparisons. These recommendations and the collective results of this study thus provide a framework for comparing hydroclimate variability within paleoclimatic, observational, and modeled data

REPLY TO COMMENTS OF NOLAN AND COOK

We appreciate and would like to respond to the comments made by Nolan and Cook. We make three points in our reply, which are directed at clarifying our position and responding to a few of Nolan and Cook\u27s assumptions. First, as a matter of clarification, the Palmer Drought Severity Index (PDSI) is not a measure of summer rainfall as Nolan and Cook (2010) imply. Instead it is a theoretical measure of soil-moisture, the value of which evolves over several months in response to fluxes of precipitation, evaporation, and runoff. In practice, PDSI is usually heavily weighted toward the precipitation side of the soil water balance. Second, with regard to archaeological facts, the known developmental trajectory of the greater Cahokia region is not quite as Nolan and Cook describe it. In their discussion of the timing of the precocious development of social complexity in the American Bottom, Nolan and Cook refer to the American Bottom during both the Edelhardt (A.D. 1000-1050) and Lohmann phases (A.D. 1050-1100) as being relatively wet but not the wettest areas in the Upper Mississippi and Ohio Valleys. However, Cahokia\u27s big bang dates to the end of the Edelhardt phase and continued through the Lohmann phase, so the overall climatic state of the Edelhardt phase is not at issue. Nolan and Cook\u27s Figure 2 supports the concept that the Cahokia area was extremely wet during the Lohmann phase as opposed to the Edelhardt phase

Forced and unforced variability of twentieth century North American droughts and pluvials

Research on the forcing of drought and pluvial events over North America is dominated by general circulation model experiments that often have operational limitations (e.g., computational expense, ability to simulate relevant processes, etc). We use a statistically based modeling approach to investigate sea surface temperature (SST) forcing of the twentieth century pluvial (1905-1917) and drought (1932-1939, 1948-1957, 1998-2002) events. A principal component (PC) analysis of Palmer Drought Severity Index (PDSI) from the North American Drought Atlas separates the drought variability into five leading modes accounting for 62% of the underlying variance. Over the full period spanning these events (1900-2005), the first three PCs significantly correlate with SSTs in the equatorial Pacific (PC 1), North Pacific (PC 2), and North Atlantic (PC 3), with spatial patterns (as defined by the empirical orthogonal functions) consistent with our understanding of North American drought responses to SST forcing. We use a large ensemble statistical modeling approach to determine how successfully we can reproduce these drought/pluvial events using these three modes of variability. Using Pacific forcing only (PCs 1-2), we are able to reproduce the 1948-1957 drought and 1905-1917 pluvial above a 95% random noise threshold in over 90% of the ensemble members; the addition of Atlantic forcing (PCs 1-2-3) provides only marginal improvement. For the 1998-2002 drought, Pacific forcing reproduces the drought above noise in over 65% of the ensemble members, with the addition of Atlantic forcing increasing the number passing to over 80%. The severity of the drought, however, is underestimated in the ensemble median, suggesting this drought intensity can only be achieved through internal variability or other processes. Pacific only forcing does a poor job of reproducing the 1932-1939 drought pattern in the ensemble median, and less than one third of ensemble members exceed the noise threshold (28%). Inclusion of Atlantic forcing improves the ensemble median drought pattern and nearly doubles the number of ensemble members passing the noise threshold (52%). Even with the inclusion of Atlantic forcing, the intensity of the simulated 1932-1939 drought is muted, and the drought itself extends too far into the southwest and southern Great Plains. To an even greater extent than the 1998-2002 drought, these results suggest much of the variance in the 1932-1939 drought is dependent on processes other than SST forcing. This study highlights the importance of internal noise and non SST processes for hydroclimatic variability over North America, complementing existing research using general circulation models

Inter-decadal climate variability in the Southern Hemisphere: evidence from Tasmanian tree rings over the past three millennia

EXTRACT (SEE PDF FOR FULL ABSTRACT): The characterization of inter-decadal climate variability in the Southern Hemisphere is severely constrained by the shortness of the instrumental climate records. To help relieve this constraint, we have developed and analyzed a reconstruction of warm-season (November-April) temperatures from Tasmanian tree rings that now extends back to 800 BC. A detailed analysis of this reconstruction in the time and frequency domains indicates that much of the inter-decadal variability is principally confined to four frequency bands with mean periods of 31, 57, 77, and 200 years. ... In so doing, we show how a future greenhouse warming signal over Tasmania could be masked by these natural oscillations unless they are taken into account

In Memoriam: Keith R. Briffa, 1952–2017

Keith R. Briffa was one of the most influential palaeoclimatologists of the last 30 years. His primary research interests lay in Late-Holocene climate change with a geographical emphasis on northern Eurasia. His greatest impact was in the field of dendroclimatology, a field that he helped to shape. His contributions have been seminal to the development of sound methods for tree-ring analysis and in their proper application to allow the interpretation of climate variability from tree rings. This led to the development of many important records that allow us to understand natural climate variability on timescales from years to millennia and to set recent climatic trends in their historical context

A Tree-Ring-Based Reconstruction of Delaware River Basin Streamflow Using Hierarchical Bayesian Regression

A hierarchical Bayesian regression model is presented for reconstructing the average summer streamflow at five gauges in the Delaware River basin using eight regional tree-ring chronologies. The model provides estimates of the posterior probability distribution of each reconstructed streamflow series considering parameter uncertainty. The vectors of regression coefficients are modeled as draws from a common multivariate normal distribution with unknown parameters estimated as part of the analysis. This leads to a multilevel structure. The covariance structure of the streamflow residuals across sites is explicitly modeled. The resulting partial pooling of information across multiple stations leads to a reduction in parameter uncertainty. The effect of no pooling and full pooling of station information, as end points of the method, is explored. The no-pooling model considers independent estimation of the regression coefficients for each streamflow gauge with respect to each tree-ring chronology. The full-pooling model considers that the same regression coefficients apply across all streamflow sites for a particular tree-ring chronology. The cross-site correlation of residuals is modeled in all cases. Performance on metrics typically used by tree-ring reconstruction experts, such as reduction of error, coefficient of efficiency, and coverage rates under credible intervals is comparable to, or better, for the partial-pooling model relative to the no-pooling model, and streamflow estimation uncertainty is reduced. Long record simulations from reconstructions are used to develop estimates of the probability of duration and severity of droughts in the region. Analysis of monotonic trends in the reconstructed drought events do not reject the null hypothesis of no trend at the 90% significance over 1754–2000

A reconstruction of global hydroclimate and dynamical variables over the Common Era

Hydroclimate extremes critically affect human and natural systems, but there remain many unanswered questions about their causes and how to interpret their dynamics in the past and in climate change projections. These uncertainties are due, in part, to the lack of long-term, spatially resolved hydroclimate reconstructions and information on the underlying physical drivers for many regions. Here we present the first global reconstructions of hydroclimate and associated climate dynamical variables over the past two thousand years. We use a data assimilation approach tailored to reconstruct hydroclimate that optimally combines 2,978 paleoclimate proxy-data time series with the physical constraints of an atmosphere—ocean climate model. The global reconstructions are annually or seasonally resolved and include two spatiotemporal drought indices, near-surface air temperature, an index of North Atlantic variability, the location of the intertropical convergence zone, and monthly Niño indices. This database, called the Paleo Hydrodynamics Data Assimilation product (PHYDA), will provide a critical new platform for investigating the causes of past climate variability and extremes, while informing interpretations of future hydroclimate projections

A comparison of times series approaches for dendroecological reconstructions of past canopy disturbance events

Time series analysis can identify outliers in tree-ring widths that may not only indicate past disturbances, but may also estimate the subsequent effects of these disturbances on tree growth. Finding a way to isolate these disturbance signals from tree-ring time series could have broad applications in forest ecology and management. Time series outliers may be expressed as pulse, step, or trend interventions, but few dendroecological studies have explored how well these different types of interventions express the response of tree-ring widths to a canopy disturbance resulting in a release event. This study addresses that question by comparing two different time series approaches for detecting release events: a pulse intervention approach and a new combined step and trend (CST) intervention approach. These methods are tested against tree-ring collections with known historical canopy disturbance events: northern red oaks (Quercus rubra) in New York released by the chestnut blight during the early 1920s, eastern hemlocks (Tsuga canadensis) in Pennsylvania released by adjacent selective logging in 1910, and northern red oaks and chestnut oaks (Quercus montana) in West Virginia released through an experimental thinning in 1982. Clusters of CST interventions, but not pulse interventions, were detected for all three collections during and immediately after the known disturbance events, showing that a CST intervention approach consistently reconstructs these release events. In addition, a CST intervention approach isolated canopy disturbance signals from tree-ring widths as disturbance-growth indices. Detrending disturbances from tree-ring widths provides an alternative approach to reconstruct climate in closed-canopy forests; however, just as importantly, disturbance-growth indices created through this method can also reconstruct changes in tree growth rates, biomass, or carbon resulting from a past disturbance event or forest thinning