Landslide susceptibility models show the potential of landslide occurrence at a location. These models are pivotal for reducing losses associated with landslides (Godt et al., 2022). In this data release, we include susceptibility results from the associated manuscript by Woodard and Mirus (2025). This manuscript shows how a morphometric model can create consistent and effective susceptibility models over large regions (> 100 km2) by analyzing the terrain’s topography. The model assumes that areas with high relative slope and hillslope area in comparison to the rest of the terrain are more susceptible to landsliding. As the model’s only input is elevation data, it mitigates the data biases common in the data-driven statistical methods (e.g., machine learning) generally used over these scales. We compare the morphometric model outputs to a parsimonious national susceptibility map and logistic regression machine learning models. The national susceptibility map is available in Belair et al., (2024). The two logistic regression models are trained on the landslide data available in the Willamette Valley Hydrologic Unit Code (HUC) 4 watershed (DOGAMI, 2024). To account for the effects of the sampling ratio of event to non-event data points, we create two logistic regression models. The first uses a 1:1 sampling ratio of landslide to non-landslide points and the second uses all the data within the training data which results in a 1:33 sampling ratio. Environmental datasets requisite for the logistic regression models are all derived from the three-dimensional elevation program (3DEP) (U.S. Geological Survey, 2019a) preprocessed within the National Hydrography Dataset (U.S. Geological Survey, 2019b). The morphometric model was derived using only the 3DEP dataset without any input of where landslides have occurred. All model outputs are shown with slope units. This data release includes the following files: 1) logistic regression results with 1:1 sampling ratio over Willamette Valley HUC4 watershed (1709) (Logistic_1709_1.zip); 2) logistic regression results with 1:33 sampling ratio over Willamette Valley HUC4 watershed (1709) (Logistic_1709_All.zip); 3) morphometric results with uniform weights over the Willamette Valley HUC4 watershed (1709) (Morph_Uniform_1709.zip); 4) morphometric results with area weights over the 1701 HUC 4 watershed (Morph_Area_1701.zip); 5) morphometric results with area weights over the 1702 HUC 4 watershed (Morph_Area_1702.zip); 6) morphometric results with area weights over the 1703 HUC 4 watershed (Morph_Area_1703.zip); 7) morphometric results with area weights over the 1704 HUC 4 watershed (Morph_Area_1704.zip); 8) morphometric results with area weights over the 1705 HUC 4 watershed (Morph_Area_1705.zip); 9) morphometric results with area weights over the 1706 HUC 4 watershed (Morph_Area_1706.zip); 10) morphometric results with area weights over the 1707 HUC 4 watershed (Morph_Area_1707.zip); 11) morphometric results with area weights over the 1708 HUC 4 watershed (Morph_Area_1708.zip); 12) morphometric results with area weights over the 1709 HUC 4 watershed (Morph_Area_1709.zip); 13) morphometric results with area weights over the 1710 HUC 4 watershed (Morph_Area_1710.zip); 14) morphometric results with area weights over the 1711 HUC 4 watershed (Morph_Area_1711.zip); 15) morphometric results with area weights over the 1712 HUC 4 watershed (Morph_Area_1712.zip). 16) shape file field descriptors (Field_Descriptors.txt) Each zip-file contains the vector shapefiles of interest which can be extracted using most archiver software. References Cited DOGAMI. (2024). SLIDO (Version 4.5) [Data set]. https://pubs.oregon.gov/dogami/SLIDO/4.5/SLIDO_Release_4p5_wMetadata.gdb.zip. Gina M Belair, Jeanne M Jones, Sabrina N Martinez, Benjamin B Mirus, & Nathan J Wood. (2024). Slope-Relief Threshold Landslide Susceptibility Models for the United States and Puerto Rico [Data Release]. U.S. Geological Survey.

