Geology Reports⌕ Search

Geology topics

Research about British Columbia, California, Oregon, Washington

Source-linked reports with geographic coverage including British Columbia, California, Oregon, Washington.

At least 19 recordsLinked to original sources
USGS2026

3D broadband ground motion simulations of Cascadia megathrust earthquakes: Full-margin scenarios incorporating variable magnitudes and frontal thrust splay faulting

Great megathrust earthquakes (∼ M w 8–9) that rupture along the Cascadia Subduction Zone pose significant shaking and tsunami hazards to the U.S. Pacific Northwest. Although these large earthquakes are documented in the geologic record, none have been recorded with modern instrumentation, complicating regional seismic hazard assessments. To improve constraints on seismic hazard estimates for future earthquakes, we develop 3D broadband (0–10 Hz) ground‐motion simulations for 36 full‐margin megathrust earthquake scenarios through the Cascadia Coastlines and Peoples Hazards Research Hub. These simulations expand on the M9 Project ( Frankel et al. , 2018 ; Wirth, Frankel, et al. , 2018 ) by incorporating a broader range of rupture parameters—variable magnitudes, splay faulting, and geodetically informed slip distributions. Benchmarking these simulations against the Next Generation Attenuation Subduction (NGA‐Sub) ground‐motion models shows reasonable fits for non‐basin sites, although simulations predict stronger shaking at short (0–2 s) and long (6–10 s) periods. NGA‐Sub basin terms improve fits to synthetic shaking within the Seattle Basin but underpredict amplitudes by a factor of 2 at periods of 2–5 s. We also evaluate how rupture characteristics influence regional inland shaking variability and basin response in the Puget Lowlands. Larger magnitude scenarios generally cause stronger inland shaking due to farther inland extents of slip and larger strong‐motion‐generating fault patches but show little difference in basin amplification within the Puget Lowlands compared to smaller magnitude events. Along‐dip hypocenter location strongly influences overall inland shaking, with offshore (updip) hypocenters generating on average ∼45% higher shaking at long periods (>6 s) compared to onshore (downdip) hypocenters due to rupture directivity effects. Offshore splay faults minimally affect inland shaking but alter offshore deformation important for coupled ground‐shaking and tsunami scenarios. Finally, we show that rupture directivity increases amplification within the Seattle Basin through the generation and interaction of higher order mode surface waves, an effect of 3D basin geometry that cannot be captured using 1D approximations.

British Columbia, California, Oregon, Washington↗
USGS2025

Three-dimensional seismic velocity model for the Cascadia Subduction Zone with shallow soils and topography, version 1.7

The U.S. Geological Survey’s seismic velocity model for the Cascadia Subduction Zone provides P- and S-wave velocity ( V P and V S , respectively) information from 40.2° to 50.0° N. latitude and −129.0° to −121.0° W. longitude, and is used to support a variety of research topics, including three-dimensional (3D) earthquake simulations and seismic hazard assessment in the Pacific Northwest. This report describes an update to the previous version (v) 1.6 of the 3D seismic velocity model for the Cascadia Subduction Zone. This new model (herein referred to as v1.7) contains more detailed near-surface structure for improved earthquake ground motion modeling. Updated features include the addition of a new shallow soil velocity model in the top few hundred meters and the option of adding user-specified topography. Although v1.6 of the Cascadia seismic velocity model has a minimum V S of 600 meters per second (m/s), the new model (v1.7) has a minimum V S of approximately 40 m/s. Overall, this update will allow for more accurate ground motion estimates from 3D simulations of scenario earthquakes in the Cascadia Subduction Zone region.

British Columbia, California, Oregon, Washington↗
USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗
USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2024

Late-Quaternary surface displacements on accretionary wedge splay faults in the Cascadia Subduction Zone: Implications for megathrust rupture

