Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

A first chironomid-based summer temperature reconstruction (13–5 ka BP) around 49°N in inland Europe compared with local lake development

Temperature reconstructions for the end of the Pleistocene and the first half of the Holocene based on biotic proxies are rare for inland Europe around 49°N. We analysed a 7 m long sequence of lake deposits in the Vihorlat Mts in eastern Slovakia (820 m a.s.l.). Chironomid head capsules were used to reconstruct mean July temperature (TJuly), other proxies (diatoms, green algae, pollen, geochemistry) were used to reconstruct local environmental changes that might have affected the climate reconstruction, such as epilimnetic total phosphorus concentrations (TP), lake level changes and development of surrounding vegetation. During the Younger Dryas (YD), temperature fluctuated between 7 and 11 °C, with distinct, decadal to centennial scale variations, that agree with other palaeoclimate records in Europe such as δ18O content in stalagmites or Greenland ice cores. The results indicate that the site was somewhat colder than expected from the general south-to-north YD temperature gradient within Europe, possibly because of north-facing exposition. The warmer phases of the YD were characterised by low water level or even complete desiccation of the lake (12,200-12,400 cal yr BP). At the Late-Glacial/Holocene transition TJuly steeply increased from from 11 to 15.5 °C (11,700-11,400 cal yr BP) - the highest TJuly for entire sequence. This rapid climate change was reflected by all proxies as a compositional change and increasing species diversity. The open woodlands of Pinus, Betula, Larix and Picea were replaced by broad-leaved temperate forests dominated by Betula, later by Ulmus and finally by Corylus (ca 9700 cal yr BP). At the same time, input of eroded coarse-grained material into the lake decreased and organic matter (LOI) and biogenic silica increased. The Early-Holocene climate was rather stable till 8700 cal yr BP, with temporary decrease in TJuly around 11,200 cal yr BP. The lake was productive with a well-developed littoral, as indicated by both diatoms and chironomids. A distinct decline of TJuly to 10 °C between 8700 and 8000 cal yr BP was associated with decreasing chironomid diversity and increasing climate moistening indicated by pollen. Tychoplanktonic and phosphorus-demanding diatoms increased which might be explained by hydrological and land-cover changes. Later, a gradual warming started after 7000 cal yr BP and representation of macrophytes, periphytic diatoms and littoral chironomids increased. Our results suggest that the Holocene thermal maximum was taking place unusually early in the Holocene at our study site, but its timing might be affected by topography and mesoclimate. We further demonstrated that temperature changes had coincided with variations in local hydrology

A chrysophyte stomatocyst-based reconstruction of cold-season air temperature from Alpine Lake Silvaplana (AD 1500–2003); methods and concepts for quantitative inferences

Relatively little is known about past cold-season temperature variability in high-Alpine regions because of a lack of natural cold-season temperature proxies as well as under-representation of high-altitude sites in meteorological, early-instrumental and documentary data sources. Recent studies have shown that chrysophyte stomatocysts, or simply cysts (sub-fossil algal remains of Chrysophyceae and Synurophyceae), are among the very few natural proxies that can be used to reconstruct cold-season temperatures. This study presents a quantitative, high-resolution (5-year), cold-season (Oct–May) temperature reconstruction based on sub-fossil chrysophyte stomatocysts in the annually laminated (varved) sediments of high-Alpine Lake Silvaplana, SE Switzerland (1,789 m a.s.l.), since AD 1500. We first explore the method used to translate an ecologically meaningful variable based on a biological proxy into a simple climate variable. A transfer function was applied to reconstruct the ‘date of spring mixing’ from cyst assemblages. Next, statistical regression models were tested to convert the reconstructed ‘dates of spring mixing’ into cold-season surface air temperatures with associated errors. The strengths and weaknesses of this approach are thoroughly tested. One much-debated, basic assumption for reconstructions (‘stationarity’), which states that only the environmental variable of interest has influenced cyst assemblages and the influence of confounding variables is negligible over time, is addressed in detail. Our inferences show that past cold-season air-temperature fluctuations were substantial and larger than those of other temperature reconstructions for Europe and the Alpine region. Interestingly, in this study, recent cold-season temperatures only just exceed those of previous, multi-decadal warm phases since AD 1500. These findings highlight the importance of local studies to assess natural climate variability at high altitudes.

