The future of the past—an earth system framework for high resolution paleoclimatology: editorial essay

High-resolution paleoclimatology is the study of climate variability and change on interannual to multi-century time scales. Its primary focus is the past few millennia, a period lacking major shifts in external climate forcing and earth system configuration. Large arrays of proxy climate records derived from natural archives have been used to reconstruct aspects of climate in recent centuries. The main approaches used have been empirical and statistical, albeit informed by prior knowledge both of the physics of the climate, and of the processes imprinting climate information in the natural archives. We propose a new direction, in which emerging tools are used to formalize the combination of process knowledge and proxy climate records to better illuminate past climate variability on these time scales of great relevance to human concerns.     

Pairwise comparisons to reconstruct mean temperature in the Arctic Atlantic Region over the last 2,000 years

Existing multi-proxy climate reconstruction methods assume the suitably transformed proxy time series are linearly related to the target climate variable, which is likely a simplifying assumption for many proxy records. Furthermore, with a single exception, these methods face problems with varying temporal resolutions of the proxy data. Here we introduce a new reconstruction method that uses the ordering of all pairs of proxy observations within each record to arrive at a consensus time series that best agrees with all proxy records. The resulting unitless composite time series is subsequently calibrated to the instrumental record to provide an estimate of past climate. By considering only pairwise comparisons, this method, which we call PaiCo, facilitates the inclusion of records with differing temporal resolutions, and relaxes the assumption of linearity to the more general assumption of a monotonically increasing relationship between each proxy series and the target climate variable. We apply PaiCo to a newly assembled collection of high-quality proxy data to reconstruct the mean temperature of the Northernmost Atlantic region, which we call Arctic Atlantic, over the last 2,000 years. The Arctic Atlantic is a dynamically important region known to feature substantial temperature variability over recent millennia, and PaiCo allows for a more thorough investigation of the Arctic Atlantic regional climate as we include a diverse array of terrestrial and marine proxies with annual to multidecadal temporal resolutions. Comparisons of the PaiCo reconstruction to recent reconstructions covering larger areas indicate greater climatic variability in the Arctic Atlantic than for the Arctic as a whole. The Arctic Atlantic reconstruction features temperatures during the Roman Warm Period and Medieval Climate Anomaly that are comparable or even warmer than those of the twentieth century, and coldest temperatures in the middle of the nineteenth century, just prior to the onset of the recent warming trend.     

200?years of European temperature variability: insights from and tests of the proxy surrogate reconstruction analog method

Spatially resolved climate reconstructions are commonly derived from long instrumental series and proxy data via linear regression based approaches that use the main modes of the climate system. Such reconstructions have been shown to underestimate climate variability and are based upon the assumption that the main modes of climate variability are stationary back in time. Climate models simulate physically consistent climate fields but cannot be taken to represent the real past climate trajectory because of their necessarily simplified scope and chaotic internal variability. Here, we present sensitivity tests of, and a 200-year temperature reconstruction from, the PSR (Proxy Surrogate Reconstruction) method. This method simultaneously capitalizes on the individual strengths of instrumental/proxy data based reconstructions and model simulations by selecting the model states (analogs) that are most similar with proxy/instrumental data available at specific places and specific moments of time. Sensitivity experiments reveal an optimal PSR configuration and indicate that 6,500 simulation years of existing climate models provide a sufficient pool of possible analogs to skillfully reconstruct monthly European temperature fields during the past 200?years. Reconstruction verification based upon only seven instrumental stations indicates potential for extensions back in time using sparse proxy data. Additionally the PSR method allows evaluation of single time series, in this case the homogeneity of instrumental series, by identifying inconsistencies with the reconstructed climate field. We present an updated European temperature reconstruction including newly homogenized instrumental records performed with the computationally efficient PSR method that proves to capture the total variance of the target.     

Late Holocene Asian summer monsoon dynamics from small but complex networks of paleoclimate data

