Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geochemical Perspectives Letters”

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.

185 records · Page 5Linked to original sources
USGS2026

Landslide volume estimation from seismic waveform features using a global catalog of seismogenic events

Seismic waves generated by large, rapid landslides encode information about the source and can be analyzed rapidly following an event. Even remote landslides can pose hazards to downstream communities, so rapid detection and characterization using existing seismic monitoring networks could be beneficial. In this study, we expand on past regionally limited work by presenting a globally applicable method for estimating landslide volume from seismic features that could be integrated into future landslide seismic monitoring frameworks. We train the model using multivariable linear regression and five seismic features derived from recordings of 129 landslide events of a range of styles and locations with independently estimated volumes. We present two preferred models, one that combines long‐period (LP) and high‐frequency (HF) features and one for use on smaller landslides without observed long‐period signals. We find that our best‐performing model, applicable to landslides larger than 100,000 m 3 with signals containing observable long‐period energy, requires only two features: LP (20–100 s) absolute maximum amplitude and HF (1–5 Hz) rise time (time between the signal start and the maximum envelope amplitude) and has an R 2 score of 0.79. This model predicted volumes within one order of magnitude for 55 out of 58 events. We find that this combination mitigates the trade‐off between mass and acceleration and other variations of landslide style that limit methods based on amplitude alone.

Seismological Research Letters↗
USGS2025

Long‐period ground motions from dynamic rupture simulations of large earthquakes on the creeping Hayward–Calaveras–Rodgers Creek fault system

he Hayward, Calaveras, and Rodgers Creek faults in the San Francisco Bay region of California have a high probability of producing a large earthquake in the next decades. Although these faults creep, the creep is insufficient to keep up with their relatively rapid slip rates on their deepest sections, so they have been storing tectonic strain since their last large earthquakes, with the Hayward’s and Rodgers Creek’s more than 150 yr ago. We do not know what the next large Hayward–Calaveras–Rodgers Creek earthquakes will look like or how strongly they will shake the San Francisco Bay region. Harris et al. (2021) used the 3D dynamic (spontaneous) rupture method to simulate large earthquakes on these creeping faults. In this article, we examine the resulting simulated long‐period ( T > 1 s) ground shaking from 0 to 50 km distance, for earthquakes nucleating on the Hayward fault and earthquakes nucleating on the Rodgers Creek fault. We compare these simulated long‐period ground motions with the Boore et al. (2014) well‐established empirically based ground‐motion model suitable for the slowest material velocity in our 3D velocity structure. We find that the simulated long‐period ground motions from the creeping‐fault earthquake scenarios produce a reasonable agreement with the empirical expectations if frictional cohesion is included only where it is appropriate.

California↗
USGS2026

Preface to the focus section on intraplate earthquakes

More than a half century after plate tectonics provided an overarching framework to explain earthquakes along active plate boundaries, numerous theories have been proposed to explain where, why, and how often earthquakes occur well away from active plate boundaries, but a paradigm remains elusive. Even the classification of earthquakes away from active plate boundaries as "intraplate" raises issues, with potentially important distinctions between Stable Continental Regions and more actively deforming regions including passive margins and failed rifts. Some of the largest known intraplate earthquakes themselves remain enigmatic, having occurred before the modern instrumental era. Hazard assessments are often data-limited: low fault-slip rates relative to landscape modification rates result in poor discoverability of fault sources, challenging the characterization of source zones and earthquake recurrence; the completeness and homogenization of instrumental earthquake catalogs using uncertain magnitude conversions can lead to uncertainties in earthquake recurrence; and, limited strong-motion observations for large-magnitude events at near-source distances leads to uncertainties in the selection and development of ground-motion models for seismic hazard studies. Data from recent intraplate earthquakes around the world—from the moment magnitude M 7.7 2001 Bhuj, India, earthquake 25 years ago to the 2024 M 4.8 Tewksbury, New Jersey earthquake—have yielded both new insights and new questions. The papers in this special focus discuss many of the long-standing challenges involved with intraplate earthquake investigations and provide a snapshot of the state of the art with current research to advance our understanding.

