Search USGS⌕ Search

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.

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

At least 919 records · Page 51Linked to original sources
USGS2013

Spatial variability of "Did You Feel It?" intensity data: insights into sampling biases in historical earthquake intensity distributions

Recent parallel development of improved quantitative methods to analyze intensity distributions for historical earthquakes and of web‐based systems for collecting intensity data for modern earthquakes provides an opportunity to reconsider not only important individual historical earthquakes but also the overall characterization of intensity distributions for historical events. The focus of this study is a comparison between intensity distributions of historical earthquakes with those from modern earthquakes for which intensities have been determined by the U.S. Geological Survey “Did You Feel It?” (DYFI) website (see Data and Resources). As an example of a historical earthquake, I focus initially on the 1843 Marked Tree, Arkansas, event. Its magnitude has been previously estimated as 6.0–6.2. I first reevaluate the macroseismic effects of this earthquake, assigning intensities using a traditional approach, and estimate a preferred magnitude of 5.4. Modified Mercalli intensity (MMI) values for the Marked Tree earthquake are higher, on average, than those from the 2011 >Mw 5.8 Mineral, Virginia, earthquake for distances ≤500 km but comparable or lower on average at larger distances, with a smaller overall felt extent. Intensity distributions for other moderate historical earthquakes reveal similar discrepancies; the discrepancy is also even more pronounced using earlier published intensities for the 1843 earthquake. I discuss several hypotheses to explain the discrepancies, including the possibility that intensity values associated with historical earthquakes are commonly inflated due to reporting/sampling biases. A detailed consideration of the DYFI intensity distribution for the Mineral earthquake illustrates how reporting and sampling biases can account for historical earthquake intensity biases as high as two intensity units and for the qualitative difference in intensity distance decays for modern versus historical events. Thus, intensity maps for historical earthquakes tend to imply more widespread damage patterns than are revealed by intensity distributions of modern earthquakes of comparable magnitude. However, intensity accounts of historical earthquakes often include fragmentary accounts suggesting long‐period shaking effects that will likely not be captured fully in historical intensity distributions.

Bulletin of the Seismological Society of America↗
USGS2002

Relationship of the 1999 Hector Mine and 1992 Landers fault ruptures to offsets on neogene faults and distribution of late Cenozoic basins in the eastern California shear zone

This report examines the Hector Mine and Landers earthquakes in the broader context of faults and fault-related basins of the eastern California shear zone (ECSZ). We compile new estimates of total strike-slip offset (horizontal separation) at nearly 30 fault sites based on offset magnetic anomaly pairs. We also present a map of the depth to pre-Cenozoic basement rock (thickness of basin-filling late Cenozoic deposits) for the region, based on an inversion of gravity and geologic data. Our estimates of total long-term strike-slip offsets on faults that slipped during the 1999 Hector Mine (3.4 km), and the 1992 Landers earthquakes (3.1 ? to 4.6 km) fall within the 3- to 5-km range of total strike-slip offset proposed for most faults of the western ECSZ. Faults having offsets as great as 20 km are present in the eastern part of the ECSZ. Although the Landers rupture followed sections of a number of faults that had been mapped as independent structures, the similarity in total strike-slip offset associated with these faults is compatible with one of the following hypotheses: (1) the Landers multistrand rupture is a typical event for this linked fault system or (2) this complex rupture path has acted as a coherent entity when viewed over some characteristic multiearthquake cycle. The second hypothesis implies that, for each cycle, slip associated with smaller earthquakes on individual fault segments integrates to a uniform slip over the length of the linked faults. With one exception, the region surrounding the Hector Mine and Landers ruptures is devoid of deep late Cenozoic basins. In particular, no deep basins are found immediately north of the Pinto Mountain fault, a place where a number of kinematic models for development of the ECSZ have predicted basins. In contrast, some basins exist near Barstow and along the eastern part of the ECSZ, where the model of Dokka et al. (1998) predicts basins.

Bulletin of the Seismological Society of America↗
USGS1998

Faulting apparently related to the 1994 Northridge, California, earthquake and possible co-seismic origin of surface cracks in Potrero Canyon, Los Angeles County, California