Internal variability of the Asian monsoon system and the relationship amongst its sub-systems, the Indian and East Asian Summer Monsoon, are not sufficiently understood to predict its responses to a future warming climate. Past environmental variability is recorded in Palaeoclimate proxy data. In the Asian monsoon domain many records are available, e.g. from stalagmites, tree-rings or sediment cores. They have to be interpreted in the context of each other, but visual comparison is insufficient. Heterogeneous growth rates lead to uneven temporal sampling. Therefore, computing correlation values is difficult because standard methods require co-eval observation times, and sampling-dependent bias effects may occur. Climate networks are tools to extract system dynamics from observed time series, and to investigate Earth system dynamics in a spatio-temporal context. We establish paleoclimate networks to compare paleoclimate records within a spatially extended domain. Our approach is based on adapted linear and nonlinear association measures that are more efficient than interpolation-based measures in the presence of inter-sampling time variability. Based on this new method we investigate Asian Summer Monsoon dynamics for the late Holocene, focusing on the Medieval Warm Period (MWP), the Little Ice Age (LIA), and the recent period of warming in East Asia. We find a strong Indian Summer Monsoon (ISM) influence on the East Asian Summer Monsoon during the MWP. During the cold LIA, the ISM circulation was weaker and did not extend as far east. The most recent period of warming yields network results that could indicate a currently ongoing transition phase towards a stronger ISM penetration into China. We find that we could not have come to these conclusions using visual comparison of the data and conclude that paleoclimate networks have great potential to study the variability of climate subsystems in space and time.     

Tree-ring inferred annual mean temperature variations on the southeastern Tibetan Plateau during the last millennium and their relationships with the Atlantic Multidecadal Oscillation

We present two tree-ring chronologies for the southeastern Tibetan Plateau (TP), established by applying the signal-free regional curve standardization and standard dendrochronological methodologies to a set of ring-width series of Tibetan juniper. The relationship between tree growth and climate shows that temperature variability in the previous year is the primary factor controlling tree growth at the upper portion of the forest belt. Accordingly, we developed a mean annual temperature reconstruction covering the period A.D. 984–2009 and explaining 50 % of the instrumental variance. The spatial correlation patterns suggest that our temperature reconstruction is a reasonable proxy for temperature change over the TP. At long time scales, the temperature reconstruction shows similar warm-cold patterns to those in temperature records from other regions of the TP, indicating that decadal and multidecadal temperature variations were generally synchronous across the TP during the past millennium. The periods 1140–1350 and 1600–1800 were common warm and cold episodes over the TP, respectively. Comparison of our reconstruction with four Northern Hemisphere (NH) temperature series indicates that temperature changes on the southeastern TP have generally followed the NH temperature patterns during the past millennium. Our results also suggest that temperature variability over the TP is affected by the Atlantic Multidecadal Oscillation (AMO), with the warm (cool) phases of the AMO associated with above-average (below-average) temperatures over the TP.     

Propagating decadal sea surface temperature signal identified in modern proxy records of the tropical Pacific

Analysis of 86 years of multiple modern coral δ¹⁸O records in the tropical Pacific reveals a basin-scale decadal pattern of variability. Although coral δ¹⁸O records the effects of both temperature and seawater δ¹⁸O variability due to salinity effects, in practice, most of the records used here agree well with observations of sea surface temperature on longer timescales. These coral proxy records reveal strong variability near a 12-year period. Their relative phasing suggests a signal propagating from the southwestern subtropical Pacific to other regions. The results are consistent with recent studies based on instrumental data and with coupled climate model studies, in which advection of thermal anomalies leads to El Niño/Southern Oscillation-like variability on decadal timescales. Additionally, there is evidence for a significant shift in many of the time series, along with a decrease in the decadal variability, occurring in the early 1940s. Finally, the proxy records indicate the presence of strong teleconnections between the eastern tropical Pacific and high latitude climate.     

Temperature and Precipitation Changes in China During the Holocene