Seismological Research Letters↗
USGS2025

What 25+ years of "Did You Feel It" intensities tell us about shaking in California

“When will the Big One happen?” is a question that people often have for earthquake scientists. But while waiting for the “Big One” to occur, people will usually experience frightening or damaging shaking from multiple relatively smaller‐magnitude earthquakes. Given this context, it raises the question: “Where does most of the damage come from?” Could smaller, yet more frequent, earthquakes account for the majority of reported impactful shaking? To explore this question, we consider reports of earthquake damage and felt shaking experiences from a catalog of community‐collected intensity values from the U.S. Geological Survey’s “Did You Feel It?” system. Comparing these intensities to expectations from a ground‐motion model, we find that earthquakes of magnitudes smaller than expected are responsible for most reported intensities of community decimal intensities (CDI) 4.5 and above (moderate and higher shaking intensity levels). (Here “expected value” is meant in its mathematical sense of the mean or equivalently median shaking intensity.) We also present a regional analysis of observed earthquake shaking for specific areas, detailing the maximum intensity experienced within a predetermined area. We identify several instances of M < 4.5 events that generated maximum intensities of CDI > 5 in regions in California surrounding the cities of Eureka, San Francisco, Los Angeles, and San Diego. Our results motivate the need to include smaller‐magnitude earthquakes in communications about earthquake hazard and risk reduction.

California↗
USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗
USGS2026

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗
USGS2026

Corundum discovered by SuperCam and the Perseverance rover at Jezero crater, Mars

Mars is primarily composed of mafic mineral assemblages and their alteration products, but small, scattered rocks strewn across the landscape offer clues to greater petrological diversity. While traversing the Jezero crater rim, the Perseverance rover encountered several plagioclase-rich light-toned float rocks. SuperCam identified the distinctive signature of corundum (α-Al 2 O 3 ) in these rocks using time-resolved luminescence spectroscopy. Two strong peaks (692.7 and 694.1 nm) with millisecond lifetimes, and additional supporting lines, are consistent with Cr 3+ substitution for Al 3+ in corundum. Corundum forms in Al-rich, Si-depleted environments through magmatic or metamorphic processes. Given the rocks' small size, association with plagioclase, and location on the crater rim, we interpret the most plausible formation scenario to be impact induced metamorphism at the interface of a felsic and a mafic/ultramafic member with the likely action of fluids at some stage, although other possibilities are not excluded.

Geophysical Research Letters↗
USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗
USGS2025

Combining scanning electron microscopy, X-ray diffraction, and X-ray fluorescence to characterize shear zones at the Pogo gold deposit, Alaska

This study employs a multi-method analytical approach to characterize the mineralogical, geochemical, and textural properties of fault rocks from the Pogo gold mine in the Yukon-Tanana Upland, central Alaska. Specifically, we examine cataclasites, to document the structural and geochemical evolution of shear zones and their associations with gold mineralization. To investigate the shear zone, we integrate portable X-ray fluorescence (pXRF), scanning electron microscopy-based automated mineralogy (SEM-AM), X-ray diffraction (XRD), and high-resolution micro-X-ray fluorescence (micro-XRF) mapping. These methods collectively provide insights into bulk and trace element chemistry, mineralogical composition, and deformation-related textures across multiple scales. Handheld pXRF enables rapid geochemical screening, guiding SEM-AM and XRD analyses to ensure consistent mineralogical interpretation. X-ray diffraction identifies and quantifies crystalline phases, while SEM-AM produces high-resolution mineral maps, revealing mineral abundances, grain-scale textures, and gold associations. Micro-XRF mapping further refines our understanding by showing visual trace element distributions at sub-millimetre resolution. By integrating these techniques, we improve our understanding of the nature and geochemistry of Pogo shear zones, their role in gold mineralization, and support metallurgical processing strategies. This approach enhances exploration models and resource characterization for structurally complex gold deposits.