Apparent southward-dipping, reverse-fault zones are imaged to depths of about 1.5 km beneath Potrero Canyon, Los Angeles County, California. Based on their orientation and projection to the surface, we suggest that the imaged fault zones are extensions of the Oak Ridge fault. Geologic mapping by others and correlations with seismicity studies suggest that the Oak Ridge fault is the causative fault of the 17 January 1994 Northridge earthquake (Northridge fault). Our seismically imaged faults may be among several faults that collectively comprise the Northridge thrust fault system. Unusually strong shaking in Potrero Canyon during the Northridge earthquake may have resulted from focusing of seismic energy or co-seismic movement along existing, related shallow-depth faults. The strong shaking produced ground-surface cracks and sand blows distributed along the length of the canyon. Seismic reflection and refraction images show that shallow-depth faults may underlie some of the observed surface cracks. The relationship between observed surface cracks and imaged faults indicates that some of the surface cracks may have developed from nontectonic alluvial movement, but others may be fault related. Immediately beneath the surface cracks, P-wave velocities are unusually low (<400 m/sec), and there are velocity anomalies consistent with a seismic reflection image of shallow faulting to depths of at least 100 m. On the basis of velocity data, we suggest that unconsolidated soils (<800 m/sec) extend to depths of about 15 to 20 m beneath our datum (<25 m below ground surface). The underlying rocks range in velocity from about 1000 to 5000 m/sec in the upper 100 m. This study illustrates the utility of high-resolution seismic imaging in assessing local and regional seismic hazards.

California↗
USGS2010

The influence of maximum magnitude on seismic-hazard estimates in the Central and Eastern United States

I analyze the sensitivity of seismic-hazard estimates in the central and eastern United States (CEUS) to maximum magnitude (m max ) by exercising the U.S. Geological Survey (USGS) probabilistic hazard model with several m max alternatives. Seismicity-based sources control the hazard in most of the CEUS, but data seldom provide an objective basis for estimating m max . The USGS uses preferred m max values of moment magnitude 7.0 and 7.5 for the CEUS craton and extended margin, respectively, derived from data in stable continental regions worldwide. Other approaches, for example analysis of local seismicity or judgment about a source's seismogenic potential, often lead to much smaller m max . Alternative models span the m max ranges from the 1980s Electric Power Research Institute/Seismicity Owners Group (EPRI/SOG) analysis. Results are presented as haz-ard ratios relative to the USGS national seismic hazard maps. One alternative model specifies m max equal to moment magnitude 5.0 and 5.5 for the craton and margin, respectively, similar to EPRI/SOG for some sources. For 2% probability of exceedance in 50 years (about 0.0004 annual probability), the strong m max truncation produces hazard ratios equal to 0.35-0.60 for 0.2-sec spectral acceleration, and 0.15-0.35 for 1.0-sec spectral acceleration. Hazard-controlling earthquakes interact with m max in complex ways. There is a relatively weak dependence on probability level: hazardratios increase 0-15% for 0.002 annual exceedance probability and decrease 5-25% for 0.00001 annual exceedance probability. Although differences at some sites are tempered when faults are added, m max clearly accounts for some of the discrepancies that are seen in comparisons between USGS-based and EPRI/SOG-based hazard results.

Bulletin of the Seismological Society of America↗
USGS1985

A seismic-refraction profile across the San Andreas, Sargent, and Calaveras faults, west-central California

In 1981, the United States Geological Survey recorded a seismic-refraction profile across the southern Santa Cruz Mountains in west-central California to examine the shallow velocity structure of this seismogenic region. This 40-km-long profile, which consisted of three shotpoints, extended northeastward from near Watsonville, California, to Coyote Lake, crossing the San Andreas, Sargent, and Calaveras faults. This entire region is characterized by a highly heterogeneous upper crust. West of Watsonville, 1 km of alluvium with a velocity of 2.12 km/sec overlies a basement with a velocity 5.45 km/sec. The abrupt deepening of basement by 1.5 km just east of Watsonville at a subsurface fault suggests that the Zayante fault to the north and the Vergeles fault to the south are connected. The Tertiary sediment at the San Andreas fault is 2.5 km thick and has a velocity of 3.34 km/sec. The San Andreas fault is not marked by any distinctive seismic velocity features, whereas a 1-km-wide low velocity zone is found at both the Sargent and Calaveras faults. East of the Sargent fault, the basement surface forms a broad anticlinal structure, with velocities ranging from 5.4 to 6.0 km/sec. From the anticlinal crest, basement dips to the east beneath the Santa Clara Valley and reaches a maximum depth of 1 km on the east side; the overlying alluvium has a velocity of 2.7 km/sec. At the crest of the basement anticlinal structure, a vertical low-velocity zone coincident with exposed serpentine provides strong evidence of faulting.

