Close-off Depth Research Articles

Horizontal ice flow impacts the firn structure of Greenland's percolation zone

Abstract. One-dimensional simulations of firn evolution neglect horizontal advection from ice flow, which transports the firn column across climate gradients as it is buried by accumulation. Using a suite of model runs, we demonstrate the impacts of horizontal advection on the development of firn density, temperature, and the stratigraphy of melt features through the Greenland ice sheet percolation zone. The simulations isolate processes in synthetic runs and investigate four specific transects and an ice core site. Relative to one-dimensional simulations, the horizontal advection process tends to increase the pore close-off depth, reduce the heat content, and decrease the frequency of melt features with depth by emplacing firn sourced from higher locations under increasingly warm and melt-affected surface conditions. Preservation of the advected pore space and cold content is strongly dependent upon the depth of meltwater infiltration. Horizontal ice flow interacts with topography, climate gradients, and meltwater infiltration to influence the evolution of the firn column structure; the interaction between these variables modulates the impact of horizontal advection on firn at locations around Greenland. Pore close-off and firn temperature are mainly impacted in the lowermost 20–30 km of the percolation zone, which may be relevant to migration of the lower percolation zone. Relatively high in the percolation zone, however, the stratigraphy of melt features can have an advection-derived component that should not be conflated with changing climate.

Modeling Dry-Snow Densification without Abrupt Transition

An empirical model for the densification of dry snow has been calibrated using strain-rate data from Pine Island Glacier basin, Antarctica. The model provides for a smooth transition between Stage 1 and Stage 2 densification, and leads to an analytical expression for density as a function of depth. It introduces two new parameters with a simple physical basis: transition density ρ T and a scaling factor, M, which controls the extent of the transition zone. The standard (Herron and Langway) parameterization is used for strain rates away from the transition zone. Calibration, though tentative, produces best parameter values of ρ T = 580 kg m − 3 and M = 7 for the region. Using these values, the transition model produces better simulations of snow profiles from Pine Island Glacier basin than the well-established Herron and Langway and Ligtenberg models, both of which postulate abrupt transition. Simulation of density profiles from other sites using M = 7 produces the best values of ρ T = 550 kg m − 3 for a high accumulation site and 530 kg m − 3 for a low accumulation site, suggesting that transition density may vary with climatic conditions. The variation of bubble close-off depth and depth-integrated porosity with mean annual accumulation predicted by the transition model is similar to that predicted by the Simonsen model tuned for Greenland.

Pore morphology of polar firn around closure revealed by X-ray tomography

Abstract. Understanding the slow densification process of polar firn into ice is essential in order to constrain the age difference between the ice matrix and entrapped gases. The progressive microstructure evolution of the firn column with depth leads to pore closure and gas entrapment. Air transport models in the firn usually include a closed porosity profile based on available data. Pycnometry or melting–refreezing techniques have been used to obtain the ratio of closed to total porosity and air content in closed pores, respectively. X-ray-computed tomography is complementary to these methods, as it enables one to obtain the full pore network in 3-D. This study takes advantage of this nondestructive technique to discuss the morphological evolution of pores on four different Antarctic sites. The computation of refined geometrical parameters for the very cold polar sites Dome C and Lock In (the two Antarctic plateau sites studied here) provides new information that could be used in further studies. The comparison of these two sites shows a more tortuous pore network at Lock In than at Dome C, which should result in older gas ages in deep firn at Lock In. A comprehensive estimation of the different errors related to X-ray tomography and to the sample variability has been performed. The procedure described here may be used as a guideline for further experimental characterization of firn samples. We show that the closed-to-total porosity ratio, which is classically used for the detection of pore closure, is strongly affected by the sample size, the image reconstruction, and spatial heterogeneities. In this work, we introduce an alternative parameter, the connectivity index, which is practically independent of sample size and image acquisition conditions, and that accurately predicts the close-off depth and density. Its strength also lies in its simple computation, without any assumption of the pore status (open or close). The close-off prediction is obtained for Dome C and Lock In, without any further numerical simulations on images (e.g., by permeability or diffusivity calculations).

Post-bubble close-off fractionation of gases in polar firn and ice cores: effects of accumulation rate on permeation through overloading pressure