The chironomid-temperature relationship: expression in nature and palaeoenvironmental implications

Fossils of chironomid larvae (non-biting midges) preserved in lake sediments are well-established palaeotemperature indicators which, with the aid of numerical chironomid-based inference models (transfer functions), can provide quantitative estimates of past temperature change. This approach to temperature reconstruction relies on the strong relationship between air and lake surface water temperature and the distribution of individual chironomid taxa (species, species groups, genera) that has been observed in different climate regions (arctic, subarctic, temperate and tropical) in both the Northern and Southern hemisphere. A major complicating factor for the use of chironomids for palaeoclimate reconstruction which increases the uncertainty associated with chironomid-based temperature estimates is that the exact nature of the mechanism responsible for the strong relationship between temperature and chironomid assemblages in lakes remains uncertain. While a number of authors have provided state of the art overviews of fossil chironomid palaeoecology and the use of chironomids for temperature reconstruction, few have focused on examining the ecological basis for this approach. Here, we review the nature of the relationship between chironomids and temperature based on the available ecological evidence. After discussing many of the surveys describing the distribution of chironomid taxa in lake surface sediments in relation to temperature, we also examine evidence from laboratory and field studies exploring the effects of temperature on chironomid physiology, life cycles and behaviour. We show that, even though a direct influence of water temperature on chironomid development, growth and survival is well described, chironomid palaeoclimatology is presently faced with the paradoxical situation that the relationship between chironomid distribution and temperature seems strongest in relatively deep, thermally stratified lakes in temperate and subarctic regions in which the benthic chironomid fauna lives largely decoupled from the direct influence of air and surface water temperature. This finding suggests that indirect effects of temperature on physical and chemical characteristics of lakes play an important role in determining the distribution of lake-living chironomid larvae. However, we also demonstrate that no single indirect mechanism has been identified that can explain the strong relationship between chironomid distribution and temperature in all regions and datasets presently available. This observation contrasts with the previously published hypothesis that climatic effects on lake nutrient status and productivity may be largely responsible for the apparent correlation between chironomid assemblage distribution and temperature. We conclude our review by summarizing the implications of our findings for chironomid-based palaeoclimatology and by pointing towards further avenues of research necessary to improve our mechanistic understanding of the chironomid-temperature relationship.

