Characterizing the geothermal system at Surprise Valley (SV), northeastern California, is important for determining the sustainability of the energy resource, and mitigating hazards associated with hydrothermal eruptions that last occurred in 1951. Previous geochemical studies of the area attempted to reconcile different hot spring compositions on the western and eastern sides of the valley using scenarios of dilution, equilibration at low temperatures, surface evaporation, and differences in rock type along flow paths. These models were primarily supported using classical geothermometry methods, and generally assumed that fluids in the Lake City mud volcano area on the western side of the valley best reflect the composition of a deep geothermal fluid. In this contribution, we address controls on hot spring compositions using a different suite of geochemical tools, including optimized multicomponent geochemistry (GeoT) models, hot spring fluid major and trace element measurements, mineralogical observations, and stable isotope measurements of hot spring fluids and precipitated carbonates. We synthesize the results into a conceptual geochemical model of the Surprise Valley geothermal system, and show that high-temperature (quartz, Na/K, Na/K/Ca) classical geothermometers fail to predict maximum subsurface temperatures because fluids re-equilibrated at progressively lower temperatures during outflow, including in the Lake City area. We propose a model where hot spring fluids originate as a mixture between a deep thermal brine and modern meteoric fluids, with a seasonally variable mixing ratio. The deep brine has deuterium values at least 3 to 4‰ lighter than any known groundwater or high-elevation snow previously measured in and adjacent to SV, suggesting it was recharged during the Pleistocene when meteoric fluids had lower deuterium values. The deuterium values and compositional characteristics of the deep brine have only been identified in thermal springs and groundwater samples collected in proximity to structures that transmit thermal fluids, suggesting the brine may be thermal in nature. On the western side of the valley at the Lake City mud volcano, the deep brine-meteoric water mixture subsequently boils in the shallow subsurface, precipitates calcite, and re-equilibrates at about 130 °C. On the eastern side of the valley, meteoric fluid mixes to a greater extent with the deep brine, cools conductively without boiling, and the composition is modified as dissolved elements are sequestered by secondary minerals that form along the cooling and outflow path at temperatures <130 °C. Re-equilibration of geothermal fluids at lower temperatures during outflow explains why subsurface temperature estimates based on classical geothermometry methods are highly variable, and fail to agree with temperature estimates based on dissolved sulfate-oxygen isotopes and results of classical and multicomponent geothermometry applied to reconstructed deep well fluids. The proposed model is compatible with the idea suggested by others that thermal fluids on the western and eastern side of the valley have a common source, and supports the hypothesis that low temperature re-equilibration during west to east flow is the major control on hot spring fluid compositions, rather than dilution, evaporation, or differences in rock type.

The northernmost segment of the Gorda mid-ocean ridge is the site of a small-volume eruption in 1996 and the persistent off-axis Sea Cliff hydrothermal vent field. To better understand the geologic setting and formation of these features, 1-m resolution bathymetric mapping using autonomous underwater vehicles was completed in 2016. The mapped region covers 35 km2 and 15.6 km of the volcanic axis from south of the 1996 lava flows, and a cross section for ∼4.5 km perpendicular to the axis, that extends beyond the Sea Cliff hydrothermal vent field. A proposed 1996 flow ∼7 km south of previously mapped flows is an artifact from a poor pre-eruption survey. The 1996 flows consist of three discrete steep hummocky mounds of pillows and syneruptive talus. The Sea Cliff hydrothermal field is located a few km north of the narrowest, shallowest section of the ridge segment, 2.6 km east of the center of the neovolcanic zone, and ∼370 m above the average depth of the axial graben on the largest offset ridge-parallel fault. No evidence supports the prior hypothesis that the site is located where two fault systems intersect. The axial graben is asymmetrical with larger fault offsets on the east side. The ridge axis below the hydrothermal field and to the south toward the 1996 flows is constructed dominantly of hummocky flows of pillow basalt, many unusually steep-sided, with syneruptive talus at the base of their steep slopes. Three channelized flows ponded between steep hummocky flows, and then partially drained. Some low-eruption-rate hummocky flows and high-eruption-rate channelized flows have nearly identical compositions, supporting the idea that eruption rates on mid-ocean ridges vary because of different dike widths. Four volcanic structures with volumes between 0.18 and 0.25 km3 occur in the axial graben south of the 1996 flows. Two are flat-topped cones, another is a 1.5-km diameter inflated hummocky flow with 7 pit craters that demonstrate that the flow had a molten interior during growth. The fourth voluminous structure is a steep ridge with abundant syneruptive talus on its lower slopes. The North Gorda segment is an end-member, structurally and volcanically, compared with other Pacific intermediate-rate spreading ridges.

The region surrounding the Modoc Plateau, encompassing parts of northeastern California, southern Oregon, and northwestern Nevada, lies at an intersection between two tectonic provinces; the Basin and Range province and the Cascade volcanic arc. Both of these provinces have substantial geothermal resource base and resource potential. Geothermal systems with evidence of magmatic heat, associated with Cascade arc magmatism, typify the western side of the region. Systems on the eastern side of the region appear to be fault controlled with heat derived from high crustal heat flow, both of which are typical of the Basin and Range. As it has the potential to host Cascade arc-type geothermal resources, Basin and Range-type geothermal resources, and/or resources with characteristics of both provinces, and because there is relatively little current development, the Modoc Plateau region represents an intriguing potential for undiscovered geothermal resources. It remains unclear however, what specific set(s) of characteristics are diagnostic of Modoc-type geothermal systems and how or if those characteristics are distinct from Basin and Range-type or Cascade arc-type geothermal systems. In order to evaluate the potential for undiscovered geothermal resources in the Modoc area, we integrate a wide variety of existing data in order to evaluate geothermal resource potential and exploration risk utilizing ‘play-fairway’ analysis. We consider that the requisite parameters for hydrothermal circulation are: 1) heat that is sufficient to drive circulation, and 2) permeability that is sufficient to allow for fluid circulation in the subsurface. We synthesize data that indicate the extent and distribution of these parameters throughout the Modoc region. ‘Fuzzy logic’ is used to incorporate expert opinion into the utility of each dataset as an indicator of either heat or permeability, and thus geothermal favorability. The results identify several geothermal prospects, areas that are highly favorable for the occurrence of both heat and permeability. These are also areas where there is sufficient data coverage, quality, and consistency that the exploration risk is relatively low. These unknown, undeveloped, and under-developed prospects are well-suited for continued exploration efforts. The results also indicate to what degree the two ‘play-types,’ i.e. Cascade arc-type or Basin and Range-type, apply to each of the geothermal prospects, a useful guide in exploration efforts.