SEARCH · USGS Science

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5Linked to original sources

USGS2014

Widespread gas hydrate instability on the upper U.S. Beaufort margin

The most climate-sensitive methane hydrate deposits occur on upper continental slopes at depths close to the minimum pressure and maximum temperature for gas hydrate stability. At these water depths, small perturbations in intermediate ocean water temperatures can lead to gas hydrate dissociation. The Arctic Ocean has experienced more dramatic warming than lower latitudes, but observational data have not been used to study the interplay between upper slope gas hydrates and warming ocean waters. Here we use (a) legacy seismic data that constrain upper slope gas hydrate distributions on the U.S. Beaufort Sea margin, (b) Alaskan North Slope borehole data and offshore thermal gradients determined from gas hydrate stability zone thickness to infer regional heat flow, and (c) 1088 direct measurements to characterize multidecadal intermediate ocean warming in the U.S. Beaufort Sea. Combining these data with a three-dimensional thermal model shows that the observed gas hydrate stability zone is too deep by 100 to 250 m. The disparity can be partially attributed to several processes, but the most important is the reequilibration (thinning) of gas hydrates in response to significant (~0.5°C at 2σ certainty) warming of intermediate ocean temperatures over 39 years in a depth range that brackets the upper slope extent of the gas hydrate stability zone. Even in the absence of additional ocean warming, 0.44 to 2.2 Gt of methane could be released from reequilibrating gas hydrates into the sediments underlying an area of ~5–7.5 × 10 3 km 2 on the U.S. Beaufort Sea upper slope during the next century.

Alaska↗

USGS2014

Attenuation and scattering tomography of the deep plumbing system of Mount St. Helens

We present a combined 3-D P wave attenuation, 2-D S coda attenuation, and 3-D S coda scattering tomography model of fluid pathways, feeding systems, and sediments below Mount St. Helens (MSH) volcano between depths of 0 and 18 km. High-scattering and high-attenuation shallow anomalies are indicative of magma and fluid-rich zones within and below the volcanic edifice down to 6 km depth, where a high-scattering body outlines the top of deeper aseismic velocity anomalies. Both the volcanic edifice and these structures induce a combination of strong scattering and attenuation on any seismic wavefield, particularly those recorded on the northern and eastern flanks of the volcanic cone. North of the cone between depths of 0 and 10 km, a low-velocity, high-scattering, and high-attenuation north-south trending trough is attributed to thick piles of Tertiary marine sediments within the St. Helens Seismic Zone. A laterally extended 3-D scattering contrast at depths of 10 to 14 km is related to the boundary between upper and lower crust and caused in our interpretation by the large-scale interaction of the Siletz terrane with the Cascade arc crust. This contrast presents a low-scattering, 4–6 km 2 “hole” under the northeastern flank of the volcano. We infer that this section represents the main path of magma ascent from depths greater than 6 km at MSH, with a small north-east shift in the lower plumbing system of the volcano. We conclude that combinations of different nonstandard tomographic methods, leading toward full-waveform tomography, represent the future of seismic volcano imaging.

Washington↗

USGS2014

Eruptions at Lone Star geyser, Yellowstone National Park, USA: 2. Constraints on subsurface dynamics

We use seismic, tilt, lidar, thermal, and gravity data from 32 consecutive eruption cycles of Lone Star geyser in Yellowstone National Park to identify key subsurface processes throughout the geyser's eruption cycle. Previously, we described measurements and analyses associated with the geyser's erupting jet dynamics. Here we show that seismicity is dominated by hydrothermal tremor (~5–40 Hz) attributed to the nucleation and/or collapse of vapor bubbles. Water discharge during eruption preplay triggers high-amplitude tremor pulses from a back azimuth aligned with the geyser cone, but during the rest of the eruption cycle it is shifted to the east-northeast. Moreover, ~4 min period ground surface displacements recur every 26 ± 8 min and are uncorrelated with the eruption cycle. Based on these observations, we conclude that (1) the dynamical behavior of the geyser is controlled by the thermo-mechanical coupling between the geyser conduit and a laterally offset reservoir periodically filled with a highly compressible two-phase mixture, (2) liquid and steam slugs periodically ascend into the shallow crust near the geyser system inducing detectable deformation, (3) eruptions occur when the pressure decrease associated with overflow from geyser conduit during preplay triggers an unstable feedback between vapor generation (cavitation) and mass discharge, and (4) flow choking at a constriction in the conduit arrests the runaway process and increases the saturated vapor pressure in the reservoir by a factor of ~10 during eruptions.