Landslide susceptibility models show the potential of landslide occurrence at a location. These models are pivotal for reducing losses associated with landslides (Godt et al., 2022). In this data release, we include susceptibility results from the associated manuscript by Woodard and Mirus (2025). This manuscript shows how a morphometric model can create consistent and effective susceptibility models over large regions (> 100 km2) by analyzing the terrain’s topography. The model assumes that areas with high relative slope and hillslope area in comparison to the rest of the terrain are more susceptible to landsliding. As the model’s only input is elevation data, it mitigates the data biases common in the data-driven statistical methods (e.g., machine learning) generally used over these scales. We compare the morphometric model outputs to a parsimonious national susceptibility map and logistic regression machine learning models. The national susceptibility map is available in Belair et al., (2024). The two logistic regression models are trained on the landslide data available in the Willamette Valley Hydrologic Unit Code (HUC) 4 watershed (DOGAMI, 2024). To account for the effects of the sampling ratio of event to non-event data points, we create two logistic regression models. The first uses a 1:1 sampling ratio of landslide to non-landslide points and the second uses all the data within the training data which results in a 1:33 sampling ratio. Environmental datasets requisite for the logistic regression models are all derived from the three-dimensional elevation program (3DEP) (U.S. Geological Survey, 2019a) preprocessed within the National Hydrography Dataset (U.S. Geological Survey, 2019b). The morphometric model was derived using only the 3DEP dataset without any input of where landslides have occurred. All model outputs are shown with slope units. This data release includes the following files: 1) logistic regression results with 1:1 sampling ratio over Willamette Valley HUC4 watershed (1709) (Logistic_1709_1.zip); 2) logistic regression results with 1:33 sampling ratio over Willamette Valley HUC4 watershed (1709) (Logistic_1709_All.zip); 3) morphometric results with uniform weights over the Willamette Valley HUC4 watershed (1709) (Morph_Uniform_1709.zip); 4) morphometric results with area weights over the 1701 HUC 4 watershed (Morph_Area_1701.zip); 5) morphometric results with area weights over the 1702 HUC 4 watershed (Morph_Area_1702.zip); 6) morphometric results with area weights over the 1703 HUC 4 watershed (Morph_Area_1703.zip); 7) morphometric results with area weights over the 1704 HUC 4 watershed (Morph_Area_1704.zip); 8) morphometric results with area weights over the 1705 HUC 4 watershed (Morph_Area_1705.zip); 9) morphometric results with area weights over the 1706 HUC 4 watershed (Morph_Area_1706.zip); 10) morphometric results with area weights over the 1707 HUC 4 watershed (Morph_Area_1707.zip); 11) morphometric results with area weights over the 1708 HUC 4 watershed (Morph_Area_1708.zip); 12) morphometric results with area weights over the 1709 HUC 4 watershed (Morph_Area_1709.zip); 13) morphometric results with area weights over the 1710 HUC 4 watershed (Morph_Area_1710.zip); 14) morphometric results with area weights over the 1711 HUC 4 watershed (Morph_Area_1711.zip); 15) morphometric results with area weights over the 1712 HUC 4 watershed (Morph_Area_1712.zip). 16) shape file field descriptors (Field_Descriptors.txt) Each zip-file contains the vector shapefiles of interest which can be extracted using most archiver software. References Cited DOGAMI. (2024). SLIDO (Version 4.5) [Data set]. https://pubs.oregon.gov/dogami/SLIDO/4.5/SLIDO_Release_4p5_wMetadata.gdb.zip. Gina M Belair, Jeanne M Jones, Sabrina N Martinez, Benjamin B Mirus, & Nathan J Wood. (2024). Slope-Relief Threshold Landslide Susceptibility Models for the United States and Puerto Rico [Data Release]. U.S. Geological Survey.

The Susceptibility to Deep-Seated Landslides map covers the entire state of California and was originally published in May of 2011 as CGS Map Sheet 58. It made use of several data layers of varying scales and formats, such as Landslide Inventory, Geology, Rock Strength, and Slope. For the statewide analysis of landslide susceptibility, the methodology of Wilson and Keefer (1985) was used in combining the rock strength and slope data layers as implemented by Ponti, el al. (2008) to create classes of landslide susceptibility (0 to 10, low to high). These classes express the generalization that on very low slopes, landslide susceptibility is low even in weak materials, and that landslide susceptibility increases with slope and in weak rocks.

