SEARCH · Search USGS

Results for “Bulletin of the Geological Society of America”

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 901 records · Page 50Linked to original sources

USGS2003

Salton Trough regional deformation estimated from combined trilateration and survey-mode GPS data

The Salton Trough in southeastern California, United States, has one of the highest seismicity and deformation rates in southern California, including 20 earthquakes M 6 or larger since 1892. From 1972 through 1987, the U.S. Geological Survey (USGS) measured a 41-station trilateration network in this region. We remeasured 37 of the USGS baselines using survey - mode Global Positioning System methods from 1995 through 1999. We estimate the Salton Trough deformation field over a nearly 30-year period through combined analysis of baseline length time series from these two datasets. Our primary result is that strain accumulation has been steady over our observation span, at a resolution of about 0.05 μstrain/yr at 95% confidence, with no evidence for significant long-term strain transients despite the occurrence of seven large regional earthquakes during our observation period. Similar to earlier studies, we find that the regional strain field is consistent with 0.5 ± 0.03 μstrain/yr total engineering shear strain along an axis oriented 311.6° ± 23° east of north, approximately parallel to the strike of the major regional faults, the San Andreas and San Jacinto (all uncertainties in the text and tables are standard deviations unless otherwise noted). We also find that (1) the shear strain rate near the San Jacinto fault is at least as high as it is near the San Andreas fault, (2) the areal dilatation near the southeastern Salton Sea is significant, and (3) one station near the southeastern Salton Sea moved anomalously during the period 1987.95-1995.11.

California↗

USGS2006

Modeling and validation of a 3D velocity structure for the Santa Clara Valley, California, for seismic-wave simulations

A 3D seismic velocity and attenuation model is developed for Santa Clara Valley, California, and its surrounding uplands to predict ground motions from scenario earthquakes. The model is developed using a variety of geologic and geophysical data. Our starting point is a 3D geologic model developed primarily from geologic mapping and gravity and magnetic surveys. An initial velocity model is constructed by using seismic velocities from boreholes, reflection/refraction lines, and spatial autocorrelation microtremor surveys. This model is further refined and the seismic attenuation is estimated through waveform modeling of weak motions from small local events and strong-ground motion from the 1989 Loma Prieta earthquake. Waveforms are calculated to an upper frequency of 1 Hz using a parallelized finite-difference code that utilizes two regions with a factor of 3 difference in grid spacing to reduce memory requirements. Cenozoic basins trap and strongly amplify ground motions. This effect is particularly strong in the Evergreen Basin on the northeastern side of the Santa Clara Valley, where the steeply dipping Silver Creek fault forms the southwestern boundary of the basin. In comparison, the Cupertino Basin on the southwestern side of the valley has a more moderate response, which is attributed to a greater age and velocity of the Cenozoic fill. Surface waves play a major role in the ground motion of sedimentary basins, and they are seen to strongly develop along the western margins of the Santa Clara Valley for our simulation of the Loma Prieta earthquake.

California↗

USGS2023

Applications of nonergodic site response models to ShakeAlert case studies in the Los Angeles area

In this study, we explore whether the Parker and Baltay (2022) site response models for the Los Angeles (LA) basin region can improve ground‐motion forecasts in the U.S. Geological Survey ShakeAlert earthquake early warning system (hereafter ShakeAlert). We implement the peak ground acceleration and peak ground velocity site response models of Parker and Baltay (2022) in ShakeAlert via the earthquake information to ground‐motion (hereafter eqinfo2GM) module, which predicts ground motions from the estimated earthquake parameters of magnitude, rupture length, and location. The nonergodic site response models for the greater LA area were developed using ground motions from 414 M 3–7.3 earthquakes in southern California. We test nonergodic ground‐motion forecasts for five earthquakes in the LA area: the 1994 M 6.7 Northridge earthquake, the 2008 M 5.4 Chino Hills earthquake, the 2019 M 7.1 Ridgecrest earthquake, the 2020 M 4.5 South El Monte earthquake, and a synthetic M 7.8 earthquake on the southern San Andreas fault from the ShakeOut scenario, which was the basis of a statewide emergency response exercise. From the test results, we find that with the nonergodic site response applied, ShakeAlert not only alerts larger areas but can also result in longer warning times in LA region. In addition, the modified Mercalli intensity (MMI) ground‐motion predictions generated by the ShakeAlert eqinfo2GM module are improved in accuracy when compared with the corresponding ShakeMap ground‐truth MMI when the nonergodic site response model is applied.

