SeanPaul M. La Selle Finn Løvholt Steven J. Gibbons Boe J. Derosier Daniel S. Brothers 20251220 text data release DOI:10.5066/P14CB2SN Pacific Coastal and Marine Science Center, Santa Cruz, California U.S. Geological Survey Suggested Citation: La Selle, S.M., Løvholt, F., Gibbons, S. J., Derosier, B.J., Brothers, D.S., 2025, Sublacustrine landslide and tsunami models from Lake Quinault, Washington: U.S. Geological Survey data release, https://doi.org/10.5066/P14CB2SN. https://doi.org/10.5066/P14CB2SN This USGS data release provides model setup files to simulate a hypothetical sublacustrine landslide in Lake Quinault and resulting tsunami waves. The BingClaw model (Kim and others, 2019; Kim and others, 2025) is used to model slope failure and GeoClaw (Berger and LeVeque, 2023; Clawpack Development Team, 2025) is used to model tsunami propagation and inundation. A topobathymetric digital terrain model (DTM) used in both models is provided and was merged from several elevation data sources (National Centers for Environmental Information, 2023; Washington Geological Survey, 2012; 2018; 2019; 2020). Zip files are provided for each model that can be used to reproduce the simulations. An example set of modeled outputs are also provided in each zip file. Lake Quinault, on the Olympic Peninsula (Washington, U.S.), is located above the Cascadia megathrust and near several mapped upper crustal faults. Past earthquakes on these faults have likely triggered slope failures in Lake Quinault (Leithold and others, 2018). The models in this data release were run to assess the tsunamigenic potential of a hypothetical sublacustrine landslide in Lake Quinault. These data are intended for science researchers, students, policy makers, and the general public. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. 20250730 20251204 ground condition at time data were collected Planned As needed -123.99420 -123.77104 47.56284 47.41273 USGS Metadata Identifier USGS:50ca5ad5-e621-427a-a91b-14f2990279f4 Data Categories for Marine Planning Bathymetry and Elevation ISO 19115 Topic Category geoscientificInformation USGS Thesaurus mathematical simulation earthquakes landslides tsunamis Marine Realms Information Bank (MRIB) keywords modeling hydrodynamics None U.S. Geological Survey USGS Coastal and Marine Hazards and Resources Program CMHRP Pacific Coastal and Marine Science Center PCMSC BingClaw GeoClaw Geographic Names Information System (GNIS) Lake Quinault Grays Harbor County Olympic Peninsula State of Washington No access constraints USGS-authored or produced data and information are in the public domain from the U.S. Government and are freely redistributable with proper metadata and source attribution. Please recognize and acknowledge the U.S. Geological Survey as the originator(s) of the dataset and in products derived from these data. This information is not intended for navigation purposes. PCMSC Science Data Coordinator U.S. Geological Survey, Pacific Coastal and Marine Science Center mailing and physical

2885 Mission Street

