Comparison of climate field reconstruction techniques: Application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.

Establishing the skill of climate field reconstruction techniques for precipitation with pseudoproxy experiments

This study aims at assessing the skill of several climate field reconstruction techniques (CFR) to reconstruct past precipitation over continental Europe and the Mediterranean at seasonal time scales over the last two millennia from proxy records. A number of pseudoproxy experiments are performed within the virtual reality ofa regional paleoclimate simulation at 45 km resolution to analyse different aspects of reconstruction skill. Canonical Correlation Analysis (CCA), two versions of an Analog Method (AM) and Bayesian hierarchical modeling (BHM) are applied to reconstruct precipitation from a synthetic network of pseudoproxies that are contaminated with various types of noise. The skill of the derived reconstructions is assessed through comparison with precipitation simulated by the regional climate model. Unlike BHM, CCA systematically underestimates the variance. The AM can be adjusted to overcome this shortcoming, presenting an intermediate behaviour between the two aforementioned techniques. However, a trade-off between reconstruction-target correlations and reconstructed variance is the drawback of all CFR techniques. CCA (BHM) presents the largest (lowest) skill in preserving the temporal evolution, whereas the AM can be tuned to reproduce better correlation at the expense of losing variance. While BHM has been shown to perform well for temperatures, it relies heavily on prescribed spatial correlation lengths. While this assumption is valid for temperature, it is hardly warranted for precipitation. In general, none of the methods outperforms the other. All experiments agree that a dense and regularly distributed proxy network is required to reconstruct precipitation accurately, reflecting its high spatial and temporal variability. This is especially true in summer, when a specifically short de-correlation distance from the proxy location is caused by localised summertime convective precipitation events.

Circulation dynamics and its influence on European and Mediterranean January-April climate over the past half millennium: results and insights from instrumental data, documentary evidence and coupled climate models

We use long instrumental temperature series together with available field reconstructions of sea-level pressure (SLP) and three-dimensional climate model simulations to analyze relations between temperature anomalies and atmospheric circulation patterns over much of Europe and the Mediterranean for the late winter/early spring (January–April, JFMA) season. A Canonical Correlation Analysis (CCA) investigates interannual to interdecadal covariability between a new gridded SLP field reconstruction and seven long instrumental temperature series covering the past 250 years. We then present and discuss prominent atmospheric circulation patterns related to anomalous warm and cold JFMA conditions within different European areas spanning the period 1760–2007. Next, using a data assimilation technique, we link gridded SLP data with a climate model (EC-Bilt-Clio) for a better dynamical understanding of the relationship between large scale circulation and European climate. We thus present an alternative approach to reconstruct climate for the pre-instrumental period based on the assimilated model simulations. Furthermore, we present an independent method to extend the dynamic circulation analysis for anomalously cold European JFMA conditions back to the sixteenth century. To this end, we use documentary records that are spatially representative for the long instrumental records and derive, through modern analogs, large-scale SLP, surface temperature and precipitation fields. The skill of the analog method is tested in the virtual world of two three-dimensional climate simulations (ECHO-G and HadCM3). This endeavor offers new possibilities to both constrain climate model into a reconstruction mode (through the assimilation approach) and to better asses documentary data in a quantitative way.

Southern Hemisphere high-resolution palaeoclimate records of the last 2000 years

This study presents a comprehensive assessment of high-resolution Southern Hemisphere (SH) paleoarchives covering the last 2000 years. We identified 174 monthly to annually resolved climate proxy (tree ring, coral, ice core, documentary, speleothem and sedimentary) records from the Hemisphere. We assess the interannual and decadal sensitivity of each proxy record to large-scale circulation indices from the Pacific, Indian and Southern Ocean regions over the twentieth century. We then analyse the potential of this newly expanded palaeoclimate network to collectively represent predictands (sea surface temperature, sea level pressure, surface air temperature and precipitation) commonly used in climate reconstructions. The key dynamical centres-of-action of the equatorial Indo-Pacific are well captured by the palaeoclimate network, indicating that there is considerable reconstruction potential in this region, particularly in the post AD 1600 period when a number of long coral records are available. Current spatiotemporal gaps in data coverage and regions where significant potential for future proxy collection exists are discussed. We then highlight the need for new and extended records from key dynamical regions of the Southern Hemisphere. Although large-scale climate field reconstructions for the SH are in their infancy, we report that excellent progress in the development of regional proxies now makes plausible estimates of continental- to hemispheric-scale climate variations possible.

Alpine lacustrine varved record reveals summer temperature as main control of glacier fluctuations over the past 2250 years

Glacier fluctuations are a key indicator of changing climate. Their reconstruction beyond historical times unravels glacier variability and its forcing factors on long time scales, which can considerably improve our understanding of the climate–glacier relationship. Here, we present a 2250-year-long reconstruction of particle-mass accumulation rates recorded in the lacustrine sediments of Lake Trüebsee (Central Swiss Alps) that are directly related to glacier extent, thus reflecting a continuous record of fluctuations of the upstream-located Titlis Glacier. Mass accumulation rate values show strong centennial to multi-centennial fluctuations and reveal 12 well-pronounced periods of enhanced values corresponding to times of maximum extent of the neighboring Lower Grindelwald Glacier. This result supports previous studies of proglacial lake sediments that documented high mass accumulation rate values during glacier advances. The strong variability in the Lake Trüebsee mass accumulation rate record thus represents a highly sensitive paleoclimatic archive, which mirrors rapid and pronounced feedbacks of Titlis Glacier to climatic changes over the past 2250years. The comparison of our data with independent paleo-temperature reconstructions from tree rings suggests that variations in mean summer temperature were the primary driving factor of fluctuations of Titlis Glacier. Also, advances of Titlis Glacier occurred during the grand solar minima (Dalton, Maunder, Spörer, Wolf) of the last millennium. This relation of glacier extent with summer temperature reveals strong evidence that the mass balance of this Alpine glacier is primarily controlled by the intensity of glacier melting during summer.

Glacier mass balance reconstruction by sublimation induced enrichment of chemical species on Cerro Tapado (Chilean Andes)

A 36 m long ice core down to bedrock from the Cerro Tapado glacier (5536 m a.s.l, 30°08' S, 69°55' W) was analyzed to reconstruct past climatic conditions for Northern Chile. Because of the marked seasonality in the precipitation (short wet winter and extended dry summer periods) in this region, major snow ablation and related post-depositional processes occur on the glacier surface during summer periods. They include predominantly sublimation and dry deposition. Assuming that, like measured during the field campaign, the enrichment of chloride was always related to sublimation, the chemical record along the ice core may be applied to reconstruct the history of such secondary processes linked to the past climatic conditions over northern Chile. For the time period 1962–1999, a mean annual net accumulation of 316 mm water equivalent (weq) and 327 mm weq loss by sublimation was deduced by this method. This corresponds to an initial total annual accumulation of 539 mm weq. The annual variability of the accumulation and sublimation is related with the Southern Oscillation Index (SOI): higher net-accumulation during El-Niño years and more sublimation during La Niña years. The deepest part of the ice record shows a time discontinuity; with an ice body deposited under different climatic conditions: 290 mm higher precipitation but with reduced seasonal distribution (+470 mm in winter and –180 mm in summer) and –3°C lower mean annual temperature. Unfortunately, its age is unknown. The comparison with regional proxy data however let us conclude that the glacier buildup did most likely occur after the dry mid-Holocene.