USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research Solid Earth”

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 541 records · Page 30Linked to original sources
USGS2002

A Holocene paleosecular variation from 14C-dated volcanic rocks in western North America

A paleosecular variation (PSV) curve for western North America is presented on the basis of 94 virtual geomagnetic poles (VGPs) from dated volcanic rocks sampled at 446 sites. Approximately 60% of the paleomagnetic database has been previously published. A curve defined by “spherical smoothed splines” is fitted to the VGPs, ranked by the quality of the age determinations, where the data density is highest between 3690 and −30 years before present (B.P.) (A.D. 1950), between 7800 and 7050 years B.P., and between 14,060 and 12,700 years B.P. The younger segments of the curve derived from volcanic rocks are similar but less complex than other high-resolution PSV curves derived from lacustrine sediments, particularly the record at Fish Lake, Oregon. The PSV record from lava flows (PSVL), however, is perhaps more reliable in its general shape and chronology because of the higher fidelity of volcanic rocks as magnetic field recorders and because of the greater density of 14 C dates. The new PSVL record provides a partial Holocene master curve for western North America and will be of particular value in dating geological and archeological materials using paleomagnetic directions.

Arizona, California, Colorado, Idaho, New Mexico, ↗
USGS2003

The crustal structure from the Altai Mountains to the Altyn Tagh fault, northwest China

[1] We present a new crustal section across northwest China based on a seismic refraction profile and geologic mapping. The 1100‐km‐long section crosses the southern margin of the Chinese Altai Mountains, Junggar Accretional Belt and eastern Junggar basin, easternmost Tianshan Mountains, and easternmost Tarim basin. The crustal velocity structure and Poisson's ratio (σ), which provide a constraint on crustal composition, were determined from P and S wave data. Despite the complex geology, the crustal thickness along the entire profile is nearly uniform at 50 km. The thickest crust (56 km) occurs at the northern end of the profile beneath the Altai Mountains and the thinnest (46 km) crust is beneath the Junggar basin. Beneath surficial sediments, the crust is found to have three layers with P wave velocities ( V p ) of 6.0–6.3, 6.3–6.6, and 6.9–7.0 km/s, respectively. The southern half of the profile, including the eastern Tianshan Mountains and eastern margin of the Tarim basin, shows low P wave velocities and σ = 0.25 to a depth of 30 km, which suggests a quartz‐rich, granitic upper crustal composition. The northern half of the profile below the Altai Mountains and Junggar Accretional Belt has a higher Poisson's ratio of σ = 0.26–0.27 to a depth of 30 km, indicative of an intermediate crustal composition. The entire 1100‐km‐long profile is underlain by a 15–30 km thick high velocity (6.9–7.0 km/s; σ = 0.26–0.28) lower‐crustal layer that we interpret to have a bulk composition of mafic granulite. At the southern end of the profile, a 5‐km‐thick midcrustal low‐velocity layer ( V p = 5.9 km/s, σ = 0.25) underlies the Tianshan and the region to the south, and may be indicative of a near‐horizontal detachment interface. P n velocities are ∼7.7–7.8 km/s between the Tianshan and the Junggar basin, and ∼7.9–8.0 km/s below the Altai Mountains and eastern margin of the Tarim basin. We interpret the consistent three‐layer stratification of the crust to indicate that the crust has undergone partial melting and differentiation after Paleozoic terrane accretion. The thickness (50 km) of the crust appears to be related to compression resulting from the Indo‐Asian collision.

Altai Mountains↗
USGS2002

Magnetostratigraphy, paleomagnetic correlation, and deformation of pleistocene deposits in the south central Puget Lowland, Washington

Paleomagnetic results from Pleistocene sedimentary deposits in the central Puget Lowland indicate that the region has experienced widespread deformation within the last 780 kyr. Three oriented samples were collected from unaltered fine-grained sediments mostly at sea level to determine the magnetostratigraphy at 83 sites. Of these, 47 have normal, 18 have reversed, and 18 have transitional (8 localities) polarities. Records of reversed- to normal-polarity transitions of the geomagnetic field were found in thick sections of silt near the eastern end of the Tacoma Narrows Bridge, and again at Wingehaven Park near the northern tip of Vashon Island. The transitional horizons, probably related to the Bruhnes-Matuyama reversal, apparently fall between previously dated Pleistocene sediments at the Puyallup Valley type section (all reversed-polarity) to the south and the Whidbey Island type section (all normal-polarity) to the north. The samples, in general, are of sufficient quality to record paleosecular variation (PSV) of the geomagnetic field, and a statistical technique is used to correlate horizons with significant agreement in their paleomagnetic directions. Our data are consistent with the broad structures of the Seattle uplift inferred at depth from seismic reflection, gravity, and aeromagnetic profiles, but the magnitude of vertical adjustments is greatly subdued in the Pleistocene deposits.