California↗

USGS2015

Basement and regional structure along strike of the Queen Charlotte Fault in the context of modern and historical earthquake ruptures

The Queen Charlotte fault (QCF) is a dextral transform system located offshore of southeastern Alaska and western Canada, accommodating &sim;4.4  cm/yr of relative motion between the Pacific and North American plates. Oblique convergence along the fault increases southward, and how this convergence is accommodated is still debated. Using seismic reflection data, we interpret offshore basement structure, faulting, and stratigraphy to provide a geological context for two recent earthquakes, an M w 7.5 strike‐slip event near Craig, Alaska, and an M w 7.8 thrust event near Haida Gwaii, Canada. We map downwarped Pacific oceanic crust near 54° N, between the two rupture zones. Observed downwarping decreases north and south of 54° N, parallel to the strike of the QCF. Bending of the Pacific plate here may have initiated with increased convergence rates due to a plate motion change at &sim;6  Ma. Tectonic reconstruction implies convergence‐driven Pacific plate flexure, beginning at 6 Ma south of a 10° bend the QCF (which is currently at 53.2° N) and lasting until the plate translated past the bend by &sim;2  Ma. Normal‐faulted approximately late Miocene sediment above the deep flexural depression at 54° N, topped by relatively undeformed Pleistocene and younger sediment, supports this model. Aftershocks of the Haida Gwaii event indicate a normal‐faulting stress regime, suggesting present‐day plate flexure and underthrusting, which is also consistent with reconstruction of past conditions. We thus favor a Pacific plate underthrusting model to initiate flexure and accommodation space for sediment loading. In addition, mapped structures indicate two possible fault segment boundaries along the QCF at 53.2° N and at 56° N.

Alaska↗

USGS1998

Global teleseismic earthquake relocation with improved travel times and procedures for depth determination

We relocate nearly 100,000 events that occurred during the period 1964 to 1995 and are well-constrained teleseismically by arrival-time data reported to the International Seismological Centre ( ISC ) and to the U.S. Geological Survey's National Earthquake Information Center ( NEIC ). Hypocenter determination is significantly improved by using, in addition to regional and teleseismic P and S phases, the arrival times of PKiKP, PKPdf , and the teleseismic depth phases pP, pwP , and sP in the relocation procedure. A global probability model developed for later-arriving phases is used to independently identify the depth phases. The relocations are compared to hypocenters reported in the ISC and NEIC catalogs and by other sources. Differences in our epicenters with respect to ISC and NEIC estimates are generally small and regionally systematic due to the combined effects of the observing station network and plate geometry regionally, differences in upper mantle travel times between the reference earth models used, and the use of later-arriving phases. Focal depths are improved substantially over most other independent estimates, demonstrating (for example) how regional structures such as downgoing slabs can severely bias depth estimation when only regional and teleseismic P arrivals are used to determine the hypocenter. The new data base, which is complete to about Mw 5.2 and includes all events for which moment-tensor solutions are available, has immediate application to high-resolution definition of Wadati-Benioff Zones ( WBZ s) worldwide, regional and global tomographic imaging, and other studies of earth structure.

Bulletin of the Seismological Society of America↗

USGS2020

Surface displacement distributions for the July 2019 Ridgecrest, California earthquake ruptures

