New constraints on the gas age-ice age difference along the EPICA ice cores, 0-50 kyr

Gas is trapped in polar ice sheets at ~50–120 m below the surface and is therefore younger than the surrounding ice. Firn densification models are used to evaluate this ice age-gas age difference (Δage) in the past. However, such models need to be validated by data, in particular for periods colder than present day on the East Antarctic plateau. Here we bring new constraints to test a firn densification model applied to the EPICA Dome C (EDC) site for the last 50 kyr, by linking the EDC ice core to the EPICA Dronning Maud Land (EDML) ice core, both in the ice phase (using volcanic horizons) and in the gas phase (using rapid methane variations). We also use the structured 10Be peak, occurring 41 kyr before present (BP) and due to the low geomagnetic field associated with the Laschamp event, to experimentally estimate the Δage during this event. Our results seem to reveal an overestimate of the Δage by the firn densification model during the last glacial period at EDC. Tests with different accumulation rates and temperature scenarios do not entirely resolve this discrepancy. Although the exact reasons for the Δage overestimate at the two EPICA sites remain unknown at this stage, we conclude that current densification model simulations have deficits under glacial climatic conditions. Whatever the cause of the Δage overestimate, our finding suggests that the phase relationship between CO2 and EDC temperature previously inferred for the start of the last deglaciation (lag of CO2 by 800±600 yr) seems to be overestimated.

Synchronisation of the EDML and EDC ice cores for the last 52 kyr by volcanic signature matching

A common time scale for the EPICA ice cores from Dome C (EDC) and Dronning Maud Land (EDML) has been established. Since the EDML core was not drilled on a dome, the development of the EDML1 time scale for the EPICA ice core drilled in Dronning Maud Land was based on the creation of a detailed stratigraphic link between EDML and EDC, which was dated by a simpler 1D ice-flow model. The synchronisation between the two EPICA ice cores was done through the identification of several common volcanic signatures. This paper describes the rigorous method, using the signature of volcanic sulfate, which was employed for the last 52 kyr of the record. We estimated the discrepancies between the modelled EDC and EDML glaciological age scales during the studied period, by evaluating the ratio R of the apparent duration of temporal intervals between pairs of isochrones. On average R ranges between 0.8 and 1.2 corresponding to an uncertainty of up to 20% in the estimate of the time duration in at least one of the two ice cores. Significant deviations of R up to 1.4–1.5 are observed between 18 and 28 kyr before present (BP), where present is defined as 1950. At this stage our approach does not allow us unequivocally to find out which of the models is affected by errors, but assuming that the thinning function at both sites and accumulation history at Dome C (which was drilled on a dome) are correct, this anomaly can be ascribed to a complex spatial accumulation variability (which may be different in the past compared to the present day) upstream of the EDML core.

Iron fluxes to Talos Dome, Antarctica, over the past 200 kyr

Atmospheric fluxes of iron (Fe) over the past 200 kyr are reported for the coastal Antarctic Talos Dome ice core, based on acid leachable Fe concentrations. Fluxes of Fe to Talos Dome were consistently greater than those at Dome C, with the greatest difference observed during interglacial climates. We observe different Fe flux trends at Dome C and Talos Dome during the deglaciation and early Holocene, attributed to a combination of deglacial activation of dust sources local to Talos Dome and the reorganisation of atmospheric transport pathways with the retreat of the Ross Sea ice shelf. This supports similar findings based on dust particle sizes and fluxes and Rare Earth Element fluxes. We show that Ca and Fe should not be used as quantitative proxies for mineral dust, as they all demonstrate different deglacial trends at Talos Dome and Dome C. Considering that a 20 ppmv decrease in atmospheric CO2 at the coldest part of the last glacial maximum occurs contemporaneously with the period of greatest Fe and dust flux to Antarctica, we confirm that the maximum contribution of aeolian dust deposition to Southern Ocean sequestration of atmospheric CO2 is approximately 20 ppmv.�

NGRIP CH4 concentration from 120 to 10 kyr before present and its relation to a δ15N temperature reconstruction from the same ice core

During the last glacial cycle, Greenland temperature showed many rapid temperature variations, the so-called Dansgaard–Oeschger (DO) events. The past atmospheric methane concentration closely followed these temperature variations, which implies that the warmings recorded in Greenland were probably hemispheric in extent. Here we substantially extend and complete the North Greenland Ice Core Project (NGRIP) methane record from the Preboreal Holocene (PB) back to the end of the last interglacial period with a mean time resolution of 54 yr. We relate the amplitudes of the methane increases associated with DO events to the amplitudes of the local Greenland NGRIP temperature increases derived from stable nitrogen isotope (δ15N) measurements, which have been performed along the same ice core (Kindler et al., 2014). We find the ratio to oscillate between 5 parts per billion (ppb) per °C and 18 ppb °C−1 with the approximate frequency of the precessional cycle. A remarkably high ratio of 25.5 ppb °C−1 is reached during the transition from the Younger Dryas (YD) to the PB. Analysis of the timing of the fast methane and temperature increases reveals significant lags of the methane increases relative to NGRIP temperature for DO events 5, 9, 10, 11, 13, 15, 19, and 20. These events generally have small methane increase rates and we hypothesize that the lag is caused by pronounced northward displacement of the source regions from stadial to interstadial. We further show that the relative interpolar concentration difference (rIPD) of methane is about 4.5% for the stadials between DO events 18 and 20, which is in the same order as in the stadials before and after DO event 2 around the Last Glacial Maximum. The rIPD of methane remains relatively stable throughout the full last glacial, with a tendency for elevated values during interstadial compared to stadial periods.

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.

