A paleoclimate record with tephrochronological age control for the last glacial-interglacial cycle from Lake Ohrid, Albania and Macedonia

Abstract Lake Ohrid is probably of Pliocene age, and the oldest extant lake in Europe. In this study climatic and environmental changes during the last glacial-interglacial cycle are reconstructed using lithological, sedimentological, geochemical and physical proxy analysis of a 15-m-long sediment succession from Lake Ohrid. A chronological framework is derived from tephrochronology and radiocarbon dating, which yields a basal age of ca. 136 ka. The succession is not continuous, however, with a hiatus between ca. 97.6 and 81.7 ka. Sediment accumulation in course of the last climatic cycle is controlled by the complex interaction of a variety of climate-controlled parameters and their impact on catchment dynamics, limnology, and hydrology of the lake. Warm interglacial and cold glacial climate conditions can be clearly distinguished from organic matter, calcite, clastic detritus and lithostratigraphic data. During interglacial periods, short-term fluctuations are recorded by abrupt variations in organic matter and calcite content, indicating climatically-induced changes in lake productivity and hydrology. During glacial periods, high variability in the contents of coarse silt to fine sand sized clastic matter is probably a function of climatically-induced changes in catchment dynamics and wind activity. In some instances tephra layers provide potential stratigraphic markers for short-lived climate perturbations. Given their widespread distribution in sites across the region, tephra analysis has the potential to provide insight into variation in the impact of climate and environmental change across the Mediterranean.

Greenland accumulation and its connection to the large-scale atmospheric circulation in ERA-Interim and paleoclimate simulations

Changes in Greenland accumulation and the stability in the relationship between accumulation variability and large-scale circulation are assessed by performing time-slice simulations for the present day, the preindustrial era, the early Holocene, and the Last Glacial Maximum (LGM) with a comprehensive climate model. The stability issue is an important prerequisite for reconstructions of Northern Hemisphere atmospheric circulation variability based on accumulation or precipitation proxy records from Greenland ice cores. The analysis reveals that the relationship between accumulation variability and large-scale circulation undergoes a significant seasonal cycle. As the contributions of the individual seasons to the annual signal change, annual mean accumulation variability is not necessarily related to the same atmospheric circulation patterns during the different climate states. Interestingly, within a season, local Greenland accumulation variability is indeed linked to a consistent circulation pattern, which is observed for all studied climate periods, even for the LGM. Hence, it would be possible to deduce a reliable reconstruction of seasonal atmospheric variability (e.g., for North Atlantic winters) if an accumulation or precipitation proxy were available that resolves single seasons. We further show that the simulated impacts of orbital forcing and changes in the ice sheet topography on Greenland accumulation exhibit strong spatial differences, emphasizing that accumulation records from different ice core sites regarding both interannual and long-term (centennial to millennial) variability cannot be expected to look alike since they include a distinct local signature. The only uniform signal to external forcing is the strong decrease in Greenland accumulation during glacial (LGM) conditions and an increase associated with the recent rise in greenhouse gas concentrations.

Paleoclimate data-model comparison and the role of climate forcings over the past 1500 years

The past 1500 years provide a valuable opportunity to study the response of the climate system to external forcings. However, the integration of paleoclimate proxies with climate modeling is critical to improving the understanding of climate dynamics. In this paper, a climate system model and proxy records are therefore used to study the role of natural and anthropogenic forcings in driving the global climate. The inverse and forward approaches to paleoclimate data–model comparison are applied, and sources of uncertainty are identified and discussed. In the first of two case studies, the climate model simulations are compared with multiproxy temperature reconstructions. Robust solar and volcanic signals are detected in Southern Hemisphere temperatures, with a possible volcanic signal detected in the Northern Hemisphere. The anthropogenic signal dominates during the industrial period. It is also found that seasonal and geographical biases may cause multiproxy reconstructions to overestimate the magnitude of the long-term preindustrial cooling trend. In the second case study, the model simulations are compared with a coral δ18O record from the central Pacific Ocean. It is found that greenhouse gases, solar irradiance, and volcanic eruptions all influence the mean state of the central Pacific, but there is no evidence that natural or anthropogenic forcings have any systematic impact on El Niño–Southern Oscillation. The proxy climate relationship is found to change over time, challenging the assumption of stationarity that underlies the interpretation of paleoclimate proxies. These case studies demonstrate the value of paleoclimate data–model comparison but also highlight the limitations of current techniques and demonstrate the need to develop alternative approaches.