Late Holocene air temperature variability reconstructed from the sediments of Laguna Escondida, Patagonia, Chile (45 degrees 30 ' S)

Climate and environmental reconstructions from natural archives are important for the interpretation of current climatic change. Few quantitative high-resolution reconstructions exist for South America which is the only land mass extending from the tropics to the southern high latitudes at 56°S. We analyzed sediment cores from two adjacent lakes in Northern Chilean Patagonia, Lago Castor (45°36′S, 71°47′W) and Laguna Escondida (45°31′S, 71°49′W). Radiometric dating (210Pb, 137Cs, 14C-AMS) suggests that the cores reach back to c. 900 BC (Laguna Escondida) and c. 1900 BC (Lago Castor). Both lakes show similarities and reproducibility in sedimentation rate changes and tephra layer deposition. We found eight macroscopic tephras (0.2–5.5 cm thick) dated at 1950 BC, 1700 BC, at 300 BC, 50 BC, 90 AD, 160 AD, 400 AD and at 900 AD. These can be used as regional time-synchronous stratigraphic markers. The two thickest tephras represent known well-dated explosive eruptions of Hudson volcano around 1950 and 300 BC. Biogenic silica flux revealed in both lakes a climate signal and correlation with annual temperature reanalysis data (calibration 1900–2006 AD; Lago Castor r = 0.37; Laguna Escondida r = 0.42, seven years filtered data). We used a linear inverse regression plus scaling model for calibration and leave-one-out cross-validation (RMSEv = 0.56 °C) to reconstruct sub decadal-scale temperature variability for Laguna Escondida back to AD 400. The lower part of the core from Laguna Escondida prior to AD 400 and the core of Lago Castor are strongly influenced by primary and secondary tephras and, therefore, not used for the temperature reconstruction. The temperature reconstruction from Laguna Escondida shows cold conditions in the 5th century (relative to the 20th century mean), warmer temperatures from AD 600 to AD 1150 and colder temperatures from AD 1200 to AD 1450. From AD 1450 to AD 1700 our reconstruction shows a period with stronger variability and on average higher values than the 20th century mean. Until AD 1900 the temperature values decrease but stay slightly above the 20th century mean. Most of the centennial-scale features are reproduced in the few other natural climate archives in the region. The early onset of cool conditions from c. AD 1200 onward seems to be confirmed for this region.

Mineralogy-based quantitative precipitation and temperature reconstructions from annually laminated lake sediments (Swiss Alps) since AD 1580

[1] We present quantitative autumn, summer and annual precipitation and summer temperature reconstructions from proglacial annually laminated Lake Silvaplana, eastern Swiss Alps back to AD 1580. We used X-ray diffraction peak intensity ratios of minerals in the sediment layers (quartz qz, plagioclase pl, amphibole am, mica mi) that are diagnostic for different source areas and hydro-meteorological transport processes in the catchment. XRD data were calibrated with meteorological data (AD 1800/1864–1950) and revealed significant correlations: mi/pl with SON precipitation (r = 0.56, p < 0.05) and MJJAS precipitation (r = 0.66, p < 0.01); qz/mi with MJJAS temperature (r = −0.72, p < 0.01)and qz/am with annual precipitation (r = −0.54, p < 0.05). Geological catchment settings and hydro-meteorological processes provide deterministic explanations for the correlations. Our summer temperature reconstruction reproduces the typical features of past climate variability known from independent data sets. The precipitation reconstructions show a LIA climate moister than today. Exceptionally wet periods in our reconstruction coincide with regional glacier advances.

Inter-hemispheric temperature variability over the past millennium

The Earth’s climate system is driven by a complex interplay of internal chaotic dynamics and natural and anthropogenic external forcing. Recent instrumental data have shown a remarkable degree of asynchronicity between Northern Hemisphere and Southern Hemisphere temperature fluctuations, thereby questioning the relative importance of internal versus external drivers of past as well as future climate variability1, 2, 3. However, large-scale temperature reconstructions for the past millennium have focused on the Northern Hemisphere4, 5, limiting empirical assessments of inter-hemispheric variability on multi-decadal to centennial timescales. Here, we introduce a new millennial ensemble reconstruction of annually resolved temperature variations for the Southern Hemisphere based on an unprecedented network of terrestrial and oceanic palaeoclimate proxy records. In conjunction with an independent Northern Hemisphere temperature reconstruction ensemble5, this record reveals an extended cold period (1594–1677) in both hemispheres but no globally coherent warm phase during the pre-industrial (1000–1850) era. The current (post-1974) warm phase is the only period of the past millennium where both hemispheres are likely to have experienced contemporaneous warm extremes. Our analysis of inter-hemispheric temperature variability in an ensemble of climate model simulations for the past millennium suggests that models tend to overemphasize Northern Hemisphere–Southern Hemisphere synchronicity by underestimating the role of internal ocean–atmosphere dynamics, particularly in the ocean-dominated Southern Hemisphere. Our results imply that climate system predictability on decadal to century timescales may be lower than expected based on assessments of external climate forcing and Northern Hemisphere temperature variations5, 6 alone.

A 600 years warm-season temperature record from varved sediments of Lago Plomo, Northern Patagonia, Chile (47°S)

High-resolution records of calibrated proxy data for the past millennium are fundamental to place current changes into the context of pre-industrial natural forced and unforced variability. Although the need for regional spatially-explicit comprehensive reconstructions is widely recognized, the proxy data sources are still scarce, particularly for the Southern Hemisphere and especially for South America. We present a 600-year long warm season temperature record from varved sediments of Lago Plomo, a proglacial lake of the Northern Patagonian Ice field in Southern Chile (46°59′S, 72°52′W, 203 m a.s.l.). The thickness of the bright summer sediment layer relative to the dark winter layer (measured as total brightness; % reflectance 400–730 nm) is calibrated against warm season SONDJF temperature (1900–2009; r = 0.58, p(aut) = 0.056, RE = 0.52; CE = 0.15, RMSEP = 0.28 °C; five-year triangular filtered data). In Lago Plomo, warm summer temperatures lead to enhanced glacier melt and suspended sediment transport, which results in thicker light summer layers and to brighter sediments. Although Patagonia shows pronounced regional differences in decadal temperature trends and variability, the 600 years temperature reconstruction from Lago Plomo compares favourably with other regional/continental temperature records, but also emphasizes significant regional differences for which no data and information existed so far. These regional differences seem to be real as they are also reflected in modern climate data sets (1900–2010). The reconstruction shows pronounced subdecadal – multidecadal variability with cold phases during parts of the Little Ice Age (16th and 18th centuries) and in the beginning of the 20th century. The most prominent warm phase is the 19th century which is as warm as the second half of the 20th century. The exceptional summer warmth AD 1780–1810 is also found in other archives of Northern Patagonia and Central Chile. Our record shows the delayed 20th century warming in the Southern Hemisphere. The comparison between winter precipitation and summer temperature (inter-seasonal coupling) from Lago Plomo reveals alternating phases with parallel and contrasting decadal trends of winter precipitation and summer temperature (positive and negative running correlations Rwinter PP; summer TT). This observation from the sediment proxy data is also confirmed by two sets of reanalysis data for the 20th century. Reanalysis data show that phases with negative correlations between winter precipitation and summer temperature (e.g., dry winters and warm summers) at Lago Plomo are characteristic for periods when circumpolar Westerly flow is displaced southward and enhanced around 60°S.

A multi-proxy approach for revealing recent climatic changes in the Russian Altai

For the first time we present a multi-proxy data set for the Russian Altai, consisting of Siberian larch tree-ring width (TRW), latewood density (MXD), δ13C and δ18O in cellulose chronologies obtained for the period 1779–2007 and cell wall thickness (CWT) for 1900–2008. All of these parameters agree well between each other in the high-frequency variability, while the low-frequency climate information shows systematic differences. The correlation analysis with temperature and precipitation data from the closest weather station and gridded data revealed that annual TRW, MXD, CWT, and δ13C data contain a strong summer temperature signal, while δ18O in cellulose represents a mixed summer and winter temperature and precipitation signal. The temperature and precipitation reconstructions from the Belukha ice core and Teletskoe lake sediments were used to investigate the correspondence of different independent proxies. Low frequency patterns in TRW and δ13C chronologies are consistent with temperature reconstructions from nearby Belukha ice core and Teletskoe lake sediments showing a pronounced warming trend in the last century. Their combination could be used for the regional temperature reconstruction. The long-term δ18O trend agrees with the precipitation reconstruction from the Teletskoe lake sediment indicating more humid conditions during the twentieth century. Therefore, these two proxies could be combined for the precipitation reconstruction.

Impacts of changing climate and land use on vegetation dynamics in a Mediterranean ecosystem: insights from paleoecology and dynamic modeling

Forests near the Mediterranean coast have been shaped by millennia of human disturbance. Consequently, ecological studies relying on modern observations or historical records may have difficulty assessing natural vegetation dynamics under current and future climate. We combined a sedimentary pollen record from Lago di Massacciucoli, Tuscany, Italy with simulations from the LandClim dynamic vegetation model to determine what vegetation preceded intense human disturbance, how past changes in vegetation relate to fire and browsing, and the potential of an extinct vegetation type under present climate. We simulated vegetation dynamics near Lago di Massaciucoli for the last 7,000 years using a local chironomid-inferred temperature reconstruction with combinations of three fire regimes (small infrequent, large infrequent, small frequent) and three browsing intensities (no browsing, light browsing, and moderate browsing), and compared model output to pollen data. Simulations with low disturbance support pollen-inferred evidence for a mixed forest dominated by Quercus ilex (a Mediterranean species) and Abies alba (a montane species). Whereas pollen data record the collapse of A. alba after 6000 cal yr bp, simulated populations expanded with declining summer temperatures during the late Holocene. Simulations with increased fire and browsing are consistent with evidence for expansion by deciduous species after A. alba collapsed. According to our combined paleo-environmental and modeling evidence, mixed Q. ilex and A. alba forests remain possible with current climate and limited disturbance, and provide a viable management objective for ecosystems near the Mediterranean coast and in regions that are expected to experience a mediterranean-type climate in the future.

Climatic and environmental changes during the Weichselian Lateglacial Interstadial in the Weerterbos region, the Netherlands

Lake sediment records from the Weerterbos region, in the southern Netherlands, were studied to reconstruct summer temperature and environmental changes during the Weichselian Lateglacial Interstadial. A sediment core obtained from a small lacustrine basin was analysed for multiple proxies, including lithological changes, oxygen isotopes of bulk carbonates, pollen and chironomids. It was found that the oxygen isotope record differed strongly from the other proxies. Based on a comparison with three additional lake sediment records from the same region, it emerged that the oxygen isotope records were strongly affected by local environmental conditions, impeding the distinction of a regional palaeoclimate signal. The chironomid-inferred July air temperature reconstruction produced inferred interstadial temperatures ranging between ∼15° and 18°C, largely consistent with previously published results from the northern part of the Netherlands. A temporary regressive phase in the pollen record, which can be tentatively correlated with the Older Dryas, preceded the expansion of birch woodland. Despite differences between the four pollen records from the Weerterbos region, a comparable regressive vegetation phase that was possibly the result of a shift to drier conditions could be discerned in all of the profiles. In addition, a temporary temperature decline of ∼1.5°C was inferred from the chironomid record during this regressive phase. The multi-proxy approach used here enabled a direct comparison of inferred changes in temperature, vegetation and environmental conditions at an individual site, while the multi-site approach provided insight into the factors influencing the pollen and isotope records from these small-scale depressions.

Late Holocene summer temperatures in the central Andes reconstructed from the sediments of high-elevatin Laguna Chepical, Chile (32° S)

High-resolution reconstructions of climate variability that cover the past millennia are necessary to improve the understanding of natural and anthropogenic climate change across the globe. Although numerous records are available for the mid- and high-latitudes of the Northern Hemisphere, global assessments are still compromised by the scarcity of data from the Southern Hemisphere. This is particularly the case for the tropical and subtropical areas. In addition, high elevation sites in the South American Andes may provide insight into the vertical structure of climate change in the mid-troposphere. This study presents a 3000 yr-long austral summer (November to February) temperature reconstruction derived from the 210Pb- and 14C-dated organic sediments of Laguna Chepical (32°16' S, 70°30' W, 3050 m a.s.l.), a high-elevation glacial lake in the subtropical Andes of central Chile. Scanning reflectance spectroscopy in the visible light range provided the spectral index R570/R630, which reflects the clay mineral content in lake sediments. For the calibration period (AD 1901–2006), the R570/R630 data were regressed against monthly meteorological reanalysis data, showing that this proxy was strongly and significantly correlated with mean summer (NDJF) temperatures (R3 yr = −0.63, padj = 0.01). This calibration model was used to make a quantitative temperature reconstruction back to 1000 BC. The reconstruction (with a model error RMSEPboot of 0.33 °C) shows that the warmest decades of the past 3000 yr occurred during the calibration period. The 19th century (end of the Little Ice Age (LIA)) was cool. The prominent warmth reconstructed for the 18th century, which was also observed in other records from this area, seems systematic for subtropical and southern South America but remains difficult to explain. Except for this warm period, the LIA was generally characterized by cool summers. Back to AD 1400, the results from this study compare remarkably well to low altitude records from the Chilean Central Valley and southern South America. However, the reconstruction from Laguna Chepical does not show a warm Medieval Climate Anomaly during the 12–13th century, which is consistent with records from tropical South America. The Chepical record also indicates substantial cooling prior to 800 BC. This coincides with well-known regional as well as global glacier advances which have been attributed to a grand solar minimum. This study thus provides insight into the climatic drivers and temperature patterns in a region for which currently very few data are available. It also shows that since ca. AD 1400, long-term temperature patterns were generally similar at low and high altitudes in central Chile.

A model-data comparison of Holocene timberline changes in the Swiss Alps reveals past and future drivers of mountain forest dynamics

Mountain vegetation is strongly affected by temperature and is expected to shift upwards with climate change. Dynamic vegetation models are often used to assess the impact of climate on vegetation and model output can be compared with paleobotanical data as a reality check. Recent paleoecological studies have revealed regional variation in the upward shift of timberlines in the Northern and Central European Alps in response to rapid warming at the Younger Dryas/Preboreal transition ca. 11700years ago, probably caused by a climatic gradient across the Alps. This contrasts with previous studies that successfully simulated the early Holocene afforestation in the (warmer) Central Alps with a chironomid-inferred temperature reconstruction from the (colder) Northern Alps. We use LandClim, a dynamic landscape vegetation model to simulate mountain forests under different temperature, soil and precipitation scenarios around Iffigsee (2065m a.s.l.) a lake in the Northwestern Swiss Alps, and compare the model output with the paleobotanical records. The model clearly overestimates the upward shift of timberline in a climate scenario that applies chironomid-inferred July-temperature anomalies to all months. However, forest establishment at 9800 cal. BP at Iffigsee is successfully simulated with lower moisture availability and monthly temperatures corrected for stronger seasonality during the early Holocene. The model-data comparison reveals a contraction in the realized niche of Abies alba due to the prominent role of anthropogenic disturbance after ca. 5000 cal. BP, which has important implications for species distribution models (SDMs) that rely on equilibrium with climate and niche stability. Under future climate projections, LandClim indicates a rapid upward shift of mountain vegetation belts by ca. 500m and treeline positions of ca. 2500m a.s.l. by the end of this century. Resulting biodiversity losses in the alpine vegetation belt might be mitigated with low-impact pastoralism to preserve species-rich alpine meadows.

Stacking of discontinuous regional palaeoclimate records: Chironomid-based summer temperatures from the Alpine region

Since multi-site reconstructions are less affected by site-specific climatic effects and artefacts, regional palaeotemperature reconstructions based on a number of sites can provide more robust estimates of centennial- to millennial-scale temperature trends than individual, site-specific records. Furthermore, reconstructions based on multiple records are necessary for developing continuous climate records over time scales longer than covered by individual sequences. Here, we present a procedure for developing such reconstructions based on relatively short (centuries to millennia), discontinuously sampled records as are typically developed when using biotic proxies in lake sediments for temperature reconstruction. The approach includes an altitudinal correction of temperatures, an interpolation of individual records to equal time intervals, a stacking procedure for sections of the interval of interest that have the same records available, as well as a splicing procedure to link the individual stacked records into a continuous reconstruction. Variations in the final, stacked and spliced reconstruction are driven by variations in the individual records, whereas the absolute temperature values are determined by the stacked segment based on the largest number of records. With numerical simulations based on the NGRIP δ18O record, we demonstrate that the interpolation and stacking procedure provides an approximation of a smoothed palaeoclimate record if based on a sufficient number of discontinuously sampled records. Finally, we provide an example of a stacked and spliced palaeotemperature reconstruction 15000–90 calibrated 14C yr BP based on six chironomid records from the northern and central Swiss Alps and eastern France to discuss the potential and limitations of this approach.