The scientific background information for Emergency Assessment of Post-Fire Debris-Flow Hazards.

The preliminary hazard assessment relies upon empirical models to estimate the likelihood and volume of debris flows for selected basins in response to a design storm. Beginning in 2016, the empirical models are based upon historical debris–flow occurrence and magnitude data, rainfall storm conditions, terrain and soils information, and burn–severity data from recently burned areas (Staley et al., 2016, 2017). Hazard assessments for fires prior to May 2016 relied upon an earlier set of empirical methods described here.

The post-fire debris-flow likelihood, volume, and combined hazards are estimated at both the drainage-basin scale and in a spatially distributed manner along the drainage network within each basin. The characteristics of basins affected by the fire were calculated using a geographic information system (GIS). Debris-flow likelihood and volume were estimated for each basin outlet as well as along the upstream drainage networks (pixels where the contributing area is greater than or equal to 0.02 km^{2}). Independent variable values were calculated for each pixel along the drainage network and summarized at the stream segment scale to obtain estimates of debris-flow likelihood and volume.

Required input data include a differenced normalized burn ratio image (dNBR), and field-validated soil burn severity data. If field-validated soil burn severity data are unavailable, it is possible to use a Burned Area Reflectance Classification (BARC4) image in substitution. Please note that the use of BARC4 data increases the degree of uncertainty in the model estimates.

**Likelihood Model**

The likelihood of a debris-flow in response to a given peak 15-minute rainfall intensity are based upon a logistic regression approach, which combines the following two equations:

(1) P = e^{x} / (1 + e^{x})

**Where**

- P is the probability of debris-flow occurrence in fractional form, and
- e
^{x}is the exponential function where e represents the mathematical constant 2.718.

For recently burned areas in the western United States, equation 2 is used to calculate x:

(2) x = -3.63 + (0.41 × X1R) + (0.67 × X2R) + (0.7 × X3R)

**Where**

- X1R is the proportion of upslope area classified as high or moderate soil burn severity and with gradients ≥ 23°, multiplied by the peak 15-minute rainfall accumulation of the design storm (in millimeters [mm]),
- X2R is the average differenced normalized burn ratio (dNBR) of the upslope area, multiplied by the peak 15-minute rainfall accumulation of the design storm (in millimeters [mm]),
- X3R is the soil KF-Factor (Schwartz and Alexander, 1995) of the upslope area, multiplied by the peak 15-minute rainfall accumulation of the design storm (in millimeters [mm]).

Likelihood values predicted by the equation potentially range from 0 (least likely) to 1 (most likely). The predicted likelihood values are assigned to 1 of 5 equal interval classes for cartographic display and are represented as a percentage likelihood (rather than a ratio).

**Volume Model**

Debris-flow volumes both at the basin outlet and along the drainage network are predicted using a multiple linear regression model (Gartner and others, 2014). The multiple linear regression models are used to estimate the volume (V, in m3) of material that could issue from a point along the drainage network in response to a storm of a given rainfall intensity.

Potential debris-flow volume is calculated with equation 4:

(3) ln(V) = 4.22 + (0.13 × sqrt(ElevRange)) + (0.36 × ln(HM_{km})) + (0.39 × sqrt(i15))

**Where**

- ElevRange is the range (maximum elevation–minimum elevation) of elevation values within the upstream watershed (in meters),
- HM
_{km}is the area upstream of the calculation point that was burned at high or moderate severity (in km^{2}), and - i15 is the spatially-averaged peak 15-min rainfall intensity for the design storm in the upstream watershed (in mm/h).

Volume estimates were classified in order of magnitude scale ranges 0–1,000 m^{3}; 1,000–10,000 m^{3}; 10,000–100,000 m^{3}; and greater than 100,000 m^{3} for cartographic display.

**Combined Hazard**

Debris-flow hazards from a given basin can be considered as the combination of both probability and volume. For example, in a given setting, the most hazardous basins will show both a high probability of occurrence and a large estimated volume of material. Slightly less hazardous would be basins that show a combination of either relatively low probabilities and larger volume estimates or high probabilities and smaller volume estimates. The lowest relative hazard would be for basins that show both low probabilities and the smallest volumes.

