scholarly journals Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ<sup>15</sup>N measurements

2013 ◽  
Vol 9 (3) ◽  
pp. 983-999 ◽  
Author(s):  
E. Capron ◽  
A. Landais ◽  
D. Buiron ◽  
A. Cauquoin ◽  
J. Chappellaz ◽  
...  

Abstract. 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.

2012 ◽  
Vol 8 (6) ◽  
pp. 6051-6091 ◽  
Author(s):  
E. Capron ◽  
A. Landais ◽  
D. Buiron ◽  
A. Cauquoin ◽  
J. Chappellaz ◽  
...  

Abstract. Correct estimate of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice cores studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: output of a firn densification model and measurements of δ15N of N2 in air trapped in ice core. 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 δ15N measurements performed from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rate and temperature conditions. While firn densification simulations are able to correctly represent most of the δ15N trends over the last deglaciation measured in the EDC, BI, TALDICE and EDML ice cores, they systematically fail to capture BI and EDML δ15N glacial levels, a mismatch previously seen for Central East Antarctic ice cores. Using empirical constraints of the EDML gas-ice depth offset during the Laschamp event (~ 41 ka), we can rule out the existence of a large convective zone as the explanation of the glacial firn model-δ15N data mismatch for this site. The good match between modelled and measured δ15N at TALDICE as well as the lack of any clear correlation between insoluble dust concentration in snow and δ15N records in the different ice cores suggest that past changes in loads of impurities are not the only main driver of glacial-interglacial changes in firn lock-in depth. We conclude that firn densification dynamics may instead be driven mostly by accumulation rate changes. The mismatch between modelled and measured δ15N may be due to inaccurate reconstruction of past accumulation rate or underestimated influence of accumulation rate in firnification models.


2007 ◽  
Vol 3 (3) ◽  
pp. 527-540 ◽  
Author(s):  
L. Loulergue ◽  
F. Parrenin ◽  
T. Blunier ◽  
J.-M. Barnola ◽  
R. Spahni ◽  
...  

Abstract. 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.


2018 ◽  
Author(s):  
Amaëlle Landais ◽  
Emilie Capron ◽  
Valérie Masson-Delmotte ◽  
Samuel Toucanne ◽  
Rachael Rhodes ◽  
...  

Abstract. The last deglaciation represents the most recent example of natural global warming associated with large-scale climate changes. In addition to the long-term global temperature increase, the last deglaciation onset is punctuated by a sequence of abrupt changes in the Northern Hemisphere. Such interplay between orbital- and millennial-scale variability is widely documented in paleoclimatic records but the underlying mechanisms are not fully understood. Limitations arise from the difficulty in constraining the sequence of events between external forcing, high- and low- latitude climate and environmental changes. Greenland ice cores provide sub-decadal-scale records across the last deglaciation and contain fingerprints of climate variations occurring in different regions of the Northern Hemisphere. Here, we combine new ice d-excess and 17O-excess records, tracing changes in the mid-latitudes, with ice δ18O records of polar climate. Within Heinrich Stadial 1, we demonstrate a decoupling between climatic conditions in Greenland and those of the lower latitudes. While Greenland temperature remains mostly stable from 17.5 to 14.7 ka, significant change in the mid latitudes of northern Atlantic takes place at ~ 16.2 ka, associated with warmer and wetter conditions of Greenland moisture sources. We show that this climate modification is coincident with abrupt changes in atmospheric CO2 and CH4 concentrations recorded in an Antarctic ice core. Our coherent ice core chronological framework and comparison with other paleoclimate records suggests a mechanism involving two-step freshwater fluxes in the North Atlantic associated with a southward shift of the intertropical convergence zone.


2020 ◽  
Author(s):  
Loïc Schmidely ◽  
Lucas Silva ◽  
Christoph Nehrbass-Ahles ◽  
Juhyeong Han ◽  
Jinhwa Shin ◽  
...  