Washington↗
USGS2002

The 12 September 1999 Upper East Rift Zone dike intrusion at Kilauea Volcano, Hawaii

Deformation associated with an earthquake swarm on 12 September 1999 in the Upper East Rift Zone of Kilauea Volcano was recorded by continuous GPS receivers and by borehole tiltmeters. Analyses of campaign GPS, leveling data, and interferometric synthetic aperture radar (InSAR) data from the ERS-2 satellite also reveal significant deformation from the swarm. We interpret the swarm as resulting from a dike intrusion and model the deformation field using a constant pressure dike source. Nonlinear inversion was used to find the model that best fits the data. The optimal dike is located beneath and slightly to the west of Mauna Ulu, dips steeply toward the south, and strikes nearly east-west. It is approximately 3 by 2 km across and was driven by a pressure of ∼15 MPa. The total volume of the dike was 3.3 × 10 6 m 3 . Tilt data indicate a west to east propagation direction. Lack of premonitory inflation of Kilauea's summit suggests a passive intrusion; that is, the immediate cause of the intrusion was probably tensile failure in the shallow crust of the Upper East Rift Zone brought about by persistent deep rifting and by continued seaward sliding of Kilauea's south flank.

Hawaii↗
USGS2002

Strain accumulation and rotation in western Oregon and southwestern Washington

Velocities of 75 geodetic monuments in western Oregon and southwestern Washington extending from the coast to more than 300 km inland have been determined from GPS surveys over the interval 1992–2000. The average standard deviation in each of the horizontal velocity components is ∼1 mm yr −1 . The observed velocity field is approximated by a combination of rigid rotation (Euler vector relative to interior North America: 43.40°N ± 0.14°, 119.33°W ± 0.28°, and 0.822 ± 0.057° Myr −1 clockwise; quoted uncertainties are standard deviations), uniform regional strain rate (ε EE = −7.4 ± 1.8, ε EN = −3.4 ± 1.0, and ε NN = −5.0 ± 0.8 nstrain yr −1 , extension reckoned positive), and a dislocation model representing subduction of the Juan de Fuca plate beneath North America. Subduction south of 44.5°N was represented by a 40‐km‐wide locked thrust and subduction north of 44.5°N by a 75‐km‐wide locked thrust.

Oregon, Washington↗
USGS2002

High-resolution image of Calaveras Fault seismicity

By measuring relative earthquake arrival times using waveform cross correlation and locating earthquakes using the double difference technique, we are able to reduce hypocentral errors by 1 to 2 orders of magnitude over routine locations for nearly 8000 events along a 35-km section of the Calaveras Fault. This represents ∼92% of all seismicity since 1984 and includes the rupture zone of the M 6.2 1984 Morgan Hill, California, earthquake. The relocated seismicity forms highly organized structures that were previously obscured by location errors. There are abundant repeating earthquake sequences as well as linear clusters of earthquakes. Large voids in seismicity appear with dimensions of kilometers that have been aseismic over the 30-year time interval, suggesting that these portions of the fault are either locked or creeping. The area of greatest slip in the Morgan Hill main shock coincides with the most prominent of these voids, suggesting that this part of the fault may be locked between large earthquakes. We find that the Calaveras Fault at depth is extremely thin, with an average upper bound on fault zone width of 75 m. Given the location error, however, this width is not resolvably different from zero. The relocations reveal active secondary faults, which we use to solve for the stress field in the immediate vicinity of the Calaveras Fault. We find that the maximum compressive stress is at a high angle, only 13° from the fault normal, supporting previous interpretations that this fault is weak.

California↗
USGS2002

Global Omori law decay of triggered earthquakes: Large aftershocks outside the classical aftershock zone