California↗
USGS2026

Detecting earthquakes in noisy real-time GNSS data with deep learning for improved PGD magnitude estimation

To disseminate accurate and useful warnings, earthquake early warning (EEW) systems must quickly determine the size and location of an earthquake to estimate expected shaking. Traditional seismic‐based algorithms tend to underestimate the true magnitudes of large earthquakes, a phenomenon known as magnitude saturation. This limitation motivated the recent inclusion of Global Navigation Satellite Systems (GNSS) data into the U.S. Geological Survey’s ShakeAlert EEW system with the Geodetic First Approximation of Size and Time (GFAST) algorithm because GNSS data do not saturate with large ground motions. However, the noise levels of GNSS data are very high compared with traditional seismic data, which obscures P ‐wave arrivals and can result in less accurate magnitude estimations if displacement amplitudes are low, such as for lower magnitude earthquakes or large source–station distances. In this study, we develop a deep‐learning model that detects earthquakes in GNSS data and use the Ridgecrest, California, earthquake sequence as a case study to demonstrate how the model could act as a filter to reduce the amount of low‐quality data that enters an algorithm like GFAST. To preserve our limited real earthquake data for model inference, we generated a training dataset composed of >700,000 synthetic displacement waveforms. We combined the synthetic waveforms with real‐time GNSS noise to produce realistically noisy training waveforms and then tested our model on additional synthetic data and performed inference using the real data that were held back. We discuss the performance of our trained model on both the unseen synthetic data and real inference data. Our model can be used to selectively filter only high‐quality data where an earthquake signal is observed for input into an algorithm like GFAST (outperforming a simple signal‐to‐noise ratio–based filter) to reduce the error in GFAST’s real‐time earthquake magnitude estimations.

California↗
USGS2020

Basin amplification effects in the Puget Lowland, Washington from strong motion recordings and 3D simulations

Sedimentary basins in the Puget Sound region, Washington State, increase ground‐motion intensity and duration of shaking during local earthquakes. We analyze Pacific Northwest Seismic Network and U.S. Geological Survey strong‐motion recordings of five local earthquakes ( M 3.9–6.8), including the 2001 Nisqually earthquake, to characterize sedimentary basin effects within the Seattle and Tacoma basins. We observe basin‐edge generated surface waves at sites within the Seattle basin for most ray paths that cross the Seattle fault zone. We also note previously undocumented basin‐edge surface waves in the Tacoma basin during one of the local earthquakes. To place quantitative constraints on basin amplification, we determine amplification factors by computing the spectral ratios of inside‐basin sites to outside‐basin sites at 1, 2, 3, and 5 s periods. Ground shaking is amplified in the Seattle basin for all the earthquakes analyzed and for a subset of events in the Tacoma basin. We find that the largest amplification factors in the Seattle basin are produced by a shallow earthquake located to the southwest of the basin. Our observation suggests that future shallow crustal and megathrust earthquakes rupturing west of the Puget Lowland will produce greater amplification within the Seattle basin than has been seen for intraslab events. We also perform ground‐motion simulations using a finite‐difference method to validate a 3D Cascadia velocity model (CVM) by comparing properties of observed and synthetic waveforms up to a frequency of 1 Hz. Basin‐edge effects are well reproduced in the Seattle basin, but are less well resolved in the Tacoma basin. Continued study of basin effects in the Tacoma basin would improve the CVM.

Washington↗
USGS2014

A method and example of seismically imaging near‐surface fault zones in geologically complex areas using Vp, Vs, and their ratios

