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.

Pinus nigra (European black pine) as the dominant species of the last glacial pinewoods in south-western to central Iberia: a morphological study of modern and fossil pollen

Aim Our aim was to discriminate different species of Pinus via pollen analysis in order to assess the responses of particular pine species to orbital and millennial-scale climate changes, particularly during the last glacial period. Location Modern pollen grains were collected from current pine populations along transects from the Pyrenees to southern Iberia and the Balearic Islands. Fossil pine pollen was recovered from the south-western Iberian margin core MD95-2042. Methods We measured a set of morphological traits of modern pollen from the Iberian pine species Pinus nigra, P. sylvestris, P. halepensis, P. pinea and P. pinaster and of fossil pine pollen from selected samples of the last glacial period and the early to mid-Holocene. Classification and regression tree (CART) analysis was used to establish a model from the modern dataset that discriminates pollen from the different pine species and allows identification of fossil pine pollen at the species level. Results The CART model was effective in separating pollen of P. nigra and P. sylvestris from that of the Mediterranean pine group (P. halepensis, P. pinea and P. pinaster). The pollen of Pinus nigra diverged from that of P. sylvestris by having a more flattened corpus. Predictions using this model suggested that fossil pine pollen is mainly from P. nigra in all the samples analysed. Pinus sylvestris was more abundant in samples from Greenland stadials than Heinrich stadials, whereas Mediterranean pines increased in samples from Greenland interstadials and during the early to mid-Holocene. Main conclusions Morphological parameters can be successfully used to increase the taxonomic resolution of fossil pine pollen at the species level for the highland pines (P. nigra and P. sylvestris) and at the group of species level for the Mediterranean pines. Our study indicates that P. nigra was the dominant component of the last glacial south-western/central Iberian pinewoods, although the species composition of these woodlands varied in response to abrupt climate changes.

Chemical compositions of solid particles present in the Greenland NEEM ice core over the last 110,000 years

This study reports the chemical composition of particles present along Greenland’s North Greenland Eemian Ice Drilling (NEEM) ice core, back to 110,000 years before present. Insoluble and soluble particles larger than 0.45 μm were extracted from the ice core by ice sublimation, and their chemical composition was analyzed using scanning electron microscope and energy dispersive X-ray spectroscopy and micro-Raman spectroscopy. We show that the dominant insoluble components are silicates, whereas NaCl, Na₂SO₄, CaSO ₄, and CaCO₃ represent major soluble salts. For the first time, particles of CaMg(CO₃)₂ and Ca(NO₃)₂ 4H₂O are identified in a Greenland ice core. The chemical speciation of salts varies with past climatic conditions. Whereas the fraction of Na salts (NaCl + Na₂SO₄) exceeds that of Ca salts (CaSO₄+ CaCO₃) during the Holocene (0.6–11.7 kyr B.P.), the two fractions are similar during the Bølling-Allerød period (12.9–14.6 kyr B.P.). During cold climate such as over the Younger Dryas (12.0–12.6 kyr B.P.) and the Last Glacial Maximum (15.0–26.9 kyr B.P.), the fraction of Ca salts exceeds that of Na salts, showing that the most abundant ion generally controls the salt budget in each period. High-resolution analyses reveal changing particle compositions: those in Holocene ice show seasonal changes, and those in LGM ice show a difference between cloudy bands and clear layers, which again can be largely explained by the availability of ionic components in the atmospheric aerosol body of air masses reaching Greenland.

Cool Glacial Temperatures and Changes in Moisture Source Recorded in Oman Groundwaters

Concentrations of atmospheric noble gases (neon, argon, krypton, and xenon) dissolved in groundwaters from northern Oman indicate that the average ground temperature during the Late Pleistocene (15,000 to 24,000 years before present) was 6.5° ± 0.6°C lower than that of today. Stable oxygen and hydrogen isotopic groundwater data show that the origin of atmospheric water vapor changed from a primarily southern, Indian Ocean source during the Late Pleistocene to a dominantly northern, Mediterranean source today. The reduced northern water vapor source is consistent with a drier Last Glacial Maximum through much of northern Africa and Arabia.

Covariation of deep Southern Ocean oxygenation and atmospheric CO2 through the last ice age