[1] Triggered earthquakes can be large, damaging, and lethal as evidenced by the1999 shocks in Turkey and the 2001 earthquakes in El Salvador. In this study, earthquakes with M s ≥ 7.0 from the Harvard centroid moment tensor (CMT) catalog are modeled as dislocations to calculate shear stress changes on subsequent earthquake rupture planes near enough to be affected. About 61% of earthquakes that occurred near (defined as having shear stress change ∣Δτ∣ ≥ 0.01 MPa) the M s ≥ 7.0 shocks are associated with calculated shear stress increases, while ∼39% are associated with shear stress decreases. If earthquakes associated with calculated shear stress increases are interpreted as triggered, then such events make up at least 8% of the CMT catalog. Globally, these triggered earthquakes obey an Omori law rate decay that lasts between ∼7–11 years after the main shock. Earthquakes associated with calculated shear stress increases occur at higher rates than background up to 240 km away from the main shock centroid. Omori's law is one of the few time‐predictable patterns evident in the global occurrence of earthquakes. If large triggered earthquakes habitually obey Omori's law, then their hazard can be more readily assessed. The characteristic rate change with time and spatial distribution can be used to rapidly assess the likelihood of triggered earthquakes following events of M s ≥ 7.0. I show an example application to the M = 7.7 13 January 2001 El Salvador earthquake where use of global statistics appears to provide a better rapid hazard estimate than Coulomb stress change calculations.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Temporal evolution of a hydrothermal system in Kusatsu-Shirane Volcano, Japan, inferred from the complex frequencies of long-period events

We present a detailed description of temporal variations in the complex frequencies of long-period (LP) events observed at Kusatsu-Shirane Volcano. Using the Sompi method, we analyze 35 LP events that occurred during the period from August 1992 through January 1993. The observed temporal variations in the complex frequencies can be divided into three periods. During the first period the dominant frequency rapidly decreases from 5 to 1 Hz, and Q of the dominant spectral peak remains roughly constant with an average value near 100. During the second period the dominant frequency gradually increases up to 3 Hz, and Q gradually decreases from 160 to 30. During the third period the dominant frequency increases more rapidly from 3 to 5 Hz, and Q shows an abrupt increase at the beginning of this period and then remains roughly constant with an average value near 100. Such temporal variations can be consistently explained by the dynamic response of a hydrothermal crack to a magmatic heat pulse. During the first period, crack growth occurs in response to the overall pressure increase in the hydrothermal system caused by the heat pulse. Once crack formation is complete, heat gradually changes the fluid in the crack from a wet misty gas to a dry gas during the second period. As heating of the hydrothermal system gradually subsides, the overall pressure in this system starts to decrease, causing the collapse of the crack during the third period.

Kusatsu-Shirane Volcano↗
USGS2002

Crustal structure and relocated earthquakes in the Puget Lowland, Washington, from high-resolution seismic tomography

The availability of regional earthquake data from the Pacific Northwest Seismograph Network (PNSN), together with active source data from the Seismic Hazards Investigation in Puget Sound (SHIPS) seismic experiments, has allowed us to construct a new high-resolution 3-D, P wave velocity model of the crust to a depth of about 30 km in the central Puget Lowland. In our method, earthquake hypocenters and velocity model are jointly coupled in a fully nonlinear tomographic inversion. Active source data constrain the upper 10–15 km of the model, and earthquakes constrain the deepest portion of the model. A number of sedimentary basins are imaged, including the previously unrecognized Muckleshoot basin, and the previously incompletely defined Possession and Sequim basins. Various features of the shallow crust are imaged in detail and their structural transitions to the mid and lower crust are revealed. These include the Tacoma basin and fault zone, the Seattle basin and fault zone, the Seattle and Port Ludlow velocity highs, the Port Townsend basin, the Kingston Arch, and the Crescent basement, which is arched beneath the Lowland from its surface exposure in the eastern Olympics. Strong lateral velocity gradients, consistent with the existence of previously inferred faults, are observed, bounding the southern Port Townsend basin, the western edge of the Seattle basin beneath Dabob Bay, and portions of the Port Ludlow velocity high and the Tacoma basin. Significant velocity gradients are not observed across the southern Whidbey Island fault, the Lofall fault, or along most of the inferred location of the Hood Canal fault. Using improved earthquake locations resulting from our inversion, we determined focal mechanisms for a number of the best recorded earthquakes in the data set, revealing a complex pattern of deformation dominated by general arc-parallel regional tectonic compression. Most earthquakes occur in the basement rocks inferred to be the lower Tertiary Crescent formation. The sedimentary basins and the eastern part of the Olympic subduction complex are largely devoid of earthquakes. Clear association of hypocenters and focal mechanisms with previously mapped or proposed faults is difficult; however, seismicity, structure, and focal mechanisms associated with the Seattle fault zone suggest a possible high-angle mode of deformation with the north side up. We suggest that this deformation may be driven by isostatic readjustment of the Seattle basin.

