Search for: All records

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 the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally. 

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. 

Uplift of the overriding plate at a subduction zone denotes interseismic strain accumulation, which is subsequently released during a megathrust earthquake. Although most interseismic strain is thought to be released elastically, observations of uplifted coastal regions at subduction zones worldwide indicate that some strain may result in permanent uplift. The Grays Harbor and Willapa Bay (Washington, USA) coastal region of the Cascadia subduction zone hosts flights of marine terraces testifying to late Pleistocene rock uplift. Our new detailed mapping of the marine terraces recognizes nine new units, including estuarine and fluvial sediments. Luminescence dating, relative age based on soil maturity and terrace elevation, and an evaluation of previous ages from fossil shells collectively constrain the probable ages of three estuarine units to sea-level high stands during Marine Isotope Stages 5a, 5c, and 5e. We estimate an average uplift rate of 0.4 ± 0.1 mm/yr for the terraced estuarine units, consistent with other Pleistocene uplift and incision rates in Cascadia. When compared with observed interseismic vertical deformation, these rates suggest that about one-tenth of interseismic strain may become permanent. The values are permissible within the uncertainties of uplift based on regional estimates of interseismic vertical strain rates and of coseismic subsidence. 

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.