Santa Cruz CA 95060-5792 USA 831-427-4747 pcmsc_data@usgs.gov quinault_mccormick_preview_image.png a) Red-relief image of topobathymetric digital terrain model at McCormick Creek in Lake Quinault. b) BingClaw model output of slide thickness in meters. c) Resulting tsunami modeled in GeoClaw. PNG Ubuntu 22.04.5 LTS (GNU/Linux 5.15.153.1-microsoft-standard-WSL2x86_64) Kim, J. Løvholt, F. Issler, D. Forsberg, C.F. 2019 Kim, J., Løvholt, F., Issler, D., and Forsberg, C. F. (2019). Landslide material control on tsunami genesis - The Storegga Slide and tsunami (8,100 years BP). Journal of Geophysical Research: Oceans, 124(6), 3607-3627. https://doi.org/10.1029/2018JC014893 Kim, J. Gibbons, S.J. Løvholt, F. 2025 Kim, J., Gibbons, S. J., and Løvholt, F. (2025). BingCLAW_5.6.1 (Version v1.0.3) [Computer software]. Accessed August, 2025. https://github.com/norwegian-geotechnical-institute/BingCLAW_5.6.1 https://github.com/norwegian-geotechnical-institute/BingCLAW_5.6.1 Berger, M.J. LeVeque, R.J. 2023 Berger, M. J., and LeVeque, R. J. (2023). Implicit adaptive mesh refinement for dispersive tsunami propagation. Society for Industrial Applied Mathematics Journal on Scientific Computing, 46(2), B554-B578. https://doi.org/10.1137/23M1585210 Clawpack Development Team 2025 Clawpack Development Team (2024), Clawpack Version 5.11.00 [Computer software], Accessed August, 2025. http://www.clawpack.org. https://zenodo.org/records/13376470 Leithold, E.L. Wegmann, K. W. Bohnenstiehl, D. R. Smith, S. G. Noren, A. O'Grady, R. 2018 Leithold, E. L., Wegmann, K. W., Bohnenstiehl, D. R., Smith, S. G., Noren, A., and O'Grady, R. (2018). Slope failures within and upstream of Lake Quinault, Washington, as uneven responses to Holocene earthquakes along the Cascadia subduction zone. Quaternary Research, 89(1), 178-200. https://doi.org/10.1017/qua.2017.96 No formal attribute accuracy tests were conducted. No formal logical accuracy tests were conducted. Dataset is considered complete for the information presented, as described in the abstract. Users are advised to read the metadata for each part of this data release carefully for additional details. A formal accuracy assessment of the horizontal positional information in the data set has either not been conducted or is not applicable. A formal accuracy assessment of the vertical positional information in the data set has either not been conducted or is not applicable. National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Information (NCEI) 2023 Raster Digital Dataset online NOAA https://www.ngdc.noaa.gov/nos/D00001-D02000/D00278.html online database 202508 Date data were accessed NCEI, 2023 bathymetry data Washington Geological Survey 2012 Raster Digital Dataset online Washington Geological Survey https://lidarportal.dnr.wa.gov online database 20250730 Date data were accessed Washington Geological Survey, 2012 lidar data, originally contracted by the Quinault Indian Nation and Puget Sound LiDAR Consortium Washington Geological Survey 2018 Raster Digital Dataset online Washington Geological Survey https://lidarportal.dnr.wa.gov online database 20250730 Date data were accessed Washington Geological Survey, 2018 lidar data, originally contracted by the USDA Forest Service, Region 6 Washington Geological Survey 2019 Raster Digital Dataset online Washington Geological Survey https://lidarportal.dnr.wa.gov online database 20250730 Date data were accessed Washington Geological Survey, 2019 lidar data, originally contracted by the U.S. Geological Survey Washington Geological Survey 2020 Raster Digital Dataset online Washington Geological Survey https://lidarportal.dnr.wa.gov online database 20250730 Date data were accessed Washington Geological Survey, 2020 lidar data, originally contracted by NASA Airborne Snow Observatory. A topobathymetric digital terrain model (DTM) was generated by merging several topographic light detection and ranging (LiDAR) datasets and a multibeam bathymetric dataset. Bathymetry was collected in 2023 by the NOAA Navigation Response Team Seattle (NCEI, 2023), and topographic LiDAR data were accessed through the Washington LiDAR portal (Washington Geological Survey, 2012, 2018, 2019, and 2020). If needed, datasets were reprojected into NAD83/UTM Zone 10 N [EPSG:26910]. Topographic LiDAR data were provided in the NAVD88 vertical datum. Bathymetric soundings were collected in a vertical datum based on averaged water level data at the outlet of Lake Quinault from 2016-2021 at USGS streamgage 12039500 (NCEI, 2023). The bathymetric vertical datum was adjusted to NAVD88 by applying a vertical offset of +56.317 meters (m), based on a modeled orthometric height provided in the hydrographic survey report (NCEI, 2023). Hydroflattened "water surfaces" were removed from LiDAR datasets overlapping Lake Quinault, resulting in a ring of no data values between the bathymetric and topographic data, representing shallow water depths where the multibeam survey could not reach. A bilinear interpolation was used to fill in the missing data using the SciPy package in Python. Artifacts from interpolation are apparent in the merged DTM, especially along the broad shallows of the northeastern Lake Quinault shoreline. Merged datasets were regridded to a final grid resolution of 2.74 m, based on the resolution of the coarsest LiDAR dataset used (Washington Geological Survey, 2020). 20250730 A portion of the topobathymetric DTM was extracted in the vicinity of McCormick Creek along the northern shore of Lake Quinault in order to run the BingClaw landslide model. The geomorphic expression of slide runout from a previous slope failure are visible in the bathymetric data. A DTM representing a possible pre-slide lake-bed surface was generated by tracing the footprint of the observed runout debris, removing this region from the topobathymetric DTM, and performing a spline interpolation to fill in the gap. A minimum volume for the observed runout debris of 210,000 cubic meters was estimated by differencing the smoothed pre-slide DTM with the modern DTM. A raster of initial slide thickness was created by tracing the upslope region of the observed runout and filling this region with an initial guess of 5 m of sediment that tapers to 0 m thickness along the edges of the traced region. Including the modern delta, the initial modeled slide volume provided in the BingClaw example is 260,000 cubic meters. 20251117 A suite of BingClaw models was run with the modeled slide thickness inputs, varying rheologic input parameters until the modeled runout roughly matched the observed perimeter of runout debris. In the example provided in this data release, an initial yield stress of 5,000 Pascals, residual yield stress of 500 Pascals, and remolding coefficient (gamma) of 0.005 was used. The model was run using a spatial resolution of 10 m, and model results were written to output files every 2 seconds over 120 seconds of simulation time. The resulting model outputs of evolving slide thickness were used to force the GeoClaw tsunami model. BingClaw outputs were converted to GeoClaw inputs using a python script. Simple plots of slide thickness were generated using python scripts (setplot.py) provided with the BingClaw model. 20251118 GeoClaw models were run using the BingClaw model outputs as initial water level conditions based on the method of Kim and others (2019). For each output timestep in the BingClaw model, the change in water depth due to the change in landslide thickness is used as a change in the GeoClaw water surface elevations, assuming hydrostatic pressure. This coupling scheme does not directly account for longitudinal momentum transfer from the slide to the generated tsunami waves. Two types of hydrodynamic tsunami propagation models were run: a non-dispersive wave model using the nonlinear shallow water equations (swe) and a dispersive wave model using the Serre-Green-Naghdi Boussinesq-type equations (bouss). Both models use the shallow water equations for inundation and in modeled water depths less than 5 m. Both models were run with adaptive mesh refinement (AMR) at refinement level resolutions of 1000, 200, 40, 10, and 5 m. The parameter for sea level in GeoClaw was set to 57 m to represent the lake water surface elevation. Model figures and animations were plotted using the included "tsunami_plots.ipynb" Python Jupyter Notebooks. 20251118 Files needed to run BingClaw and GeoClaw models along with example modeled outputs, were compiled into compressed archives for distribution. A supplemental information file "readme.md " describes in greater detail the file structure of the included compressed archived, and limited instructions on usage. 20251212 Horizontal coordinates are in the UTM Zone 10 N coordinate system for all arc-ascii format model files (.asc and .tt3). Elevation values in these format are adjusted to the NAVD88 vertical datum, in meters. The raster properties below are for the topobathymetric DTM Raster Grid Cell 6014 6200 1 Universal Transverse Mercator 10 0.9996 -123.0 0.0 500000.0 0.0 row and column 2.74 2.74 Meter D_North_American_1983 GRS_1980 6378137.0 298.257222101 North American Vertical Datum of 1988 0.001 meters Explicit elevation coordinate included with horizontal coordinates Compressed .zip archives contain files of various types and formats. Model input files compatible with BingClaw v5.6.1 and GeoClaw v5.11.0 are provided in the zip archives "quinault_landslide_model_bingclaw.zip" and "quinault_tsunami_model_geoclaw.zip". An example set of modeled outputs are also provided in each zip file. For additional information, see the File Structure section of the supplemental file "readme.md". Source code for BingClaw and GeoClaw can be installed using instructions in Kim and others, 2025 and Clawpack Development Team, 2025. The entity and attribute information were generated by the individual and/or agency identified as the originator of the data set. Please review the rest of the metadata record for additional details and information. U.S. Geological Survey - CMGDS Mailing and Physical Address