We review here proxy records of temperature and precipitation in China during the Holocene, especially the last two millennia. The quality of proxy data, methodology of reconstruction, and uncertainties in reconstruction were emphasized in comparing different temperature and precipitation reconstruction and clarifying temporal and spatial patterns of temperature and precipitation during the Holocene. The Holocene climate was generally warm and wet. The warmest period occurred in 9.6-6.2 cal ka BP, whereas a period of maximum monsoon precipitation started at about 11.0 cal ka BP and lasted until about 8.0-5.0 cal ka BP. There were a series of millennial-scale cold or dry events superimposed on the general trend of climate changes. During past two millennia, a warming trend in the 20th century was clearly detected, but the warming magnitude was smaller than the maximum level of the Medieval Warm Period and the Middle Holocene. Cold conditions occurred over the whole of China during the Little Ice Age （AD 1400-AD 1900）, but the warming of the Medieval Warm Period （AD 900-AD 1300） was not distinct in China, especially west China. The spatial pattern of precipitation showed significant regional differences in China, especially east China. The modern warm period has lasted 20 years from 1987 to 2006. Bi-decadal oscillation in precipitation variability was apparent over China during the 20th century. Solar activity and volcanic eruptions both were major forcings governing the climate variability during the last millennium.     

Marine proxy evidence linking decadal North Pacific and Atlantic climate

Decadal- to multidecadal variability in the extra-tropical North Pacific is evident in 20th century instrumental records and has significant impacts on Northern Hemisphere climate and marine ecosystems. Several studies have discussed a potential linkage between North Pacific and Atlantic climate on various time scales. On decadal time scales no relationship could be confirmed, potentially due to sparse instrumental observations before 1950. Proxy data are limited and no multi-centennial high-resolution marine geochemical proxy records are available from the subarctic North Pacific. Here we present an annually-resolved record (1818–1967) of Mg/Ca variations from a North Pacific/Bering Sea coralline alga that extends our knowledge in this region beyond available data. It shows for the first time a statistically significant link between decadal fluctuations in sea-level pressure in the North Pacific and North Atlantic. The record is a lagged proxy for decadal-scale variations of the Aleutian Low. It is significantly related to regional sea surface temperature and the North Atlantic Oscillation (NAO) index in late boreal winter on these time scales. Our data show that on decadal time scales a weaker Aleutian Low precedes a negative NAO by several years. This atmospheric link can explain the coherence of decadal North Pacific and Atlantic Multidecadal Variability, as suggested by earlier studies using climate models and limited instrumental data.     

Modelling mid-Holocene tropical climate and ENSO variability: towards constraining predictions of future change with palaeo-data

Palaeoclimate simulations provide an opportunity for climate model evaluation as well as having a potential role in assigning relative likelihood to different ensemble members in probabilistic climate change prediction, supplementing constraints provided by the instrumental record. Here we take some initial steps towards such an approach by performing ensemble experiments with the Hadley Centre HadCM3 model under pre-industrial and mid-Holocene (6,000 years before present) forcing conditions. We examine the changes in both mean tropical climate and El Niño-Southern Oscillation (ENSO) variability, as palaeoclimate records suggest that ENSO amplitude was reduced in the mid-Holocene. Experiments are performed with perturbations to physical parameters in the atmosphere–surface component of the model, and with different implementations of heat and freshwater flux adjustments. Heat flux adjustments are required to stabilise model versions in which perturbations cause a net radiative imbalance. While we find broad agreement between different model versions in terms of changes in mean climate in the mid-Holocene, a detailed and quantitative comparison with the geographically-sparse palaeo-record is limited by systematic model biases. In the simulations without seasonally-varying flux adjustments there are modest reductions in ENSO amplitude of the order of 10–15%, lower than the range of reductions inferred from coral proxy records. We examine the mechanisms for these changes, and discuss the implications for the design of future ensemble experiments to formally quantify uncertainty in climate change predictions using palaeoclimate simulations.     

The influence of regional circulation patterns on wet and dry mineral dust and sea salt deposition over Greenland

