Julien Emile-Geay

North American Droughts of the Last Millennium from a Gridded Network of Tree-Ring Data

Drought is the most economically expensive recurring natural disaster to strike North America in modern times. Recently available gridded drought reconstructions have been developed for most of North America from a network of drought-sensitive tree-ring chronologies, many of which span the last 1000 yr. These reconstructions enable the authors to put the famous droughts of the instrumental record (i.e., the 1930s Dust Bowl and the 1950s Southwest droughts) into the context of 1000 yr of natural drought variability on the continent. We can now, with this remarkable new record, examine the severity, persistence, spatial signatures, and frequencies of drought variability over the past milllennium, and how these have changed with time. The gridded drought reconstructions reveal the existence of successive “megadroughts,” unprecedented in persistence (20‐40 yr), yet similar in year-to-year severity and spatial distribution to the major droughts experienced in today’s North America. These megadroughts occurred during a 400-yr-long period in the early to middle second millennium A.D., with a climate varying as today’s, but around a drier mean. The implication is that the mechanism forcing persistent drought in the West and the Plains in the instrumental era is analagous to that underlying the megadroughts of the medieval period. The leading spatial mode of drought variability in the recontructions resembles the North American ENSO pattern: widespread drought across the United States, centered on the Southwest, with a hint of the opposite phase in the Pacific Northwest. Recently, climate models forced by the observed history of tropical Pacific SSTs have been able to successfully simulate all of the major North American droughts of the last 150 yr. In each case, cool “La Nina‐like” conditions in the tropical Pacific are consistent with North American drought. With ENSO showing a pronounced signal in the gridded drought recontructions of the last millennium, both in terms of its link to the leading spatial mode, and the leading time scales of drought variability (revealed by multitaper spectral analysis and wavelet analysis), it is postulated that, as for the modern day, the medieval megadroughts were forced by protracted La Nina‐like tropical Pacific SSTs. Further evidence for this comes from the global hydroclimatic “footprint” of the medieval era revealed by existing paleoclimatic archives from the tropical Pacific and ENSO-sensitive tropical and extratropical land regions. In general, this global pattern matches that observed for modern-day persistent North American drought, whereby a La Nina‐like tropical Pacific is accompanied by hemispheric, and in the midlatitudes, zonal, symmetry of hydroclimatic anomalies.

Volcanoes and ENSO over the Past Millennium

Abstract The controversial claim that El Nino events might be partially caused by radiative forcing due to volcanic aerosols is reassessed. Building on the work of Mann et al., estimates of volcanic forcing over the past millennium and a climate model of intermediate complexity are used to draw a diagram of El Nino likelihood as a function of the intensity of volcanic forcing. It is shown that in the context of this model, only eruptions larger than that of Mt. Pinatubo (1991, peak dimming of about 3.7 W m−2) can shift the likelihood and amplitude of an El Nino event above the level of the model’s internal variability. Explosive volcanism cannot be said to trigger El Nino events per se, but it is found to raise their likelihood by 50% on average, also favoring higher amplitudes. This reconciles, on one hand, the demonstration by Adams et al. of a statistical relationship between explosive volcanism and El Nino and, on the other hand, the ability to predict El Nino events of the last 148 yr without knowled...

Estimating Central Equatorial Pacific SST Variability over the Past Millennium. Part II: Reconstructions and Implications

AbstractReducing the uncertainties surrounding the impacts of anthropogenic climate change requires vetting general circulation models (GCMs) against long records of past natural climate variability. This is particularly challenging in the tropical Pacific Ocean, where short, sparse instrumental data preclude GCM validation on multidecadal to centennial time scales. This two-part paper demonstrates the application of two statistical methodologies to a network of accurately dated tropical climate records to reconstruct sea surface temperature (SST) variability in the Nino-3.4 region over the past millennium. While Part I described the methods and established their validity and limitations, this paper presents several reconstructions of Nino-3.4, analyzes their sensitivity to procedural choices and input data, and compares them to climate forcing time series and previously published tropical Pacific SST reconstructions. The reconstructions herein show remarkably similar behavior at decadal to multidecadal s...

Estimating Central Equatorial Pacific SST Variability over the Past Millennium. Part I: Methodology and Validation

AbstractConstraining the low-frequency (LF) behavior of general circulation models (GCMs) requires reliable observational estimates of LF variability. This two-part paper presents multiproxy reconstructions of Nino-3.4 sea surface temperature over the last millennium, applying two techniques [composite plus scale (CPS) and hybrid regularized expectation maximization (RegEM) truncated total least squares (TTLS)] to a network of tropical, high-resolution proxy records. This first part presents the data and methodology before evaluating their predictive skill using frozen network analysis (FNA) and pseudoproxy experiments. The FNA results suggest that about half of the Nino-3.4 variance can be reconstructed back to A.D. 1000, but they show little LF skill during certain intervals. More variance can be reconstructed in the interannual band where climate signals are strongest, but this band is affected by dating uncertainties (which are not formally addressed here). The CPS reliably estimates interannual varia...

PRYSM: An open‐source framework for PRoxY System Modeling, with applications to oxygen‐isotope systems

Paleoclimate observations constitute the only constraint on climate behavior prior to the instrumental era. However, such observations only provide indirect (proxy) constraints on physical variables. Proxy system models aim to improve the interpretation of such observations and better quantify their inherent uncertainties. However, existing models are currently scattered in the literature, making their integration difficult. Here, we present a comprehensive modeling framework for proxy systems, named PRYSM. For this initial iteration, we focus on water-isotope based climate proxies in ice cores, corals, tree ring cellulose, and speleothem calcite. We review modeling approaches for each proxy class, and pair them with an isotope-enabled climate simulation to illustrate the new scientific insights that may be gained from this framework. Applications include parameter sensitivity analysis, the quantification of archive-specific processes on the recorded climate signal, and the quantification of how chronological uncertainties affect signal detection, demonstrating the utility of PRYSM for a broad array of climate studies.

