1UJF-Grenoble 1/CNRS, Grenoble Image Parole Signal Automatique (GIPSA-lab), UMR5216, B.P. 46, 38402 St Martin d'Hères, France
2UJF-Grenoble 1/CNRS, Laboratoire de Glaciologie et Géophysique de l'Environnement (LGGE) UMR5183, Grenoble, 38041, France
3School of Environmental Sciences, University of East Anglia, Norwich NR4 7TJ, UK
4National Institute of Polar Research, 10-3 Midorichou, Tachikawa, Tokyo 190-8518, Japan
5Institute of Biogeosciences, Japan Agency for Marine-Earth Science and Technology, 2-15 Natsushima-cho, Yokosuka 237-0061, Japan
6Laboratoire des Sciences du Climat et de L'Environnement, IPSL/CEA-CNRS-UVSQ, 91191 Gif-sur-Yvette, France
7British Antarctic Survey, Natural Environment Research Council, High Cross, Madingley Road, Cambridge CB3 0ET, UK
8NOAA Earth System Research Laboratory, Boulder, Colorado, USA
9Commonwealth Scientific and Industrial Research Organisation, Marine and Atmospheric Research, PMB 1, Aspendale, Vic. 3195, Australia
10Centre for Ice and Climate, Niels Bohr Institute, University of Copenhagen, Juliane Maries, VEJ 30, 2100 Copenhagen Ø, Denmark
Received: 17 Jul 2011 – Published in Atmos. Chem. Phys. Discuss.: 16 Aug 2011
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.
Revised: 04 Oct 2012 – Accepted: 20 Oct 2012 – Published: 04 Dec 2012
Citation: Witrant, E., Martinerie, P., Hogan, C., Laube, J. C., Kawamura, K., Capron, E., Montzka, S. A., Dlugokencky, E. J., Etheridge, D., Blunier, T., and Sturges, W. T.: A new multi-gas constrained model of trace gas non-homogeneous transport in firn: evaluation and behaviour at eleven polar sites, Atmos. Chem. Phys., 12, 11465-11483, doi:10.5194/acp-12-11465-2012, 2012.