Search USGS⌕ Search

SEARCH · Search USGS

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 577 records · Page 32Linked to original sources
USGS1975

Preliminary report on the reconnaissance engineering geology of the Yakutat area, Alaska, with emphasis on evaluation of earthquake and other geologic hazards

Yakutat, situated about 225 miles northwest of Juneau, Alaska, near the shores of the Gulf of Alaska, has a setting that calls for superlatives. Within the Yakutat region are some of the tallest mountains, some of the heaviest snowfalls, and the largest glacier in North America. Between the abrupt mountain front and the Gulf of Alaska lies a very gently sloping plain of outwash derived from repeated cycles of advance followed by melt and retreat of glaciers during the Quaternary Period. The latest melting probably took place 500 to 600 years ago. Yakutat is built upon the moderately steep moraine that is the product of melting of one of these glaciers. Near Yakutat, surficial deposits are as much as 700 feet thick and probably overlie siltstone, sandstone, and mudstone. The eight general categories of mapped surficial deposits include artificial fill, organic, eolian, beach, delta-estuarine, alluvial, and outwash deposits, and deposits of the outer Yakutat Bay moraine complex.

Alaska↗
USGS2024

Using open-science workflow tools to produce SCEC CyberShake physics-based probabilistic seismic hazard models

The Statewide (formerly Southern) California Earthquake Center (SCEC) conducts multidisciplinary earthquake system science research that aims to develop predictive models of earthquake processes, and to produce accurate seismic hazard information that can improve societal preparedness and resiliency to earthquake hazards. As part of this program, SCEC has developed the CyberShake platform, which calculates physics-based probabilistic seismic hazard analysis (PSHA) models for regions with high-quality seismic velocity and fault models. The CyberShake platform implements a sophisticated computational workflow that includes over 15 individual codes written by 6 developers. These codes are heterogeneous, ranging from short-running high-throughput serial CPU codes to large, long-running, parallel GPU codes. Additionally, CyberShake simulation campaigns are computationally extensive, typically producing tens of terabytes of meaningful scientific data and metadata over several months of around-the-clock execution on leadership-class supercomputers. To meet the needs of the CyberShake platform, we have developed an extreme-scale workflow stack, including the Pegasus Workflow Management System, HTCondor, Globus, and custom tools. We present this workflow software stack and identify how the CyberShake platform and supporting tools enable us to meet a variety of challenges that come with large-scale simulations, such as automated remote job submission, data management, and verification and validation. This platform enabled us to perform our most recent simulation campaign, CyberShake Study 22.12, from December 2022 to April 2023. During this time, our workflow tools executed approximately 32,000 jobs, and used up to 73% of the Summit system at Oak Ridge Leadership Computing Facility. Our workflow tools managed about 2.5 PB of total temporary and output data, and automatically staged 19 million output files totaling 74 TB back to archival storage on the University of Southern California's Center for Advanced Research Computing systems, including file-based relational data and large binary files to efficiently store millions of simulated seismograms. CyberShake extreme-scale workflows have generated simulation-based probabilistic seismic hazard models that are being used by seismological, engineering, and governmental communities.

Frontiers Earth Science Journal↗
USGS2004

Apparent stress, fault maturity and seismic hazard for normal-fault earthquakes at subduction zones

