Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Local Rates of Ice-Sheet Thickness Change in Greenland

The rate of ice-sheet thickness change is calculated for 10 sites in Greenland by comparing measured values of ice vertical velocity and snow-accumulation rate. Vertical velocities are derived from repeat surveys of markers using precision global positioning system techniques, and accumulation rates are determined from stratigraphic analysis of firn cores. The results apply to time-scales covered by the firn-core records, which in most cases are a few decades. A spectrum of thickness-change rates is obtained, ranging from substantial thinning to slow thickening. The sites where ice-sheet thinning is indicated are located near the ice-sheet margin or in outlet glacier catchments. Interior and high-elevation sites are predominantly in balance or thickening slowly. Uncertainties in the rates of thickness change are dominated by errors in the determination of accumulation rates. The results of this work are broadly comparable with regional estimates of mass balance obtained from the analysis of catchment input vs discharge.

Stratigraphic Continuity in 400 Mhz Short-Pulse Radar Profiles of Firn in West Antarctica

We track dated firn horizons within 400 MHz short-pulse radar profiles to find the continuous extent over which they can be used as historical benchmarks to study past accumulation rates in West Antarctica. The 30-40 cm pulse resolution compares with the accumulation rates of most areas. We tracked a particular set that varied from 30 to 90 m in depth over a distance of 600 km. The main limitations to continuity are fading at depth, pinching associated with accumulation rate differences within hills and valleys, and artificial fading caused by stacking along dips. The latter two may be overcome through multi-kilometer distances by matching the relative amplitude and spacing of several close horizons, along with their pulse forms and phases. Modeling of reflections from thin layers suggests that the - 37 to - 50 dB range of reflectivity and the pulse waveforms we observed are caused by the numerous thin ice layers observed in core stratigraphy. Constructive interference between reflections from these close, high-density layers can explain the maintenance of reflective strength throughout the depth of the firn despite the effects of compaction. The continuity suggests that these layers formed throughout West Antarctica and possibly into East Antarctica as well.

A 700-Year Record of Atmospheric Circulation Developed from the Law Dome Ice Core, East Antarctica

A 700-year, high-resolution, multivariate ice core record from Dome Summit South (DSS) (66degrees46'S, 112degrees48'E; 1370 m), Law Dome, is used to investigate sea level pressure (SLP) variability in the region of East Antarctica. Empirical orthogonal function (EOF) analysis reveals that the first EOF (LDEOF1) of the combined glaciochemical, oxygen isotope ratio, and accumulation rate record from DSS represents most of the variability in sea salt seen in the record. LDEOF1 is positively correlated (at least 95% confidence level) to instrumental June mean SLP across most of East Antarctica. Over the last 700 years, LDEOF1 levels at Law Dome were the highest during the nineteenth century, suggesting an increase in intensification of winter circulation during this period. The Law Dome DSS oxygen isotope ratio series also indicates that the nineteenth century had the coldest winters of any century in the record. In contrast, LDEOF1 levels were the lowest at Law Dome during the eighteenth century, suggesting a significant shift in the patterns and/or intensity of East Antarctic atmospheric circulation between the eighteenth and the nineteenth centuries. The LDEOF1 sea salt record is characterized by significant decadal-scale variability with a strong 25-year periodic structure.

Influence of Equatorial Diatom Processes on Si Deposition and Atmospheric CO(2) Cycles at Glacial/Interglacial Timescales

The causes of the glacial cycle remain unknown, although the primary driver is changes in atmospheric CO(2), likely controlled by the biological pump and biogeochemical cycles. The two most important regions of the ocean for exchange of CO(2) with the atmosphere are the equatorial Pacific and the Southern Ocean ( SO), the former a net source and the latter a net sink under present conditions. The equatorial Pacific has been shown to be a Si(OH)(4)-limited ecosystem, a consequence of the low source Si(OH)(4) concentrations in upwelled water that has its origin in the SO. This teleconnection for nutrients between the two regions suggests an oscillatory relationship that may influence or control glacial cycles. Opal mass accumulation rate (MAR) data and delta(15)N measurements in equatorial cores are interpreted with predictions from a one- dimensional Si(OH)(4)-limited ecosystem model (CoSINE) for the equatorial Pacific. The results suggest that equatorial Pacific surface CO(2) processes are in opposite phase to that of the global atmosphere, providing a negative feedback to the glacial cycle. This negative feedback is implemented through the effect of the SO on the equatorial Si(OH)(4) supply. An alternative hypothesis, that the whole ocean becomes Si(OH)(4) poor during cooling periods, is suggested by low opal MAR in cores from both equatorial and Antarctic regions, perhaps as a result of low river input. terminations in this scenario would result from blooms of coccolithophorids triggered by low Si(OH)(4) concentrations.

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.

Alpine lacustrine varved record reveals summer temperature as main control of glacier fluctuations over the past 2250 years

