Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research B: 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 37 records · Page 2Linked to original sources

The strength and rheology of methane clathrate hydrate

Methane clathrate hydrate (structure I) is found to be very strong, based on laboratory triaxial deformation experiments we have carried out on samples of synthetic, high-purity, polycrystalline material. Samples were deformed in compressional creep tests (i.e., constant applied stress, ??), at conditions of confining pressure P = 50 and 100 MPa, strain rate 4.5 ?? 10-8 ??? ?? ??? 4.3 ?? 10-4 s-1, temperature 260 ??? T ??? 287 K, and internal methane pressure 10 ??? PCH4 ??? 15 MPa. At steady state, typically reached in a few percent strain, methane hydrate exhibited strength that was far higher than expected on the basis of published work. In terms of the standard high-temperature creep law, ?? = A??ne-(E*+PV*)/RT the rheology is described by the constants A = 108.55 MPa-n s-1, n = 2.2, E* = 90,000 J mol-1, and V* = 19 cm3 mol-1. For comparison at temperatures just below the ice point, methane hydrate at a given strain rate is over 20 times stronger than ice, and the contrast increases at lower temperatures. The possible occurrence of syntectonic dissociation of methane hydrate to methane plus free water in these experiments suggests that the high strength measured here may be only a lower bound. On Earth, high strength in hydrate-bearing formations implies higher energy release upon decomposition and subsequent failure. In the outer solar system, if Titan has a 100-km-thick near-surface layer of high-strength, low-thermal conductivity methane hydrate as has been suggested, its interior is likely to be considerably warmer than previously expected.

Journal of Geophysical Research B: Solid Earth

A dislocation model of strain accumulation and release at a subduction zone

Strain accumulation and release at a subduction zone are attributed to stick slip on the main thrust zone and steady aseismic slip on the remainder of the plate interface. This process can be described as a superposition of steady state subduction and a repetitive cycle of slip on the main thrust zone, consisting of steady normal slip at the plate convergence rate plus occasional thrust events that recover the accumulated normal slip. Because steady state subduction does not contribute to the deformation at the free surface, deformation observed there is completely equivalent to that produced by the slip cycle alone. The response to that slip is simply the response of a particular earth model to embedded dislocations. For a purely elastic earth model, the deformation cycle consists of a coseismic offset followed by a linear‐in‐time recovery to the initial value during the interval between earthquakes. For an elastic‐viscoelastic earth model (elastic lithosphere over a viscoelastic asthenosphere), the postearthquake recovery is not linear in time. Records of local uplift as a function of time indicate that the long‐term postseismic recovery is approximately linear, suggesting that elastic earth models are adequate to describe the deformation cycle. However, the deformation predicted for a simple elastic half‐space earth model does not reproduce the deformation observed along the subduction zones in Japan at all well if stick slip is restricted to the main thrust zone. As recognized earlier by Shimazaki, Seno, and Kato, the uplift profiles could be explained if stick slip were postulated to extend along the plate interface beyond the main thrust zone to a depth of perhaps 100 km, but independent evidence suggests that stick slip at such depths is unlikely.

Journal of Geophysical Research B: Solid Earth

Evaluation of temporally correlated noise in global navigation satellite system time series: Geodetic monument performance

Estimates of background noise of Global Positioning System‐derived time series of positions for 740 sites in the western United States are examined. These data consist of daily epochs of three components of displacements that are at least 9.75 years long within the interval between 2000 and 2018. We find that these time series have significant temporal correlations that could be represented as a combination of white, flicker, random‐walk, and band‐pass filtered noise. From this noise model, two other metrics are computed: the root‐mean‐square of seasonal noise, that is, the integrated power spectrum between 0.5 and 2 cycles per year, and the standard error in position rate for a 10‐year‐long time series. These two metrics are used to evaluate potential correlations with different geographic regions and with different methods of construction of monuments used to attach the Global Positioning System antenna to the Earth's surface. The sites with the lowest noise, both in terms of rate error and seasonal root‐mean‐square, are located in semiarid regions east of the rain shadow provided by the Cascade and Sierra Nevada mountain ranges. In addition, according to statistical rank tests, monuments known as drilled‐braced monuments perform 30% to 50% better than other monument types (buildings, boreholes, piers, etc.) in terms of having smaller rate errors and lower seasonal noise.