Depositional modes and lake-level variability at Lake Towuti, Indonesia, during the past ~29 kyr BP

Lake Towuti (2.5°S, 121.5°E) is a long-lived, tectonic lake located on the Island of Sulawesi, Indonesia, and in the center of the Indo-Pacific warm pool (IPWP). Lake Towuti is connected with upstream lakes Matano and Mahalona through the Mahalona River, which constitutes the largest inlet to the lake. The Mahalona River Delta is prograding into Lake Towuti’s deep northern basin thus exerting significant control on depositional processes in the basin. We combine high-resolution seismic reflection and sedimentological datasets from a 19.8-m-long sediment piston core from the distal edge of this delta to characterize fluctuations in deltaic sedimentation during the past ~29 kyr BP and their relation to climatic change. Our datasets reveal that, in the present, sedimentation is strongly influenced by deposition of laterally transported sediments sourced from the Mahalona River Delta. Variations in the amount of laterally transported sediments, as expressed by coarse fraction amounts in pelagic muds and turbidite recurrence rates and cumulative thicknesses, are primarily a function of lake-level induced delta slope instability and delta progradation into the basin. We infer lowest lake-levels between ~29 and 16, a gradual lake level rise between ~16 and 11, and high lake-levels between ~11 and 0 kyr BP. Periods of highest turbidite deposition, ~26 to 24 and ~18 to 16 kyr BP coincide with Heinrich events 2 and 1, respectively. Our lake-level reconstruction therefore supports previous observations based on geochemical hydroclimate proxies of a very dry last glacial and a wet Holocene in the region, and provides new evidence of millennial-scale variations in moisture balance in the IPWP.

Paleohydraulic reconstruction of a 40 kyr-old terrace sequence implies that water discharge was larger than today

The evolution of landscapes crucially depends on the climate history. This is particularly evident in South America where landscape responses to orbital climate shifts have been well documented. However, while most studies have focused on inferring temperature variations from paleoclimate proxy data, estimates of water budget changes have been complicated because of a lack of adequate physical information. Here, we present a methodology and related results, which allowed us to extract water discharge values from the sedimentary record of the 40 Ka-old fluvial terrace deposits in the Pisco valley, western Peru. In particular, this valley hosts a Quaternary cut-and-fill succession that we used, in combination with beryllium-10 (10Be)-based sediment flux, gauging records, channel geometries and grain size measurements, to quantitatively assess sediment and water discharge values c. 40 Ka ago in relation to present-day conditions. We compare these discharge estimates to the discharge regime of the modern Pisco River and find that the water discharge of the paleo-Pisco River, during the Minchin pluvial period c. 40 Ka ago, was c. 7–8 times greater than the modern Pisco River if considering the mean and the maximum water discharge. In addition, the calculations show that inferred water discharge estimates are mainly dependent on channel gradients and grain size values, and to a lesser extent on channel width measures. Finally, we found that the c. 40 Ka-old Minchin terrace material was poorer sorted than the modern deposits, which might reflect that sediment transport during the past period was characterized by a larger divergence from equal mobility compared to the modern situation. In summary, the differences in grain size distribution and inferred water discharge estimates between the modern and the paleo-Pisco River suggests that the 40 Ka-old Minchin period was characterized by a wetter climate and more powerful flood events.

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.

Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed ~300 ppm at any point during the last 800 000 years, East Antarctica was at least ~3–4 °C warmer than preindustrial (CO2~280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (~3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as "Warmer than Present Transients" (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations preceding WPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming.

Rapid parallel adaptive radiations from a single hybridogenic ancestral population

The Alpine lake whitefish (Coregonus lavaretus) species complex is a classic example of a recent radiation, associated with colonization of the Alpine lakes following the glacial retreat (less than 15 kyr BP). They have formed a unique array of endemic lake flocks, each with one to six described sympatric species differing in morphology, diet and reproductive ecology. Here, we present a genomic investigation of the relationships between and within the lake flocks. Comparing the signal between over 1000 AFLP loci and mitochondrial control region sequence data, we use phylogenetic tree-based and population genetic methods to reconstruct the phylogenetic history of the group and to delineate the principal centres of genetic diversity within the radiation. We find significant cytonuclear discordance showing that the genomically monophyletic Alpine whitefish clade arose from a hybrid swarm of at least two glacial refugial lineages. Within this radiation, we find seven extant genetic clusters centred on seven lake systems. Most interestingly, we find evidence of sympatric speciation within and parallel evolution of equivalent phenotypes among these lake systems. However, we also find the genetic signature of human-mediated gene flow and diversity loss within many lakes, highlighting the fragility of recent radiations.