No single mechanism can account for the full amplitude of past atmospheric carbon dioxide (CO2) concentration variability over glacial–interglacial cycles. A build-up of carbon in the deep ocean has been shown to have occurred during the Last Glacial Maximum. However, the mechanisms responsible for the release of the deeply sequestered carbon to the atmosphere at deglaciation, and the relative importance of deep ocean sequestration in regulating millennial-timescale variations in atmospheric CO2 concentration before the Last Glacial Maximum, have remained unclear. Here we present sedimentary redox-sensitive trace metal records from the Antarctic Zone of the Southern Ocean that provide a reconstruction of transient changes in deep ocean oxygenation and, by inference, respired carbon storage throughout the last glacial cycle. Our data suggest that respired carbon was removed from the abyssal Southern Ocean during the Northern Hemisphere cold phases of the deglaciation, when atmospheric CO2 concentration increased rapidly, reflecting—at least in part— a combination of dwindling iron fertilization by dust and enhanced deep ocean ventilation. Furthermore, our records show that the observed covariation between atmospheric CO2 concentration and abyssal Southern Ocean oxygenation was maintained throughout most of the past 80,000 years. This suggests that on millennial timescales deep ocean circulation and iron fertilization in the Southern Ocean played a consistent role in modifying atmospheric CO2 concentration.

Plant-wax, alkenone and major element analyses of ODP Site 659 for the Last Glacial cycle and associated coretop data

We present a hydrologic reconstruction of the Sahara-Sahel transition, covering the complete last glacial cycle (130 ka), based on a combination of plant-wax-specific hydrogen (dD) and carbon isotopes (d13C). The dD and d13C signatures of long-chain n-alkanes from ODP Site 659 off NW Africa reveal a significant anti-correlation. Complementary to published pollen data, we infer that this plant-wax signal reflects sensitive responses of the vegetation cover to precipitation changes in the Sahel region, as well as varying contributions from biomes north of the Sahara (C3 domain) by North-East Trade Winds (NETW). During arid phases, especially the northern parts of the Sahel likely experienced crucial water stress, which resulted in a pronounced contraction of the vegetation cover, thus reducing the amount of C4 plant waxes from the region. The increase in NETW strength during dry periods further promoted a more pronounced C3-plant-wax signal derived from the North African C3 plant domain. During humid periods, the C4-dominated Sahelian environments spread northward into the Saharan realm, in association with lower NETW inputs of C3 plant waxes. Arid-humid cycles deduced from plant-wax dD are in accordance with concomitant changes in weathering intensity reflected in varying major element distributions. Environmental shifts are generally linked to periods with large fluctuations in Northern Hemisphere summer insolation. During Marine Isotope Stages 2 and 3, when insolation variability was low, coupling of the hydrologic regime to alkenone-based estimates of NE Atlantic sea-surface temperatures becomes apparent.

Changes in bottom water corrosiveness and carbonate ion concentrations in the sub-Antarctic Atlantic during the last glacial period

The glacial climate system transitioned rapidly between cold (stadial) and warm (interstadial) conditions in the Northern Hemisphere. This variability, referred to as Dansgaard-Oeschger variability, is widely believed to arise from perturbations of the Atlantic Meridional Overturning Circulation. Evidence for such changes during the longer Heinrich stadials has been identified, but direct evidence for overturning circulation changes during Dansgaard-Oeschger events has proven elusive. Here we reconstruct bottom water [CO3]2- variability from B/Ca ratios of benthic foraminifera and indicators of sedimentary dissolution, and use these reconstructions to infer the flow of northern-sourced deep water to the deep central sub-Antarctic Atlantic Ocean. We find that nearly every Dansgaard-Oeschger interstadial is accompanied by a rapid incursion of North Atlantic Deep Water into the deep South Atlantic. Based on these results and transient climate model simulations, we conclude that North Atlantic stadial-interstadial climate variability was associated with significant Atlantic overturning circulation changes that were rapidly transmitted across the Atlantic. However, by demonstrating the persistent role of Atlantic overturning circulation changes in past abrupt climate variability, our reconstructions of carbonate chemistry further indicate that the carbon cycle response to abrupt climate change was not a simple function of North Atlantic overturning.

Mean sortable silt values from Holocene and glacial Scotia Sea sediments

The Antarctic Circumpolar Current is key to the mixing and ventilation of the world's oceans. This current flows from west to east between about 45° and 70° S connecting the Atlantic, Pacific and Indian oceans, and is driven by westerly winds and buoyancy forcing. High levels of productivity in the current regulate atmospheric CO2 concentrations. Reconstructions of the current during the last glacial period suggest that flow speeds were faster or similar to present, and it is uncertain whether the strength and position of the westerly winds changed. Here we reconstruct Antarctic Circumpolar Current bottom speeds through the constricting Drake Passage and Scotia Sea during the Last Glacial Maximum and Holocene based on the mean grain size of sortable silt from a suite of sediment cores. We find essentially no change in bottom flow speeds through the region, and, given that the momentum imparted by winds, and modulated by sea-ice cover, is balanced by the interaction of these flows with the seabed, this argues against substantial changes in wind stress. However, glacial flow speeds in the sea-ice zone south of 56° S were significantly slower than present, whereas flow in the north was faster, but not significantly so. We suggest that slower flow over the rough topography south of 56° S may have reduced diapycnal mixing in this region during the last glacial period, possibly reducing the diapycnal contribution to the Southern Ocean overturning circulation.

