Delineating Spatial Patterns in the Yellowstone Hydrothermal System using Geothermometry

Abstract:

Yellowstone National Park is unmatched with regard to its quantity of active hydrothermal features. Origins of thermal waters in its geyser basins have been traced to mixing of a deep parent water with meteoric waters in shallow local reservoirs (Fournier, 1989). A mineral-solution equilibrium model was developed to calculate water-rock chemical re-equilibration temperatures in these shallow reservoirs. We use the GeoT program, which uses water composition data as input to calculate saturation indices of selected minerals; the “best-clustering” minerals are then statistically determined to infer reservoir temperatures (Spycher et al., 2013). We develop the method using water composition data from Heart Lake Geyser Basin (HLGB), for which both chemical and isotopic geothermometers predict a reservoir water temperature of 205°C ± 10°C (Lowenstern et al., 2012), and minerals found in drill cores in Yellowstone’s geyser basins. We test the model for sensitivity to major element composition, pH, Total Inorganic Carbon (TIC) and selected minerals to optimize model parameters. Calculated temperatures are most accurate at pH values below 9.0, and closely match the equilibrium saturation indices of quartz, stilbite, microcline, and albite. The model is optimized with a TIC concentration that is consistent with the mass of diffuse CO2 flux in HLGB (Lowenstern et al., 2012). We then use water compositions from other thermal basins in Yellowstone in search of spatial variations in reservoir temperatures. We then compare the calculated temperatures with various SiO2 and cation geothermometers.