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.

Diffusive equilibration of N2, O2 and CO2 mixing ratios in a 1.5-million-years-old ice core

In the framework of the International Partnerships in Ice Core Sciences, one of the most important targets is to retrieve an Antarctic ice core that extends over the last 1.5 million years (i.e. an ice core that enters the climate era when glacial–interglacial cycles followed the obliquity cycles of the earth). In such an ice core the annual layers of the oldest ice would be thinned by a factor of about 100 and the climatic information of a 10 000 yr interval would be contained in less than 1 m of ice. The gas record in such an Antarctic ice core can potentially reveal the role of greenhouse gas forcing on these 40 000 yr cycles. However, besides the extreme thinning of the annual layers, also the long residence time of the trapped air in the ice and the relatively high ice temperatures near the bedrock favour diffusive exchanges. To investigate the changes in the O2 / N2 ratio, as well as the trapped CO2 concentrations, we modelled the diffusive exchange of the trapped gases O2, N2 and CO2 along the vertical axis. However, the boundary conditions of a potential drilling site are not yet well constrained and the uncertainties in the permeation coefficients of the air constituents in the ice are large. In our simulations, we have set the drill site ice thickness at 2700 m and the bedrock ice temperature at 5–10 K below the ice pressure melting point. Using these conditions and including all further uncertainties associated with the drill site and the permeation coefficients, the results suggest that in the oldest ice the precessional variations in the O2 / N2 ratio will be damped by 50–100%, whereas CO2 concentration changes associated with glacial–interglacial variations will likely be conserved (simulated damping 5%). If the precessional O2 / N2 signal will have disappeared completely in this future ice core, orbital tuning of the ice-core age scale will be limited.

High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core

The Greenland NEEM (North Greenland Eemian Ice Drilling) operation in 2010 provided the first opportunity to combine trace-gas measurements by laser spectroscopic instruments and continuous-flow analysis along a freshly drilled ice core in a field-based setting. We present the resulting atmospheric methane (CH4) record covering the time period from 107.7 to 9.5 ka b2k (thousand years before 2000 AD). Companion discrete CH4 measurements are required to transfer the laser spectroscopic data from a relative to an absolute scale. However, even on a relative scale, the high-resolution CH4 data set significantly improves our knowledge of past atmospheric methane concentration changes. New significant sub-millennial-scale features appear during interstadials and stadials, generally associated with similar changes in water isotopic ratios of the ice, a proxy for local temperature. In addition to the midpoint of Dansgaard–Oeschger (D/O) CH4 transitions usually used for cross-dating, sharp definition of the start and end of these events brings precise depth markers (with ±20 cm uncertainty) for further cross-dating with other palaeo- or ice core records, e.g. speleothems. The method also provides an estimate of CH4 rates of change. The onsets of D/O events in the methane signal show a more rapid rate of change than their endings. The rate of CH4 increase associated with the onsets of D/O events progressively declines from 1.7 to 0.6 ppbv yr−1 in the course of marine isotope stage 3. The largest observed rate of increase takes place at the onset of D/O event #21 and reaches 2.5 ppbv yr−1.

High-resolution mineral dust and sea ice proxy records from the Talos Dome ice core

In this study we report on new non-sea salt calcium (nssCa2+, mineral dust proxy) and sea salt sodium (ssNa+, sea ice proxy) records along the East Antarctic Talos Dome deep ice core in centennial resolution reaching back 150 thousand years (ka) before present. During glacial conditions nssCa2+ fluxes in Talos Dome are strongly related to temperature as has been observed before in other deep Antarctic ice core records, and has been associated with synchronous changes in the main source region (southern South America) during climate variations in the last glacial. However, during warmer climate conditions Talos Dome mineral dust input is clearly elevated compared to other records mainly due to the contribution of additional local dust sources in the Ross Sea area. Based on a simple transport model, we compare nssCa2+ fluxes of different East Antarctic ice cores. From this multi-site comparison we conclude that changes in transport efficiency or atmospheric lifetime of dust particles do have a minor effect compared to source strength changes on the large-scale concentration changes observed in Antarctic ice cores during climate variations of the past 150 ka. Our transport model applied on ice core data is further validated by climate model data. The availability of multiple East Antarctic nssCa2+ records also allows for a revision of a former estimate on the atmospheric CO2 sensitivity to reduced dust induced iron fertilisation in the Southern Ocean during the transition from the Last Glacial Maximum to the Holocene (T1). While a former estimate based on the EPICA Dome C (EDC) record only suggested 20 ppm, we find that reduced dust induced iron fertilisation in the Southern Ocean may be responsible for up to 40 ppm of the total atmospheric CO2 increase during T1. During the last interglacial, ssNa+ levels of EDC and EPICA Dronning Maud Land (EDML) are only half of the Holocene levels, in line with higher temperatures during that period, indicating much reduced sea ice extent in the Atlantic as well as the Indian Ocean sector of the Southern Ocean. In contrast, Holocene ssNa+ flux in Talos Dome is about the same as during the last interglacial, indicating that there was similar ice cover present in the Ross Sea area during MIS 5.5 as during the Holocene.

