Assessing locations susceptible to shallow landslide initiation during prolonged intense rainfall in the Lares, Utuado, and Naranjito municipalities of Puerto Rico

Hurricane Maria induced about 70 000 landslides throughout Puerto Rico, USA, including thousands each in three municipalities situated in Puerto Rico's rugged Cordillera Central range. By combining a nonlinear soil-depth model, presumed wettest-case pore pressures, and quasi-three-dimensional (3D) slope-stability analysis, we developed a landslide susceptibility map that has very good performance and continuous susceptibility zones having smooth, buffered boundaries. Our landslide susceptibility map enables assessment of potential ground-failure locations and their use as landslide sources in a companion assessment of inundation and debris-flow runout. The quasi-3D factor of safety, inline-formulaF₃, showed strong inverse correlation to landslide density (high density at low inline-formulaF₃). Area under the curve (AUC) of true positive rate (TPR) versus false positive rate (FPR) indicated success of inline-formulaF₃ in identifying head-scarp points (AUC inline-formula= 0.84) and source-area polygons (0.85 inline-formula≤ AUC inline-formula≤ 0.88). The susceptibility zones enclose specific percentages of observed landslides. Thus, zone boundaries use successive inline-formulaF₃ levels for increasing TPR of landslide head-scarp points, with zones bounded by inline-formulaF₃ at TPR inline-formula= 0.75, very high; inline-formulaF₃ at TPR inline-formula= 0.90, high; and the remainder moderate to low. The very high susceptibility zone, with 118 landslides kminline-formula⁻², covered 23 % of the three municipalities. The high zone (51 landslides kminline-formula⁻²) covered another 10 %.