USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13Linked to original sources
USGS1997

Stressing of the New Madrid seismic zone by a lower crust detachment fault

A new mechanical model for the cause of the New Madrid seismic zone in the central United States is analyzed. The model contains a subhorizontal detachment fault which is assumed to be near the domed top surface of locally thickened anomalous lower crust (“rift pillow”). Regional horizontal compression induces slip on the fault, and the slip creates a stress concentration in the upper crust above the rift pillow dome. In the coseismic stage of the model earthquake cycle, where the three largest magnitude 7–8 earthquakes in 1811–1812 are represented by a single model mainshock on a vertical northeast trending fault, the model mainshock has a moment equivalent to a magnitude 8 event. During the interseismic stage, corresponding to the present time, slip on the detachment fault exerts a right-lateral shear stress on the locked vertical fault whose failure produces the model mainshock. The sense of shear is generally consistent with the overall sense of slip of 1811–1812 and later earthquakes. Predicted rates of horizontal strain at the ground surface are about 10 −7 year −1 and are comparable to some observed rates. The model implies that rift pillow geometry is a significant influence on the maximum possible earthquake magnitude.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas fault

The Parkfield, California, segment of the San Andreas fault is transitional in character between the creeping segment of the fault to the northwest and the locked Carrizo Plain segment to the southeast. The rate of shallow fault slip decreases from 25–30 mm/yr northwest of the epicenter of the 1966 Parkfield earthquake to zero at the southeastern end of the 1966 rupture zone. Data from a network of trilateration lines spanning the San Andreas fault near Parkfield and extending to the Pacific coast near San Luis Obispo shed light on the rate of fault slip at depth since the 1966 earthquake. In this study, average rates of line length change and shallow fault slip were inverted to determine the slip rate at depth on the Parkfield fault segment. The fault is taken to be a vertical surface with unknown distribution of strike-slip displacement in an elastic half-space. A striking result of the inversions is that all solutions providing acceptable fits to the data exhibit a locked zone essentially coincident with the rupture surface of the 1966 Parkfield earthquake. The data require that the locked zone extend nearly as far north as the 1966 epicenter; however, the vertical extent of the locked zone is not well resolved. Over much of the Parkfield segment the fault is slipping faster at the earth's surface than it is at seismogenic depths. In order to fit the trilateration measurements it is necessary to include a component of contraction normal to the trend of the San Andreas. The inversion results suggest a spatially uniform normal strain of −0.06 μstrain/yr. The orientation of the contraction is compatible with geologic and seismic evidence of active folding and reverse faulting in the region. The magnitude of the contraction is consistent with convergence rates inferred from global plate motion models.

Journal of Geophysical Research Solid Earth↗
USGS2003

Three-dimensional crustal structure of Long Valley caldera, California, and evidence for the migration of CO2 under Mammoth Mountain

A temporary network of 69 three-component seismic stations captured a major seismic sequence in Long Valley caldera in 1997. We performed a tomographic inversion for crustal structure beneath a 28 km ?? 16 km area encompassing part of the resurgent dome, the south moat, and Mammoth Mountain. Resolution of crustal structure beneath the center of the study volume was good down to ???3 km below sea level (???5 km below the surface). Relatively high wave speeds are associated with the Bishop Tuff and lower wave speeds characterize debris in the surrounding moat. A low-Vp/Vs anomaly extending from near the surface to ???1 km below sea level beneath Mammoth Mountain may represent a CO2 reservoir that is supplying CO2-rich springs, venting at the surface, and killing trees. We investigated temporal variations in structure beneath Mammoth Mountain by differencing our results with tomographic images obtained using data from 1989/1990. Significant changes in both Vp and Vs were consistent with the migration of CO2 into the upper 2 km or so beneath Mammoth Mountain and its depletion in peripheral volumes that correlate with surface venting areas. Repeat tomography is capable of detecting the migration of gas beneath active silicic volcanoes and may thus provide a useful volcano monitoring tool.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Tidal calibration of Plate Boundary Observatory borehole strainmeters: Roles of vertical and shear coupling

