The impact of different glacial boundary conditions on atmospheric dynamics and precipitation in the North Atlantic region

Using a highly resolved atmospheric general circulation model, the impact of different glacial boundary conditions on precipitation and atmospheric dynamics in the North Atlantic region is investigated. Six 30-yr time slice experiments of the Last Glacial Maximum at 21 thousand years before the present (ka BP) and of a less pronounced glacial state – the Middle Weichselian (65 ka BP) – are compared to analyse the sensitivity to changes in the ice sheet distribution, in the radiative forcing and in the prescribed time-varying sea surface temperature and sea ice, which are taken from a lower-resolved, but fully coupled atmosphere-ocean general circulation model. The strongest differences are found for simulations with different heights of the Laurentide ice sheet. A high surface elevation of the Laurentide ice sheet leads to a southward displacement of the jet stream and the storm track in the North Atlantic region. These changes in the atmospheric dynamics generate a band of increased precipitation in the mid-latitudes across the Atlantic to southern Europe in winter, while the precipitation pattern in summer is only marginally affected. The impact of the radiative forcing differences between the two glacial periods and of the prescribed time-varying sea surface temperatures and sea ice are of second order importance compared to the one of the Laurentide ice sheet. They affect the atmospheric dynamics and precipitation in a similar but less pronounced manner compared with the topographic changes.

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.

Response of subfossil Cladocera in Gerzensee (Swiss Plateau) to early Late Glacial environmental change

Sub-fossil Cladocera were studied in a core from Gerzensee (Swiss Plateau) for the late-glacial periods of Oldest Dryas, Bølling, and Allerød. Cladocera assemblages were dominated by cold-tolerant littoral taxa Chydorus sphaericus, Acroperus harpae, Alonella nana, Alona affinis, and Alonella excisa. The rapid warming at the beginning of the Bølling (GI-1e) ca. 14,650 yr before present (BP: before AD 1950) was indicated by an abrupt 2‰ shift in carbonate δ18O and a clear change in pollen assemblages. Cladocera assemblages, in contrast, changed more gradually. C. sphaericus and A. harpae are the most cold-tolerant, and their abundance was highest in the earliest part of the record. Only 150–200 years after the beginning of the Bølling warming we observed an increase in less cold-tolerant A. excisa and A. affinis. The establishment of Alona guttata, A. guttata var. tuberculata, and Pleuroxus unicatus was delayed by ca. 350, 770, and 800 years respectively after the onset of the Bølling. The development of the Cladocera assemblages suggests increasing water temperatures during the Bølling/Allerød, which agrees with the interpretation by von Grafenstein et al. (2013-this issue) that decreasing δ18O values in carbonates in this period reflect increasing summer water temperatures at the sediment–water interface. Other processes also affected the Cladocera community, including the development and diversification of aquatic vegetation favourable for Cladocera. The record is clearly dominated by Chydoridae, as expected for a littoral core. Yet, the planktonic Eubosmina-group occurred throughout the core, with the exception of a period at ca. 13,760–13,420 yr BP. Lake levels reconstructed for this period are relatively low, indicating that the littoral location might have become too shallow for Eubosmina in that period.

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.

Centennial mineral dust variability in high-resolution ice core data from Dome C, Antarctica

Ice core data from Antarctica provide detailed insights into the characteristics of past climate, atmospheric circulation, as well as changes in the aerosol load of the atmosphere. We present high-resolution records of soluble calcium (Ca2+), non-sea-salt soluble calcium (nssCa2+), and particulate mineral dust aerosol from the East Antarctic Plateau at a depth resolution of 1 cm, spanning the past 800 000 years. Despite the fact that all three parameters are largely dust-derived, the ratio of nssCa2+ to particulate dust is dependent on the particulate dust concentration itself. We used principal component analysis to extract the joint climatic signal and produce a common high-resolution record of dust flux. This new record is used to identify Antarctic warming events during the past eight glacial periods. The phasing of dust flux and CO2 changes during glacial-interglacial transitions reveals that iron fertilization of the Southern Ocean during the past nine glacial terminations was not the dominant factor in the deglacial rise of CO2 concentrations. Rapid changes in dust flux during glacial terminations and Antarctic warming events point to a rapid response of the southern westerly wind belt in the region of southern South American dust sources on changing climate conditions. The clear lead of these dust changes on temperature rise suggests that an atmospheric reorganization occurred in the Southern Hemisphere before the Southern Ocean warmed significantly.