Washington↗
USGS2002

Prehistoric ruptures of the Gurvan Bulag fault, Gobi Altay, Mongolia

The 1957 Gobi Altay M8.3 earthquake in southern Mongolia was associated with the simultaneous rupture of several faults, including the Gurvan Bulag reverse fault, which is located about 25 km south of the main strike-slip Bogd fault. Our study of paleoseismic excavations across the Gurvan Bulag fault suggests that the penultimate surface rupture occurred after 6.0 ka, most likely between 2.6 and 4.4 ka, and a possible earlier rupture occurred after 7.3 ka. Our interpretation of the stratigraphic relations in one of the exposures suggests that at least five earthquakes have generated surface rupture of the Gurvan Bulag fault since the abandonment of an ancient alluvial fan surface. Luminescence dating of sediment associated with this surface indicates that it formed either 26.6 ± 2.1 ka or 16.1 ± 2.0 ka. These data imply that the recurrence intervals for surface faulting on the Gurvan Bulag and Bogd faults are similar, on the order of several thousands of years, but that the penultimate surface ruptures of the two faults did not occur during the same earthquake.

Gobi Altay↗
USGS2002

Steady subsidence of Medicine Lake volcano, northern California, revealed by repeated leveling surveys

Leveling surveys of a 193‐km circuit across Medicine Lake volcano (MLV) in 1954 and 1989 show that the summit area subsided by as much as 302 ± 30 mm (−8.6 ± 0.9 mm/yr) with respect to a datum point near Bartle, California, 40 km to the southwest. This result corrects an error in the earlier analysis of the same data by Dzurisin et al. [1991] , who reported the subsidence rate as −11.1 ± 1.2 mm/yr. The subsidence pattern extends across the entire volcano, with a surface area of nearly 2000 km 2 . Two areas of localized subsidence by as much as 20 cm can be attributed to shallow normal faulting near the volcano's periphery. Surveys of an east–west traverse across Lava Beds National Monument on the north flank of the volcano in 1990 and of a 23‐km traverse across the summit area in 1999 show that subsidence continued at essentially the same rate during 1989–1999 as 1954–1989. Volcano‐wide subsidence can be explained by either a point source of volume loss (Mogi) or a contracting horizontal rectangular dislocation (sill) at a depth of 10–11 km. Volume loss rate estimates range from 0.0013 to 0.0032 km 3 /yr, depending mostly on the source depth estimate and source type. Based on first‐order quantitative considerations, we can rule out that the observed subsidence is due to volume loss from magma withdrawal, thermal contraction, or crystallizing magma at depth. Instead, we attribute the subsidence and faulting to: (1) gravitational loading of thermally weakened crust by the mass of the volcano and associated intrusive rocks, and (2) thinning of locally weakened crust by Basin and Range deformation. The measured subsidence rate exceeds long‐term estimates from drill hole data, suggesting that over long timescales, steady subsidence and episodic uplift caused by magmatic intrusions counteract each other to produce the lower net subsidence rate.

California↗
USGS2002

Preeruptive inflation and surface interferometric coherence characteristics revealed by satellite radar interferometry at Makushin Volcano, Alaska: 1993-2000

