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

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.

2012 ◽  
Vol 12 (23) ◽  
pp. 11465-11483 ◽  
Author(s):  
E. Witrant ◽  
P. Martinerie ◽  
C. Hogan ◽  
J. C. Laube ◽  
K. Kawamura ◽  
...  

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.


2014 ◽  
Vol 10 (4) ◽  
pp. 3585-3616 ◽  
Author(s):  
L. Bazin ◽  
B. Lemieux-Dudon ◽  
A. Landais ◽  
M. Guillevic ◽  
P. Kindler ◽  
...  

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.


2021 ◽  
Author(s):  
Yuzhen Yan ◽  
Nicole E. Spaulding ◽  
Michael L. Bender ◽  
Edward J. Brook ◽  
John A. Higgins ◽  
...  

Abstract. The S27 ice core, drilled in the Allan Hills Blue Ice Area of East Antarctica, is located in Southern Victoria Land ~80 km away from the present-day northern edge of the Ross Ice Shelf. Here, we utilize the reconstructed accumulation rate of S27 covering the Last Interglacial (LIG) period between 129 and 116 thousand years before present (ka) to infer moisture transport into the region. The accumulation rate is based on the ice age-gas age differences calculated from the ice chronology, which is constrained by the stable water isotopes of the ice, and an improved gas chronology based on measurements of oxygen isotopes of O2 in the trapped gases. The peak accumulation rate in S27 occurred at 128.2 ka, near the peak LIG warming in Antarctica. Even the most conservative estimate yields a six-fold increase in the accumulation rate in the LIG, whereas other Antarctic ice cores are typically characterized by a glacial-interglacial difference of a factor of two to three. While part of the increase in S27 accumulation rates must originate from changes in the large-scale atmospheric circulation, additional mechanisms are needed to explain the large changes. We hypothesize that the exceptionally high snow accumulation recorded in S27 reflects open-ocean conditions in the Ross Sea, created by reduced sea ice extent and increased polynya size, and perhaps by a southward retreat of the Ross Ice Shelf relative to its present-day position near the onset of LIG. The proposed ice shelf retreat would also be compatible with a sea-level high stand around 129 ka significantly sourced from West Antarctica. The peak in S27 accumulation rates is transient, suggesting that if the Ross Ice Shelf had indeed retreated during the early LIG, it would have re-advanced by 125 ka.


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.


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.


1994 ◽  
Vol 20 ◽  
pp. 121-128 ◽  
Author(s):  
W. Graf ◽  
H. Moser ◽  
O. Reinwarth ◽  
J. Kipfstuhl ◽  
H. Oerter ◽  
...  

The accumulation and distribution of the2H content of near-surface layers in the eastern part of the Ronne Ice Shelf were determined from 16 firn cores drilled to about 10 m depth during the Filchner IIIa and IV campaigns in 1990 and 1992, respectively. The cores were dated stratigraphically by seasonal δ2H variations in the firn. In addition,3H and high-resolution chemical profiles were used to assist in dating. Both the accumulation rate and the stable-isotope content decrease with increasing distance from the ice edge: the δ2H values range from about 195‰ at the ice edge to -25‰ at BAS sites 5 and 6, south of Henry Ice Rise, and the accumulation rates from about 210 to 90 kgm-2a-1. The δ2H values of the near-surface firn and the 10 m firn temperatures (Θ) at individual sites are very well correlated: dδ2H/dΘ = (10.3 ± 0.6)‰K-1; r = 0.97.The δ2H profiles of the two ice cores BI3 and BI5 drilled in 1990 and 1992 to 215 and 320 m depth, respectively, reflect the gradual depletion in2H in the firn upstream of the drill sites. Comparison with the surface data indicates that the ice above 142 m in core BIS and above 137 m in core BI3 was deposited on the ice shelf, whereas the deeper ice, down to 152.8 m depth, most probably originated from the margin of the Antarctic ice sheet.


2017 ◽  
Vol 13 (11) ◽  
pp. 1685-1693 ◽  
Author(s):  
Christoph Florian Schaller ◽  
Johannes Freitag ◽  
Olaf Eisen

Abstract. In order to interpret the paleoclimatic record stored in the air enclosed in polar ice cores, it is crucial to understand the fundamental lock-in process. Within the porous firn, bubbles are sealed continuously until the respective horizontal layer reaches a critical porosity. Present-day firn air models use a postulated temperature dependence of this value as the only parameter to adjust to the surrounding conditions of individual sites. However, no direct measurements of the firn microstructure could confirm these assumptions. Here we show that the critical porosity is a climate-independent constant by providing an extensive data set of micrometer-resolution 3-D X-ray computer tomographic measurements for ice cores representing different extremes of the temperature and accumulation ranges. We demonstrate why indirect measurements suggest a climatic dependence and substantiate our observations by applying percolation theory as a theoretical framework for bubble trapping. The incorporation of our results significantly influences the dating of trace gas records, changing gas-age–ice-age differences by up to more than 1000 years. This may further help resolve inconsistencies, such as differences between East Antarctic δ15N records (as a proxy for firn height) and model results. We expect our findings to be the basis for improved firn air and densification models, leading to lower dating uncertainties. The reduced coupling of proxies and surrounding conditions may allow for more sophisticated reinterpretations of trace gas records in terms of paleoclimatic changes and will benefit the development of new proxies, such as the air content as a marker of local insolation.


