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 55 records · Page 3Linked to original sources
USGS2021

Geodetic constraints on a 25-year magmatic inflation episode near Three Sisters, central Oregon

Crustal inflation near the Three Sisters volcanic center documented since the mid-1990s has persisted for more than two decades. We update past analyses of the event through 2020 by simultaneously inverting InSAR interferograms, GPS time series, and leveling data for time-dependent volcanic deformation source parameters. We explore several source models to estimate how the deformation rate varied through time and to identify parameters that can reproduce measured deformation. Our preferred model is a Mogi source 4.1 km below sea level (5.9 km below the surface) about 5 km west of the summit of South Sister. Inflation started in late 1995 or 1996; the rate increased rapidly during 1998–1999, and peaked in late 1999, resulting in maximum surface uplift of about 30 cm by mid-2020. Since 2000, the inflation rate generally declined exponentially with a time constant of about 6 years. Two source inflation scenarios fit the data equally well. In the first, the crust surrounding the source is elastic and the net source-volume increase, which we attribute to persistent magma input, has been about 49 × 10 6 m 3 . The second scenario adds a viscoelastic shell surrounding the Mogi source. In that case, an injection of about 21 × 10 6 m 3 of magma prior to 2000, followed by continuing relaxation of the viscoelastic shell, can account for most of the observed surface deformation. In both scenarios, modeling reveals quasiperiodic increases in the inflation rate (pulses) with a recurrence interval of 3–4 years, both before and after 2000.

Oregon↗
USGS2021

Cataloging tectonic tremor energy radiation in the Cascadia subduction zone

For the past ∼12 years the Pacific Northwest Seismic Network has been automatically detecting and locating tectonic tremor across the Cascadia subduction zone, resulting in a catalog of more than 500,000 tremor epicenters to date, which has served as a valuable resource for tremor and slip research. This manuscript presents an updated methodology for routine tremor detection in Cascadia and a new catalog of over 180,000 tremor epicenters including amplitudes detected along the subduction zone margin from 2017 to 2021. The events are detected via cross-correlation of continuous vertical envelope data of 128 stations from northern California to northern Vancouver Island. The modified approach results in less scatter and a 55% increase in detected epicenters than previously observed, as well as a newly identified tremor source offset updip from the main tremor and slip region at the southern edge of the subduction zone. Radiated seismic energy in the 1.5–5 Hz band is used to assign epicenters an energy magnitude ( M eL ), which is calibrated to the M L of local earthquakes. Southern Cascadia is most active, but the highest tremor energy rates occur in northern Cascadia. Tremor in central Cascadia is systematically weaker and less frequent. Individual epicenter magnitudes range from ∼0.5–2 and spatiotemporally cluster into 1,060 swarms with cumulative M eL ranging from ∼0.8 to 3.7. The swarms reflect underlying slow slip events and occur with an earthquake-like energy distribution with a b value ∼1. Tremor epicenters, however, follow a tapered Gutenberg-Richter distribution with high b values, suggesting individual tremor bursts and their constituent low-frequency earthquakes are fault-dimension limited.

California, Oregon, Washington↗
USGS2022

Possible effects of multiphase methane evolution during a glacial cycle on underpressure development in sedimentary basins: An analysis with application to the northeast Michigan Basin

In low-permeability systems, groundwater may be accompanied by separate-phase fluids, and measured pore water pressures may deviate from those expected in steady-state, single-phase systems. These same systems may be of interest for storage of nuclear waste in Deep Geologic Repositories. Therefore, it is important to understand the relationship between the presence of a separate phase and anomalous pressure development. At the Bruce site in Southern Ontario, a significant underpressure was observed, and there is evidence for the presence of gas-phase methane in situ. This study used a one-dimensional (vertical) numerical model of the subsurface down to a depth of 844 m beneath the Bruce site to evaluate possible effects of hydromechanical coupling with multiphase flow on pressure evolution during glacial loading and unloading. The simulated pressure conditions were affected strongly by the amount of methane initially present in the system, and the maximum simulated underpressure varied nonmonotonically with increasing initial methane content. When the initial methane content was below the solubility limit, exsolution led to underpressures that briefly exceeded those that formed in the single-phase case. At intermediate initial methane contents (sufficient to produce an immobile gas phase), the gas phase dampened the hydromechanical effects of the glacial cycle. At large initial methane contents (when a mobile gas phase was present), gas migration caused a large decrease in relative liquid permeability, which further contributed to underpressure development in the pore water. Multiple scenarios that spanned a range of initial methane contents yielded underpressures like those observed at the Bruce site.