&lt;p&gt; Small air inclusions in ice cores represent a direct archive of past atmospheric compositions, allowing us to measure the concentration of the three most potent non-condensable Greenhouse Gases (GHG) CO&lt;sub&gt;2&lt;/sub&gt;, CH&lt;sub&gt;4&lt;/sub&gt; and N&lt;sub&gt;2&lt;/sub&gt;O as far back as 800,000 years before present (kyr BP). These records demonstrate that transitions from glacial to interglacial conditions are accompanied by a substantial net increase of CO&lt;sub&gt;2&lt;/sub&gt;, CH&lt;sub&gt;4&lt;/sub&gt; and N&lt;sub&gt;2&lt;/sub&gt;O in the atmosphere (L&amp;#252;thi et al. 2008, Loulergue et al. 2008, Schilt et al. 2010). A sound understanding of the interplay between the reorganization of the climate system and the perturbation of GHG inventories during glacial terminations is partly limited by the temporal resolution of the records derived from ice cores. In fact, with the exception of the last deglaciation (23-9 kyr BP) centennial-scale GHG variability remained uncaptured for precedings glacial terminations.&lt;/p&gt;&lt;p&gt;In this work, we exploit the exceptionally long temporal coverage of the EPICA Dome C (EDC) ice core to reconstruct, for the first time, centennial-scale fluctuations of CH&lt;sub&gt;4&lt;/sub&gt; mole fractions from 145 to 125 kyr BP, encompassing the entire penultimate deglaciation (138-128 kyr BP). With a temporal resolution of ~100 years, our new record is now unveiling all climate-driven signals enclosed into the EDC ice core, exploiting the maximum resolution possible at Dome C (). This offers us the opportunity to study the timing and rates of change of CH&lt;sub&gt;4&lt;/sub&gt; in unprecedented details.&lt;/p&gt;&lt;p&gt;Preliminary analysis reveals that the deglacial CH&lt;sub&gt;4 &lt;/sub&gt;rise is a superimposition of gradual millennial-scale increases (~0.01-0.02 ppb/year) and abrupt and partly intermittent centennial-scale events (~80-200 ppb in less than a millennium). We will investigate processes modulating the observed changes in the CH&lt;sub&gt;4&lt;/sub&gt; cycle, compare the structure of our record with the CH&lt;sub&gt;4&lt;/sub&gt; profile of the last deglaciation (Marcott, 2014) and contrast it with the EDC CO&lt;sub&gt;2&lt;/sub&gt; and N&lt;sub&gt;2&lt;/sub&gt;O records over the penultimate glacial termination now available in similar resolution.&lt;/p&gt;


2011 ◽  
Vol 7 (1) ◽  
pp. 1-16 ◽  
Author(s):  
D. Buiron ◽  
J. Chappellaz ◽  
B. Stenni ◽  
M. Frezzotti ◽  
M. Baumgartner ◽  
...  

Abstract. A new deep ice core drilling program, TALDICE, has been successfully handled by a European team at Talos Dome, in the Ross Sea sector of East Antarctica, down to 1620 m depth. Using stratigraphic markers and a new inverse method, we produce the first official chronology of the ice core, called TALDICE-1. We show that it notably improves an a priori chronology resulting from a one-dimensional ice flow model. It is in agreement with a posteriori controls of the resulting accumulation rate and thinning function along the core. An absolute uncertainty of only 300 yr is obtained over the course of the last deglaciation. This uncertainty remains lower than 600 yr over Marine Isotope Stage 3, back to 50 kyr BP. The phasing of the TALDICE ice core climate record with respect to the central East Antarctic plateau and Greenland records can thus be determined with a precision allowing for a discussion of the mechanisms at work at sub-millennial time scales.


2017 ◽  
Author(s):  
Léa Gest ◽  
Frédéric Parrenin ◽  
Jai Chowdhry Beeman ◽  
Dominique Raynaud ◽  
Tyler J. Fudge ◽  
...  

Abstract. To understand causal relationships in past climate variations, it is essential to have accurate chronologies of paleoclimate records. The last deglaciation, which occurred from 18 000 to 11 000 years ago, is especially interesting, since it is the most recent large climatic variation of global extent. Ice cores in Antarctica provide important paleoclimate proxies, such as regional temperature and global atmospheric CO2. However, temperature is recorded in the ice while CO2 is recorded in the enclosed air bubbles. The ages of the former and of the latter are different since air is trapped at 50–120 m below the surface. It is therefore necessary to correct for this air-ice shift to accurately infer the sequence of events. Here we accurately determine the phasing between East Antarctic temperature and atmospheric CO2 variations during the last deglacial warming based on Antarctic ice core records. We build a stack of East Antarctic temperature variations by averaging the records from 4 ice cores (EPICA Dome C, Dome Fuji, EPICA Dronning Maud Land and Talos Dome), all accurately synchronized by volcanic event matching. We place this stack onto the WAIS Divide WD2014 age scale by synchronizing EPICA Dome C and WAIS Divide using volcanic event matching, which allows comparison with the high resolution CO2 record from WAIS Divide. Since WAIS Divide is a high accumulation site, its air age scale, which has previously been determined by firn modeling, is more robust. Finally, we assess the CO2/Antarctic temperature phasing by determining four periods when their trends change abruptly. We find that at the onset of the last deglaciation and at the onset of the Antarctic Cold Reversal (ACR) period CO2 and Antarctic temperature are synchronous within a range of 210 years. Then CO2 slightly leads by 165 ± 116 years at the end of the Antarctic Cold Reversal (ACR) period. Finally, Antarctic temperature significantly leads by 406 ± 200 years at the onset of the Holocene period. Our results further support the hypothesis of no convective zone at EPICA Dome C during the last deglaciation and the use of nitrogen-15 to infer the height of the diffusive zone. Future climate and carbon cycle modeling works should take into account this robust phasing constraint.


2012 ◽  
Vol 8 (2) ◽  
pp. 1089-1131 ◽  
Author(s):  
F. Parrenin ◽  
S. Barker ◽  
T. Blunier ◽  
J. Chappellaz ◽  
J. Jouzel ◽  
...  

Abstract. We compare a variety of methods for estimating the gas/ice depth offset (Δdepth) at EPICA Dome C (EDC, East Antarctica). (1) Purely based on modelling efforts, Δdepth can be estimated combining a firn densification with an ice flow model. Observations allow direct and indirect estimate of Δdepth. (2) The diffusive column height can be estimated from δ15N and converted to Δdepth using an ice flow model and assumptions about past average firn density and thickness of the convective zone. (3) Ice and gas synchronisation of the EDC ice core to the GRIP, EDML and TALDICE ice cores shifts the ice/gas offset problem into higher accumulation ice cores where it can be more accurately evaluated. (4) Finally, the bipolar seesaw hypothesis allows us to synchronise the ice isotopic record with the gas CH4 record, the later being taken as a proxy of Greenland temperature. The bipolar seesaw antiphase relationship is generally supported by the ice-gas cross synchronisation between EDC and the GRIP, EDML and TALDICE ice cores, which provide support for method 4. Applying the bipolar seesaw hypothesis to the deeper section of the EDC core confirms that the ice flow is complex and can help improving our reconstruction of the thinning function and thus of the EDC age scale. We confirm that method 1 overestimates the glacial Δdepth at EDC and we suggested that it is due to an overestimation of the glacial Close Off Depth by the firn densification model. In contrast we find that the glaciological models probably underestimate the Δdepth during termination II. Finally, we show that method 2 based on 15N data produces for the last deglaciation a Δdepth estimate which is in good agreement with methods 3 and 4.


2004 ◽  
Vol 39 ◽  
pp. 339-345 ◽  
Author(s):  
Marzena Kaczmarska ◽  
Elisabeth Isaksson ◽  
Lars Karlöf ◽  
Jan-Gunnar Winther ◽  
Jack Kohler ◽  
...  

AbstractA 100 m long ice core was retrieved from the coastal area of Dronning Maud Land (DML), Antarctica, in the 2000/01 austral summer. The core was dated to AD 1737 by identification of volcanic horizons in dielectrical profiling and electrical conductivity measurement records in combination with seasonal layer counting from high-resolution oxygen isotope (δ18O) data. A mean long-term accumulation rate of 0.29 ma–1w.e. was derived from the high-resolution δ18O record as well as accumulation rates during periods in between the identified volcanic horizons. A statistically significant decrease in accumulation was found from about 1920 to the present. A comparison with other coastal ice cores from DML suggests that this is a regional pattern.


2017 ◽  
Vol 13 (7) ◽  
pp. 833-853 ◽  
Author(s):  
Camille Bréant ◽  
Patricia Martinerie ◽  
Anaïs Orsi ◽  
Laurent Arnaud ◽  
Amaëlle Landais

Abstract. The transformation of snow into ice is a complex phenomenon that is difficult to model. Depending on surface temperature and accumulation rate, it may take several decades to millennia for air to be entrapped in ice. The air is thus always younger than the surrounding ice. The resulting gas–ice age difference is essential to documenting the phasing between CO2 and temperature changes, especially during deglaciations. The air trapping depth can be inferred in the past using a firn densification model, or using δ15N of air measured in ice cores. All firn densification models applied to deglaciations show a large disagreement with δ15N measurements at several sites in East Antarctica, predicting larger firn thickness during the Last Glacial Maximum, whereas δ15N suggests a reduced firn thickness compared to the Holocene. Here we present modifications of the LGGE firn densification model, which significantly reduce the model–data mismatch for the gas trapping depth evolution over the last deglaciation at the coldest sites in East Antarctica (Vostok, Dome C), while preserving the good agreement between measured and modelled modern firn density profiles. In particular, we introduce a dependency of the creep factor on temperature and impurities in the firn densification rate calculation. The temperature influence intends to reflect the dominance of different mechanisms for firn compaction at different temperatures. We show that both the new temperature parameterization and the influence of impurities contribute to the increased agreement between modelled and measured δ15N evolution during the last deglaciation at sites with low temperature and low accumulation rate, such as Dome C or Vostok. We find that a very low sensitivity of the densification rate to temperature has to be used in the coldest conditions. The inclusion of impurity effects improves the agreement between modelled and measured δ15N at cold East Antarctic sites during the last deglaciation, but deteriorates the agreement between modelled and measured δ15N evolution at Greenland and Antarctic sites with high accumulation unless threshold effects are taken into account. We thus do not provide a definite solution to the firnification at very cold Antarctic sites but propose potential pathways for future studies.


2020 ◽  
Vol 20 (9) ◽  
pp. 5861-5885 ◽  
Author(s):  
V. Holly L. Winton ◽  
Alison Ming ◽  
Nicolas Caillon ◽  
Lisa Hauge ◽  
Anna E. Jones ◽  
...  

Abstract. The nitrogen stable isotopic composition in nitrate (δ15N-NO3-) measured in ice cores from low-snow-accumulation regions in East Antarctica has the potential to provide constraints on past ultraviolet (UV) radiation and thereby total column ozone (TCO) due to the sensitivity of nitrate (NO3-) photolysis to UV radiation. However, understanding the transfer of reactive nitrogen at the air–snow interface in polar regions is paramount for the interpretation of ice core records of δ15N-NO3- and NO3- mass concentrations. As NO3- undergoes a number of post-depositional processes before it is archived in ice cores, site-specific observations of δ15N-NO3- and air–snow transfer modelling are necessary to understand and quantify the complex photochemical processes at play. As part of the Isotopic Constraints on Past Ozone Layer Thickness in Polar Ice (ISOL-ICE) project, we report new measurements of NO3- mass concentration and δ15N-NO3- in the atmosphere, skin layer (operationally defined as the top 5 mm of the snowpack), and snow pit depth profiles at Kohnen Station, Dronning Maud Land (DML), Antarctica. We compare the results to previous studies and new data, presented here, from Dome C on the East Antarctic Plateau. Additionally, we apply the conceptual 1D model of TRansfer of Atmospheric Nitrate Stable Isotopes To the Snow (TRANSITS) to assess the impact of NO3- recycling on δ15N-NO3- and NO3- mass concentrations archived in snow and firn. We find clear evidence of NO3- photolysis at DML and confirmation of previous theoretical, field, and laboratory studies that UV photolysis is driving NO3- recycling and redistribution at DML. Firstly, strong denitrification of the snowpack is observed through the δ15N-NO3- signature, which evolves from the enriched snowpack (−3 ‰ to 100 ‰), to the skin layer (−20 ‰ to 3 ‰), to the depleted atmosphere (−50 ‰ to −20 ‰), corresponding to mass loss of NO3- from the snowpack. Based on the TRANSITS model, we find that NO3- is recycled two times, on average, before it is archived in the snowpack below 15 cm and within 0.75 years (i.e. below the photic zone). Mean annual archived δ15N-NO3- and NO3- mass concentration values are 50 ‰ and 60 ng g−1, respectively, at the DML site. We report an e-folding depth (light attenuation) of 2–5 cm for the DML site, which is considerably lower than Dome C. A reduced photolytic loss of NO3- at DML results in less enrichment of δ15N-NO3- than at Dome C mainly due to the shallower e-folding depth but also due to the higher snow accumulation rate based on TRANSITS-modelled sensitivities. Even at a relatively low snow accumulation rate of 6 cm yr−1 (water equivalent; w.e.), the snow accumulation rate at DML is great enough to preserve the seasonal cycle of NO3- mass concentration and δ15N-NO3-, in contrast to Dome C where the depth profiles are smoothed due to longer exposure of surface snow layers to incoming UV radiation before burial. TRANSITS sensitivity analysis of δ15N-NO3- at DML highlights that the dominant factors controlling the archived δ15N-NO3- signature are the e-folding depth and snow accumulation rate, with a smaller role from changes in the snowfall timing and TCO. Mean TRANSITS model sensitivities of archived δ15N-NO3- at the DML site are 100 ‰ for an e-folding depth change of 8 cm, 110 ‰ for an annual snow accumulation rate change of 8.5 cm yr−1 w.e., 10 ‰ for a change in the dominant snow deposition season between winter and summer, and 10 ‰ for a TCO change of 100 DU (Dobson units). Here we set the framework for the interpretation of a 1000-year ice core record of δ15N-NO3- from DML. Ice core δ15N-NO3- records at DML will be less sensitive to changes in UV than at Dome C; however the higher snow accumulation rate and more accurate dating at DML allows for higher-resolution δ15N-NO3- records.


Sign in / Sign up

Export Citation Format

Share Document