Annually resolved ice core records from different regions over the Greenland ice sheet (GrIS) are used to investigate the spatial and temporal variability of calcium (Ca²⁺, mainly from mineral dust) and sodium (Na⁺, mainly from sea salt) deposition. Cores of high common inter-annual variability are grouped with an EOF analysis, resulting in regionally representative Ca²⁺ and Na⁺ records for northeastern and central Greenland. Utilizing a regression and validation method with ERA-40 reanalysis data, these common records are associated with distinct regional atmospheric circulation patterns over the North American Arctic, Greenland, and Central to Northern Europe. These patterns are interpreted in terms of transport and deposition of the impurities. In the northeastern part of the GrIS sea salt records reflect the intrusion of marine air masses from southeasterly flow. A large fraction of the Ca²⁺ variability in this region is connected to a circulation pattern suggesting transport from the west and dry deposition. This pattern is consistent with the current understanding of a predominantly Asian source of the dust deposited over the GrIS. However, our results also indicate that a significant fraction of the inter-annual dust variability in NE and Central Greenland is determined by the frequency and intensity of wet deposition during the season of high atmospheric dust loading, rather than representing the variability of the Asian dust source and/or long-range transport to Greenland. The variances in the regional proxy records explained by the streamfunction patterns are high enough to permit reconstructions of the corresponding regional deposition regimes and the associated circulation patterns.     

Intra-seasonal rainfall and piped water revenue variability in rural Africa

Rainfall patterns influence water usage and revenue from user payments in rural Africa. We explore these dynamics by examining monthly rainfall against 4,888 records of rural piped water revenue in Ghana, Rwanda, and Uganda and quantifying revenue changes over 635 transitions between dry and wet seasons.Results show operators experience revenue variability at regional and intra-seasonal scales. Revenues fall by an average of 30 percent during the wettest months of the year in climate regimes with consistent wet season rainfall. However, seasonally stable revenues are observed in areas where consecutive dry days are common during the wet season, potentially reflecting a dependency on reliable services. We also find changes in tariff level, waterpoint connection type, and payment approach do not consistently prevent or increase seasonal revenue variability.Local revenue generation underpins delivery of drinking water services. Where rainfall patterns remain consistent, piped water operators can expect to encounter seasonal revenue reductions regardless of whether services are provided on or off premises and of how services are paid for. Revenue projections that assume consistent volumetric demand year-round may lead to shortfalls that threaten sustainability and undermine the case for future investment. Intra-seasonal rainfall analysis can enhance rural piped water revenue planning by offering localised insight into demand dynamics and revealing where climate variability may increase dependency on reliable services.     

Climatic fluctuations on the century time scale: A review of high-resolution proxy data and possible mechanisms

We review the century time scale climatic variability that is observed in high-resolution proxy data records covering the past 10 000 yr. Cyclic variations with time scales ranging from 50 to 400 yr occur in oxygen isotope ratios derived from ice cores, tree-ring index series, pollen records and sea-ice extents. Century time scale cycles can also be identified in some biological and historical records and in long-term instrumental observations. In order to appreciate the century scale cycles in the context of climatic variability in general, a brief survey of all climatic time scales is presented.The traditional interpretation that decadal-to-century scale fluctuations in the climate system are externally forced, e.g. by variations in solar properties, is questioned. A different mechanism for these fluctuations is proposed on the basis of recent findings of numerical models of the ocean's thermohaline circulation. The results indicate that this oceanic circulation exhibits natural variability on the century time scale which produces oscillations in the ocean-to-atmosphere heat flux. Although global in extent, these fluctuations are largest in the Atlantic Ocean.     

Functional clustering for Italian climate zones identification

This work presents a functional clustering procedure applied to meteorological time series. Our proposal combines time series interpolation with smoothing penalized B-spline and the partitioning around medoids clustering algorithm. Our final goal is to obtain homogeneous climate zones of Italy. We compare this approach to standard methods based on a combination of principal component analysis and Cluster Analysis (CA) and we discuss it in relation to other functional clustering approaches based on Fourier analysis and CA. We show that a functional approach is simpler than the standard methods from a methodological and interpretability point of view. Indeed, it becomes natural to find a clear connection between mathematical results and physical variability mechanisms. We discuss how the choice of the basis expansion (splines, Fourier) affects the analysis and propose some comments on their use. The basis for classification is formed by monthly values of temperature and precipitation recorded during the period 1971–2000 over 95 and 94 Italian monitoring stations, respectively. An assessment based on climatic patterns is presented to prove the consistency of the clustering and a comparison of results obtained with different methods is used to judge the functional data approach.     

Spatial gradient of temperature and potential vegetation feedback across Europe during the late Quaternary

