Numerical modeling of cold magmatic CO 2 flux measurements for the exploration of hidden geothermal systems
The most accepted conceptual model to explain surface degassing of cold magmatic CO 2 in volcanic-geothermal systems involves the presence of a gas reservoir. Here, numerical simulations using the TOUGH2-ECO2N V2.0 package are performed to get quantitative insights into how cold CO 2 soil flux measurements are related to reservoir and fluid properties. While the modeling is based on flux data measured at a specific geothermal site, the Acoculco caldera (Mexico), some general insights have been gained. Both the CO 2 fluxes at the surface and the depth at which CO 2 exsolves are highly sensitive to the dissolved CO 2 content of the deep fluid. If CO 2 mainly exsolves above the reservoir within a fracture zone, the surface CO 2 fluxes are not sensitive to the reservoir size but depend on the CO 2 dissolved content and the rock permeability. For gas exsolution below the top of the reservoir, surface CO 2 fluxes also depend on the gas saturation of the deep fluid as well as the reservoir size. The absence of thermal anomalies at the surface is mainly a consequence of the low enthalpy of CO 2 . The heat carried by CO 2 is efficiently cooled down by heat conduction and to a certain extent by isoenthalpic volume expansion depending on the temperature gradient. Thermal anomalies occur at higher CO 2 fluxes (>37,000 g m -2 d -1 ) when the heat flux of the rising CO 2 is not balanced anymore. Lastly, specific results are obtained for the Acoculco area (reservoir depth, CO 2 dissolved content, and gas saturation state).