The determination of near‐surface (vadose zone and slightly below) fault locations and geometries is important because assessment of ground rupture, strong shaking, geologic slip rates, and rupture histories occurs at shallow depths. However, seismic imaging of fault zones at shallow depths can be difficult due to near‐surface complexities, such as weathering, groundwater saturation, massive (nonlayered) rocks, and vertically layered strata. Combined P ‐ and S ‐wave seismic‐refraction tomography data can overcome many of the near‐surface, fault‐zone seismic‐imaging problems because of differences in the responses of elastic (bulk and shear) moduli of P and S waves to shallow‐depth, fault‐zone properties. We show that high‐resolution refraction tomography images of P ‐ to S ‐wave velocity ratios ( V P / V S ) can reliably identify near‐surface faults. We demonstrate this method using tomography images of the San Andreas fault (SAF) surface‐rupture zone associated with the 18 April 1906 ∼ M 7.9 San Francisco earthquake on the San Francisco peninsula in California. There, the SAF cuts through Franciscan mélange, which consists of an incoherent assemblage of greywacke, chert, greenstone, and serpentinite. A near‐vertical zone (∼75° northeast dip) of high P ‐wave velocities (up to 3000 m/s), low S ‐wave velocities (∼150–600 m/s), high V P / V S ratios (4–8.8), and high Poisson’s ratios (0.44–0.49) characterizes the main surface‐rupture zone to a depth of about 20 m and is consistent with nearby trench observations. We suggest that the combined V P / V S imaging approach can reliably identify most near‐surface fault zones in locations where many other seismic methods cannot be applied.

Bulletin of the Seismological Society of America↗
USGS2000

Accounting for site effects in probabilistic seismic hazard analyses of southern California: Overview of the SCEC Phase III Report

This article presents an overview of the Southern California Earthquake Center (SCEC) Phase-III effort to determine the extent to which probabilistic seismic hazard analysis (PSHA) can be improved by accounting for site effects. The contributions made in this endeavor are represented in the various articles that compose this special issue of BSSA. Given the somewhat arbitrary nature of the site-effect distinction, it must be carefully defined in any given context. With respect to PSHA, we define the site effect as the response, relative to an attenuation relationship, averaged over all damaging earthquakes in the region. A diligent effort has been made to identify any attributes that predispose a site to greater or lower levels of shaking. The most detailed maps of Quaternary geology are not found to be helpful; either they are overly detailed in terms of distinguishing different amplification factors or present southern California strong-motion observations are inadequate to reveal their superiority. A map based on the average shear-wave velocity in the upper 30 m, however, is found to delineate significantly different amplification factors. A correlation of amplification with basin depth is also found to be significant, implying up to a factor of two difference between the shallowest and deepest parts of the Los Angeles basin. In fact, for peak acceleration the basin-depth correction is more influential than the 30-m shear-wave velocity. Questions remain, however, as to whether basin depth is a proxy for some other site attribute. In spite of these significant and important site effects, the standard deviation of an attenuation relationship (the prediction error) is not significantly reduced by making such corrections. That is, given the influence of basin-edge-induced waves, subsurface focusing, and scattering in general, any model that attempts to predict ground motion with only a few parameters will have a substantial intrinsic variability. Our best hope for reducing such uncertainties is via waveform modeling based on first principals of physics. Finally, questions remain with respect to the overall reliability of attenuation relationships at large magnitudes and short distances. Current discrepancies between viable models produce up to a factor of 3 difference among predicted 10% in 50-yr exceedance levels, part of which results from the uncertain influence of sediment nonlinearity.

California↗
USGS2025

Seismic moment and local magnitude scales in Ridgecrest, CA from the SCEC/USGS Community Stress Drop Validation Study

We illustrate the systematic difference between moment magnitude and local magnitude caused by underlying earthquake source physics, using seismic moments submitted to the Statewide California Earthquake Center/United States Geological Survey Community Stress Drop Validation Study 2019 Ridgecrest data set. While the relationship between seismic moment and moment magnitude ( M or M w ) of log 10 ( M 0 ) ~ 1.5* M is uniformly valid for all earthquake sizes by definition (Hanks and Kanamori, 1979), the relationship between local magnitude M L and moment is itself magnitude dependent. For moderate events, ~3< M < ~6, M and M L are coincident; for earthquakes smaller than ~3, M L ~ 1.0 log 10 M 0 (Hanks and Boore, 1984). This is a physical consequence of the corner frequency fc becoming larger than the upper frequency of observation and implies that M L and M differ systematically by a factor of 1.5 for these small events. While this idea is not new, we propose a new, continuous relationship between local magnitude and moment, for magnitudes 2 to 6 which extrapolates to smaller and larger magnitudes, applicable to southern California specific to the Ridgecrest region. We make use of the plethora of seismic moments as submitted by many participants of the Community Stress Drop study, compared to the Southern California Seismic Network (SCSN) catalog magnitudes. Overall, the seismic moments in the Community Study recover moment magnitude well, so we use our new M L - M 0 to convert M L to M , refining the SCSN operational M Lr scale. This systematic difference of 50% in slope between local and moment magnitude at small magnitudes has implications for spectral stress drop estimates, earthquake ground motion modeling, as well as other magnitude scales and earthquake occurrence statistics.