Alaska↗
USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗
USGS2025

Evidence for offset of Cretaceous plutons by the Tintina fault in eastern Alaska: Implications for regional metallogeny

Cretaceous magmatism in eastern interior Alaska is voluminous, but temporally and spatially diverse – suggestive of varying sources and drivers. More than 150 new U-Pb zircon and more than 500 geochemical analyses of Cretaceous plutonic units allow for the grouping of distinct plutonic suites. Magmatism was continuous from 120-66 Ma but can be grouped into temporally distinct pulses from ca. 115-100 Ma, 100-90 Ma, and 75-66 Ma. Geochemical diversity occurs during each pulse, further distinguishing multiple suites. Diverse metallogenic epochs are strongly correlated to pluton chemistry. Mineralization is largely absent prior to 108 Ma. From 108-100 Ma, plutonism is coeval with sparse, but notable Au-quartz veins with variable Bi, As, W, and Mo. From 100-90 Ma, intrusion-related mineralization zones from Au-Cu(-Bi) and U-Th in the northwest to central Au- Bi-As-Te(-W), and Mo-W to the southeast. Porphyry style Cu-Mo(-Au) occurrences occur with the latest Cretaceous plutons emplaced from 75-66 Ma . Restoration of ~450 km of dextral movement on the Tintina fault and comparison of metallogenic and geochemical characteristics of Alaska plutons suggest 100-90 Ma plutons may be the continuation of the metallogenically significant Tombstone, Mayo, and Tungsten suites from the Yukon.

Alaska↗
USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2025

Searching for seismic precursors - The Barry Landslide hazard clean up

The Barry Landslide, located in Barry Arm of Prince William Sound, Alaska, poses a major hazard due to its steep, unstable slopes and the potential for a massive landslide-generated tsunami. With an estimated volume of 500–700 million cubic meters, the Barry Landslide could trigger highly destructive waves. In this study, we focus on seismic signals from the Barry Landslide, which are critical for providing timely tsunami warnings. Since the summer of 2020, the region has been instrumented to monitor the landslide, but the seismic record is complicated by the presence of nearby glaciers and frequent regional earthquakes. Among these signals, we analyze a specific class of short-duration, high-frequency seismic events that exhibit strong seasonal variability, increasing in rate from late summer to mid-winter before ceasing abruptly in late winter or early spring. Our analysis suggests that the source of these signals is likely near or beneath Cascade Glacier, adjacent to the landslide, rather than within the landslide mass itself. We apply detection algorithms to construct a time history for this signal type, which we then compare with environmental factors like precipitation, temperature, and slope displacement data from ground-based radar and remote sensing. Correlations indicate that these seismic events may be driven by seasonal hydrological changes, particularly the freeze-up of subglacial water pathways. While these events are not directly linked to landslide motion, they serve as indirect markers of subsurface hydrological conditions that influence slope stability. Our findings highlight the complex interplay between glaciers, groundwater, and landslide dynamics, emphasizing the need for multi-parameter monitoring to assess evolving geohazards in the region.

Alaska↗
USGS2026

Making many out of one: Synthetic geologic deformation model distributions for use in USGS NSHM25‐PRVI Puerto Rico-U.S. Virgin Island update

A key use‐case of geologic slip rates is within deformation models used in probabilistic seismic hazard analyses. Field‐derived geologic slip rates have formed the cornerstone of deformation models in such applications for decades. Recent advancements in seismic hazard analyses have expanded the use of faults for which geologic slip rates are not well constrained using categorical slip rate estimates. Because of these advancements, application of a geologic deformation model for use in 2025 U.S. Geological Survey National Seismic Hazard Model Puerto Rico‐U.S. Virgin Islands (NSHM25‐PRVI) proved challenging due to: (1) a lack of field‐based geologic slip rates, and (2) a lack of epistemic uncertainty distributions within a broad range of estimated slip rates. Preliminary versions of the NSHM25‐PRVI model sampled these slip rate bins in a coincident manner along preferred and extreme value branches, which yielded untenable correlations in mean hazard results. To minimize the influence of correlated uncertainties amid these challenges, we develop a synthetic epistemic uncertainty distribution for deformation rate on each crustal fault. Each fault has a weighting schema across four possible distribution shapes: uniform, normal, triangular favoring local minima, and triangular favoring local maxima. The synthetic distributions are then sampled several times for each logic tree branch. The results provide a more realistic distribution of rates across the study region as compared with using correlated extrema sampling. This exploration of our method in a small region like PRVI can pave the way for larger‐scale, more complicated applications (e.g., western United States).

Puerto Rico, US Virgin Islands↗
USGS2026

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗
USGS2026

Earthquake ground-motion model adjustments for the San Francisco Bay area

We develop adjustments to ergodic ground‐motion models (GMMs) to improve their performance in the San Francisco Bay Area (SFBA). GMMs are widely used in hazard assessments to estimate characteristics of ground shaking based on known properties of the source, path, and site. Such models are often developed using datasets containing records from various regions, resulting in models that represent median ground‐motion behavior, which may not adequately represent ground motions within subregions. This is true for the SFBA, where ground motions attenuate more rapidly with distance than in many other parts of California that dominate GMM databases. To support improved seismic hazard estimates in the SFBA, we calculate regional constants and anelastic attenuation coefficient adjustments relative to two commonly used ergodic GMMs: BSSA14 ( Boore et al. , 2014 ) and ASK14 ( Abrahamson et al. , 2014 ). These adjustments are obtained for a suite of ground‐motion intensity measures (peak ground acceleration, peak ground velocity, and 5%‐damped pseudospectral acceleration at oscillator periods ranging from 0.075 to 10 s) using mixed‐effects regression. Use of the regionally adjusted models reduces the overall bias by up to 0.5 natural log units for BSSA14 and up to 0.6 natural log units for ASK14. We demonstrate one application of our attenuation adjustments and their implications in an earthquake early warning case study of the 2014 M 6.0 South Napa earthquake. The predicted extent of shaking using the adjusted models better matches observed shaking at large source‐to‐site distances, especially for lower shaking intensities, thus potentially reducing overalerting. We encourage the use of our model adjustments when ergodic models are considered for seismic hazard studies in the SFBA.

California↗
USGS2025

Reflections on a trio of North American earthquakes in 1925

In 1925, three moderately large damaging earthquakes occurred in North America over four months: the 28 February (local time; LT) M 6.2 Charlevoix, 27 June (LT) M 6.6 Montana, and 29 June M 6.5 Santa Barbara earthquakes. The centennial anniversaries of these events motivated this retrospective consideration focused on the ground motions generated by the three events, including a reconsideration of early intensity assignments for the Montana earthquake. At the time, these three earthquakes appeared to support the arguments of some geologists who downplayed the severity of seismic hazard in southern California relative to other parts of the country. Some of the arguments advanced at that time, for example that Los Angeles “has the least to fear from ‘Acts of God’ of any city under the American flag,” ( Hill, 1928 ) sound naïve if not laughable now, but a comparison of well‐constrained shaking distributions for the three earthquakes reveals the dramatic difference in wave propagation efficiency in western versus eastern North America (ENAM), which leads to moderate ENAM events being felt to much larger distances. At M 6.2, the 1925 Charlevoix earthquake was a notably large event in ENAM. This earthquake was the largest event in eastern Canada since 1870 and caused damage in the epicentral region in addition to towns as far away as 200 km, with felt shaking extending over 1000 km. In contrast, felt shaking from the Santa Barbara earthquake barely extended beyond ∼200 km. Compiling published intensity distributions for larger ENAM earthquakes, we show that perceptible earthquake shaking is not uncommon in ENAM over century time scales, but experience with weakly felt shaking may incline people to downplay potential earthquake risk.

California, Montana, Quebec↗