SEARCH · USGS Science

Results for “Geophysical Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25Linked to original sources

USGS2012

Reflections from midcrustal rocks within the Mesozoic subduction complex near the eastern Aleutian Trench

Seismic reflection data collected in 1973 by Western Geophysical Company show that highly reflective rocks make up the midcrust of the convergent margin adjacent to the eastern Aleutian Trench. These rocks form an arch that strikes obliquely across the strongly expressed northeast-southwest structural grain of exposed Mesozoic rocks. In an earlier report we proposed that the deep events mark the location of either the presently subducting plate or the top of underplated rocks. However, the short radius of curvature of this arch precludes the possibility that the deep events indicate the igneous oceanic crust. Instead, the deep reflections could be from underplated rocks that have been arched by the imbrication or underplating of strata below the reflective rocks. High-amplitude water layer and interbed multiples prevent precise connection of shallow and deep geology. For example, the Border Ranges fault, a suture between major tectonostratigraphic terranes, is not revealed on any of our seismic sections. We speculate, however, that one band of reflections that rises toward but does not reach the surface is from the Eagle River thrust fault, which separates Late Cretaceous melange from deformed turbidite sequences of the same age.

Journal of Geophysical Research Solid Earth↗

USGS2012

Geophysical constraints on Washington convergent margin structure

Gravity and magnetic maps of western Washington reveal the lateral structure and fabric of the Washington Coast Range, Puget Basin, and southern Washington Cascade Range. The magnetic and gravity maps show large amplitude positive anomalies associated with the shallow but largely buried section of Washington Coast Range mafic rocks which are separated by negative anomalies over deep sedimentary basins. The positive anomalies indicate that the Coast Range mafic basement extends farther east than previously thought, at least as far east as the longitude of Seattle. Linear and steep gravity and magnetic gradients indicate many unmapped, often buried faults in the Washington Coast Range Province. Magnetic highs are also associated with mapped batholiths in the Cascade arc. Several magnetic highs observed east of the Coast Range rocks and west of these batholiths may be associated with buried Tertiary plutons or ophiolites. Two-dimensional gravity and magnetic modeling constrained with geological and other geophysical data indicate that the Coast Range Province rocks are about 1 km thick at the coast, thickening to as much as 30 km near their postulated eastern edge. A maximum boundary on the average density of the upper 15–20 km of the rocks that compose the Coast Range Province of 2920 kg/m 3 was established by the modeling, suggesting a composition largely of basalt and gabbro with little interbedded sediments. Under these rocks may be mantle or a subduction complex composed of dense mafic, ultramafic, and sedimentary rocks like that proposed to underlie Vancouver Island. Previous gravity models of the Washington margin include lower densities for the proposed subduction complex than for Vancouver Island, suggesting a lower component of mafic and ultramafic rocks than the rocks underlying Vancouver Island. However, my Washington model requires that the proposed subduction complex be more dense than the trench sediments and, therefore, that material denser than sediments be incorporated within it. The absence of continental mantle and the modeled wedge shape of the Coast Range Province upper crust suggest that erosion of the bottom of the overriding plate by subduction processes may have occurred.

Washington↗

USGS2012

Thermal regime of the State 2-14 well, Salton Sea Scientific Drilling Project

Temperature logs were made repeatedly during breaks in drilling and both during and after flow tests in the Salton Sea Scientific Drilling Project well (State 2–14). The purpose of these logs was to assist in identifying zones of fluid loss or grain and to characterize reservoir temperatures. At the conclusion of the active phase of the project, a series of logs was begun in an attempt to establish the equilibrium temperature profile. Initially, we were able to log to depths below 3 km, but beginning in late May of 1986, it was impossible to log below about 1.8 km owing to casing failure. Our best estimates of formation temperature below 1.8 km are 305° ± 5°C at 1890 m and 355° ± 10°C at 3170 m. For the upper 1.8 km the latest temperature log (October 24, 1986), using a digital “slickline” (heat-shielded downhole recording) device, was within a few degrees Celsius of equilibrium, as confirmed by a more recent log (July 31, 1987) to a depth of ∼ 1 km. As in most other wells in the Salton Sea geothermal field, there is an impermeable, thermally conductive “cap” on the hydrothermal system; this cap extends to a depth of more than 900 m at the State 2–14 well. Thermal conductivities of 19 samples of drill cuttings from this interval were measured at room temperature. The conductivity values were corrected for in situ porosity as determined from geophysical logs and for the effects of elevated temperature. Thermal gradients decrease from about 250 mK m −1 (same as degrees Celsius per kilometer) in the upper few hundred meters to just below 200 mK m −1 near the base of the conductive cap. Using one interpretation, thermal conductivities increase with depth (mainly because of decreasing porosity), resulting in component heat flows that agree reasonably well with the mean of about 450 m W m −2 . This value agrees well with heat flow data from shallow wells within the Salton Sea geothermal field. A second interpretation, in which measured temperature coefficients of quartz- and carbonate-rich rocks are used to correct thermal conductivity, results in lower mean conductivities that are roughly constant with depth and, consequently, systematically decreasing heat flux averaging about 350 mW m −2 below 300 m. This interpretation is consistent with the inference (from fluid inclusion studies) that the rocks in this part of the field were once several tens of degrees Celsius hotter than they are now. The age of this possible disturbance is estimated at a few thousand years.