Because splay faults branch at a steep dip angle from the plate-boundary décollement in an accretionary wedge, their coseismic displacement can potentially result in larger tsunamis with distinct characteristics compared to megathrust-only fault ruptures, posing an enhanced hazard to coastal communities. Elsewhere, there is evidence of coseismic slip on splay faults during many of the largest subduction zone earthquakes, but our understanding of potentially active splay faults and their hazards at the Cascadia subduction zone remains limited. To identify the most recently active splay faults at Cascadia, we conduct stratigraphic and structural interpretations of near-surface deformation in the outer accretionary wedge for the ~400 km along-strike length of the landward vergence zone. We analyze recently acquired high-frequency sparker seismic data and crustal-scale multi-channel seismic data to examine the record of deformation in shallow slope basins and the upper ~1 km of the surrounding accreted sediments and to investigate linkages to deeper décollement structure. We present a new fault map for widest, most completely locked portion of Cascadia from 45 to 48°N latitude, which documents the distribution of faults that show clear evidence of recent late Quaternary activity. We find widespread evidence for active splay faulting up to 30 km landward of the deformation front, in what we define as the active domain, and diminished fault activity landward outside of this zone. The abundance of surface-deforming splay faults in the active outer wedge domain suggests Cascadia megathrust events may commonly host distributed shallow rupture on multiple splay faults located within 30 km of the deformation front.

British Columbia, California, Oregon, Washington↗
USGS2022

Fluid transport and storage in the Cascadia forearc influenced by overriding plate lithology

Subduction of hydrated oceanic lithosphere can carry water deep into the Earth, with consequences for a range of tectonic and magmatic processes. Most of the fluid is released in the forearc where it plays a critical role in controlling the mechanical properties and seismic behaviour of the subduction megathrust. Here we present results from three-dimensional inversions of data from nearly 400 long-period magnetotelluric sites, including 64 offshore, to provide insights into the distribution of fluids in the forearc of the Cascadia subduction zone. We constrain the geometry of the electrically resistive Siletz terrane, a thickened section of oceanic crust accreted to North America in the Eocene, and the conductive accretionary complex underthrust along the margin. We find that fluids accumulate over timescales exceeding 1 My above the plate in metasedimentary units, while the mafic rocks of Siletzia remain dry. Fluid concentrations tend to peak at slab depths of 17.5 and 30 km, suggesting control by metamorphic processes, but also concentrate around the edges of Siletzia, suggesting that this mafic block is impermeable, with dehydration fluids escaping up-dip along the megathrust. Our results demonstrate that the lithology of the overriding crust can play a critical role in controlling fluid transport in a subduction zone.

British Columbia, California, Oregon, Washington↗
USGS2022

The nature of active magma reservoirs and storage underneath Cascade volcanoes

Volcanic activity is dictated by crustal pathways and reservoirs through which magma ascends and collects. The Cascades are a natural laboratory to study the variable influence of the crust on erupted magmas and subterranean reservoirs. The interaction of ascending magmas with the variable subduction geometry and tectonics of the overriding North American plate have given rise to a diversity of magmatic storage conditions. The confluence of geochemical and geophysical investigations emphasizes that most magmatic systems in the Cascades have been built at multiple levels in the crust as determined by tectonics, pre-existing structure, and magmatic flux from the mantle.

British Columbia, California, Oregon, Washington↗
USGS2021

Imaging the next Cascadia earthquake: Optimal design for a seafloor GNSS- A network

The Cascadia subduction zone in the Pacific Northwest of the United States of America capable of producing magnitude ∼9 earthquakes, likely often accompanied by tsunamis. An outstanding question in this region is the degree and spatial extent of interseismic strain accumulation on the subduction megathrust. Seafloor geodetic methods combining GNSS and underwater acoustic ranging (GNSS-A) are capable of imaging this strain accumulation on the offshore portion of the subduction zone and therefore anticipating the potential size and rupture pattern of a future earthquake. However, the high cost of seafloor geodesy means that only a limited number of stations may be deployed and monitored. To facilitate expansion of current geodetic networks offshore, we develop a quantitative recommendation of optimal locations for future seafloor geodetic observations, based on the amount of new information provided by that observation. The optimal network depends on the problem that one is trying to solve with those observations (mapping subduction locking rates, coupling rates, constraining total moment rate, etc.), and on a number of modelling and data uncertainty assumptions. In particular, data uncertainty assumptions will change over time, as more position observations reduce velocity uncertainties. We find that near-trench observations on the megathrust hangingwall, distributed along-strike, consistently provide significant reduction in differential entropy over a large suite of assumptions, and that a well-placed seafloor observation can provide up to ∼30 times the information gain of the most optimal onshore observation.