We combined the results of the probability and the volume maps following the methods of Cannon and others (2010). A rank of 1 to 5 (with 5 being the highest) is assigned to each of the probability classes, and a rank of 1 to 4 is assigned to each of the volume classes. The ranks of the probability and volume classes are then added together to produce a map of the combined relative hazard ranking for each basin (with 9 being the highest combined hazard).

**References**

- Schwartz, G.E., and Alexander, R.B., 1995, Soils data for the conterminous United States derived from the NRCS State Soil Geographic (STATSGO) Database: U.S. Geological Survey Open-File Report 95–449, accessed July 2013, article

Below are publications associated with this project.

#### Empirical models for predicting volumes of sediment deposited by debris flows and sediment-laden floods in the transverse ranges of southern California

#### Prediction of spatially explicit rainfall intensity–duration thresholds for post-fire debris-flow generation in the western United States

#### Updated logistic regression equations for the calculation of post-fire debris-flow likelihood in the western United States

#### Probability and volume of potential postwildfire debris flows in the 2012 Waldo Canyon Burn Area near Colorado Springs, Colorado

#### Predicting the probability and volume of postwildfire debris flows in the intermountain western United States

- Overview
The scientific background information for Emergency Assessment of Post-Fire Debris-Flow Hazards.

The preliminary hazard assessment relies upon empirical models to estimate the likelihood and volume of debris flows for selected basins in response to a design storm. Beginning in 2016, the empirical models are based upon historical debris–flow occurrence and magnitude data, rainfall storm conditions, terrain and soils information, and burn–severity data from recently burned areas (Staley et al., 2016, 2017). Hazard assessments for fires prior to May 2016 relied upon an earlier set of empirical methods described here.

The post-fire debris-flow likelihood, volume, and combined hazards are estimated at both the drainage-basin scale and in a spatially distributed manner along the drainage network within each basin. The characteristics of basins affected by the fire were calculated using a geographic information system (GIS). Debris-flow likelihood and volume were estimated for each basin outlet as well as along the upstream drainage networks (pixels where the contributing area is greater than or equal to 0.02 km

^{2}). Independent variable values were calculated for each pixel along the drainage network and summarized at the stream segment scale to obtain estimates of debris-flow likelihood and volume.Required input data include a differenced normalized burn ratio image (dNBR), and field-validated soil burn severity data. If field-validated soil burn severity data are unavailable, it is possible to use a Burned Area Reflectance Classification (BARC4) image in substitution. Please note that the use of BARC4 data increases the degree of uncertainty in the model estimates.

**Likelihood Model**The likelihood of a debris-flow in response to a given peak 15-minute rainfall intensity are based upon a logistic regression approach, which combines the following two equations:

(1) P = e

^{x}/ (1 + e^{x})**Where**- P is the probability of debris-flow occurrence in fractional form, and
- e
^{x}is the exponential function where e represents the mathematical constant 2.718.

For recently burned areas in the western United States, equation 2 is used to calculate x:

(2) x = -3.63 + (0.41 × X1R) + (0.67 × X2R) + (0.7 × X3R)

**Where**- X1R is the proportion of upslope area classified as high or moderate soil burn severity and with gradients ≥ 23°, multiplied by the peak 15-minute rainfall accumulation of the design storm (in millimeters [mm]),
- X2R is the average differenced normalized burn ratio (dNBR) of the upslope area, multiplied by the peak 15-minute rainfall accumulation of the design storm (in millimeters [mm]),
- X3R is the soil KF-Factor (Schwartz and Alexander, 1995) of the upslope area, multiplied by the peak 15-minute rainfall accumulation of the design storm (in millimeters [mm]).

Likelihood values predicted by the equation potentially range from 0 (least likely) to 1 (most likely). The predicted likelihood values are assigned to 1 of 5 equal interval classes for cartographic display and are represented as a percentage likelihood (rather than a ratio).

