Effective coefficient of diffusion and permeability of firn at Dome C and Lock In, Antarctica, and of various snow types - estimates over the 100-850kgm-3 density range

Abstract

Modeling air transport through the entire firn column of polar ice sheets is needed to interpret climate archives. To this end, different regressions have been proposed in the past to estimate the effective coefficient of diffusion and permeability of firn. These regressions are often valid for specific depth or porosity ranges only. Also, they constitute a source of uncertainty as evaluations have been limited by the lack of reliable data of firn transport properties. To contribute with a new dataset, this study presents the effective coefficient of diffusion and the permeability at Dome C and Lock In, Antarctica, from the near-surface to the close-off (23 to 133m depth). Also, microstructure is characterized based on density, specific surface area, closed porosity ratio, connectivity index, and structural anisotropy through the correlation lengths. All properties were estimated based on pore-scale computations from 3D tomographic images of firn samples. The normalized diffusion coefficient ranges from 1.9×10-1 to 8.3×10-5, and permeability ranges from 1.2×10-9 to 1.1 ×10-12m2, for densities between 565 and 888kgm-3. No or little anisotropy is reported. Next, we investigate the relationship of the transport properties with density over the firn density range (550-850kgm-3), as well as over the entire density range encountered in the ice sheets (100-850kgm-3), by extending the datasets with transport properties of alpine and artificial snow from previous studies. Classical analytical models and regressions from literature are evaluated against the estimates from pore-scale simulations. For firn, good agreements are found for permeability and the diffusion coefficient with two existing regressions of the literature based on open porosity despite the rather different site conditions (Greenland). Over the entire 100-850kgm-3 density range, permeability is accurately reproduced by the Carman-Kozeny and self-consistent (spherical bi-composite) models when expressed in terms of a rescaled porosity, φresCombining double low line(φ-φoff)/(1-φoff), to account for pore closure, where φoff is the close-off porosity. For the normalized diffusion coefficient, none of the evaluated formulas were satisfactory, so we propose a new regression based on the rescaled porosity that reads D/DairCombining double low line(φres)1.61.