2885 Mission Street

Santa Cruz CA 95060 USA 831-427-4747 pcmsc_data@usgs.gov Compressed .zip archives contain files of various types and formats. Model input files compatible with BingClaw v5.6.1 and GeoClaw v5.11.0 are provided in the zip archives "quinault_landslide_model_bingclaw.zip" and "quinault_tsunami_model_geoclaw.zip". Source code for BingClaw and GeoClaw can be installed using instructions in Kim and others, 2025 and Clawpack Development Team, 2025 Unless otherwise stated, all data, metadata and related materials are considered to satisfy the quality standards relative to the purpose for which the data were collected. Although these data and associated metadata have been reviewed for accuracy and completeness and approved for release by the U.S. Geological Survey (USGS), no warranty expressed or implied is made regarding the display or utility of the data on any other system or for general or scientific purposes, nor shall the act of distribution constitute any such warranty. ASCII, plain text BingClaw v5.6.1 and GeoClaw v5.11.0 zipped file folders containing model input and output files for BingClaw v5.6.1 (quinault_landslide_model_bingclaw.zip) and GeoClaw v5.11.0 (quinault_tsunami_model_geoclaw.zip). Each zip file contains a folder named "mccormick_slide". Within the "mccormick_slide" folders are relevant model input files such as topobathymetric DTMs and slide thickness. One example run is provided for BingClaw, in the folder "run01_i5000_r500_g005". Two example tsunami model runs using non-dispersive (swe) and dispersive (bouss) equations are provided for GeoClaw along with modeled outputs. See supplemental files "readme.md" or "readme.pdf" for more details. zip 2202.93 https://doi.org/10.5066/P14CB2SN Data may be accessed and downloaded via the Internet by using the Network_Resource_Name link and scrolling to the Simulation Data section. None 20251220 PCMSC Science Data Coordinator U.S. Geological Survey, Pacific Coastal and Marine Science Center mailing and physical

2885 Mission Street

Santa Cruz CA 95060 831-427-4747 pcmsc_data@usgs.gov Content Standard for Digital Geospatial Metadata FGDC-STD-001-1998