British Columbia, California, Oregon, Washington↗
USGS2021

The PLUM earthquake early warning algorithm: A retrospective case study of West Coast, USA, data

The PLUM (Propagation of Local Undamped Motion) earthquake early warning (EEW) algorithm differs from typical source-based EEW algorithms as it predicts shaking directly from observed shaking without first deriving earthquake source information (e.g., magnitude and epicenter). Here, we determine optimal PLUM event detection thresholds for U.S. West Coast earthquakes using two data sets: 558 M3.5+ earthquakes (California, Oregon, Washington; 2012–2017) and the ShakeAlert test suite of historic and problematic signals (1999–2015). PLUM computes Modified Mercalli Intensity ( I MMI ) using velocity and acceleration data, leveraging co-located sensors to avoid problematic signals. An event detection is issued when the observed I MMI exceeds a given threshold(s). We find a two-station detection method using I MMI trigger thresholds of 4.0 and 3.0 for the first and second stations, respectively, is optimal for detecting M4.5+ earthquakes. PLUM detected 79 events in the 2012–2017 data set, reporting (not including telemetry or alert dissemination) detection times on par, and sometimes faster than current EEW methods (mean 8 s; median 6 s). As expected, detection times were slower for the older 1999–2015 earthquakes ( N = 21; mean 11 s; median 6 s) when station coverage was sparser. Of the 31 PLUM detected M5+ events (10 2012–2017; 21 1999–2015), theoretically 20 (∼65%) could provide timely warnings. PLUM issued no false detections and avoided issuing detections for all calibration/anomalous signals, regional and teleseismic events. We conclude PLUM can successfully identify I MMI 4+ shaking from local earthquakes and could complement and enhance EEW in the U.S.

British Columbia, California, Oregon, Washington↗
USGS2021

The productivity of Cascadia aftershock sequences

This study addresses questions about the productivity of Cascadia mainshock–aftershock sequences using earthquake catalogs produced by the Geological Survey of Canada and the Pacific Northwest Seismic Network. Questions concern the likelihood that future moderate to large intermediate depth intraslab earthquakes in Cascadia would have as few detectable aftershocks as those documented since 1949. More broadly, for Cascadia, we consider if aftershock productivities vary spatially, if they are outliers among global subduction zones, and if they are consistent with a physical model in which aftershocks are clock‐advanced versions of tectonically driven background seismicity. A practical motivation for this study is to assess the likely accuracy of aftershock forecasts based on productivities derived from global data that are now being issued routinely by the U.S. Geological Survey. For this reason, we estimated productivity following the identical procedures used in those forecasts and described in Page et al. (2016) . Results indicate that in Cascadia we can say that the next intermediate depth intraslab earthquake will likely have just a few detectable aftershocks and that aftershock productivity appears to be an outlier among global subduction zones, with rates that on average are lower by more than half, except for mainshocks in the upper plate. Our results are consistent with a clock‐advance model; productivities may be related to the proximity of mainshocks to a population of seismogenic fault patches and correlate with background seismicity rates. The latter and a clear correlation between productivities with mainshock depth indicate that both factors may have predictive value for aftershock forecasting.

British Columbia, California, Oregon, Washington↗
USGS2020

Minimal stratigraphic evidence for coseismic coastal subsidence during 2000 yr of megathrust earthquakes at the central Cascadia subduction zone

