skip to main content


Title: Do Multi‐Model Ensembles Improve Reconstruction Skill in Paleoclimate Data Assimilation?
Abstract

Reconstructing past climates remains a difficult task because pre‐instrumental observational networks are composed of geographically sparse and noisy paleoclimate proxy records that require statistical techniques to inform complete climate fields. Traditionally, instrumental or climate model statistical relationships are used to spread information from proxy measurements to other locations and to other climate variables. Here ensembles drawn from single climate models and from combinations of multiple climate models are used to reconstruct temperature variability over the last millennium in idealized experiments. We find that reconstructions derived from multi‐model ensembles produce lower error than reconstructions from single‐model ensembles when reconstructing independent model and instrumental data. Specifically, we find the largest decreases in error over regions far from proxy locations that are often associated with large uncertainties in model physics, such as mid‐ and high‐latitude ocean and sea‐ice regions. Furthermore, we find that multi‐model ensemble reconstructions outperform single‐model reconstructions that use covariance localization. We propose that multi‐model ensembles could be used to improve paleoclimate reconstructions in time periods beyond the last millennium and for climate variables other than air temperature, such as drought metrics or sea ice variables.

 
more » « less
Award ID(s):
1702423
NSF-PAR ID:
10374358
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Earth and Space Science
Volume:
8
Issue:
4
ISSN:
2333-5084
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract We use theNorthern Hemisphere Tree-RingNetwork Development (NTREND) tree-ring database to examine the effects of using a small, highly-sensitive proxy network for paleotemperature data assimilation over the last millennium. We first evaluate our methods using pseudo-proxy experiments. These indicate that spatial assimilations using this network are skillful in the extratropical Northern Hemisphere and improve on previous NTREND reconstructions based on Point-by-Point regression. We also find our method is sensitive to climate model biases when the number of sites becomes small. Based on these experiments, we then assimilate the real NTREND network. To quantify model prior uncertainty, we produce 10 separate reconstructions, each assimilating a different climate model. These reconstructions are most dissimilar prior to 1100 CE, when the network becomes sparse, but show greater consistency as the network grows. Temporal variability is also underestimated before 1100 CE. Our assimilation method produces spatial uncertainty estimates and these identify treeline North America and eastern Siberia as regions that would most benefit from development of new millennial-length temperature-sensitive tree-ring records. We compare our multi-model mean reconstruction to five existing paleo-temperature products to examine the range of reconstructed responses to radiative forcing. We find substantial differences in the spatial patterns and magnitudes of reconstructed responses to volcanic eruptions and in the transition between the Medieval epoch and Little Ice Age. These extant uncertainties call for the development of a paleoclimate reconstruction intercomparison framework for systematically examining the consequences of proxy network composition and reconstruction methodology and for continued expansion of tree-ring proxy networks. 
    more » « less
  2. Abstract. Reconstructions of past temperature and precipitation are fundamental to modeling the Greenland Ice Sheet and assessing its sensitivity to climate. Paleoclimate information is sourced from proxy records and climate-model simulations; however, the former are spatially incomplete while the latter are sensitive to model dynamics and boundary conditions. Efforts to combine these sources of information to reconstruct spatial patterns of Greenland climate over glacial–interglacial cycles have been limited by assumptions of fixed spatial patterns and a restricted use of proxy data. We avoid these limitations by using paleoclimate data assimilation to create independent reconstructions of mean-annual temperature and precipitation for the last 20 000 years. Our method uses oxygen isotope ratios of ice and accumulation rates from long ice-core records and extends this information to all locations across Greenland using spatial relationships derived from a transient climate-model simulation. Standard evaluation metrics for this method show that our results capture climate at locations without ice-core records. Our results differ from previous work in the reconstructed spatial pattern of temperature change during abrupt climate transitions; this indicates a need for additional proxy data and additional transient climate-model simulations. We investigate the relationship between precipitation and temperature, finding that it is frequency dependent and spatially variable, suggesting that thermodynamic scaling methods commonly used in ice-sheet modeling are overly simplistic. Our results demonstrate that paleoclimate data assimilation is a useful tool for reconstructing the spatial and temporal patterns of past climate on timescales relevant to ice sheets. 
    more » « less
  3. null (Ed.)
    Abstract. The Last Millennium Reanalysis (LMR) utilizes an ensemble methodology to assimilate paleoclimate data for the production of annually resolved climate field reconstructions of the Common Era. Two key elements are the focus of this work: the set of assimilated proxy records and the forward models that map climate variables to proxy measurements. Results based on an updated proxy database and seasonal regression-based forward models are compared to the LMR prototype, which was based on a smaller set of proxy records and simpler proxy models formulated as univariate linear regressions against annual temperature. Validation against various instrumental-era gridded analyses shows that the new reconstructions of surface air temperature and 500 hPa geopotential height are significantly improved (from 10 % to more than 100 %), while improvements in reconstruction of the Palmer Drought Severity Index are more modest. Additional experiments designed to isolate the sources of improvement reveal the importance of the updated proxy records, including coral records for improving tropical reconstructions, and tree-ring density records for temperature reconstructions, particularly in high northern latitudes. Proxy forward models that account for seasonal responses, and dependence on both temperature and moisture for tree-ring width, also contribute to improvements in reconstructed thermodynamic and hydroclimate variables in midlatitudes. The variability of temperature at multidecadal to centennial scales is also shown to be sensitive to the set of assimilated proxies, especially to the inclusion of primarily moisture-sensitive tree-ring-width records. 
    more » « less
  4. Waitt RB, Thackray GD (Ed.)
    Mountain glacier moraine sequences and their chronologies allow us to evaluate the timing and climate conditions that underpin changes in the equilibrium line altitudes (ELAs), which can provide valuable information on the paleoclimatology of understudied regions such as tropical East Africa. However, moraine sequences are inherently discontinuous, and the precise climate conditions that they represent can be ambiguous due to the sensitivity of mountain glaciers to temperature, precipitation, and other environmental variables. Here, we used a two-dimensional (2-D) iceflow and mass-balance model to simulate glacier extents and ELAs in the Rwenzori Mountains in East Africa over the past 31,000 yr (31 k.y.), including the Last Glacial Maximum (LGM), late glacial period, and the Holocene Epoch. We drove the glacier model with two independent, continuous temperature reconstructions to simulate possible glacier length changes through time. Model input paleoclimate values came from branched glycerol dialkyl glycerol tetraether (brGDGT) temperature reconstructions from alpine lakes on Mount Kenya for the last ~31 k.y., and precipitation reconstructions for the LGM came from various East African locations. We then compared the simulated fluctuations with the positions and ages (where known) of the Rwenzori moraines. The simulated glacier extents reached within 1.1 km of the dated LGM moraines in one valley (93% of the full LGM extent) when forced by the brGDGT temperature reconstructions (maximum cooling of 6.1 °C) and a decrease in precipitation (-10% than modern amounts). These simulations suggest that the Rwenzori glaciers required a cooling of at least 6.1 °C to reach the dated LGM moraines. Based on the model output, we predict an age of 12–11 ka for moraines located halfway between the LGM and modern glacier extents. We also predict ice-free conditions in the Rwenzori Mountains for most of the early to middle Holocene, followed by a late Holocene glacier readvance within the last 2000 yr. 
    more » « less
  5. Abstract

    Large uncertainties exist in climate model projections of the Asian summer monsoon (ASM). The El Niño‐Southern Oscillation (ENSO) is an important modulator of the ASM, but the ENSO‐ASM teleconnection is not stationary. Furthermore, teleconnections between ENSO and the East Asian versus South Asian subcomponents of the ASM exhibit distinct characteristics. Therefore, understanding the variability of the ENSO‐ASM teleconnection is critical for anticipating future variations in ASM intensity. To this end, we here use paleoclimate records to extend temporal coverage beyond the instrumental era by millennia. Recently, data assimilation techniques have been applied for the last millennium, which facilitates physically consistent, globally gridded climate reconstructions informed by paleoclimate observations. We use these novel data assimilation products to investigate variations in the ENSO‐ASM relationship over the last 1,000 years. We find that correlations between ENSO and ASM indices are mostly negative in the last millennium, suggesting that strong ASM years are often associated with La Niña events. During periods of weak correlations between ENSO and the East Asian summer monsoon, we observe an El Niño‐like sea surface temperature (SST) pattern in the Pacific. Additionally, SST patterns associated with periods of weak correlations between ENSO and South Asian summer monsoon rainfall are not consistent among data assimilation products. This underscores the importance of developing more precipitation‐sensitive paleoclimate proxies in the Indian subcontinental realm over the last millennium. Our study serves as a baseline for future appraisals of paleoclimate assimilation products and an example of informing our understanding of decadal‐scale ENSO‐ASM teleconnection variability using paleoclimate data sets.

     
    more » « less