Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands

This paper is the maritime and sub–Antarctic contribution to the Scientific Committee for Antarctic Research (SCAR) Past Antarctic Ice Sheet Dynamics (PAIS) community Antarctic Ice Sheet reconstruction. The overarching aim for all sectors of Antarctica was to reconstruct the Last Glacial Maximum (LGM) ice sheet extent and thickness, and map the subsequent deglaciation in a series of 5000 year time slices. However, our review of the literature found surprisingly few high quality chronological constraints on changing glacier extents on these timescales in the maritime and sub–Antarctic sector. Therefore, in this paper we focus on an assessment of the terrestrial and offshore evidence for the LGM ice extent, establishing minimum ages for the onset of deglaciation, and separating evidence of deglaciation from LGM limits from those associated with later Holocene glacier fluctuations. Evidence included geomorphological descriptions of glacial landscapes, radiocarbon dated basal peat and lake sediment deposits, cosmogenic isotope ages of glacial features and molecular biological data. We propose a classification of the glacial history of the maritime and sub–Antarctic islands based on this assembled evidence. These include: (Type I) islands which accumulated little or no LGM ice; (Type II) islands with a limited LGM ice extent but evidence of extensive earlier continental shelf glaciations; (Type III) seamounts and volcanoes unlikely to have accumulated significant LGM ice cover; (Type IV) islands on shallow shelves with both terrestrial and submarine evidence of LGM (and/or earlier) ice expansion; (Type V) Islands north of the Antarctic Polar Front with terrestrial evidence of LGM ice expansion; and (Type VI) islands with no data. Finally, we review the climatological and geomorphological settings that separate the glaciological history of the islands within this classification scheme.

Peak Last Glacial weathering intensity on the North American continent recorded by the authigenic Hf isotope composition of North Atlantic deep-sea sediments

We have retrieved radiogenic hafnium (Hf) isotope compositions (ɛHf) from authigenic Fe–Mn oxyhydroxides of deep northwest Atlantic sediments deposited over the past 26 ka to investigate the oceanic evidence of changes in dissolved weathering inputs from NE America during the last deglaciation. The extraction of seawater-derived Hf isotopic compositions from Fe–Mn oxyhydroxides is not a standard procedure. Comparisons between the Al/Hf ratios and Hf isotopic compositions of the chemically extracted authigenic phase on the one hand, and those of the corresponding detrital fractions on the other, provide evidence that the composition of past seawater has been reliably obtained for most sampled depths with our leaching procedures. This is endorsed most strongly by data for a sediment core from 4250 m water depth at the deeper Blake Ridge, for which consistent replicates were produced throughout. The Hf isotopic composition of the most recent sample in this core also closely matches that of nearby present day central North Atlantic seawater. Comparison with previously published seawater Nd and Pb isotope compositions obtained on the same cores shows that both Hf and Pb were released incongruently during incipient chemical weathering, but responded differently to the deglacial retreat of the Laurentide Ice Sheet. Hafnium was released more congruently during peak glacial conditions of the Last Glacial Maximum (LGM) and changed to typical incongruent interglacial ɛHf signatures either during or shortly after the LGM. This indicates that some zircon-derived Hf was released to seawater during the LGM. Conversely, there is no clear evidence for an increase in the influence of weathering of Lu-rich mineral phases during deglaciation, possibly since relatively unradiogenic Hf contributions from feldspar weathering were superimposed. While the authigenic Pb isotope signal in the same marine sediment samples traced peak chemical weathering rates on continental North America during the transition to the Holocene a similar incongruent excursion is notably absent in the Hf isotope record. The early change towards more radiogenic ɛHf in relation to the LGM may provide direct evidence for the transition from a cold-based to a warm-based Laurentide Ice Sheet on the Atlantic sector of North America.

RADIX: a minimal-resources rapid-access drilling system

Determining the expected age at a potential ice-core drilling site on a polar ice sheet generally depends on a combination of information from remote-sensing methods, estimates of current accumulation and modelling. This poses irreducible uncertainties in retrieving an undisturbed ice core of the desired age. Although recently perfected radar techniques will improve the picture of the ice sheet below future drilling sites, rapid prospective drillings could further increase the success of deep drilling projects. Here we design and explore a drilling system for a minimum-size rapid-access hole. The advantages of a small hole are the low demand for drilling fluid, low overall weight of the equipment, fast installing and de-installing and low costs. We show that, in theory, drilling of a 20mm hole to a depth of 3000m is possible in ~4 days. First concepts have been realized and verified in the field. Both the drill cuttings and the hole itself can be used to characterize the properties of the ice sheet and its potential to provide a trustworthy palaeo-record. A candidate drilling site could be explored in ~2 weeks, which would enable the characterization of several sites in one summer season.

