scholarly journals Optimisation of glaciological parameters for ice core chronology by implementing counted layers between identified depth levels

2014 ◽  
Vol 10 (4) ◽  
pp. 3585-3616 ◽  
Author(s):  
L. Bazin ◽  
B. Lemieux-Dudon ◽  
A. Landais ◽  
M. Guillevic ◽  
P. Kindler ◽  
...  
Keyword(s):  
Bayesian Approach ◽  
Accumulation Rate ◽  
Ice Core ◽  
Ice Cores ◽  
Age Difference ◽  
Sensitivity Tests ◽  
New Type ◽  
Lock In ◽  
Ice Core Dating

Abstract. A~recent coherent chronology has been built for 4 Antarctic ice cores and the NorthGRIP (NGRIP) Greenland ice core (Antarctic Ice Core Chronology 2012, AICC2012) using a bayesian approach for ice core dating (Datice). When building the AICC2012 chronology, and in order to prevent any confusion with official ice cores chronology, it has been imposed that the AICC2012 chronology for NGRIP should respect exactly the GICC05 chronology based on layer counting. However, such a strong tuning did not satisfy the hypothesis of independence of background parameters and observations for the NGRIP core as required by Datice. We present here the implementation in Datice of a new type of markers that is better suited to constraints deduced from layer counting: the markers of age-difference. Using this type of markers for NGRIP in a 5 cores dating exercise with Datice, we have performed several sensitivity tests and show that the new ice core chronologies obtained with these new markers do not differ by more than 400 years from AICC2012 for Antarctic ice cores and by more than 130 years from GICC05 for NGRIP over the last 60 000 years. With this new parameterization, the accumulation rate and lock-in depth associated with NGRIP are more coherent with independent estimates than those obtained in AICC2012. While these new chronologies should not be used yet as new ice core chronologies, the improved methodology presented here should be considered in the next coherent ice core dating exercise.

scholarly journals Implementation of counted layers for coherent ice core chronology

2015 ◽  
Vol 11 (6) ◽  
pp. 959-978 ◽  
Author(s):  
B. Lemieux-Dudon ◽  
L. Bazin ◽  
A. Landais ◽  
H. Toyé Mahamadou Kele ◽  
M. Guillevic ◽  
...  
Keyword(s):  
Bayesian Approach ◽  
Ice Core ◽  
Ice Cores ◽  
Global Error ◽  
A Posteriori ◽  
Sensitivity Tests ◽  
New Type ◽  
Ice Core Dating

Abstract. A recent coherent chronology has been built for four Antarctic ice cores and the NorthGRIP (NGRIP) Greenland ice core (Antarctic Ice Core Chronology 2012, AICC2012) using a Bayesian approach for ice core dating (Datice). When building the AICC2012 chronology, and in order to prevent any confusion with official ice core chronology, the AICC2012 chronology for NGRIP was forced to fit exactly the GICC05 chronology based on layer counting. However, such a strong tuning did not satisfy the hypothesis of independence of background parameters and observations for the NGRIP core, as required by Datice. We present here the implementation in Datice of a new type of markers that is better suited for constraints deduced from layer counting: the duration constraints. Estimating the global error on chronology due to such markers is not straightforward and implies some assumption on the correlation between individual counting errors for each interval of duration. We validate this new methodological implementation by conducting twin experiments and a posteriori diagnostics on the NGRIP ice core. Several sensitivity tests on marker sampling and correlation between counting errors were performed to provide some guidelines when using such a method for future dating experiments. Finally, using these markers for NGRIP in a five-core dating exercise with Datice leads to new chronologies that do not differ by more than 410 years from AICC2012 for Antarctic ice cores and 150 years from GICC05 for NGRIP over the last 60 000 years.

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 ◽  
...  
Keyword(s):  
Accumulation Rate ◽  
Ice Core ◽  
Ice Cores ◽  
Detailed Comparison ◽  
Last Deglaciation ◽  
Depth Study ◽  
Dronning Maud Land ◽  
The Last Deglaciation ◽  
Correct Estimation ◽  
Lock In

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.

scholarly journals A first chronology for the NEEM ice core