Journal of Geophysical Research B: Solid Earth

Slip rates and spatially variable creep on faults of the northern San Andreas system inferred through Bayesian inversion of Global Positioning System data

Fault creep, depending on its rate and spatial extent, is thought to reduce earthquake hazard by releasing tectonic strain aseismically. We use Bayesian inversion and a newly expanded GPS data set to infer the deep slip rates below assigned locking depths on the San Andreas, Maacama, and Bartlett Springs Faults of Northern California and, for the latter two, the spatially variable interseismic creep rate above the locking depth. We estimate deep slip rates of 21.5 ± 0.5, 13.1 ± 0.8, and 7.5 ± 0.7 mm/yr below 16 km, 9 km, and 13 km on the San Andreas, Maacama, and Bartlett Springs Faults, respectively. We infer that on average the Bartlett Springs fault creeps from the Earth's surface to 13 km depth, and below 5 km the creep rate approaches the deep slip rate. This implies that microseismicity may extend below the locking depth; however, we cannot rule out the presence of locked patches in the seismogenic zone that could generate moderate earthquakes. Our estimated Maacama creep rate, while comparable to the inferred deep slip rate at the Earth's surface, decreases with depth, implying a slip deficit exists. The Maacama deep slip rate estimate, 13.1 mm/yr, exceeds long-term geologic slip rate estimates, perhaps due to distributed off-fault strain or the presence of multiple active fault strands. While our creep rate estimates are relatively insensitive to choice of model locking depth, insufficient independent information regarding locking depths is a source of epistemic uncertainty that impacts deep slip rate estimates.

California

Plume-driven plumbing and crustal formation in Iceland

Through combination of surface wave and body wave constraints we derive a three-dimensional (3-D) crustal S velocity model and Moho map for Iceland. It reveals a vast plumbing system feeding mantle plume melt into upper crustal magma chambers where crustal formation takes place. The method is based on the partitioned waveform inversion to which we add additional observations. Love waves from six local events recorded on the HOTSPOT- SIL networks are fitted, S n travel times from the same events measured, previous observations of crustal thickness are added, and all three sets of constraints simultaneously inverted for our 3-D model. In the upper crust (0–15 km) an elongated low-velocity region extends along the length of the Northern, Eastern and Western Neovolcanic Zones. The lowest velocities (−7%) are found at 5–10 km below the two most active volcanic complexes: Hekla and Bárdarbunga-Grímsvötn. In the lower crust (>15 km) the low-velocity region can be represented as a vertical cylinder beneath central Iceland. The low-velocity structure is interpreted as the thermal halo of pipe work which connects the region of melt generation in the uppermost mantle beneath central Iceland to active volcanoes along the neovolcanic zones. Crustal thickness in Iceland varies from 15–20 km beneath the Reykjanes Peninsula, Krafla and the extinct Snæfellsnes rift zone, to 46 km beneath central Iceland. The average crustal thickness is 29 km. The variations in thickness can be explained in terms of the temporal variation in plume productivity over the last ∼20 Myr, the Snæfellsnes rift zone being active during a minimum in plume productivity. Variations in crustal thickness do not depart significantly from an isostatically predicted crustal thickness. The best fit linear isostatic relation implies an average density jump of 4% across the Moho. Rare earth element inversions of basalt compositions on Iceland suggest a melt thickness (i.e., crustal thickness) of 15–20 km, given passive upwelling. The observed crustal thickness of up to 46 km implies active fluxing of source material through the melt zone by the mantle plume at up to 3 times the passive rate.