Two-phase change in CO2, Antarctic temperature and global climate during Termination II

The end of the Last Glacial Maximum (Termination I), roughly 20 thousand years ago (ka), was marked by cooling in the Northern Hemisphere, a weakening of the Asian monsoon, a rise in atmospheric CO2 concentrations and warming over Antarctica. The sequence of events associated with the previous glacial–interglacial transition (Termination II), roughly 136 ka, is less well constrained. Here we present high-resolution records of atmospheric CO2 concentrations and isotopic composition of N2—an atmospheric temperature proxy—from air bubbles in the EPICA Dome C ice core that span Termination II. We find that atmospheric CO2 concentrations and Antarctic temperature started increasing in phase around 136 ka, but in a second phase of Termination II, from 130.5 to 129 ka, the rise in atmospheric CO2 concentrations lagged that of Antarctic temperature unequivocally. We suggest that during this second phase, the intensification of the low-latitude hydrological cycle resulted in the development of a CO2 sink, which counteracted the CO2 outgassing from the Southern Hemisphere oceans over this period.

Eemian and Holocene sea-surface conditions in the southern Black Sea: organic-walled dinoflagellate cyst record from core 22-GC3

In order to compare the sea-surface conditions in the Black Sea during the Holocene and Eemian, sapropelic parts of marine core 22-GC3 (42°13.53′N/36°29.55′E, 838 m water depth) were studied for organic-walled dinoflagellate cyst content. The record shows a change from freshwater/brackish assemblages (Pyxidinopsis psilata, Spiniferites cruciformis, and Caspidinium rugosum) to more marine assemblages (Lingulodinium machaerophorum and Spiniferites ramosus complex) during each interglacial, due to the inflow of saline Mediterranean water. The lacustrine–marine transitions in 22-GC3 occurred at ~ 8.3 cal kyr BP during the early Holocene and ~ 128 kyr BP during the early Eemian, slightly later compared to the onset of interglacial conditions on the adjacent land. Dinoflagellate cyst assemblages reveal higher sea-surface salinity (~ 28–30) (e.g. Spiniferites pachydermus, Bitectatodinium tepikiense, and Spiniferites mirabilis) around ~ 126.5–121 kyr BP in comparison to the Holocene (~ 15–20) as well as relatively high sea-surface temperature (e.g. Tuberculodinium vancampoae, S. pachydermus, and S. mirabilis) especially at ~ 127.6–125.3 kyr BP. Establishment of high sea-surface salinity during the Eemian correlates very well with reconstructed relatively high global sea-level and is explained as a combined effect of increased Mediterranean supply and high temperatures at the beginning of the last interglacial. The observed changes in the dinocyst record highlight the importance of nutrients for the composition of the Eemian and Holocene dinocyst assemblages.

Solar forcing of North Atlantic surface temperature and salinity over the past millennium