California↗
USGS2023

The USGS 2023 Conterminous U.S. time‐independent earthquake rupture forecast

We present the 2023 U.S. Geological Survey time‐independent earthquake rupture forecast for the conterminous United States, which gives authoritative estimates of the magnitude, location, and time‐averaged frequency of potentially damaging earthquakes throughout the region. In addition to updating virtually all model components, a major focus has been to provide a better representation of epistemic uncertainties. For example, we have improved the representation of multifault ruptures, both in terms of allowing more and less fault connectivity than in the previous models, and in sweeping over a broader range of viable models. An unprecedented level of diagnostic information has been provided for assessing the model, and the development was overseen by a 19‐member participatory review panel. Although we believe the new model embodies significant improvements and represents the best available science, we also discuss potential model limitations, including the applicability of logic tree branch weights with respect different types of hazard and risk metrics. Future improvements are also discussed, with deformation model enhancements being particularly worthy of pursuit, as well as better representation of sampling errors in the gridded seismicity components. We also plan to add time‐dependent components, and assess implications with a wider range of hazard and risk metrics.

Bulletin of the Seismological Society of America↗
USGS2019

Shallow structure and geomorphology along the offshore northern San Andreas Fault, Tomales Point to Fort Ross, California

We mapped a poorly documented 35-km-long section of the northern San Andreas fault zone (NSAF) between Tomales Point and Fort Ross, California. Mapping is largely based on high-resolution seismic-reflection profiles (38 fault crossings), multibeam bathymetry, and onshore geology. NSAF strike in this section is nearly parallel to plate motion, characterized by a slight (~2°) northerly (transtensional) bend in the south between Tomales Bay and the Bodega isthmus, and a northwesterly (transpressional) ~5° bend in the north between the Bodega isthmus and Fort Ross. The southern transtensional bend is the northern part of the now-submerged, linear, ~50-km-long and 1- to 2-km-wide, "Tomales-Bodega valley." The valley floor is cut by a complex zone of subparallel, variably continuous fault strands and the deformed valley fill is an inferred mix of late Quaternary marine and nonmarine strata. In the northern part of this elongate valley, Holocene fault offset occurred on two fault strands about 740 m apart. The northern transpressional bend is characterized by narrow, elongate, asymmetric basins containing as much as 56 m of inferred latest Pleistocene to Holocene sediment. Between Bodega Head and Fort Ross, the gently dipping (~0.8°) shelf includes two large (4.8 and 5.9 km2) zones of sediment failure that we speculatively correlate with the 1906 San Francisco NSAF earthquake. Similar sediment-failure zones should be common along offshore reaches of the NSAF and other nearshore fault zones, but have apparent limited preservation potential. Onland geomorphic impacts of the mainly offshore NSAF include: (1) Northward upwarping of uplifted marine terraces in the transpressional zone north of Bodega Bay; and (2) Blocking of littoral sediment transport by uplifts on the west flank of the NSAF at Bodega Head and Tomales Point, resulting in rapidly accreting beaches and large coastal sand dune complexes.

California↗
USGS2025

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California↗
USGS2004

Interpretation of the Seattle Uplift, Washington, as a passive-roof duplex