A multicomponent borehole strainmeter directly measures changes in the diameter of its cylindrical housing at several azimuths. To transform these measurements to formation strains requires a calibration matrix, which must be estimated by analyzing the installed strainmeter's response to known strains. Typically, theoretical calculations of Earth tidal strains serve as the known strains. This paper carries out such an analysis for 12 Plate Boundary Observatory (PBO) borehole strainmeters, postulating that each of the strainmeters' four gauges responds ("couples") to all three horizontal components of the formation strain tensor, as well as to vertical strain. Orientation corrections are also estimated. The fourth extensometer in each PBO strainmeter provides redundant information used to reduce the chance that coupling coefficients could be misleadingly fit to inappropriate theoretical tides. Satisfactory fits between observed and theoretically calculated tides were obtained for three PBO strainmeters in California, where the calculated tides are corroborated by other instrumentation, as well as for six strainmeters in Oregon and Washington, where no other instruments have ever recorded Earth tidal strain. Several strainmeters have unexpectedly large coupling coefficients for vertical strain, which increases the strainmeter's response to atmospheric pressure. Vertical coupling diminishes, or even changes the sign of, the apparent response to areal strain caused by Earth tides or deep Earth processes because near the free surface, vertical strains are opposite in sign to areal strain. Vertical coupling does not impair the shear strain response, however. PBO borehole strainmeters can provide calibrated shear strain time series of transient strain associated with tectonic or magmatic processes.

Journal of Geophysical Research↗
USGS2010

Inelastic off-fault response and three-dimensional dynamics of earthquake rupture on a strike-slip fault

Large dynamic stress off the fault incurs an inelastic response and energy loss, which contributes to the fracture energy, limiting the rupture and slip velocity. Using an explicit finite element method, we model three-dimensional dynamic ruptures on a vertical strike-slip fault in a homogeneous half-space. The material is subjected to a pressure-dependent Drucker-Prager yield criterion. Initial stresses in the medium increase linearly with depth. Our simulations show that the inelastic response is confined narrowly to the fault at depth. There the inelastic strain is induced by large dynamic stresses associated with the rupture front that overcome the effect of the high confining pressure. The inelastic zone increases in size as it nears the surface. For material with low cohesion (~5 MPa) the inelastic zone broadens dramatically near the surface, forming a "flowerlike" structure. The near-surface inelastic strain occurs in both the extensional and the compressional regimes of the fault, induced by seismic waves ahead of the rupture front under a low confining pressure. When cohesion is large (~10 MPa), the inelastic strain is significantly reduced near the surface and confined mostly to depth. Cohesion, however, affects the inelastic zone at depth less significantly. The induced shear microcracks show diverse orientations near the surface, owing to the low confining pressure, but exhibit mostly horizontal slip at depth. The inferred rupture-induced anisotropy at depth has the fast wave direction along the direction of the maximum compressive stress.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Evidence for water influx from a caldera lake during the explosive hydromagmatic eruption of 1790, Kilauea volcano, Hawaii

In 1790 a major hydromagmatic eruption at the summit of Kilauea volcano, Hawaii, deposited up to 10 m of pyroclastic fall and surge deposits and killed several dozen Hawaiian natives who were crossing the island. Previous studies have hypothesized that the explosivity of this eruption was due to the influx of groundwater into the conduit and mixing of the groundwater with ascending magma. This study proposes that surface water, not groundwater, was the agent responsible for the explosiveness of the eruption. That is, a lake or pond may have existed in the caldera in 1790 and explosions may have taken place when magma ascended into the lake from below. That assertion is based on two lines of evidence: (1) high vesicularity (averaging 73% of more than 3000 lapilli) and high vesicle number density (10 5 –10 7 cm −3 melt) of pumice clasts suggest that some phases of the eruption involved vigorous, sustained magma ascent; and (2) numerical calculations suggest that under most circumstances, hydrostatic pressure would not be sufficient to drive water into the eruptive conduit during vigorous magma ascent unless the water table were above the ground surface. These results are supported by historical data on the rate of infilling of the caldera floor during the early 1800s. When extrapolated back to 1790, they suggest that the caldera floor was below the water table.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Three-dimensional surface deformation derived from airborne interferometric UAVSAR: Application to the Slumgullion Landslide