Stability of the Atlantic meridional overturning circulation: A model intercomparison.

The evolution of the Atlantic Meridional Overturning Circulation (MOC) in 30 models of varying complexity is examined under four distinct Representative Concentration Pathways. The models include 25 Atmosphere-Ocean General Circulation Models (AOGCMs) or Earth System Models (ESMs) that submitted simulations in support of the 5th phase of the Coupled Model Intercomparison Project (CMIP5) and 5 Earth System Models of Intermediate Complexity (EMICs). While none of the models incorporated the additional effects of ice sheet melting, they all projected very similar behaviour during the 21st century. Over this period the strength of MOC reduced by a best estimate of 22% (18%–25%; 5%–95% confidence limits) for RCP2.6, 26% (23%–30%) for RCP4.5, 29% (23%–35%) for RCP6.0 and 40% (36%–44%) for RCP8.5. Two of the models eventually realized a slow shutdown of the MOC under RCP8.5, although no model exhibited an abrupt change of the MOC. Through analysis of the freshwater flux across 30°–32°S into the Atlantic, it was found that 40% of the CMIP5 models were in a bistable regime of the MOC for the duration of their RCP integrations. The results support previous assessments that it is very unlikely that the MOC will undergo an abrupt change to an off state as a consequence of global warming.

North Atlantic Eddy-Driven Jet in interglacial and glacial winter climates

The atmospheric westerly flow in the North Atlantic (NA) sector is dominated by atmospheric waves or eddies generating via momentum flux convergence, the so-called eddy-driven jet. The position of this jet is variable and shows for the present-day winter climate three preferred latitudinal states: a northern, central, and southernposition in the NA. Here, the authors analyze the behavior of the eddy-driven jet under different glacial and interglacial boundary conditions using atmosphere–land-only simulations with the CCSM4 climate model. As state-of-the-art climate models tend to underestimate the trimodality of the jet latitude, the authors apply a bias correction and successfully extract the trimodal behavior of the jet within CCSM4. The analysis shows that during interglacial times (i.e., the early Holocene and the Eemian) the preferred jet positions are rather stable and the observed multimodality is the typical interglacial character of the jet. During glacial times, the jet is strongly enhanced, its position is shifted southward, and the trimodal behavior vanishes. This is mainly due to the presence of the Laurentide ice sheet (LIS). The LIS enhances stationary waves downstream, thereby accelerating and displacing the NA eddy-driven jet by anomalous stationary momentum flux convergence. Additionally, changes in the transient eddy activity caused by topography changes as well as other glacial boundary conditions lead to an acceleration of the westerly winds over the southern NA at the expenseof more northernareas. Consequently, bothstationaryand transient eddiesfoster the southward shift of the NA eddy-driven jet during glacial winter times.

Links between atmospheric carbon dioxide, the land carbon reservoir and climate over the past millennium

The stability of terrestrial carbon reservoirs is thought to be closely linked to variations in climate 1, but the magnitude of carbon–climate feedbacks has proved dificult to constrain for both modern 2–4 and millennial 5–13 timescales. Reconstructions of atmospheric CO2 concentrations for the past thousand years have shown fluctuations on multidecadal to centennial timescales 5–7, but the causes of these fluctuations are unclear. Here we report high-resolution carbon isotope measurements of CO2 trapped within the ice of the West Antarctic Ice Sheet Divide ice core for the past 1,000 years. We use a deconvolution approach 14 to show that changes in terrestrial organic carbon stores best explain the observed multidecadal variations in the 13 C of CO2 and in CO2 concentrations from 755 to 1850 CE. If significant long-term carbon emissions came from pre-industrial anthropogenic land-use changes over this interval, the emissions must have been offset by a natural terrestrial sink for 13 C-depleted carbon, such as peatlands. We find that on multidecadal timescales, carbon cycle changes seem to vary with reconstructed regional climate changes. We conclude that climate variability could be an important control of fluctuations in land carbon storage on these timescales.

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.

Modern solar maximum forced late twentieth century Greenland cooling