Journal of Geophysical Research - Solid Earth↗
USGS2022

Detection of aseismic slip and poroelastic reservoir deformation at the North Brawley Geothermal Field from 2009 to 2019

The North Brawley Geothermal Field, located within the Brawley Seismic Zone of Southern California, presents a case study for understanding seismic hazards linked to fluid injection and geothermal energy extraction. An earthquake swarm near the geothermal field in 2012 included two earthquakes with magnitudes greater than 5 and was potentially preceded by a years-long aseismic slip transient. To better understand ground deformation around the geothermal field, including its evolution with time and its physical mechanisms, we analyze deformation before, during, and after the swarm using ground- and satellite-based geodetic techniques between 2009 and 2019. We integrate observations from GNSS, Sentinel-1, TerraSAR-X, UAVSAR, and leveling surveys into a single deformation history. Modeling of this new collection of observations at the North Brawley Geothermal Field provides evidence for 80% more pre-swarm aseismic slip than previously recognized from 2009 to 2012. During the 2012 Brawley swarm, our geodetic slip inversions closely match the results of seismic waveform inversions from the swarm events. After the 2012 swarm, surface deformation is dominated by poroelastic deformation of a shallow fluid reservoir at <1 km depth rather than fault slip. The deformation history and seismicity catalogs at North Brawley suggest a cessation of fault-related slip during the ∼7 years after the 2012 earthquake swarm.

California↗
USGS2022

Quantifying large-scale surface change using SAR amplitude images: Crater morphology changes during the 2019-2020 Shishaldin Volcano eruption

Morphological processes often induce meter-scale elevation changes. When a volcano erupts, tracking such processes provides insights into the style and evolution of eruptive activity and related hazards. Compared to optical remote-sensing products, synthetic aperture radar (SAR) observes surface change during inclement weather and at night. Differential SAR interferometry estimates phase change between SAR acquisitions and is commonly applied to quantify deformation. However, large deformation or other coherence loss can limit its use. We develop a new approach applicable when repeated digital elevation models (DEMs) cannot be otherwise retrieved. Assuming an isotropic radar cross-section, we estimate meter-scale vertical morphological change directly from SAR amplitude images via an optimization method that utilizes a high-quality DEM. We verify our implementation through simulation of a collapse feature that we modulate onto topography. We simulate radar effects and recover the simulated collapse. To validate our method, we estimate elevation changes from TerraSAR-X stripmap images for the 2011–2012 eruption of Mount Cleveland. Our results reproduce those from two previous studies; one that used the same dataset, and another based on thermal satellite data. By applying this method to the 2019–2020 eruption of Shishaldin Volcano, Alaska, we generate elevation change time series from dozens of co-registered TerraSAR-X high-resolution spotlight images. Our results quantify previously unresolved cone growth in November 2019, collapses associated with explosions in December–January, and further changes in crater elevations into spring 2020. This method can be used to track meter-scale morphology changes for ongoing eruptions with low latency as SAR imagery becomes available.

Alaska↗
USGS2022

Magnetotelluric investigations of the Kīlauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kīlauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kīlauea and the East and Southwest Rift zones, which has implications for understanding Kīlauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kīlauea. Melt fractions between 0.096 and 0.117 are calculated for the Kīlauea and Puʻuʻōʻō low resistivity anomalies, respectively.

Hawai'i↗
USGS2022

Stress heterogeneity as a driver of aseismic slip during the 2011 Prague, Oklahoma aftershock sequence