2013 ◽  
Vol 9 (3) ◽  
pp. 2967-3013 ◽  
Author(s):  
S. O. Rasmussen ◽  
P. Abbott ◽  
T. Blunier ◽  
A. Bourne ◽  
E. Brook ◽  
...  
Keyword(s):  
Accumulation Rate ◽  
Ice Core ◽  
Conductivity Measurement ◽  
Ice Cores ◽  
Monte Carlo Analysis ◽  
Heat Diffusion ◽  
Ice Age ◽  
Age Difference ◽  
Electrical Conductivity Measurement ◽  
Volcanic Origin

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.

scholarly journals A new multi-gas constrained model of trace gas non-homogeneous transport in firn: evaluation and behavior at eleven polar sites

2011 ◽  
Vol 11 (8) ◽  
pp. 23029-23080 ◽  
Author(s):  
E. Witrant ◽  
P. Martinerie ◽  
C. Hogan ◽  
J. C. Laube ◽  
K. Kawamura ◽  
...  
Keyword(s):  
Accumulation Rate ◽  
Trace Gases ◽  
Ice Cores ◽  
Gas Diffusion ◽  
Snow Accumulation ◽  
Least Square ◽  
Trace Gas ◽  
Age Difference ◽  
Depth Gradients ◽  
Lock In

Abstract. Insoluble trace gases are trapped in polar ice at the firn-ice transition, at approximately 50 to 100 m below the surface, depending primarily on the site temperature and snow accumulation. Due to the different time scales for snow accumulation versus diffusion of gases through the snowpack, age differences between gases and the ice in which they are "trapped" can be large; e.g. several thousand years in central Antarctica (a low snow accumulation area). Models of trace gas diffusion in polar firn are used to relate firn air and ice core records of trace gases to their atmospheric history. We propose a new diffusion model based on the following contributions. First, the airflow transport model is revised in a poromechanics framework with specific emphasis on the non-homogeneous properties (convective layer, depth-dependent diffusivity and lock-in zone) and an almost-stagnant behavior described by Darcy's law (gravity effect). We then derive a non-linear least square multi-gas optimization scheme to calculate the effective firn diffusivity (automatic diffusivity tuning). The improvements associated with the additional constraints gained by the multi-gas approach are investigated (up to eleven gases for a single site are included in the optimization process). The model is applied to measured data from four Arctic (Devon Island, NEEM, North GRIP, Summit) and seven Antarctic (DE08, Berkner Island, Siple Dome, Dronning Maud Land, South Pole, Dome C, Vostok) sites and the depth-dependent diffusivity profiles are calculated. Among these different sites, a relationship between an increasing thickness of the lock-in zone defined from the isotopic composition of molecular nitrogen in firn air (denoted δ15N) and the snow accumulation rate is obtained, in accordance with observations. It is associated with reduced diffusivity depth-gradients in deep firn, which decreases gas density depth-gradients, at high accumulation rate sites. This has implications for the understanding of δ15N of N2 records in ice cores, in relation with past variations of the snow accumulation rate. Although the extent of layering is clearly a primary control on the thickness of the lock-in zone, our new approach that allows calculation of an estimated lock-in depth may lead to a better constraint on the age difference between the ice and entrapped gases.

scholarly journals A first chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core

2013 ◽  
Vol 9 (6) ◽  
pp. 2713-2730 ◽  
Author(s):  
S. O. Rasmussen ◽  
P. M. Abbott ◽  
T. Blunier ◽  
A. J. Bourne ◽  
E. Brook ◽  
...  
Keyword(s):  
Accumulation Rate ◽  
Ice Core ◽  
Conductivity Measurement ◽  
Ice Cores ◽  
Heat Diffusion ◽  
Ice Age ◽  
Age Difference ◽  
Volcanic Origin ◽  
The North ◽  
North Greenland

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.

scholarly journals IceChrono v1: a probabilistic model to compute a common and optimal chronology for several ice cores

2014 ◽  
Vol 7 (5) ◽  
pp. 6811-6844
Author(s):  
F. Parrenin
Keyword(s):  
Probabilistic Model ◽  
Accumulation Rate ◽  
Ice Core ◽  
Ice Cores ◽  
Air Bubbles ◽  
Ice Flow ◽  
Orbital Parameters ◽  
Information Depth ◽  
Dating Methods ◽  
Lock In

