skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on April 28, 2026

Title: Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise
Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1 % floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km² (low), 160 km² (medium), or 300 km² (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km² (low), 240 km² (medium), or 370 km² (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from earthquake-cycle and climate-driven sea-level rise, and provide critical insights for tectonically active coastlines globally.  more » « less
Award ID(s):
2325312
PAR ID:
10586375
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences of the United States of America
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences of the United States of America
ISSN:
0027-8424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km2(low), 160 km2(medium), or 300 km2(high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km2(low), 240 km2(medium), or 370 km2(high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally. 
    more » « less
  2. Abstract Stratigraphic, lithologic, foraminiferal, and radiocarbon analyses indicate that at least four abrupt mud-over-peat contacts are recorded across three sites (Jacoby Creek, McDaniel Creek, and Mad River Slough) in northern Humboldt Bay, California, USA (∼44.8°N, −124.2°W). The stratigraphy records subsidence during past megathrust earthquakes at the southern Cascadia subduction zone ∼40 km north of the Mendocino Triple Junction. Maximum and minimum radiocarbon ages on plant macrofossils from above and below laterally extensive (>6 km) contacts suggest regional synchroneity of subsidence. The shallowest contact has radiocarbon ages that are consistent with the most recent great earthquake at Cascadia, which occurred at 250 cal yr B.P. (1700 CE). Using Bchron and OxCal software, we model ages for the three older contacts of ca. 875 cal yr B.P., ca. 1120 cal yr B.P., and ca. 1620 cal yr B.P. For each of the four earthquakes, we analyze foraminifera across representative mud-over-peat contacts selected from McDaniel Creek. Changes in fossil foraminiferal assemblages across all four contacts reveal sudden relative sea-level (RSL) rise (land subsidence) with submergence lasting from decades to centuries. To estimate subsidence during each earthquake, we reconstructed RSL rise across the contacts using the fossil foraminiferal assemblages in a Bayesian transfer function. The coseismic subsidence estimates are 0.85 ± 0.46 m for the 1700 CE earthquake, 0.42 ± 0.37 m for the ca. 875 cal yr B.P. earthquake, 0.79 ± 0.47 m for the ca. 1120 cal yr B.P. earthquake, and ≥0.93 m for the ca. 1620 cal yr B.P. earthquake. The subsidence estimate for the ca. 1620 cal yr B.P. earthquake is a minimum because the pre-subsidence paleoenvironment likely was above the upper limit of foraminiferal habitation. The subsidence estimate for the ca. 875 cal yr B.P. earthquake is less than (<50%) the subsidence estimates for other contacts and suggests that subsidence magnitude varied over the past four earthquake cycles in southern Cascadia. 
    more » « less
  3. Abstract The risk of compound coastal flooding in the San Francisco Bay Area is increasing due to climate change yet remains relatively underexplored. Using a novel hybrid statistical-dynamical downscaling approach, this study investigates the impacts of climate change induced sea-level rise and higher river discharge on the magnitude and frequency of flooding events as well as the relative importance of various forcing drivers to compound flooding within the Bay. Results reveal that rare occurrences of flooding under the present-day climate are projected to occur once every few hundred years under climate change with relatively low sea-level rise (0.5 m) but would become annual events under climate change with high sea-level rise (1.0 to 1.5 m). Results also show that extreme water levels that are presently dominated by tides will be dominated by sea-level rise in most locations of the Bay in the future. The dominance of river discharge to the non-tidal and non-sea-level rise driven water level signal in the North Bay is expected to extend ~15 km further seaward under extreme climate change. These findings are critical for informing climate adaptation and coastal resilience planning in San Francisco Bay. 
    more » « less
  4. Abstract Globally, coastal communities experience flood hazards that are projected to worsen from climate change and sea level rise. The 100-year floodplain or record flood are commonly used to identify risk areas for planning purposes. Remote communities often lack measured flood elevations and require innovative approaches to estimate flood elevations. This study employs observation-based methods to estimate the record flood elevation in Alaska communities and compares results to elevation models, infrastructure locations, and sea level rise projections. In 46 analyzed communities, 22% of structures are located within the record floodplain. With sea level rise projections, this estimate increases to 30–37% of structures by 2100 if structures remain in the same location. Flood exposure is highest in western Alaska. Sea level rise projections suggest northern Alaska will see similar flood exposure levels by 2100 as currently experienced in western Alaska. This evaluation of record flood height, category, and history can be incorporated into hazard planning documents, providing more context for coastal flood exposure than previously existed for Alaska. This basic flood exposure method is transferable to other areas with similar mapping challenges. Identifying current and projected hazardous zones is essential to avoid unintentional development in floodplains and improve long-term safety. 
    more » « less
  5. Abstract 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 14C 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. 
    more » « less