After the last glacial maximum (LGM; 18,000 years. bp) plants migrated onto areas made free by the retreating ice sheets. Pollen records document changes in ecosystems and are an important proxy for the reconstruction of climate parameters. Here, we quantify changes in January temperature (Tjan) over the past 14,000 years using 216 European pollen records. We used functional principal component analysis to evaluate the patterns of change over this period of time. Between 14 and 12 ka, Tjan records show comparable changes across Europe. After this time period, northwest and southeast Europe experienced differing trends in Tjan, trends that became pronounced by 9 ka. The spatial gradient of Tjan records is well-correlated with the overall change from steppe to increasingly diverse forest ecosystems. In particular, differential incorporation of conifer vs. deciduous species into forest ecosystems, based on differential colonization out of refugia, likely resulted in albedo gradients over Europe that affected regional Tjan.     

Evaluation of proxy-based millennial reconstruction methods

A range of existing statistical approaches for reconstructing historical temperature variations from proxy data are compared using both climate model data and real-world paleoclimate proxy data. We also propose a new method for reconstruction that is based on a state-space time series model and Kalman filter algorithm. The state-space modelling approach and the recently developed RegEM method generally perform better than their competitors when reconstructing interannual variations in Northern Hemispheric mean surface air temperature. On the other hand, a variety of methods are seen to perform well when reconstructing surface air temperature variability on decadal time scales. An advantage of the new method is that it can incorporate additional, non-temperature, information into the reconstruction, such as the estimated response to external forcing, thereby permitting a simultaneous reconstruction and detection analysis as well as future projection. An application of these extensions is also demonstrated in the paper.     

Support for global climate reorganization during the “Medieval Climate Anomaly”

Widely distributed proxy records indicate that the Medieval Climate Anomaly (MCA; ~900–1350 AD) was characterized by coherent shifts in large-scale Northern Hemisphere atmospheric circulation patterns. Although cooler sea surface temperatures in the central and eastern equatorial Pacific can explain some aspects of medieval circulation changes, they are not sufficient to account for other notable features, including widespread aridity through the Eurasian sub-tropics, stronger winter westerlies across the North Atlantic and Western Europe, and shifts in monsoon rainfall patterns across Africa and South Asia. We present results from a full-physics coupled climate model showing that a slight warming of the tropical Indian and western Pacific Oceans relative to the other tropical ocean basins can induce a broad range of the medieval circulation and climate changes indicated by proxy data, including many of those not explained by a cooler tropical Pacific alone. Important aspects of the results resemble those from previous simulations examining the climatic response to the rapid Indian Ocean warming during the late twentieth century, and to results from climate warming simulations—especially in indicating an expansion of the Northern Hemisphere Hadley circulation. Notably, the pattern of tropical Indo-Pacific sea surface temperature (SST) change responsible for producing the proxy-model similarity in our results agrees well with MCA-LIA SST differences obtained in a recent proxy-based climate field reconstruction. Though much remains unclear, our results indicate that the MCA was characterized by an enhanced zonal Indo-Pacific SST gradient with resulting changes in Northern Hemisphere tropical and extra-tropical circulation patterns and hydroclimate regimes, linkages that may explain the coherent regional climate shifts indicated by proxy records from across the planet. The findings provide new perspectives on the nature and possible causes of the MCA—a remarkable, yet incompletely understood episode of Late Holocene climatic change.     

Temperature and Precipitation Changes in China During the Holocene

We review here proxy records of temperatare and precipitation in China during the Holocene,especially the last two millennia.The quality of proxy data,methodology of reconstruction,and uncertainties in reconstruction were emphasized in comparing different temperatare and precipitation reconstruction and clarilying temporal and spatial patterns of temperature and precipitation during the Holocene.The Holocene climate was generally warm and wet.The warmest period occurred in 9.6-6.2 cal ka BP,whereas a period of maximum monsoon precipitation started at about 11.0 cal ka BP and lasted until about 8.O-5.0 cal ka BP.There were a series of millennial-scale cold or dry events superimposed on the general trend of climate changes.During past two millennia,a warming trend in the 20th century was clearly detected,but the warming magnitude was smaller than the maximum level of the Medieval Warm Period and the Middle Holocene.Cold conditions occurred over the whole of China during the Little Ice Age (AD 1400-AD 1900),but the warming of the Medieval Warm Period(AD 900-AD 1300)was not distinct in China,especially west China.The spatial pattern of precipitation showed significant regional differences in China,especially east China.The modern warm period has lasted 20、years from 1987 to 2006.Bi-decadal oscillation in precipitation variability was apparent over China during the 20th century. Solar activity and volcanic eruptions both were major forcings governing the climate variability during the last millennium.     

