Geology Reports⌕ Search

Geology topics

Robert C. Witter

Publications and source records attributed to Robert C. Witter.

At least 19 recordsLinked to original sources
USGS2026

The Cascadia Region Earthquake Science Center (CRESCENT): Advancing understanding of Cascadia's earthquake hazards

The Cascadia Subduction Zone (CSZ) is capable of producing large megathrust earthquakes and tsunamis, posing significant risks to the Pacific Northwest. This paper examines the seismic characteristics of the CSZ, including its offshore coupling, sparse seismicity, and a complex paleoseismic record. While recent advances in geophysical methods and data collection have improved our understanding of subduction zone processes, important uncertainties remain regarding rupture dynamics, ground shaking, and tsunami impacts. The Cascadia Region Earthquake Science Center (CRESCENT) has been established to address these challenges through coordinated research, education, and community engagement. By integrating scientific studies, fostering collaboration, and leveraging new technologies, CRESCENT seeks to improve hazard assessment and resilience in the region. This paper outlines the primary challenges and opportunities in advancing subduction zone science and enhancing societal preparedness for future CSZ events.

California, Oregon, Washington↗
USGS2026

Subaerially exposed Iceberg Lake sediments: An exceptional record of historical subaqueous earthquake disturbance at the eastern edge of the Alaskan-Aleutian subduction zone

Paleoseismic records are limited in the Yakutat Terrane (eastern edge of the Alaskan-Aleutian subduction zone) due to the extensive ice cover that hinders traditional methods such as trenching of the faults, but lacustrine sediments offer an alternative archive. We investigated lakebed sediments exposed after recent outburst floods (1999 CE) at Iceberg Lake, a glacier-dammed basin whose stratigraphy was revealed after the lake drained and partially eroded. We logged outcrops across the basin and sampled their sediments. Between annually laminated background deposits, we identified sediment gravity flow beds and in-situ soft-sediment deformation structures (convolutions, sand blows, and fractures) interpreted to be earthquake-induced. Our age model links some of the uppermost depositional and deformational events to the 1958, 1964 and 1979 CE earthquakes. These results demonstrate that Iceberg Lake was a sensitive recorder of seismic shaking and its sediments hold strong potential for producing a paleoseismic record for the northern Yakutat Terrane.