The abrupt Northern Hemispheric warming at the end of the twentieth century has been attributed to an enhanced greenhouse effect. Yet Greenland and surrounding subpolar North Atlantic remained anomalously cold in 1970s to early 1990s. Here we reconstructed robust Greenland temperature records (North Greenland Ice Core Project and Greenland Ice Sheet Project 2) over the past 2100 years using argon and nitrogen isotopes in air trapped within ice cores and show that this cold anomaly was part of a recursive pattern of antiphase Greenland temperature responses to solar variability with a possible multidecadal lag. We hypothesize that high solar activity during the modern solar maximum (approximately 1950s–1980s) resulted in a cooling over Greenland and surrounding subpolar North Atlantic through the slowdown of Atlantic Meridional Overturning Circulation with atmospheric feedback processes.

Mid-Pliocene palynomorph record of the North Atlantic

The mid-Pliocene was an episode of prolonged global warmth and strong North Atlantic thermohaline circulation, interrupted briefly at circa 3.30 Ma by a global cooling event corresponding to marine isotope stage (MIS) M2. Paleoceanographic changes in the eastern North Atlantic have been reconstructed between circa 3.35 and 3.24 Ma at Deep Sea Drilling Project Site 610 and Integrated Ocean Drilling Program Site 1308. Mg/Ca ratios and d18O from Globigerina bulloides are used to reconstruct the temperature and relative salinity of surface waters, and dinoflagellate cyst assemblages are used to assess variability in the North Atlantic Current (NAC). Our sea surface temperature data indicate warm waters at both sites before and after MIS M2 but a cooling of ~2-3°C during MIS M2. A dinoflagellate cyst assemblage overturn marked by a decline in Operculodinium centrocarpum reflects a southward shift or slowdown of the NAC between circa 3.330 and 3.283 Ma, reducing northward heat transport 23-35 ka before the global ice volume maximum of MIS M2. This will have established conditions that ultimately allowed the Greenland ice sheet to expand, leading to the global cooling event at MIS M2. Comparison with an ice-rafted debris record excludes fresh water input via icebergs in the northeast Atlantic as a cause of NAC decline. The mechanism causing the temporary disruption of the NAC may be related to a brief reopening of the Panamanian Gateway at about this time.

Dinoflagellate cyst assemblages, oxygen isotope records of foraminifera, and alkenone and Mg/Ca-based SST estimates for the time interval 3.40-3.18 Ma (Late Pliocene) from five North Atlantic and Caribbean sediment cores

The early Late Pliocene (3.6 to ~3.0 million years ago) is the last extended interval in Earth's history when atmospheric CO2 concentrations were comparable to today's and global climate was warmer. Yet a severe global glaciation during marine isotope stage (MIS) M2 interrupted this phase of global warmth ~3.30 million years ago, and is seen as a premature attempt of the climate system to establish an ice-age world. Our geochemical and palynological records from five marine sediment cores along a Caribbean to eastern North Atlantic transect show that increased Pacific-to-Atlantic flow via the Central American Seaway weakened the North Atlantic Current (NAC) and attendant northward heat transport prior to MIS M2. The consequent cooling of the northern high latitude oceans permitted expansion of the Greenland ice sheet during MIS M2, despite near-modern atmospheric CO2 concentrations. Before and after MIS M2, heat transport via the NAC was crucial in maintaining warm climates comparable to those predicted for the end of this century.

Stable isotope analysis and temperature reconstruction data from DSDP Hole 94-609B and ODP Hole 162-984B

