GSA Annual Meeting in Seattle, Washington, USA - 2017

Paper No. 58-1
Presentation Time: 1:30 PM

THE TIMING AND DEFORMATION OF CASCADIA SUBDUCTION ZONE EARTHQUAKES FROM COASTAL MARSH RECORDS OF OREGON (Invited Presentation)


ENGELHART, Simon E., Department of Geosciences, University of Rhode Island, Kingston, RI 02881, CAHILL, Niamh, Department of Biostatistics and Epidemiology, University of Massachusetts Amherst, Amherst, MA 01003, HAWKES, Andrea, Earth and Ocean Sciences, University of North Carolina Wilmington, 601 S. College Road DeLoach Hall, Wilmington, NC 28403, HORTON, Benjamin P., Sea Level Research, Department of Marine and Coastal Science, Rutgers University, 71 Dudley Road, New Brunswick, NJ 08901, KEMP, Andrew C., Department of Earth and Ocean Sciences, Tufts University, Medford, MA 02155, MILKER, Yvonne, Department of Earth Sciences, Center for Earth System Research and Sustainability, Institut for Geology, Hamburg, 20146, Germany, NELSON, Alan R., Geologic Hazards Science Center, U.S.G.S., Golden, CO 80401, PADGETT, J. Scott, Department of Geosciences, University of Rhode Island, Woodward Hall, 9 E Alumni Avenue, Kingston, RI 02881, WANG, Kelin, School of Earth and Ocean Sciences, University of Victoria, Victoria, BC V8P 5C2, Canada and WITTER, Robert C., Alaska Science Center, U.S.G.S., Anchorage, AK 99508, engelhart@uri.edu

Determining the timing of Cascadia subduction zone earthquakes and their associated magnitude of coseismic deformation remains a challenge for paleoseismology. Here, we use (1) new coastal records; (2) microfossil-based paleogeodetic techniques (errors < 0.3m); and (3) high-resolution AMS radiocarbon ages (errors < ±20 years) to better understand past ruptures and, therefore, the seismic hazard posed by the Cascadia subduction zone.

We develop a new regional foraminiferal-based Bayesian Transfer Function (BTF) for inferring the magnitude of coseismic subsidence from these records. The BTF method has the advantage of incorporating information readily available from secondary proxies (e.g., sedimentology) and allows non-unimodal taxa-elevation relationships. The use of secondary priors compensates for the effect of mixing across coseismic contacts. Application of the validated BTF to records pertaining to the AD 1700 rupture suggests greater coseismic subsidence in northern Oregon and southern Washington than previously estimated. However, these results still indicate heterogeneous rupture in this event and a magnitude near 9 as inferred by Wang et al. (2013 JGR).

To extend this analysis to earlier earthquakes we must produce robust correlations between earthquakes at individual sites. To achieve this, we compile a new dataset of high-quality AMS radiocarbon dates on earthquakes recorded in sediment beneath Oregon coastal marshes. We employ two age models Bchron and Oxcal. Model results from both approaches are consistent when maximum and minimum ages are available and closely-spaced in age, and are proximal to the sedimentary contacts inferred to record coseismic subsidence and postseismic recovery. When these conditions are not met, the age estimated for an earthquake by Bchron and Oxcal may differ. This is most likely due to the additional sedimentation rate information used by Bchron to refine earthquake timings. These differences are most pronounced for an earthquake(s) occurring ~800 cal yrs BP with Bchron suggesting multiple earthquakes closely spaced in time in Oregon. In contrast, Oxcal results allow for both single and multiple-rupture scenarios. Both Oxcal and Bchron agree that more than one earthquake is recorded in Oregon between ~950 and 1300 cal yrs BP.