Post-fire Debris Flow Initiation Model
We collected post-fire data on rainfall and hydrogeomorphic outcomes from fourteen burned regions in Arizona and New Mexico. This dataset includes 163 observations from 200 watersheds in the first year after the fire, totaling 3,144 records for specific events. Rainfall intensities were averaged over 10, 15, 30, and 60 minutes. Following Staley et al. (2017), we derived features using the product of rainfall with terrain characteristics (RT), fire metrics (RF), and soil properties (RS).


We developed models using logistic regression (LR), linear discriminant analysis (LDA), random forest (RF), and XGBoost (XGB), resulting in 3,200 two-feature and 16,000 three-feature models. We chose two two-feature LR models, SWA and SWB, for their interpretability and user-friendliness, offering comparable or better performance than complex models.
The SWA model incorporates rainfall, terrain, and fire metrics based on 15-minute averaged rainfall intensities, mean slope, and the fraction of watershed area burned at different severity levels. In contrast, the SWB model uses the same rainfall and slope metrics but incorporates the mean watershed differenced normalized burn ratio (dNBR) divided by 1000.

We calculated rainfall intensity–duration (ID) thresholds for a 15-minute duration using the SWA and SWB models. For comparison, we also used the M1 model (Staley et al., 2017), which is widely utilized for assessing postfire debris-flow risks. The study included the 651 South watershed from the 2017 Pinal Fire in Arizona and Watershed D from the 2020 Tadpole Fire in New Mexico, while Watershed CON1 was withheld from the training/testing dataset and used only for calculating intensity–duration (ID) thresholds.
