Interpretation of ice core trace gas records depends on an accurate understanding of the processes that smooth the atmospheric signal in the firn. Much work has been done to understand the processes affecting air transport in the open pores of the firn, but a paucity of data from air trapped in bubbles in the firn-ice transition region has limited the ability to constrain the effect of bubble closure processes. Here we present high-resolution measurements of firn density, methane concentrations, nitrogen isotopes, and total air content that show layering in the firn-ice transition region at the WAIS Divide ice core site. Using the notion that bubble trapping is a stochastic process, we derive a new parameterization for closed porosity that incorporates the effects of layering in a steady-state firn modeling approach. We include the process of bubble trapping into an open-porosity firn air transport model, and obtain a good fit to the firn core data. We find that layering broadens the depth range over which bubbles are trapped, widens the gas-age distribution of air in closed bubbles, reduces the mean gas-age of air in closed bubbles, and introduces stratigraphic irregularities in the gas-age scale that have a peak-to-peak variability of ~10 years at WAIS-Divide. For a more complete understanding of gas occlusion and its impact on ice core records we suggest that this experiment be repeated at sites climatically different from WAIS Divide, for example on the East Antarctic plateau.