AbstractBarometric pumping plays a crucial role in the release of trace gases from fractured porous media to the atmosphere, and it requires a rigorous and complete modeling in order to go beyond the approximate schemes available in the literature. Therefore, a coupled set of convection and convection‐diffusion equations for a slightly compressible fluid in unsteady conditions should be solved. The numerical methodology is presented, and it is applied to conditions close to the ones of the Roselend Natural Laboratory (France). The precision of the code is assessed and the mechanism of barometric pumping is explained. The usual schematization by simple vertical fractures is shown to be only qualitative. Finally, barometric pumping is shown to be efficient in a narrow range of parameter values; its efficiency is a decreasing function of the matrix porosity and of the fracture density.