Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

Follow this topic with RSS

At least 829 records · Page 46Linked to original sources

USGS1988

The seismic radiation from composite models of faulting

The failure of an asperity, i.e., the dynamic rupture of a small fault area with finite stress drop surrounded by a broken or weak fault area which has no stress drop but which slips after the asperity fails, is proposed as a model for the rupture process of a subevent in a composite earthquake. The rupture area of the composite earthquake surrounding the subevent is modeled by the weak fault area surrounding the asperity in the subevent model. The resulting seismic moment of the subevent is proportional to the stress drop and the rupture area of the subevent, as well as the radius of the composite earthaquake. By setting the stress drops of the asperity models equal to the dynamic stress drops of the subevents, the composite earthquake can be modeled as the sum of a set of subevents which cover the rupture area of the composite earthquake. The scaling of the high- and low-frequency radiation from composite earthquakes composed of asperities is commensurate with generally observed spectral scaling laws, in contrast to composite earthquakes composed of cracks, or smaller earthquakes. A simple filtering strategy is proposed for filtering the waveforms radiated by cracks to approximate the waveforms radiated by asperities. The P and S waves radiated by an ML = 5.2 earthquake which occurred on 9 May 1983, at Coalinga, California, are simulated using the P and S waves radiated by an ML = 3.6 aftershock. The aftershock waveforms are first filtered to approximate the radiation from asperities with the appropriate rupture areas, and then the waveforms of 12 asperity subevents are summed together to simulate the waveforms and spectra of the composite earthquake.

Bulletin of the Seismological Society of America↗

USGS2011

Coulomb 3.3 Graphic-rich deformation and stress-change software for earthquake, tectonic, and volcano research and teaching-user guide

Coulomb is intended both for publication-directed research and for college and graduate school classroom instruction. We believe that one learns best when one can see the most and can explore alternatives quickly. So the principal feature of Coulomb is ease of input, rapid interactive modification, and intuitive visualization of the results. The program has menus and check-items, and dialogue boxes to ease operation. The internal graphics are suitable for publication, and can be easily imported into Illustrator, GMT, Google Earth, or Flash for further enhancements. Coulomb is designed to let one calculate static displacements, strains, and stresses at any depth caused by fault slip, magmatic intrusion, or dike expansion/contraction. One can calculate static displacements (on a surface or at GPS stations), strains, and stresses caused by fault slip, magmatic intrusion, or dike expansion. Problems, such as how an earthquake promotes or inhibits failure on nearby faults, or how fault slip or dike expansion will compress a nearby magma chamber, are germane to Coulomb. Geologic deformation associated with strike-slip faults, normal faults, or fault-bend folds is also a useful application. Calculations are made in an elastic halfspace with uniform isotropic elastic properties following Okada (1992)-see list of key papers in section 1.9. There is substantial evidence to support the hypothesis that faults interact by the transfer of stress, both on the time scales of earthquake sequences and aftershocks and on longer time scales associated with the interevent time of the largest shocks that occur in a given region. There is also evidence that faults and magmatic systems interact as well, and that static stress changes influence intrusions and eruptions. Processes not included in Coulomb are also important, such as dynamic stresses, pore-fluid diffusion, and viscoelastic rebound. Further, elastic stiffness differences between basins and crustal layering modify the stresses in comparison to the elastic halfspace implemented in Coulomb. Nevertheless, we believe that a simple tool that permits exploration of a key component of earthquake interaction has great value for understanding and discovery. To download the software (a MATLAB application) and tutorial files, please go to http://www.coulombstress.org

Open-File Report↗

USGS2005

Three-dimensional geologic map of the Hayward fault, northern California: Correlation of rock unites with variations in seismicity, creep rate, and fault dip