Oxygen isotope ratios (18O/16O) of hemicellulose-derived sugar biomarkers in plants, soils and sediments as paleoclimate proxy II: insight from a climate transect study

The oxygen isotopic composition of precipitation (δ18Oprec) is well known to be a valuable (paleo-)climate proxy. Paleosols and sediments and hemicelluloses therein have the potential to serve as archives recording the isotopic composition of paleoprecipitation. In a companion paper (Zech et al., 2014) we investigated δ18Ohemicellulose values of plants grown under different climatic conditions in a climate chamber experiment. Here we present results of compound-specific δ18O analyses of arabinose, fucose and xylose extracted from modern topsoils (n = 56) along a large humid-arid climate transect in Argentina in order to answer the question whether hemicellulose biomarkers in soils reflect δ18Oprec. The results from the field replications indicate that the homogeneity of topsoils with regard to δ18Ohemicellulose is very high for most of the 20 sampling sites. Standard deviations for the field replications are 1.5‰, 2.2‰ and 1.7‰, for arabinose, fucose and xylose, respectively. Furthermore, all three hemicellulose biomarkers reveal systematic and similar trends along the climate gradient. However, the δ18Ohemicellulose values (mean of the three sugars) do not correlate positively with δ18Oprec (r = −0.54, p < 0.014, n = 20). By using a Péclet-modified Craig-Gordon (PMCG) model it can be shown that the δ18Ohemicellulose values correlate highly significantly with modeled δ18Oleaf water values (r = 0.81, p < 0.001, n = 20). This finding suggests that hemicellulose biomarkers in (paleo-)soils do not simply reflect δ18Oprec but rather δ18Oprec altered by evaporative 18O enrichment of leaf water due to evapotranspiration. According to the modeling results, evaporative 18O enrichment of leaf water is relatively low (∼10‰) in the humid northern part of the Argentinian transect and much higher (up to 19‰) in the arid middle and southern part of the transect. Model sensitivity tests corroborate that changes in relative air humidity exert a dominant control on evaporative 18O enrichment of leaf water and thus δ18Ohemicellulose, whereas the effect of temperature changes is of minor importance. While oxygen exchange and degradation effects seem to be negligible, further factors needing consideration when interpreting δ18Ohemicellulose values obtained from (paleo-)soils are evaporative 18O enrichment of soil water, seasonality effects, wind effects and in case of abundant stem/root-derived organic matter input a partial loss of the evaporative 18O enrichment of leaf water. Overall, our results prove that compound-specific δ18O analyses of hemicellulose biomarkers in soils and sediments are a promising tool for paleoclimate research. However, disentangling the two major factors influencing δ18Ohemicellulose, namely δ18Oprec and relative air humidity controlled evaporative 18O enrichment of leaf water, is challenging based on δ18O analyses alone.

Patterns of millennial variability over the last 500 ka

Millennial variability is a robust feature of many paleoclimate records, at least throughout the last several glacial cycles. Here we use the mean signal from Antarctic climate events 1 to 4 to probe the EPICA Dome C temperature proxy reconstruction through the last 500 ka for similar millennial-scale events. We find that clusters of millennial events occurred in a regular fashion over half of the time during this with a mean recurrence interval of 21 kyr. We find that there is no consistent link between ice-rafted debris deposition and millennial variability. Instead we speculate that changes in the zonality of atmospheric circulation over the North Atlantic form a viable alternative to freshwater release from icebergs as a trigger for millennial variability. We suggest that millennial changes in the zonality of atmospheric circulation over the North Atlantic are linked to precession via sea-ice feedbacks and that this relationship is modified by the presence of the large, Northern Hemisphere ice sheets during glacial periods.

Fingerprints of changes in the terrestrial carbon cycle in response to large reorganizations in ocean circulation