Glacier fluctuations are a key indicator of changing climate. Their reconstruction beyond historical times unravels glacier variability and its forcing factors on long time scales, which can considerably improve our understanding of the climate–glacier relationship. Here, we present a 2250-year-long reconstruction of particle-mass accumulation rates recorded in the lacustrine sediments of Lake Trüebsee (Central Swiss Alps) that are directly related to glacier extent, thus reflecting a continuous record of fluctuations of the upstream-located Titlis Glacier. Mass accumulation rate values show strong centennial to multi-centennial fluctuations and reveal 12 well-pronounced periods of enhanced values corresponding to times of maximum extent of the neighboring Lower Grindelwald Glacier. This result supports previous studies of proglacial lake sediments that documented high mass accumulation rate values during glacier advances. The strong variability in the Lake Trüebsee mass accumulation rate record thus represents a highly sensitive paleoclimatic archive, which mirrors rapid and pronounced feedbacks of Titlis Glacier to climatic changes over the past 2250years. The comparison of our data with independent paleo-temperature reconstructions from tree rings suggests that variations in mean summer temperature were the primary driving factor of fluctuations of Titlis Glacier. Also, advances of Titlis Glacier occurred during the grand solar minima (Dalton, Maunder, Spörer, Wolf) of the last millennium. This relation of glacier extent with summer temperature reveals strong evidence that the mass balance of this Alpine glacier is primarily controlled by the intensity of glacier melting during summer.

Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

Global pulses of organic carbon burial in deep-sea sediments during glacial maxima

The burial of organic carbon in marine sediments removes carbon dioxide from the ocean–atmosphere pool, provides energy to the deep biosphere, and on geological timescales drives the oxygenation of the atmosphere. Here we quantify natural variations in the burial of organic carbon in deep-sea sediments over the last glacial cycle. Using a new data compilation of hundreds of sediment cores, we show that the accumulation rate of organic carbon in the deep sea was consistently higher (50%) during glacial maxima than during interglacials. The spatial pattern and temporal progression of the changes suggest that enhanced nutrient supply to parts of the surface ocean contributed to the glacial burial pulses, with likely additional contributions from more efficient transfer of organic matter to the deep sea and better preservation of organic matter due to reduced oxygen exposure. These results demonstrate a pronounced climate sensitivity for this global carbon cycle sink.

Weather and Growing Season Summary

Precipitation for 2011 was less than the longterm climate average. Early in the year, precipitation lagged behind normal, but then tracked close to the normal accumulation rate from mid-April through mid-August. After that time, precipitation amounts greatly lagged behind normal, and the year ended almost 7 in. behind the long-term average. (Figure 1). Overall, 2011 will be remembered for good moisture early, but ending the season with almost no rainfall.

(Table 1) Summary of general lithology in IODP Exp302 cores

The Cenozoic ice-rafted debris (IRD) history of the central Arctic is reconstructed utilizing the terrigenous coarse sand fraction in IODP 302 cores from 0 to 273 meters composite depth. This Holocene - middle Eocene quantitative record of terrigenous sand accumulation on the Lomonosov Ridge, along with qualitative information on grain texture and composition, confirms the interpretation that ice initiation (sea ice and glacial ice) occurred ~46 Ma in the Arctic, and provides a long-term pattern of Arctic ice expansion and decay since the middle Eocene. IRD mass accumulation rates range from 0 to 0.13 g/cm2/ka in the middle Eocene and from 0 to 0.36 g/cm2/ka in the Neogene. IRD mass accumulation rate (MAR) maxima in the Miocene and Pliocene cooccur with either glacial initiation or intensification in the sub-Arctic. The 46.25 Ma IRD onset in the central Arctic slightly precedes the earliest evidence of ice in the Antarctic, and compares in timing with a >1000 ppm decrease in atmospheric concentrations of CO2. The decline of pCO2 in the middle Eocene may have driven both poles across the temperature threshold that enabled the nucleation of glaciers on land and partial freezing of the surface Arctic Ocean, especially during times of low insolation.

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.

Sedimentological and clay mineralogical analysis of ODP sites in the central Fram Strait

This study presents the results of high-resolution sedimentological and clay mineralogical investigations on sediments from ODP Sites 908A and 909AlC located in the central Fram Strait. The objective was to reconstruct the paleoclimate and paleoceanography of the high northern latitudes since the middle Miocene. The sediments are characterised in particular by a distinctive input of ice-rafted material, which most probably occurs since 6 Ma and very likely since 15 Ma. A change in the source area at 1 1.2 Ma is clearly marked by variations within clay mineral composition and increasing accumulation rates. This is interpreted as a result of an increase in water mass exchange through the Fram Strait. A further period of increasing exchange between 4-3 Ma is identified by granulometric investigations and points to a synchronous intensification of deep water production in the North Atlantic during this time interval. A comparison of the components of coarse and clay fraction clearly shows that both are not delivered by the Same transport process. The input of the clay fraction can be related to transport mechanisms through sea ice and glaciers and very likely also through oceanic currents. A reconstruction of source areas for clay minerals is possible only with some restrictions. High smectite contents in middle and late Miocene sediments indicate a background signal produced by soil formation together with sediment input, possibly originating from the Greenland- Scotland Ridge. The applicability of clay mineral distribution as a climate proxy for the high northern latitudes can be confirmed. Based on a comparison of sediments from Site 909C, characterised by the smectite/illite and chlorite ratio, with regional and global climatic records (oxygen isotopes), a middle Miocene cooling phase between 14.8-14.6 Ma can be proposed. A further cooling phase between 10-9 Ma clearly shows similarities in its Progress toward drastic decrease in carbonate sedimentation and preservation in the eastern equatorial Pacific. The modification in sea water and atmosphere chemistry may represent a possible link due to the built-up of equatorial carbonate platforms. Between 4.8-4.6 Ma clay mineral distribution indicates a distinct cooling trend in the Fram Strait region. This is not accompanied by relevant glaciation, which would otherwise be indicated by the coarse fraction. The intensification of glaciation in the northern hemisphere is distinctly documented by a rapid increase of illite and chlorite starting from 3.3 Ma, which corresponds to oxygen isotope data trends from North Atlantic.