Biogeochemical variability during the past 3.6 million years recorded by FTIR spectroscopy in the sediment record of Lake El'gygytgyn, Far East Russian Arctic

Abstract. A number of studies have shown that Fourier transform infrared spectroscopy (FTIRS) can be applied to quantitatively assess lacustrine sediment constituents. In this study, we developed calibration models based on FTIRS for the quantitative determination of biogenic silica (BSi; n = 420; gradient: 0.9–56.5 %), total organic carbon (TOC; n = 309; gradient: 0–2.9 %), and total inorganic carbon (TIC; n = 152; gradient: 0–0.4 %) in a 318 m-long sediment record with a basal age of 3.6 million years from Lake El’gygytgyn, Far East Russian Arctic. The developed partial least squares (PLS) regression models yield high cross-validated (CV) R2 CV = 0.86–0.91 and low root mean square error of crossvalidation (RMSECV) (3.1–7.0% of the gradient for the different properties). By applying these models to 6771 samples from the entire sediment record, we obtained detailed insight into bioproductivity variations in Lake El’gygytgyn throughout the middle to late Pliocene and Quaternary. High accumulation rates of BSi indicate a productivity maximum during the middle Pliocene (3.6–3.3 Ma), followed by gradually decreasing rates during the late Pliocene and Quaternary. The average BSi accumulation during the middle Pliocene was �3 times higher than maximum accumulation rates during the past 1.5 million years. The indicated progressive deterioration of environmental and climatic conditions in the Siberian Arctic starting at ca. 3.3 Ma is consistent with the first occurrence of glacial periods and the finally complete establishment of glacial–interglacial cycles during the Quaternary.

Phylogeography of silver pheasant (Lophura nycthemera L.) across China: aggregate effects of refugia, introgression, and riverine barriers

The role of Pleistocene glacial cycles in forming the contemporary genetic structure of organisms has been well studied in China with a particular focus on the Tibetan Plateau. However, China has a complex topography and diversity of local climates, and how glacial cycles may have shaped the subtropical and tropical biota of the region remains mostly unaddressed. To investigate the factors that affected the phylogeography and population history of a widely distributed and nondeciduous forest species, we analysed morphological characters, mitochondrial DNA sequences and nuclear microsatellite loci in the Silver Pheasant (Lophura nycthemera). In a pattern generally consistent with phenotypic clusters, but not nominal subspecies, deeply divergent mitochondrial lineages restricted to different geographic regions were detected. Coalescent simulations indicated that the time of main divergence events corresponded to major glacial periods in the Pleistocene and gene flow was only partially lowered by drainage barriers between some populations. Intraspecific cytonuclear discordance was revealed in mitochondrial lineages from Hainan Island and the Sichuan Basin with evidence of nuclear gene flow from neighbouring populations into the latter. Unexpectedly, hybridization was revealed in Yingjiang between the Silver Pheasant and Kalij Pheasant (Lophura leucomelanos) with wide genetic introgression at both the mtDNA and nuclear levels. Our results highlight a novel phylogeographic pattern in a subtropical area generated from the combined effects of climate oscillation, partial drainage barriers and interspecific hybridization. Cytonuclear discordance combined with morphological differentiation implies that complex historical factors shaped the divergence process in this biodiversity hot spot area of southern China.

