| Popis: |
Modeling air transport through the entire ice sheet column is needed to interpret climate archives. To this end, different regressions were proposed to estimate the effective coefficient of diffusion and permeability of firn. Such regressions are often valid for specific depth or porosity ranges and were little evaluated as data of these properties are scarce. 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 133 m 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 on 3D tomographic images of firn samples. 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−12 m2, for densities between 565 and 888 kg m−3. No or little anisotropy is reported. Next, we investigate the relationship of the transport properties with density over the firn density range as well as over the entire density range encountered in ice sheets by including snow data. Classical analytical models and regressions from literature are evaluated. For firn (550–850 kg m−3), good agreements are found for permeability and diffusion coefficient with the regressions based on the open porosity of Freitag et al. (2002) and Adolph and Albert (2014), despite the rather different site conditions (Greenland). Over the entire 100–850 kg m−3 density range, permeability is accurately reproduced by the Carman-Kozeny and Self-Consistent (spherical bi-composite) model when expressed in terms of a rescaled porosity ϕres = (ϕ–ϕoff) / (1–ϕoff) to account for pore closure, with ϕoff 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/Dair = (ϕres)1.61. |