Wyoming↗

USGS2014

The impact of static stress change, dynamic stress change, and the background stress on aftershock focal mechanisms

The focal mechanisms of earthquakes in Southern California before and after four M  ≥ 6.7 main shocks provide insight into how fault systems respond to stress and changes in stress. The main shock static stress changes have two observed impacts on the seismicity: changing the focal mechanisms in a given location to favor those aligned with the static stress change and changing the spatial distribution of seismicity to favor locations where the static stress change aligns with the background stress. The aftershock focal mechanisms are significantly aligned with the static stress changes for absolute stress changes of ≥ 0.02 MPa, for up to ~20 years following the main shock. The dynamic stress changes have similar, although smaller, effects on the local focal mechanisms and the spatial seismicity distribution. Dynamic stress effects are best observed at long periods (30–60 s) and for metrics based on repeated stress cycling in the same direction. This implies that dynamic triggering operates, at least in part, through cyclic shear stress loading in the direction of fault slip. The background stress also strongly controls both the preshock and aftershock mechanisms. While most aftershock mechanisms are well oriented in the background stress field, 10% of aftershocks are identified as poorly oriented outliers, which may indicate limited heterogeneity in the postmain shock stress field. The fault plane orientations of the outliers are well oriented in the background stress, while their slip directions are not, implying that the background stress restricts the distribution of available fault planes.

California↗

USGS2015

Postearthquake relaxation evidence for laterally variable viscoelastic structure and water content in the Southern California mantle

I reexamine the lower crust and mantle relaxation following two large events in the Mojave Desert: the 1992 M 7.3 Landers and 1999 M 7.1 Hector Mine, California, earthquakes. Time series from continuous GPS sites out to 300 km from the ruptures are used to constrain models of postseismic relaxation. Crustal motions in the Mojave Desert region are elevated above background for several years following each event. To account for broadscale relaxation of the lower crust and mantle, the Burgers body model is employed, involving transient and steady state viscosities. Joint afterslip/postseismic relaxation modeling of the GPS time series up to one decade following the Hector Mine earthquake reveals a significant rheological contrast between a northwest trending “southwest domain” (that envelopes the San Andreas fault system and western Mojave Desert) and an adjacent “northeast domain” (that envelopes the Landers and Hector Mine rupture areas in the central Mojave Desert). The steady state viscosity of the northeast domain mantle asthenosphere is inferred to be ∼4 times greater than that of the southwest domain. This pattern is counter to that expected for regional heat flow, which is higher in the northeast domain, but it is explicable by means of a nonlinear rheology that includes dependence on both strain rate and water concentration. I infer that the southwest domain mantle has a relatively low steady state viscosity because of its high strain rate and water content. The relatively low mantle water content of the northeast domain is interpreted to result from the continual extraction of water through igneous and volcanic activity over the past ∼20 Myr.

California↗

USGS2014

Ice sheet load cycling and fluid underpressures in the Eastern Michigan Basin, Ontario, Canada

Strong fluid underpressures have been detected in Paleozoic strata in the eastern Michigan Basin, with hydraulic heads reaching ~400 m below land surface (~4 MPa underpressure) and ~200 m below sea level in strata where unusually low permeabilities (~10 −20 –10 −23  m 2 ) were measured in situ. Multiple glaciations, including three with as much as 3 km of ice cover at the site in the last 120 ka, suggest a causal link with the underpressures. We examined this possibility using a one-dimensional groundwater flow model incorporating mechanical loading from both ice weight and lithospheric flexure. Because hydrologic and mechanical changes during glaciation are not well characterized and subsurface properties are imperfectly known, the model was used inversely to estimate flexural loads and loosely constrained permeabilities by matching observed pressures. Acceptable matches were obtained for a surprisingly wide range of scenarios with permeabilities close to measured values and plausible flexural loads. Matches were not obtained when too many parameters were preselected, or when permeabilities were constrained to be significantly larger than measured values. In successful model runs groundwater expulsion under glacial-mechanical loads caused the underpressuring, and flexural loads were important if aquifer and sub-glacial pressures were significantly elevated during glaciation. Simulated fluid pressures in the low-permeability strata fluctuated by 30–40 MPa during glacial cycles but resulted in advective transport of only tens of meters or less. Although other mechanisms cannot be ruled out, we conclude that glacial-mechanical forcing of a water-saturated system can explain the observed underpressures.