We interpret seismic lines and a wide variety of other geological and geophysical data to suggest that the Seattle uplift is a passive-roof duplex. A passive-roof duplex is bounded top and bottom by thrust faults with opposite senses of vergence that form a triangle zone at the leading edge of the advancing thrust sheet. In passive-roof duplexes the roof thrust slips only when the floor thrust ruptures. The Seattle fault is a south-dipping reverse fault forming the leading edge of the Seattle uplift, a 40-km-wide fold-and-thrust belt. The recently discovered, north-dipping Tacoma reverse fault is interpreted as a back thrust on the trailing edge of the belt, making the belt doubly vergent. Floor thrusts in the Seattle and Tacoma fault zones, imaged as discontinuous reflections, are interpreted as blind faults that flatten updip into bedding plane thrusts. Shallow monoclines in both the Seattle and Tacoma basins are interpreted to overlie the leading edges of thrust-bounded wedge tips advancing into the basins. Across the Seattle uplift, seismic lines image several shallow, short-wavelength folds exhibiting Quaternary or late Quaternary growth. From reflector truncation, several north-dipping thrust faults (splay thrusts) are inferred to core these shallow folds and to splay upward from a shallow roof thrust. Some of these shallow splay thrusts ruptured to the surface in the late Holocene. Ages from offset soils in trenches across the fault scarps and from abruptly raised shorelines indicate that the splay, roof, and floor thrusts of the Seattle and Tacoma faults ruptured about 1100 years ago.

Washington↗
USGS2008

Observation and prediction of dynamic ground strains, tilts, and torsions caused by the Mw 6.0 2004 Parkfield, California, earthquake and aftershocks, derived from UPSAR array observations

The 28 September 2004 Parkfield, California, earthquake (Mw 6.0) and four aftershocks (Mw 4.7-5.1) were recorded on 12 accelerograph stations of the U.S. Geological Survey Parkfield seismic array (UPSAR), an array of three-component accelerographs occupying an area of about 1 km2 located 8.8 km from the San Andreas fault. Peak horizontal acceleration and velocity at UPSAR during the mainshock were 0.45g and 27 cm/sec, respectively. We determined both time-varying and peak values of ground dilatations, shear strains, torsions, tilts, torsion rates, and tilt rates by applying a time-dependent geodetic analysis to the observed array displacement time series. Array-derived dilatations agree fairly well with point measurements made on high sample rate recordings of the Parkfield-area dilatometers (Johnston et al., 2006). Torsion Fourier amplitude spectra agree well with ground velocity spectra, as expected for propagating plane waves. A simple predictive relation, using the predicted peak velocity from the Boore-Atkinson ground-motion prediction relation (Boore and Atkinson, 2007) scaled by a phase velocity of 1 km/sec, predicts observed peak Parkfield and Chi-Chi rotations (Huang, 2003) well. However, rotation rates measured during Mw 5 Ito, Japan, events observed on a gyro sensor (Takeo, 1998) are factors of 5-60 greater than those predicted by our predictive relation. This discrepancy might be caused by a scale dependence in rotation, with rotations measured over a short baseline exceeding those measured over long baselines. An alternative hypothesis is that events having significant non-double-couple mechanisms, like the Ito events, radiate much stronger rotations than double-couple events. If this is true, then rotational observations might provide an important source of new information for monitoring seismicity in volcanic areas.

Bulletin of the Seismological Society of America↗
USGS1988

Teleseismic and near-field analysis of the Nahanni earthquakes in the Northwest Territories, Canada

The analysis of the Nahanni earthquakes of October 5, 1985 (MS 6.6), and December 23, 1985 (MS 6.9), will have important implications for the assessment of seismic hazards in intraplate environments. To maximize the information available to seismic engineers, broadband data recorded teleseismically are analyzed jointly with strong-motion data recorded in the near field. The time-domain analysis of teleseismic data yields the source mechanisms, depths, and complexities of rupture of each earthquake. Both earthquakes occurred as shallow thrusts with centroid depths (6 to 7 km) and shallowly dipping fault planes that correspond well with the aftershock distributions obtained from a local survey run by the Canadian Geological Survey. The shallow nodal plane for the October 5 earthquake dips 30° to the WSW, while the shallow nodal planes of the subevents for the December 23 earthquake dip an average of 23° to the WSW. The October 5 earthquake has an impulsive initial rupture, followed by a weak subevent of longer duration but smaller moment release. The December 23 earthquake exhibits more complexity, being comprised of three subevents of similar size. The subevent delays derived from the teleseismic analysis are used to help interpret arrivals in records of ground velocity recorded in the near field of the December 23 earthquake. The rupture geometries inferred from the joint near- and far-field analysis suggest that the rupture processes were unusually complicated and that the 2g peak that occurs late in one of the near-field records could be a localized phenomenon. Spectral analyses of the teleseismic P waves yield the following source parameters for the October 5 and December 23 earthquakes, respectively: the seismic moments are 1.2 and 1.8 × 1026 dyne-cm, the radiated energies are 1.8 and 2.8 × 1021 dyne-cm, and the dynamic stress drops are 65 and 50 bar. The acceleration source spectra of both earthquakes exhibit an intermediate slope (| üα(ω) | ∝ω) from 0.03 to 0.3 Hz, suggesting that the earthquakes represent the failure of asperities. Extrapolating the teleseismic P-wave spectra to estimate the near-field S-wave spectra yields good fits to the acceleration spectra from two strong motion records, but underestimates the spectra from a third strong motion record with the strongest, but possibly localized, accelerations.