,

Earthquake-triggered landslides can significantly contribute to human and economic losses during and immediately following earthquakes, but data on the runout behavior of such ground failures is limited. Hazard assessment of coseismic landslide risk can vary dramatically depending on landslide mobility and runout extent, which makes modeling of such behavior imperative. Predictive and empirical models require comprehensive datasets with diverse climatic, topographic, and geologic factors. We present an openly accessible global dataset of coseismic landslide runout lengths, produced from an automated method for estimating runout length from existing landslide inventories. This tool was developed and validated using manually measured runout lengths of 1,726 landslides from five global earthquake-induced landslide inventories spanning a variety of terrains and geologic settings. The resultant database contains 73,665 measured and estimated runout lengths of coseismic landslides from 23 global earthquakes derived from the USGS’s open repository of earthquake-triggered ground-failure inventories on ScienceBase (https://doi.org/10.5066/F7H70DB4, v4.0, Schmitt et al., 2022). We present separate data files for each inventory, reporting area and predicted or measured runout lengths of individual landslides.

Landslide susceptibility maps are essential tools in infrastructure planning, hazard mitigation, and risk reduction. Susceptibility maps trained in one area have been found to be unreliable when applied to different areas (Woodard et al., 2023). This limitation leads to the need for a national map that is higher resolution and rigorous, but simple enough to be applied to diverse terrains and landslide types. The susceptibility maps presented here cover the conterminous United States (CONUS), Alaska (AK), Hawaii (HI), and Puerto Rico (PR) with a resolution of 90-m. Other United States (U.S.) territories were not considered due to insufficient landslide and digital elevation data. We also provide information on the proportion of susceptible terrain as well as the density (landslides per square kilometer) of documented landslides within susceptible terrain for each U.S. county. To generate the susceptibility maps we used 1/3 arc-second digital elevation models (DEMs) (U.S. Geological Survey, 2019) to calculate slope and 100-m relief, 613,724 unique landslides from our national landslide inventory compilation (Belair et al., 2022) to train the models and compute U.S. county aggregated susceptibility information, and high-performance computing resources to train the models (Falgout and Gordon, 2023). We present two slope-relief threshold models: (1) a linear regression model weighted by landslide density of each ecoregion (Wiken et al., 2011), and (2) a quantile nonlinear regression model fitted to the 10th quantile of the data. We (1) removed extraneous landslide data, (2) averaged 50 model runs, and then (3) down-sampled the maps from 10-m to 90-m resolution to account for uncertainty in the DEM and landslide position. The nonlinear model (n10) performs better under most topographic conditions and optimally balances our priorities of capturing observed landslides (98.9%) while minimizing area covered by susceptible terrain (44.6%). The weighted linear model (lw) captures slightly fewer landslides (98.8%) and has slightly less susceptible terrain (43.1%). The values of both maps represent the number of susceptible 10-m cells within each 90-m cell after down-sampling and can range from 0 to 81. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration (or cell value) capture the majority of landslides, thus representing higher susceptibility areas. The susceptibility maps were then used to determine the total area of landslide susceptible terrain (square kilometers) for each U.S. county. The national landslide inventory compilation was used to determine the number of documented landslides within susceptible terrain for each county. This information was then used to calculate the proportion of susceptible terrain and the density of documented landslides within susceptible terrain for each county in the United States. This information is provided in tabular format, with columns corresponding to the information discussed above, and each row corresponding to a U.S. county. Further information about this analysis can be found in an interpretive publication (Mirus et al., 2024). This data release includes: (1) weighted linear susceptibility maps (lw_susc.zip), (2) quantile nonlinear susceptibility maps (n10_susc.zip), (3) landslide data used to develop the models (landslides.csv), (4) county aggregated susceptibility information (county_analysis.csv), (5) readme and analysis files, and (6) metadata. References Cited Belair, G. M., Jones, E. S., Slaughter, S. L., and Mirus, B. B., 2022, Landslide Inventories across the United States version 2: U.S. Geological Survey data release, https://doi.org/10.5066/P9FZUX6N Falgout, J. T., and Gordon, J., 2023, USGS Advanced Research Computing, USGS Yeti Supercomputer: U.S. Geological Survey, https://doi.org/10.5066/F7D798MJ Mirus, B. B., Belair, G. M., Wood, N. J., Jones, J. M., and Martinez, S. M., 2024, Parsimonious high-resolution

Landslide susceptibility maps are essential tools in infrastructure planning, hazard mitigation, and risk reduction. Susceptibility maps trained in one area have been found to be unreliable when applied to different areas (Woodard et al., 2023). This limitation leads to the need for a national map that is higher resolution and rigorous, but simple enough to be applied to diverse terrains and landslide types. The susceptibility maps presented here cover the conterminous United States (CONUS), Alaska (AK), Hawaii (HI), and Puerto Rico (PR) with a resolution of 90-m. Other United States (U.S.) territories were not considered due to insufficient landslide and digital elevation data. We also provide information on the proportion of susceptible terrain as well as the density (landslides per square kilometer) of documented landslides within susceptible terrain for each U.S. county. To generate the susceptibility maps we used 1/3 arc-second digital elevation models (DEMs) (U.S. Geological Survey, 2019) to calculate slope and 100-m relief, 613,724 unique landslides from our national landslide inventory compilation (Belair et al., 2022) to train the models and compute U.S. county aggregated susceptibility information, and high-performance computing resources to train the models (Falgout and Gordon, 2023). We present two slope-relief threshold models: (1) a linear regression model weighted by landslide density of each ecoregion (Wiken et al., 2011), and (2) a quantile nonlinear regression model fitted to the 10th quantile of the data. We (1) removed extraneous landslide data, (2) averaged 50 model runs, and then (3) down-sampled the maps from 10-m to 90-m resolution to account for uncertainty in the DEM and landslide position. The nonlinear model (n10) performs better under most topographic conditions and optimally balances our priorities of capturing observed landslides (98.9%) while minimizing area covered by susceptible terrain (44.6%). The weighted linear model (lw) captures slightly fewer landslides (98.8%) and has slightly less susceptible terrain (43.1%). The values of both maps represent the number of susceptible 10-m cells within each 90-m cell after down-sampling and can range from 0 to 81. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration (or cell value) capture the majority of landslides, thus representing higher susceptibility areas. The susceptibility maps were then used to determine the total area of landslide susceptible terrain (square kilometers) for each U.S. county. The national landslide inventory compilation was used to determine the number of documented landslides within susceptible terrain for each county. This information was then used to calculate the proportion of susceptible terrain and the density of documented landslides within susceptible terrain for each county in the United States. This information is provided in tabular format, with columns corresponding to the information discussed above, and each row corresponding to a U.S. county. Further information about this analysis can be found in an interpretive publication (Mirus et al., 2024). This data release includes: (1) weighted linear susceptibility maps (lw_susc.zip), (2) quantile nonlinear susceptibility maps (n10_susc.zip), (3) landslide data used to develop the models (landslides.csv), (4) county aggregated susceptibility information (county_analysis.csv), (5) readme and analysis files, and (6) metadata. References Cited Belair, G. M., Jones, E. S., Slaughter, S. L., and Mirus, B. B., 2022, Landslide Inventories across the United States version 2: U.S. Geological Survey data release, https://doi.org/10.5066/P9FZUX6N Falgout, J. T., and Gordon, J., 2023, USGS Advanced Research Computing, USGS Yeti Supercomputer: U.S. Geological Survey, https://doi.org/10.5066/F7D798MJ Mirus, B. B., Belair, G. M., Wood, N. J., Jones, J. M., and Martinez, S. M., 2024, Parsimonious high-resolution