scholarly journals Modelling the firn thickness evolution during the last deglaciation: constrains on sensitivity to temperature and impurities

2016 ◽  
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 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 that the surrounding ice. The resulting gas-ice age difference is essential to document 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 in several sites of East Antarctica, predicting larger firn thickness during the Last Glacial Maximum, whereas δ15N suggests a reduced firn thickness compared to the Holocene. We present here modifications of the LGGE firn densification model, which significantly reduce the model-data mismatch for the gas trapping depth evolution over the last deglaciation, while preserving the good agreement between measured and modelled modern firn density profiles. In particular, we introduce a dependency of the activation energy to temperature and impurities in the firn densification rate calculation. The temperature influence reflects the existence 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. However, the inclusion of impurities effects deteriorates the agreement between modelled and measured δ15N evolution in Greenland and Antarctic sites with high accumulation.

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.


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.


1982 ◽  
Vol 3 ◽  
pp. 17-22 ◽  
Author(s):  
J. P. Benoist ◽  
J. Jouzel ◽  
C. Lorius ◽  
L. Merlivat ◽  
M. Pourchet

Data on climatic changes over thousands of years is needed for a better understanding of the shorter term variations which are of interest to man. For this purpose we measured the isotope composition (δD‰) of two adjacent ice cores drilled in the Dome C area. The time scale was established using the remarkably constant mean annual accumulation rate (37 kg m−2) determined by various techniques. The detailed isotope records were smoothed to filter out the δ value fluctuations not directly related to local temperature changes. With respect to conditions over the last 2.5 ka, the combined smoothed δ curve indicates a cooler climate from about 1800 to 1200 AD and a slightly warmer period from about 1200 to 700 AD. These periods may well correspond to the suggested world-wide Little Ice Age and medieval warm phase. Using the present δD‰/T°C measured at the surface, the maximum amplitude for these two periods, after smoothing with a low pass filter of 512 a, is approximately -0.35 and +0.3°C, respectively.


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.


2013 ◽  
Vol 9 (3) ◽  
pp. 2967-3013 ◽  
Author(s):  
S. O. Rasmussen ◽  
P. Abbott ◽  
T. Blunier ◽  
A. Bourne ◽  
E. Brook ◽  
...  

Abstract. A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the Electrical Conductivity Measurement (ECM) and Dielectrical Profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide additional horizons used for the time scale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard–Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age–gas age difference (Δage) has been reconstructed using a coupled firn densification–heat diffusion model. Temperature and accumulation inputs to the Δage model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from δ15N of nitrogen and high-resolution methane data during the abrupt onsets of interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, and based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed, providing the necessary basis for further analysis of the records from NEEM.


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.


2013 ◽  
Vol 9 (6) ◽  
pp. 2713-2730 ◽  
Author(s):  
S. O. Rasmussen ◽  
P. M. Abbott ◽  
T. Blunier ◽  
A. J. Bourne ◽  
E. Brook ◽  
...  

Abstract. A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard–Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age–gas age difference (Δage) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Δage model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from δ15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.


2020 ◽  
Vol 16 (5) ◽  
pp. 1691-1713 ◽  
Author(s):  
James E. Lee ◽  
Edward J. Brook ◽  
Nancy A. N. Bertler ◽  
Christo Buizert ◽  
Troy Baisden ◽  
...  

Abstract. In 2013 an ice core was recovered from Roosevelt Island, an ice dome between two submarine troughs carved by paleo-ice-streams in the Ross Sea, Antarctica. The ice core is part of the Roosevelt Island Climate Evolution (RICE) project and provides new information about the past configuration of the West Antarctic Ice Sheet (WAIS) and its retreat during the last deglaciation. In this work we present the RICE17 chronology, which establishes the depth–age relationship for the top 754 m of the 763 m core. RICE17 is a composite chronology combining annual layer interpretations for 0–343 m (Winstrup et al., 2019) with new estimates for gas and ice ages based on synchronization of CH4 and δ18Oatm records to corresponding records from the WAIS Divide ice core and by modeling of the gas age–ice age difference. Novel aspects of this work include the following: (1) an automated algorithm for multiproxy stratigraphic synchronization of high-resolution gas records; (2) synchronization using centennial-scale variations in methane for pre-anthropogenic time periods (60–720 m, 1971 CE to 30 ka), a strategy applicable for future ice cores; and (3) the observation of a continuous climate record back to ∼65 ka providing evidence that the Roosevelt Island Ice Dome was a constant feature throughout the last glacial period.


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 (2) ◽  
pp. 435-467 ◽  
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 are not well tested on low accumulation and cold sites of the East Antarctic plateau, especially for periods with different climatic conditions. 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 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 and Δdepth during this event. It allows us to evaluate the model and to link together climatic archives from EDC and EDML to NorthGRIP (Greenland). Our results 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. Our finding suggests that the phase relationship between CO2 and EDC temperature inferred at the start of the last deglaciation (lag of CO2 by 800±600 yr) is overestimated and that the CO2 increase could well have been in phase or slightly leading the temperature increase at EDC.


Sign in / Sign up

Export Citation Format

Share Document