Journal of Geophysical Research Solid Earth↗

USGS2003

Three-dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications

Using P and S arrival times from 4625 local and regional earthquakes recorded at 174 seismic stations and associated geophysical investigations, this paper presents a three‐dimensional crustal and upper mantle velocity structure of southwestern China (21°–34°N, 97°–105°E). Southwestern China lies in the transition zone between the uplifted Tibetan plateau to the west and the Yangtze continental platform to the east. In the upper crust a positive velocity anomaly exists in the Sichuan Basin, whereas a large‐scale negative velocity anomaly exists in the western Sichuan Plateau, consistent with the upper crustal structure under the southern Tibetan plateau. The boundary between these two anomaly zones is the Longmen Shan Fault. The negative velocity anomalies at 50‐km depth in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with temperature and composition variations in the upper mantle. The Red River Fault is the boundary between the positive and negative velocity anomalies at 50‐km depth. The overall features of the crustal and the upper mantle structures in southwestern China are a low average velocity, large crustal thickness variations, the existence of a high‐conductivity layer in the crust or/and upper mantle, and a high heat flow value. All these features are closely related to the collision between the Indian and the Asian plates.

Journal of Geophysical Research B: Solid Earth↗

USGS2025

Crustal structure along and surrounding the Hawaiian Islands: Volcanic island construction across scarred oceanic crust

The Hawaiian Ridge is a classic example of an intraplate volcanic island chain emplaced on oceanic lithosphere. We seek to constrain both the deformation from island loading around the Hawaiian Ridge and the influence of the oceanic lithosphere, including the Molokaʻi fracture zone (MFZ), on off-axis volcanic emplacement. To examine these processes, we conducted a marine geophysical experiment in 2018 that included the acquisition of eight multichannel seismic reflection lines and used a 6,600 in 3 tuned air gun array and an ultra-long hydrophone streamer cable towed behind R/V Marcus G. Langseth across and around the southern Hawaiian Ridge. We image both the top of igneous oceanic crust and the Moho, and we observe significant variations in regional crustal structure and thickness variations, Moho characteristics, the locations of the buried MFZ, and the flexure of the Pacific oceanic lithosphere under the Hawaiian Ridge. We observe up to ∼4.5 km of deflection of the igneous oceanic crust in response to the volcanic load with sediment thickness increasing to ∼3–3.4 km near the ridge but not correlating with the deflection. A systematic difference in igneous oceanic crustal thickness is found north (average ∼5.2 km) and south (average ∼5.9 km), of the MFZ that also correspond to changes in Moho characteristics. The MFZ itself is associated with the largest crustal thickness variations (∼3.7–7.6 km). Magma ascent through these different crustal regions may account for some variations of magmatic flux to the surface along the Hawaiian Ridge.

Hawaii↗

USGS2001

Rock types of South Pole-Aitken basin and extent of basaltic volcanism

The enormous pre-Nectarian South Pole-Aitken (SPA) basin represents a geophysically and compositionally unique region on the Moon. We present and analyze the mineralogical diversity across this basin and discuss the implications for basin evolution. Rock types are derived from Clementine multispectral data based on diagnostic characteristics of ferrous absorptions in fresh materials. Individual areas are characterized as noritic (dominated by low-Ca pyroxene), gabbroic/basaltic (dominated by high-Ca pyroxene), feldspathic (<3-6% FeO), and olivine-gabbro (dominated by high-Ca pyroxene and olivine). The anorthositic crust has effectively been removed from the interior of the basin. The style of volcanism within the basin extends over several 100 Myr and includes mare basalt and pyroclastic deposits. Several areas of ancient (pre-Orientale) volcanism, or cryptomaria, have also been identified. The nonmare mafic lithology that occurs across the basin is shown to be noritic in composition and is pervasive laterally and vertically. We interpret this to represent impact melt/breccia deposits derived from the lower crust. A few localized areas are identified within the basin that contain more diverse lithologies (gabbro, olivine-gabbro), some of which may represent material from the deepest part of the lower crust and perhaps uppermost mantle involved in the SPA event.