In order to better understand mechanisms of active faults, we studied relationships between fault behavior and rock units along the Hayward fault using a three-dimensional geologic map. The three-dimensional map-constructed from hypocenters, potential field data, and surface map data-provided a geologic map of each fault surface, showing rock units on either side of the fault truncated by the fault. The two fault-surface maps were superimposed to create a rock-rock juxtaposition map. The three maps were compared with seismicity, including aseismic patches, surface creep, and fault dip along the fault, by using visuallization software to explore three-dimensional relationships. Fault behavior appears to be correlated to the fault-surface maps, but not to the rock-rock juxtaposition map, suggesting that properties of individual wall-rock units, including rock strength, play an important role in fault behavior. Although preliminary, these results suggest that any attempt to understand the detailed distribution of earthquakes or creep along a fault should include consideration of the rock types that abut the fault surface, including the incorporation of observations of physical properties of the rock bodies that intersect the fault at depth. ?? 2005 Geological Society of America.

Geology↗

USGS1975

A review of recently active faults in Taiwan

Six faults associated with five large earthquakes produced surface displacements ranging from 1 to 3 m in the period 1906 through 1951. Four of the ruptures occurred in the western coastal plain and foothills, and two occurred in the Longitudinal Valley of eastern Taiwan. Maps are included showing the locations and dimensions of the displacements. The published geological literature probably would not lead one to infer the existence of a fault along most of the 1906 rupture, except for descriptions of the rupture itself. Over most of its length the 1935 rupture on the Chihhu fault is parallel to but more than 0.5 km from nearby faults shown on geologic maps published in 1969 and 1971; only about 1.5 km of its 15 km length coincides with a mapped fault. The coastal plain part of the Tuntzuchio fault which ruptured in 1935 is apparently not revealed by landforms, and only suggested by other data. Part of the 1946 Hsinhua faulting coincides with a fault identified in the subsurface by seismic work but surface indications of the fault are obscure. The 1951 Meilun faulting occurred along a conspicuous pre-1951 scarp and the 1951 Yuli faulting occurred near or in line with pre-1951 scarps. More than 40 faults which, according to the published literature, have had Pleistocene or later movement are shown on a small-scale map. Most of these faults are in the densely-populated western part of Taiwan. The map and text calls attention to faults that may be active and therefore may be significant in planning important structures. Equivocal evidence suggestive of fault creep was found on the Yuli fault and the Hsinhua fault. Fault creep was not found at several places examined along the 1906 fault trace. Tectonic uplift has occurred in Taiwan in the last 10,000 years and application of eustatic sea level curves to published radiocarbon dates shows that the minimum rate of uplift is considerably different in different parts of the island. Incomplete data indicate that the rate is high near Hualien, where an uplift of at least 0.6 m and probably more than 1 m occurred in the 1951 earthquake, and near and south of the 1946 faulting. Sudden uplifts can have serious consequences for installations near the shore. Investigation of this process, study of recently active faults, and continuing study of seismicity are necessary parts of a practical earthquake-hazard reduction program.

Open-File Report↗

USGS2013

The northwest trending north Boquerón Bay-Punta Montalva Fault Zone; A through going active fault system in southwestern Puerto Rico

The North Boquerón Bay–Punta Montalva fault zone has been mapped crossing the Lajas Valley in southwest Puerto Rico. Identification of the fault was based upon detailed analysis of geophysical data, satellite images, and field mapping. The fault zone consists of a series of Cretaceous bedrock faults that reactivated and deformed Miocene limestone and Quaternary alluvial fan sediments. The fault zone is seismically active (local magnitude greater than 5.0) with numerous locally felt earthquakes. Focal mechanism solutions suggest strain partitioning with predominantly east–west left-lateral displacements with small normal faults striking mostly toward the northeast. Northeast-trending fractures and normal faults can be found in intermittent streams that cut through the Quaternary alluvial fan deposits along the southern margin of the Lajas Valley, an east–west-trending 30-km-long fault-controlled depression. Areas of preferred erosion within the alluvial fan trend toward the west-northwest parallel to the onland projection of the North Boquerón Bay fault. The North Boquerón Bay fault aligns with the Punta Montalva fault southeast of the Lajas Valley. Both faults show strong southward tilting of Miocene strata. On the western end, the Northern Boquerón Bay fault is covered with flat-lying Holocene sediments, whereas at the southern end the Punta Montalva fault shows left-lateral displacement of stream drainage on the order of a few hundred meters.