The interaction of aseismic and seismic slip before and after an earthquake is fundamental for both earthquake nucleation and postseismic stress relaxation. However, it can be difficult to determine where and when aseismic slip occurs within the seismogenic zone because geodetic techniques are limited to detecting moderate to large slip amplitudes or long duration small slip amplitudes. Here, we use repeating earthquakes (earthquakes that re-rupture the same fault patch) as a proxy for aseismic slip during the 2011 Prague, Oklahoma earthquake sequence. We find that aseismic slip in the Prague earthquake sequence occurs both within the granitic basement and the overlying sedimentary rocks. The repeating earthquakes show that patches of aseismic slip are mostly located at fault intersections. These fault intersections hosted possible mainshock slip, abundant aftershocks, and afterslip. We estimate that ∼40% of the aftershocks are driven by afterslip. We interpret that aseismic slip occurs at fault intersections where stress heterogeneity creates patches of lower stress that are stable within a nonsteady state, rate-state framework.

Oklahoma↗
USGS2022

Off-fault deformation in regions of complex fault geometries: the 2013, Mw7.7, Baluchistan rupture (Pakistan)

Observations of recent earthquake surface ruptures show that ground deformations include a localized component occurring on faults, and an off-fault component affecting the surrounding medium. This second component is also referred to as off-fault deformation (OFD). The localized component generally occurs on complex networks of faults that connect at depth onto a unique fault plane, whereas OFD consists of distributed fracturing and diffuse deformation of the bulk volume, and occurs over scales of hundreds of meters to kilometers around the faults. High-resolution optical image correlation presents a unique potential to characterize the complexity of the surface displacements, including on-fault displacements and OFDs. In this study, we used sub-pixel correlation of 0.5-m resolution optical images to measure the surface displacement field with a <20 cm accuracy for a 30-km long section of the 2013 M w 7.7 Baluchistan, Pakistan, rupture. Our results document significant variability in the fault displacements, associated with large proportions of OFD in regions of fault geometrical complexity. Conversely, in regions where the fault geometry is simple, surface deformation is entirely accommodated by the primary faults with 0% OFD. When combining the localized deformation on faults with the OFD, we show that the total surface displacement budget is constant along the strike of the rupture, despite strong variations observed in the rupture geometry. Based on this analysis, we propose an idealized scenario of earthquake surface deformation as a function of the rupture geometrical variations.

Journal of Geophysical Research, Solid Earth↗
USGS2023

Spatial and temporal analysis of geologic slip rates, Cucamonga Fault, California, USA: Implications for along-strike applications and multi-fault rupture

To constrain fault processes and hazard, fault slip rates may be extrapolated over different fault lengths or time intervals. Here, we investigate slip rates for the Cucamonga Fault (CF). The CF is located at the junction of the Transverse Range fault system with the San Andreas and San Jacinto Faults, and it is hypothesized to connect with these faults, promoting the propagation of large, multi-fault earthquakes. Previous work has shown that CF displacements on late Quaternary alluvial fan surfaces are highly variable along strike. We present two new 10 Be surface exposure ages from depth profiles on the alluvial fans. Slip rates are consistent with a rate of 1.4 ± 0.3 m/kyr over time intervals of ∼20, ∼30, and ∼40 kyr. If the CF participates in multi-fault ruptures, then these earthquakes occur either rarely or with sufficient regularity to maintain apparently steady rates over multiple intervals. We also explore along-strike fault displacement variability using a calibrated morphological model. The model successfully reproduces scarp profiles and indicates that fault displacement variability can be explained in part by scarp age but not uplift rate. We infer that both erosion by ephemeral gullying and distributed deformation contribute to fault displacement variability, although both are difficult to detect confidently without excavations across the scarp. These investigations show that better characterization of cumulative-slip variability along strike may improve accuracy and precision of slip rates. Slip rates that do not consider epistemic uncertainties may not be suitable for extrapolation over longer fault sections.

California↗
USGS2022