Bulletin of the Seismological Society of America↗
USGS2015

Earthquake geology of the Bulnay Fault (Mongolia)

The Bulnay earthquake of July 23, 1905 (M w 8.3-8.5), in north-central Mongolia, is one of the world's largest recorded intracontinental earthquakes and one of four great earthquakes that occurred in the region during the 20th century. The 375-km-long surface rupture of the left-lateral, strike-slip, N095°E trending Bulnay Fault associated with this earthquake is remarkable for its pronounced expression across the landscape and for the size of features produced by previous earthquakes. Our field observations suggest that in many areas the width and geometry of the rupture zone is the result of repeated earthquakes; however, in those areas where it is possible to determine that the geomorphic features are the result of the 1905 surface rupture alone, the size of the features produced by this single earthquake are singular in comparison to most other historical strike-slip surface ruptures worldwide. Along the 80 km stretch, between 97.18°E and 98.33°E, the fault zone is characterized by several meters width and the mean left-lateral 1905 offset is 8.9 ± 0.6 m with two measured cumulative offsets that are twice the 1905 slip. These observations suggest that the displacement produced during the penultimate event was similar to the 1905 slip. Morphotectonic analyses carried out at three sites along the eastern part of the Bulnay fault, allow us to estimate a mean horizontal slip rate of 3.1 ± 1.7 mm/yr over the Late Pleistocene-Holocene period. In parallel, paleoseismological investigations show evidence for two earthquakes prior to the 1905 event with recurrence intervals of ~2700-4000 years.

Bulnay Fault↗
USGS2023

3-D wave propagation simulations of Mw 6.5+ earthquakes on the Tacoma Fault, Washington state, considering the effects of topography, a geotechnical gradient, and a fault damage zone

We simulate shaking in Tacoma, Washington, and surrounding areas from M w 6.5 and 7.0 earthquakes on the Tacoma fault. Ground motions are directly modeled up to 2.5 Hz using kinematic, finite‐fault sources; a 3D seismic velocity model considering regional geology; and a model mesh with 30 m sampling at the ground surface. In addition, we explore how adjustments to the seismic velocity model affect predicted shaking over a range of periods. These adjustments include the addition of a region‐specific geotechnical gradient, surface topography, and a fault damage zone. We find that the simulated shaking tends to be near estimates from empirical ground‐motion models (GMMs). However, long‐period ( T = 5.0 s) shaking within the Tacoma basin is typically underpredicted by the GMMs. The fit between simulated and GMM‐derived short‐period ( T = 0.5 s) shaking is significantly improved with the addition of the geotechnical gradient. From comparing different M w 6.5 earthquake scenarios, we also find that the response of the Tacoma basin is sensitive to the azimuth of incoming seismic waves. In adding surface topography to the simulation, we find that average ground motion is similar to that produced from the nontopography model. However, shaking is often amplified at topographic highs and deamplified at topographic lows, and the wavefield undergoes extensive scattering. Adding a fault damage zone has the effect of amplifying short‐period shaking adjacent to the fault, while reducing far‐field shaking. Intermediate‐period shaking is amplified within the Tacoma basin, likely due to enhanced surface‐wave generation attributable to the fault damage zone waveguide. When applied in the same model, the topography and fault damage zone adjustments often enhance or reduce the effects of one another, adding further complexity to the wavefield. These results emphasize the importance of improving near‐surface velocity model resolution as waveform simulations progress toward higher frequencies.

Washington↗