Improving accuracy and precision of ice core δD(CH4) analyses using methane pre-pyrolysis and hydrogen post-pyrolysis trapping and subsequent chromatographic separation

Firn and polar ice cores offer the only direct palaeoatmospheric archive. Analyses of past greenhouse gas concentrations and their isotopic compositions in air bubbles in the ice can help to constrain changes in global biogeochemical cycles in the past. For the analysis of the hydrogen isotopic composition of methane (δD(CH4) or δ2H(CH4)) 0.5 to 1.5 kg of ice was hitherto used. Here we present a method to improve precision and reduce the sample amount for δD(CH4) measurements in (ice core) air. Pre-concentrated methane is focused in front of a high temperature oven (pre-pyrolysis trapping), and molecular hydrogen formed by pyrolysis is trapped afterwards (post-pyrolysis trapping), both on a carbon-PLOT capillary at −196 °C. Argon, oxygen, nitrogen, carbon monoxide, unpyrolysed methane and krypton are trapped together with H2 and must be separated using a second short, cooled chromatographic column to ensure accurate results. Pre- and post-pyrolysis trapping largely removes the isotopic fractionation induced during chromatographic separation and results in a narrow peak in the mass spectrometer. Air standards can be measured with a precision better than 1‰. For polar ice samples from glacial periods, we estimate a precision of 2.3‰ for 350 g of ice (or roughly 30 mL – at standard temperature and pressure (STP) – of air) with 350 ppb of methane. This corresponds to recent tropospheric air samples (about 1900 ppb CH4) of about 6 mL (STP) or about 500 pmol of pure CH4.

Sedimentological processes and environmental variability at Lake Ohrid (Macedonia, Albania) between 637 ka and the present

Lake Ohrid (Macedonia, Albania) is thought to be more than 1.2 million years old and host more than 300 endemic species. As a target of the International Continental scientific Drilling Program (ICDP), a successful deep drilling campaign was carried out within the scope of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project in 2013. Here, we present lithological, sedimentological, and (bio-)geochemical data from the upper 247.8 m composite depth of the overall 569 m long DEEP site sediment succession from the central part of the lake. According to an age model, which is based on 11 tephra layers (first-order tie points) and on tuning of bio-geochemical proxy data to orbital parameters (second-order tie points), the analyzed sediment sequence covers the last 637 kyr. The DEEP site sediment succession consists of hemipelagic sediments, which are interspersed by several tephra layers and infrequent, thin (< 5 cm) mass wasting deposits. The hemipelagic sediments can be classified into three different lithotypes. Lithotype 1 and 2 deposits comprise calcareous and slightly calcareous silty clay and are predominantly attributed to interglacial periods with high primary productivity in the lake during summer and reduced mixing during winter. The data suggest that high ion and nutrient concentrations in the lake water promoted calcite precipitation and diatom growth in the epilimnion during MIS15, 13, and 5. Following a strong primary productivity, highest interglacial temperatures can be reported for marine isotope stages (MIS) 11 and 5, whereas MIS15, 13, 9, and 7 were comparably cooler. Lithotype 3 deposits consist of clastic, silty clayey material and predominantly represent glacial periods with low primary productivity during summer and longer and intensified mixing during winter. The data imply that the most severe glacial conditions at Lake Ohrid persisted during MIS16, 12, 10, and 6, whereas somewhat warmer temperatures can be inferred for MIS14, 8, 4, and 2. Interglacial-like conditions occurred during parts of MIS14 and 8.

Northern Mediterranean climate since the Middle Pleistocene: a 637 ka stable isotope record from Lake Ohrid (Albania/Macedonia)