Identification of linear relationships from noisy data using errors-in-variables models??relevance for reconstruction of past climate from tree-ring and other proxy information

Reliable paleoclimate reconstructions are needed to assess if the recent climatic changes are unusual compared to pre-industrial climate variability. Here, we focus on one important problem in climate reconstructions: Transfer functions relating proxies (predictors) and target climatic quantities (predictands) can be seriously biased when predictand and predictor noise is not adequately accounted for, resulting in biased amplitudes of reconstructed climatic time series. We argue for errors-in-variables models (EVM) for unbiased identification of linear structural relationships between noisy proxies and target climatic quantities by (1) introducing underlying statistical concepts and (2) demonstrating the potential biases of using the EVM approach, the most commonly used direct ordinary least squares (OLS) regression, inverse OLS regression, or the reduced major axis method (??variance matching??) with a simulation example of artificial noise-disturbed sinusoidal time series. We then develop an alternative strategy for paleoclimate reconstruction from tree-ring and other proxy data, explicitly accounting for the identified problem.     

Reconstructed cool- and warm-season precipitation over the tribal lands of northeastern Arizona

For over a decade, the Hopi Tribe and Navajo Nation of northeastern Arizona have suffered the effects of persistent drought conditions. Severe dry spells have critically impacted natural ecosystems, water resources, and regional livelihoods including dryland farming and ranching. Drought planning and resource management efforts in the region are based largely on the instrumental climate record, which contains a limited number of severe, sustained droughts. In this study, a new network of moisture-sensitive tree-ring chronologies provides the basis for evaluating the longer-term temporal variability of precipitation in the Four Corners region. By analyzing the earlywood and latewood components within each annual tree ring, we are able to generate separate, centuries-long reconstructions of both cool- (October-April) and warm-season (July-August) precipitation. These proxy records offer new insights into seasonal drought characteristics and indicate that the instrumental record fails to adequately represent precipitation variability over the past 400 years. Through the use of two different analysis techniques, we identify multiyear and decadal-scale drought events more severe than any in the modern era. Furthermore, the reconstructions suggest that many of the historically significant droughts of the past (e.g., 17th century Puebloan drought) were not merely winter phenomena, but persisted through the summer season as well. By comparing these proxy records with historical documents, we are able to independently validate the reconstructions and better understand the socioeconomic and environmental significance of past climate anomalies on the tribal lands of northeastern Arizona.     

Climate change at the ecosystem scale: a 50-year record in New Hampshire

Observing the full range of climate change impacts at the local scale is difficult. Predicted rates of change are often small relative to interannual variability, and few locations have sufficiently comprehensive long-term records of environmental variables to enable researchers to observe the fine-scale patterns that may be important to understanding the influence of climate change on biological systems at the taxon, community, and ecosystem levels. We examined a 50-year meteorological and hydrological record from the Hubbard Brook Experimental Forest (HBEF) in New Hampshire, an intensively monitored Long-Term Ecological Research site. Of the examined climate metrics, trends in temperature were the most significant (ranging from 0.7 to 1.3 °C increase over 40–50 year records at 4 temperature stations), while analysis of precipitation and hydrologic data yielded mixed results. Regional records show generally similar trends over the same time period, though longer-term (70–102 year) trends are less dramatic. Taken together, the results from HBEF and the regional records indicate that the climate has warmed detectably over 50 years, with important consequences for hydrological processes. Understanding effects on ecosystems will require a diversity of metrics and concurrent ecological observations at a range of sites, as well as a recognition that ecosystems have existed in a directionally changing climate for decades, and are not necessarily in equilibrium with the current climate.