Catalogs of microseismicity are routinely compiled at geothermal reservoirs and provide valuable insights into reservoir structure and fluid movement. Hypocentral locations are typically used to infer the orientations of structures and constrain the extent of the permeable reservoir. However, frequency-magnitude distributions may contain additional, and underused, information about the distribution of pressure. Here, we present a four-year catalog of seismicity for the Rotokawa geothermal field in the central Taupō Volcanic Zone, New Zealand starting two years after the commissioning of the 140 MWe Nga Awa Purua power station. Using waveform-correlation-based signal detection we double the size of the previous earthquake catalog, refine the location and orientation of two reservoir faults and identify a new structure. We find the rate of seismicity to be insensitive to major changes in injection strategy during the study period, including the injectivity decline and shift of injection away from the dominant injector, RK24. We also map the spatial distribution of the earthquake frequency-magnitude distribution, or b-value, and show that it increases from ∼1.0 to ∼1.5 with increasing depth below the reservoir. As has been proposed at other reservoirs, we infer that these spatial variations reflect the distribution of pressure in the reservoir, where areas of high b-value correspond to areas of high pore-fluid pressure and a broad distribution of activated fractures. This analysis is not routinely conducted by geothermal operators but shows promise for using earthquake b-value as an additional tool for reservoir monitoring and management.

We investigate the relation between geothermal field production and fracture density and orientation in the Ngatamariki and Rotokawa geothermal fields, located in the Taupo Volcanic Zone, New Zealand using shear wave splitting (SWS). We determine the SWS parameters for 17 702 microseismic events across 38 stations spanning close to 4 yr from 2012 to 2015. We compare the strength of anisotropy to changes in field production and injection. We also compare the orientation of the anisotropy to in situ and regional measurements of maximum horizontal stress orientation. (SHmax). Due to the volume of unique events (approximately 160 000), shear wave phases are picked automatically. We carry out automatic SWS measurements using the Multiple Filter Automatic Splitting Technique (MFAST). The SWS measurements are interpreted in the context of stress aligned microcracks. Outside both fields and within Ngatamariki, fast polarizations align with the NE-SW regional orientation of SHmax. Within Rotokawa a greater complexity is observed, with polarizations tending toward N-S. We observe increases in per cent anisotropy coinciding with the start of production/injection inNgatamariki and then a later correlated drop in per cent anisotropy and vP/vS ratios in southern Ngatamariki as injection is shifted to the north. This relationship is consistent with pore fluid pressure within the reservoir being affected by local changes in production and injection of geothermal fluids causing cracks to open and close in response.

Fluid injection into the Earth's crust can induce seismic events that cause damage to local infrastructure but also offer valuable insight into seismogenesis. The factors that influence the magnitude, location, and number of induced events remain poorly understood but include injection flow rate and pressure as well as reservoir temperature and permeability. The relationship between injection parameters and injection-induced seismicity in high-temperature, high-permeability reservoirs has not been extensively studied. Here we focus on the Ngatamariki geothermal field in the central Taupō Volcanic Zone, New Zealand, where three stimulation/injection tests have occurred since 2012. We present a catalog of seismicity from 2012 to 2015 created using a matched-filter detection technique. We analyze the stress state in the reservoir during the injection tests from first motion-derived focal mechanisms, yielding an average direction of maximum horizontal compressive stress (SHmax) consistent with the regional NE-SW trend. However, there is significant variation in the direction of maximum compressive stress (σ1), which may reflect geological differences between wells. We use the ratio of injection flow rate to overpressure, referred to as injectivity index, as a proxy for near-well permeability and compare changes in injectivity index to spatiotemporal characteristics of seismicity accompanying each test. Observed increases in injectivity index are generally poorly correlated with seismicity, suggesting that the locations of microearthquakes are not coincident with the zone of stimulation (i.e., increased permeability). Our findings augment a growing body of work suggesting that aseismic opening or slip, rather than seismic shear, is the active process driving well stimulation in many environments.

Temperature and fluid pressure conditions control rock deformation and mineralization on geological faults, and hence the distribution of earthquakes. Typical intraplate continental crust has hydrostatic fluid pressure and a near-surface thermal gradient of 31 ± 15 degrees Celsius per kilometre. At temperatures above 300-450 degrees Celsius, usually found at depths greater than 10-15 kilometres, the intra-crystalline plasticity of quartz and feldspar relieves stress by aseismic creep and earthquakes are infrequent. Hydrothermal conditions control the stability of mineral phases and hence frictional-mechanical processes associated with earthquake rupture cycles, but there are few temperature and fluid pressure data from active plate-bounding faults. Here we report results from a borehole drilled into the upper part of the Alpine Fault, which is late in its cycle of stress accumulation and expected to rupture in a magnitude 8 earthquake in the coming decades. The borehole (depth 893 metres) revealed a pore fluid pressure gradient exceeding 9 ± 1 per cent above hydrostatic levels and an average geothermal gradient of 125 ± 55 degrees Celsius per kilometre within the hanging wall of the fault. These extreme hydrothermal conditions result from rapid fault movement, which transports rock and heat from depth, and topographically driven fluid movement that concentrates heat into valleys. Shear heating may occur within the fault but is not required to explain our observations. Our data and models show that highly anomalous fluid pressure and temperature gradients in the upper part of the seismogenic zone can be created by positive feedbacks between processes of fault slip, rock fracturing and alteration, and landscape development at plate-bounding faults.