Lake Ohrid (Macedonia/Albania) is an ancient lake with unique biodiversity and a site of global significance for investigating the influence of climate, geological, and tectonic events on the generation of endemic populations. Here, we present oxygen (δ18O) and carbon (δ13C) isotope data from carbonate over the upper 243 m of a composite core profile recovered as part of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project. The investigated sediment succession covers the past ca. 637 ka. Previous studies on short cores from the lake (up to 15 m, < 140 ka) have indicated the total inorganic carbon (TIC) content of sediments to be highly sensitive to climate change over the last glacial–interglacial cycle. Sediments corresponding to warmer periods contain abundant endogenic calcite; however, an overall low TIC content in glacial sediments is punctuated by discrete bands of early diagenetic authigenic siderite. Isotope measurements on endogenic calcite (δ18Oc and δ13Cc) reveal variations both between and within interglacials that suggest the lake has been subject to palaeoenvironmental change on orbital and millennial timescales. We also measured isotope ratios from authigenic siderite (δ18Os and δ13Cs) and, with the oxygen isotope composition of calcite and siderite, reconstruct δ18O of lake water (δ18Olw) over the last 637 ka. Interglacials have higher δ18Olw values when compared to glacial periods most likely due to changes in evaporation, summer temperature, the proportion of winter precipitation (snowfall), and inflow from adjacent Lake Prespa. The isotope stratigraphy suggests Lake Ohrid experienced a period of general stability from marine isotope stage (MIS) 15 to MIS 13, highlighting MIS 14 as a particularly warm glacial. Climate conditions became progressively wetter during MIS 11 and MIS 9. Interglacial periods after MIS 9 are characterised by increasingly evaporated and drier conditions through MIS 7, MIS 5, and the Holocene. Our results provide new evidence for long-term climate change in the northern Mediterranean region, which will form the basis to better understand the influence of major environmental events on biological evolution within Lake Ohrid.

The Southern Hemisphere at glacial terminations: insights from the Dome C ice core

The many different proxy records from the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core allow for the first time a comparison of nine glacial terminations in great detail. Despite the fact that all terminations cover the transition from a glacial maximum into an interglacial, there are large differences between single terminations. For some terminations, Antarctic temperature increased only moderately, while for others, the amplitude of change at the termination was much larger. For the different terminations, the rate of change in temperature is more similar than the magnitude or duration of change. These temperature changes were accompanied by vast changes in dust and sea salt deposition all over Antarctica. Here we investigate the phasing between a South American dust proxy (non-sea-salt calcium flux, nssCa2+), a sea ice proxy (sea salt sodium flux, ssNa+) and a proxy for Antarctic temperature (deuterium, δD). In particular, we look into whether a similar sequence of events applies to all terminations, despite their different characteristics. All proxies are derived from the EPICA Dome C ice core, resulting in a relative dating uncertainty between the proxies of less than 20 years. At the start of the terminations, the temperature (δD) increase and dust (nssCa2+ flux) decrease start synchronously. The sea ice proxy (ssNa+ flux), however, only changes once the temperature has reached a particular threshold, approximately 5°C below present day temperatures (corresponding to a δD value of −420‰). This reflects to a large extent the limited sensitivity of the sea ice proxy during very cold periods with large sea ice extent. At terminations where this threshold is not reached (TVI, TVIII), ssNa+ flux shows no changes. Above this threshold, the sea ice proxy is closely coupled to the Antarctic temperature, and interglacial levels are reached at the same time for both ssNa+ and δD. On the other hand, once another threshold at approximately 2°C below present day temperature is passed (corresponding to a δD value of −402‰), nssCa2+ flux has reached interglacial levels and does not change any more, despite further warming. This threshold behaviour most likely results from a combination of changes to the threshold friction velocity for dust entrainment and to the distribution of surface wind speeds in the dust source region.

A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy

