Biogeochemical reactions and vadose zone transport, in particular gas phase transport, are inherently coupled processes. To explore feedback mechanisms between these processes in a quantitative manner, multicomponent gas diffusion and advection are implemented into an existing reactive transport model that includes a full suite of geochemical reactions. Multicomponent gas diffusion is described based on the Dusty Gas Model, which provides the most generally applicable description for gas diffusion. Gas advection is described by Darcy's Law, which in the current formulation, is directly substituted into the transport equations. The model is used to investigate the interactions between geochemical reactions and transport processes with an emphasis to quantify reaction-induced gas migration in the vadose zone. Simulations of pyrite oxidation in mine tailings, gas attenuation in partially saturated landfill soil covers, and methane production and oxidation in aquifers contaminated by organic compounds demonstrate how biogeochemical reactions drive diffusive and advective transport of reactive and non-reactive gases. Pyrite oxidation in mine tailings causes a pressure reduction in the reaction zone and drives advective gas flow into the sediment column, enhancing the oxidation process. Release of carbon dioxide by carbonate mineral dissolution partly offsets pressure reduction, and illustrates the role of water-rock interaction on gas transport. Microbially mediated methane oxidation in landfill covers reduces the existing upward pressure gradient, thereby decreasing the contribution of advective methane emissions to the atmosphere and enhancing the net flux of atmospheric oxygen into the soil column. At an oil spill site, both generation of CH4 in the methanogenic zone and oxidation of CH4 in the methanotrophic zone contribute to drive advective and diffusive fluxes. The model confirmed that non-reactive gases tend to accumulate in zones of gas consumption and become depleted in zones of gas production. In most cases, the model was able to quantify existing conceptual models, but also proved useful to identify data gaps, sensitivity, and inconsistencies in conceptual models. The formulation of the model is general and can be applied to other vadose zone systems in which reaction-induced gas transport is of importance.