**Volume Model**Debris-flow volumes both at the basin outlet and along the drainage network are predicted using a multiple linear regression model (Gartner and others, 2014). The multiple linear regression models are used to estimate the volume (V, in m3) of material that could issue from a point along the drainage network in response to a storm of a given rainfall intensity.

Potential debris-flow volume is calculated with equation 4:

(3) ln(V) = 4.22 + (0.13 × sqrt(ElevRange)) + (0.36 × ln(HM

_{km})) + (0.39 × sqrt(i15))**Where**- ElevRange is the range (maximum elevation–minimum elevation) of elevation values within the upstream watershed (in meters),
- HM
_{km}is the area upstream of the calculation point that was burned at high or moderate severity (in km^{2}), and - i15 is the spatially-averaged peak 15-min rainfall intensity for the design storm in the upstream watershed (in mm/h).

Volume estimates were classified in order of magnitude scale ranges 0–1,000 m

^{3}; 1,000–10,000 m^{3}; 10,000–100,000 m^{3}; and greater than 100,000 m^{3}for cartographic display.**Combined Hazard**Debris-flow hazards from a given basin can be considered as the combination of both probability and volume. For example, in a given setting, the most hazardous basins will show both a high probability of occurrence and a large estimated volume of material. Slightly less hazardous would be basins that show a combination of either relatively low probabilities and larger volume estimates or high probabilities and smaller volume estimates. The lowest relative hazard would be for basins that show both low probabilities and the smallest volumes.

We combined the results of the probability and the volume maps following the methods of Cannon and others (2010). A rank of 1 to 5 (with 5 being the highest) is assigned to each of the probability classes, and a rank of 1 to 4 is assigned to each of the volume classes. The ranks of the probability and volume classes are then added together to produce a map of the combined relative hazard ranking for each basin (with 9 being the highest combined hazard).

**References**- Schwartz, G.E., and Alexander, R.B., 1995, Soils data for the conterminous United States derived from the NRCS State Soil Geographic (STATSGO) Database: U.S. Geological Survey Open-File Report 95–449, accessed July 2013, article

- Publications
Below are publications associated with this project.

#### Empirical models for predicting volumes of sediment deposited by debris flows and sediment-laden floods in the transverse ranges of southern California

Debris flows and sediment-laden floods in the Transverse Ranges of southern California pose severe hazards to nearby communities and infrastructure. Frequent wildfires denude hillslopes and increase the likelihood of these hazardous events. Debris-retention basins protect communities and infrastructure from the impacts of debris flows and sediment-laden floods and also provide critical data for vo#### Prediction of spatially explicit rainfall intensity–duration thresholds for post-fire debris-flow generation in the western United States

Early warning of post-fire debris-flow occurrence during intense rainfall has traditionally relied upon a library of regionally specific empirical rainfall intensity–duration thresholds. Development of this library and the calculation of rainfall intensity-duration thresholds often require several years of monitoring local rainfall and hydrologic response to rainstorms, a time-consuming approach w#### Updated logistic regression equations for the calculation of post-fire debris-flow likelihood in the western United States

Wildfire can significantly alter the hydrologic response of a watershed to the extent that even modest rainstorms can generate dangerous flash floods and debris flows. To reduce public exposure to hazard, the U.S. Geological Survey produces post-fire debris-flow hazard assessments for select fires in the western United States. We use publicly available geospatial data describing basin morphology,#### Probability and volume of potential postwildfire debris flows in the 2012 Waldo Canyon Burn Area near Colorado Springs, Colorado

This report presents a preliminary emergency assessment of the debris-flow hazards from drainage basins burned by the 2012 Waldo Canyon fire near Colorado Springs in El Paso County, Colorado. Empirical models derived from statistical evaluation of data collected from recently burned basins throughout the intermountain western United States were used to estimate the probability of debris-flow occur#### Predicting the probability and volume of postwildfire debris flows in the intermountain western United States

Empirical models to estimate the probability of occurrence and volume of postwildfire debris flows can be quickly implemented in a geographic information system (GIS) to generate debris-flow hazard maps either before or immediately following wildfires. Models that can be used to calculate the probability of debris-flow production from individual drainage basins in response to a given storm were de