ALBERT

Description: Coastal stratigraphy from the Pacific Northwest of the United States contains evidence of sudden subsidence during ruptures of the Cascadia subduction zone. Transfer functions (empirical relationships between assemblages and elevation) can convert microfossil data into coastal subsidence estimates. Coseismic deformation models use the subsidence values to constrain earthquake magnitudes. To test the response of foraminifera, the accuracy of the transfer function method, and the presence of a pre-seismic signal, we simulated a great earthquake near Coos Bay, Oregon, by transplanting a bed of modern high salt-marsh sediment into the tidal flat, an elevation change that mimics a coseismic subsidence of 0.64 m. The transplanted bed was quickly buried by mud; after 12 mo and 5 yr, we sampled it for foraminifera. Reconstruction of the simulated coseismic subsidence using our transfer function was 0.61 m, nearly identical to the actual elevation change. Our transplant experiment, and additional analyses spanning the A.D. 1700 earthquake contact at the nearby Coquille River 15 km to the south, show that sediment mixing may explain assemblage changes previously interpreted as evidence of pre-seismic land-level change in Cascadia and elsewhere.

Description: The Sallys Bend swamp and marsh area on the central Oregon coast onshore of the Cascadia subduction zone contains a sequence of buried coastal wetland soils that extends back ~4500 yr B.P. The upper 10 of the 12 soils are represented in multiple cores. Each soil is abruptly overlain by a sandy deposit and then, in most cases, by greater than 10 cm of mud. For eight of the 10 buried soils, times of soil burial are constrained through radiocarbon ages on fine, delicate detritus from the top of the buried soil; for two of the buried soils, diatom and foraminifera data constrain paleoenvironment at the time of soil burial. We infer that each buried soil represents a Cascadia subduction zone earthquake because the soils are laterally extensive and abruptly overlain by sandy deposits and mud. Preservation of coseismically buried soils occurred from 4500 yr ago until ~500–600 yr ago, after which preservation was compromised by cessation of gradual relative sea-level rise, which in turn precluded drowning of marsh soils during instances of coseismic subsidence. Based on grain-size and microfossil data, sandy deposits overlying buried soils accumulated immediately after a subduction zone earthquake, during tsunami incursion into Sallys Bend. The possibility that the sandy deposits were sourced directly from landslides triggered upstream in the Yaquina River basin by seismic shaking was discounted based on sedimentologic, microfossil, and depositional site characteristics of the sandy deposits, which were inconsistent with a fluvial origin. Biostratigraphic analyses of sediment above two buried soils—in the case of two earthquakes, one occurring shortly after 1541–1708 cal. yr B.P. and the other occurring shortly after 3227–3444 cal. yr B.P.—provide estimates that coseismic subsidence was a minimum of 0.4 m. The average recurrence interval of subduction zone earthquakes is 420–580 yr, based on an ~3750–4050-yr-long record and seven to nine interearthquake intervals. The comparison of the Yaquina Bay earthquake record to similar records at other Cascadia coastal sites helps to define potential patterns of rupture for different earthquakes, although inherent uncertainty in dating precludes definitive statements about rupture length during earthquakes. We infer that in the first half of the last millennia, the northern Oregon part of the subduction zone had a different rupture history than the southern Oregon part of the subduction zone, and we also infer that at ca. 1.6 ka, two earthquakes closely spaced in time together ruptured a length of the megathrust that extends at least from southwestern Washington to southern Oregon.

Description: Characterizations of tsunami hazards along the Cascadia subduction zone hinge on uncertainties in megathrust rupture models used for simulating tsunami inundation. To explore these uncertainties, we constructed 15 megathrust earthquake scenarios using rupture models that supply the initial conditions for tsunami simulations at Bandon, Oregon. Tsunami inundation varies with the amount and distribution of fault slip assigned to rupture models, including models where slip is partitioned to a splay fault in the accretionary wedge and models that vary the updip limit of slip on a buried fault. Constraints on fault slip come from onshore and offshore paleoseismological evidence. We rank each rupture model using a logic tree that evaluates a model’s consistency with geological and geophysical data. The scenarios provide inputs to a hydrodynamic model, SELFE, used to simulate tsunami generation, propagation, and inundation on unstructured grids with 〈5–15 m resolution in coastal areas. Tsunami simulations delineate the likelihood that Cascadia tsunamis will exceed mapped inundation lines. Maximum wave elevations at the shoreline varied from ~4 m to 25 m for earthquakes with 9–44 m slip and M w 8.7–9.2. Simulated tsunami inundation agrees with sparse deposits left by the A.D. 1700 and older tsunamis. Tsunami simulations for large (22–30 m slip) and medium (14–19 m slip) splay fault scenarios encompass 80%–95% of all inundation scenarios and provide reasonable guidelines for land-use planning and coastal development. The maximum tsunami inundation simulated for the greatest splay fault scenario (36–44 m slip) can help to guide development of local tsunami evacuation zones.