Puerto Rico↗

USGS2006

Some Key Features of the Strong-Motion Data from the M 6.0 Parkfield, California, Earthquake of 28 September 2004

The 2004 Parkfield, California, earthquake was recorded by an extensive set of strong-motion instruments well positioned to record details of the motion in the near-fault region, where there has previously been very little recorded data. The strong-motion measurements obtained are highly varied, with significant variations occurring over only a few kilometers. The peak accelerations in the near fault region range from 0.13g to over 1.8g (one of the highest acceleration recorded to date, exceeding the capacity of the recording instrument The largest accelerations occurred near the northwest end of the inferred rupture zone. These motions are consistent with directivity for a fault rupturing from the hypocenter near Gold Hill toward the northwest. However, accelerations up to 0.8g were also observed in the opposite direction, at the south end of the Cholame Valley near Highway 41, consistent with bilateral rupture, with rupture southeast of the hypocenter. Several stations near and over the rupturing fault recorded relatively weak motions, consistent with seemingly paradoxical observations of low shaking damage near strike-slip faults. This event had more ground-motion observations within 10 km of the fault than many other earthquakes combined. At moderate distances peak horizontal ground acceleration (PGA) values dropped off more rapidly with distance than standard relationships. At close-in distance the wide variation of PGA suggests a distance-dependent sigma may be important to consider. The near-fault ground-motion variation is greater than that assumed in ShakeMap interpolations, based on the existing set of observed data. Higher density of stations near faults may be the only means in the near future to reduce uncertainty in the interpolations. Outside of the near-fault zone the variance is closer to that assumed. This set of data provides the first case where near-fault radiation has been observed at an adequate number of stations around the fault to allow detailed study of the fault-normal and fault-parallel motion and the near-field S-wave radiation. The fault-normal motions are significant, but they are not large at the central part of the fault, away from the ends. The fault-normal and fault-parallel motions drop off quite rapidly with distance from the fault. Analysis of directivity indicates increased values of peak velocity in the rupture direction. No such dependence is observed in the peak acceleration, except for stations close to the strike of the fault near and beyond the ends of the faulting.

Bulletin of the Seismological Society of America↗

USGS2003

A shallow fault-zone structure illuminated by trapped waves in the Karadere-Duzce branch of the North Anatolian Fault, western Turkey

We discuss the subsurface structure of the Karadere-Duzce branch of the North Anatolian Fault based on analysis of a large seismic data set recorded by a local PASSCAL network in the 6 months following the Mw = 7.4 1999 Izmit earthquake. Seismograms observed at stations located in the immediate vicinity of the rupture zone show motion amplification and long-period oscillations in both P- and S-wave trains that do not exist in nearby off-fault stations. Examination of thousands of waveforms reveals that these characteristics are commonly generated by events that are well outside the fault zone. The anomalous features in fault-zone seismograms produced by events not necessarily in the fault may be referred to generally as fault-zone-related site effects. The oscillatory shear wave trains after the direct S arrival in these seismograms are analysed as trapped waves propagating in a low-velocity fault-zone layer. The time difference between the S arrival and trapped waves group does not grow systematically with increasing source-receiver separation along the fault. These observations imply that the trapping of seismic energy in the Karadere-Duzce rupture zone is generated by a shallow fault-zone layer. Traveltime analysis and synthetic waveform modelling indicate that the depth of the trapping structure is approximately 3-4 km. The synthetic waveform modelling indicates further that the shallow trapping structure has effective waveguide properties consisting of thickness of the order of 100 m, a velocity decrease relative to the surrounding rock of approximately 50 per cent and an S-wave quality factor of 10-15. The results are supported by large 2-D and 3-D parameter space studies and are compatible with recent analyses of trapped waves in a number of other faults and rupture zones. The inferred shallow trapping structure is likely to be a common structural element of fault zones and may correspond to the top part of a flower-type structure. The motion amplification associated with fault-zone-related site effects increases the seismic shaking hazard near fault-zone structures. The effect may be significant since the volume of sources capable of generating motion amplification in shallow trapping structures is large.

Geophysical Journal International↗

USGS2010

Shallow subsurface structure of the Wasatch fault, Provo segment, Utah, from integrated compressional and shear-wave seismic reflection profiles with implications for fault structure and development

Integrated vibroseis compressional and experimental hammer-source, shear-wave, seismic reflection profiles across the Provo segment of the Wasatch fault zone in Utah reveal near-surface and shallow bedrock structures caused by geologically recent deformation. Combining information from the seismic surveys, geologic mapping, terrain analysis, and previous seismic first-arrival modeling provides a well-constrained cross section of the upper ~500 m of the subsurface. Faults are mapped from the surface, through shallow, poorly consolidated deltaic sediments, and cutting through a rigid bedrock surface. The new seismic data are used to test hypotheses on changing fault orientation with depth, the number of subsidiary faults within the fault zone and the width of the fault zone, and the utility of integrating separate elastic methods to provide information on a complex structural zone. Although previous surface mapping has indicated only a few faults, the seismic section shows a wider and more complex deformation zone with both synthetic and antithetic normal faults. Our study demonstrates the usefulness of a combined shallow and deeper penetrating geophysical survey, integrated with detailed geologic mapping to constrain subsurface fault structure. Due to the complexity of the fault zone, accurate seismic velocity information is essential and was obtained from a first-break tomography model. The new constraints on fault geometry can be used to refine estimates of vertical versus lateral tectonic movements and to improve seismic hazard assessment along the Wasatch fault through an urban area. We suggest that earthquake-hazard assessments made without seismic reflection imaging may be biased by the previous mapping of too few faults. ?? 2010 Geological Society of America.

Geological Society of America Bulletin↗

USGS2024

Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?

The Rose Canyon fault is the southern extension of the larger Newport–Inglewood–Rose Canyon fault system, which represents a major structural boundary in the Inner Continental Borderland (ICB) offshore of southern California. Ten to fifteen percent of total plate boundary motion in southern California is thought to be accommodated by the faults of the ICB, but the exact distribution of slip is uncertain. With an onshore segment, the Rose Canyon fault offers an opportunity to measure the slip rate using traditional geodetic methods. In this study, we use Global Positioning System (GPS) surface velocities from a combined campaign and continuous GPS network to constrain elastic models of the Rose Canyon fault. We then compare the observed surface velocities with proposed conceptual models of regional fault connections that facilitate the transfer of slip into the Rose Canyon fault to assess how well the observations are explained by the models. The results of elastic half‐space models suggest that the Rose Canyon fault may be slipping toward the higher end of geologic estimates, with the preferred model indicating a slip rate of 2.4 ± 0.5 mm/yr. Although limited in terms of near‐fault benchmarks, we find an improved model fit using an asymmetrical elastic half‐space model and a higher slip rate, suggesting a potential rheological contrast across the Rose Canyon fault, similar to observations from the northern Newport–Inglewood fault segments. Observed GPS surface velocities, background seismicity, and gravity anomalies south of San Diego Bay point toward a more easterly trace for the Rose Canyon fault, suggesting a possible connection with the San Miguel–Vallecitos fault system. Such a connection could increase the potential rupture lengths of future earthquakes and have important consequences for regional seismic hazards.

California↗

USGS2014

Paleoearthquakes at Frazier Mountain, California delimit extent and frequency of past San Andreas Fault ruptures along 1857 trace

Large earthquakes are infrequent along a single fault, and therefore historic, well-characterized earthquakes exert a strong influence on fault behavior models. This is true of the 1857 Fort Tejon earthquake (estimated M7.7–7.9) on the southern San Andreas Fault (SSAF), but an outstanding question is whether the 330 km long rupture was typical. New paleoseismic data for six to seven ground-rupturing earthquakes on the Big Bend of the SSAF restrict the pattern of possible ruptures on the 1857 stretch of the fault. In conjunction with existing sites, we show that over the last ~650 years, at least 75% of the surface ruptures are shorter than the 1857 earthquake, with estimated rupture lengths of 100 to <300 km. These results suggest that the 1857 rupture was unusual, perhaps leading to the long open interval, and that a return to pre-1857 behavior would increase the rate of M7.3–M7.7 earthquakes.

California↗

USGS1987

Nucleation and triggering of earthquake slip: Effect of periodic stresses

Results of stability analyses for spring and slider systems, with state variable constitutive properties, are applied to slip on embedded fault patches. Unstable slip may nucleate only if the slipping patch exceeds some minimum size. Subsequent to the onset of instability the earthquake slip may propagate well beyond the patch. It is proposed that the seismicity of a volume of the earth's crust is determined by the distribution of initial conditions on the population of fault patches that nucleate earthquake slip, and the loading history acting upon the volume. Patches with constitutive properties inferred from laboratory experiments are characterized by an interval of self-driven accelerating slip prior to instability, if initial stress exceeds a minimum threshold. This delayed instability of the patches provides an explanation for the occurrence of aftershocks and foreshocks including decay of earthquake rates by time −1 . A population of patches subjected to loading with a periodic component results in periodic variation of the rate of occurrence of instabilities. The change of the rate of seismicity for a sinusoidal load is proportional to the amplitude of the periodic stress component and inversely proportional to both the normal stress acting on the fault patches and the constitutive parameter, A 1 , that controls the direct velocity dependence of fault slip. Values of A 1 representative of laboratory experiments indicate that in a homogeneous crust, correlation of earthquake rates with earth tides should not be detectable at normal stresses in excess of about 8 MPa. Correlation of earthquakes with tides at higher normal stresses can be explained if there exist inhomogeneities that locally amplify the magnitude of the tidal stresses. Such amplification might occur near magma chambers or other soft inclusions in the crust and possibly near the ends of creeping fault segments if the creep or afterslip rates vary in response to tides. Observations of seismicity rate variations associated with seasonal fluctuations of reservoir levels appear to be consistent with the model.

Tectonophysics↗

USGS2007

Deep faulting and structural reactivation beneath the southern Illinois basin

The investigation of deep fault structure and seismogenesis within "stable" continental interiors has been hindered by the paucity of detailed subsurface information and by low levels of seismicity. Outstanding seismotectonic questions for these areas include whether pre-existing structures govern the release of seismic energy as earthquakes, can reactivation of such structures be recognized, and to what extent have Precambrian basement structures exerted long-lived controls on the development of overlying Phanerozoic features. The southern portion of the Illinois basin provides a premier area in which to study the relation between contemporary seismicity and pre-existing structures due to the frequency of seismic events, the concentration of available geophysical data, and the wealth of borehole information. We have integrated the study of this information in order to create a 2.5-dimensional picture of the earth for local seismogenic depths (0-15 km) for a study area of moderate 20th century earthquake activity. The area is located along the western flanks of two of the major structures within the Illinois basin, the Wabash Valley fault system (WVFS) and the La Salle anticlinal belt (LSA). The results of reprocessing seismic reflection profiles, combined with earthquake hypocenter parameters, suggest three distinct seismotectonic environments in the upper crust. First, we have delineated a fault pattern that appears to correspond to the steep nodal plane of a strike-slip mechanism event (1974.04.03; mb = 4.7). The fault pattern is interpreted to be a deeply buried rift zone or zone of intense normal faulting underpinning a major Paleozoic depocenter of the Illinois basin (Fairfield basin). Second, a similar event (1987.06.10; mb = 5.2) and its well-located aftershocks define a narrow zone of deformation that occurs along and parallel to the frontal thrust of the LSA. Third, the hypocenter of the largest event in the study area (1968.11.09; mb = 5.5) may be spatially associated with a prominent zone of dipping middle crustal reflections, just west of the WVFS, which have been interpreted as a deeply buried blind thrust. The proposed correlation of pre-existing structures with earthquakes having consistently oriented structural parameters supports the reactivation of old deformation zones by contemporary stresses as previously proposed by earlier workers. However, the degree to which deformation has propagated upward from Precambrian basement into the Paleozoic rocks varied significantly even over a small study area. The societal value of associating an earthquake with a specific pre-existing deformation zone in the seismogenic crust is to improve the assessment of seismic hazard or to assess the integrity of a stratigraphic formation, being considered as a target for natural gas storage or carbon sequestration. ?? 2007 Elsevier B.V. All rights reserved.

Precambrian Research↗

USGS2020

Dynamic rupture simulations of the M6.4 and M7.1 July 2019 Ridgecrest, California earthquakes

The largest earthquakes of the 2019 Ridgecrest, California, sequence were a M 6.4 left‐lateral rupture followed 34 hr later by a M 7.1 on a perpendicular right‐lateral fault. We use dynamic rupture modeling to address the questions of why the first earthquake did not propagate through the right‐lateral fault in one larger event, whether stress changes from the M 6.4 were necessary for the M 7.1 to occur, and how the Ridgecrest earthquakes affected the nearby Garlock Fault. We find that dynamic clamping and shear stress reduction confined surface rupture in the M 6.4 to the left‐lateral fault. We also find that stress changes from the M 6.4 were not necessary to allow a M 7.1 on the right‐lateral fault but that they affected the slip and likely accelerated the timing of the M 7.1. Lastly, we find that the Ridgecrest earthquakes may have brought the central Garlock Fault closer to failure.

California↗

USGS2021

Physics‐based evaluation of the maximum magnitude of potential earthquakes induced by the Hutubi (China) underground gas storage

The world’s largest underground gas storage facility in Hutubi (HUGS), China, is a unique case where cyclic gas injection‐extraction induced both seismicity and ground deformation. To assess the potential for future induced seismicity, we develop a framework physically based on a well‐constrained hydro‐geomechanical model and on fully coupled poroelastic simulations. We first interpret the spatiotemporal distribution and focal mechanisms of induced earthquakes and use these to estimate the magnitude and location of the largest potential earthquake. The sharp increase in seismicity was controlled by poroelastic loading on secondary southwest‐dipping thrust faults with spatial scales too small to be resolved by 3D seismic surveys. Both operational and local geological factors affect the seismic productivity at the HUGS site, distinguishing it from most cases of seismicity induced by wastewater disposal and hydraulic fracturing. We then conduct slip tendency analyses for major faults imaged by the seismic data, including the largest reservoir‐bounding Hutubi fault hydraulically connected to injection wells. The reactivation potentials of these imaged faults are estimated to be extremely low. Accordingly, future seismicity would most likely occur on failure‐prone secondary faults in regions with positive stress perturbation due to poroelastic loading. The maximum magnitude likely depends on the spatial scales of the secondary faults. As the occurrence of detected earthquakes is spatially and temporally consistent with the simulated evolution of Coulomb stress perturbation, the location of the largest potential earthquake probably depends on the sizes of the poroelastic stressing regions.

Hutubi↗

USGS2012

A teleseismic analysis of the New Brunswick earthquake of January 9, 1982

The analysis of the New Brunswick earthquake of January 9, 1982, has important implications for the evaluation of seismic hazards in eastern North America. Although moderate in size ( m b 5.7), it was well-recorded teleseismically. Source characteristics of this earthquake have been determined from analysis of data that were digitally recorded by the Global Digital Seismograph Network. From broadband displacement and velocity records of P waves, we have obtained a dynamic description of the rupture process as well as conventional static properties of the source. The depth of the hypocenter is estimated to be 9 km from depth phases. The focal mechanism determined from the broadband data corresponds to predominantly thrust faulting. From the variation in the waveforms the direction of slip is inferred to be updip on a west dipping NNE striking fault plane. The steep dip of the inferred fault plane suggests that the earthquake occurred on a preexisting fault that was at one time a normal fault. From an inversion of bodywave pulse durations, the estimated rupture length is 5.5 km. Average properties of the rupture process were examined by a moment tensor analysis of long-period P and SH body waves. The long-period moment of this earthquake was 5.3 × 10 24 dyne cm. The static and dynamic stress drops are 41 and 65 bars, respectively, similar to those of many earthquakes with similar moment in regions that are more seismically active. The joint epicenter determination algorithm was used to locate, relative to the mainshock, the three teleseismically recorded aftershocks that occurred through March 31, 1982. The relocated hypocenters of the aftershocks are significantly different from each other and from that of the mainshock; they provide additional support for the source dimensions inferred from the waveform analysis.

Journal of Geophysical Research Solid Earth↗

USGS1999

Analysis of the tsunami generated by the MW 7.8 1906 San Francisco earthquake

We examine possible sources of a small tsunami produced by the 1906 San Francisco earthquake, recorded at a single tide gauge station situated at the opening to San Francisco Bay. Coseismic vertical displacement fields were calculated using elastic dislocation theory for geodetically constrained horizontal slip along a variety of offshore fault geometries. Propagation of the ensuing tsunami was calculated using a shallow-water hydrodynamic model that takes into account the effects of bottom friction. The observed amplitude and negative pulse of the first arrival are shown to be inconsistent with small vertical displacements (∼4–6 cm) arising from pure horizontal slip along a continuous right bend in the San Andreas fault offshore. The primary source region of the tsunami was most likely a recently recognized 3 km right step in the San Andreas fault that is also the probable epicentral region for the 1906 earthquake. Tsunami models that include the 3 km right step with pure horizontal slip match the arrival time of the tsunami, but underestimate the amplitude of the negative first-arrival pulse. Both the amplitude and time of the first arrival are adequately matched by using a rupture geometry similar to that defined for the 1995 M w (moment magnitude) 6.9 Kobe earthquake: i.e., fault segments dipping toward each other within the stepover region (83° dip, intersecting at 10 km depth) and a small component of slip in the dip direction (rake = −172°). Analysis of the tsunami provides confirming evidence that the 1906 San Francisco earthquake initiated at a right step in a right-lateral fault and propagated bilaterally, suggesting a rupture initiation mechanism similar to that for the 1995 Kobe earthquake.

Geology↗

USGS2006

Geophysical setting of the 2000 ML 5.2 Yountville, California, earthquake: Implications for seismic Hazard in Napa Valley, California

The epicenter of the 2000 ML 5.2 Yountville earthquake was located 5 km west of the surface trace of the West Napa fault, as defined by Helley and Herd (1977). On the basis of the re-examination of geologic data and the analysis of potential field data, the earthquake occurred on a strand of the West Napa fault, the main basin-bounding fault along the west side of Napa Valley. Linear aeromagnetic anomalies and a prominent gravity gradient extend the length of the fault to the latitude of Calistoga, suggesting that this fault may be capable of larger-magnitude earthquakes. Gravity data indicate an ???2-km-deep basin centered on the town of Napa, where damage was concentrated during the Yountville earthquake. It most likely played a minor role in enhancing shaking during this event but may lead to enhanced shaking caused by wave trapping during a larger-magnitude earthquake.

Bulletin of the Seismological Society of America↗

USGS2024

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America↗