Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research--Solid Earth”

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

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

At least 73 records · Page 4Linked to original sources
USGS2023

Quantifying site effects and their influence on earthquake source parameter estimations using a dense array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

Oklahoma↗
USGS2024

Exploring the dynamic interactions between the Southern San Andreas Fault and a normal fault under the Salton Sea

We investigate the dynamic interactions between the Southern San Andreas Fault (SSAF) and a proximal normal fault (NF) beneath the Salton Sea in southern California. The NF, positioned near the SSAF terminus at Bombay Beach, exhibits 11–15 displacement events across 14 stratigraphic sequences, with a range of 0.2–1.4 m of vertical offset since ∼2–3 ka. Notably, four of these events may align temporally with SSAF earthquakes, raising questions about the possible interplay between the two faults. Utilizing dynamic rupture models, we analyze the coseismic interactions between the SSAF and NF, addressing under what conditions the SSAF induces slip on the NF. Our findings reveal that a suite of SSAF ruptures, particularly those propagating from north to south, can trigger slip on the normal fault and replicate observed vertical offsets. If the SSAF extends beneath the Salton Sea, earthquakes originating south of the NF intersection are less likely to trigger normal fault slip, although we cannot exclude this possibility. Some SSAF ruptures do not trigger discernible slip on the NF, rendering such events undetectable in the stratigraphic record. Our research contributes toward discussions regarding the seismic hazard in southern California, shedding light on the interplay between the SSAF and NF.

California↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Earthquake cycle mechanics during caldera collapse: Simulating the 2018 Kīlauea eruption

In multiple observed caldera-forming eruptions, the rock overlying a draining magma reservoir dropped downward along ring faults in sequences of discrete collapse earthquakes. These sequences are analogous to tectonic earthquake cycles and provide opportunities to examine fault mechanics and collapse eruption dynamics over multiple events. Collapse earthquake cycles have been studied with zero-dimensional slider-block models, but these do not account for the complicated interplay between fluid and elastic dynamics or for factors such as the heterogeneous fault properties and non-vertical ring fault geometries often inferred at volcanoes. We present two-dimensional axisymmetric mafic piston-like collapse earthquake cycle models that include rate-and-state friction, fully-dynamic elasticity, and compressible viscous fluid magma flow. We demonstrate that collapse earthquake intervals and magnitudes are highly sensitive to inertial effects, evolving stress fields, fault geometry, and depth-varying fault friction. Given the consistent earthquake cycles observed in most eruptions, this suggests that ring faults can quickly stabilize and often become nearly vertical at depth. We use the well-monitored 2018 collapse sequence at Kı̄lauea as a case study. Our model can produce many features of Kı̄lauea seismic and geodetic observations, except for a significant amount of interseismic slip, which cannot be readily explained with simple rate-and-state friction parameterizations.

Hawaii↗
USGS2024

Subduction zone geometry modulates the megathrust earthquake cycle: Magnitude, recurrence, and variability

Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably-dipping planar and variably-curved nonplanar megathrusts using the volumetric, high-order accurate code tandem to account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly-dipping megathrusts to host larger earthquakes than steeply-dipping ones. Under elevated pore pressure and less strongly velocity-weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly-dipping and sharply-curved megathrusts host multi-period supercycles of slow-to-fast, small-to-large slip events under higher effective stresses and more strongly velocity-weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second-order effects from structurally-derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co- and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Brittle regime slip partitioned damage and deformation mechanisms along the eastern Denali fault zone in southwestern, Yukon

