The Cascadia Subduction Zone (Figure 1, CSZ) is a 1,200 km plate boundary that extends from the Mendocino Triple Junction in Northern California to the tectonically complex region surrounding the Explorer Plate offshore of Northern Vancouver Island (Braunmiller & Nábělek, 2002;Savard et al., 2020). It accommodates roughly 4 cm of annual convergence between the Juan de Fuca and North American plates (DeMets et al., 2010). It is well established that the CSZ hosts magnitude 9+ earthquakes and is the dominant source of seismic hazard in the region (Petersen et al., 2020).The most recent great earthquake occurred in 1700, this is supported by a variety of geologic observations summarized in Walton & Staisch, et al. (2021). Across the margin there are instances of buried marsh soils and sand over mud contacts interpreted to represent sudden coseismic subsidence and deposition soon after caused by a tsunami. Microfossil analyses at these coastal marshes allow quantification of the subsidence (Figure 2). Similarly, there is evidence of deep-sea turbidite in cores, these are inferred to be deposited simultaneously at locations across the plate boundary and are thought to be seismically-induced. Historical documents in Japan reveal that a tsunami occurred in Honshu at sites as far as ∼800 km apart without an accompanying earthquake. From these documents Atwater et al. (2005) and Satake et al. ( 2003) deduced the amplitudes of the tsunami there and hydrodynamic modeling found that an earthquake in the M8.7-M9.2 range rupturing most of the CSZ could simultaneously explain the North-American and Japanese observations. Despite this knowledge, quantifying seismic hazards has been challenging because recurrence intervals are long, and no large events have occurred in the instrumental record. Geodetic locking models (Li et al., 2018;Schmalzle et al., 2014) can provide estimates of fault moment rates but do not resolve the offshore region clearly. Paleoseismic inferences of magnitudes and recurrences of events (Goldfinger et al., 2012(Goldfinger et al., , 2017Nelson et al., 2021) can be used as well. This information can be combined with suitable ground motion models (GMMs; e.g., Abrahamson et al., 2016), which estimate the severity of ground motions for an earthquake as a function of magnitude, source dimensions, distance, and local geology to make inferences about likely shaking and its probability of occurring.