Simulation of water–rock interaction in the Yellowstone geothermal system using TOUGHREACT

The Yellowstone geothermal system provides an ideal opportunity to test the ability of reactive transport models to simulate the chemical and hydrological effects ofwater–rock interaction. Previous studies of the Yellowstone geothermal system have characterized water–rock interaction through analysis of rocks and fluids obtained from both surface and downhole samples. Fluid chemistry, rock mineralogy, permeability, porosity, and thermal data obtained from the Y-8 borehole in Upper Geyser Basin were used to constrain a series of reactive transport simulations of the Yellowstone geothermal system using TOUGHREACT. Three distinct stratigraphic units were encountered in the 153.4m deep Y-8 drill core: volcaniclastic sandstone, perlitic rhyolitic lava, and nonwelded pumiceous tuff. The main alteration phases identified in the Y-8 core samples include clay minerals, zeolites, silica polymorphs, adularia, and calcite. Temperatures observed in the Y-8 borehole increase with depth from sub-boiling conditions at the surface to a maximum of 169.8 ◦C at a depth of 104.1 m, with near-isothermal conditions persisting down to the well bottom. 1-D models of the Y-8 core hole were constructed to simulate the observed alteration mineral assemblage given the initial rock mineralogy and observed fluid chemistry and temperatures. Preliminary simulations involving the perlitic rhyolitic lava unit are consistent with the observed alteration of rhyolitic glass to form celadonite. © 2004 CNR. Published by Elsevier Ltd. All rights reserved