Abstract. Gases in ice cores are invaluable archives of past environmental changes (e.g., the past atmosphere). However, gas fractionation processes after bubble closure in the firn are poorly understood, although increasing evidence indicates preferential leakages of smaller molecules (e.g., neon, oxygen, and argon) from the closed bubbles through the ice matrix. These fractionation processes are believed to be responsible for the observed millennial δO2/N2 variations in ice cores, linking ice core chronologies with orbital parameters. In this study, we investigated high-resolution δAr/N2 of the GISP2 (Greenland Ice Sheet Project 2), NGRIP (North Greenland Ice Core Project), and Dome Fuji ice cores for the past few thousand years. We find that δAr/N2 at multidecadal resolution on the "gas-age scale" in the GISP2 ice core has a significant negative correlation with accumulation rate and a positive correlation with air contents over the past 6000 years, indicating that changes in overloading pressure induced δAr/N2 fractionation in the firn. Furthermore, the GISP2 temperature and accumulation rate for the last 4000 years have nearly equal effects on δAr/N2 with sensitivities of 0.72 ± 0.1 ‰ °C−1 and −0.58 ± 0.09 ‰ (0.01 m ice year−1)−1, respectively. To understand the fractionation processes, we applied a permeation model for two different processes of bubble pressure build-up in the firn, "pressure sensitive process" (e.g., microbubbles: 0.3–3 % of air contents) with a greater sensitivity to overloading pressures and "normal bubble process". The model indicates that δAr/N2 in the bubbles under the pressure sensitive process are negatively correlated with the accumulation rate due to changes in overloading pressure. On the other hand, the normal bubbles experience only limited depletion (< 0.5 ‰) in the firn. Colder temperatures in the firn induce more depletion in δAr/N2 through thicker firn. The pressure sensitive bubbles are so depleted in δAr/N2 at the bubble close-off depth that they dominate the total δAr/N2 changes in spite of their smaller air contents. The model also indicates that δAr/N2 of ice cores should have experienced several per mil of depletion during the storage 14–18 years after coring. Further understanding of the δAr/N2 fractionation processes in the firn, combined with nitrogen and argon isotope data, may lead to a new proxy for the past temperature and accumulation rate.

Steady-state size distribution of air bubbles in polar ice

Between the close-off depth and the bubble-to-hydrate transition zone in polar ice sheets, the geometrical properties of air bubbles, such as number concentration and size of bubbles, are mainly controlled by firn temperature and ice accumulation rate prevailing during the snow to ice transformation [3], and by the bubble compression in the course of bubbly ice densification [13]. This implies that the data on the bubble properties can be used for reconstruction of the past climate change. On the basis of our earlier studies of bubbly ice densification and the new measurements of air bubbles in the Antarctic ice cores, we have developed a theory of bubble evolution in polar ice and propose an inverse procedure for bubble size conversion to specified conditions at the close-off depth. Both outcomes of the research contribute to elaboration of the new paleoclimatological tool based on the bubble properties.

Assessment of diffusive isotopic fractionation in polar firn, and application to ice core trace gas records

During rapid variations of the atmospheric mixing ratio of a trace gas, diffusive transport in the porous firn layer atop ice sheets and glaciers alters the isotopic composition of that gas relative to the overlying atmosphere. Records of past atmospheric trace gas isotopic composition from ice cores and firn need to be corrected for this diffusive fractionation artifact.We present a novel, semi-empirical method to accurately estimate the magnitude of the diffusive fractionation in the ice core record. Our method (1) consists of a relatively simple analytical calculation; (2) requires only commonly available ice core data; (3) is not subject to the uncertainties inherent to estimating the accumulation rate, temperature, close-off depth and depth–diffusivity relationship back in time; (4) does not require knowledge of the true atmospheric variations, but uses the smoothed records obtained from ice cores; (5) arguably gives more accurate results than a combined firn densification–firn air transport modeling study would.We apply the method to records of CH4, CO2 and N2O mixing ratios, and we find that the correction is particularly important for δ13C–CH4. We apply the correction to δ13C–CH4 records over the last glacial termination and the 8.2ka event. In both cases the diffusive signal exceeds the analytical precision of the data, and has a significant impact on the observed isotopic trends. For the 8.2ka event the corrected data show an isotopic enrichment in δ13C–CH4 for the duration of the event, consistent with reduced wetland emissions.

Dating a 109.9 m ice core from Dome A (East Antarctica) with volcanic records and a firn densification model