Ontario↗

USGS2015

Near-surface versus fault zone damage following the 1999 Chi-Chi earthquake: Observation and simulation of repeating earthquakes

We observe crustal damage and its subsequent recovery caused by the 1999  M 7.6 Chi-Chi earthquake in central Taiwan. Analysis of repeating earthquakes in Hualien region, ~70 km east of the Chi-Chi earthquake, shows a remarkable change in wave propagation beginning in the year 2000, revealing damage within the fault zone and distributed across the near surface. We use moving window cross correlation to identify a dramatic decrease in the waveform similarity and delays in the S wave coda. The maximum delay is up to 59 ms, corresponding to a 7.6% velocity decrease averaged over the wave propagation path. The waveform changes on either side of the fault are distinct. They occur in different parts of the waveforms, affect different frequencies, and the size of the velocity reductions is different. Using a finite difference method, we simulate the effect of postseismic changes in the wavefield by introducing S wave velocity anomaly in the fault zone and near the surface. The models that best fit the observations point to pervasive damage in the near surface and deep, along-fault damage at the time of the Chi-Chi earthquake. The footwall stations show the combined effect of near-surface and the fault zone damage, where the velocity reduction (2–7%) is twofold to threefold greater than the fault zone damage observed in the hanging wall stations. The physical models obtained here allow us to monitor the temporal evolution and recovering process of the Chi-Chi fault zone damage.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Self-similar rupture implied by scaling properties of volcanic earthquakes occurring during the 2004-2008 eruption of Mount St. Helens, Washington

We analyze a group of 6073 low-frequency earthquakes recorded during a week-long temporary deployment of broadband seismometers at distances of less than 3 km from the crater at Mount St. Helens in September of 2006. We estimate the seismic moment ( M 0 ) and spectral corner frequency ( f 0 ) using a spectral ratio approach for events with a high signal-to-noise (SNR) ratio that have a cross-correlation coefficient of 0.8 or greater with at least five other events. A cluster analysis of cross-correlation values indicates that the group of 421 events meeting the SNR and cross-correlation criteria forms eight event families that exhibit largely self-similar scaling. We estimate the M 0 and f 0 values of the 421 events and calculate their static stress drop and scaled energy ( E R / M 0 ) values. The estimated values suggest self-similar scaling within families, as well as between five of eight families (i.e., and constant). We speculate that differences in scaled energy values for the two families with variable scaling may result from a lack of resolution in the velocity model. The observation of self-similar scaling is the first of its kind for such a large group of low-frequency volcanic tectonic events occurring during a single active dome extrusion eruption.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Seismic source dynamics of gas-piston activity at Kı̄lauea Volcano, Hawai‘i

Since 2008, eruptive activity at the summit of Kı̄lauea Volcano, Hawai‘i has been confined to the new Overlook pit crater within the Halema‘uma‘u Crater. Among the broad range of magmatic processes observed in the new pit are recurring episodes of gas pistoning. The gas-piston activity is accompanied by seismic signals that are recorded by a broadband network deployed in the summit caldera. We use raw data recorded with this network to model the source mechanism of representative gas-piston events in a sequence that occurred on 20–25 August 2011 during a gentle inflation of the Kı̄lauea summit. 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. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–10,000 s. Most of the seismic wavefield produced by gas-pistoning originates in a source region ∼1 km below the eastern perimeter of the Halema‘uma‘u pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks featuring an east striking crack (dike) dipping 80°to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each gas-piston event is marked by a similar rapid inflation lasting a few minutes, trailed by a slower deflation ramp extending up to 15 min, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes accompanying the growth and collapse of a layer of foam at the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes : foam thickness (10–50 m), foam cell diameter (0.04–0.10 m), and gas-injection velocity (0.01–0.06 m s −1 ). Gas-piston activity occurs in a narrow pipe with diameter of 6 m connecting the Halema‘uma‘u pit crater to the subjacent dike system. The height of the magma column is estimated at ∼104 m at the start of the sequence based on the period of very long period (VLP) oscillations accompanying the onset of the gas-piston signal. Based on the change in the period of VLP oscillations and tilt evidence, the height of the magma column is inferred to have risen by up to ∼23 m by the end of the 5 day long sequence. A penny-shaped crack model of the dike geometry yields effective diameters of ∼1.2–2.9 km for the east dike and 0.7 km for the north dike. The shallower north dike segment is embedded in a relatively weak medium, compatible with expected mechanical properties in the hydrothermal environment of this dike.