Assessing El Niño Southern Oscillation variability during the past millennium

evaluate the relative strength/frequency of El Nino and La Nina events. In contrast to previous studies, we use robust and resistant statistics to quantify the spread and symmetry of the d 18 O distributions; an approach motivated by the relatively small sample size and the presence of outliers. Furthermore, we use a pseudo‐proxy approach to investigate the effects of the different paleo‐environmental factors on the statistics of the d 18 O distributions, which could bias the paleo‐ENSO reconstruction. We find no systematic difference in the magnitude/strength of ENSO during the Northern Hemisphere MCA or LIA. However, our results suggest that ENSO during the MCA was skewed toward stronger/more frequent La Nina than El Nino, an observation consistent with the medieval megadroughts documented from sites in western North America.

A global multiproxy database for temperature reconstructions of the Common Era

Reproducible climate reconstructions of the Common Era (1 CE to present) are key to placing industrial-era warming into the context of natural climatic variability. Here we present a community-sourced database of temperature-sensitive proxy records from the PAGES2k initiative. The database gathers 692 records from 648 locations, including all continental regions and major ocean basins. The records are from trees, ice, sediment, corals, speleothems, documentary evidence, and other archives. They range in length from 50 to 2000 years, with a median of 547 years, while temporal resolution ranges from biweekly to centennial. Nearly half of the proxy time series are significantly correlated with HadCRUT4.2 surface temperature over the period 1850–2014. Global temperature composites show a remarkable degree of coherence between high- and low-resolution archives, with broadly similar patterns across archive types, terrestrial versus marine locations, and screening criteria. The database is suited to investigations of global and regional temperature variability over the Common Era, and is shared in the Linked Paleo Data (LiPD) format, including serializations in Matlab, R and Python.

Multidecadal rainfall variability in South Pacific Convergence Zone as revealed by stalagmite geochemistry

Pacifi c decadal variability (PDV) causes widespread, persistent fl uctuations that affect climate, water resources, and fi sheries throughout the Pacifi c basin, yet the magnitude, frequency, and causes of PDV remain poorly constrained. Here we present an absolutely dated, subannually resolved, 446 yr stable oxygen isotope (δ 18 O) cave record of rainfall variability in Vanuatu (southern Pacifi c Ocean), a location that has a climate heavily infl uenced by the South Pacifi c Convergence Zone (SPCZ). The δ 18 O-based proxy rainfall record is dominated by changes in stalagmite δ 18 O that are large (~1‰), quasi-periodic (~50 yr period), and generally abrupt (within 5‐10 yr). These isotopic changes imply abrupt rainfall changes of as much as ~1.8 m per wet season, changes that can be ~2.5◊ larger than the 1976 C.E. shift in rainfall amount associated with a PDV phase switch. The Vanuatu record also shares little commonality with previously documented changes in the Intertropical Convergence Zone during the Little Ice Age or solar forcing. We conclude that multidecadal SPCZ variability is likely of an endogenous nature. Large, spontaneous, and low-frequency changes in SPCZ rainfall during the past 500 yr have important implications for the relative magnitude of natural PDV possible in the coming century.

Statistical paleoclimate reconstructions via Markov random fields

Understanding centennial scale climate variability requires data sets that are accurate, long, continuous and of broad spatial coverage. Since instrumental measurements are generally only available after 1850, temperature fields must be reconstructed using paleoclimate archives, known as proxies. Various climate field reconstructions (CFR) methods have been proposed to relate past temperature to such proxy networks. In this work, we propose a new CFR method, called GraphEM, based on Gaussian Markov random fields embedded within an EM algorithm. Gaussian Markov random fields provide a natural and flexible framework for modeling high-dimensional spatial fields. At the same time, they provide the parameter reduction necessary for obtaining precise and well-conditioned estimates of the covariance structure, even in the sample-starved setting common in paleoclimate applications. In this paper, we propose and compare the performance of different methods to estimate the graphical structure of climate fields, and demonstrate how the GraphEM algorithm can be used to reconstruct past climate variations. The performance of GraphEM is compared to the widely used CFR method RegEM with regularization via truncated total least squares, using synthetic data. Our results show that GraphEM can yield significant improvements, with uniform gains over space, and far better risk properties. We demonstrate that the spatial structure of temperature fields can be well estimated by graphs where each neighbor is only connected to a few geographically close neighbors, and that the increase in performance is directly related to recovering the underlying sparsity in the covariance of the spatial field. Our work demonstrates how significant improvements can be made in climate reconstruction methods by better modeling the covariance structure of the climate field.

Fragility of reconstructed temperature patterns over the Common Era: Implications for model evaluation

Climate field reconstructions(CFRs) enable spatially resolved estimates of past climates, providing important insights about climate variability over the Common Era. In particular, a reconstructed “La Nina-like” pattern during the transition from the Medieval Climate Anomaly (MCA) to the Little Ice Age has been widely tied to medieval droughts in southwest North America. This pattern is now used as a key benchmark for global climate model simulations of the last millennium, which have yet to reproduce it. Here we test the patterns robustness by using four different CFR methods and two proxy networks. With the older network, we find the reconstructed patterns to be highly method-dependent, with the La Nina-like pattern not reproduced by two of the CFR methodologies. With the updated proxy network, a globally uniform MCA emerges with all methods, in agreement with simulations from the Paleoclimate Modelling Intercomparison Project Phase 3 ensemble. Our results caution against drawing dynamical interpretations from a single CFR and affirm the importance of developing CFRs through improved statistical methodology and community-driven proxy syntheses.