There were several centennial-scale fluctuations in the climate and oceanography of the North Atlantic region over the past 1,000 years, including a period of relative cooling from about AD 1450 to 1850 known as the Little Ice Age1. These variations may be linked to changes in solar irradiance, amplified through feedbacks including the Atlantic meridional overturning circulation2. Changes in the return limb of the Atlantic meridional overturning circulation are reflected in water properties at the base of the mixed layer south of Iceland. Here we reconstruct thermocline temperature and salinity in this region from AD 818 to 1780 using paired δ18O and Mg/Ca ratio measurements of foraminifer shells from a subdecadally resolved marine sediment core. The reconstructed centennial-scale variations in hydrography correlate with variability in total solar irradiance. We find a similar correlation in a simulation of climate over the past 1,000 years. We infer that the hydrographic changes probably reflect variability in the strength of the subpolar gyre associated with changes in atmospheric circulation. Specifically, in the simulation, low solar irradiance promotes the development of frequent and persistent atmospheric blocking events, in which a quasi-stationary high-pressure system in the eastern North Atlantic modifies the flow of the westerly winds. We conclude that this process could have contributed to the consecutive cold winters documented in Europe during the Little Ice Age.

A High-Resolution Record of Atmospheric Dust Composition and Variability Since Ad 1650 from a Mount Everest Ice Core

A Mount Everest ice core analyzed at high resolution for major and trace elements (Sr, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Tb, Dy, Ho, Er, Tm, Yb, Lu, Bi, U, Tl, Al, S, Ca, Ti, V, Cr, Mn, Fe, Co) and spanning the period A. D. 1650- 2002 is used to investigate the sources of and variations in atmospheric dust through time. The chemical composition of dust varies seasonally, and peak dust concentrations occur during the winter-spring months. Significant correlations between the Everest dust record and dust observations at stations suggest that the Everest record is representative of regional variations in atmospheric dust loading. Back-trajectory analysis in addition to a significant correlation of Everest dust concentrations and the Total Ozone Mapping Spectrometer (TOMS) aerosol index indicates that the dominant winter sources of dust are the Arabian Peninsula, Thar Desert, and northern Sahara. Factors that contribute to dust generation at the surface include soil moisture and temperature, and the long-range transport of dust aerosols appears to be sensitive to the strength of 500-mb zonal winds. There are periods of high dust concentration throughout the 350-yr Mount Everest dust record; however, there is an increase in these periods since the early 1800s. The record was examined for recent increases in dust emissions associated with anthropogenic activities, but no recent dust variations can be conclusively attributed to anthropogenic inputs of dust.

Twentieth Century Increase of Atmospheric Ammonia Recorded in Mount Everest Ice Core

An NH4+ record covering the period A.D. 1845-1997 was reconstructed using an 80.4 m ice core from East Rongbuk Glacier at an elevation of 6450 m on the northern slope of Mount Everest. Variations in NH4+ are characterized by a dramatic increase since the 1950s. The highest NH4+ concentrations occur in the 1980s. They are about twofold more than those in the first half of twentieth century. Empirical orthogonal function (EOF) analysis on the eight major ion (Na+,K+,Mg2+,NH4+,Ca2+,NO3-,SO42- and Cl-) series from this core indicates that NH4+ is loaded mainly on EOF3 (60% of NH4+ variance), suggesting that NH4+ has a unique signature. Instrumental sea level pressure (SLP) and regional temperatures are used to explore the relationship between NH4+ variations and both atmospheric circulation and natural source strength over Asia. Higher NH4+ concentrations are associated with an enhanced winter Mongolian High and a deepened summer Mongolian Low. A positive relationship also exists between NH4+ concentrations and regional temperature changes of the GIS Box 36 (Indian subcontinent), indicating that an increase in temperature may contribute to the strengthening of natural ammonia emissions (e. g., from plants and soils). A close positive correlation between NH4+ and acidic species (SO42- plus NO3-) concentrations suggests that a portion of the increase in NH4+ concentrations could be contributed by enhanced atmospheric acidification. Anthropogenic ammonia emissions from enhanced agricultural activities and energy consumption over Asia in concert with population increase since the 1950s appear also to be a significant factor in the dramatic increase of NH4+ concentrations during the last few decades.

Modeling Englacial Radar Attenuation at Siple Dome, West Antarctica, Using Ice Chemistry and Temperature Data