Due to their outstanding resolution and well-constrained chronologies, Greenland ice-core records provide a master record of past climatic changes throughout the Last Interglacial–Glacial cycle in the North Atlantic region. As part of the INTIMATE (INTegration of Ice-core, MArine and TErrestrial records) project, protocols have been proposed to ensure consistent and robust correlation between different records of past climate. A key element of these protocols has been the formal definition and ordinal numbering of the sequence of Greenland Stadials (GS) and Greenland Interstadials (GI) within the most recent glacial period. The GS and GI periods are the Greenland expressions of the characteristic Dansgaard–Oeschger events that represent cold and warm phases of the North Atlantic region, respectively. We present here a more detailed and extended GS/GI template for the whole of the Last Glacial period. It is based on a synchronization of the NGRIP, GRIP, and GISP2 ice-core records that allows the parallel analysis of all three records on a common time scale. The boundaries of the GS and GI periods are defined based on a combination of stable-oxygen isotope ratios of the ice (δ18O, reflecting mainly local temperature) and calcium ion concentrations (reflecting mainly atmospheric dust loading) measured in the ice. The data not only resolve the well-known sequence of Dansgaard–Oeschger events that were first defined and numbered in the ice-core records more than two decades ago, but also better resolve a number of short-lived climatic oscillations, some defined here for the first time. Using this revised scheme, we propose a consistent approach for discriminating and naming all the significant abrupt climatic events of the Last Glacial period that are represented in the Greenland ice records. The final product constitutes an extended and better resolved Greenland stratotype sequence, against which other proxy records can be compared and correlated. It also provides a more secure basis for investigating the dynamics and fundamental causes of these climatic perturbations.

Millennial changes in North American wildfire and soil activity over the last glacial cycle

Climate changes in the North Atlantic region during the last glacial cycle were dominated by the slow waxing and waning of the North American ice sheet as well as by intermittent, millennial-scale Dansgaard–Oeschger climate oscillations. However, prior to the last deglaciation, the responses of North American vegetation and biomass burning to these climate variations are uncertain. Ammonium in Greenland ice cores, a product from North American soil emissions and biomass burning events, can help to fill this gap. Here we use continuous, high-resolution measurements of ammonium concentrations between 110,000 to 10,000 years ago from the Greenland NGRIP and GRIP ice cores to reconstruct North American wildfire activity and soil ammonium emissions. We find that on orbital timescales soil emissions increased under warmer climate conditions when vegetation expanded northwards into previously ice-covered areas. For millennial-scale interstadial warm periods during Marine Isotope Stage 3, the fire recurrence rate increased in parallel to the rapid warmings, whereas soil emissions rose more slowly, reflecting slow ice shrinkage and delayed ecosystem changes. We conclude that sudden warming events had little impact on soil ammonium emissions and ammonium transport to Greenland, but did result in a substantial increase in the frequency of North American wildfires.

A Late-glacial and Holocene record of climatic change from a Swiss peat humification profile

Colorimetric measurements of alkaline extracts from two Swiss peat cores have provided a complete 14500-year-long record of peat humification, a proxy of effective precipitation. Peat from the cold Younger Dryas (11050–9550 cal. bc) was well preserved despite low levels of precipitation. A particularly dry period, peaking at c. 7100 cal. bc, is indicated by well-decomposed peat. Peat from c. 6750–4250 cal. bc shows a low degree of decomposition, indicating a wet bog surface despite relatively warm temperatures and therefore indicating high levels of precipitation. A sharp transition to higher levels of decomposition c. 4450–3750 cal. bc indicates a major transition to a drier bog surface. Subsequently, peat humification generally decreases towards the end of the deeper profile (c. cal. ad 1050), indicating a gradual return to wetter conditions. This gradual decrease is punctuated by periods of particularly low humification which appear to be due to shifts to higher levels of effective precipitation from c. 2500 to 1350 cal. bc, c. 1050 to 550 cal. bc, centered around 150 cal. bc, and from c. cal. ad 550 onwards. Anthropogenic influences appear to have affected peat humification at the site at least since the Middle Ages. This study indicates that humification in colder regions/time periods could be more affected by temperature than precipitation and vice versa.