Rare bedrock exposures of the eastern Denali fault zone in southwestern Yukon allow for the measurement, sampling, and analyses of brittle regime fault slip data and deformation mechanisms to explore relations to far field, oblique plate motions. Host rock lithologies and associated slip surfaces show episodic damage zone‐related deformation and calcite ± hematite ± chlorite related hydrothermal fluid flow. This regional scale network of asymmetric fault damage is spatially and kinematically linked to a discrete and narrow fault core. Fault network observations, orientations, slip data, and strain inversions document a slip partitioned strike‐slip fault system with locally and mutually overprinting strike‐, oblique‐, and dip‐slip components. Microstructural analyses reveal crystal plastic and co‐seismic brittle deformation mechanisms active in a narrow range of upper crustal temperature, pressure, fluid, and chemical conditions. The net damage related slip is not exclusively formed by a single kinematic system, but rather a fully partitioned, time integrated system likely operative for much of the fault's brittle regime evolution temporally constrained by previously published thermochronometric data. Although the fault slip data was collected from outcrop‐scale exposures at sites tens of kilometers apart, results show remarkable correlation between fault kinematics and plate motions along the ∼580 km long eastern Denali fault segment. End member, subhorizontal, northeast directed reverse and north directed dextral strike slip fault strain axes closely reflect relative plate motion interactions over at least the last 30 m.y. and act as a proxy for far‐field stresses compatible with the kinematics of the damage zone network.

Alaska↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗
USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗
USGS2025

New paleomagnetic constraints on the eruption timing, stratigraphy, and post-emplacement deformation of the Picture Gorge Basalt within the Columbia River Basalt Group

The Picture Gorge Basalt (PGB) is part of the Miocene Columbia River Basalt Group (CRBG). The PGB, which outcrops in eastern Oregon, is considered coincident in time with the voluminous Grande Ronde Basalt. New radiometric ages have expanded the age‐range of the PGB, suggesting it began erupting prior to the Steens Basalt to its south and continued until after cessation of the Grande Ronde Basalt eruptions, an interval of 1.5 Ma. However, the existing paleomagnetism of the PGB implies this eruption timeline is an overestimate. To reconcile the radiometric and paleomagnetic timescales for the PGB, we conducted a paleomagnetic study on sections of the PGB to construct a detailed, high‐quality magnetostratigraphy. Our data indicate the stratigraphically lowest lava flows in the PGB are of reversed polarity, revealing a new paleomagnetic transition with the PGB and a reversed (R)–normal (N)–reversed (R) sequence. This suggests one of two timeline possibilities for PGB volcanism: (a) eruptions began and during through CRBG polarity chrons R0–N0–R1, penecontemporaneous with Steens Basalt, or (b) eruptions began and persisted during CRBG polarity chrons R1–N1–R2. Our work supports a longer interval of PGB volcanism than was suggested by previous paleomagnetic data but is at odds with the suggestion that PGB eruptions lasted through the entire main CRBG. We favor a scenario wherein PGB eruptions begin with R0 and continue into the R1 paleomagnetic interval. The paleomagnetic results also record a ∼18° vertical‐axis rotation of east‐central Oregon after ∼16 Ma with respect to stable North America.

Oregon↗
USGS2025

Rupture continuity through intermittent pauses in Cascadia slow slip events

Cascadia slow slip events (SSEs) are often envisioned as smooth, continuous ruptures, progressively activating tremor asperities as they propagate. Macroscopically, geodetic inversions and spatiotemporal maps of tremor epicenters show steady, uniform migration. In detail tremor is more chaotic and discontinuous. Larger long-term SSEs observed in daily geodetic solutions are inferred to exhibit intermittent pauses that reflect temporary re-locking of the fault, but this temporal resolution limits tests for similar re-locking on shorter timescales. We use temporal measurements of the areal growth and radiated energy of tremor clusters to investigate SSE intermittence. We find that ruptures mirror tremor pauses. Areal growth rate, however, does not reset, and removing the pauses results in smoother and more similar growth measurements among all SSEs. The rupture similarity occurs regardless of size or location and hints at an underlying uniformity and lack of predeterminism in eventual SSE size. Epicentral uncertainty precludes quantifying early rupture stages, but for larger events areal growth follows a power-law and slows with increasing size. Temporal correlations in tremor energy with inferred SSE propagation velocities and tremor rates suggest its use as a proxy for slip velocity. We find that tremor energy is tidally modulated at daily and sub-daily frequencies, and this modulation is continuous through pauses, suggesting a memory of slip state is sustained through them. We argue these pauses reflect unsteady propagation of the slip front, marked by rapid re- and un-locking, and excluding them removes rupture complexity to reveal a diffusive-like slip process and underlying universality in growth.

