Regional‐scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics‐based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three‐dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad‐scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km2) using 3D, limit‐equilibrium (LE)‐based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back‐calculated from inventoried landslides using 3D LE‐based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA. Plain Language Summary: Landslides are potentially destructive natural hazards that may impact topography, ecology, and important infrastructure. Often, previously triggered landslides are used to better understand where landslides may occur in the future, but in some regions, previous landslides may be poorly documented or characterized. In lieu of these data, models capturing the physics tied to landsliding may be paired with digital topography to predict landslide potential across landscapes. However, existing models are limited in their ability to efficiently capture the realistic geometry of three‐dimensional (3D) landslides across large landscapes. This study presents RegionGrow3D, a new model that identifies slope instabilities throughout large regions, and then grows those instabilities into landslides of previously unknown geometry by balancing forces within the underlying soil. The model is applied to a region in the Oregon Coast Range, USA to develop empirical relationships between landslide area, volume, and frequency for a range of rainfall scenarios and to determine the amount of rainfall required to trigger landslides on specific terrain features. Furthermore, because soil parameters are highly uncertain at this spatial scale (≥1,000 km2), a 3D landslide model is used to forensically assess soil strength from previously documented landslides in the area. Key Points: RegionGrow3D is a regional‐scale three‐dimensional slope stability model that grows discrete landslide volumes to achieve force equilibriumRegionGrow3D is parameterized using a range of input parameters to quantify the size, location, and likelihood of landslide initiationModeled landslides are used to quantify geomorphic scaling laws and rainfall thresholds that consider model input uncertainty [ABSTRACT FROM AUTHOR]