The behavior of apparent stress for normal-fault earthquakes at subduction zones is derived by examining the apparent stress (?? a = ??Es/Mo, where E s is radiated energy and Mo is seismic moment) of all globally distributed shallow (depth, ?? < 70 km) earthquakes with normal-fault mechanisms that occurred in or near subduction zones between 1987 and 2001 for which Es and M o are available. Accurately determined hypocentres from the Engdahl-Hilst-Buland (EHB) catalogue establish the fine detail of the Wadati-Benioff zone. In many cases, we can relate trends in apparent stress to specific features within the subduction zone and compare these with trends for interplate-thrust earthquakes in the same subduction zones. There are two depth ranges over which ??a, for normal-fault earthquakes shows maxima. The highest and most anomalous values of ??a are found in the deeper depth range from 35 to 70 km. The high-??a, events (up to 5 MPa) are characteristically intraslab and in proximity to zones of intense deformation such as a sharp slab bend or where opposing slabs collide. High-??a events in the region of the shallower maximum (hypocentres between 10-35 km and ??a > 1 MPa) are also generally intraslab, but occur where the lithosphere has just begun subduction beneath the overriding plate. They usually occur in cold slabs near trenches where the direction of plate motion across the trench is oblique to the trench axis, or where there are local contortions or geometrical complexities of the plate boundary. Lower ??a (< 1 MPa) is associated with events occurring at the outer rise (OR) complex (between the OR and the trench axis), as well as with intracrustal events occurring just landward of the trench. The average apparent stress of intraslab-normal-fault earthquakes is considerably higher than the average apparent stress of interplate-thrust-fault earthquakes. In turn, the average ?? a of strike-slip earthquakes in intraoceanic environments is considerably higher than that of intraslab-normal-fault earthquakes. The variation of average ??a with focal mechanism and tectonic regime suggests that the level of ?? a is related to fault maturity. Lower stress drops are needed to rupture mature faults such as those found at plate interfaces that have been smoothed by large cumulative displacements (from hundreds to thousands of kilometres). In contrast, immature faults, such as those on which intraslab-normal-fault earthquakes generally occur, are found in cold and intact lithosphere in which total fault displacement has been much less (from hundreds of metres to a few kilometres). Also, faults on which high ??a oceanic strike-slip earthquakes occur are predominantly intraplate or at evolving ends of transforms. At subduction zones, earthquakes occurring on immature faults are likely to be more hazardous as they tend to generate higher amounts of radiated energy per unit of moment than earthquakes occurring on mature faults. We have identified earthquake pairs in which an interplate-thrust and an intraslab-normal earthquake occurred remarkably close in space and time. The intraslab-normal member of each pair radiated anomalously high amounts of energy compared to its thrust-fault counterpart. These intraslab earthquakes probably ruptured intact slab mantle and are dramatic examples in which Mc (an energy magnitude) is shown to be a far better estimate of the potential for earthquake damage than Mw. This discovery may help explain why loss of life as a result of intraslab earthquakes was greater in the 20th century in Latin America than the fatalities associated with interplate-thrust events that represented much higher total moment release. ?? 2004 RAS.

Geophysical Journal International↗
USGS2018

Variabilities in probabilistic seismic hazard maps for natural and induced seismicity in the central and eastern United States

Probabilistic seismic hazard analysis (PSHA) characterizes ground-motion hazard from earthquakes. Typically, the time horizon of a PSHA forecast is long, but in response to induced seismicity related to hydrocarbon development, the USGS developed one-year PSHA models. In this paper, we present a display of the variability in USGS hazard curves due to epistemic uncertainty in its informed submodel using a simple bootstrapping approach. We find that variability is highest in low-seismicity areas. On the other hand, areas of high seismic hazard, such as the New Madrid seismic zone or Oklahoma, exhibit relatively lower variability simply because of more available data and a better understanding of the seismicity. Comparing areas of high hazard, New Madrid, which has a history of large naturally occurring earthquakes, has lower forecast variability than Oklahoma, where the hazard is driven mainly by suspected induced earthquakes since 2009. Overall, the mean hazard obtained from bootstrapping is close to the published model, and variability increased in the 2017 one-year model relative to the 2016 model. Comparing the relative variations caused by individual logic-tree branches, we find that the highest hazard variation (as measured by the 95% confidence interval of bootstrapping samples) in the final model is associated with different ground-motion models and maximum magnitudes used in the logic tree, while the variability due to the smoothing distance is minimal. It should be pointed out that this study is not looking at the uncertainty in the hazard in general, but only as it is represented in the USGS one-year models.

The Leading Edge↗
USGS2011

Seismic hazard maps for Haiti

We have produced probabilistic seismic hazard maps of Haiti for peak ground acceleration and response spectral accelerations that include the hazard from the major crustal faults, subduction zones, and background earthquakes. The hazard from the Enriquillo-Plantain Garden, Septentrional, and Matheux-Neiba fault zones was estimated using fault slip rates determined from GPS measurements. The hazard from the subduction zones along the northern and southeastern coasts of Hispaniola was calculated from slip rates derived from GPS data and the overall plate motion. Hazard maps were made for a firm-rock site condition and for a grid of shallow shear-wave velocities estimated from topographic slope. The maps show substantial hazard throughout Haiti, with the highest hazard in Haiti along the Enriquillo-Plantain Garden and Septentrional fault zones. The Matheux-Neiba Fault exhibits high hazard in the maps for 2% probability of exceedance in 50 years, although its slip rate is poorly constrained.

Earthquake Spectra↗
USGS2025

A soil velocity model for improved ground motion simulations in the U. S. Pacific Northwest

