Search for: All records

We use online data assimilation to combine information from a linear inverse model of coupled atmosphere‐ocean dynamics with proxy records to create a new annual‐resolution reconstruction of atmosphere and ocean fields over the last millennium. Instrumental validation of reconstructed sea‐surface temperature and 0–700 m ocean heat content shows broad regions of positive spatial correlations, and high correlations (∼0.6–0.9) for global averages and indices of large‐scale modes of atmospheric variability. Compared to previous reconstructions, the online reconstructions show global and hemispheric averages with little‐to‐no millennial‐scale trend and global‐mean temperatures ∼0.25–0.5 K cooler during early periods (1000–1400 C.E.). The spatial anomaly differences of average temperature between an early (1000–1250 C.E.) and later (1400–1700 C.E.) period show warm anomalies over high‐latitude Europe and cool tropical conditions in partial agreement with previous assessments. The addition of online data assimilation, which provides dynamical memory to climate proxy information, is shown to be crucial for adequately characterizing decadal‐to‐centennial‐scale variability of 0–700 m ocean heat content. Furthermore, the climate forecasts provide model‐based physical constraints for atmosphere–ocean interaction, which become increasingly important during early periods when less proxy information is available for assimilation.

 

Paleoclimate reconstruction relies on estimates of spatiotemporal relationships among climate quantities to interpolate between proxy data. This work quantifies how structural uncertainties in those relationships translate to uncertainties in reconstructions of past climate. We develop and apply a data assimilation uncertainty quantification approach to paleoclimate networks and observational uncertainties representative of data for the last millennium. We find that structural uncertainties arising from uncertain spatial covariance relationships typically contribute 10% of the total uncertainty in reconstructed temperature variability at small (∼200 km), continental, and hemispheric length scales, with larger errors (50% or larger) in regions where long‐range climate covariances are least certain. These structural uncertainties contribute far more to errors in uncertainty quantification, sometimes by a factor of 5 or higher. Accounting for and reducing uncertainties in climate model dynamics and resulting covariance relationships will improve paleoclimate reconstruction accuracy.

 

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.

 

Despite the importance of interdecadal climate variability, we have a limited understanding of which geographic regions are associated with global temperature variability at these timescales. The instrumental record tends to be too short to develop sample statistics to study interdecadal climate variability, and Coupled Model Intercomparison Project, Phase 5 (CMIP5) climate models tend to disagree about which locations most strongly influence global mean interdecadal temperature variability. Here we use a new paleoclimate data assimilation product, the Last Millennium Reanalysis (LMR), to examine where local variability is associated with global mean temperature variability at interdecadal timescales. The LMR framework uses an ensemble Kalman filter data assimilation approach to combine the latest paleoclimate data and state‐of‐the‐art model data to generate annually resolved field reconstructions of surface temperature, which allow us to explore the timing and dynamics of preinstrumental climate variability in new ways. The LMR consistently shows that the middle‐ to high‐latitude north Pacific and the high‐latitude North Atlantic tend to lead global temperature variability on interdecadal timescales. These findings have important implications for understanding the dynamics of low‐frequency climate variability in the preindustrial era.

 

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. 

Abstract Scientific understanding of low-frequency tropical Pacific variability, especially responses to perturbations in radiative forcing, suffers from short observational records, sparse proxy networks, and bias in model simulations. Here, we combine the strengths of proxies and models through coral-based paleoclimate data assimilation. We combine coral archives ( δ 18 O, Sr/Ca) with the dynamics, spatial teleconnections, and intervariable relationships of the CMIP5/PMIP3 Past1000 experiments using the Last Millennium Reanalysis data assimilation framework. This analysis creates skillful reconstructions of tropical Pacific temperatures over the observational era. However, during the period of intense volcanism in the early nineteenth century, southwestern Pacific corals produce El Niño–Southern Oscillation (ENSO) reconstructions that are of opposite sign from those from eastern Pacific corals and tree ring records. We systematically evaluate the source of this discrepancy using 1) single-proxy experiments, 2) varied proxy system models (PSMs), and 3) diverse covariance patterns from the Past1000 simulations. We find that individual proxy records and coral PSMs do not significantly contribute to the discrepancy. However, following major eruptions, the southwestern Pacific corals locally record more persistent cold anomalies than found in the Past1000 experiments and canonical ENSO teleconnections to the southwest Pacific strongly control the reconstruction response. Furthermore, using covariance patterns independent of ENSO yields reconstructions consistent with coral archives across the Pacific. These results show that model bias can strongly affect how proxy information is processed in paleoclimate data assimilation. As we illustrate here, model bias influences the magnitude and persistence of the response of the tropical Pacific to volcanic eruptions. 

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.