Journal of Geophysical Research B: Solid Earth

Volcanic plume height measured by seismic waves based on a mechanical model

In August 2008 an unmonitored, largely unstudied Aleutian volcano, Kasatochi, erupted catastrophically. Here we use seismic data to infer the height of large eruptive columns such as those of Kasatochi based on a combination of existing fluid and solid mechanical models. In so doing, we propose a connection between a common, observable, short-period seismic wave amplitude to the physics of an eruptive column. To construct a combined model, we estimate the mass ejection rate of material from the vent on the basis of the plume height, assuming that the height is controlled by thermal buoyancy for a continuous plume. Using the estimated mass ejection rate, we then derive the equivalent vertical force on the Earth through a momentum balance. Finally, we calculate the far-field surface waves resulting from the vertical force. The model performs well for recent eruptions of Kasatochi and Augustine volcanoes if v , the velocity of material exiting the vent, is 120-230 m s -1 . The consistency between the seismically inferred and measured plume heights indicates that in these cases the far-field ~1 s seismic energy radiated by fluctuating flow in the volcanic jet during the eruption is a useful indicator of overall mass ejection rates. Thus, use of the model holds promise for characterizing eruptions and evaluating ash hazards to aircraft in real time on the basis of far-field short-period seismic data. This study emphasizes the need for better measurements of eruptive plume heights and a more detailed understanding of the full spectrum of seismic energy radiated coeruptively.

Journal of Geophysical Research B: Solid Earth

Triggering and modulation of geyser eruptions in Yellowstone National Park by earthquakes, earth tides, and weather

We analyze intervals between eruptions (IBEs) data acquired between 2001 and 2011 at Daisy and Old Faithful geysers in Yellowstone National Park. We focus our statistical analysis on the response of these geysers to stress perturbations from within the solid earth (earthquakes and earth tides) and from weather (air pressure and temperature, precipitation, and wind). We conclude that (1) the IBEs of these geysers are insensitive to periodic stresses induced by solid earth tides and barometric pressure variations; (2) Daisy (pool geyser) IBEs lengthen by evaporation and heat loss in response to large wind storms and cold air; and (3) Old Faithful (cone geyser) IBEs are not modulated by air temperature and pressure variations, wind, and precipitation, suggesting that the subsurface water column is decoupled from the atmosphere. Dynamic stress changes of 0.1−0.2 MPa resulting from the 2002 M-7.9 Denali, Alaska, earthquake surface waves caused a statistically significant shortening of Daisy geyser's IBEs. Stresses induced by other large global earthquakes during the study period were at least an order of magnitude smaller. In contrast, dynamic stresses of >0.5 MPa from three large regional earthquakes in 1959, 1975, and 1983 caused lengthening of Old Faithful's IBEs. We infer that most subannual geyser IBE variability is dominated by internal processes and interaction with other geysers. The results of this study provide quantitative bounds on the sensitivity of hydrothermal systems to external stress perturbations and have implications for studying the triggering and modulation of volcanic eruptions by external forces.

Wyoming

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

Seismic source mechanism of degassing bursts at Kilauea volcano, Hawaii: Results from waveform inversion in the 10–50 s band