In order to provide surface geodetic measurements with “landslide-wide” spatial coverage, we develop and validate a method for the characterization of 3-D surface deformation using the unique capabilities of the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne repeat-pass radar interferometry system. We apply our method at the well-studied Slumgullion Landslide, which is 3.9 km long and moves persistently at rates up to ∼2 cm/day. A comparison with concurrent GPS measurements validates this method and shows that it provides reliable and accurate 3-D surface deformation measurements. The UAVSAR-derived vector velocity field measurements accurately capture the sharp boundaries defining previously identified kinematic units and geomorphic domains within the landslide. We acquired data across the landslide during spring and summer and identify that the landslide moves more slowly during summer except at its head, presumably in response to spatiotemporal variations in snowmelt infiltration. In order to constrain the mechanics controlling landslide motion from surface velocity measurements, we present an inversion framework for the extraction of slide thickness and basal geometry from dense 3-D surface velocity fields. We find that the average depth of the Slumgullion Landslide is 7.5 m, several meters less than previous depth estimates. We show that by considering a viscoplastic rheology, we can derive tighter theoretical bounds on the rheological parameter relating mean horizontal flow rate to surface velocity. Using inclinometer data for slow-moving, clay-rich landslides across the globe, we find a consistent value for the rheological parameter of 0.85 ± 0.08.

Colorado↗
USGS1998

Seismotectonics of the Loma Prieta, California, region determined from three-dimensional Vp, Vp/Vs, and seismicity

Three-dimensional V p and V p / V s velocity models for the Loma Prieta region were developed from the inversion of local travel time data (21,925 P arrivals and 1,116 S arrivals) from earthquakes, refraction shots, and blasts recorded on 1700 stations from the Northern California Seismic Network and numerous portable seismograph deployments. The velocity and density models and microearthquake hypocenters reveal a complex structure that includes a San Andreas fault extending to the base of the seismogenic layer. A body with high V p extends the length of the rupture and fills the 5 km wide volume between the Loma Prieta mainshock rupture and the San Andreas and Sargent faults. We suggest that this body controls both the pattern of background seismicity on the San Andreas and Sargent faults and the extent of rupture during the mainshock, thus explaining how the background seismicity outlined the along-strike and depth extent of the mainshock rupture on a different fault plane 5 km away. New aftershock focal mechanisms, based on three-dimensional ray tracing through the velocity model, support a heterogeneous postseismic stress field and can not resolve a uniform fault normal compression. The subvertical (or steeply dipping) San Andreas fault and the fault surfaces that ruptured in the 1989 Loma Prieta earthquake are both parts of the San Andreas fault zone and this section of the fault zone does not have a single type of characteristic event.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Crustal velocity field near the big bend of California's San Andreas fault

We use geodetic data spanning the 1920–1992 interval to estimate the horizontal velocity field near the big bend segment of California's San Andreas fault (SAF). More specifically, we estimate a horizontal velocity vector for each node of a two-dimensional grid that has a 15-min-by-15-min mesh and that extends between latitudes 34.0°N and 36.0°N and longitudes 117.5°W and 120.5°W. For this estimation process, we apply bilinear interpolation to transfer crustal deformation information from geodetic sites to the grid nodes. The data include over a half century of triangulation measurements, over two decades of repeated electronic distance measurements, a decade of repeated very long baseline interferometry measurements, and several years of Global Positioning System measurements. Magnitudes for our estimated velocity vectors have formal standard errors ranging from 0.7 to 6.8 mm/yr. Our derived velocity field shows that (1) relative motion associated with the SAF exceeds 30 mm/yr and is distributed on the Earth's surface across a band (>100 km wide) that is roughly centered on this fault; (2) when velocities are expressed relative to a fixed North America plate, the motion within our primary study region has a mean orientation of N44°W ± 2° and the surface trace of the SAF is congruent in shape to nearby contours of constant speed yet this trace is oriented between 5° and 10° counterclockwise relative to these contours; and (3) large strain rates (shear rates > 150 nrad/yr and/or areal dilatation rates < −150 nstr/yr) exist near the Garlock fault, near the White Wolf fault, and in the Ventura basin.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Dislocation pileup as a representation of strain accumulation on a strike-slip fault