Lithology and microfossil biostratigraphy beneath the marshes of a central Oregon estuary limit geophysical models of Cascadia megathrust rupture during successive earthquakes by ruling out >0.5 m of coseismic coastal subsidence for the past 2000 yr. Although the stratigraphy in cores and outcrops includes as many as 12 peat-mud contacts, like those commonly inferred to record subsidence during megathrust earthquakes, mapping, qualitative diatom analysis, foraminiferal transfer function analysis, and 14 C dating of the contacts failed to confirm that any contacts formed through subsidence during great earthquakes. Based on the youngest peat-mud contact’s distinctness, >400 m distribution, ∼0.6 m depth, and overlying probable tsunami deposit, we attribute it to the great 1700 CE Cascadia earthquake and(or) its accompanying tsunami. Minimal changes in diatom assemblages from below the contact to above its probable tsunami deposit suggest that the lower of several foraminiferal transfer function reconstructions of coseismic subsidence across the contact (0.1–0.5 m) is most accurate. The more limited stratigraphic extent and minimal changes in lithology, foraminifera, and(or) diatom assemblages across the other 11 peat-mud contacts are insufficient to distinguish them from contacts formed through small, gradual, or localized changes in tide levels during river floods, storm surges, and gradual sea-level rise. Although no data preclude any contacts from being synchronous with a megathrust earthquake, the evidence is equally consistent with all contacts recording relative sea-level changes below the ∼0.5 m detection threshold for distinguishing coseismic from nonseismic changes.

British Columbia, California, Oregon, Washington↗
USGS2020

Systematic characterization of morphotectonic variability along the Cascadia convergent margin: Implications for shallow megathrust behavior and tsunami hazards

Studies of recent destructive megathrust earth­quakes and tsunamis along subduction margins in Japan, Sumatra, and Chile have linked forearc mor­phology and structure to megathrust behavior. This connection is based on the idea that spatial varia­tions in the frictional behavior of the megathrust influence the tectono-morphological evolution of the upper plate. Here we present a comprehen­sive examination of the tectonic geomorphology, outer wedge taper, and structural vergence along the marine forearc of the Cascadia subduction zone (offshore northwestern North America). The goal is to better understand geologic controls on outer wedge strength and segmentation at spatial scales equivalent to rupture lengths of large earthquakes (≥M 6.7), and to examine potential linkages with shallow megathrust behavior. We use cross-margin profiles, spaced 25 km apart, to characterize along-strike variation in outer wedge width, steepness, and structural vergence (measured between the toe and the outer arc high). The width of the outer wedge varies between 17 and 93 km, and the steepness ranges from 0.9° to 6.5°. Hierarchical cluster analysis of outer wedge width and steepness reveals four distinct regions that also display unique patterns of structural ver­gence and shape of the wedge: Vancouver Island, British Columbia, Canada (average width, linear wedge, seaward and mixed vergence); Washington, USA (higher width, concave wedge, landward and mixed vergence); northern and central Oregon, USA (average width, linear and convex wedge, mixed and seaward vergence); and southern Oregon and northern California, USA (lower width, convex wedge, seaward and mixed vergence). Variabil­ity in outer wedge morphology and structure is broadly associated with along-strike megathrust segmentation inferred from differences in oceanic asthenospheric velocities, patterns of episodic tremor and slow slip, GPS models of plate locking, and the distribution of seismicity near the plate interface. In more detail, our results appear to delin­eate the extent, geometry, and lithology of dynamic and static backstops along the margin. Varying backstop configurations along the Cascadia mar­gin are interpreted to represent material-strength contrasts within the wedge that appear to regulate the along- and across-strike taper and structural vergence in the outer wedge. We argue that the morphotectonic variability in the outer wedge may reflect spatial variations in shallow megathrust behavior occurring over roughly the last few million years. Comparing outer wedge taper along the Cascadia margin to a global compilation suggests that observations in the global catalog are not accurately representing the range of hetero­geneity within individual margins and highlights the need for detailed margin-wide morphotectonic analyses of subduction zones worldwide.

British Columbia, California, Oregon, Washington↗
USGS2020

The potential of using dynamic strains in earthquake early warning applications

We investigate the potential of using borehole strainmeter data from the Network of the Americas (NOTA) and the U.S. Geological Survey networks to estimate earthquake moment magnitudes for earthquake early warning (EEW) applications. We derive an empirical equation relating peak dynamic strain, earthquake moment magnitude, and hypocentral distance, and investigate the effects of different types of instrument calibration on model misfit. We find that raw (uncalibrated) strains fit the model as accurately as calibrated strains. We test the model by estimating moment magnitudes of the largest two earthquakes in the July 2019 Ridgecrest earthquake sequence—the M 6.4 foreshock and the M 7.1 mainshock—using two strainmeters located within ∼50 km of the rupture. In both the cases, the magnitude based on the dynamic strain component is within ∼0.1–0.4 magnitude units of the catalog moment magnitude. We then compare the temporal evolution of our strain‐derived magnitudes for the largest two Ridgecrest events to the real‐time performance of the ShakeAlert EEW System (SAS). The final magnitudes from NOTA borehole strainmeters are close to SAS real‐time estimates for the M 6.4 foreshock, and significantly more accurate for the M 7.1 mainshock.