The Greenland ice sheet is accepted as a key factor controlling the Quaternary glacial scenario. However, the origin and mechanisms of major Arctic glaciation starting at 3.15 Ma and culminating at 2.74 Ma are still controversial. For this phase of intense cooling Ravelo et al. proposed a complex gradual forcing mechanism. In contrast, our new submillennial-scale paleoceanographic records from the Pliocene North Atlantic suggest a far more precise timing and forcing for the initiation of northern hemisphere glaciation (NHG), since it was linked to a 2-3 °C surface water warming during warm stages from 2.95 to 2.82 Ma. These records support previous models, claiming that the final closure of the Panama Isthmus (3.0- ~2.5 Ma induced an increased poleward salt and heat transport. Associated strengthening of North Atlantic Thermohaline Circulation and in turn, an intensified moisture supply to northern high latitudes resulted in the build-up of NHG, finally culminating in the great, irreversible climate crash at marine isotope stage G6 (2.74 Ma). In summary, there was a two-step threshold mechanism that marked the onset of NHG with glacial-to-interglacial cycles quasi-persistent until today.

Age determination and sea-surface reconstruction of sediment cores from the Norwegian Sea

Polycystine radiolarians are used to reconstruct summer sea surface temperatures (SSSTs) for the Late Pleistocene-Holocene (600-13,400 14C years BP) in the Norwegian Sea. At 13,200 14C years BP, the SSST was close to the average Holocene SSST (~12°C). It then gradually dropped to 7.1°C in the Younger Dryas. Near the Younger Dryas-Holocene transition (~10,000 14C years BP), the SSST increased 5°C in about 530 years. Four abrupt cooling events, with temperature drops of up to 2.1°C, are recognized during the Holocene: at 9340, 7100 ("8200 calendar years event"), 6400 and 1650 14C years BP. Radiolarian SSSTs and the isotopic signal from the GISP2 ice core are strongly coupled, stressing the importance of the Norwegian Sea as a mediator of heat/precipitation exchange between the North Atlantic, the atmosphere, and the Greenland ice sheet. Radiolarian and diatom-derived SSSTs display similarities, with the former not showing the recently reported Holocene cooling trend.

Diatom assemblages, HBIs and TEXL86 in sediment core NBP99-03_10

The West Antarctic ice sheet is particularly sensitive to global warming and its evolution and impact on global climate over the next few decades remains difficult to predict. In this context, investigating past sea ice conditions around Antarctica is of primary importance. Here, we document changes in sea ice presence, upper water column temperatures (0-200 m) and primary productivity over the last 9000 yr BP (before present) in the western Antarctic Peninsula (WAP) margin from a sedimentary core collected in the Palmer Deep Basin. Employing a multi-proxy approach, based on the combination of two biomarkers proxies (highly branched isoprenoid (HBI) alkenes for sea ice and TEXL86 for temperature) and micropaleontological data (diatom assemblages), we derived new Holocene records of sea ice conditions and upper water column temperatures. The early Holocene (9000-7000 yr BP) was characterized by a cooling phase with a short sea ice season. During the mid-Holocene (~7000-3800 yr BP), local climate evolved towards slightly colder conditions and a prominent extension of the sea ice season occurred, promoting a favorable environment for intensive diatom growth. The late Holocene (the last ~2100 yr) was characterized by warmer temperatures and increased sea ice presence, accompanied by reduced local primary productivity, likely in response to a shorter growing season compared to the early or mid-Holocene. The gradual increase in annual sea ice duration over the last 7000 yr might have been influenced by decreasing mean annual and spring insolation, despite increasing summer insolation. We postulate that, in addition to precessional changes in insolation, seasonal variability, via changes in the strength of the circumpolar Westerlies and upwelling activity, was further amplified by the increasing frequency/amplitude of the El Nino-Southern Oscillation (ENSO). However, between 3800 and 2100 yr BP, the lack of correlation between ENSO and climate variability in the WAP suggests that other climatic factors might have been more important in controlling WAP climate at this time.

Grain size analysis on CRP sediment cores from the Ross Sea, Antarctica

The upper 1200 m of pre-Pliocene sediment recovered by Cape Roberts Project (CRP) drilling off the Victoria Land coast of Antarctica between 1997-1999 has been subdivided into 54 unconformity-bound stratigraphic sequences, spanning the period c. 32 to 17 Ma. The sequences are recognised on the basis of the cyclical vertical stacking of their constituent lithofacies, which are enclosed by erosion surfaces produced during the grounding of the advancing ice margin onto the sea floor. Each sequence represents deposition in a range of offshore shelf to coastal glacimarine sedimentary environments during oscillations in the ice margin across the Western Ross Sea shelf, and coeval fluctuations in water depth. This paper applies spectral analysis techniques to depth- and time-series of sediment grain size (500 samples) for intervals of the core with adequate chronological data. Time series analysis of 0.5-l.0m-spaced grainsize data spanning sequences 9-11 (CRP-2/2A) and sequences 1-7 (CRP-3) suggests that the length of individual sequences correspond to Milankovitch frequencies, probably 41 k.y., but possibly as low as 100 k.y. Higher frequency periodic components at 23 k.y. (orbital precession) and 15-10 k.y. (sub-orbital) are recognised at the intrasequence-scale, and may represent climatic cycles akin to the ice rafting episodes described in the North Atlantic Ocean during the Quaternary. The cyclicity recorded by glacimarine sequences in CRP core provides direct evidence from the periphery of Antarctica for orbital oscillations in the size of the Oligocene-Early Miocene East Antarctic Ice Sheet.