The conventional model of strain accumulation on a vertical transform fault is a discrete screw dislocation in an elastic half-space with the Burgers vector of the dislocation increasing at the rate of relative plate motion. It would be more realistic to replace that discrete dislocation by a dislocation distribution, presumably a pileup in which the individual dislocations are in equilibrium. The length of the pileup depends upon the applied stress and the amount of slip that has occurred at depth. I argue here that the dislocation pileup (the transition on the fault from no slip to slip at the full plate rate) occupies a substantial portion of the lithosphere thickness. A discrete dislocation at an adjustable depth can reproduce the surface deformation profile predicted by a pileup so closely that it will be difficult to distinguish between the two models. The locking depth (dislocation depth) of that discrete dislocation approximation is substantially (???30%) larger than that (depth to top of the pileup) in the pileup model. Thus, in inverting surface deformation data using the discrete dislocation model, the locking depth in the model should not be interpreted as the true locking depth. Although dislocation pileup models should provide a good explanation of the surface deformation near the fault trace, that explanation may not be adequate at greater distances from the fault trace because approximating the expected horizontally distributed deformation at subcrustal depths by uniform slip concentrated on the fault is not justified.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Vertical coseismic offsets from differential high-resolution stereogrammetric DSMs: The 2013 Baluchistan, Pakistan earthquake

The recent proliferation of high-resolution (< 3-m spatial resolution) digital topography datasets opens a spectrum of geodetic applications in differential topography, including the quantification of coseismic vertical displacement fields. Most investigations of coseismic vertical displacements to date rely, in part, on pre- or post-event lidar surveys that are intractable or non-existent in many locales. Stereogrammetric digital surface models (DSMs) derived from high-resolution satellite optical imagery provide a new avenue for the retrieval of spatially-dense vertical coseismic displacements on a global scale. In this study, we generated 2-m resolution pre- and post-seismic DSMs from satellite optical imagery spanning the 2013 Mw7.7 Baluchistan strike-slip earthquake that occurred on the Hoshab fault in southern Pakistan. We applied the Iterative Closest Point algorithm to the DSMs to quantify the coseismic vertical displacement field at a spatial resolution of 10-30 m and to generate 3D coseismic strain tensors. We found that across-fault vertical offsets alternated between uplift and subsidence and varied between ~1-3 m in a non-systematic manner along the Hoshab fault. We show that the pre-existing topography and near-fault geomorphology are variably consistent and inconsistent with the displacement kinematics of the 2013 earthquake, and we argue that these relationships highlight varied slip sense history along the Hoshab fault. Notably, topography along the southern extents of the Hoshab fault requires different surface displacement kinematics than occurred in the 2013 earthquake, suggesting that the Hoshab fault accommodates varying senses of slip (bimodal slip) through time.

Baluchistan Province↗
USGS1996

Rheology of water ices V and VI

We have measured the mechanical strength (σ) of pure water ices V and VI under steady state deformation conditions. Constant displacement rate compressional tests were conducted in a gas apparatus at confining pressures from 400 < P < 800 MPa, temperatures from 209 < T < 270 K, and strain rates 7 × 10 −7 < < 7 × 10 −4 s −1 . Most of the results fit to an empirical flow law of the form = A σ n exp {−( E * + PV *)/ RT }, where the four material constants A , n , E *, and V * are (for in inverse seconds and P and σ in megapascals) 10 23.0 , 6.0, 136 kJ/mol, and 29 cm 3 /mol, respectively, for ice V; and 10 6.7 , 4.5, 66 kJ/mol, and 11 cm 3 /mol, respectively, for ice VI. Ice VI may weaken to a mechanism of higher E * at T > 250 K. Ices V and VI are thus rheologically distinct but by coincidence have approximately the same strength under the conditions chosen for these experiments. To avoid misidentification, these tests are therefore accompanied by careful observations of the occurrences and characteristics of phase changes. One sample each of ice V and VI was quenched at pressure to metastably retain the high-pressure phase and the acquired deformation microstructures; X ray diffraction analysis of these samples confirmed the phase identification. Surface replicas of the deformed and quenched samples suggest that ice V probably deforms largely by dislocation creep, while ice VI deforms by a more complicated process involving substantial grain size reduction through recrystallization.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Seismic-geodynamic constraints on three-dimensional structure, vertical flow, and heat transfer in the mantle