Hawaii↗

USGS2015

Direct measurement of asperity contact growth in quartz at hydrothermal conditions

Earthquake recurrence requires interseismic fault restrengthening which results from solid state deformation in room-temperature friction and indentation experiments. In contrast exhumed fault zones show solution-transport processes such as pressure solution and contact overgrowths influence fault zone properties . In the absence of fluid flow, overgrowths are driven by gradients in surface curvature where material is dissolved, diffuses, and precipitates at the contact without convergence normal to the contact. To determine the rate of overgrowth for quartz, we conducted single contact experiments in an externally heated pressure vessel. Convergence was continuously monitored using reflected-light interferometry through a long-working-distance microscope. Contact normal force was constant with an initial effective normal stress of 1.7 MPa, temperature was between 350 and 530{degree sign}C, and water pressure was constant at 150 MPa. Two control experiments were conducted: one dry at 425{degree sign}C and one bi-material (sapphire) at 425{degree sign}C and 150 MPa water pressure. No contact growth or convergence was observed in the controls. For wet single-phase contacts, growth was initially rapid and then decreased with time. No convergence was observed. Fluid inclusions indicate that the contact is not uniformly wetted. The contact is bounded by small regions of high aperture, reflecting local free-face dissolution as the source for the overgrowth. The apparent activation energy is ~125 kJ/mol. Extrapolation predicts rates of contact area increase orders of magnitude faster than in dry, room-temperature and hydrothermal friction experiments, suggesting that natural strength recovery near the base of the seismogenic zone could be dominated by contact overgrowth.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

The Palos Verdes Fault offshore southern California: late Pleistocene to present tectonic geomorphology, seascape evolution and slip rate estimate based on AUV and ROV surveys

The Palos Verdes Fault (PVF) is one of few active faults in Southern California that crosses the shoreline and can be studied using both terrestrial and subaqueous methodologies. To characterize the near-seafloor fault morphology, tectonic influences on continental slope sedimentary processes and late Pleistocene to present slip rate, a grid of high-resolution multibeam bathymetric data, and chirp subbottom profiles were acquired with an autonomous underwater vehicle (AUV) along the main trace of PVF in water depths between 250 and 600 m. Radiocarbon dates were obtained from vibracores collected using a remotely operated vehicle (ROV) and ship-based gravity cores. The PVF is expressed as a well-defined seafloor lineation marked by subtle along-strike bends. Right-stepping transtensional bends exert first-order control on sediment flow dynamics and the spatial distribution of Holocene depocenters; deformed strata within a small pull-apart basin record punctuated growth faulting associated with at least three Holocene surface ruptures. An upper (shallower) landslide scarp, a buried sedimentary mound, and a deeper scarp have been right-laterally offset across the PVF by 55 ± 5, 52 ± 4 , and 39 ± 8 m, respectively. The ages of the upper scarp and buried mound are approximately 31 ka; the age of the deeper scarp is bracketed to 17–24 ka. These three piercing points bracket the late Pleistocene to present slip rate to 1.3–2.8 mm/yr and provide a best estimate of 1.6–1.9 mm/yr. The deformation observed along the PVF is characteristic of strike-slip faulting and accounts for 20–30% of the total right-lateral slip budget accommodated offshore Southern California.

California↗

USGS2015

Holocene geologic slip rate for the Banning strand of the southern San Andreas Fault, southern California