Near-surface seismic velocity structure may significantly impact the intensity, duration, and frequency content of ground shaking during an earthquake. In this study, we compile 649 shear wave velocity (Vs) profiles throughout the U.S. Pacific Northwest and southern British Columbia (PNW) and use these measured profiles to develop a representative soil velocity model for four major Holocene soil provinces: Puget Lowlands, Willamette Valley, fill and alluvium, and `other' soils. The resulting soil velocity model shows good agreement to measured data for a wide range of site conditions, with variability between different geologic domains reflecting fundamental differences in depositional environments. We then show that using this regional soil velocity model in simulations of the 2001 M6.8 Nisqually, Washington earthquake improves the fit to observed high-frequency (≥ 0.5 Hz) ground motions in the Puget Sound region compared to simulations that do not incorporate shallow (≤ 200 m) seismic velocity structure. Overall, this work shows that incorporating localized soil velocity profiles into seismic velocity models is important for accurately estimating high-frequency ground motion and regional seismic hazard in earthquake simulations. Future earthquake simulations and hazard studies in the PNW could incorporate these soil velocity profiles to capture the region's distinct site response characteristics.

Washington↗
USGS2012

Evaluating fault rupture hazard for strike-slip earthquakes

We present fault displacement data, regressions, and a methodology to calculate in both a probabilistic and deterministic framework the fault rupture hazard for strike-slip faults. To assess this hazard we consider: (1) the size of the earthquake and probability that it will rupture to the surface, (2) the rate of all potential earthquakes on the fault (3) the distance of the site along and from the mapped fault, (4) the complexity of the fault and quality of the fault mapping, (5) the size of the structure that will be placed at the site, and (6) the potential and size of displacements along or near the fault. Probabilistic fault rupture hazard analysis should be an important consideration in design of structures or lifelines that are located within about 50m of well-mapped active faults.

Conference Paper↗
USGS2024

Earthquake scenario development in conjunction with the 2023 USGS National Seismic Hazard Model

We present earthquake scenarios developed to accompany the release of the 2023 update to the US Geological Survey National Seismic Hazard Model (NSHM). Scenarios can serve a range of local and regional needs, from developing proactive-targeted mitigation strategies for minimizing impending risk to aiding emergency management planning. These deterministic scenarios can also be used to communicate seismic hazard and risk to audiences who are not well versed in methods, such as probabilistic seismic hazard analyses. Specifically, we discuss the scenarios developed, challenges, and lessons learned in the development process, and how this work aided the development of the 2023 NSHM itself. In total, 28 scenarios were developed for Hawaii, Utah, Alaska, and Virginia considering the 2023 NSHM science, past scenario efforts, and input from local experts and stakeholders. Finally, we investigate how NSHM modeling decisions can change estimated impacts to Utah and Hawaii in more detail showing, for example, that a shallower dip of the Wasatch fault under Salt Lake City can increase predicted ground-motion intensities and therefore estimated losses and deaths.

Earthquake Spectra↗
USGS2010

Locations and magnitudes of historical earthquakes in the Sierra of Ecuador (1587–1996)

The whole territory of Ecuador is exposed to seismic hazard. Great earthquakes can occur in the subduction zone (e.g. Esmeraldas, 1906, M w 8.8), whereas lower magnitude but shallower and potentially more destructive earthquakes can occur in the highlands. This study focuses on the historical crustal earthquakes of the Andean Cordillera. Several large cities are located in the Interandean Valley, among them Quito, the capital (∼2.5 millions inhabitants). A total population of ∼6 millions inhabitants currently live in the highlands, raising the seismic risk. At present, precise instrumental data for the Ecuadorian territory is not available for periods earlier than 1990 (beginning date of the revised instrumental Ecuadorian seismic catalogue); therefore historical data are of utmost importance for assessing seismic hazard. In this study, the Bakun & Wentworth method is applied in order to determine magnitudes, locations, and associated uncertainties for historical earthquakes of the Sierra over the period 1587–1976. An intensity-magnitude equation is derived from the four most reliable instrumental earthquakes ( M w between 5.3 and 7.1). Intensity data available per historical earthquake vary between 10 (Quito, 1587, Intensity ≥VI) and 117 (Riobamba, 1797, Intensity ≥III). The bootstrap resampling technique is coupled to the B&W method for deriving geographical confidence contours for the intensity centre depending on the data set of each earthquake, as well as confidence intervals for the magnitude. The extension of the area delineating the intensity centre location at the 67 per cent confidence level (±1σ) depends on the amount of intensity data, on their internal coherence, on the number of intensity degrees available, and on their spatial distribution. Special attention is dedicated to the few earthquakes described by intensities reaching IX, X and XI degrees. Twenty-five events are studied, and nineteen new epicentral locations are obtained, yielding equivalent moment magnitudes between 5.0 and 7.6. Large earthquakes seem to be related to strike slip faults between the North Andean Block and stable South America to the east, while moderate earthquakes ( M w ≤ 6) seem to be associated with to thrust faults located on the western internal slopes of the Interandean Valley.

Geophysical Journal International↗
USGS2018

Evaluation of ground motion models for USGS seismic hazard forecasts: Induced and tectonic earthquakes in the Central and Eastern U.S.

Ground motion model (GMM) selection and weighting introduces a significant source of uncertainty in United States Geological Survey (USGS) seismic hazard models. The increase in moderate moment magnitude induced earthquakes (Mw 4 to 5.8) in Oklahoma and Kansas since 2009, due to increased wastewater injection related to oil and gas production (Keranen et al., 2013; 2014; Weingarten et al., 2015; McNamara et al., 2015a), provides useful near-source (< 40 km) instrumental ground-motion observations for comparisons between central and eastern US (CEUS) induced (Rennolet et al., 2017) and tectonic (Goulet et al., 2014) earthquakes. In this study, we evaluate over 50 GMMs using two well-established probabilistic scoring methods: log likelihood (LLH) (Scherbaum et al., 2004; 2009) and multivariate LLH (MLLH) (Mak et al., 2017). The LLH approach compares the mean and standard deviation (σ) of the observed and modeled ground motions. The MLLH approach advances the LLH method by considering the variability (φ,τ) of multiple correlated variables namely intra- (within) and inter- (between) event residuals. For the probabilistic scoring GMM evaluation methods (LLH, MLLH), we compute horizontal component peak ground acceleration (PGA) and 1s period pseudo spectral acceleration (PSA1.0) total residuals using GMM software (nshmp-haz) recently implemented by the USGS National Seismic Hazard Model Project (NSHMP). We observe from LLH and MLLH scores that: 1) newer GMMs with lower standard deviations (σ,φ,τ) score better than older GMMs with higher published uncertainty; 2) 2014 CEUS GMMs score better for CEUS tectonic earthquakes than induced earthquakes; 3) NGA-West2, G17 and A15 GMMs score well for CEUS induced earthquake ground motions; and 4) NGA-East GMMs score well for CEUS tectonic earthquake ground motions. We also use the LLH and MLLH scores to evaluate GMM weights applied in past USGS seismic hazard forecasts and to inform weighting of GMMs in future seismic hazard forecasts.

Bulletin of the Seismological Society of America↗
USGS2012

Deep-Sea Turbidites as Guides to Holocene Earthquake History at the Cascadia Subduction Zone&mdash;Alternative Views for a Seismic-Hazard Workshop

This report reviews the geological basis for some recent estimates of earthquake hazards in the Cascadia region between southern British Columbia and northern California. The largest earthquakes to which the region is prone are in the range of magnitude 8-9. The source of these great earthquakes is the fault down which the oceanic Juan de Fuca Plate is being subducted or thrust beneath the North American Plate. Geologic evidence for their occurrence includes sedimentary deposits that have been observed in cores from deep-sea channels and fans. Earthquakes can initiate subaqueous slumps or slides that generate turbidity currents and which produce the sedimentary deposits known as turbidites. The hazard estimates reviewed in this report are derived mainly from deep-sea turbidites that have been interpreted as proxy records of great Cascadia earthquakes. The estimates were first published in 2008. Most of the evidence for them is contained in a monograph now in press. We have reviewed a small part of this evidence, chiefly from Cascadia Channel and its tributaries, all of which head offshore the Pacific coast of Washington State. According to the recent estimates, the Cascadia plate boundary ruptured along its full length in 19 or 20 earthquakes of magnitude 9 in the past 10,000 years; its northern third broke during these giant earthquakes only, and southern segments produced at least 20 additional, lesser earthquakes of Holocene age. The turbidite case for full-length ruptures depends on stratigraphic evidence for simultaneous shaking at the heads of multiple submarine canyons. The simultaneity has been inferred primarily from turbidite counts above a stratigraphic datum, sandy beds likened to strong-motion records, and radiocarbon ages adjusted for turbidity-current erosion. In alternatives proposed here, this turbidite evidence for simultaneous shaking is less sensitive to earthquake size and frequency than previously thought. Turbidites far below a channel confluence, instead of representing the merged flows from two tributaries, monitor the dominant tributary only. Sandy beds low in the turbidites, instead of matching from channel to channel, permit divergent stratigraphic correlations; and rather than approximating strong-motion seismograms, the sandy beds more likely record processes internal to the generation and transformation of subaqueous mass movements. The age adjustments, instead of supporting other evidence that all the northern ruptures were long, are uncertain enough to accord with variation in rupture mode, and this variation improves agreement with onshore paleoseismology. Many of the turbidites counted as evidence for frequent earthquakes on the southern Cascadia plate boundary may instead reflect nearness to steep slopes. This report is meant to aid in the updating of national maps of seismic hazards in Canada and the United States. It offers three main conclusions for consideration at a U.S. hazard-map workshop slated for March 21-22, 2012: If giant earthquakes are the norm for the plate boundary offshore southern Washington, the strongest paleoseismic evidence for this rupture mode is the average earthquake-recurrence interval of about 500 years that is evidenced both offshore in lower Cascadia Channel and onshore at estuaries of southern Washington and northernmost Oregon. The plate boundary offshore southern British Columbia and northern Washington may be capable of producing great earthquakes at an average interval as short as 300 years that is evidenced mainly onshore. Review of more of the turbidite evidence now in press may clarify implications for the hazard maps. Further work on the deep-sea turbidites could target sedimentary processes and chronological uncertainties that may affect the turbidites' sensitivity to fault-rupture lengths and recurrence rates.

Washington;Oregon;California↗
USGS2000

Geologic Map and Map Database of the Oakland Metropolitan Area, Alameda, Contra Costa, and San Francisco Counties, California

Introduction This report contains a new geologic map at 1:50,000 scale, derived from a set of geologic map databases containing information at a resolution associated with 1:24,000 scale, and a new description of geologic map units and structural relationships in the mapped area. The map database represents the integration of previously published reports and new geologic mapping and field checking by the author (see Sources of Data index map on the map sheet or the Arc-Info coverage pi-so and the textfile pi-so.txt). The descriptive text (below) contains new ideas about the Hayward fault and other faults in the East Bay fault system, as well as new ideas about the geologic units and their relations. These new data are released in digital form in conjunction with the Federal Emergency Management Agency Project Impact in Oakland. The goal of Project Impact is to use geologic information in land-use and emergency services planning to reduce the losses occurring during earthquakes, landslides, and other hazardous geologic events. The USGS, California Division of Mines and Geology, FEMA, California Office of Emergency Services, and City of Oakland participated in the cooperative project. The geologic data in this report were provided in pre-release form to other Project Impact scientists, and served as one of the basic data layers for the analysis of hazard related to earthquake shaking, liquifaction, earthquake induced landsliding, and rainfall induced landsliding. The publication of these data provides an opportunity for regional planners, local, state, and federal agencies, teachers, consultants, and others outside Project Impact who are interested in geologic data to have the new data long before a traditional paper map could be published. Because the database contains information about both the bedrock and surficial deposits, it has practical applications in the study of groundwater and engineering of hillside materials, as well as the study of geologic hazards and the academic research on the geologic history and development of the region.

Miscellaneous Field Studies Map↗
USGS2023

Rapid shallow megathrust afterslip from the 2021 M8.2 Chignik, Alaska earthquake revealed by seafloor geodesy

The shallower portions of subduction zone megathrust faults host Earth’s most hazardous tsunamigenic earthquakes, yet understanding how and when they slip remains elusive because of challenges making seafloor observations. We performed Global Navigation Satellite System Acoustic seafloor geodetic surveys before and ~2.5 months after the 29 July 2021 M w (moment magnitude) 8.2 Chignik, Alaska, earthquake and determine ~1.4 meters cumulative co- and post-seismic horizontal displacement ~60 kilometers from the megathrust front. Only for the 2011 M w 9 Tohoku event have closer subduction zone earthquake displacements been observed. We estimate ~2 to 3 meters of megathrust afterslip shallower than 20 kilometers, a portion of the megathrust on which both inter- and co-seismic slip likely had occurred previously. Our analysis demonstrates that by 2.5 months, shallower and deeper moment had effectively equilibrated on the megathrust, suggesting that its tsunamigenic potential remains no more elevated than before the earthquake.

Alaska↗
USGS2005

Hazard calculations for the WGCEP-2002 earthquake forecast using OpenSHA and distributed object technologies

We present seismic-hazard calculations for what is arguably the most sophisticated earthquake forecast ever developed—the model by the 2002 Working Group on California Earthquake Probabilities (2003), or WGCEP-2002 as referred to hereafter. These calculations have been made possible by developments in both OpenSHA ( Field et al. , 2003 ) and the Information Technology Research (ITR) Collaboration of the Southern California Earthquake Center (SCEC) ( Jordan et al. , 2003 ). In particular, we demonstrate the use of a freely available, platform-independent, and graphical-user-interface-based application for computing hazard curves. This application utilizes distributed-object technologies, meaning the code running on the client's computer communicates over the Internet with code hosted elsewhere on server computers. This interoperability, which is invisible to the user, makes the application fast and executable without heavy installation requirements. Instructions on how to reproduce the calculations presented here are available under “Publications” at the Web site http://www.OpenSHA.org/ .

Seismological Research Letters↗
USGS2016

Dense lower crust elevates long-term earthquake rates in the New Madrid seismic zone

Knowledge of the local state of stress is critical in appraising intraplate seismic hazard. Inverting earthquake moment tensors, we demonstrate that principal stress directions in the New Madrid seismic zone (NMSZ) differ significantly from those in the surrounding region. Faults in the NMSZ that are incompatible with slip in the regional stress field are favorably oriented relative to local stress. We jointly analyze seismic velocity, gravity, and topography to develop a 3-D crustal and upper mantle density model, revealing uniquely dense lower crust beneath the NMSZ. Finite element simulations then estimate the stress tensor due to gravitational body forces, which sums with regional stress. The anomalous lower crust both elevates gravity-derived stress at seismogenic depths in the NMSZ and rotates it to interfere more constructively with far-field stress, producing a regionally maximal deviatoric stress coincident with the highest concentration of modern seismicity. Moreover, predicted principal stress directions mirror variations (observed independently in moment tensors) at the NMSZ and across the region.

Arkansas, Missouri, Tennessee↗
USGS2006

Slip on the San Andreas fault at Parkfield, California, over two earthquake cycles, and the implications for seismic hazard

Parkfield, California, which experienced M 6.0 earthquakes in 1934, 1966, and 2004, is one of the few locales for which geodetic observations span multiple earthquake cycles. We undertake a comprehensive study of deformation over the most recent earthquake cycle and explore the results in the context of geodetic data collected prior to the 1966 event. Through joint inversion of the variety of Parkfield geodetic measurements (trilateration, two-color laser, and Global Positioning System), including previously unpublished two-color data, we estimate the spatial distribution of slip and slip rate along the San Andreas using a fault geometry based on precisely relocated seismicity. Although the three most recent Parkfield earthquakes appear complementary in their along-strike distributions of slip, they do not produce uniform strain release along strike over multiple seismic cycles. Since the 1934 earthquake, more than 1 m of slip deficit has accumulated on portions of the fault that slipped in the 1966 and 2004 earthquakes, and an average of 2 m of slip deficit exists on the 33 km of the fault southeast of Gold Hill to be released in a future, perhaps larger, earthquake. It appears that the fault is capable of partially releasing stored strain in moderate earthquakes, maintaining a disequilibrium through multiple earthquake cycles. This complicates the application of simple earthquake recurrence models that assume only the strain accumulated since the most recent event is relevant to the size or timing of an upcoming earthquake. Our findings further emphasize that accumulated slip deficit is not sufficient for earthquake nucleation.

Bulletin of the Seismological Society of America↗
USGS2016

Seismic imaging beneath an InSAR anomaly in eastern Washington State: Shallow faulting associated with an earthquake swarm in a low-hazard area

In 2001, a rare swarm of small, shallow earthquakes beneath the city of Spokane, Washington, caused ground shaking as well as audible booms over a five‐month period. Subsequent Interferometric Synthetic Aperture Radar (InSAR) data analysis revealed an area of surface uplift in the vicinity of the earthquake swarm. To investigate the potential faults that may have caused both the earthquakes and the topographic uplift, we collected ∼3 km of high‐resolution seismic‐reflection profiles to image the upper‐source region of the swarm. The two profiles reveal a complex deformational pattern within Quaternary alluvial, fluvial, and flood deposits, underlain by Tertiary basalts and basin sediments. At least 100 m of arching on a basalt surface in the upper 500 m is interpreted from both the seismic profiles and magnetic modeling. Two west‐dipping faults deform Quaternary sediments and project to the surface near the location of the Spokane fault defined from modeling of the InSAR data.

Washington↗