Joint inversions of seismic and geodynamic data are carried out in which we simultaneously constrain global-scale seismic heterogeneity in the mantle as well as the amplitude of vertical mantle flow across the 670 km seismic discontinuity. These inversions reveal the existence of a family of three-dimensional (3-D) mantle models that satisfy the data while at the same time yielding predictions of layered mantle flow. The new 3-D mantle models we obtain demonstrate that the buoyancy forces due to the undulations of the 670 km phase-change boundary strongly inhibit the vertical flow between the upper and lower mantle. The strong stabilizing effect of the 670 km topography also has an important impact on the predicted dynamic topography of the Earth's solid surface and on the surface gravity anomalies. The new 3-D models that predict strongly or partially layered mantle flow provide essentially identical fits to the global seismic data as previous models that have, until now, predicted only whole-mantle flow. The convective vertical transport of heat across the mantle predicted on the basis of the new 3-D models shows that the heat flow is a minimum at 1000 km depth. This suggests the presence at this depth of a globally defined horizon across which the pattern of lateral heterogeneity changes rapidly.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Postseismic relaxation associated with transient creep rheology

[1] Perfettini and Avouac (2004) postulated that both the aftershock rate (assumed proportional to the local stressing rate) and the postseismic relaxation are driven by the loading imposed by postseismic slip on the brittle creep fault zone (BCFZ), the downdip extension of the fault zone below the coseismic rupture. I explore the consequences of that hypothesis for a long, strike-slip fault in the case where the BCFZ rheology is compatible with ordinary transient creep (creep strain proportional to log e (1 + t / τ 2 )). Because the important relaxation occurs near the bottom of the coseismic rupture, I calculate the postearthquake response with a model in which the BCFZ is represented by a viscoelastic half-space below the coseismic rupture. I find that both the predicted postseismic relaxation and the cumulative number of aftershocks can be approximated by the same temporal dependence N MO ( t ) = a MO (1−(1 + t / τ ) 1− p )/( p − 1), where t is the time after the earthquake and a MO , τ , and p are the constants chosen to fit either data set. Notice that d N MO ( t )/d t = ( a MO / τ )/(1 + t / τ ) p is the modified Omori law used to describe the rate of aftershock occurrence. Thus, the modified Omori law can be understood as a consequence of the Perfettini–Avouac hypothesis (aftershocks driven by slip on the BCFZ) and a BCFZ rheology compatible with ordinary transient creep. Moreover, the temporal dependence N MO ( t ) has been shown to fit postseismic surface deformation following at least 9 earthquakes. I also show that the conventional, one-dimensional, spring-block model of a BFCZ with a rheology compatible with ordinary transient creep leads to the same temporal dependence ( N MO ( t )).

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Magma intrusion and volatile ascent beneath Norris Geyser Basin, Yellowstone National Park

Recent activity has provided new insights into the causes of surface deformation in and around the Yellowstone Caldera, a topic that has been debated since the discovery of caldera floor uplift more than four decades ago. An episode of unusually rapid uplift (>15 cm/yr) centered near Norris Geyser Basin along the north caldera rim began in late 2013 and continued until a M w 4.9 earthquake on 30 March 2014; thereafter, uplift abruptly switched to subsidence. Uplift at rates of several centimeters per year resumed in 2016 and continued at least through the end of 2018. Modeling of Global Positioning System and interferometric synthetic aperture radar data suggests an evolving process of deep magma intrusion during 1996–2001 followed by volatile ascent and accumulation at shallow levels, perhaps as shallow as a few hundred meters depth. The depth of shallow volatile accumulation appears to have shallowed from the 2014 to the 2016 deformation episode, and frequent eruptions of Steamboat Geyser since March 2018 are likely a surface manifestation of this ongoing process. Hydrothermal explosion features are prominent in the Norris Geyser Basin area, and the apparent shallow nature of the volatile accumulation implies an increased risk of hydrothermal explosions.