Journal of Geophysical Research E: Planets↗

USGS2007

Hydrogeologic controls imposed by mechanical stratigraphy in layered rocks of the Chateauguay River Basin, a U.S.-Canada transborder aquifer

[1] The Châteauguay River Basin delineates a transborder watershed with roughly half of its surface area located in northern New York State and half in southern Québec Province, Canada. As part of a multidisciplinary study designed to characterize the hydrogeologic properties of this basin, geophysical logs were obtained in 12 wells strategically located to penetrate the four major sedimentary rock formations that constitute the regional aquifers. The layered rocks were classified according to their elastic properties into three primary units: soft sandstone, hard sandstone, and dolostone. Downhole measurements were analyzed to identify fracture patterns associated with each unit and to evaluate their role in controlling groundwater flow. Fracture networks are composed of orthogonal sets of laterally extensive, subhorizontal bedding plane partings and bed-delimited, subvertical joints with spacings that are consistent with rock mechanics principles and stress models. The vertical distribution of transmissive zones is confined to a few select bedding plane fractures, with soft sandstone having the fewest (one per 70-m depth) and hard sandstone the most (five per 70-m depth). Bed-normal permeability is examined using a probabilistic model that considers the lengths of flow paths winding along joints and bedding plane fractures. Soft sandstone has the smallest bed-normal permeability primarily because of its wide, geomechanically undersaturated joint spacing. Results indicate that the three formations have similar values of bulk transmissivity, within roughly an order of magnitude, but that each rock unit has its own unique system of groundwater flow paths that constitute that transmissivity.

Chateauguay River Basin↗

USGS2012

Elastic moduli of rock glasses under pressure to 8 kilobars and geophysical implications

Shear and longitudinal velocities were measured by the ultrasonic phase comparison method as a function of pressure to 8 kbar on synthetic glasses of basalt, andesite, rhyolite, and quartz composition and on natural obsidian. Velocities of most of the glasses decrease anomalously with pressure, but increasingly more-normal behavior occurs with decrease in SiO 2 content. The pressure derivatives of rigidity and bulk modulus increase linearly, from −3.39 to −0.26 and from −5.91 to +2.09, respectively, with decrease in SiO 2 content from 100 to 49%. The change from negative to positive in the pressure derivatives of both moduli and observed at Poisson's ratio of about 0.25 is consistent with the Smyth model for the anomalous elastic behavior of glass. If the temperature in the upper mantle is about 1500°C, tholeiitic basalt would be molten in accordance with the partial melt explanation for the low-velocity zone; at 1300°C and below, basalt would be in the glassy state, especially if more felsic than tholeiite. At a temperature of 1370°C and at 30 kbar, reasonable values for the upper mantle at 100 km depth, the basalt glass of this study would have a viscosity of about 10 13 P. On the basis of the theory of viscoelasticity the glass would support shear wave propagation at frequencies above 0.1 Hz. Under such conditions of PT, 10 to 30% basalt glass in a matrix either of eclogite or olivine would reduce the seismic velocities by 3 to 9% and could also account for the values observed in the low-velocity zone.

Journal of Geophysical Research Solid Earth↗

USGS2012

Fracture characterization by means of attenuation and generation of tube waves in fractured crystalline rock at Mirror Lake, New Hampshire

Results are presented from experiments carried out in conjunction with the U. S. Geological Survey at the Hubbard Brook Experimental Forest near Mirror Lake, New Hampshire. The study focuses on our ability to obtain orientation and transmissivity estimates of naturally occurring fractures. The collected data set includes a four-offset hydrophone vertical seismic profile, full waveform acoustic logs at 5, 15, and 34 kHz, borehole televiewer, temperature, resistivity, and self-potential logs, and borehole-to-borehole pump test data. Borehole televiewer and other geophysical logs indicate that permeable fractures intersect the Mirror Lake boreholes at numerous depths, but less than half of these fractures appear to have significant permeability beyond the annulus of drilling disturbance on the basis of acoustic waveform log analysis. The vertical seismic profiling (VSP) data indicate a single major permeable fracture near a depth of 44 m, corresponding to one of the most permeable fractures identified in the acoustic waveform log analysis. VSP data also indicate a somewhat less permeable fracture at 220 m and possible fractures at depths of 103 and 135 m; all correspond to major permeable fractures in the acoustic waveform data set. Pump test data confirm the presence of a hydraulic connection between the Mirror Lake boreholes through a shallow dipping zone of permeability at 44 m in depth. Effective fracture apertures calculated from modeled transmissivities correspond to those estimated for the largest fractures indicated on acoustic waveform logs but are over an order of magnitude larger than effective apertures calculated from tube waves in the VSP data set. This discrepancy is attributed to the effect of fracture stiffness. A new model is presented to account for the mechanical strength of asperities in resisting fracture closure during the passage of seismic waves during the generation of VSPs.

