Simulating discontinuous gas flow and methane reactive transport in the saturated zone