2019 ◽  
Vol 65 (1) ◽  
pp. 46-62 ◽  
Author(s):  
A. A. Ekaykin ◽  
D. O. Vladimirova ◽  
N. A. Tebenkova ◽  
E. V. Brovkov ◽  
A. N. Veres ◽  
...  

The knowledge of the spatial distribution of the snow accumulation rate and isotopic composition in different scales, from local to continental, over the Antarctic Ice Sheet is critically important for the interpretation of the paleoclimate data obtained from deep ice cores, for correct assessment of the ice sheet mass balance, etc. With this in mind, we have synthesized geodetic, glaciological and geochemical data collected in the vicinity of central Antarctic Vostok station in 1970–2017 in order to shed light on the processes governing the spatial distribution of snow isotopic composition and accumulation rate in the spatial scale from 100 to 1000 m. First, we have discovered that snow surface height and snow accumulation rate field are strongly affected by the influence of the logistic convoy route annually operating between Russian Antarctic stations Vostok and Progress. This influence is detectable up to 1 km leeward from the route. At the same time the isotopic composition of the upper 10 cm of the snow does not show any anomalies in the vicinity of the route. This is an unexpected result, because large anomalies of the ice sheet surface (e.g., megadunes) are known to affect the snow isotopic composition. Second, in the undisturbed part of the snow surface near Vostok station we have discovered quasi-periodic (with the wavelength of about 400 m) low-amplitude variations of the surface height that are covariant with the corresponding waves in snow accumulation and isotopic composition. We suggest that spatial variability of the snow isotopic composition is due to the different ratio of summer and winter precipitation deposited in different locations, as evident from a strong negative correlation between δD and dxs parameters. The results of this study may explain the nature of the low-frequency noise (with the time-scale from decades to centuries) observed in the climate records obtained from shallow and deep ice cores in central Antarctica.


2015 ◽  
Vol 9 (4) ◽  
pp. 1481-1504 ◽  
Author(s):  
V. Masson-Delmotte ◽  
H. C. Steen-Larsen ◽  
P. Ortega ◽  
D. Swingedouw ◽  
T. Popp ◽  
...  

Abstract. Combined records of snow accumulation rate, δ18O and deuterium excess were produced from several shallow ice cores and snow pits at NEEM (North Greenland Eemian Ice Drilling), covering the period from 1724 to 2007. They are used to investigate recent climate variability and characterise the isotope–temperature relationship. We find that NEEM records are only weakly affected by inter-annual changes in the North Atlantic Oscillation. Decadal δ18O and accumulation variability is related to North Atlantic sea surface temperature and is enhanced at the beginning of the 19th century. No long-term trend is observed in the accumulation record. By contrast, NEEM δ18O shows multidecadal increasing trends in the late 19th century and since the 1980s. The strongest annual positive δ18O values are recorded at NEEM in 1928 and 2010, while maximum accumulation occurs in 1933. The last decade is the most enriched in δ18O (warmest), while the 11-year periods with the strongest depletion (coldest) are depicted at NEEM in 1815–1825 and 1836–1846, which are also the driest 11-year periods. The NEEM accumulation and δ18O records are strongly correlated with outputs from atmospheric models, nudged to atmospheric reanalyses. Best performance is observed for ERA reanalyses. Gridded temperature reconstructions, instrumental data and model outputs at NEEM are used to estimate the multidecadal accumulation–temperature and δ18O–temperature relationships for the strong warming period in 1979–2007. The accumulation sensitivity to temperature is estimated at 11 ± 2 % °C−1 and the δ18O–temperature slope at 1.1 ± 0.2 ‰ °C−1, about twice as large as previously used to estimate last interglacial temperature change from the bottom part of the NEEM deep ice core.


2004 ◽  
Vol 39 ◽  
pp. 214-218 ◽  
Author(s):  
Gordon S. Hamilton

AbstractSnow-accumulation rates are known to be sensitive to local changes in ice-sheet surface slope because of the effect of katabatic winds. These topographic effects can be preserved in ice cores that are collected at non-ice-divide locations. The trajectory of an ice-core site at South Pole is reconstructed using measurements of ice-sheet motion to show that snow was probably deposited at places of different surface slope during the past 1000 years. Recent accumulation rates, derived from shallow firn cores, vary along this trajectory according to surface topography, so that on a relatively steep flank mean annual accumulation is ∼18% smaller than on a nearby topographic depression. These modern accumulation rates are used to reinterpret the cause of accumulation rate variability with time in the long ice-core record as an ice-dynamics effect and not a climate-change signal. The results highlight the importance of conducting ancillary ice-dynamics measurements as part of ice-coring programs so that topographic effects can be deconvolved from potential climate signals.


Sign in / Sign up

Export Citation Format

Share Document