British Columbia, California, Oregon, Washington↗
USGS2014

Slip rate and tremor genesis in Cascadia

At many plate boundaries, conditions in the transition zone between seismogenic and stable slip produce slow earthquakes. In the Cascadia subduction zone, these events are consistently observed as slow, aseismic slip on the plate interface accompanied by persistent tectonic tremor. However, not all slow slip at other plate boundaries coincides spatially and temporally with tremor, leaving the physics of tremor genesis poorly understood. Here we analyze seismic, geodetic, and strainmeter data in Cascadia to observe for the first time a large, tremor-generating slow earthquake change from tremor-genic to silent and back again. The tremor falls silent at reduced slip speeds when the migrating slip front pauses as it loads the stronger adjacent fault segment to failure. The finding suggests that rheology and slip-speed-regulated stressing rate control tremor genesis, and the same section of fault can slip both with and without detectable tremor, limiting tremor's use as a proxy for slip.

British Columbia, California, Oregon, Washington↗
USGS2012

Turbidite event history—Methods and implications for Holocene paleoseismicity of the Cascadia subduction zone

Turbidite systems along the continental margin of Cascadia Basin from Vancouver Island, Canada, to Cape Mendocino, California, United States, have been investigated with swath bathymetry; newly collected and archive piston, gravity, kasten, and box cores; and accelerator mass spectrometry radiocarbon dates. The purpose of this study is to test the applicability of the Holocene turbidite record as a paleoseismic record for the Cascadia subduction zone. The Cascadia Basin is an ideal place to develop a turbidite paleoseismologic method and to record paleoearthquakes because (1) a single subduction-zone fault underlies the Cascadia submarine-canyon systems; (2) multiple tributary canyons and a variety of turbidite systems and sedimentary sources exist to use in tests of synchronous turbidite triggering; (3) the Cascadia trench is completely sediment filled, allowing channel systems to trend seaward across the abyssal plain, rather than merging in the trench; (4) the continental shelf is wide, favoring disconnection of Holocene river systems from their largely Pleistocene canyons; and (5) excellent stratigraphic datums, including the Mazama ash and distinguishable sedimentological and faunal changes near the Pleistocene-Holocene boundary, are present for correlating events and anchoring the temporal framework. Multiple tributaries to Cascadia Channel with 50- to 150-km spacing, and a wide variety of other turbidite systems with different sedimentary sources contain 13 post-Mazama-ash and 19 Holocene turbidites. Likely correlative sequences are found in Cascadia Channel, Juan de Fuca Channel off Washington, and Hydrate Ridge slope basin and Astoria Fan off northern and central Oregon. A probable correlative sequence of turbidites is also found in cores on Rogue Apron off southern Oregon. The Hydrate Ridge and Rogue Apron cores also include 12-22 interspersed thinner turbidite beds respectively. We use 14 C dates, relative-dating tests at channel confluences, and stratigraphic correlation of turbidites to determine whether turbidites deposited in separate channel systems are correlative - triggered by a common event. In most cases, these tests can separate earthquake-triggered turbidity currents from other possible sources. The 10,000-year turbidite record along the Cascadia margin passes several tests for synchronous triggering and correlates well with the shorter onshore paleoseismic record. The synchroneity of a 10,000-year turbidite-event record for 500 km along the northern half of the Cascadia subduction zone is best explained by paleoseismic triggering by great earthquakes. Similarly, we find a likely synchronous record in southern Cascadia, including correlated additional events along the southern margin. We examine the applicability of other regional triggers, such as storm waves, storm surges, hyperpycnal flows, and teletsunami, specifically for the Cascadia margin. The average age of the oldest turbidite emplacement event in the 10-0-ka series is 9,800±~210 cal yr B.P. and the youngest is 270±~120 cal yr B.P., indistinguishable from the A.D. 1700 (250 cal yr B.P.) Cascadia earthquake. The northern events define a great earthquake recurrence of ~500-530 years. The recurrence times and averages are supported by the thickness of hemipelagic sediment deposited between turbidite beds. The southern Oregon and northern California margins represent at least three segments that include all of the northern ruptures, as well as ~22 thinner turbidites of restricted latitude range that are correlated between multiple sites. At least two northern California sites, Trinidad and Eel Canyon/pools, record additional turbidites, which may be a mix of earthquake and sedimentologically or storm-triggered events, particularly during the early Holocene when a close connection existed between these canyons and associated river systems. The combined stratigraphic correlations, hemipelagic analysis, and 14 C framework suggest that the Cascadia margin has three rupture modes: (1) 19-20 full-length or nearly full length ruptures; (2) three or four ruptures comprising the southern 50-70 percent of the margin; and (3) 18-20 smaller southern-margin ruptures during the past 10 k.y., with the possibility of additional southern-margin events that are presently uncorrelated. The shorter rupture extents and thinner turbidites of the southern margin correspond well with spatial extents interpreted from the limited onshore paleoseismic record, supporting margin segmentation of southern Cascadia. The sequence of 41 events defines an average recurrence period for the southern Cascadia margin of ~240 years during the past 10 k.y. Time-independent probabilities for segmented ruptures range from 7-12 percent in 50 years for full or nearly full margin ruptures to ~21 percent in 50 years for a southern-margin rupture. Time-dependent probabilities are similar for northern margin events at ~7-12 percent and 37-42 percent in 50 years for the southern margin. Failure analysis suggests that by the year 2060, Cascadia will have exceeded ~27 percent of Holocene recurrence intervals for the northern margin and 85 percent of recurrence intervals for the southern margin. The long earthquake record established in Cascadia allows tests of recurrence models rarely possible elsewhere. Turbidite mass per event along the Cascadia margin reveals a consistent record for many of the Cascadia turbidites. We infer that larger turbidites likely represent larger earthquakes. Mass per event and magnitude estimates also correlate modestly with following time intervals for each event, suggesting that Cascadia full or nearly full margin ruptures weakly support a time-predictable model of recurrence. The long paleoseismic record also suggests a pattern of clustered earthquakes that includes four or five cycles of two to five earthquakes during the past 10 k.y., separated by unusually long intervals. We suggest that the pattern of long time intervals and longer ruptures for the northern and central margins may be a function of high sediment supply on the incoming plate, smoothing asperities, and potential barriers. The smaller southern Cascadia segments correspond to thinner incoming sediment sections and potentially greater interaction between lower-plate and upper-plate heterogeneities. The Cascadia Basin turbidite record establishes new paleoseismic techniques utilizing marine turbidite-event stratigraphy during sea-level highstands. These techniques can be applied in other specific settings worldwide, where an extensive fault traverses a continental margin that has several active turbidite systems.

British Columbia, California, Oregon, Washington↗
USGS2011

A Geo-referenced 3D model of the Juan de Fuca Slab and associated seismicity

We present a Geographic Information System (GIS) of a new 3-dimensional (3D) model of the subducted Juan de Fuca Plate beneath western North America and associated seismicity of the Cascadia subduction system. The geo-referenced 3D model was constructed from weighted control points that integrate depth information from hypocenter locations and regional seismic velocity studies. We used the 3D model to differentiate earthquakes that occur above the Juan de Fuca Plate surface from earthquakes that occur below the plate surface. This GIS project of the Cascadia subduction system supersedes the one previously published by McCrory and others (2006). Our new slab model updates the model with new constraints. The most significant updates to the model include: (1) weighted control points to incorporate spatial uncertainty, (2) an additional gridded slab surface based on the Generic Mapping Tools (GMT) Surface program which constructs surfaces based on splines in tension (see expanded description below), (3) double-differenced hypocenter locations in northern California to better constrain slab location there, and (4) revised slab shape based on new hypocenter profiles that incorporate routine depth uncertainties as well as data from new seismic-reflection and seismic-refraction studies. We also provide a 3D fly-through animation of the model for use as a visualization tool.

British Columbia, California, Oregon, Washington↗