Wyoming↗
USGS1996

Coral ages and island subsidence, Hilo drill hole

A 25.8-m-thick sedimentary section containing coral fragments occurs directly below a surface lava flow (the ∼1340 year old Panaewa lava flow) at the Hilo drill hole. Ten coral samples from this section dated by accelerator mass spectrometry (AMS) radiocarbon and five by thermal infrared multispectral scanner (TIMS) 230 Th/U methods show good agreement. The calcareous unit is 9790 years old at the bottom and 1690 years old at the top and was deposited in a shallow lagoon behind an actively growing reef. This sedimentary unit is underlain by a 34-m-thick lava flow which in turn overlies a thin volcaniclastic silt with coral fragments that yield a single 14 C date of 10,340 years. The age-depth relations of the dated samples can be compared with proposed eustatic sea level curves after allowance for island subsidence is taken. Island subsidence averages 2.2 mm/yr for the last 47 years based on measurements from a tide gage near the drill hole or 2.5–2.6 mm/yr for the last 500,000 years based on the ages and depths of a series of drowned coral reefs offshore from west Hawaii. The age-depth measurements of coral fragments are more consistent with eustatic sea levels as determined by coral dating at Barbados and Albrolhos Islands than those based on oxygen isotopic data from deep sea cores. The Panaewa lava flow entered a lagoon underlain by coral debris and covered the drill site with 30.9 m of lava of which 11 m was above sea level. This surface has now subsided to 4.2 m above sea level, but it demonstrates how a modern lava flow entering Hilo Bay would not only change the coastline but could extensively modify the offshore shelf.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Evolution of melt-vapor surface tension in silicic volcanic systems: Experiments with hydrous melts

We evaluate the melt‐vapor surface tension (σ) of natural, water‐saturated dacite melt at 200 MPa, 950–1055°C, and 4.8–5.7 wt % H 2 O. We experimentally determine the critical supersaturation pressure for bubble nucleation as a function of dissolved water and then solve for σ at those conditions using classical nucleation theory. The solutions obtained give dacite melt‐vapor surface tensions that vary inversely with dissolved water from 0.042 (±0.003) J m −2 at 5.7 wt % H 2 O to 0.060 (±0.007) J m −2 at 5.2 wt % H 2 O to 0.073 (±0.003) J m −2 at 4.8 wt % H 2 O. Combining our dacite results with data from published hydrous haplogranite and high‐silica rhyolite experiments reveals that melt‐vapor surface tension also varies inversely with the concentration of mafic melt components (e.g., CaO, FeO total , MgO). We develop a thermodynamic context for these observations in which melt‐vapor surface tension is represented by a balance of work terms controlled by melt structure. Overall, our results suggest that cooling, crystallization, and vapor exsolution cause systematic changes in σ that should be considered in dynamic modeling of magmatic processes.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Geodetic analysis of reservoir depletion at the Geyser steam field in northern California

Reservoir depletion at the Geysers from 1974 to 1977 is evident in measured changes in gravity, surface strain, and pore pressure drainage. The drainage area increased about 20%, the maximum gravity decrease was about −120 μGal, and the maximum elevation change was about 6 cm during this period. Since the net mass withdrawal is known, it may be combined with the gravity change to estimate a drainage volume. The maximum drainage volume is 25 km 3 . Because the depth of the caprock and extent of the pore pressure drainage is known, this volume limit implies that no significant mass withdrawal occurs below a depth of 4 km. The ratios of surface elevation changes to horizontal contraction coinciding with the drainage area imply an equant drainage geometry, assuming that fluid production produces negative dilatation. Using the same cylindrical geometry used to model the gravity, negative dilatational strain rates of 4 to 5 × 10 −5 /yr were found to produce the observed surface displacements. The likelihood of boiling in the system to produce steam combined with apparent large bulk moduli (from seismic velocities) and small pore pressure declines suggest that most of the reservoir contraction is due to cooling.

Journal of Geophysical Research Solid Earth↗