The radar reflectivity of an ice-sheet bed is a primary measurement for discriminating between thawed and frozen beds. Uncertainty in englacial radar attenuation and its spatial variation introduces corresponding uncertainty in estimates of basal reflectivity. Radar attenuation is proportional to ice conductivity, which depends on the concentrations of acid and sea-salt chloride and the temperature of the ice. We synthesize published conductivity measurements to specify an ice-conductivity model and find that some of the dielectric properties of ice at radar frequencies are not yet well constrained. Using depth profiles of ice-core chemistry and borehole temperature and an average of the experimental values for the dielectric properties, we calculate an attenuation rate profile for Siple Dome, West Antarctica. The depth-averaged modeled attenuation rate at Siple Dome (20.0 +/- 5.7 dB km(-1)) is somewhat lower than the value derived from radar profiles (25.3 +/- 1.1 dB km(-1)). Pending more experimental data on the dielectric properties of ice, we can match the modeled and radar-derived attenuation rates by an adjustment to the value for the pure ice conductivity that is within the range of reported values. Alternatively, using the pure ice dielectric properties derived from the most extensive single data set, the modeled depth-averaged attenuation rate is 24.0 +/- 2.2 dB km(-1). This work shows how to calculate englacial radar attenuation using ice chemistry and temperature data and establishes a basis for mapping spatial variations in radar attenuation across an ice sheet.

Oxygen Isotopic and Soluble Ionic Composition of a Shallow Firn Core, Inilchek Glacier, Central Tien Shan