The current (March 2008 to February 2009) summit eruptive activity at Kilauea Volcano is characterized by explosive degassing bursts accompanied by very long period (VLP) seismic signals. We model the source mechanisms of VLP signals in the 10–50 s band using data recorded for 15 bursts with a 10‐station broadband network deployed in the summit caldera. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. The VLP signals associated with the bursts originate in a source region ∼1 km below the eastern perimeter of Halemaumau pit crater. The observed waveforms are well explained by the combination of a volumetric component and a vertical single force component. For the volumetric component, several source geometries are obtained which equally explain the observed waveforms. These geometries include (1) a pipe dipping 64° to the northeast; (2) two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east; (3) a pipe dipping 58° to the northeast, intersecting a crack dipping 48° to the west–southwest; and (4) a pipe dipping 57° to the northeast, intersecting a pipe dipping 58° to the west–southwest. Using the dual‐crack model as reference, the largest volume change obtained among the 15 bursts is ∼24,400 m 3 , and the maximum amplitude (peak to peak) of the force is ∼20 GN. Each burst is marked by a similar sequence of deflation and inflation, trailed by decaying oscillations of the volumetric source. The vertical force is initially upward, synchronous with source deflation, then downward, synchronous with source reinflation, followed by oscillations with polarity opposite to the volumetric oscillations. This combination of force and volume change is attributed to pressure and momentum changes induced during a fluid dynamic source mechanism involving the ascent, expansion, and burst of a large slug of gas within the upper ∼150 m of the magma conduit. As the slug expands upon approach to the surface and more liquid becomes wall supported by viscous shear forces, the pressure below the slug decreases, inducing conduit deflation and an upward force on the Earth. The final rapid slug expansion and burst stimulate VLP and LP oscillations of the conduit system, which slowly decay due to viscous dissipation and elastic radiation. Consideration of the fluid dynamic arguments leads us to prefer the dual‐crack VLP source model as it is the only candidate model capable of producing plausible values of length scales and pressure changes. The magnitudes of the vertical forces observed in the 15 bursts appear consistent with slug masses of 10 4 to 10 6 kg.

Hawaii

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth

Strain accumulation across the Prince William Sound asperity, Southcentral Alaska

The surface velocities predicted by the conventional subduction model are compared to velocities measured in a GPS array (surveyed in 1993, 1995, 1997, 2000, and 2004) spanning the Prince William Sound asperity. The observed velocities in the comparison have been corrected to remove the contributions from postseismic (1964 Alaska earthquake) mantle relaxation. Except at the most seaward monument (located on Middleton Island at the seaward edge of the continental shelf, just 50 km landward of the deformation front in the Aleutian Trench), the corrected velocities qualitatively agree with those predicted by an improved, two-dimensional, back slip, subduction model in which the locked megathrust coincides with the plate interface identified by seismic refraction surveys, and the back slip rate is equal to the plate convergence rate. A better fit to the corrected velocities is furnished by either a back slip rate 20% greater than the plate convergence rate or a 30% shallower megathrust. The shallow megathrust in the latter fit may be an artifact of the uniform half-space Earth model used in the inversion. Backslip at the plate convergence rate on the megathrust mapped by refraction surveys would fit the data as well if the rigidity of the underthrust plate was twice that of the overlying plate, a rigidity contrast higher than expected. The anomalous motion at Middleton Island is attributed to continuous slip at near the plate convergence rate on a postulated, listric fault that splays off the megathrust at depth of about 12 km and outcrops on the continental slope south-southeast of Middleton Island.

Journal of Geophysical Research B: Solid Earth

The Steens Mountain (Oregon) geomagnetic polarity transition: 3. Its regional significance

Study of the variations of direction and intensity of the geomagnetic field as recorded by the Miocene lava flows on Steens Mountain, southeastern Oregon, has resulted in a detailed description of total field behavior during a reversal in polarity. In addition to information about the polarity reversal itself, the detailed paleomagnetic record includes several thousand years of geomagnetic history preceding and following the polarity transition at 15.5 Ma. In order to test the feasibility of using this record as a means of correlation in this part of the western United States, comparisons are made of reconnaissance and previously published paleomagnetic records obtained from what has been thought to be the Steens Basalt or rocks of equivalent age. Despite the fact that many of these earlier studies were not done in detail and were not intended for correlation purposes, convincing similarities among some of the records are evident. The Steens Basalt paleomagnetic record does, indeed, have potential as a correlation tool during this time of widespread basaltic volcanism. Additionally, paleomagnetic data from flows that were sampled in detail yield a middle Miocene paleomagnetic pole at 88.3°N, 209.0° (α 95 = 6.3°) for the High Lava Plains of Oregon. This pole position is statistically indistinguishable from the earth's rotational axis and implies that no tectonic rotation of this region has occurred since these lava flows were erupted. Data from selected sites within the coeval part of the Columbia River Basalt Group yield a paleomagnetic pole at 88.7°N, 171.6°E (α 95 = 4.0°). The Columbia River Basalt Group pole is statistically indistinguishable from either the rotational axis or from the High Lava Plains pole. These findings indicate no post‐20 Ma differential rotation between south‐eastern Washington and south‐central Oregon, in contrast to previous interpretations.