California, Oregon, Washington↗
USGS2026

Role of volcano-tectonic interactions during early-phase magma-assisted continental rifting

We study volcanic-tectonic interactions during the early phase of rifting in the Natron Rift, a magmatic segment of the East African Rift System. We use the open-source crustal deformation modeling code, PyLith, to model deflation during the 2007–2008 explosive eruptions and inflation during 2016–2023. We compare predictions of surface displacements resulting from magma reservoir inflation from 3D models with and without topography and homogeneous and heterogeneous material properties against Global Navigation Satellite System displacements from 2016 to 2023. We prefer the model with heterogeneous material properties and topography for fitting the observed displacements. Including topography improves the fit to GNSS data and modulates stress changes from shallow magma sources. Our numerical modeling of Coulomb stress changes suggests that deflation of a magma reservoir under the active Ol Doinyo Lengai volcano could trigger oblique normal slip on the Natron Fault with Coulomb stress increases exceeding ∼0.1 MPa on the volcano-adjacent fault segment. Fault slip on the Natron Fault is less likely during inflation of the magma reservoir because normal compressive stress increases on the fault. We also demonstrate that Coulomb stress changes are sensitive to the inclusion of topography in our model, particularly because the depth of the magma reservoir is similar to the amplitude of the topographic relief.

Ol Doinyo Lenga↗
USGS2012

Compositions, growth mechanisms, and temporal relations of hydrothermal sulfide‐sulfate‐silica chimneys at the northern Cleft segment, Juan de Fuca Ridge

Three active hydrothermal vents forming sulfide mounds and chimneys (Monolith, Fountain, and Pipe Organ) and more widely distributed inactive chimneys are spatially related to a system of discontinuous fissures and young sheet flow lavas at the northern Cleft segment, Juan de Fuca Ridge. The formation of zoned tubular Curich chimneys (type I) on the Monolith sulfide mound is related to focused flow of high‐temperature (to 328°C) fluid. Bulbous chimneys (type II or “beehives”) at the Monolith and Fountain vents are products of diffuse high‐temperature (to 315°C) discharge. A broader zone of vigorous mixing between the hydrothermal fluid and seawater results in quench crystallization of anhydrite‐rich shells. Columnar Zn‐sulfide‐rich chimneys with narrow channelways (type III) are constructed where focused and relatively low‐temperature (261°C) fluid vents directly from the basalt substrate. The bulk chemistry (low Cu; high Pb, Ag, and SiO 2 contents), mineralogy (pyrite‐marcasite‐wurtzite‐amorphous silica‐anglesite), colloform and filamentous textures, and oxygen isotope characteristics of inactive (type IV) chimneys indicate a low‐temperature (<250°C) origin involving diffuse and sluggish flow patterns and conductive cooling. Seafloor observations and 210 Pb data indicate that (1) type IV chimneys are products of an earlier period of hydrothermal activity that ended no more than 60 years ago but prior to the sheet flow eruption, (2) the high‐temperature Monolith and Fountain vents are manifestations of the same heating event (shallow emplacement of magma) that led to the sheet flow eruption and recent megaplumes, and (3) the Pipe Organ Vent is in a very youthful stage of development and chimney deposition postdates the sheet flow eruption.

California, Oregon, Washington↗
USGS2012

Surface seismic and electrical methods to detect fluids related to faulting