A 109.9 m ice core was extracted at a location about 300 m away from the Dome A summit (80°00′S, 77°21″E) by the Chinese team of the International Trans-Antarctic Science Expedition (ITASE) during the 21st Chinese National Antarctica Research Expedition (CHINARE) in January 2005. Two independent methods were used for dating the ice core, volcanic event markers shown by prominent non-sea-salt sulfate (nss-SO4 2−) and the Herron and Langway (H-L) firn densification model. Six prominent volcanic events (Agung 1963 AD, Tambora 1815 AD, Kuwae 1453 AD, Unknown 1259 AD, Taupo 186 AD and Pinatubo 1050 BC) were identified by comparison with other Antarctic ice cores. Based on the mean accumulation rates between adjacent events, we estimate the age at the firn pore close-off depth (102 m) was 3516±100 a BP. This is the oldest close-off age ever reported from the Antarctic and the Greenland ice sheets. Calculations using the H-L model show that the age at the same depth is 3581±100 a BP. The two dating techniques differ by 65 years, or ∼1.8% of the record. We calculated the bottom age of the ice core as 4009±150 a BP using the volcanic dating method and 4115±150 a BP using the H-L model method.

Preliminary results of the close-off depth and the stable isotopic records along a 109.91 m ice core from Dome A, Antarctica

A 109.91 m ice core was recovered from Dome A (or Dome Argus), the highest ice feature in Antarctica, during the 2004/05 austral summer by the 21st Chinese National Antarctic Research Expedition (CHINARE-21). Both methane profile along the core and firn densification model calculation suggest that the close-off depth is at about 102.0 m with an ice age about 4200 a. Stable isotopes (δ18O and δD) of the chips samples produced during each run of ice core drilling at Dome A, together with those of the other cores recovered from the eastern inland Antarctica, suggest a relative stable climate with a temperature fluctuation amplitude about ±0.6°C at the eastern inland Antarctica during the late Holocene. The average d-excess (or d =δ D-8δ18O) of 17.1‰ along the Dome A core is probably the highest among the Antarctic inland ice cores, which may be resulted from the kinetic fractionation during the snow formation under an oversaturation condition. Moreover, the increasing trend of d-excess during the late Holocene reflects mainly the migration of the water source area for precipitation at Dome A towards low latitudes. This paper presents the first results of a shallow ice core recovered from the unexplored highest area of the Antarctic ice sheet, providing a background for the proposing deep ice core drilling at Dome A.

Properties of a marine ice layer under the Amery Ice Shelf, East Antarctica

AbstractThe Amery Ice Shelf, East Antarctica, undergoes high basal melt rates near the southern limit of its grounding line where 80% of the ice melts within 240 km of becoming afloat. A considerable portion of this later refreezes downstream as marine ice. This produces a marine ice layer up to 200 m thick in the northwest sector of the ice shelf concentrated in a pair of longitudinal bands that extend some 200 km all the way to the calving front. We drilled through the eastern marine ice band at two locations 70 km apart on the same flowline. We determine an average accretion rate of marine ice of 1.1 ± 0.2 m a−1, at a reference density of 920 kg m−3 between borehole sites, and infer a similar average rate of 1.3 ± 0.2 m a−1 upstream. The deeper marine ice was permeable enough that a hydraulic connection was made whilst the drill was still 70–100 m above the ice-shelf base. Below this marine close-off depth, borehole video imagery showed permeable ice with water-filled cavities and individual ice platelets fused together, while the upper marine ice was impermeable with small brine-cell inclusions. We infer that the uppermost portion of the permeable ice becomes impermeable with the passage of time and as more marine ice is accreted on the base of the shelf. We estimate an average closure rate of 0.3 m a−1 between the borehole sites; upstream the average closure rate is faster at 0.9 m a−1. We estimate an average porosity of the total marine ice layer of 14–20%, such that the deeper ice must have even higher values. High permeability implies that sea water can move relatively freely through the material, and we propose that where such marine ice exists this renders deep parts of the ice shelf particularly vulnerable to changes in ocean properties.

Simple relations for the close-off depth and age in dry-snow densification

AbstractA physical model for the snow/firn densification process (Salamatin and others, 2006) and Martinerie and others’ (1992, 1994) correlation for the firn density at the pore closure are employed to perform a scale analysis and computational experiments in order to deduce simplified relations for the close-off depth and ice age in quasi-stationary ice formation conditions. The critical snow density at which ice-grain rearrangement stops is used to take into account variability of snow structures subjected to densification. The results obtained are validated on a representative set of ice-core data from 22 sites which covers wide ranges of present-day temperatures and ice accumulation rates. A simple analytical approximation for the density–depth profile is proposed.

Evidence for molecular size dependent gas fractionation in firn air derived from noble gases, oxygen, and nitrogen measurements