Abstract. Polar ice cores provides exceptional archives of past environmental conditions. Dating ice and air bubbles/hydrates in ice cores is complicated since it involves different dating methods: modeling of the sedimentation process (accumulation of snow at surface, densification of snow into ice with air trapping and ice flow), use of dated horizons by comparison to other well dated targets (other dated paleo-archives or calculated variations of Earth's orbital parameters), use of dated depth intervals, use of Δdepth information (depth shift between synchronous events in the ice matrix and its air/hydrate content), use of stratigraphic links in between ice cores (ice-ice, air-air or mix ice-air links). Here I propose IceChrono v1, a new probabilistic model to combine these different kinds of chronological information to obtain a common and optimized chronology for several ice cores, as well as its confidence interval. It is based on the inversion of three quantities: the surface accumulation rate, the Lock-In Depth (LID) of air bubbles and the vertical thinning function. IceChrono is similar in scope to the Datice model, but has differences on the mathematical, numerical and programming point of views. I apply IceChrono on two dating experiments. The first one is similar to the AICC2012 experiment and I find similar results than Datice within a few centuries, which is a confirmation of both IceChrono and Datice codes. The second experiment involves only the Berkner ice core in Antarctica and I produce the first dating of this ice core. IceChrono v1 is freely available under the GPL v3 open source license.

scholarly journals A new multi-gas constrained model of trace gas non-homogeneous transport in firn: evaluation and behaviour at eleven polar sites

2012 ◽  
Vol 12 (23) ◽  
pp. 11465-11483 ◽  
Author(s):  
E. Witrant ◽  
P. Martinerie ◽  
C. Hogan ◽  
J. C. Laube ◽  
K. Kawamura ◽  
...  
Keyword(s):  
Transport Model ◽  
Accumulation Rate ◽  
Molecular Nitrogen ◽  
Trace Gases ◽  
Ice Cores ◽  
Snow Accumulation ◽  
Least Square ◽  
Trace Gas ◽  
Age Difference ◽  
Lock In

Abstract. Insoluble trace gases are trapped in polar ice at the firn-ice transition, at approximately 50 to 100 m below the surface, depending primarily on the site temperature and snow accumulation. Models of trace gas transport in polar firn are used to relate firn air and ice core records of trace gases to their atmospheric history. We propose a new model based on the following contributions. First, the firn air transport model is revised in a poromechanics framework with emphasis on the non-homogeneous properties and the treatment of gravitational settling. We then derive a nonlinear least square multi-gas optimisation scheme to calculate the effective firn diffusivity (automatic diffusivity tuning). The improvements gained by the multi-gas approach are investigated (up to ten gases for a single site are included in the optimisation process). We apply the model to four Arctic (Devon Island, NEEM, North GRIP, Summit) and seven Antarctic (DE08, Berkner Island, Siple Dome, Dronning Maud Land, South Pole, Dome C, Vostok) sites and calculate their respective depth-dependent diffusivity profiles. Among these different sites, a relationship is inferred between the snow accumulation rate and an increasing thickness of the lock-in zone defined from the isotopic composition of molecular nitrogen in firn air (denoted δ15N). It is associated with a reduced diffusivity value and an increased ratio of advective to diffusive flux in deep firn, which is particularly important at high accumulation rate sites. This has implications for the understanding of δ15N of N2 records in ice cores, in relation with past variations of the snow accumulation rate. As the snow accumulation rate is clearly a primary control on the thickness of the lock-in zone, our new approach that allows for the estimation of the lock-in zone width as a function of accumulation may lead to a better constraint on the age difference between the ice and entrapped gases.

scholarly journals IceChrono1: a probabilistic model to compute a common and optimal chronology for several ice cores

2015 ◽  
Vol 8 (5) ◽  
pp. 1473-1492 ◽  
Author(s):  
F. Parrenin ◽  
L. Bazin ◽  
E. Capron ◽  
A. Landais ◽  
B. Lemieux-Dudon ◽  
...  
Keyword(s):  
Probabilistic Model ◽  
Accumulation Rate ◽  
Ice Core ◽  
Ice Cores ◽  
Complex Problem ◽  
Polar Ice ◽  
Ice Flow ◽  
Marquardt Algorithm ◽  
Lock In

Abstract. Polar ice cores provide exceptional archives of past environmental conditions. The dating of ice cores and the estimation of the age-scale uncertainty are essential to interpret the climate and environmental records that they contain. It is, however, a complex problem which involves different methods. Here, we present IceChrono1, a new probabilistic model integrating various sources of chronological information to produce a common and optimized chronology for several ice cores, as well as its uncertainty. IceChrono1 is based on the inversion of three quantities: the surface accumulation rate, the lock-in depth (LID) of air bubbles and the thinning function. The chronological information integrated into the model are models of the sedimentation process (accumulation of snow, densification of snow into ice and air trapping, ice flow), ice- and air-dated horizons, ice and air depth intervals with known durations, Δdepth observations (depth shift between synchronous events recorded in the ice and in the air) and finally air and ice stratigraphic links in between ice cores. The optimization is formulated as a least squares problem, implying that all densities of probabilities are assumed to be Gaussian. It is numerically solved using the Levenberg–Marquardt algorithm and a numerical evaluation of the model's Jacobian. IceChrono follows an approach similar to that of the Datice model which was recently used to produce the AICC2012 (Antarctic ice core chronology) for four Antarctic ice cores and one Greenland ice core. IceChrono1 provides improvements and simplifications with respect to Datice from the mathematical, numerical and programming point of views. The capabilities of IceChrono1 are demonstrated on a case study similar to the AICC2012 dating experiment. We find results similar to those of Datice, within a few centuries, which is a confirmation of both IceChrono1 and Datice codes. We also test new functionalities with respect to the original version of Datice: observations as ice intervals with known durations, correlated observations, observations as air intervals with known durations and observations as mixed ice–air stratigraphic links. IceChrono1 is freely available under the General Public License v3 open source license.

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

2012 ◽  
Vol 8 (6) ◽  
pp. 6051-6091 ◽  
Author(s):  
E. Capron ◽  
A. Landais ◽  
D. Buiron ◽  
A. Cauquoin ◽  
J. Chappellaz ◽  
...  
Keyword(s):  
Accumulation Rate ◽  
Ice Core ◽  
Ice Cores ◽  
Clear Correlation ◽  
Last Deglaciation ◽  
Main Driver ◽  
Dronning Maud Land ◽  
James Ross Island ◽  
The Last Deglaciation ◽  
Lock In

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.

scholarly journals Impact of physical properties and accumulation rate on pore close-off in layered firn

2013 ◽  
Vol 7 (3) ◽  
pp. 2533-2566
Author(s):  
S. A. Gregory ◽  
M. R. Albert ◽  
I. Baker
Keyword(s):  
Pore Structure ◽  
Open Porosity ◽  
Accumulation Rate ◽  
Ice Core ◽  
Ice Age ◽  
Coarse Grained ◽  
Age Difference ◽  
Dense Layer ◽  
High Accumulation ◽  
Lock In

Abstract. Investigations into the physical characteristics of deep firn near the lock-in zone through pore close-off are needed to improve understanding of ice core records of past atmospheric concentrations. Specifically, the permeability and microstructure profiles of the firn through the diffusive column influence the entrapment of air into bubbles and thus the ice age-gas age difference. The purpose of this study is to examine the nature of pore closure processes at two polar sites with very different local temperatures and accumulation rates. Density, permeability, and microstructure measurements were made on firn cores from WAIS Divide in West Antarctica and Megadunes in East Antarctica. We found that the open pore structure plays a more important role than density in predicting gas transport properties, through the porous firn matrix. For both WAIS Divide and Megadunes, fine grained layers experience close-off shallower in the firn column than do coarse grained layers, regardless of which grain sized layer is the more dense layer at depth. Pore close-off occurs at an open porosity that is accumulation rate dependent. Low accumulation sites, with coarse grains, close-off at lower open porosities (< 10%) than the open porosity (> 10%) of high accumulation sites with finer grains. The depth and length of the lock-in zone is primarily dependent upon accumulation rate and microstructural variability due to differences in grain size and pore structure, rather than the density variability of the layers.

Sign in / Sign up

Export Citation Format

Share Document