New Hampshire↗

USGS2013

Bayesian inversion of data from effusive volcanic eruptions using physics-based models: Application to Mount St. Helens 2004--2008

Physics-based models of volcanic eruptions can directly link magmatic processes with diverse, time-varying geophysical observations, and when used in an inverse procedure make it possible to bring all available information to bear on estimating properties of the volcanic system. We develop a technique for inverting geodetic, extrusive flux, and other types of data using a physics-based model of an effusive silicic volcanic eruption to estimate the geometry, pressure, depth, and volatile content of a magma chamber, and properties of the conduit linking the chamber to the surface. A Bayesian inverse formulation makes it possible to easily incorporate independent information into the inversion, such as petrologic estimates of melt water content, and yields probabilistic estimates for model parameters and other properties of the volcano. Probability distributions are sampled using a Markov-Chain Monte Carlo algorithm. We apply the technique using GPS and extrusion data from the 2004–2008 eruption of Mount St. Helens. In contrast to more traditional inversions such as those involving geodetic data alone in combination with kinematic forward models, this technique is able to provide constraint on properties of the magma, including its volatile content, and on the absolute volume and pressure of the magma chamber. Results suggest a large chamber of >40 km 3 with a centroid depth of 11–18 km and a dissolved water content at the top of the chamber of 2.6–4.9 wt%.

Washington↗

USGS2017

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel↗

USGS2020

Combined seismic and geodetic analysis before, during and after the 2018 Mt. Etna eruption

In December 2018, Etna volcano experienced one of the largest episodes of unrest since the installation of geophysical monitoring networks in 1970. The unrest culminated in a short eruption with a small volume of lava erupted, a significant seismic crisis and deformation of the entire volcanic edifice of magnitude never recorded before at Mount Etna. Here we describe the evolution of the 2018 eruptive cycle from the analysis of seismic and geodetic data collected in the months preceding, during, and following the intrusion. We model the space-time evolution of high-rate deformation data starting from the active source previously identified from deformation data and the propagation of seismicity in a 3-D velocity model. The intrusion model suggests emplacement of two dikes: a smaller dike located beneath the eruptive fissure and a second, deeper dike between 1 and 5 km below sea level that opened ~2 m. The rise and eruption of magma from the shallower dike did not interrupt the pressurization of a long-lasting deeper reservoir (~6 km) that induced continuous inflation and intense deformation of the eastern flank. Shortly after the intrusion, on 26 December 2018, a M L 4.8 earthquake occurred near Pisano, destroying buildings and roads in two villages. We propose a time-dependent intrusion model that supports the hypothesis of the inflation inducing flank deformation and that this process has been active since September 2018.

Mt. Etna↗

USGS2019

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

USGS2012

Grain-size segregation and levee formation in geophysical mass flows

Data from large-scale debris-flow experiments are combined with modeling of particle-size segregation to explain the formation of lateral levees enriched in coarse grains. The experimental flows consisted of 10 m 3 of water-saturated sand and gravel, which traveled ∼80 m down a steeply inclined flume before forming an elongated leveed deposit 10 m long on a nearly horizontal runout surface. We measured the surface velocity field and observed the sequence of deposition by seeding tracers onto the flow surface and tracking them in video footage. Levees formed by progressive downslope accretion approximately 3.5 m behind the flow front, which advanced steadily at ∼2 m s −1 during most of the runout. Segregation was measured by placing ∼600 coarse tracer pebbles on the bed, which, when entrained into the flow, segregated upwards at ∼6–7.5 cm s −1 . When excavated from the deposit these were distributed in a horseshoe-shaped pattern that became increasingly elevated closer to the deposit termination. Although there was clear evidence for inverse grading during the flow, transect sampling revealed that the resulting leveed deposit was strongly graded laterally, with only weak vertical grading. We construct an empirical, three-dimensional velocity field resembling the experimental observations, and use this with a particle-size segregation model to predict the segregation and transport of material through the flow. We infer that coarse material segregates to the flow surface and is transported to the flow front by shear. Within the flow head, coarse material is overridden, then recirculates in spiral trajectories due to size-segregation, before being advected to the flow edges and deposited to form coarse-particle-enriched levees.