Northwest directed slip from the southern San Andreas Fault is transferred to the Mission Creek, Banning, and Garnet Hill fault strands in the northwestern Coachella Valley. How slip is partitioned between these three faults is critical to southern California seismic hazard estimates but is poorly understood. In this paper, we report the first slip rate measured for the Banning fault strand. We constrain the depositional age of an alluvial fan offset 25 ± 5 m from its source by the Banning strand to between 5.1 ± 0.4 ka (95% confidence interval (CI)) and 6.4 + 3.7/−2.1 ka (95% CI) using U-series dating of pedogenic carbonate clast coatings and 10 Be cosmogenic nuclide exposure dating of surface clasts. We calculate a Holocene geologic slip rate for the Banning strand of 3.9 + 2.3/−1.6 mm/yr (median, 95% CI) to 4.9 + 1.0/−0.9 mm/yr (median, 95% CI). This rate represents only 25–35% of the total slip accommodated by this section of the southern San Andreas Fault, suggesting a model in which slip is less concentrated on the Banning strand than previously thought. In rejecting the possibility that the Banning strand is the dominant structure, our results highlight an even greater need for slip rate and paleoseismic measurements along faults in the northwestern Coachella Valley in order to test the validity of current earthquake hazard models. In addition, our comparison of ages measured with U-series and 10 Be exposure dating demonstrates the importance of using multiple geochronometers when estimating the depositional age of alluvial landforms.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Crustal deformation in the New Madrid seismic zone and the role of postseismic processes

Global Navigation Satellite System data across the New Madrid seismic zone (NMSZ) in the central United States over the period from 2000 through 2014 are analyzed and modeled with several deformation mechanisms including the following: (1) creep on subsurface dislocations, (2) postseismic frictional afterslip and viscoelastic relaxation from the 1811–1812 and 1450 earthquakes in the NMSZ, and (3) regional strain. In agreement with previous studies, a dislocation creeping at about 4 mm/yr between 12 and 20 km depth along the downdip extension of the Reelfoot fault reproduces the observations well. We find that a dynamic model of postseismic frictional afterslip from the 1450 and February 1812 Reelfoot fault events can explain this creep. Kinematic and dynamic models involving the Cottonwood Grove fault provide minimal predictive power. This is likely due to the smaller size of the December 1811 event on the Cottonwood Grove fault and a distribution of stations better suited to constrain localized strain across the Reelfoot fault. Regional compressive strain across the NMSZ is found to be less than 3 × 10 −9 /yr. If much of the present-day surface deformation results from afterslip, it is likely that many of the earthquakes we see today in the NMSZ are aftershocks from the 1811–1812 New Madrid earthquakes. Despite this conclusion, our results are consistent with observations and models of intraplate earthquake clustering. Given this and the recent paleoseismic history of the region, we suggest that seismic hazard is likely to remain significant.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Changes in seismic velocity during the first 14 months of the 2004–2008 eruption of Mount St. Helens, Washington

Mount St. Helens began erupting in late 2004 following an 18 year quiescence. Swarms of repeating earthquakes accompanied the extrusion of a mostly solid dacite dome over the next 4 years. In some cases the waveforms from these earthquakes evolved slowly, likely reflecting changes in the properties of the volcano that affect seismic wave propagation. We use coda-wave interferometry to quantify small changes in seismic velocity structure (usually <1%) between two similar earthquakes and employed waveforms from several hundred families of repeating earthquakes together to create a continuous function of velocity change observed at permanent stations operated within 20 km of the volcano. The high rate of earthquakes allowed tracking of velocity changes on an hourly time scale. Changes in velocity were largest near the newly extruding dome and likely related to shallow deformation as magma first worked its way to the surface. We found strong correlation between velocity changes and the inverse of real-time seismic amplitude measurements during the first 3 weeks of activity, suggesting that fluctuations of pressure in the shallow subsurface may have driven both seismicity and velocity changes. Velocity changes during the remainder of the eruption likely result from a complex interplay of multiple effects and are not well explained by any single factor alone, highlighting the need for complementary geophysical data when interpreting velocity changes.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Effective stress, friction and deep crustal faulting