Oxygen isotopic and soluble ionic measurements made on snow-pit (2 in depth) and firn-core (12.4 m depth samples recovered from the accumulation zone 5100 m) of Inilchek glacier 43degrees N, 79degrees E) provide information on recent (1992-98) climatic and environmental conditions in the central Tien Shan region of central Asia. The combined 14.4 m snow-pit/firn-core profile lies within the firn zone, arid contains only one observed melt feature (10 m temperature = - 12 degreesC), Although some post-depositional attenuation of the sub-seasonal delta(18)O record is possible, annual cycles are apparent throughout the isotope profile. We therefore use the preserved delta(18)O record to establish a depth/age scale for the core. Mean delta(18)O values for the entire core and for summer periods are consistent with delta(18)O/temperature observations, and suggest the delta(18)O record provides a means to reconstruct past changes in summer surface temperature at the site. Major-ion (Na(+), K(+), Mg(2+), Ca(2+), NH(4)(+), Cl(-), NO(3)(-), SO(4)(2-)) data from the core demonstrate the dominant influence of dust deposition on the soluble chemistry at the site, arid indicate significant interannual variability in atmospheric-dust loading during the 1900s. Anthropogenic impacts oil NH(4)(+) concentrations are observed at the site, and suggest a summer increase in atmospheric NH(4)(+) that may be related to regional agricultural (nitrogen-rich fertilizer use activities.

Atmospheric Soluble Dust Records from a Tibetan Ice Core: Possible Climate Proxies and Teleconnection With the Pacific Decadal Oscillation

In autumn 2005, a joint expedition between the University of Maine and the Institute of Tibetan Plateau Research recovered three ice cores from Guoqu Glacier (33 degrees 34'37.80 '' N, 91 degrees 10'35.3 '' E, 5720 m above sea level) on the northern side of Mt. Geladaindong, central Tibetan Plateau. Isotopes ( delta(18)O), major soluble ions (Na(+), K(+), Mg(2+), Ca(2+), Cl(-), NO(3)(-), SO(4)(2-)), and radionuclide (beta-activity) measurements from one of the cores revealed a 70-year record (1935-2005). Statistical analysis of major ion time series suggests that atmospheric soluble dust species dominate the chemical signature and that background dust levels conceal marine ion species deposition. The soluble dust time series have interspecies relations and common structure (empirical orthogonal function (EOF) 1), suggesting a similar soluble dust source or transport route. Annual and seasonal correlations between the EOF 1 time series and National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis climate variables (1948-2004) suggest that the Mt. Geladaindong ice core record provides a proxy for local and regional surface pressure. An approximately threefold decrease of soluble dust concentrations in the middle to late 1970s, accompanied by regional increases in pressure and temperature and decreases in wind velocity, coincides with the major 1976-1977 shift of the Pacific Decadal Oscillation (PDO) from a negative to a positive state. This is the first ice core evidence of a potential teleconnection between central Asian atmospheric soluble dust loading and the PDO. Analysis of temporally longer ice cores from Mt. Geladaindong may enhance understanding of the relationship between the PDO and central Asian atmospheric circulation and subsequent atmospheric soluble dust loading.

Association Between Atmospheric Circulation Patterns and Firn-Ice Core Records from the Inilchek Glacierized Area, Central Tien Shan, Asia

Glacioclimatological research in the central Tien Shan was performed in the summers of 1998 and 1999 on the South Inilchek Glacier at 5100 - 5460 m. A 14.36 m firn-ice core and snow samples were collected and used for stratigraphic, isotopic, and chemical analyses. The firn-ice core and snow records were related to snow pit measurements at an event scale and to meteorological data and synoptic indices of atmospheric circulation at annual and seasonal scales. Linear relationships between the seasonal air temperature and seasonal isotopic composition in accumulated precipitation were established. Changes in the delta(18)O air temperature relationship, in major ion concentration and in the ratios between chemical species, were used to identify different sources of moisture and investigate changes in atmospheric circulation patterns. Precipitation over the central Tien Shan is characterized by the lowest ionic content among the Tien Shan glaciers and indicates its mainly marine origin. In seasons of minimum precipitation, autumn and winter, water vapor was derived from the arid and semiarid regions in central Eurasia and contributed annual maximal solute content to snow accumulation in Tien Shan. The lowest content of major ions was observed in spring and summer layers, which represent maximum seasonal accumulation when moisture originates over the Atlantic Ocean and Mediterranean and Black Seas.

A Record of Atmospheric Co2 During the Last 40,000 Years from the Siple Dome, Antarctica Ice Core

We have measured the CO2 concentration of air occluded during the last 40,000 years in the deep Siple Dome A ( hereafter Siple Dome) ice core, Antarctica. The general trend of CO2 concentration from Siple Dome ice follows the temperature inferred from the isotopic composition of the ice and is mostly in agreement with other Antarctic ice core CO2 records. CO2 rose initially at similar to 17.5 kyr B. P. ( thousand years before 1950), decreased slowly during the Antarctic Cold Reversal, rose during the Younger Dryas, fell to a local minimum at around 8 kyr B. P., and rose continuously since then. The CO2 concentration never reached steady state during the Holocene, as also found in the Taylor Dome and EPICA Dome C ( hereafter Dome C) records. During the last glacial termination, a lag of CO2 versus Siple Dome isotopic temperature is probable. The Siple Dome CO2 concentrations during the last glacial termination and in the Holocene are at certain times greater than in other Antarctic ice cores by up to 20 ppm (mumol CO2/mol air). While in situ production of CO2 is one possible cause of the sporadic elevated levels, the mechanism leading to the enrichment is not yet clear.

Dependence of Eemian Greenland temperature reconstructions on the ice sheet topography

The influence of a reduced Greenland Ice Sheet (GrIS) on Greenland's surface climate during the Eemian interglacial is studied using a set of simulations with different GrIS realizations performed with a comprehensive climate model. We find a distinct impact of changes in the GrIS topography on Greenland's surface air temperatures (SAT) even when correcting for changes in surface elevation, which influences SAT through the lapse rate effect. The resulting lapse-rate-corrected SAT anomalies are thermodynamically driven by changes in the local surface energy balance rather than dynamically caused through anomalous advection of warm/cold air masses. The large-scale circulation is indeed very stable among all sensitivity experiments and the Northern Hemisphere (NH) flow pattern does not depend on Greenland's topography in the Eemian. In contrast, Greenland's surface energy balance is clearly influenced by changes in the GrIS topography and this impact is seasonally diverse. In winter, the variable reacting strongest to changes in the topography is the sensible heat flux (SHF). The reason is its dependence on surface winds, which themselves are controlled to a large extent by the shape of the GrIS. Hence, regions where a receding GrIS causes higher surface wind velocities also experience anomalous warming through SHF. Vice-versa, regions that become flat and ice-free are characterized by low wind speeds, low SHF, and anomalous low winter temperatures. In summer, we find surface warming induced by a decrease in surface albedo in deglaciated areas and regions which experience surface melting. The Eemian temperature records derived from Greenland proxies, thus, likely include a temperature signal arising from changes in the GrIS topography. For the Eemian ice found in the NEEM core, our model suggests that up to 3.1 °C of the annual mean Eemian warming can be attributed to these topography-related processes and hence is not necessarily linked to large-scale climate variations.