Isotope–geochemical characterization and geothermometrical modeling of Uttarakhand geothermal field, India

Abstract

Uttarakhand geothermal area, located in the central belt of the Himalayan geothermal province, is one of the important high temperature geothermal fields in India. In this study, the chemical characteristics of the thermal waters are investigated to identify the main geochemical processes affecting the composition of thermal waters during its ascent toward the surface as well as to determine the subsurface temperature of the feeding reservoir. The thermal waters are mainly Ca–Mg–HCO3 type with moderate silica and TDS concentrations. Mineral saturation states calculated from PHREEQC geochemical code indicate that thermal waters are supersaturated with respect to calcite, dolomite, aragonite, chalcedony, quartz (SI > 0), and undersaturated with respect to gypsum, anhydrite, and amorphous silica (SI < 0). XRD study of the spring deposit samples fairly corroborates the predicted mineral saturation state of the thermal waters. Stable isotopes (δ18O, δ2H) data confirm the meteoric origin of the thermal waters with no oxygen-18 shift. The mixing phenomenon between thermal water with shallow ground water is substantiated using tritium (3H) and chemical data. The extent of dilution is quantified using tritium content of thermal springs and non-thermal waters. Classical geothermometers, mixing model, and multicomponent fluid geothermometry modeling (GeoT) have been applied to estimate the subsurface reservoir temperature. Among different classical geothermometers, only quartz geothermometer provide somewhat reliable estimation (96–140 °C) of the reservoir temperature. GeoT modeling results suggest that thermal waters have attained simultaneous equilibrium with respect to minerals like calcite, quartz, chalcedony, brucite, tridymite, cristobalite, talc, at the temperature 130 ± 5 °C which is in good agreement with the result obtained from the mixing model.