Studies of crustal faulting and rock friction invariably assume the effective normal stress that determines fault shear resistance during frictional sliding is the applied normal stress minus the pore pressure. Here we propose an expression for the effective stress coefficient α f at temperatures and stresses near the brittle-ductile transition (BDT) that depends on the percentage of solid-solid contact area across the fault. α f varies with depth and is only near 1 when the yield strength of asperity contacts greatly exceeds the applied normal stress. For a vertical strike-slip quartz fault zone at hydrostatic pore pressure and assuming 1 mm and 1 km shear zone widths for friction and ductile shear, respectively, the BDT is at ~13 km. α f near 1 is restricted to depths where the shear zone is narrow. Below the BDT α f = 0 is due to a dramatically decreased strain rate. Under these circumstances friction cannot be reactivated below the BDT by increasing the pore pressure alone and requires localization. If pore pressure increases and the fault localizes back to 1 mm, then brittle behavior can occur to a depth of around 35 km. The interdependencies among effective stress, contact-scale strain rate, and pore pressure allow estimates of the conditions necessary for deep low-frequency seismicity seen on the San Andreas near Parkfield and in some subduction zones. Among the implications are that shear in the region separating shallow earthquakes and deep low-frequency seismicity is distributed and that the deeper zone involves both elevated pore fluid pressure and localization.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Resolving bathymetry from airborne gravity along Greenland fjords

Recent glacier mass loss in Greenland has been attributed to encroaching warming waters, but knowledge of fjord bathymetry is required to investigate this mechanism. The bathymetry in many Greenland fjords is unmapped and difficult to measure. From 2010 to 2012, National Aeronautics and Space Administration's Operation IceBridge collected a unique set of airborne gravity, magnetic, radar, and lidar data along the major outlet glaciers and fjords in Greenland. We applied a consistent technique using the IceBridge gravity data to create 90 bathymetric profiles along 54 Greenland fjords. We also used this technique to recover subice topography where warm or crevassed ice prevents the radar system from imaging the bed. Here we discuss our methodology, basic assumptions and error analysis. We present the new bathymetry data and discuss observations in six major regions of Greenland covered by IceBridge. The gravity models provide a total of 1950 line kilometers of bathymetry, 875 line kilometers of subice topography, and 12 new grounding line depths.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Tomographic Rayleigh-wave group velocities in the Central Valley, California centered on the Sacramento/San Joaquin Delta

If shaking from a local or regional earthquake in the San Francisco Bay region were to rupture levees in the Sacramento/San Joaquin Delta then brackish water from San Francisco Bay would contaminate the water in the Delta: the source of fresh water for about half of California. As a prelude to a full shear-wave velocity model that can be used in computer simulations and further seismic hazard analysis, we report on the use of ambient noise tomography to build a fundamental-mode, Rayleigh-wave group velocity model for the region around the Sacramento/San Joaquin Delta in the western Central Valley, California. Recordings from the vertical component of about 31 stations were processed to compute the spatial distribution of Rayleigh wave group velocities. Complex coherency between pairs of stations were stacked over 8 months to more than a year. Dispersion curves were determined from 4 to about 18 seconds. We calculated average group velocities for each period and inverted for deviations from the average for a matrix of cells that covered the study area. Smoothing using the first difference is applied. Cells of the model were about 5.6 km in either dimension. Checkerboard tests of resolution, which is dependent on station density, suggest that the resolving ability of the array is reasonably good within the middle of the array with resolution between 0.2 and 0.4 degrees. Overall, low velocities in the middle of each image reflect the deeper sedimentary syncline in the Central Valley. In detail, the model shows several centers of low velocity that may be associated with gross geologic features such as faulting along the western margin of the Central Valley, oil and gas reservoirs, and large cross cutting features like the Stockton arch. At shorter periods around 5.5s, the model’s western boundary between low and high velocities closely follows regional fault geometry and the edge of a residual isostatic gravity low. In the eastern part of the valley, the boundaries of the low velocity zone and gravity anomaly are better aligned at longer periods (around 10.5s) suggesting that the eastern edge of the gravity low is associated with deeper structure. There is a strong correspondence between a low in gravity near the Kirby Hills fault and low velocities from the ambient noise tomography. At longer periods, higher velocities creep in from the east and narrow the overall dimension defined by the lower velocities. Overall, there is a strong correspondence between the shape and location of low velocities in the Rayleigh wave velocity images, and geological and geophysical features.

California↗