Mineral dust deposition in interglacial and glacial climate conditions: model results

The global aerosol/climate model ECHAM5-HAM is used in order to investigate the dust cycle for four interglacial and one glacial climate conditions. The 20-year time-slices are the pre-industrial control (CTRL), mid-Holocene (6000 years BP), last glacial inception (115000 years BP), Eemian (126000 years BP) and Last Glacial Maximum (LGM) (21000 years BP) time intervals. The study is focused on the Antarctic region. The model is able to reproduce the magnitude order of dust deposition globally for the pre-industial and LGM climates. Correlation coefficient of the natural logarithm of the observed and modeled values is 0.78 for the CTRL and 0.81 for the LGM. For the pre-industrial simulation the model overestimates observed values in Antarctica by a factor of about 2-3 due to overestimation of the Australian dust source and too high wet deposition in the Antarctica interior. In the LGM, the model underestimates dust deposition in eastern Antarctica by a factor of about 4-5 due to underestimation of the South American dust source. More records are needed to validate dust deposition for the past interglacial time-slices. The modeled results show that dust deposition in Antarctica in the past interglacial time-slices is higher than in the CTRL simulation. The largest increase of dust deposition in Antarctica is simulated for the LGM, showing about 10-fold increase compared to CTRL.

Glacial-interglacial record of ODP Site 181-1119 in the Canterbury Bright, Southwest Pacific Ocean

Ocean Drilling Program (ODP) Site 1119 is located at water depth 395 m near the subtropical front (STF; here represented by the Southland Front), just downslope from the shelf edge of eastern South Island, New Zealand. The upper 86.19 metres composite depth (mcd) of Site 1119 sediment was deposited at an average sedimentation rate of 34 cm/kyr during Marine Isotope Stages (MIS) 1-8 (0-252 ka), and is underlain across a ~25 kyr intra-MIS 8 unconformity by MIS 8.5-11 (277-367 ka) and older sediment deposited at ~14 cm/kyr. A time scale is assigned to Site 1119 using radiocarbon dates for the period back to ~39 ka, and, prior to then, by matching its climatic record with that of the Vostok ice core, which it closely resembles. Four palaeoceanographic proxy measures for surface water masses vary together with the sandy-muddy, glacial-interglacial (G/I) cyclicity at the site. Interglacial intervals are characterised by heavy delta13C, high colour reflectance (a proxy for carbonate content), low Q-ray (a proxy for clay content) and light delta18O; conversely, glacial intervals exhibit light delta13C, low reflectance, high Q-ray and heavy delta18O signatures. Early interglacial intervals are represented by silty clays with 10-105-cm-thick beds of sharp-based (Chondrites-burrowed), shelly, graded, fine sand. The sands are rich in foraminifera, and were deposited distant from the shoreline under the influence of longitudinal flow in relatively deep water. Glacial intervals comprise mostly micaceous silty clay, though with some thin (2-10 cm thick) sands present also at peak cold periods, and contain the cold-water scallop Zygochlamys delicatula. Interglacial sandy intervals are characterised by relatively low sedimentation rates of 5-32 cm/kyr; cold climate intervals MIS 10, 6 and 2 have successively higher sedimentation rates of 45, 69 and 140 cm/kyr. Counter-intuitively,and forced by the bathymetric control of a laterally-moving shoreline during G/I and I/G transitions, the 1119 core records a southeasterly (seaward) movement of the STF during early glacial periods, accompanied by the incursion of subtropical water (STW) above the site, and northwesterly (landward) movement during late glacial and interglacial times, resulting in a dominant influence then of subantarctic surface water (SAW). The history of passage of these different water masses at the site is clearly delineated by their characteristic delta13C values. The intervals of thin, graded sands-muds which occur within MIS 2-3, 6, 7.4 and 10 indicate the onset at times of peak cold of intermittent bottom currents caused by strengthened and expanded frontal flows along the STF, which at such times lay near Site 1119 in close proximity to seaward-encroaching subantarctic waters within the Bounty gyre. In common with other nearby Southern Hemisphere records, the cold period which represents the last glacial maximum lasted between ~23-18 ka at Site 1119, during which time the STF and Subantarctic Front (SAF) probably merged into a single intense frontal zone around the head of the adjacent Bounty Trough.