Search USGSSearch

USGS · 70199863

Wetland stratigraphic evidence for variable megathrust earthquake rupture modes at the Cascadia subduction zone

Abstract

Although widespread agreement that the Cascadia subduction zone produces great earthquakes of magnitude 8 to 9 was reached decades ago, debate continues about the rupture lengths, magnitudes, and frequency of megathrust earthquakes recorded by wetland stratigraphy fringing Cascadia’s estuaries. Correlation of such coastal earthquake evidence along the subduction zone has largely relied on relative position in a stratigraphic sequence and maximum-limiting 14C ages with errors of decades to hundreds of years. Offshore, a 10,000-year record of turbidites in marine cores is interpreted as an archive of strong shaking from great earthquakes, with an average frequency of about 500 years in northern Cascadia versus 200-300 years in southern Oregon and northern California. Onshore, fewer events marked by sharp (<3 mm) peat-mud (mud-over-peat) contacts in tidal wetland stratigraphic sequences have been widely inferred to record sudden relative sea-level rise due to coseismic subsidence during megathrust earthquakes: 4-7 sharp subsidence contacts in 3500 years at estuaries in northern Oregon and southern Washington (500-800 year average recurrence), and 9-12 sharp subsidence contacts in over 6000 years in sequences in central and southern Oregon (500-900 year average recurrence). Improved understanding of the onshore and offshore records is critical to the assessment of earthquake hazard in western North America and of tsunami hazard in the Pacific basin. However, because dating the turbidite record is inherently much less precise than are age models for subsidence events in the most thoroughly studied tidal wetland sequences, accurate reconstruction of the times of Cascadia’s great earthquakes depends on the ages from the onshore record. Although methods to reduce uncertainty in the limits of resolution of tidal stratigraphy for recording earthquakes of a particular magnitude, and ways to distinguish earthquake subsidence stratigraphic contacts from non-seismic contacts, have been discussed for decades (e.g., Nelson, 1992; Atwater, 1992; Darienzo et al., 1994; Nelson et al., 1996a; Atwater and Hemphill-Haley, 1997; Witter et al., 2001; Kelsey et al., 2002; Nelson et al., 2006; Graehl et al., 2014; Milker et al., 2016), consensus about the threshold of resolution (minimum identifiable evidence of an earthquake) of tidal stratigraphy and, therefore, the completeness of Cascadia’s coastal record of great earthquakes, remains elusive. Although the most distinct, widespread contacts likely record close to a meter of coastal subsidence during the greatest megathrust earthquakes (e.g., M8.8-M9), other contacts may record <0.5 m of subsidence onshore of patches of low stress release on the megathrust during great earthquakes, during lesser megathrust earthquakes (e.g., M8.2-8.6), or from localized subsidence near upper-plate faults that slip during or independently of megathrust earthquakes (Nelson et al., 1996b; Wang et al., 2013; Kemp et al., 2018). At the Siuslaw River estuary in central Oregon (lat. 43.97°) a stratigraphy of 9-12 peat-mud contacts, similar to those described from many Cascadia estuaries, may record a greater number of earthquakes during the past 2000 years than at any other of the tens of tidal wetland sites to the north and south. Here, as well as at many tens of other Cascadia tidal wetland sites, peat-mud contacts mark the tops of couplets of tidal flat and low marsh mud gradually shoaling upward into middle and high marsh peat. Along core transects across an 800-m-wide, island marsh in the Siuslaw River, we traced the 9 most continuous of 12-15 peaty beds dating from the past 2000 years for 250-500 m, but we had difficulty correlating the 3-6 intervening beds >50-100 m. We attribute the sharper, more extensive upper contacts on peaty beds—two capped by sandy beds probably deposited by tsunamis—to sudden coseismic subsidence of middle and high marshes, but origins for other upper and lower contacts boundi

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Alan R. Nelson, Robert C. Witter, Simon Englehart, Andrea Hawkers, Benjamin P. Horton. 2018. Wetland stratigraphic evidence for variable megathrust earthquake rupture modes at the Cascadia subduction zone. https://pubs.usgs.gov/publication/70199863

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Time-to-depth conversion of seismic-reflection data from eastern Lake Superior and implications for the eastern arm of the Midcontinent Rift

Seismic-reflection data were acquired in the mid 1980s along several lines across eastern Lake Superior by industry and the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE) (Fig. 1). The lines form part of a larger network of crossing lines over the entire lake, which can be used to develop three-dimensional geologic models of the Mesoproterozoic Midcontinent Rift that lies below. To better interpret these lines, we developed velocity models to convert seismic reflections versus two-way travel time (TWTT) to reflections versus depth. In addition, the velocity models themselves provide insights into the structure of the Midcontinent Rift by recognizing common velocity ranges for certain rock types (Grauch, 2023).

eastern Lake Superior

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona

Early Miocene volcanic rocks and associated tectonics, Lava Hills and southern Bristol Mountains, California

Volcanic rocks of latest Oligocene to early Miocene age form an east-west belt across part of the central eastern Mojave Desert from the Whipple Mountains on the east to the Rosamond Hills on the west. We term this the central belt because it is separated from northern and southern belts by swaths with no volcanic rocks. Limited geochronologic data indicate that much of the belt is latest Oligocene and early Miocene in age, about 24 to 19 Ma, a finding that is consistent with these rocks being overlain by the 18.8 Ma Peach Spring Tuff in many places. We describe Miocene geology in a central area of the belt, in the Lava Hills, southern Bristol Mountains, and southern Old Dad Mountains. Sedimentary basins formed coeval with early andesite to rhyolite volcanism, progressing from fluvial and lacustrine tuffaceous sandstone to volcanic lavas, tuffs, and breccias, indicating that early basins formed proximal to volcanic edifices. Higher strata are fluvial and lacustrine with lavas punctuating the sequence. Although basins may partly have been formed within topographic lows bounded by volcanic domes, plateaus, and stratovolcanoes, consistent stratigraphic sections over wide areas indicate that tectonic basin evolution affected broad areas. The volcanic section is capped by local basalt flows and the regional Peach Spring Tuff. Limited data on normal faults support interpretations of early extensional basin development caused by northeast-southwest oriented stretching. Later extension caused stratal rotations, tilting early deposits down to the southwest. This tilted and subsequently beveled basin architecture was overlain by the youngest volcanic deposits, primarily rhyolite and basalt. The Peach Spring Tuff, 18.8 Ma, lies within this upper unit. Similar stratigraphic and structural relations are exposed in the nearby Marble Mountains and Van Winkle Mountain sections, reinforcing that a broad area underwent similar volcanism and tectonism. In our study area the upper unit is only very gently tilted except near dextral strike-slip faults of the eastern California shear zone. These late Miocene to Recent faults are represented as four main fault zones spaced about 5 km apart, representing distributed shear north of the Bristol Lake basin.

California