Surface rupture in the 2019 Ridgecrest, California, earthquake sequence occurred along two orthogonal cross faults and includes dominantly left‐lateral and northeast‐striking rupture in the M w "> M w 6.4 foreshock and dominantly right‐lateral and northwest‐striking rupture in the M w "> M w 7.1 mainshock. We present > 650 "> > 650 field‐based, surface‐displacement observations for these ruptures and synthesize our results into cumulative along‐strike displacement distributions. Using these data, we calculate displacement gradients and compare our results with historical strike‐slip ruptures in the eastern California shear zone. For the M w "> M w 6.4 rupture, we report 96 displacements measured along 18 km of northeast‐striking rupture. Cumulative displacement curves for the rupture yield a mean left‐lateral displacement of 0.3–0.5 m and maximum of 0.7–1.6 m. Net mean vertical displacement based on the difference of down‐to‐the‐west (DTW) and down‐to‐the‐east (DTE) displacement curves is close to zero (0.02 m DTW). The M w "> M w 6.4 displacement distribution shows that the majority of displacement occurred southwest of the intersection with the M w "> M w 7.1 rupture. The M w "> M w 7.1 rupture is northwest‐striking and 50 km long based on 576 field measurements. Displacement curves indicate a mean right‐lateral displacement of 1.2–1.7 m and a maximum of 4.3–7.0 m. Net vertical displacement in the rupture averages 0.3 m DTW. The M w "> M w 7.1 displacement distributions demonstrate that maximum displacement occurred along a 12‐km‐long portion of the fault near the M w "> M w 7.1 epicenter, releasing 66% of the geologically based seismic moment along 24% of the total rupture length. Using our displacement distributions, we calculate kilometer‐scale displacement gradients for the M w "> M w 7.1 rupture. The steepest gradients ( ⁠ ∼ 1 – 3    m / km "> ∼ 1 – 3 m / km ⁠ ) flank the 12‐km‐long region of maximum displacement. In contrast, gradients for the 1992 M w "> M w 7.3 Landers and 1999 M w "> M w 7.1 Hector Mine earthquakes are < 0.6    m / km "> < 0.6 m / km . Our displacement distributions are important for understanding the influence of cross‐fault rupture on M w "> M w 6.4 and 7.1 rupture length and displacement and will facilitate comparisons with distributions generated remotely and at broader scales.

California↗

USGS2006

San Andreas fault geometry in the Parkfield, California, region

In map view, aftershocks of the 2004 Parkfield earthquake lie along a line that forms a straighter connection between San Andreas fault segments north and south of the Parkfield reach than does the mapped trace of the fault itself. A straightedge laid on a geologic map of Central California reveals a ???50-km-long asymmetric northeastward warp in the Parkfield reach of the fault. The warp tapers gradually as it joins the straight, creeping segment of the San Andreas to the north-west, but bends abruptly across Cholame Valley at its southeast end to join the straight, locked segment that last ruptured in 1857. We speculate that the San Andreas fault surface near Parkfield has been deflected in its upper ???6 km by nonelastic behavior of upper crustal rock units. These units and the fault surface itself are warped during periods between large 1857-type earthquakes by the presence of the 1857-locked segment to the south, which buttresses intermittent coseismic and continuous aseismic slip on the Parkfield reach. Because of nonelastic behavior, the warping is not completely undone when an 1857-type event occurs, and the upper portion of the three-dimensional fault surface is slowly ratcheted into an increasingly prominent bulge. Ultimately, the fault surface probably becomes too deformed for strike-slip motion, and a new, more vertical connection to the Earth's surface takes over, perhaps along the Southwest Fracture Zone. When this happens a wedge of material currently west of the main trace will be stranded on the east side of the new main trace.

Bulletin of the Seismological Society of America↗

USGS2022

Rupture scenarios for the 3 June 1770 Haiti earthquake

The 2010 M 7.0 Haiti earthquake provided the impetus to reconsider historical earthquakes in Hispaniola ( Bakun et al. , 2012 ). That earthquake also shed new light on complex fault systems along Haiti’s southern peninsula ( Douilly et al. , 2013 ; Saint Fleur et al. , 2015 ). Recently, the 2021 M 7.2 Nippes earthquake ( Calais et al. , 2022 ; Douilly et al. , 2022 ), and a recent study reconsidering the 1860 sequence ( Martin et al. , 2022 ) further underscored the complexity of fault systems and large earthquake ruptures along the peninsula. Motivated by these studies and recent geological investigations ( Prentice et al. , 2010 ; Saint Fleur et al. , 2020 ), we reconsider the 3 June 1770 Haiti earthquake to explore the conventional assumption that it was the last major ( M ≥7.5) earthquake along the Enriquillo–Plantain Garden fault (EPGF). Accounts provide compelling evidence for substantial liquefaction in the Cul‐de‐Sac plain, one or more likely landslide‐driven tsunami in Gonaïves Bay, and extensive landsliding that created at least three documented landslide dams. We consider three end‐member rupture scenarios that are consistent with available constraints: two scenarios with M 7.7 and rupture lengths of 150–170 km, and one scenario with a ∼90 km rupture and M 7.5. Absent future work to identify and date paleoevents along the southern peninsula, none of these scenarios can be ruled out. Our preferred rupture model extends from the Miragoâne pull‐apart to near la Selle mountain, with a rupture length of 127 km, M 7.6, and a high stress drop. Rupture could have been on the EPGF or on an oblique thrust fault associated with overthrusting of the Massif de la Selle. The results do support the conclusion that the 1770 earthquake was the last major earthquake in southern Haiti, with a magnitude upward of M 7.5 and significantly more severe shaking in southern Haiti than during the 2010 earthquake.

Bulletin of the Seismological Society of America↗

USGS2004

Revisiting the 23 February 1892 Laguna Salada earthquake

According to some compilations, the Laguna Salada, Baja California, earthquake of 23 February 1892 ranks among the largest earthquakes in California and Baja California in historic times. Although surface rupture was not documented at the time of the earthquake, recent geologic investigations have identified and mapped a rupture on the Laguna Salada fault that can be associated with high probability with the 1892 event (Mueller and Rockwell, 1995). The only intensity-based magnitude estimate for the earthquake, M 7.8, was made by Strand (1980) based on an interpretation of macroseismic effects and a comparison of isoseismal areas with those from instrumentally recorded earthquakes. In this study we reinterpret original accounts of the Laguna Salada earthquake. We assign modified Mercalli intensity (MMI) values in keeping with current practice, focusing on objective descriptions of damage rather than subjective human response and not assigning MMI values to effects that are now known to be poor indicators of shaking level, such as liquefaction and rockfalls. The reinterpreted isoseismal contours and the estimated magnitude are both significantly smaller than those obtained earlier. Using the method of Bakun and Wentworth (1997) we obtain a magnitude estimate of M 7.2 and an optimal epicenter less than 15 km from the center of the mapped Laguna Salada rupture. The isoseismal contours are elongated toward the northwest, which is qualitatively consistent with a directivity effect, assuming that the fault ruptured from southeast to northwest. We suggest that the elongation may also thus reflect wave propagation effects, with more efficient propagation of crustal surface (Lg) waves in the direction of the overall regional tectonic fabric.

Bulletin of the Seismological Society of America↗

USGS2000

Liquefaction evidence for the strength of ground motions resulting from Late Holocene Cascadia subduction earthquakes, with emphasis on the event of 1700 A.D.

During the past decade, paleoseismic studies done by many researchers in the coastal regions of the Pacific Northwest have shown that regional downdropping and subsequent tsunami inundation occurred in response to a major earthquake along the Cascadia subduction zone. This earthquake occurred almost certainly in 1700 A.D., and is believed by many to have been of M 8.5-9 or perhaps larger. In order to characterize the severity of ground motions from this earthquake, we report on a field search and analysis of seismically induced liquefaction features. The search was conducted chiefly along the banks of islands in the lowermost Columbia River of Oregon and Washington and in stream banks along smaller rivers throughout southwestern Washington. To a lesser extent, the investigation included rivers in central Oregon. Numerous small- to moderate-sized liquefaction features from the earthquake of 1700 A.D. were found in some regions, but there was a notable lack of liquefaction features in others. The regional distribution of liquefaction features is evaluated as a function of geologic and geotechnical factors in different field settings near the coast. Our use of widely different field settings, each in which we independently assess the strength of shaking and arrive at the same conclusion, enhances the credibility of our interpretations. Our regional inventory of liquefaction features and preliminary geotechnical analysis of liquefaction potential provide substantial evidence for only moderate levels of ground shaking in coastal Washington and Oregon during the subduction earthquake of 1700 A.D. Additionally, it appears that a similar conclusion can be reached for an earlier subduction earthquake that occurred within the past 1100 years, which also has been characterized by others as being M 8 or greater. On the basis of more limited data for older events collected in our regional study, it appears that seismic shaking has been no stronger throughout Holocene time. Our interpreted levels of shaking are considerably lower than current estimates in the technical literature that use theoretical and statistical models to predict ground motions of subduction earthquakes in the Cascadia region. Because of the influence of estimated ground motions from Cascadia subduction-zone earthquakes on seismic hazard evaluations, more paleoliquefaction and geotechnical field studies are needed to definitively bracket the strength of shaking. With further work, it should be possible to extend the record of seismic shaking through much of Holocene time in large portions of Washington and Oregon.

Oregon, Washington↗

USGS2020

Regionally Optimized Background Earthquake Rates from ETAS (ROBERE) for probabilistic seismic hazard assessment

We use an epidemic‐type aftershock sequence (ETAS) based approach to develop a regionally optimized background earthquake rates from ETAS (ROBERE) method for probabilistic seismic hazard assessment. ROBERE fits parameters to the full seismicity catalog for a region with maximum‐likelihood estimation, including uncertainty. It then averages the earthquake rates over a suite of catalogs from which foreshocks and aftershocks have been removed using stochastic declustering while maintaining the same Gaussian smoothing currently used for the U.S. Geological Survey National Seismic Hazard Model (NSHM). The NSHM currently determines these rates by smoothing a single catalog from which foreshocks and aftershocks have been removed using the method of Gardner and Knopoff (1974 ; hereafter, GK74). The parameters used in GK74 were determined from subjectively identified aftershock sequences, unlike ROBERE, in which both background rate and aftershock triggering parameters are objectively fitted. A major difference between the impacts of the two methods is GK74 significantly reduces the b "> b ‐value, a critical value for seismic hazard analysis, whereas ROBERE maintains the original b "> b ‐value from the full catalog. We apply these methods to the induced seismicity in Oklahoma and Kansas and tectonic activity in the San Francisco Bay Region. Using GK74 gives lower overall earthquake rates but estimates higher hazard due to the reduction in the b "> b ‐value. ROBERE provides higher earthquake rates, at the magnitude of completeness, but lower hazard because it does not alter the b "> b ‐value. We test two other declustering methods that produce results closer to ROBERE but do not use objectively fit parameters, include uncertainty, and may not work as well in other areas. We suggest adopting ROBERE for the NSHM so that our hazard estimates are based on an objective analysis, including uncertainty, and do not depend strongly on potentially biased b "> b ‐values, which was never the goal of the existing methodology.

Bulletin of the Seismological Society of America↗

USGS2020

Earthquake magnitude and Lg Q variations between the Grenville and northern Appalachian geologic provinces of eastern Canada

This article assesses the ability of regionally specific, frequency‐dependent crustal attenuation ( ⁠ 1 / Q "> 1 / Q ⁠ ) to reduce mean magnitude discrepancies between seismic stations in the northern Appalachian and Grenville provinces (NAP and GP) of Canada. Lg Q ( f ) "> Q ( f ) is an important parameter in ground‐motion models used in probabilistic seismic hazard analysis. Discrepancies in regional magnitude estimates have long been noted to exist between stations in the two provinces for common event origins. Such discrepancies could arise from systematic site condition variations between the geologic provinces or from varying crustal attenuative properties. To evaluate the effect of frequency‐dependent anelastic attenuation, Lg Q ( f ) "> Q ( f ) on estimated magnitudes, we analyze Lg amplitudes from > 6000 "> > 6000 waveforms recorded by Grenville and northern Appalachian receivers from 420 natural earthquakes of M N "> M N magnitude 3–5.6. Waveform analysis is strictly limited to analyst‐reviewed, vertical‐component waveforms in which Lg is clearly identified, ensuring that the datasets exhibit dominant, high‐frequency energy in the Lg velocity window. Lg Q ( f ) "> Q ( f ) is found to be higher in the GP than in the northern Appalachians. In the Grenville, Q ( f ) = 761 ( ± 145 ) f 0.25 ( ± 0.014 ) "> Q ( f ) = 761 ( ± 145 ) f 0.25 ( ± 0.014 ) ⁠ , and in the northern Appalachians, attenuation is higher: Q ( f ) = 506 ( ± 172 ) f 0.33 ( ± 0.310 ) "> Q ( f ) = 506 ( ± 172 ) f 0.33 ( ± 0.310 ) . Earthquake magnitude determined using the peak amplitude of the Lg phase ( ⁠ m b L g "> m b L g ) for eastern Canada is corrected to incorporate the frequency‐dependent, regionally specific Lg Q ( f ) "> Q ( f ) determined in this study. Using the new Lg Q ( f ) "> Q ( f ) values diminishes and nearly resolves magnitude discrepancies between the provinces. Correcting regional magnitude discrepancies between provinces is critical for reliable regional seismic hazard estimates because magnitude error in a particular region could lead to increased uncertainty in seismic hazard models.

Grenville Geologic Province, Northern Appalachian ↗

USGS1985

Characteristics of the aftershock sequence of the Borah Peak, Idaho, earthquake determined from digital recordings of the events

The U.S. Geological Survey, Menlo Park, deployed and maintained a network of twelve digital instruments over the 2 weeks following the 28 October 1983 Borah Peak, Idaho, earthquake. The network recorded 45 events with M ≧ 3.0, and 6 events with M ≦ 4.0. The epicenters are located in a narrow band which parallels the trace of the surface fauiting up to the Willow Creek summit; the depths of the events range from 5 to 16 km. In the south, the distribution of hypocenters delineate a plane which dips to the southwest at 50°; to the north, the hypocenters dip steeply to the east. Composite focal mechanisms for three groups of events show normal faulting mechanisms; the mechanism of the aftershocks in the north appear rotated in both strike and dip from the aftershocks in the south. The seismic moments of the aftershocks increase with increasing hypocentral depth below 12 km. The dynamic stress drops of the events do not show any systematic variation with depth, however. Most of the events with large stress drops are clustered in the northwest limb of the aftershock distribution; the average stress drop of the southern events is 31 ± 16 bars, while the average stress drop of the events in the northwest limb is 77 ± 52 bars. This clustering of events with large stress drops marks an apparent stress concentration, possibly associated with the arrest of the main shock rupture propagation by a fracture barrier at depth.

Idaho↗

USGS2006

Rupture propagation of the 2004 Parkfield, California, earthquake from observations at the UPSAR

Using a short-baseline seismic array (U.S. Geological Survey Parkfield Dense Seismograph Array [UPSAR]) about 12 km west of the rupture initiation of the 28 September 2004 M 6.0 Parkfield, California, earthquake, we have observed the movement of the rupture front of this earthquake on the San Andreas fault. The sources of high-frequency arrivals at UPSAR, which we use to identify the rupture front, are mapped onto the San Andreas fault using their apparent velocity and back azimuth. Measurements of apparent velocity and back azimuth are calibrated using aftershocks, which have a compact source and known location. Aftershock back azimuths show considerable lateral refraction, consistent with a high-velocity ridge on the southwest side of the fault. We infer that the initial mainshock rupture velocity was approximately the Rayleigh speed (with respect to slower side of the fault), and the rupture then slowed to about 0.66?? near the town of Parkfield after 2 sec. The last well-correlated pulse, 4 sec after S, is the largest at UPSAR, and its source is near the region of large accelerations recorded by strong-motion accelerographs and close to northern extent of continuous surface fractures on the southwest fracture zone. Coincidence of sources with preshock and aftershock distributions suggests fault material properties control rupture behavior. High-frequency sources approximately correlate with the edges of asperities identified as regions of high slip derived from inversion of strong-motion waveforms.

Bulletin of the Seismological Society of America↗

USGS2004

Gravity study through the Tualatin Mountains, Oregon: Understanding crustal structure and earthquake hazards in the Portland urban area

A high-resolution gravity survey through the Tualatin Mountains (Portland Nills) west of downtown Portland exhibits evidence of faults previously identified from surface geologic and aeromagnetic mapping. The gravity survey was conducted in 1996 along the 4.5-km length of a twin-bore tunnel, then under construction and now providing light-rail service between downtown Portland and communities west of the Portland Hills. Gravitational attraction gradually increases from west to east inside the tunnel, which reflects the tunnel's location between low-density sedimentary deposits of the Tualatin basin to the west and high-density, mostly concealed Eocene basalt to the east. Superimposed on this gradient are several steplike anomalies that we interpret as evidence for faulted contacts between rocks of contrasting density. The largest of these anomalies occurs beneath Sylvan Creek, where a fault had previously been mapped inside the tunnel. Another occurs 1200 m from the west portal, at the approximate intersection of the tunnel with an aeromagnetic anomaly associated with the Sylvan fault (formerly called the Oatfield fault). Lithologic cross sections based on these gravity data show that the steplike anomalies are consistent with steeply dipping reverse faults, although strike-slip displacements also may be important. Three gravity lows correspond with topographic lows directly overhead and may reflect zones of shearing. Several moderate earthquakes (M ??? 3.5) occurred near the present-day location of the tunnel in 1991, suggesting that some of these faults or other faults in the Portland Hills fault zone are seismically active.

Bulletin of the Seismological Society of America↗

USGS2021

Improved scaling relationships for seismic moment and average slip of strike-slip earthquakes incorporating fault slip rate, fault width and stress drop

We develop a self‐consistent scaling model relating magnitude M w to surface rupture length ( ⁠L E ⁠ ), surface displacement D E ⁠ , and rupture width W E ⁠ , for strike‐slip faults. Knowledge of the long‐term fault‐slip rate S F improves magnitude estimates. Data are collected for 55 ground‐rupturing strike‐slip earthquakes that have geological estimates of L E ⁠ , D E ⁠ , and S F⁠ , and geophysical estimates of W E ⁠ . We begin with the model of Anderson et al. (2017) , which uses a closed form equation for the seismic moment of a surface‐rupturing strike‐slip fault of arbitrary aspect ratio and given stress drop, Δτ C ⁠ . Using W E estimates does not improve M w estimates. However, measurements of D E plus the relationship between Δτ C and surface slip provide an alternate approach to study W E ⁠ . A grid of plausible stress drop and width pairs were used to predict displacement and earthquake magnitude. A likelihood function was computed from within the uncertainty ranges of the corresponding observed M w and D E values. After maximizing likelihoods over earthquakes in length bins, we found the most likely values of W E for constant stress drop; these depend on the rupture length. The best‐fitting model has the surprising form W E ∝logL E —a gentle increase in width with rupture length. Residuals from this model are convincingly correlated to the fault‐slip rate and also show a weak correlation with the crustal thickness. The resulting model thus supports a constant stress drop for ruptures of all lengths, consistent with teleseismic observation. The approach can be extended to test other observable factors that might improve the predictability of magnitude from a mapped fault for seismic hazard analyses.

Bulletin of the Seismological Society of America↗

USGS2017

3D ground‐motion simulations of Mw 7 earthquakes on the Salt Lake City segment of the Wasatch fault zone: Variability of long‐period (T≥1 s) ground motions and sensitivity to kinematic rupture parameters

We examine the variability of long‐period ( T ≥1 s) earthquake ground motions from 3D simulations of M w 7 earthquakes on the Salt Lake City segment of the Wasatch fault zone, Utah, from a set of 96 rupture models with varying slip distributions, rupture speeds, slip velocities, and hypocenter locations. Earthquake ruptures were prescribed on a 3D fault representation that satisfies geologic constraints and maintained distinct strands for the Warm Springs and for the East Bench and Cottonwood faults. Response spectral accelerations (SA; 1.5–10 s; 5% damping) were measured, and average distance scaling was well fit by a simple functional form that depends on the near‐source intensity level SA 0 ( T ) and a corner distance R c :SA( R , T )=SA 0 ( T )(1+( R / R c )) −1 . Period‐dependent hanging‐wall effects manifested and increased the ground motions by factors of about 2–3, though the effects appeared partially attributable to differences in shallow site response for sites on the hanging wall and footwall of the fault. Comparisons with modern ground‐motion prediction equations (GMPEs) found that the simulated ground motions were generally consistent, except within deep sedimentary basins, where simulated ground motions were greatly underpredicted. Ground‐motion variability exhibited strong lateral variations and, at some sites, exceeded the ground‐motion variability indicated by GMPEs. The effects on the ground motions of changing the values of the five kinematic rupture parameters can largely be explained by three predominant factors: distance to high‐slip subevents, dynamic stress drop, and changes in the contributions from directivity. These results emphasize the need for further characterization of the underlying distributions and covariances of the kinematic rupture parameters used in 3D ground‐motion simulations employed in probabilistic seismic‐hazard analyses.

Utah↗

USGS1998

Seismicity alert probabilities at Parkfield, California, revisited

For a decade, the U.S. Geological Survey has used the Parkfield Earthquake Prediction Experiment scenario document to estimate the probability that earthquakes observed on the San Andreas fault near Parkfield will turn out to be foreshocks followed by the expected magnitude 6 mainshocks. During this time, we have learned much about the seismogenic process at Parkfield, about the long-term probability of the Parkfield mainshock, and about the estimation of these types of probabilities. The probabilities for potential foreshocks at Parkfield are reexamined and revised in light of these advances. As part of this process, we have confirmed both the rate of foreshocks before strike-slip earthquakes in the San Andreas physiographic province and the uniform distribution of foreshocks with magnitude proposed by earlier studies. Compared to the earlier assessment, these new estimates of the long-term probability of the Parkfield mainshock are lower, our estimate of the rate of background seismicity is higher, and we find that the assumption that foreshocks at Parkfield occur in a unique way is not statistically significant at the 95% confidence level. While the exact numbers vary depending on the assumptions that are made, the new alert probabilities are lower than previously estimated. Considering the various assumptions and the statistical uncertainties in the input parameters, we also compute a plausible range for the probabilities. The range is large, partly due to the extra knowledge that exists for the Parkfield segment, making us question the usefulness of these numbers.

California↗