Journal of Geophysical Research B: Solid Earth

Crustal deformation associated with glacial fluctuations in the eastern Chugach Mountains, Alaska

The changes of the solid Earth in south central Alaska in response to two major glacial fluctuations on different temporal and spatial scales have been estimated and we evaluated their influence on the stress state and ongoing tectonic deformation of the region. During the recent (1993–1995) Bering Glacier surge, a large transfer of ice from the Bagley Ice Field to the Bering Glacier terminus region occurred. We estimated the elastic displacement of the solid Earth due to ice mass redistribution from Global Positioning System (GPS) measurements at sites near the surging glacier. We can account for these displacements by transfer of an ice volume of about 14 km 3 from the surge reservoir area to the terminus region. We examined the background seismicity (M L ≥ 2.5) before, during, and after the surge. We found that the occurrence of small earthquakes (M L ≤ 4.0) in the surge reservoir region increased during the surge time interval possibly in response to a decrease in ice mass. This suggests that a small decrease in the vertical stress, σ 3 , could be enough to modulate the occurrence of small, shallow earthquakes in this dominantly thrust fault setting. During this century the southern Alaska coastal glaciers have been undergoing an overall decrease in volume. Based on our compilation of changes in the extent and thickness of the coastal glaciers between the Malaspina and Bering, we calculated surface displacements due to the Earth's viscoelastic response to annual thinning and to the cumulative retreat over the last 100 years. The uplift of the region due to an average annual thinning rate of 1–6 m/yr in the ablation region is 1–12 mm/yr. For our reference model with a viscosity of 5×10 19 Pa s for depths between ≈ 40 and 200 km the total viscoelastic response due to the retreat over the last century may be as much as a couple of meters within the coastal ablation zone near Icy Bay. The maximum decrease in σv between 0 and 10 km was ≈ 1.0 MPa, which is significant in relation to the stress drops in recent earthquakes (≈ 2 to 10 MPa) but small in relation to the estimated tectonic stress magnitude. Therefore the occurrence of an earthquake such as the St. Elias (1979, M S = 7.2) may have been advanced in time; however, most of the ongoing stress accumulation would be primarily due to tectonic forces.

Journal of Geophysical Research B: Solid Earth

Temporal stress changes caused by earthquakes: A review

Earthquakes can change the stress field in the Earth’s lithosphere as they relieve and redistribute stress. Earthquake-induced stress changes have been observed as temporal rotations of the principal stress axes following major earthquakes in a variety of tectonic settings. The stress changes due to the 2011 Mw9.0 Tohoku-Oki, Japan, earthquake were particularly well documented. Earthquake stress rotations can inform our understanding of earthquake physics, most notably addressing the long-standing problem of whether the Earth’s crust at plate boundaries is “strong” or “weak.” Many of the observed stress rotations, including that due to the Tohoku-Oki earthquake, indicate near-complete stress drop in the mainshock. This implies low background differential stress, on the order of earthquake stress drop, supporting the weak crust model. Earthquake stress rotations can also be used to address other important geophysical questions, such as the level of crustal stress heterogeneity and the mechanisms of postseismic stress reloading. The quantitative interpretation of stress rotations is evolving from those based on simple analytical methods to those based on more sophisticated numerical modeling that can capture the spatial-temporal complexity of the earthquake stress changes.

Journal of Geophysical Research B: Solid Earth