In the absence of drilling, surface-based geophysical methods are necessary to observe fault zones and fault zone physical properties at seismogenic depths. These in situ physical properties can then be used to infer the presence and distribution of fluids along faults, although such observations are by nature indirect and become less exact with greater depth. Multiple observations of a range of such geophysical properties as compressional and shear seismic velocity ( V p and V s ), V p /V5 ratio (related to Poisson's ratio), resistivity and attenuation in and adjacent to fault zones offer the greatest hope of making inferences of the fault zone geometry, fluids in the fault zone, and fluid reservoirs in the surrounding crust. For simple geometries, fault zone guided waves can provide information on fault zone width and velocities for faults of the order of 200 m wide. To address the question of whether a narrow fault zone can be imaged well enough at depths of seismic rupture to infer the presence of anomalously high fluid/rock ratios, we present synthetic seismic tomography and magnetotelluric examples for an ideal case of a narrow fault zone with a simple geometry, large changes in material properties, and numerous earthquakes within the fault zone. A synthetic 0.5-km wide fault zone with 20% velocity reduction is well imaged using local earthquake tomography. When sequential velocity inversions are done, the true fault width is found, even to 9 km depth, although the calculated amplitude of the velocity reduction is lower than the actual amplitude. V p / V s is as well determined as V p . Magnetotelluric imaging of a synthetic fault zone shows that a conductive fault zone can be well imaged within the upper 10 km. Further, a narrow (1 km) very low resistivity (3 ohm m) fault core can be imaged within a broad (5 km) low resistivity (10 ohm m) fault zone, illustrating that regions of a fault containing large quantities of interconnected fluids within a broader, conductive fault zone should be detectable. Thus variations in fluid content and fluid pressure can be inferred from electrical and seismic methods but there will always be uncertainty in these inferences due to the trade-off with other factors, such as intrinsic variations in porosity, mineralogy, and pore geometry. The best approach is combined modeling of varied seismic and electrical data.

Journal of Geophysical Research Solid Earth↗
USGS1998

The coseismic slip distributions of the 1940 and 1979 Imperial Valley, California, earthquakes and their implications

Geodetic arrays observed by the U.S. Coast and Geodetic Survey span the Imperial fault in southern California. For the 1940 M 7.1 Imperial Valley earthquake, a 1934–1941 triangulation network has sufficient resolution to allow inversion for the coseismic slip distribution on fault segments 5 to 25 km long extending from the surface to a depth of 9 km. The estimated right-lateral slip is 0.8 to 1.7 m on the northern 30 km of the main trace of the Imperial fault, 4.8±0.2 m on a 10-km-long segment straddling the United States - Mexico border, and 1.3±0.4 m on a southern 25-km-long segment in Mexico. Fixing this strike-slip model and inverting 1940 leveling data only for dip slip yields 0.1 m of east-side-down dip slip. The seismic moment for this model is M 0 = (3.2±0.3) ×10 19 N m. The 1979 geodetic data set, mostly elevation changes from leveling routes, has insufficient resolution for inversion. However, it is possible to use this geodetic data set and results published by others to infer that the 1940 and 1979 earthquakes may be similar on the rupture zone common to both events. Our preferred 1940 model is similar to the 1979 geodetic results of Crook [1984] on the segments where both networks have good resolution. Elevation changes from 1940 and 1979 leveling data are very similar. Thus the geodetic data corroborate the surface slip evidence of Sharp [1982b] that the 1940 and 1979 slip distributions are examples of “characteristic slip” on the northern Imperial fault.

California↗
USGS1998

Suppression of large earthquakes by stress shadows: A comparison of Coulomb and rate-and-state failure

Stress shadows generated by California's two most recent great earthquakes (1857 Fort Tejon and 1906 San Francisco) substantially modified 19th and 20th century earthquake history in the Los Angeles basin and in the San Francisco Bay area. Simple Coulomb failure calculations, which assume that earthquakes can be modeled as static dislocations in an elastic half-space, have done quite well at approximating how long the stress shadows, or relaxing effects, should last and at predicting where subsequent large earthquakes will not occur. There has, however, been at least one apparent exception to the predictions of such simple models. The 1911 M >6.0 earthquake near Morgan Hill, California, occurred at a relaxed site on the Calaveras fault. We examine how the more complex rate-and-state friction formalism based on laboratory experiments might have allowed the 1911 earthquake. Rate-and-state time-to-failure calculations are consistent with the occurrence of the 1911 event just 5 years after 1906 if the Calaveras fault was already close to failure before the effects of 1906. We also examine the likelihood that the entire 78 years of relative quiet (only four M ≥6 earthquakes) in the bay area after 1906 is consistent with rate-and-state assumptions, given that the previous 7 decades produced 18 M ≥6 earthquakes. Combinations of rate-and-state variables can be found that are consistent with this pattern of large bay area earthquakes, assuming that the rate of earthquakes in the 7 decades before 1906 would have continued had 1906 not occurred. These results demonstrate that rate-and-state offers a consistent explanation for the 78-year quiescence and the 1911 anomaly, although they do not rule out several alternate explanations.

California↗
USGS1998

Strain accumulation rates in the San Francisco Bay area, 1972–1989

Maps of the strain accumulation rate in the San Francisco Bay area have been constructed from trilateration observations extending from about 1972 until the Loma Prieta earthquake in late 1989. The observations were corrected to remove offsets imposed by shallow fault creep and by four M ∼6 earthquakes that occurred in the Bay area during that time interval. The Bay area was divided into 32 contiguous polygons, and the uniform (in both space and time) strain rates that best explain the changes in the corrected (earthquake and shallow fault creep offsets removed) distances within each polygon were calculated. In a coordinate system with the 1 axis directed N58°E and the 2 axis N32°W (perpendicular and parallel to the local tangent to the small circle drawn about the Pacific-Sierra Nevada pole of rotation) the averages of these 32 strain rates (each weighted by the area of the polygon) are ε11=9.2±7.4, ε12=−160.7±4.6, and ε22=8.2±6.2 nanostrain/yr, where extension is reckoned positive and quoted uncertainties are standard deviations. As expected from the Pacific-Sierra Nevada relative plate motion, the overall strain rate is predominantly right-lateral shear across a vertical plane striking N32°W. The net increase in the 12,225 km2 area of the trilateration network is only 212±110 m2/yr, which arises from almost equal extensions in the N32°W and N58°E directions. Within the network the strain rates vary from polygon to polygon. Significant areal dilatation rates are observed in almost ⅓ of the individual polygons and the N32°W extension rates tend to be negative to the west of the Hayward-Rodgers Creek fault trend and positive east of it. The N58E extension rate is positive in 22 out of 32 polygons, a proportion that is significantly larger than would be expected by chance if the N58°E extension rate were zero or negative. The pre-1989 strain accumulation across the eventual site of the Loma Prieta rupture involves fault normal contraction as well as right-lateral shear, consistent with the rupture mechanism.

California↗
USGS2020

Timescales and processes of methane hydrate formation and breakdown, with application to geologic systems

Gas hydrate is an ice-like form of water and low molecular weight gas stable at temperatures of roughly -10ºC to 25ºC and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one-sixth of Earth’s methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean-atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean-atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth’s near-surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some fields studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

Journal of Geophysical Research-Solid Earth↗
USGS2021

Temporal seismic velocity variations: Recovery following from the 2019 Mw 7.1 Ridgecrest, California earthquake

We investigated seismic velocity changes ( dv / v ) associated with the 2019 Ridgecrest earthquake sequence with high‐frequency autocorrelations of ambient seismic noise data. Daily autocorrelation functions were computed for the entirety of 2019 and the first quarter of 2020 for broadband stations within the region, including the temporary broadband stations installed during the aftershock deployment. Travel time shifts in the daily autocorrelation functions, relative to the mean autocorrelation waveform, were computed to produce dv / v time series, which are sensitive to the evolving material properties of the shallow crust surrounding the Ridgecrest fault zone (RFZ). A short‐term velocity drop follows the M w 7.1 earthquake at stations in the vicinity of the rupture surface, while those greater than 50 km away showed no such drop. The maximum, absolute changes in seismic velocity are proportional to the logarithm of distance from the fault rupture and to the peak dynamic strain experienced during the earthquake. Near the areas of the highest coseismic slip within the RFZ, seismic velocities recovered over 3 months. However, in the vicinity of the nearby Garlock fault, where triggered slip manifested, and north of the RFZ, seismic velocities recovered within a month. We interpret the seismic velocity changes and their recovery to be largely due to changes in the physical properties of the shallow crust, such as fault zone damage recovery caused by the earthquake rupture process and in response to the large dynamic stresses of passing seismic waves from the mainshock.

California↗