Pilot reports in January 1995 and geologic field observations from the summer of 1996 indicate that a relatively small explosive eruption of Makushin, one of the more frequently active volcanoes in the Aleutian arc of Alaska, occurred on 30 January 1995. Several independent radar interferograms that each span the time period from October 1993 to September 1995 show evidence of ???7 cm of uplift centered on the volcano's east flank, which we interpret as preeruptive inflation of a ???7-km-deep magma source (??V = 0.022 km3). Subsequent interferograms for 1995-2000, a period that included no reported eruptive activity, show no evidence of additional ground deformation. Interferometric coherence at C band is found to persist for 3 years or more on lava flow and other rocky surfaces covered with short grass and sparsely distributed tall grass and for at least 1 year on most pyroclastic deposits. On lava flow and rocky surfaces with dense tall grass and on alluvium, coherence lasts for a few months. Snow and ice surfaces lose coherence within a few days. This extended timeframe of coherence over a variety of surface materials makes C band radar interferometry an effective tool for studying volcano deformation in Alaska and other similar high-latitude regions.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Post-1906 stress recovery of the San Andreas fault system calculated from three-dimensional finite element analysis

The M = 7.8 1906 San Francisco earthquake cast a stress shadow across the San Andreas fault system, inhibiting other large earthquakes for at least 75 years. The duration of the stress shadow is a key question in San Francisco Bay area seismic hazard assessment. This study presents a three-dimensional (3-D) finite element simulation of post-1906 stress recovery. The model reproduces observed geologic slip rates on major strike-slip faults and produces surface velocity vectors comparable to geodetic measurements. Fault stressing rates calculated with the finite element model are evaluated against numbers calculated using deep dislocation slip. In the finite element model, tectonic stressing is distributed throughout the crust and upper mantle, whereas tectonic stressing calculated with dislocations is focused mostly on faults. In addition, the finite element model incorporates postseismic effects such as deep afterslip and viscoelastic relaxation in the upper mantle. More distributed stressing and postseismic effects in the finite element model lead to lower calculated tectonic stressing rates and longer stress shadow durations (17–74 years compared with 7–54 years). All models considered indicate that the 1906 stress shadow was completely erased by tectonic loading no later than 1980. However, the stress shadow still affects present-day earthquake probability. Use of stressing rate parameters calculated with the finite element model yields a 7–12% reduction in 30-year probability caused by the 1906 stress shadow as compared with calculations not incorporating interactions. The aggregate interaction-based probability on selected segments (not including the ruptured San Andreas fault) is 53–70% versus the noninteraction range of 65–77%.

California↗
USGS2002

Fault structure and kinematics of the Long Valley Caldera region, California, revealed by high-accuracy earthquake hypocenters and focal mechanism stress inversions

We have determined high-resolution hypocenters for 45,000+ earthquakes that occurred between 1980 and 2000 in the Long Valley caldera area using a double-difference earthquake location algorithm and routinely determined arrival times. The locations reveal numerous discrete fault planes in the southern caldera and adjacent Sierra Nevada block (SNB). Intracaldera faults include a series of east/west-striking right-lateral strike-slip faults beneath the caldera's south moat and a series of more northerly striking strike-slip/normal faults beneath the caldera's resurgent dome. Seismicity in the SNB south of the caldera is confined to a crustal block bounded on the west by an east-dipping oblique normal fault and on the east by the Hilton Creek fault. Two NE-striking left-lateral strike-slip faults are responsible for most seismicity within this block. To understand better the stresses driving seismicity, we performed stress inversions using focal mechanisms with 50 or more first motions. This analysis reveals that the least principal stress direction systematically rotates across the studied region, from NE to SW in the caldera's south moat to WNW-ESE in Round Valley, 25 km to the SE. Because WNW-ESE extension is characteristic of the western boundary of the Basin and Range province, caldera area stresses appear to be locally perturbed. This stress perturbation does not seem to result from magma chamber inflation but may be related to the significant (???20 km) left step in the locus of extension along the Sierra Nevada/Basin and Range province boundary. This implies that regional-scale tectonic processes are driving seismic deformation in the Long Valley caldera.

California↗
USGS2003

A revised dislocation model of interseismic deformation of the Cascadia subduction zone

