Metadata: Identification_Information: Citation: Citation_Information: Originator: SeanPaul M. La Selle Originator: Finn Løvholt Originator: Steven J. Gibbons Originator: Boe J. Derosier Originator: Daniel S. Brothers Publication_Date: 20251220 Title: Sublacustrine landslide and tsunami models from Lake Quinault, Washington Geospatial_Data_Presentation_Form: text Series_Information: Series_Name: data release Issue_Identification: DOI:10.5066/P14CB2SN Publication_Information: Publication_Place: Pacific Coastal and Marine Science Center, Santa Cruz, California Publisher: U.S. Geological Survey Other_Citation_Details: 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. Online_Linkage: https://doi.org/10.5066/P14CB2SN Description: Abstract: 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. Purpose: 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. Supplemental_Information: Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Time_Period_of_Content: Time_Period_Information: Range_of_Dates/Times: Beginning_Date: 20250730 Ending_Date: 20251204 Currentness_Reference: ground condition at time data were collected Status: Progress: Planned Maintenance_and_Update_Frequency: As needed Spatial_Domain: Bounding_Coordinates: West_Bounding_Coordinate: -123.99420 East_Bounding_Coordinate: -123.77104 North_Bounding_Coordinate: 47.56284 South_Bounding_Coordinate: 47.41273 Keywords: Theme: Theme_Keyword_Thesaurus: USGS Metadata Identifier Theme_Keyword: USGS:50ca5ad5-e621-427a-a91b-14f2990279f4 Theme: Theme_Keyword_Thesaurus: Data Categories for Marine Planning Theme_Keyword: Bathymetry and Elevation Theme: Theme_Keyword_Thesaurus: ISO 19115 Topic Category Theme_Keyword: geoscientificInformation Theme: Theme_Keyword_Thesaurus: USGS Thesaurus Theme_Keyword: mathematical simulation Theme_Keyword: earthquakes Theme_Keyword: landslides Theme_Keyword: tsunamis Theme: Theme_Keyword_Thesaurus: Marine Realms Information Bank (MRIB) keywords Theme_Keyword: modeling Theme_Keyword: hydrodynamics Theme: Theme_Keyword_Thesaurus: None Theme_Keyword: U.S. Geological Survey Theme_Keyword: USGS Theme_Keyword: Coastal and Marine Hazards and Resources Program Theme_Keyword: CMHRP Theme_Keyword: Pacific Coastal and Marine Science Center Theme_Keyword: PCMSC Theme_Keyword: BingClaw Theme_Keyword: GeoClaw Place: Place_Keyword_Thesaurus: Geographic Names Information System (GNIS) Place_Keyword: Lake Quinault Place_Keyword: Grays Harbor County Place_Keyword: Olympic Peninsula Place_Keyword: State of Washington Access_Constraints: No access constraints Use_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. Point_of_Contact: Contact_Information: Contact_Person_Primary: Contact_Person: PCMSC Science Data Coordinator Contact_Organization: U.S. Geological Survey, Pacific Coastal and Marine Science Center Contact_Address: Address_Type: mailing and physical Address: 2885 Mission Street City: Santa Cruz State_or_Province: CA Postal_Code: 95060-5792 Country: USA Contact_Voice_Telephone: 831-427-4747 Contact_Electronic_Mail_Address: pcmsc_data@usgs.gov Browse_Graphic: Browse_Graphic_File_Name: quinault_mccormick_preview_image.png Browse_Graphic_File_Description: 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. Browse_Graphic_File_Type: PNG Native_Data_Set_Environment: Ubuntu 22.04.5 LTS (GNU/Linux 5.15.153.1-microsoft-standard-WSL2x86_64) Cross_Reference: Citation_Information: Originator: Kim, J. Originator: Løvholt, F. Originator: Issler, D. Originator: Forsberg, C.F. Publication_Date: 2019 Title: Landslide material control on tsunami genesis - The Storegga Slide and tsunami (8,100 years BP). Other_Citation_Details: 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. Online_Linkage: https://doi.org/10.1029/2018JC014893 Cross_Reference: Citation_Information: Originator: Kim, J. Originator: Gibbons, S.J. Originator: Løvholt, F. Publication_Date: 2025 Title: BingClaw_5.6.1 (Version v1.0.3) [Computer software] Other_Citation_Details: 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 Online_Linkage: https://github.com/norwegian-geotechnical-institute/BingCLAW_5.6.1 Cross_Reference: Citation_Information: Originator: Berger, M.J. Originator: LeVeque, R.J. Publication_Date: 2023 Title: Implicit adaptive mesh refinement for dispersive tsunami propagation Other_Citation_Details: 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. Online_Linkage: https://doi.org/10.1137/23M1585210 Cross_Reference: Citation_Information: Originator: Clawpack Development Team Publication_Date: 2025 Title: Clawpack Version 5.12.00 [Computer software] Other_Citation_Details: Clawpack Development Team (2024), Clawpack Version 5.11.00 [Computer software], Accessed August, 2025. http://www.clawpack.org. Online_Linkage: https://zenodo.org/records/13376470 Cross_Reference: Citation_Information: Originator: Leithold, E.L. Originator: Wegmann, K. W. Originator: Bohnenstiehl, D. R. Originator: Smith, S. G. Originator: Noren, A. Originator: O'Grady, R. Publication_Date: 2018 Title: Slope failures within and upstream of Lake Quinault, Washington, as uneven responses to Holocene earthquakes along the Cascadia subduction zone. Other_Citation_Details: 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. Online_Linkage: https://doi.org/10.1017/qua.2017.96 Data_Quality_Information: Attribute_Accuracy: Attribute_Accuracy_Report: No formal attribute accuracy tests were conducted. Logical_Consistency_Report: No formal logical accuracy tests were conducted. Completeness_Report: 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. Positional_Accuracy: Horizontal_Positional_Accuracy: Horizontal_Positional_Accuracy_Report: A formal accuracy assessment of the horizontal positional information in the data set has either not been conducted or is not applicable. Vertical_Positional_Accuracy: Vertical_Positional_Accuracy_Report: A formal accuracy assessment of the vertical positional information in the data set has either not been conducted or is not applicable. Lineage: Source_Information: Source_Citation: Citation_Information: Originator: National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Information (NCEI) Publication_Date: 2023 Title: Descriptive Report for D00278 Geospatial_Data_Presentation_Form: Raster Digital Dataset Publication_Information: Publication_Place: online Publisher: NOAA Online_Linkage: https://www.ngdc.noaa.gov/nos/D00001-D02000/D00278.html Type_of_Source_Media: online database Source_Time_Period_of_Content: Time_Period_Information: Single_Date/Time: Calendar_Date: 202508 Source_Currentness_Reference: Date data were accessed Source_Citation_Abbreviation: NCEI, 2023 Source_Contribution: bathymetry data Source_Information: Source_Citation: Citation_Information: Originator: Washington Geological Survey Publication_Date: 2012 Title: Quinault River Basin 2011 project [lidar data] Geospatial_Data_Presentation_Form: Raster Digital Dataset Publication_Information: Publication_Place: online Publisher: Washington Geological Survey Online_Linkage: https://lidarportal.dnr.wa.gov Type_of_Source_Media: online database Source_Time_Period_of_Content: Time_Period_Information: Single_Date/Time: Calendar_Date: 20250730 Source_Currentness_Reference: Date data were accessed Source_Citation_Abbreviation: Washington Geological Survey, 2012 Source_Contribution: lidar data, originally contracted by the Quinault Indian Nation and Puget Sound LiDAR Consortium Source_Information: Source_Citation: Citation_Information: Originator: Washington Geological Survey Publication_Date: 2018 Title: Olympic National Forest 2017 project [lidar data] Geospatial_Data_Presentation_Form: Raster Digital Dataset Publication_Information: Publication_Place: online Publisher: Washington Geological Survey Online_Linkage: https://lidarportal.dnr.wa.gov Type_of_Source_Media: online database Source_Time_Period_of_Content: Time_Period_Information: Single_Date/Time: Calendar_Date: 20250730 Source_Currentness_Reference: Date data were accessed Source_Citation_Abbreviation: Washington Geological Survey, 2018 Source_Contribution: lidar data, originally contracted by the USDA Forest Service, Region 6 Source_Information: Source_Citation: Citation_Information: Originator: Washington Geological Survey Publication_Date: 2019 Title: Olympic Peninsula, Washington 3DEP 2019 project [lidar data] Geospatial_Data_Presentation_Form: Raster Digital Dataset Publication_Information: Publication_Place: online Publisher: Washington Geological Survey Online_Linkage: https://lidarportal.dnr.wa.gov Type_of_Source_Media: online database Source_Time_Period_of_Content: Time_Period_Information: Single_Date/Time: Calendar_Date: 20250730 Source_Currentness_Reference: Date data were accessed Source_Citation_Abbreviation: Washington Geological Survey, 2019 Source_Contribution: lidar data, originally contracted by the U.S. Geological Survey Source_Information: Source_Citation: Citation_Information: Originator: Washington Geological Survey Publication_Date: 2020 Title: Olympic Park 2014 project [lidar data] Geospatial_Data_Presentation_Form: Raster Digital Dataset Publication_Information: Publication_Place: online Publisher: Washington Geological Survey Online_Linkage: https://lidarportal.dnr.wa.gov Type_of_Source_Media: online database Source_Time_Period_of_Content: Time_Period_Information: Single_Date/Time: Calendar_Date: 20250730 Source_Currentness_Reference: Date data were accessed Source_Citation_Abbreviation: Washington Geological Survey, 2020 Source_Contribution: lidar data, originally contracted by NASA Airborne Snow Observatory. Process_Step: Process_Description: 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). Process_Date: 20250730 Process_Step: Process_Description: 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. Process_Date: 20251117 Process_Step: Process_Description: 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. Process_Date: 20251118 Process_Step: Process_Description: 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. Process_Date: 20251118 Process_Step: Process_Description: 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. Process_Date: 20251212 Spatial_Data_Organization_Information: Indirect_Spatial_Reference: 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 Direct_Spatial_Reference_Method: Raster Raster_Object_Information: Raster_Object_Type: Grid Cell Row_Count: 6014 Column_Count: 6200 Vertical_Count: 1 Spatial_Reference_Information: Horizontal_Coordinate_System_Definition: Planar: Grid_Coordinate_System: Grid_Coordinate_System_Name: Universal Transverse Mercator Universal_Transverse_Mercator: UTM_Zone_Number: 10 Transverse_Mercator: Scale_Factor_at_Central_Meridian: 0.9996 Longitude_of_Central_Meridian: -123.0 Latitude_of_Projection_Origin: 0.0 False_Easting: 500000.0 False_Northing: 0.0 Planar_Coordinate_Information: Planar_Coordinate_Encoding_Method: row and column Coordinate_Representation: Abscissa_Resolution: 2.74 Ordinate_Resolution: 2.74 Planar_Distance_Units: Meter Geodetic_Model: Horizontal_Datum_Name: D_North_American_1983 Ellipsoid_Name: GRS_1980 Semi-major_Axis: 6378137.0 Denominator_of_Flattening_Ratio: 298.257222101 Vertical_Coordinate_System_Definition: Altitude_System_Definition: Altitude_Datum_Name: North American Vertical Datum of 1988 Altitude_Resolution: 0.001 Altitude_Distance_Units: meters Altitude_Encoding_Method: Explicit elevation coordinate included with horizontal coordinates Entity_and_Attribute_Information: Overview_Description: Entity_and_Attribute_Overview: 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. Entity_and_Attribute_Detail_Citation: 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. Distribution_Information: Distributor: Contact_Information: Contact_Organization_Primary: Contact_Organization: U.S. Geological Survey - CMGDS Contact_Address: Address_Type: Mailing and Physical Address Address: 2885 Mission Street City: Santa Cruz State_or_Province: CA Postal_Code: 95060 Country: USA Contact_Voice_Telephone: 831-427-4747 Contact_Electronic_Mail_Address: pcmsc_data@usgs.gov Resource_Description: 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 Distribution_Liability: 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. Standard_Order_Process: Digital_Form: Digital_Transfer_Information: Format_Name: ASCII, plain text Format_Version_Number: BingClaw v5.6.1 and GeoClaw v5.11.0 Format_Information_Content: 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. File_Decompression_Technique: zip Transfer_Size: 2202.93 Digital_Transfer_Option: Online_Option: Computer_Contact_Information: Network_Address: Network_Resource_Name: https://doi.org/10.5066/P14CB2SN Access_Instructions: Data may be accessed and downloaded via the Internet by using the Network_Resource_Name link and scrolling to the Simulation Data section. Fees: None Metadata_Reference_Information: Metadata_Date: 20251220 Metadata_Contact: Contact_Information: Contact_Person_Primary: Contact_Person: PCMSC Science Data Coordinator Contact_Organization: U.S. Geological Survey, Pacific Coastal and Marine Science Center Contact_Address: Address_Type: mailing and physical Address: 2885 Mission Street City: Santa Cruz State_or_Province: CA Postal_Code: 95060 Contact_Voice_Telephone: 831-427-4747 Contact_Electronic_Mail_Address: pcmsc_data@usgs.gov Metadata_Standard_Name: Content Standard for Digital Geospatial Metadata Metadata_Standard_Version: FGDC-STD-001-1998