Toward next-generation lava flow forecasting: Development of a fast, physics-based lava propagation model

During effusive volcanic crises, the eruption and propagation of lava flows pose a significant hazard to nearby populations, homes, and infrastructure. Consequently, timely lava flow forecasts are a critical need for volcano observatory and emergency management operations. Previous lava flow modeling tools are typically either too slow to produce timely forecasts, or are fast, but lack critical aspects of lava physics or important forecasting outputs. In particular, the strong thermal stratification present in laminar, high-Prandtl number flows has generally been neglected. Bulk rheological changes have previously been computed from cell-averaged temperatures, assuming that the flow is thermally mixed. Here, we detail the development and initial testing of Lava2d, a new two-dimensional depth-averaged finite volume model of lava flow propagation over natural terrain which accounts for bulk rheological changes due to thermorheological stratification. We use a novel approach to energy conservation based on tracking cooling and solidifying at the flow base and at the moving flow surface, allowing for the estimation of more realistic vertical thermorheological profiles, while maintaining computational efficiency, producing very timely model runs. We validate our approach with three examples: comparison with theoretical propagation of crust-dominated lava flows, comparison with a large-scale molten basalt experiment from the Syracuse University Lava Project, and efficiency testing and comparison with the initial phase of the 1984 Mauna Loa lava flows. Our model is shown to produce rapid, realistic forecasts, making it a good candidate for operationalization in active volcanic regions such as in Hawai'i.

Journal of Geophysical Research - Solid Earth↗
USGS2022

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗
USGS2023

Using corrected and imputed polarity measurements to improve focal mechanisms in a regional earthquake catalog near the Mt. Lewis Fault Zone, California

We utilized relative polarity measurements and machine learning techniques to better resolve focal mechanisms and stress orientations considering a catalog of ∼29,000 relocated earthquakes that occurred during 1984–2021 in the southeastern San Francisco Bay Area. Earthquake focal mechanisms are commonly produced using P wave first motion polarities, which traditionally requires events to be well-recorded across a seismic network with good focal sphere coverage. We adapted recently developed approaches that are less dependent on high signal-to-noise records and exploit similar waveforms to produce relative polarity and amplitude measurements between earthquake pairs. These techniques were previously only applied on localized earthquake sequences, and we further developed these approaches so that they can be utilized for regional catalogs. We validated or corrected manually identified polarities by performing polarity consensuses using earthquake pairs. Missing and unreliable polarity measurements were imputed using iterative random forests, an unsupervised ensemble machine learning method. Relative P and S wave amplitude measurements were made between earthquakes, constraining S / P ratios for low signal-to-noise waveforms. Using these techniques, we were able to reduce focal mechanism uncertainties by an average of ∼13° and produced well-constrained focal mechanisms for ∼6 times as many earthquakes than those produced using only the traditionally derived polarities. We performed stress inversions using the focal mechanisms by grouping the focal mechanism results into a quadtree structure. Our stress results are consistent with previous work, albeit at a higher spatial resolution, and demonstrate these techniques can aid our understanding of fault structures and kinematics in more detail than was previously possible.

California↗
USGS2023

Strength recovery in quartzite is controlled by changes in friction in experiments at hydrothermal conditions up to 200°C

The rate of fault zone restrengthening between earthquakes can be influenced by both frictional and cohesive healing processes. Friction is dependent on effective normal stress while cohesion is independent of normal stress, potentially explaining—in part—the lack of depth dependence of earthquake stress drops. Although amenable to laboratory testing, few studies have systematically addressed the normal stress dependence of restrengthening rate. This is partially due to difficulty in separating relative contributions of friction and cohesion in recovery of fault strength. We present results from a series of slide-hold-slide tests on thin layers (≤10 𝜇m) of ultrafine quartz gouge that develop during shearing of initially bare-surface quartzite. Tests were conducted at 10 MPa constant pore pressure, 20–200 MPa constant effective normal stress, and temperatures of 22°–200°C. Restrengthening, defined as the difference between peak shear stress measured after resumption of sliding and steady-state sliding shear stress, increases with the log of hold duration. The 200°C healing rate, 0.014 per e-fold increase in time, is comparable to that determined from seismological observations along the Calaveras Fault, California. Construction of Mohr-Coulomb failure envelopes shows that changes in cohesion are small (<1 MPa) and independent of hold durations to 10 5 s, indicating that the increased strength is due to changes in the friction coefficient. These experimental results are inconsistent with the hypothesis that cohesive healing explains the depth independence of earthquake stress drop, but higher temperatures, longer time-scales, and more complex mineralogy could facilitate cohesive healing in natural fault systems.