CO2 and carbon cycle changes in the land, ocean and atmosphere are investigated using the comprehensive carbon cycle-climate model NCAR CSM1.4-carbon. Ensemble simulations are forced with freshwater perturbations applied at the North Atlantic and Southern Ocean deep water formation sites under pre-industrial climate conditions. As a result, the Atlantic Meridional Overturning Circulation reduces in each experiment to varying degrees. The physical climate fields show changes qualitatively in agreement with results documented in the literature, but there is a clear distinction between northern and southern perturbations. Changes in the physical variables, in turn, affect the land and ocean biogeochemical cycles and cause a reduction, or an increase, in the atmospheric CO2 concentration by up to 20 ppmv, depending on the location of the perturbation. In the case of a North Atlantic perturbation, the land biosphere reacts with a strong reduction in carbon stocks in some tropical locations and in high northern latitudes. In contrast, land carbon stocks tend to increase in response to a southern perturbation. The ocean is generally a sink of carbon although large reorganizations occur throughout various basins. The response of the land biosphere is strongest in the tropical regions due to a shift of the Intertropical Convergence Zone. The carbon fingerprints of this shift, either to the south or to the north depending on where the freshwater is applied, can be found most clearly in South America. For this reason, a compilation of various paleoclimate proxy records of Younger Dryas precipitation changes are compared with our model results. The proxy records, in general, show good agreement with the model's response to a North Atlantic freshwater perturbation.

Climate forcing reconstructions for use in the PMIP simulations of the last millennium

Simulations of climate over the Last Millennium (850–1850 CE) have been incorporated into the third phase of the Paleoclimate Modelling Intercomparison Project (PMIP3). The drivers of climate over this period are chiefly orbital, solar, volcanic, changes in land use/land cover and some variation in greenhouse gas levels. While some of these effects can be easily defined, the reconstructions of solar, volcanic and land use-related forcing are more uncertain. We describe here the approach taken in defining the scenarios used in PMIP3, document the forcing reconstructions and discuss likely implications.

Surface uplift and climate change - the geomorphic evolution of the western escarpment of the Andes of northern Chile between the Miocene and present

The Western Escarpment of the Andes at 18.30°S (Arica area, northern Chile) is a classical example for a transient state in landscape evolution. This part of the Andes is characterized by the presence of >10,000 km2 plains that formed between the Miocene and the present, and >1500 m deeply incised valleys. Although processes in these valleys scale the rates of landscape evolution, determinations of ages of incision, and more importantly, interpretations of possible controls on valley formation have been controversial. This paper uses morphometric data and observations, stratigraphic information, and estimates of sediment yields for the time interval between ca. 7.5 Ma and present to illustrate that the formation of these valleys was driven by two probably unrelated components. The first component is a phase of base-level lowering with magnitudes of∼300–500 m in the Coastal Cordillera. This period of base-level change in the Arica area, that started at ca. 7.5 Ma according to stratigraphic data, caused the trunk streams to dissect headward into the plains. The headward erosion interpretation is based on the presence of well-defined knickzones in stream profiles and the decrease in valley widths from the coast toward these knickzones. The second component is a change in paleoclimate. This interpretation is based on (1) the increase in the size of the largest alluvial boulders (from dm to m scale) with distal sources during the last 7.5 m.y., and (2) the calculated increase in minimum fluvial incision rates of ∼0.2 mm/yr between ca. 7.5 Ma and 3 Ma to ∼0.3 mm/yr subsequently. These trends suggest an increase in effective water discharge for systems sourced in the Western Cordillera (distal source). During the same time, however, valleys with headwaters in the coastal region (local source) lack any evidence of fluvial incision. This implies that the Coastal Cordillera became hyperarid sometime after 7.5 Ma. Furthermore, between 7.5 Ma and present, the sediment yields have been consistently higher in the catchments with distal sources (∼15 m/m.y.) than in the headwaters of rivers with local sources (<7 m/m.y.). The positive correlation between sediment yields and the altitude of the headwaters (distal versus local sources) seems to reflect the effect of orographic precipitation on surface erosion. It appears that base-level change in the coastal region, in combination with an increase in the orographic effect of precipitation, has controlled the topographic evolution of the northern Chilean Andes.

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments