USGS Science⌕ Search

SEARCH · USGS Science

Results for “earthquake hazards”

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 1,333 records · Page 74Linked to original sources
USGS1997

Episodic nature of earthquake activity in stable continental regions revealed by palaeoseismicity studies of Australian and North American Quaternary faults

Palaeoseismic investigations of recent faulting in stable continental regions of Australia, North America and India show that these faults typically have a long-term behaviour characterised by episodes of activity separated by quiescent intervals of at least 10 000 and commonly 100 000 years or more. Long recurrence intervals such as these are well documented by detailed studies of the faults that ruptured during the 1986 Marryat Creek, South Australia and 1988 Tennant Creek, Northern Territory earthquakes. Thus, neotectonic features associated with stable continental region faults such as scarps and grabens commonly have subtle geomorphic expression and may be poorly preserved. Many potentially hazardous faults in stable continental regions are aseismic, which is one reason why the inventory of these faults is incomplete. Although they may be currently aseismic, faults in stable continental regions that are favourably oriented for movement in the current stress field could produce damaging earthquakes, often in unexpected places. Comprehensive palaeoseismic investigations of modern and prehistoric faulting events in stable continental regions are needed to understand the long-term behaviour of these faults, and thereby, improve seismic-hazard assessments.

Australian Journal of Earth Sciences↗
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↗
USGS2017

The HayWired Earthquake Scenario

Foreword The 1906 Great San Francisco earthquake (magnitude 7.8) and the 1989 Loma Prieta earthquake (magnitude 6.9) each motivated residents of the San Francisco Bay region to build countermeasures to earthquakes into the fabric of the region. Since Loma Prieta, bay-region communities, governments, and utilities have invested tens of billions of dollars in seismic upgrades and retrofits and replacements of older buildings and infrastructure. Innovation and state-of-the-art engineering, informed by science, including novel seismic-hazard assessments, have been applied to the challenge of increasing seismic resilience throughout the bay region. However, as long as people live and work in seismically vulnerable buildings or rely on seismically vulnerable transportation and utilities, more work remains to be done. With that in mind, the U.S. Geological Survey (USGS) and its partners developed the HayWired scenario as a tool to enable further actions that can change the outcome when the next major earthquake strikes. By illuminating the likely impacts to the present-day built environment, well-constructed scenarios can and have spurred officials and citizens to take steps that change the outcomes the scenario describes, whether used to guide more realistic response and recovery exercises or to launch mitigation measures that will reduce future risk. The HayWired scenario is the latest in a series of like-minded efforts to bring a special focus onto the impacts that could occur when the Hayward Fault again ruptures through the east side of the San Francisco Bay region as it last did in 1868. Cities in the east bay along the Richmond, Oakland, and Fremont corridor would be hit hardest by earthquake ground shaking, surface fault rupture, aftershocks, and fault afterslip, but the impacts would reach throughout the bay region and far beyond. The HayWired scenario name reflects our increased reliance on the Internet and telecommunications and also alludes to the interconnectedness of infrastructure, society, and our economy. How would this earthquake scenario, striking close to Silicon Valley, impact our interconnected world in ways and at a scale we have not experienced in any previous domestic earthquake? The area of present-day Contra Costa, Alameda, and Santa Clara Counties contended with a magnitude-6.8 earthquake in 1868 on the Hayward Fault. Although sparsely populated then, about 30 people were killed and extensive property damage resulted. The question of what an earthquake like that would do today has been examined before and is now revisited in the HayWired scenario. Scientists have documented a series of prehistoric earthquakes on the Hayward Fault and are confident that the threat of a future earthquake, like that modeled in the HayWired scenario, is real and could happen at any time. The team assembled to build this scenario has brought innovative new approaches to examining the natural hazards, impacts, and consequences of such an event. Such an earthquake would also be accompanied by widespread liquefaction and landslides, which are treated in greater detail than ever before. The team also considers how the now-prototype ShakeAlert earthquake early warning system could provide useful public alerts and automatic actions. Scientific Investigations Report 2017–5013 and accompanying data releases are the products of an effort led by the USGS, but this body of work was created through the combined efforts of a large team including partners who have come together to form the HayWired Coalition (see chapter A). Use of the HayWired scenario has already begun. More than a full year of intensive partner engagement, beginning in April 2017, is being directed toward producing the most in-depth look ever at the impacts and consequences of a large earthquake on the Hayward Fault. With the HayWired scenario, our hope is to encourage and support the active ongoing engagement of the entire community of the San Francisco Bay region by providing the scientific, engineering, and economic and social science inputs for use in exercises and planning well into the future. As HayWired volumes are published, they will be made available at https://doi.org/10.3133/sir20175013 .

California↗
USGS2012

Seismic hazard exposure for the Trans-Alaska Pipeline

The discovery of oil on Alaska's North Slope and the construction of a pipeline to transport that oil across Alaska coincided with the National Environmental Policy Act of 1969 and a destructive Southern California earthquake in 1971 to cause stringent stipulations, state-of-the-art investigations, and innovative design for the pipeline. The magnitude 7.9 earthquake on the Denali fault in November 2002 was remarkably consistent with the design earthquake and fault displacement postulated for the Denali crossing of the Trans-Alaska Pipeline route. The pipeline maintained its integrity, and disaster was averted. Recent probabilistic studies to update previous hazard exposure conclusions suggest continuing pipeline integrity.

Alaska↗
USGS2025

Potential for continental scientific drilling to inform fault mechanics and earthquake science

Our understanding of fault mechanics and earthquake processes remains limited, largely due to minimal direct observations near active faults at seismogenic depths. This lack of data restricts our ability to accurately assess and mitigate both natural and human-induced seismic hazards. However, recent advancements in drilling capabilities and downhole sensing technologies offer an opportunity: the ability to observe the physical conditions within a volume near active fault zones. In this contribution, we highlight how scientific drilling can provide access to the near-fault environment, enabling measurements of the stress, temperature, fluid pressure, and rock properties at depths where ruptures initiate, propagate, and arrest. These observations are essential to refine models of earthquake nucleation and dynamic rupture, bridging gaps between laboratory experiments, numerical simulations, and surface observations. These insights can advance fundamental understanding in earthquake science but also support the development of more effective seismic hazard assessments and risk mitigation strategies.

Seismica↗
USGS2017

Reducing risk where tectonic plates collide—U.S. Geological Survey subduction zone science plan

The U.S. Geological Survey (USGS) serves the Nation by providing reliable scientific information and tools to build resilience in communities exposed to subduction zone earthquakes, tsunamis, landslides, and volcanic eruptions. Improving the application of USGS science to successfully reduce risk from these events relies on whole community efforts, with continuing partnerships among scientists and stakeholders, including researchers from universities, other government labs and private industry, land-use planners, engineers, policy-makers, emergency managers and responders, business owners, insurance providers, the media, and the general public. Motivated by recent technological advances and increased awareness of our growing vulnerability to subduction-zone hazards, the USGS is uniquely positioned to take a major step forward in the science it conducts and products it provides, building on its tradition of using long-term monitoring and research to develop effective products for hazard mitigation. This science plan provides a blueprint both for prioritizing USGS science activities and for delineating USGS interests and potential participation in subduction zone science supported by its partners. The activities in this plan address many USGS stakeholder needs: High-fidelity tools and user-tailored information that facilitate increasingly more targeted, neighborhood-scale decisions to mitigate risks more cost-effectively and ensure post-event operability. Such tools may include maps, tables, and simulated earthquake ground-motion records conveying shaking intensity and frequency. These facilitate the prioritization of retrofitting of vulnerable infrastructure; Information to guide local land-use and response planning to minimize development in likely hazardous zones (for example, databases, maps, and scenario documents to guide evacuation route planning in communities near volcanoes, along coastlines vulnerable to tsunamis, and built on landslide-prone terrain); New tools to assess the potential for cascading hazards, such as landslides, tsunamis, coastal changes, and flooding caused by earthquakes or volcanic eruptions; Geospatial models of permanent, widespread land- and sea-level changes that may occur in the immediate aftermath of great ( M ≥8.0) subduction zone earthquakes; Strong partnerships between scientists and public safety providers for effective decision making during periods of elevated hazard and risk; Accurate forecasts of far-reaching hazards (for example, ash clouds, tsunamis) to avert catastrophes and unnecessary disruptions in air and sea transportation; Aftershock forecasts to guide decisions about when and where to re-enter, repair, or rebuild buildings and infrastructure, for all types of subduction zone earthquakes.

Circular↗
USGS2021

California deepwater investigations and groundtruthing (Cal DIG) I: Fault and shallow geohazard analysis offshore Morro Bay

The California Deepwater Investigations and Groundtruthing (Cal DIG) I project focuses on the potential seafloor hazards and impacts of alternative energy infrastructure in the outer continental shelf region offshore of south-central California. This is one of three reports covering a single study area located between Monterey and Point Conception, California in federal waters outside of the State of California three nautical mile limit and in water depths of 400 to 1500 meters. The goal of this report is to provide baseline geologic interpretations of the area for the purpose of improving regional models of shallow geologic hazards and sedimentary processes. Geophysical and geological data from this project will help to address important issues associated with marine spatial planning and potential offshore infrastructure development, such as offshore floating wind turbines. Datasets covered in this report include comprehensive high-resolution sub-bottom (multi-channel and Chirp seismic reflection profiles), seafloor (bathymetry), and sampling (piston, gravity, and vibracore) data collected in 2018-2019 during a series of seven seagoing geological and geophysical surveys. Specifically, this report outlines interpretation of subsurface geologic structure from the geophysical data, details preliminary core analysis results related to fluid, gas, and sediment transport activity, provides interpretations of the current geohazards in the area, and suggests next steps for improving interpretations of geohazard processes. Specific targets of geohazard interest in the study area are geological structures such as faults and folds, seafloor pockmarks within a large field (the Big Sur pockmark field), submarine channels, and mass wasting (slope failure) features. The vast majority of faults and other structures in the study occur within sediment and rock formations we interpret to be pre-Quaternary (older than 2.58 Myr BP), and thus we interpret that these structures are unlikely to present significant current hazard to seabed infrastructure, although we note that the numerous structures mapped in the study area may have the potential to become reactivated. Similarly, we find no new evidence of Holocene (younger than 11,650 years BP) fluid or gas advection in the Big Sur pockmark field. However, such fluid and gas hazards are currently difficult to assess, as additional analyses and sampling of existing core data is needed to better understand pockmark formation processes and potential gas accumulations we have mapped in the subsurface. Mass wasting along the eastern and western edges of the Santa Lucia Bank during earthquakes, as well as sediment transport down the Lucia Chica and San Simeon channels, are among the most significant, although still likely infrequent during the Holocene, hazards to seabed stability in the study area. Further analyses of the existing cores, including radiocarbon dating, stable isotope, and compositional analyses, are again needed to better understand the timing and sources of the numerous sand deposits found throughout the study area, which may have been transported downslope due to mass wasting and/or earthquake shaking processes.

California↗
USGS2007

GIS data for the Seaside, Oregon, Tsunami Pilot Study to modernize FEMA flood hazard maps

A Tsunami Pilot Study was conducted for the area surrounding the coastal town of Seaside, Oregon, as part of the Federal Emergency Management's (FEMA) Flood Insurance Rate Map Modernization Program (Tsunami Pilot Study Working Group, 2006). The Cascadia subduction zone extends from Cape Mendocino, California, to Vancouver Island, Canada. The Seaside area was chosen because it is typical of many coastal communities subject to tsunamis generated by far- and near-field (Cascadia) earthquakes. Two goals of the pilot study were to develop probabilistic 100-year and 500-year tsunami inundation maps using Probabilistic Tsunami Hazard Analysis (PTHA) and to provide recommendations for improving tsunami hazard assessment guidelines for FEMA and state and local agencies. The study was an interagency effort by the National Oceanic and Atmospheric Administration, U.S. Geological Survey, and FEMA, in collaboration with the University of Southern California, Middle East Technical University, Portland State University, Horning Geoscience, Northwest Hydraulics Consultants, and the Oregon Department of Geological and Mineral Industries. The pilot study model data and results are published separately as a geographic information systems (GIS) data report (Wong and others, 2006). The flood maps and GIS data are briefly described here.

Oregon↗
USGS2002

Differences in attenuation among the stable continental regions

There are systematic differences in the attenuation of damaging earthquake ground motions between different stable continental regions (SCRs). Seismic intensity and weak-motion data show that the attenuation in seismic waves for eastern North America (ENA) is less than for India, Africa, Australia, and northwest Europe. If ENA ground-motion attenuation relations are used in seismic hazard models for other SCRs, as is commonly done, then the estimated ground motions and resulting hazard may be too large. If an attenuation model that averages observations from ENA and the other SCRs is used to estimate the magnitudes of large historical earthquakes in ENA, as is the case for recent estimates of M for the 1811-1812 New Madrid, Missouri and the 1886 Charleston, South Carolina events, then the magnitude estimates for these events will be too large, as will be the resulting hazard.

Geophysical Research Letters↗
USGS2022

A study on the effect of site response on California seismic hazard map assessment

Prior studies have repeatedly shown that probabilistic seismic hazard maps from several different countries predict higher shaking than that observed. Previous map assessments have not, however, considered the influence of site response on hazard. Seismologists have long acknowledged the influence of near-surface geology, in particular low-impedance sediment layers, on earthquake ground-motion at frequencies of engineering concern. Although the overall effects of site response are complex, modern ground-motion models (GMMs) account for site effects using terms based on V S30 , the time-averaged shear-wave velocity in the upper 30 m of the Earth’s surface. In this study, we consider general implications of incorporating site terms from modern GMMs using site-specific V S30 as a proxy in probabilistic seismic hazard maps for California. At the long periods (1–5 s) that affect tall buildings, site terms amplify the mapped hazard by factors of 1–3 at many sites relative to maps calculated for the standard reference soft-rock site condition, V S30 = 760 m/s. However, at the short periods of ground-motion that are the main contributors to peak ground acceleration (PGA) and thus affect smaller structures, only negligible effects occur due to nonlinear deamplification of strong ground-motion at high frequencies. Nonlinear deamplification increases as the shaking level increases. For very strong shaking, deamplification can overcome the linear amplification, yielding net deamplification. We explore the implications of these results for the evaluation of hazard maps. Because site effects do not change the maps appreciably at short periods, we can exclude site response as an explanation for why the maps overpredict historically observed shaking as captured by the California Historical Intensity Mapping Project (CHIMP) dataset. The results are expected to be generalizable to regions that are comparable to California in terms of structure and seismicity rates. In low-to-moderate-seismicity regions where the hazard reflects weaker shaking, nonlinear site response is expected to be less important for the hazard.

California↗
USGS1990

Geology and land use

Geologists' eyes are trained to find and trace such natural landmarks as flood plains, landslide scars, retreating shoreline bluffs, or surface traces of active earthquake faults. more and more often, in developing areas, we find these obvious signs of trouble being erased by urban development. A geological hazard concealed by landscaping or hosing is fully as dangerous as when it is visible. The geologic limitations for building sites of some areas can be overcome, in part, by skilled engineering and expensive construction practices. But the costs can be prohibitively high, and the solutions are not always completely effective. In "earthquake country," history has shown that costs are highest and risk factors most uncertain in a few easily recognized settings: unstable hill sloped, land at the edge of rapidly eroding sea cliffs, lowlands underlain by saturated estuarine mud of ill, and areas near faults capable of producing magnitude 7 or greater earthquakes. Safety immediately after an earthquake is also a concern in these places, for extreme damage and ground distortion may impede or prevent timely access by emergency equipment.

Earthquakes & Volcanoes (USGS)↗
USGS2015

Crustal deformation in the New Madrid seismic zone and the role of postseismic processes

Global Navigation Satellite System data across the New Madrid seismic zone (NMSZ) in the central United States over the period from 2000 through 2014 are analyzed and modeled with several deformation mechanisms including the following: (1) creep on subsurface dislocations, (2) postseismic frictional afterslip and viscoelastic relaxation from the 1811–1812 and 1450 earthquakes in the NMSZ, and (3) regional strain. In agreement with previous studies, a dislocation creeping at about 4 mm/yr between 12 and 20 km depth along the downdip extension of the Reelfoot fault reproduces the observations well. We find that a dynamic model of postseismic frictional afterslip from the 1450 and February 1812 Reelfoot fault events can explain this creep. Kinematic and dynamic models involving the Cottonwood Grove fault provide minimal predictive power. This is likely due to the smaller size of the December 1811 event on the Cottonwood Grove fault and a distribution of stations better suited to constrain localized strain across the Reelfoot fault. Regional compressive strain across the NMSZ is found to be less than 3 × 10 −9 /yr. If much of the present-day surface deformation results from afterslip, it is likely that many of the earthquakes we see today in the NMSZ are aftershocks from the 1811–1812 New Madrid earthquakes. Despite this conclusion, our results are consistent with observations and models of intraplate earthquake clustering. Given this and the recent paleoseismic history of the region, we suggest that seismic hazard is likely to remain significant.

Journal of Geophysical Research B: Solid Earth↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗
USGS2014

Bayesian historical earthquake relocation: an example from the 1909 Taipei earthquake

Locating earthquakes from the beginning of the modern instrumental period is complicated by the fact that there are few good-quality seismograms and what traveltimes do exist may be corrupted by both large phase-pick errors and clock errors. Here, we outline a Bayesian approach to simultaneous inference of not only the hypocentre location but also the clock errors at each station and the origin time of the earthquake. This methodology improves the solution for the source location and also provides an uncertainty analysis on all of the parameters included in the inversion. As an example, we applied this Bayesian approach to the well-studied 1909 M w 7 Taipei earthquake. While our epicentre location and origin time for the 1909 Taipei earthquake are consistent with earlier studies, our focal depth is significantly shallower suggesting a higher seismic hazard to the populous Taipei metropolitan area than previously supposed.

Taipei↗
USGS2019

Variability in synthetic earthquake ground motions caused by source variability and errors in wave propagation models

Numerical simulations of earthquake ground motions are used both to anticipate the effects of hypothetical earthquakes by forward simulation and to infer the behaviour of the real earthquake source ruptures by the inversion of recorded ground motions. In either application it is necessary to assume some Earth structure that is necessarily inaccurate and to use a computational method that is also inaccurate for simulating the wavefield Green's functions. We refer to these two sources of error as ‘propagation inaccuracies’, which might be considered to be epistemic. We show that the variance of the Fourier spectrum of the synthetic earthquake seismograms caused by propagation inaccuracies is related to the spatial covariance on the rupture surface of errors in the computed Green's functions, which we estimate for the case of the 2009 L'Aquila, Italy, earthquake by comparing erroneous computed Green's functions with observed L'Aquila aftershock seismograms (empirical Green's functions). We further show that the variance of the synthetic seismograms caused by the rupture variability (aleatory uncertainty) is related to the spatial covariance on the rupture surface of aleatory variations in the rupture model, and we investigate the effect of correlated variations in Green's function errors and variations in rupture models. Thus, we completely characterize the variability of synthetic earthquake seismograms induced by errors in propagation and variability in the rupture behaviour. We calculate the spectra of the variance of the ground motions of the L'Aquila main shock caused by propagation inaccuracies for two specific broad-band stations, the AQU and the FIAM stations. These variances are distressingly large, being comparable or in some cases exceeding the data amplitudes, suggesting that the best-fitting L'Aquila rupture model significantly overfits the data and might be seriously in error. If these computed variances are typical, the accuracy of many other rupture models for past earthquakes may need to be reconsidered. The results of this work might be useful in seismic hazard estimation because the variability of the computed ground motion, caused both by propagation inaccuracies and variations in the rupture model, can be computed directly, not requiring laborious consideration of multiple Earth structures.

Geophysical Journal International↗
USGS2007

Integrated geologic and geophysical studies of North American continental intraplate seismicity

The origin of earthquakes within stable continental regions has been the subject of debate over the past thirty years. Here, we examine the correlation of North American stable continental region earthquakes using five geologic and geophysical data sets: (1) a newly compiled age-province map; (2) Bouguer gravity data; (3) aeromagnetic anomalies; (4) the tectonic stress field; and (5) crustal structure as revealed by deep seismic-reflection profiles. We find that: (1) Archean-age (3.8–2.5 Ga) North American crust is essentially aseismic, whereas post-Archean (less than 2.5 Ga) crust shows no clear correlation of crustal age and earthquake frequency or moment release; (2) seismicity is correlated with continental paleorifts; and (3) seismicity is correlated with the NE-SW structural grain of the crust of eastern North America, which in turn reflects the opening and closing of the proto– and modern Atlantic Ocean. This structural grain can be discerned as clear NE-SW lineaments in the Bouguer gravity and aeromagnetic anomaly maps. Stable continental region seismicity either: (1) follows the NE-SW lineaments; (2) is aligned at right angles to these lineaments; or (3) forms clusters at what have been termed stress concentrators (e.g., igneous intrusions and intersecting faults). Seismicity levels are very low to the west of the Grenville Front (i.e., in the Archean Superior craton). The correlation of seismicity with NE-SW–oriented lineaments implies that some stable continental region seismicity is related to the accretion and rifting processes that have formed the North American continental crust during the past 2 b.y. We further evaluate this hypothesis by correlating stable continental region seismicity with recently obtained deep seismicreflection images of the Appalachian and Grenville crust of southern Canada. These images show numerous faults that penetrate deep (40 km) into the crust. An analysis of hypocentral depths for stable continental region earthquakes shows that the frequency and moment magnitude of events are nearly uniform for the entire 0–35 km depths over which crustal earthquakes extend. This is in contradiction with the hypothesis that larger events have deeper focal depths. We conclude that the deep structure of the crust, in particular the existence of deeply penetrating faults, is the controlling parameter, rather than lateral variations in temperature, rheology, or high pore pressure. The distribution of stable continental region earthquakes in eastern North America is consistent with the existence of deeply penetrating crustal faults that have been reactivated in the present stress field. We infer that future earthquakes may occur anywhere along the geophysical lineations that we have identified. This implies that seismic hazard is more widespread in central and eastern North America than indicated by the limited known historical distribution of seismicity.

Special Paper of the Geological Society of America↗
USGS2026

Deformed submarine terraces in Puget Sound, Pacific Northwest, indicate only one M >~7.5 earthquake on the Seattle fault zone in the past 11,000 yr

Submerged marine terraces in Puget Sound, deformed across the Seattle fault zone (SFZ), indicate that only one earthquake as large as M~7.5 has occurred in at least the past 11 kyr. Previous paleoseismic studies document a M~7.5 earthquake between 923–4 CE, which uplifted coastal marine terraces by as much as 8 m. We demonstrate that this earthquake was the only such event since ~11 ka by mapping and quantifying deformation of older marine terraces, now submerged in Puget Sound. The submerged terraces, attributed to a late-glacial sea-level lowstand, record both glacial isostatic rebound and tectonic deformation. Vertical offset of the ~11 ka terraces within the SFZ is comparable to that of the marine terraces uplifted in 923 CE, implying no additional large (M>~7.5) earthquake on the SFZ since ~11 ka. This result implies a longer recurrence interval than current hazard estimates, which assumes recurrence of M>7.1 events every 5 kyr. Our mapping of SFZ deformation since ~11 ka also supports fault segmentation and contiguous block uplift between the Seattle and Tacoma fault zones.

Washington↗
USGS2000

Discrepancy between earthquake rates implied by historic earthquakes and a consensus geologic source model for California

We examine the difference between expected earthquake rates inferred from the historical earthquake catalog and the geologic data that was used to develop the consensus seismic source characterization for the state of California [California Department of Conservation, Division of Mines and Geology (CDMG) and U.S. Geological Survey (USGS) Petersen et al., 1996; Frankel et al., 1996]. On average the historic earthquake catalog and the seismic source model both indicate about one M 6 or greater earthquake per year in the state of California. However, the overall earthquake rates of earthquakes with magnitudes (M) between 6 and 7 in this seismic source model are higher, by at least a factor of 2, than the mean historic earthquake rates for both southern and northern California. The earthquake rate discrepancy results from a seismic source model that includes earthquakes with characteristic (maximum) magnitudes that are primarily between M 6.4 and 7.1. Many of these faults are interpreted to accommodate high strain rates from geologic and geodetic data but have not ruptured in large earthquakes during historic time. Our sensitivity study indicates that the rate differences between magnitudes 6 and 7 can be reduced by adjusting the magnitude-frequency distribution of the source model to reflect more characteristic behavior, by decreasing the moment rate available for seismogenic slip along faults, by increasing the maximum magnitude of the earthquake on a fault, or by decreasing the maximum magnitude of the background seismicity. However, no single parameter can be adjusted, consistent with scientific consensus, to eliminate the earthquake rate discrepancy. Applying a combination of these parametric adjustments yields an alternative earthquake source model that is more compatible with the historic data. The 475-year return period hazard for peak ground and 1-sec spectral acceleration resulting from this alternative source model differs from the hazard resulting from the standard CDMG-USGS model by less than 10% across most of California but is higher (generally about 10% to 30%) within 20 km from some faults.

Bulletin of the Seismological Society of America↗