We present elemental and isotopic measurements of noble gases (He, Ne, Ar, Kr, and Xe), oxygen and nitrogen of firn air from two sites. The first set of samples was taken in 1998 at the summit of the Devon Ice Cap in the eastern part of Devon Island. The second set was taken in 2001 at NGRIP location (North Greenland). He and Ne are heavily enriched relative to Ar with respect to the atmosphere in the air near the close-off depth at around 50–70 m. The enrichment increases with depth and reaches the maximum value in the deepest samples just above the zone of impermeable ice where no free air could be extracted anymore. Similarly, elemental ratios of O 2 / N 2, O 2 / Ar and Ar / N 2 are increasing with depth. In contrast but in line with expectations, isotopic ratios of 15N / 14N, 18O / 16O, and 36Ar / 40Ar show no significant enrichment near the close-off depth. The observed isotopic ratios in the firn air column can be explained within the uncertainty ranges by the well-known processes of gravitational enrichment and thermal diffusion. To explain the elemental ratios, however, an additional fractionation process during bubble inclusion has to be considered. We implemented this additional process into our firn air model. The fractionation factors were found by fitting model profiles to the data. We found a very similar close-off fractionation behavior for the different molecules at both sites. For smaller gas species (mainly He and Ne) the fractionation factors are linearly correlated to the molecule size, whereas for diameters greater than about 3.6 Å the fractionation seems to be significantly smaller or even negligible. An explanation for this size dependent fractionation process could be gas diffusion through the ice lattice. At Devon Island the enrichment at the bottom of the firn air column is about four times higher compared to NGRIP. We explain this by lower firn diffusivity at Devon Island, most probably due to melt layers, resulting in significantly reduced back diffusion of the excess gas near the close-off depth. The results of this study considerably increase the understanding of the processes occurring during air bubble inclusion near the close-off depth in firn and can help to improve the interpretation of direct firn air measurements, as well as air bubble measurements in ice cores, which are used in numerous studies as paleo proxies.

Model calculations of the age of firn air across the Antarctic continent

Abstract. The age of firn air in Antarctica at pore close-off depth is only known for a few specific sites where firn air has been sampled for analyses. We present a model that calculates the age of firn air at pore close-off depth for the entire Antarctic continent. The model basically uses four meteorological parameters as input (surface temperature, pressure, accumulation rate and wind speed). Using parameterisations for surface snow density, pore close-off density and tortuosity, in combination with a density-depth model and data of a regional atmospheric climate model, distribution of pore close-off depth for the entire Antarctic continent is determined. The deepest pore close-off depth was found for the East Antarctic Plateau near 72° E, 82° S, at 150±15 m (2σ). A firn air diffusion model was applied to calculate the age of CO2 at pore close-off depth. The results predict that the oldest firn gas (CO2 age) is located between Dome Fuji, Dome Argos and Vostok at 43° E, 78° S being 148±23 (1σ or 38 for 2σ) years old. At this location an atmospheric trace gas record should be obtained. In this study we show that methyl chloride could be recorded with a predicted length of 125 years as an example for trace gas records at this location. The longest record currently available from firn air is derived at South Pole, being 80 years. Sensitivity tests reveal that the locations with old firn air (East Antarctic Plateau) have an estimated uncertainty (2σ) for the modelled CO2 age at pore close-off depth of 30% and of about 40% for locations with younger firn air (CO2 age typically 40 years). Comparing the modelled age of CO2 at pore close-off depth with directly determined ages at seven sites yielded a correlation coefficient of 0.90 and a slope close to 1, suggesting a high level of confidence for the modelled results in spite of considerable remaining uncertainties.

Methanesulphonic acid movement in solid ice cores

AbstractMethanesulphonic acid (MSA) is an important trace-ion constituent in ice cores, with connections to biological activity and sea-ice distribution. Post-depositional movement of MSA has been documented in firn, and this study investigates movement in solid ice by measuring variations in MSA distribution across several horizontal sections from an ice core after 14.5 years storage. The core used is from below the bubble close-off depth at Dome Summit South, Law Dome, East Antarctica. MSA concentration was studied at 3 and 0.5 cm resolution across the core widths. Its distribution was uniform through the core centres, but the outer 3 cm showed gradients in concentrations down to less than half of the central value at the core edge. This effect is consistent with diffusion to the surrounding air during its 14.5 year storage. The diffusion coefficient is calculated to be 2 ×10–13 m2 s–1, and the implications for the diffusion mechanism are discussed.