UC San Diego

Abstract. Ancient air trapped in ice core bubbles has been paramount to developing our understanding of past climate and atmospheric composition. Before air bubbles become isolated in ice, the atmospheric signal is altered in the firn column by transport processes such as advection and diffusion. However, the influence of impermeable layers and barometric pumping (driven by surface pressure variability) on firn air transport is not well understood and cannot be captured in conventional 1-dimensional firn air models. Here we present a 2-dimensional (2D) trace gas advection-diffusion-dispersion model that accounts for discontinuous horizontal layers of reduced permeability. We find that layering and barometric pumping individually yield too small a reduction in gravitational settling to match observations. In contrast, a combination of both effects more strongly suppresses gravitational fractionation. Layering locally focuses airflows in the 2D model and thus amplifies the dispersive mixing resulting from barometric pumping. Hence, the representation of both factors is needed to obtain a more natural emergence of the lock-in zone. Moreover, we find that barometric pumping in the layered 2D model does not substantially change the differential kinetic fractionation of fast and slow diffusing trace gases, which is observed in nature. This suggests that further subgrid-scale physics may be missing in the current generation of firn air transport models. However, we find robust scaling relationships between kinetic isotope fractionation of different noble gas isotope and elemental ratios. These relationships may be used to correct for kinetic fractionation in future high precision ice core studies.