Journal of Geophysical Research - Solid Earth↗
USGS2023

Accuracy of finite fault slip estimates in subduction zone regions with topographic Green's functions and seafloor geodesy

Until recently, the lack of seafloor geodetic instrumentation and the use of unrealistically simple, half-space based forward models have resulted in poor resolution of near-trench slip in subduction zone settings. Here, we use a synthetic framework to investigate the impact of topography and geodetic data distribution on coseismic slip estimates in various subduction zone settings. We calculate surface displacements in two synthetic topographic domains that have topography similar to that of Chile and Japan, respectively. We then attempt to image target slip distributions by using a Bayesian approach to solve for slip with two sets of Green's functions—one that accounts for topography and one that does not—and five sets of 50 or more observation points selected from the synthetic surface displacements. Three of these sets of observation points are entirely onland, and two include 5–10 seafloor geodetic sites. We find that the use of topographic Green's functions always improves inferred slip models, and with seafloor geodetic data, it enables an almost perfect recovery of a target slip model, even in the near-trench region. Critically, our results demonstrate that it would be impossible for non-topographic Green's functions to properly recover the true slip distribution, particularly in the near-trench region. We also perform a parameter study with approximately 4,000 slip models estimated using a least-square approach, and find that topographic Green's functions yield significantly more accurate slip models in cases where good data (well distributed and reasonably dense) are available, even in the absence of seafloor geodetic sites.

Journal of Geophysical Research: Solid Earth↗
USGS2023

Pre-eruptive outgassing and pressurization, and post-fragmentation bubble nucleation, recorded by vesicles in breadcrust bombs from vulcanian activity at Guagua Pichincha Volcano, Ecuador

Breadcrust bombs formed during Vulcanian eruptions are assumed to originate from the shallow plug or dome. Their rim to core texture reflects the competition between cooling and degassing timescales, which results in a dense crust with isolated vesicles contrasting with a highly vesicular vesicle network in the interior. Due to relatively fast quenching, the crust can shed light on pre- and syn-eruptive conditions prior to or during fragmentation, whereas the interior allows us to explore post-fragmentation vesiculation. Investigation of pre- to post-fragmentation processes in breadcrust bombs from the 1999 Vulcanian activity at Guagua Pichincha, Ecuador, via 2D and 3D textural analysis reveals a complex vesiculation history, with multiple, spatially localized nucleation and growth events. Large vesicles (Type 1), present in low number density in the crust, are interpreted as pre-eruptive bubbles formed by outgassing and collapse of a permeable bubble network during ascent or stalling in the plug. Haloes of small, syn-fragmentation vesicles (Type 2), distributed about large vesicles, are formed by pressurization and enrichment of volatiles in these haloes. The nature of the pressurization process in the plug is discussed in light of seismicity and ground deformation signals, and previous textural and chemical studies. A third population (Type 3) of post-fragmentation small vesicles appears in the interior of the bomb, and growth and coalescence of Type 2 and 3 vesicles causes the transition from isolated to interconnected bubble network in the interior. We model the evolution of viscosity, bubble growth rate, diffusion timescales, bubble radius and porosity during fragmentation and cooling. These models reveal that thermal quenching dominates in the crust whereas the interior undergoes a viscosity quench caused by degassing, and that the transition from crust to interior corresponds to the onset of percolation and development of permeability in the bubble network.

Guagua Pichincha Volcano↗
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↗