Climatic impact of glacial cycle polar motion: Coupled oscillations of ice sheet mass and rotation pole position

Precessional motion of Earth's rotation axis relative to its orbit is a well-known source of long-period climatic variation. It is less well appreciated that growth and decay of polar ice sheets perturb the symmetry of the global mass distribution enough that the geographic location of the rotation axis will change by at least 15 km and possibly as much as 100 km during a single glacial cycle. This motion of the pole will change the seasonal and latitudinal pattern of temperatures. We present calculations, based on a diurnal average energy balance, which compare the summer and winter temperature anomalies due to a 1° decrease in obliquity with those due to a 1° motion of the rotation pole toward Hudson Bay. Both effects result in peak temperature perturbations of about 1° Celsius. The obliquity change primarily influences the amplitude of the seasonal cycle, while the polar motion primarily changes the annual mean temperatures. The polar motion induced temperature anomaly is such that it will act as a powerful negative feedback on ice sheet growth. We also explore the evolution of the coupled system composed of ice sheet mass and pole position. Oscillatory solutions result from the conflicting constraints of rotational and thermal stability. A positive mass anomaly on an otherwise featureless Earth is in rotational equilibrium only at the poles or the equator. The two polar equilibria are rotationally unstable, and the equatorial equilibrium, though rotationally stable, is thermally unstable. We find that with a plausible choice for the strength of coupling between the thermal and rotational systems, relatively modest external forcing can produce significant response at periods of 10 4 –10 6 years, but it strongly attenuates polar motion at longer periods. We suggest that these coupled oscillations may contribute to the observed dominance of 100 kyr glacial cycles since the mid-Pleistocene and will tend to stabilize geographic patterns that are suitable to glaciations.

Journal of Geophysical Research B: Solid Earth

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

Pore fluid pressure, apparent friction, and Coulomb failure

Many recent studies of stress-triggered seismicity rely on a fault failure model with a single free parameter, the apparent coefficient of friction, presumed to be a material constant with possible values 0 ≤ μ′ ≤ 1. These studies may present a misleading view of fault strength and the role of pore fluid pressure in earthquake failure. The parameter μ′ is intended to incorporate the effects of both friction and pore pressure, but is a material constant only if changes in pore fluid pressure induced by changes in stress are proportional to the normal stress change across the potential failure plane. Although specific models of fault zones permit such a relation, neither is it known that fault zones within the Earth behave this way, nor is this behavior expected in all cases. In contrast, for an isotropic homogeneous poroelastic model the pore pressure changes are proportional to changes in mean stress, μ′ is not a material constant, and −∞ ≤ μ′ ≤ +∞. Analysis of the change in Coulomb failure stress for tectonically loaded reverse and strike-slip faults shows considerable differences between these two pore pressure models, suggesting that such models might be distinguished from one another using observations of triggered seismicity (e.g., aftershocks). We conclude that using the constant apparent friction model exclusively in studies of Coulomb failure stress is unwise and could lead to significant errors in estimated stress change and seismic hazard.

Journal of Geophysical Research B: Solid Earth

Mass addition at Mount St. Helens, Washington, inferred from repeated gravity surveys

Measurements of subtle changes in the Earth’s gravityfield can provideinformation on the addition/loss of mass (e.g., magma or aqueousfluids) beneath a volcano. In this study, wemeasured gravity at Mount St. Helens from 2010 to 2016 to investigate possible mass changes followingthe 2004–2008 dome-forming eruption. The raw gravity measurements were corrected for changes in themass and shape of Crater Glacier and for deformation of the volcanic edifice to obtain residual gravity values.We found positive residual gravity changes that we interpret as evidence for partial recharge of the magmareservoir that fed the 2004–2008 eruption and/or for accumulation of groundwater in one or more shallowaquifers. Most of the gravity signal can be explained by groundwater accumulation; magma recharge is notrequired. However, there is additional evidence for at least partial recharge from seismic, deformation, andgeochemical observations

Washington