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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22Linked to original sources

USGS2021

Shear-wave velocity site characterization in Oklahoma from joint inversion of multi-method surface seismic measurements: Implications for central U.S. Ground Motion Prediction

We analyze multimethod shear (SH)‐wave velocity ( ⁠ V S "> V S ) site characterization data acquired at three permanent and 25 temporary seismograph stations in Oklahoma that recorded M 4+ earthquakes within a 50 km hypocentral distance of at least one of the 2016 M 5.1 Fairview, M 5.8 Pawnee, or M 5.0 Cushing earthquakes to better constrain earthquake ground‐motion modeling in the region. We acquired active‐source seismic data for time‐averaged V S "> V S to 30 m depth ( V S "> V S V S 30 "> 30 ) at 28 seismograph stations near the Fairview, Pawnee, and Cushing epicentral areas. The SH‐wave refraction travel times coupled with Rayleigh‐ and Love‐wave phase velocity dispersion were extracted and modeled in a nonlinear least‐squares (L2) joint inversion to obtain a best‐fit 1D V S "> V S versus depth profile for each site. At a subset of sites where the preferred L2 inverse model did not optimally fit each of the Love, Rayleigh, and SH travel‐time datasets, we explore application of simulated annealing in a joint inversion to find a more global solution. V S "> V S V S 30 "> 30 values range from 262 to 807    m / s "> 807 m / s for the preferred measured (in situ) V S "> V S profiles, or National Earthquake Hazards Reduction Program (NEHRP) site class D to B, and are broadly comparable with estimates from previous data reports in the region. Site amplification estimates were calculated next from 1D SH transfer functions of the preferred V S "> V S VS profiles and then compared against observed horizontal‐to‐vertical spectral ratios (HVSRs) from nearby seismograph stations. We generally see good agreement between the predicted in situ model and the observed HVSR resonant frequencies, with nominal amplifications between 2 and 10 within the 2–15 Hz frequency band. Next, using 40 known in situ V S "> V S V S 30 "> 30 measurements in the region, we demonstrate that the in situ V S "> V S V S 30 "> 30 values improve the fit for selected suites of ground‐motion models (GMMs) for M 4+ earthquakes within a 50 km hypocentral distance when compared with proxy methods, arguing for future development of GMMs implementing in situ V S "> V S profiles.

Oklahoma↗

USGS1988

Measuring contemporary crustal motions; NASA’s Crustal Dynamics Project

In the early and middle 1970's, two new space-based geodetic techniques became available that offered unprecedented accuracy in the measurement of distances over long baselines. As described below, Satellite Laser Ranging (SLR) and Very Long Baseline Interferometry (VLBI) provided a capability to determine in relatively short periods of time the inter-site distance between two observing stations to a level of several centimeters, even if they were separated by thousands or tens of thousands of kilometers. During the 1980's the two techniques have evolved to the point where baselines can now be routinely measured to a level below one centimeter. This is a tenfold improvement in about ten years. Perhaps more important than providing the measurement capability are the consequences of this capability when applied over a few years: the motion of any one site with respect to another can be monitored at the level of better than a cm/yr. This has made possible, for the first time, direct measurement of the motion of the Earth's tectonic plates and of the deformation of its crust in active plate boundary regions. Because such measurements have the potential for providing important information on the mechanisms that drive plates and cause them to slip during earthquakes, a coordinated federal program for the application of this space technology to crustal dynamics and earthquake research was established in 1979. The participating agencies were the National Aeronautics and Space Administration (NASA), the National Oceanic and Atmospheric Administration (NOAA), the United States Geological Survey (USGS), the National Science Foundation (NSF), and the Defense Mapping Agency. NASA formed the Crustal Dynamics Project to further develop the SLR and VLBI techniques, to implement global networks of stations with the cooperation of many different countries, and to conduct measurements of plate motions and regional deformations. this project has grown into an international effort to collect information directly relevant to understanding the threat of the earthquake hazard as well as fundamental research about the crust and upper mantle. In this article we describe briefly the two space geodetic techniques and how they are used by the Crustal Dynamics Project, show some of the very exciting results that have emerged at the halfway point in the project's life, describe the availability and utilization of the data being collected, and consider what the future may hold when measurement accuracies eventually exceed even those now available and when other international groups become more heavily involved.

Earthquakes & Volcanoes (USGS)↗

USGS2014

Slip rates and spatially variable creep on faults of the northern San Andreas system inferred through Bayesian inversion of Global Positioning System data

Fault creep, depending on its rate and spatial extent, is thought to reduce earthquake hazard by releasing tectonic strain aseismically. We use Bayesian inversion and a newly expanded GPS data set to infer the deep slip rates below assigned locking depths on the San Andreas, Maacama, and Bartlett Springs Faults of Northern California and, for the latter two, the spatially variable interseismic creep rate above the locking depth. We estimate deep slip rates of 21.5 ± 0.5, 13.1 ± 0.8, and 7.5 ± 0.7 mm/yr below 16 km, 9 km, and 13 km on the San Andreas, Maacama, and Bartlett Springs Faults, respectively. We infer that on average the Bartlett Springs fault creeps from the Earth's surface to 13 km depth, and below 5 km the creep rate approaches the deep slip rate. This implies that microseismicity may extend below the locking depth; however, we cannot rule out the presence of locked patches in the seismogenic zone that could generate moderate earthquakes. Our estimated Maacama creep rate, while comparable to the inferred deep slip rate at the Earth's surface, decreases with depth, implying a slip deficit exists. The Maacama deep slip rate estimate, 13.1 mm/yr, exceeds long-term geologic slip rate estimates, perhaps due to distributed off-fault strain or the presence of multiple active fault strands. While our creep rate estimates are relatively insensitive to choice of model locking depth, insufficient independent information regarding locking depths is a source of epistemic uncertainty that impacts deep slip rate estimates.

California↗

USGS1997

Future of the US National Strong-Motion Program

These reports are presented in response to a charge of the Earthquake Hazards Reduction Program (EHRP) Program Council of the U.S. Geological Survey to 'define the future of the USGS National Strong-Motion Program (NSMP)' (Appendix A). The council requested that a 'Vision Paper' and an 'Options Document' be prepared. Each of these reports is a separate document. The 'Executive Summary' of both reports is presented here under separate cover. The content of the reports and the executive summary is organized to the extent possible according to the main charges of the EHRP Program Council. Detailed evaluation of charges and questions posed by the Council with the benefit of modern GIS tools and the recent deliberations of three national workshops have yielded substantial new insight regarding the status of strong-motion measurement in the US. These reports are intended as an objective evaluation of the nation's need and capability to record the next major earthquake at locations of most significance for public safety.

Open-File Report↗

USGS1997

Executive summary of vision and options for the future of the US National Strong-Motion Program

Open-File Report↗

USGS2014

Quantitative study of tectonic geomorphology along Haiyuan fault based on airborne LiDAR

High-precision and high-resolution topography are the fundamental data for active fault research. Light detection and ranging (LiDAR) presents a new approach to build detailed digital elevation models effectively. We take the Haiyuan fault in Gansu Province as an example of how LiDAR data may be used to improve the study of active faults and the risk assessment of related hazards. In the eastern segment of the Haiyuan fault, the Shaomayin site has been comprehensively investigated in previous research because of its exemplary tectonic topographic features. Based on unprecedented LiDAR data, the horizontal and vertical coseismic offsets at the Shaomayin site are described. The measured horizontal value is about 8.6 m, and the vertical value is about 0.8 m. Using prior dating ages sampled from the same location, we estimate the horizontal slip rate as 4.0 ± 1.0 mm/a with high confidence and define that the lower bound of the vertical slip rate is 0.4 ± 0.1 mm/a since the Holocene. LiDAR data can repeat the measurements of field work on quantifying offsets of tectonic landform features quite well. The offset landforms are visualized on an office computer workstation easily, and specialized software may be used to obtain displacement quantitatively. By combining precious chronological results, the fundamental link between fault activity and large earthquakes is better recognized, as well as the potential risk for future earthquake hazards.

Gansu↗

USGS1988

What is worse than the “big one”?

The Whittier Narrows California earthquake sequence (local magnitude, Ml=5.9 or 1 October, 1987), which caused over $358 million damage, indicates that assessments of earthquake hazards in Los Angeles metropolitan area may be underestimated. the sequence ruptured a previously unidentified thrust fault that may be part of a large system of thrust faults that extends across the entire east-west length of the northern margin of the Los Angeles basin. Peak horizontal accelerations from the main shock, which were measured at ground level and in structures, were as high as 0.6g (where g is acceleration of gravity at sea level) within 50 kilometers of the epicenter The first thought in the minds of many residents of the city of Whittier when the first shock hit them was "Is this the big one?" the San Andreas' once-in-150-years great shaker? It might as well have been for Whittier, which is 20 kilometers east of downtown Los Angeles. The ground shook harder there this month than it will when the big one does strike the distant San Andreas, which lies 50 kilometers on the other side of the mountains. And this was only a moderate, magnitude 6.1 shock. Earthquake of magnitude 7 and large 30 times more powerful, could rupture faults beneath the feet of Angelenos at any time. The loss of life and destruction could exceed that caused by the big one.

California↗

USGS2015

Holocene geologic slip rate for the Banning strand of the southern San Andreas Fault, southern California

Northwest directed slip from the southern San Andreas Fault is transferred to the Mission Creek, Banning, and Garnet Hill fault strands in the northwestern Coachella Valley. How slip is partitioned between these three faults is critical to southern California seismic hazard estimates but is poorly understood. In this paper, we report the first slip rate measured for the Banning fault strand. We constrain the depositional age of an alluvial fan offset 25 ± 5 m from its source by the Banning strand to between 5.1 ± 0.4 ka (95% confidence interval (CI)) and 6.4 + 3.7/−2.1 ka (95% CI) using U-series dating of pedogenic carbonate clast coatings and 10 Be cosmogenic nuclide exposure dating of surface clasts. We calculate a Holocene geologic slip rate for the Banning strand of 3.9 + 2.3/−1.6 mm/yr (median, 95% CI) to 4.9 + 1.0/−0.9 mm/yr (median, 95% CI). This rate represents only 25–35% of the total slip accommodated by this section of the southern San Andreas Fault, suggesting a model in which slip is less concentrated on the Banning strand than previously thought. In rejecting the possibility that the Banning strand is the dominant structure, our results highlight an even greater need for slip rate and paleoseismic measurements along faults in the northwestern Coachella Valley in order to test the validity of current earthquake hazard models. In addition, our comparison of ages measured with U-series and 10 Be exposure dating demonstrates the importance of using multiple geochronometers when estimating the depositional age of alluvial landforms.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

Ground motion in Anchorage, Alaska, from the 2002 Denali fault earthquake: Site response and displacement pulses

Data from the 2002 Denali fault earthquake recorded at 26 sites in and near Anchorage, Alaska, show a number of systematic features important in studies of site response and in constructing long-period spectra for use in earthquake engineering. The data demonstrate that National Earthquake Hazards Reduction Program (NEHRP) site classes are a useful way of grouping stations according to site amplification. In general, the sites underlain by lower shear-wave velocities have higher amplification. The amplification on NEHRP class D sites exceeds a factor of 2 relative to an average of motions on class C sites. The amplifications are period dependent. They are in rough agreement with those from previous studies, but the new data show that the amplifications extend to at least 10 sec, periods longer than considered in previous studies. At periods longer than about 14 sec, all sites have motion of similar amplitude, and the ground displacements are similar in shape, polarization, and amplitude for all stations. The displacement ground motion is dominated by a series of four pulses, which are associated with the three subevents identified in inversion studies (the first pulse is composed of P waves from the first subevent). Most of the high-frequency ground motion is associated with the S waves from subevent 1. The pulses from subevents 1 and 2, with moment releases corresponding to M 7.1 and 7.0, are similar to the pulse of displacement radiated by the M 7.1 Hector Mine earthquake. The signature from the largest subevent (M 7.6) is more subdued than those from the first two subevents. The two largest pulses produce response spectra with peaks at a period of about 15 sec. The spectral shape at long periods is in good agreement with the recent 2003 NEHRP code spectra but is in poor agreement with the shape obtained from Eurocode 8.

Alaska↗

USGS2006

High-resolution boomer seismic-reflection profiles of the shelf off southern California from cruise A-1-00-SC: Santa Monica Bay to San Diego

High-resolution boomer data were collected in the California Continental Borderland as part of the southern California Earthquake Hazards Task of the Southern California Coastal and Marine Geology Regional Investigations Project. During the period from 1997 to 2002, five data-acquisition cruises collected seismic-reflection data using several different systems from offshore Santa Barbara, California, south to the Exclusive Economic Zone boundary with Mexico. A key mission of this project was to map late Quaternary deformation in addition to improving our understanding of which offshore fault zones might have potential to damage highly populated areas of southern California. State regulations concerning the use of seismic-reflection equipment within three miles of the coastline precluded the routine gathering of high-resolution multichannel data in that swath adjacent to the coast. Boomer seismic-reflection data, however, can be obtained within the state 3-mile limit provided the operation receives authorization from the California State Lands Commission. The Geopulse boomer data accessible through this report were collected on the cruise A-1-00-SC, which was the only survey where we requested permission to work inside the 3-mile limit of the State of California. These data are critical to discovering connections between onshore and offshore faults, the overall lengths of which are related to the potential size of an earthquake that might be generated along them. The 2000 survey was designed to fill the gap between onshore data and reflection data obtained in deeper water on previous cruises as well as data anticipated from future surveys. This report includes trackline maps showing the location of the data, as well as both digital data files (SEG-Y) and images of all of the profiles.

California↗

USGS2021

Late Holocene slip rate of the Mojave section of the San Andreas Fault near Palmdale, California

The geologic slip rate on the Mojave section of the San Andreas fault is poorly constrained, despite its importance for understanding earthquake hazard, apparent discrepancies between geologic and geodetic slip rates along this fault section, and long‐term fault interactions in southern California. Here, we use surficial geologic mapping, excavations, and radiocarbon and luminescence dating to quantify the displacements and ages of late Holocene landforms offset by the fault at three sites. At the Ranch Center site, the slip rate is determined using the base of a fan marking incision and deflection of an ephemeral channel. At the adjacent Key Slide site, the margin of a landslide deposited on indigenous fire hearths provides a minimum rate. At the X‐12 site, the slip rate is determined from a channel that incised into a broad fan surface, and is deflected and beheaded by the fault. We use maximum–minimum bounds on both the displacement and age of each offset feature to calculate slip rate for each site independently. Overlap of the three independent rate ranges yields a rate of 33–39 mm/yr over the last 3 ka, under the assumption that the sites share a common history, given their proximity. Considered in sequence, site‐level epistemic uncertainties in the data permit but do not require a rate increase since ∼1200 cal B.P. Modest rate changes can be explained by aleatory variability in earthquake timing and magnitude; larger changes could suggest a shared regional variation with the Garlock and other faults. The new late Holocene slip rates are consistent with geodetic model estimates that include a viscoelastic crust and earthquake cycle effects. The geologic slip rates also provide average slip over dozens of earthquake cycles—a key constraint for long‐term earthquake rupture forecasts.

California↗

USGS1989

A satellite-based digital data system for low-frequency geophysical data

A reliable method for collection, display, and analysis of low-frequency geophysical data from isolated sites, which can be throughout North and South America and the Pacific Rim, has been developed for use with the Geostationary Operational Environmental Satellite (GOES) system. Geophysical data primarily intended for earthquake hazard and crustal deformation monitoring are digitized with either 12-bit or 16-bit resolution and transmitted every 10 min through a satellite link to a bank of UNIX-based computers in Menlo Park, California. There the data are available for analysis and display within a few seconds of their transmit time. This system provides real-time monitoring of crustal deformation parameters such as tilt, strain, fault displacement, local magnetic field, crustal geochemistry, and water levels, as well as meteorological and other parameters, along faults in California and Alaska, and in volcanic regions in the western United States, Rabaul, and other locations in the New Britain region of the South Pacific. Various mathematical, statistical, and graphical algorithms process the incoming data to detect changes in crustal deformation and fault slip that may indicate the first stages of catastrophic fault failure. Alert trigger levels based on physical models, signal resolution, and previous history have been defined for particular instrument types. Computer-driven remote paging and mail systems are used to notify appropriate personnel when alarm status is reached. The system supports continuous historical records of low-frequency geophysical data, software for extensive analysis of these data, and programs for modeling fault rupture with and without seismic radiation, as well as providing an environment for real-time attempts at earthquake prediction.

Bulletin of the Seismological Society of America↗

USGS2020

A global hybrid VS30 map with a topographic slope–based default and regional map insets

Time-averaged shear wave velocity over the upper 30 m of the earth’s surface ( V S 30 ) is a key parameter for estimating ground motion amplification as both a predictive and a diagnostic tool for earthquake hazards. The first-order approximation of V S 30 is commonly obtained through a topographic slope–based or terrain proxy due to the widely available nature of digital elevation models. However, better-constrained V S 30 maps have been developed in many regions. Such maps preferentially employ various combinations of V S 30 measurements, higher-resolution elevation models, lithologic, geologic, geomorphic, and other proxies and often utilize refined interpolation schemes. We develop a new hybrid global V S 30 map database that defaults to the global slope-based V S 30 map, but smoothly inserts regional V S 30 maps where available. In addition, we present comparisons of the default slope-based proxy maps against the new hybrid version in terms of V S 30 and amplification ratio maps, and uncertainties in assigned V S 30 values.

Earthquake Spectra↗

USGS2026

3D semantic mapping of surface geological features

Semantic mapping in 3D is fundamental to a wide range of geoscientific studies and applications, including geomorphology, hazard assessment, and environmental monitoring. However, automatically segmenting geological features from large-scale photogrammetric datasets remains a significant challenge. We present a methodology to address this gap. Using overlapping images collected over environments of interest, Structure-from-Motion (SfM) produces georeferenced point clouds and estimates camera poses. Existing large vision models, such as Segment Anything Model, segment objects in the images, generating pixel-segmentation associations. To produce pixel-point associations, we project the points back onto the camera image planes. As objects are independently segmented across multiple images with different perspectives, we develop a segmentation mosaicking algorithm to build probabilistic point-segmentation associations that combines the pixel-segmentation associations and pixel-point associations. Our methodology is validated using both synthetic data generated by Kubric and real-world UAV-SfM data. The implementation is designed to be compatible with existing SfM software, including Agisoft and OpenDroneMap, for photogrammetry mapping in geoscience studies. As a case study, we apply our method to the semantic mapping of precariously balanced rocks (PBRs), which provide upper-bound constraints on historical ground motion shaking intensity. To support object-level identification of PBRs, we additionally integrated Grounding DINO, enabling text-prompted segmentation of features of interest within UAV imagery. This case study demonstrates the effectiveness of our method in generating a 3D semantic map of PBRs, enabling spatial distribution of PBR fragility for earthquake hazard analysis.

Computers & Geosciences↗

USGS2018

The 2015 Gorkha (Nepal) Earthquake sequence: I. Source modeling and deterministic 3D ground shaking

To better quantify the relatively long period (< 0.3 Hz) shaking experienced during the 2015 Gorkha (Nepal) earthquake sequence, we study the finite rupture processes and the associated 3D ground motion of the Mw7.8 mainshock and the Mw7.2 aftershock . The 3D synthetics are then used in the broadband ground shaking in Kathmandu with a hybrid approach, summarized in a companion paper (Chen and Wei, 2017, submitted together). We determined the coseismic rupture process of the mainshock by joint inversion of InSAR/SAR, GPS (static and high-rate), strong motion and teleseismic waveforms. Our inversion for the mainshock indicates unilateral rupture towards the ESE, with an average rupture speed of 3.0 km/s and a total duration of ~ 60 s. Additionally, we find that the beginning part of the rupture (5–18 s) has about 40% longer rise time than the rest of the rupture, as well as slower rupture velocity. Our model shows two strong asperities occurring ~ 24 s and ~ 36 s after the origin and located ~ 30 km to the northwest and northeast of the Kathmandu valley, respectively. In contrast, the Mw7.2 aftershock is more compact both in time and space, as revealed by joint inversion of teleseismic body waves and InSAR data. The different rupture features between the mainshock and the aftershock could be related to difference in fault zone structure. The mainshock and aftershock ground motions in the Kathmandu valley, recorded by both strong motion and high-rate GPS stations, exhibited strong amplification around 0.2 Hz. A simplified 3D basin model, calibrated by an Mw5.2 aftershock, can match the observed waveforms reasonably well at 0.3 Hz and lower frequency. The 3D simulations indicate that the basin structure trapped the wavefield and produced an extensive ground vibration. Our study suggests that the combination of rupture characteristics and propagational complexity are required to understand the ground shaking produced by hazardous earthquakes such as the Gorkha event.

Tectonophysics↗

USGS2015

Chance findings about early holocene tidal marshes of Grays Harbor, Washington, in relation to rapidly rising seas and great subduction earthquakes

Tidal marshes commonly build upward apace with gradual rise in the level of the sea. It is expected, however, that few tidal marshes will keep up with accelerated sea-level rise later in this century. Tidal marshes have been drowned, moreover, after subsiding during earthquakes. This report tells of ancient marshes that endured rapid sea-level rise in a region that subsides during earthquakes. The soils of these marshes were unexpectedly encountered in borings for a public-works project at Grays Harbor, Washington. The borings were logged quickly and only a few of the core sections were conserved. The limited findings pose puzzles about how the ancient marshes endured and what their history implies for earthquake hazards. The borings establish that tidal marshes persisted during the early Holocene at Grays Harbor, an estuary along the Cascadia Subduction Zone of western North America. The persistent marshes are recorded by a unit of peaty mud up to 10 m thick and as much as 40 m below present sea level in the drowned valley of the Chehalis River. The unit was encountered in two areas 4 km apart that were tidal flats in the 19th century. The marshes originated less than 10,000 years ago and endured through most or all of an estimated 500–1,500 years. The borings further show that these persistent marshes eventually yielded to tidal flats, tidal channels, or both. The change is marked by sand and mud that overlie the peaty mud at a typically sharp contact. The marshes were drowned about 8,600–8,400 years ago if the sand and mud buried them promptly, or later if the sand and mud filled channels that migrated across the peaty mud. In one of the studied areas, tidal marshes became re-established locally in the early Holocene and widely in the middle Holocene, and deposits of middle Holocene marshes were overrun as recently as 1,000 years ago by a gravelly tidal channel. In the other area, tidal-flat and probably subtidal deposits make up all of the middle and late Holocene section below artificial fill; if marshes became re-established in this area after about 8,600–8,400 years ago, their deposits have been lost to erosion. The puzzles posed by these findings include: (1) How did the marshes manage to endure centuries of relative sea-level rise that likely approached 1 cm/yr on average? (2) Did the marshes also endure subsidence that accompanied great thrust earthquakes on the Cascadia Subduction Zone? (3) Was their eventual drowning triggered by a Cascadia earthquake of unusually large size, or can the drowning be explained by sea-level rise that included a jump from drainage of glacial Lake Agassiz?

Washington↗

USGS2004

Prediction of nonlinear soil effects

Mathematical models of soil nonlinearity in common use and recently developed nonlinear codes compared to investigate the range of their predictions. We consider equivalent linear formulations with and without frequency-dependent moduli and damping ratios and nonlinear formulations for total and effective stress. Average velocity profiles to 150 m depth with midrange National Earthquake Hazards Reduction Program site classifications (B, BC, C, D, and E) in the top 30 m are used to compare the response of a wide range of site conditions from rock to soft soil. Nonlinear soil models are compared using the amplification spectrum, calculated as the ratio of surface ground motion to the input motion at the base of the velocity profile. Peak input motions from 0.1 g to 0.9 g are considered. For site class B, no significant differences exist between the models considered in this article. For site classes BC and C, differences are small at low input motions (0.1 g to 0.2 g ), but become significant at higher input levels. For site classes D and E the overdamping of frequencies above about 4 Hz by the equivalent linear solution with frequency-independent parameters is apparent for the entire range of input motions considered. The equivalent linear formulation with frequency-dependent moduli and damping ratios under damps relative to the nonlinear models considered for site class C with larger input motions and most input levels for site classes D and E. At larger input motions the underdamping for site classes D and E is not as severe as the overdamping with the frequency-independent formulation, but there are still significant differences in the time domain. A nonlinear formulation is recommended for site classes D and E and for site classes BC and C with input motions greater than a few tenths of the acceleration of gravity. The type of nonlinear formulation to use is driven by considerations of the importance of water content and the availability of laboratory soils data. Our average amplification curves from a nonlinear effective stress formulation compare favorably with observed spectral amplification at class D and E sites in the Seattle area for the 2001 Nisqually earthquake.

Bulletin of the Seismological Society of America↗

USGS1997

Crustal structure along the west flank of the Cascades, western Washington

Knowledge of the crustal structure of the Washington Cascades and adjacent Puget Lowland is important to both earthquake hazards studies and geologic studies of the evolution of this tectonically active region. We present a model for crustal velocity structure derived from analysis of seismic refraction/wide-angle reflection data collected in 1991 in western Washington. The 280-km-long north-south transect skirts the west flank of the Cascades as it crosses three tectonic provinces including the Northwest Cascades Thrust System (NWCS), the Puget Lowland, and the volcanic arc of the southern Cascades. Within the NWCS, upper crustal velocities range from 4.2 to 5.7 km s-1 and are consistent with the presence of a diverse suite of Mesozoic and Paleozoic metasediments and metavolcanics. In the upper 2-3 km of the Puget Lowland velocities drop to 1.7-3.5 km s-1 and reflect the occurrence of Oligocene to recent sediments within the basin. In the southern Washington Cascades, upper crustal velocities range from 4.0 to 5.5 km s-1 and are consistent with a large volume of Tertiary sediments and volcanics. A sharp change in velocity gradient at 5-10 km marks the division between the upper and middle crust. From approximately 10 to 35 km depth the velocity field is characterized by a velocity increase from ???6.0 to 7.2 km s-1. These high velocities do not support the presence of marine sedimentary rocks at depths of 10-20 km beneath the Cascades as previously proposed on the basis of magnetotelluric data. Crustal thickness ranges from 42 to 47 km along the profile. The lowermost crust consists of a 2 to 8-km-thick transitional layer with velocities of 7.3-7.4 km s-1. The upper mantle velocity appears to be an unusually low 7.6-7.8 km s-1. When compared to velocity models from other regions, this model most closely resembles those found in active continental arcs. Distinct seismicity patterns can be associated with individual tectonic provinces along the seismic transect. In the NWCS and Puget Lowland, most of the seismicity occurs below the base of the upper crust as defined by a seismic boundary at 5-10 km depth and continues to 20-30 km depth. The region of transition between the NWCS and the Puget Lowland appears as a gap in seismicity with notably less seismic activity north of the boundary between the two. Earthquakes within the Cascades are generally shallower (0-20 km) and are dominated by events associated with the Rainier Seismic Zone.

Washington↗