Alaska↗
USGS2025

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 2 (low), 160 km 2 (medium), or 300 km 2 (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 2 (low), 240 km 2 (medium), or 370 km 2 (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.

California, Oregon, Washington↗
USGS2025

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii↗
USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗
USGS2024

An Alaska-Aleutian subduction zone interface earthquake recurrence model from geology and geodesy

We present an earthquake recurrence model for the subduction interface of the Alaska–Aleutian subduction zone from geodetic and paleoseismic data. To capture variations in rupture behavior along strike, we define fault sections based on geodetic coupling, prehistoric earthquake and tsunami recurrence, historical ruptures, and geologic and geophysical structure. In the 1964 M w 9.2 rupture area, which spans four fault sections, recurrence rates for section participation in presumed great ( M w 8.5+) events vary from ∼600 years (Prince William Sound section) to ∼380 years (Kodiak section), geodetic character varies substantially along strike, and geologic evidence indicates that rupture patches vary in space and time. Westward along the Semidi section, the recurrence of large, tsunamigenic ruptures is more frequent (∼220 years) based on geologic and geodetic data than previously assumed. The seismic potential of the Shumagin section, an area of low coupling, remains enigmatic despite a large ( M w 7.8) rupture in 2020. Prehistoric tsunami data indicate that inundation recurrence in the Fox Islands is ∼210 years. Paleoseismic data are lacking west of the Fox Islands, so rupture rates along the western 1900 km of the subduction interface to Komandorski are inferred from geodetic constraints.

Alaska↗
USGS2024

Constraining the earthquake recording threshold of intraslab earthquakes with turbidites in southcentral Alaska’s lakes and fjords

Strong ground motion from intraslab earthquakes, which do not produce primary paleoseismic evidence, may initiate gravity-driven turbidity flows in subaqueous basins. The resulting deposits (turbidites) can provide a paleoseismic proxy if the conditions that initiate these flows are known. To better constrain the initiating conditions, we use two recent intraslab earthquakes in southcentral Alaska, the M w 7.1 30 November 2018 Anchorage earthquake and the M w 7.1 24 January 2016 Iniskin earthquake, as calibration events. Through a multilake investigation, we document the occurrence, or the absence, of earthquake-generated turbidity flows from these two events. Both earthquakes are recorded by centimeter-scale turbidites that can be differentiated from climatically generated deposits, as well as other seismic sources based on deposit thickness, sedimentological properties, and deposit age. We show that a Modified Mercalli Intensity (MMI) of ∼V–V1/2 is the minimum shaking intensity required to generate localized sediment remobilization from deltaic slopes, and an MMI of ∼V1/2 is required to produce a deposit of sufficient thickness that a seismic origin can be confidently assigned. The documentation of seismically generated deposits in quick succession (∼2 years) with diagnostic features highlights the utility of using recent earthquakes as calibration events to investigate the subaqueous response to strong ground motion.

Alaska↗
USGS2024

Oblique contraction along the fastest ocean-continent transform plate boundary focuses rock uplift west of the Fairweather fault, southeast Alaska

Contraction along the Yakutat–North America plate boundary drives 4.6–9.0 mm/year Holocene rock uplift rates along Earth's fastest slipping (≥49 mm/year) ocean–continent transform fault, the Fairweather Fault. Between Icy Point and Lituya Bay, the near-vertical Fairweather fault focuses rock uplift and rapid right-lateral slip by accommodating both vertical and fault-parallel strain during oblique-slip and separate, predominantly strike-slip ruptures. Unusually high uplift rates, indicated by radiocarbon and luminescence dating, result from a 10-km-wide, asymmetric, positive flower structure along a 20°, ∼30-km-long restraining double bend in the Fairweather fault. The principal reverse fault in the flower structure, the offshore, blind Icy Point–Lituya Bay fault, ruptures no more than every 460–1040 years evidenced by uplifted Holocene shorelines. Maximum 3–5 m coseismic uplifts imply 3.1–10 m dip slip per event and earthquake magnitudes of M w 7.0–7.5. The Yakutat block collides obliquely into North America, and our model entails oblique slip on the Fairweather fault with and without corupture on the reverse fault. Oblique slip is evident by vertically offset (>25 m) fluvial and marine terraces and by the primary Fairweather fault strand that strikes >20° to the west of plate-boundary motion.

Alaska↗
USGS2024

The 2023 Alaska National Seismic Hazard Model

US Geological Survey (USGS) National Seismic Hazard Models (NSHMs) are used extensively for seismic design regulations in the United States and earthquake scenario development, as well as risk assessment and mitigation for both buildings and infrastructure. This 2023 update of the long-term, time-independent Alaska NSHM includes substantial changes to both the earthquake rupture forecast (ERF) and ground motion models (GMMs). The ERF includes numerous additions to the finite-fault model, considers two deformation models, and introduces updated declustering and smoothing algorithms in the gridded background seismicity model. For the Alaska–Aleutian subduction zone, megathrust earthquakes occur on an updated structural and segmentation model, and the moment magnitude (M) 8+ rupture and rate model include a logic tree branch that considers slip rates derived from geodetic models of interface coupling. The megathrust model considers multiple models of down-dip width, and magnitudes are computed using newly developed scaling relations. For subduction intraslab events and subduction interface events with M < 7, the 2023 update uses a smoothed seismicity model with rupture depths derived from Slab2. The 2023 model updates GMMs in all tectonic settings using the recently published Next Generation Attenuation Subduction (NGA-Sub) GMMs for subduction interface and intraslab events, and the NGA-West2 GMMs for active crustal settings. Collectively, additions and updates to the Alaska NSHM result in hazard increases across most of south-central Alaska relative to the previous model, published in 2007. These changes are primarily due to the adoption of updated rate models for the large-magnitude interface events and the NGA-Sub GMMs that have much higher aleatory variability (sigma), consistent with global observations, and that include models of epistemic uncertainty.

Alaska↗
USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗
USGS2024

Repeated coseismic uplift of coastal lagoons above the Patton Bay Splay Fault System, Montague Island, Alaska, USA

Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.

Alaska↗
USGS2024

A great tsunami earthquake component of the 1957 Aleutian Islands earthquake

The great 1957 Aleutian Islands earthquake ruptured ∼1200 km of the plate boundary along the Aleutian subduction zone and produced a destructive tsunami across Hawaiʻi. Early seismic and tsunami analyses indicated that large megathrust fault slip was concentrated in the western Aleutian Islands, but tsunami waves generated by slip in the west cannot explain the large observed runup in Hawaiʻi far to the southeast. Recently mapped 1957 geologic deposits on eastern Aleutian Islands suggest occurrence of very large nearby slip. Jointly modeling tsunami runup along the eastern Aleutian and Hawaiian Islands together with tide gauge recordings across the Pacific resolves 12-26 m shallow slip along 600 km of the eastern Aleutian Islands in addition to modest, deeper western slip inferred from seismic records. The eastern near-trench slip results in an M W 8.3-8.6 tsunami earthquake component of the M W 8.6-8.8 rupture, comparable in size to the adjacent 1946 Aleutian tsunami earthquake to the east. The reexamination of the 1957 rupture confirms the tsunami hazards posed by the eastern Aleutian subduction zone to Hawaiʻi and lays the groundwork for investigation of large prehistoric earthquakes through modeling tsunami runup inferred from stratigraphic observations to constrain their rupture processes.

Alaska↗
USGS2023

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra↗
USGS2023

Dating the penultimate great earthquake in south-central Alaska using tree-ring crossdating and radiocarbon wiggle-matching

A forest bed of tree stumps currently in the intertidal zone at Girdwood, south-central Alaska, records coseismic submergence during the penultimate great earthquake. Tree-ring samples from ten spruce stumps were crossdated to develop a 149-year-long ring-width chronology. Radiocarbon wiggle-matching found that single-ring ages from the chronology were offset 28 ± 7 years older than the IntCal20 calibration curve and that the last ring of the chronology dated as 1169 to 1189 CE (781–761 cal. yr. BP) at the 95% confidence level. Bark was observed on some stumps, six samples had the same year for the last growth ring, and so this wiggle-match date is also the best estimate of the date of the penultimate great earthquake. This date is in good agreement with a date for this event in a seismo-turbidite record from Skilak Lake but not with previous dates from Bayesian models of maximum- and minimum-limiting ages from coastal salt marshes. Reanalysis of the coastal salt marsh ages with the data grouped by area, context and material found that outer wood samples from stumps at coseismic submergence sites and a Bayesian limiting age model based on just herbaceous plant ages from Turnagain Arm and the Copper River area are both consistent with our wiggle-match date. Furthermore, coseismic emergence ages from Cape Suckling and Yakataga are older than the penultimate earthquake and so likely relate to an earlier uplift event in this eastern area. The rupture extent during the penultimate great earthquake appears to have been less than in the 1964 great earthquake and the interseismic interval between these two events was 785 ± 10 years.

Alaska↗
USGS2023

Rapid shallow megathrust afterslip from the 2021 M8.2 Chignik, Alaska earthquake revealed by seafloor geodesy

The shallower portions of subduction zone megathrust faults host Earth’s most hazardous tsunamigenic earthquakes, yet understanding how and when they slip remains elusive because of challenges making seafloor observations. We performed Global Navigation Satellite System Acoustic seafloor geodetic surveys before and ~2.5 months after the 29 July 2021 M w (moment magnitude) 8.2 Chignik, Alaska, earthquake and determine ~1.4 meters cumulative co- and post-seismic horizontal displacement ~60 kilometers from the megathrust front. Only for the 2011 M w 9 Tohoku event have closer subduction zone earthquake displacements been observed. We estimate ~2 to 3 meters of megathrust afterslip shallower than 20 kilometers, a portion of the megathrust on which both inter- and co-seismic slip likely had occurred previously. Our analysis demonstrates that by 2.5 months, shallower and deeper moment had effectively equilibrated on the megathrust, suggesting that its tsunamigenic potential remains no more elevated than before the earthquake.

Alaska↗
USGS2022

Unravelling a 2300 year long sedimentary record of megathrust and intraslab earthquakes in proglacial Skilak Lake, south-central Alaska

Seismic hazards in subduction settings typically arise from megathrust, intraslab and crustal earthquake sources. Despite the frequent occurrence of intraslab earthquakes in subduction zones and their potential threat to communities, their long-term recurrence behaviour is barely studied. Sedimentary sequences in lakes may register ground shaking from different seismic sources. This study investigates two long sediment cores (13 m and 16 m) from Skilak Lake, a proglacial lake in south-central Alaska, to evaluate whether different seismic sources leave a distinct imprint. The sedimentary record shows a continuously varved sediment sequence, occasionally interrupted by turbidites, slump deposits and tephra beds. Turbidites and slump deposits were objectively identified using a statistical outlier analysis on varve thickness. The earthquake origin of these deposits was ascertained by resemblance with deposits induced by instrumentally recorded earthquakes (for example, 1964 ce M w 9.2 megathrust and 1954 ce M w 6.4 intraslab earthquakes) and correlation with multiple coeval landslide deposits on sub-bottom profiles. The Skilak Lake record chronicles 19 earthquakes with moderate to very high confidence level in the past 1350 years. The sedimentary evidence of instrumentally-recorded intraslab and megathrust earthquakes within the past 70 years demonstrates that not only megathrust earthquakes, but also past intraslab events are recorded. Although reported seismic intensities at Skilak Lake are comparable for the 1964 ce megathrust and the 1954 ce intraslab earthquakes, the long duration and low frequency content of seismic ground motion during megathrust earthquakes facilitate the triggering of multiple, voluminous landslides and the generation of megaturbidites. In contrast, the shorter duration and higher frequency source spectrum of intraslab earthquakes may only induce surficial slope remobilization and the generation of thinner turbidites. This study demonstrates that the sedimentary record of Skilak Lake has the potential to decipher multiple seismic sources, which opens possibilities for a comprehensive seismic hazard analysis for south-central Alaska.

Alaska↗
USGS2022

Reproducibility and variability of earthquake subsidence estimates from saltmarshes of a Cascadia estuary

We examine fossil foraminiferal assemblages from 20 sediment cores to assess sudden relative sea-level (RSL) changes across three mud-over-peat contacts at three salt marshes in northern Humboldt Bay, California (~44.8°N, -124.2°W). We use a validated foraminiferal-based Bayesian transfer function to evaluate the variability of subsidence stratigraphy at a range of 30-6000 m across an estuary. We use the consistency in RSL reconstructions to support estimates of coseismic subsidence from megathrust earthquakes. To assess the variability of subsidence estimates, we analyzed: nine examples of the 1700 CE earthquake (average of 0.64 ±0.14 m subsidence; range of 0.24 ±0.27 to 1.00 ±0.44 m), five examples of the ca. 875 cal a BP earthquake (average of 0.43 ±0.16 m; range of 0.41 ±0.36 to 0.48 ±0.39 m), and six examples of the ca. 1120 cal a BP earthquake (average of 0.70±0.18 m; range of 0.47 ±0.36 to 0.80 ±0.49 m). Our subsidence estimate results suggest ~±0.3 m of within-site (intra-site) variability, which is consistent with previous research. We also identify inconsistencies between sites (inter-site) at northern Humboldt Bay greater than one-sigma uncertainties, driven by variable foraminiferal assemblages in the mud overlying the 1700 CE subsidence contact. Therefore, we recommend at least two quantitative microfossil reconstructions across the same stratigraphic sequence from different marsh sites within an estuary to account for estimate variability and provide increased confidence in vertical coseismic deformation estimates. Our results have broad implications for quantitative, microfossil-based reconstructions of coseismic subsidence at temperate coastlines globally.

California↗