CAS3D‐2, a new three‐dimensional (3‐D) dislocation model, is developed to model interseismic deformation rates at the Cascadia subduction zone. The model is considered a snapshot description of the deformation field that changes with time. The effect of northward secular motion of the central and southern Cascadia forearc sliver is subtracted to obtain the effective convergence between the subducting plate and the forearc. Horizontal deformation data, including strain rates and surface velocities from Global Positioning System (GPS) measurements, provide primary geodetic constraints, but uplift rate data from tide gauges and leveling also provide important validations for the model. A locked zone, based on the results of previous thermal models constrained by heat flow observations, is located entirely offshore beneath the continental slope. Similar to previous dislocation models, an effective zone of downdip transition from locking to full slip is used, but the slip deficit rate is assumed to decrease exponentially with downdip distance. The exponential function resolves the problem of overpredicting coastal GPS velocities and underpredicting inland velocities by previous models that used a linear downdip transition. A wide effective transition zone (ETZ) partially accounts for stress relaxation in the mantle wedge that cannot be simulated by the elastic model. The pattern of coseismic deformation is expected to be different from that of interseismic deformation at present, 300 years after the last great subduction earthquake. The downdip transition from full rupture to no slip should take place over a much narrower zone.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Undrained poroelastic response of sandstones to deviatoric stress change

Deformation of porous crustal rock through diagenesis, tectonic loading, or other processes can change pore volume and affect fluid pressure. The largest stress-induced pore pressure changes occur when fluid is trapped in pores in an “undrained” condition. We have measured the undrained poroelastic response of two sandstones to changes in mean and deviatoric stress. Pore pressure was found to respond to mean stress σ m in the usual manner: Δ p = B Δσ m ( B ranging from 0.4 to 0.7), nearly independent of the ambient deviatoric stress state. However, variations in deviatoric stress (σ d = (σ 1 − σ 3 )/2) at constant mean stress were also found to induce a reversible (elastic) pore pressure response to stress levels up to and exceeding 80% failure stress (i.e., Δ p = = const.). The coefficient η became more negative with increasing deviatoric stress level in sandstone and Ottawa sand samples. That is, η represents a dilatant response where increased deviatoric stress causes a decrease in pore pressure. The poroelastic response to deviatoric stress is explained in terms of anisotropic matrix stiffening due to closure of crack-like pore space or flattening of grain contacts at high ambient stress levels and can be important in calculations of earthquake stress transfer.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Concordant paleolatitudes for Neoproterozoic ophiolitic rocks of the Trinity Complex, Klamath Mountains, California

New paleomagnetic results from the eastern Klamath Mountains of northern California show that Neoproterozoic rocks of the Trinity ophiolitic complex and overlying Middle Devonian volcanic rocks are latitudinally concordant with cratonal North America. Combining paleomagnetic data with regional geologic and faunal evidence suggests that the Trinity Complex and related terranes of the eastern Klamath plate were linked in some fashion to the North American craton throughout that time, but that distance between them may have varied considerably. A possible model that is consistent with our paleomagnetic results and the geologic evidence is that the Trinity Complex formed and migrated parallel to paleolatitude in the basin between Laurasia and Australia–East Antarctica as the Rodinian supercontinent began to break up. It then continued to move parallel to paleolatitude at least through Middle Devonian time. Although the eastern Klamath plate served as a nucleus against which more western components of the Klamath Mountains province amalgamated, the Klamath superterrane was not accreted to North America until Early Cretaceous time.

California↗
USGS2002

Episodic thermal perturbations associated with groundwater flow: An example from Kilauea Volcano, Hawaii

Temperature measurements in deep drill holes on volcano summits or upper flanks allow a quantitative analysis of groundwater induced heat transport within the edifice. We present a new temperature-depth profile from a deep well on the summit of Kilauea Volcano, Hawaii, and analyze it in conjunction with a temperature profile measured 26 years earlier. We propose two groundwater flow models to interpret the complex temperature profiles. The first is a modified confined lateral flow model (CLFM) with a continuous flux of hydrothermal fluid. In the second, transient flow model (TFM), slow conductive cooling follows a brief, advective heating event. We carry out numerical simulations to examine the timescales associated with each of the models. Results for both models are sensitive to the initial conditions, and with realistic initial conditions it takes between 750 and 1000 simulation years for either model to match the measured temperature profiles. With somewhat hotter initial conditions, results are consistent with onset of a hydrothermal plume ∼550 years ago, coincident with initiation of caldera subsidence. We show that the TFM is consistent with other data from hydrothermal systems and laboratory experiments and perhaps is more appropriate for this highly dynamic environment. The TFM implies that volcano-hydrothermal systems may be dominated by episodic events and that thermal perturbations may persist for several thousand years after hydrothermal flow has ceased.

Hawaii↗