Journal of Geophysical Research↗

USGS2012

The hydrothermal system at Newberry Volcano, Oregon

Results of recent geological and geophysical studies at Newberry Volcano have been incorporated into conceptual and numerical models of a magma-based hydrothermal system. Numerical simulations begin with emplacement of a small magma body, the presumed source of silicic eruptions at Newberry that began about 10,000 B.P., into a thermal regime representing 100,000 years of cooling of a large underlying intrusion. Simulated flow patterns and thermal histories for three sets of hypothetical permeability values are compatible with data from four geothermal drill holes on the volcano. Meteoric recharge cools the caldera-fill deposits, but thermal water moving up a central conduit representing a permeable volcanic vent produces temperatures close to those observed in drill holes within the caldera. Meteoric recharge from the caldera moves down the flanks and creates a near-isothermal zone that extends several hundred meters below the water table, producing temperature profiles similar to those observed in drill holes on the flanks. The temperatures observed in drill holes on the flanks are not influenced by the postulated Holocene magma body. The elevated temperature gradients measured in the lower portions of these holes may be related to the cumulative effect of older intrusions. The models also indicate that meteoric recharge to the deep hydrothermal system probably originates within or near the caldera. Relatively low fluid velocities at depth suggest that at least a significant fraction of the thermal fluid may be very old.

Journal of Geophysical Research Solid Earth↗

USGS2010

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California↗

USGS2014

A continuous record of intereruption velocity change at Mount St. Helens from coda wave interferometry

In September 2004, Mount St. Helens volcano erupted after nearly 18 years of quiescence. However, it is unclear from the limited geophysical observations when or if the magma chamber replenished following the 1980–1986 eruptions in the years before the 2004–2008 extrusive eruption. We use coda wave interferometry with repeating earthquakes to measure small changes in the velocity structure of Mount St. Helens volcano that might indicate magmatic intrusion. By combining observations of relative velocity changes from many closely located earthquake sources, we solve for a continuous function of velocity changes with time. We find that seasonal effects dominate the relative velocity changes. Seismicity rates and repeating earthquake occurrence also vary seasonally; therefore, velocity changes and seismicity are likely modulated by snow loading, fluid saturation, and/or changes in groundwater level. We estimate hydrologic effects impart stress changes on the order of tens of kilopascals within the upper 4 km, resulting in annual velocity variations of 0.5 to 1%. The largest nonseasonal change is a decrease in velocity at the time of the deep M w  = 6.8 Nisqually earthquake. We find no systematic velocity changes during the most likely times of intrusions, consistent with a lack of observable surface deformation. We conclude that if replenishing intrusions occurred, they did not alter seismic velocities where this technique is sensitive due to either their small size or the finite compressibility of the magma chamber. We interpret the observed velocity changes and shallow seasonal seismicity as a response to small stress changes in a shallow, pressurized system.

Journal of Geophysical Research B: Solid Earth↗

USGS2017

Constraining the magmatic system at Mount St. Helens (2004–2008) using Bayesian inversion with physics-based models including gas escape and crystallization

Physics-based models of volcanic eruptions track conduit processes as functions of depth and time. When used in inversions, these models permit integration of diverse geological and geophysical data sets to constrain important parameters of magmatic systems. We develop a 1-D steady state conduit model for effusive eruptions including equilibrium crystallization and gas transport through the conduit and compare with the quasi-steady dome growth phase of Mount St. Helens in 2005. Viscosity increase resulting from pressure-dependent crystallization leads to a natural transition from viscous flow to frictional sliding on the conduit margin. Erupted mass flux depends strongly on wall rock and magma permeabilities due to their impact on magma density. Including both lateral and vertical gas transport reveals competing effects that produce nonmonotonic behavior in the mass flux when increasing magma permeability. Using this physics-based model in a Bayesian inversion, we link data sets from Mount St. Helens such as extrusion flux and earthquake depths with petrological data to estimate unknown model parameters, including magma chamber pressure and water content, magma permeability constants, conduit radius, and friction along the conduit walls. Even with this relatively simple model and limited data, we obtain improved constraints on important model parameters. We find that the magma chamber had low (<5wt%) total volatiles and that the magma permeability scale is well constrained at ~10-11.4 m2 to reproduce observed dome rock porosities. Compared with previous results, higher magma overpressure and lower wall friction are required to compensate for increased viscous resistance while keeping extrusion rate at the observed value.

Washington↗