Description Usage Arguments Value Author(s) References See Also Examples
Include corrections to the orginal model by Frandsen (1973), Albini (1976), and Andrews et al. (2013).
1 | ros (modeltype, w, s, delta, mx.dead, h, m, u, slope)
|
modeltype |
S(tatic), D(ynamic) |
w |
a vector or data frame of fuel load [t/ha] for fuel classes 1-hour, 10-hour, 100-hour, live herbs and live woody, respectively (5 values or columns; 0 if fuel class is absent). |
s |
a vector or data frame of surface-to-volume ratio [m2/m3] for fuel classes 1-hour, 10-hour, 100-hour, live herbs and live woody, respectively (5 values or columns; 0 if fuel class is absent). |
delta |
a value or vector of fuel bed depth [cm] |
mx.dead |
a value or vector of dead fuel moisture of extinction [percent] |
h |
a vector or data frame of heat content [kJ/kg] for fuel classes 1-hour, 10-hour, 100-hour, live herbs and live woody, respectively (5 values or columns; 0 if fuel class is absent). |
m |
a vector or data frame of percent moisture on a dry weight basis (percent) for fuel classes 1-hour, 10-hour, 100-hour, live herbs and live woody, respectively (5 values or columns; 0 if fuel class is absent). |
u |
a value or vector of midflame windspeed [km/h] |
slope |
a value or vector of site slope [percent] |
A list of values or vectors for the following variables: [1] Characteristic dead fuel moisture [percent], [2] Characteristic live fuel moisture [percent], [3] Live fuel moisture of extinction [percent], [4] Characteristic SAV [m2/m3] [5] Bulk density [kg/m3], [6] Packing ratio [dimensionless], [7] Relative packing ratio [dimensionless], [8] Dead fuel Reaction intensity [kW/m2], [9] Live fuel Reaction intensity [kW/m2], [10] Reaction intensity [kW/m2], [11] Wind factor [0-100], [12] Slope factor [0-1], [13] Heat source [kW/m2], [14] Heat sink [kJ/m3], [15] ROS [m/min].
Giorgio Vacchiano, Davide Ascoli (DISAFA, University of Torino, Italy)
Albini, F. A. (1976). Computer-based models of wildland fire behavior: A users' manual. Ogden, UT: US Department of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station.
Andrews, P. L., Cruz, M. G., and Rothermel, R. C. (2013). Examination of the wind speed limit function in the Rothermel surface fire spread model. International Journal of Wildland Fire 22 (7): 959-969. http://dx.doi.org/10.1071/WF12122.
Frandsen, W. H. (1973). Using the effective heating number as a weighting factor in Rothermel's fire spread model. Ogden, UT: US Department of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station.
Rothermel, R. C. (1972). A mathematical model for fire spread predictions in wildland fires. Research Paper INT-115. Ogden, UT: US Department of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 | # Example 1: Simulation using vectors of input values
modeltype <- "D"
w <-c (2, 1, 0.5, 3, 8)
s <- c (5600, 358, 98, 6200, 8000)
delta <- 50
mx.dead <- 30
h <- c (18622, 18622, 18622, 19500, 20000)
m <- c (7, 8, 9, 40, 60)
u <- 5
slope <- 10
ros(modeltype, w, s, delta, mx.dead, h, m, u, slope)
# Example 2: variable wind input
# Only rate of spread is reported here (i.e., element [15] of ros ( ) output)
modeltype <- "D"
w <-c (2, 1, 0.5, 3, 8)
s <- c (5600, 358, 98, 6200, 8000)
delta <- 50
mx.dead <- 30
h <- c (18622, 18622, 18622, 19500, 20000)
m <- c (7, 8, 9, 40, 60)
slope <- 10
df <- data.frame ("windspeed" = seq (3, 15, 1), ROS=NA)
for (i in 1:nrow (df)) {
df [i,2] <-
ros (modeltype, w, s, delta, mx.dead, h, m, u=df[i,1], slope) [15]
}
df
# Example 3: variable wind and slope input
# A two-entry table of rates of spread is created
modeltype <- "D"
w <-c (2, 1, 0.5, 3, 8)
s <- c (5600, 358, 98, 6200, 8000)
delta <- 50
mx.dead <- 30
h <- c (18622, 18622, 18622, 19500, 20000)
m <- c (7, 8, 9, 40, 60)
u <- seq (3, 15, 1)
slope <- seq (0, 45, 15)
df <- matrix (rep (NA, length (u) * length (slope)),
length (u),
length (slope)
)
df <- data.frame (u, df)
colnames (df) <- c ("windspeed",
paste ("slope_", as.character (slope))
)
for (i in 1:length (u)) {
for (j in 1:length (slope)) {
df [i, j+1] <- ros (
modeltype, w, s, delta, mx.dead, h, m, u[i], slope[j])[15]
}
}
df
# Example 4: prediction and validation of rate of spread
# using existing data from a field experiment
library (Rothermel)
# Observed variables
data (firexp)
m <- firexp [, 18:22]
u <- firexp [, "u"]
slope <- firexp [, "slope"]
obs <- firexp [, "ros"]
# Predict ROS using Standard Fuel Models GR5, GS3 and SH7
data (SFM_metric)
a = list ( )
models = which (rownames (SFM_metric) == "GR5" |
rownames (SFM_metric) == "GS3" |
rownames (SFM_metric) == "SH7")
for (i in 1 : length (models) ) {
modeltype <- SFM_metric [models [i], 1]
w <- SFM_metric [models [i], 2:6]
s <- SFM_metric [models [i], 7:11]
delta <- SFM_metric [models [i], "Fuel_Bed_Depth"]
mx.dead <- SFM_metric [models [i], "Mx_dead"]
h <- SFM_metric [models [i], 14:18]
a [i] <- ros (modeltype, w, s, delta, mx.dead, h,
m, u, slope)[15]}
# Plot
plot (obs, a [[1]], xlab = "Observed rate of spread (m/min)",
ylab = "Predicted rate of spread (m/min)", col = "red",
pch =19, xlim = c (0, 30), cex.lab = 1.1)
points (obs, a [[2]], pch = 19, col = "green2")
points (obs, a [[3]], pch = 19, col = "blue2")
abline (coef = c(0, 1))
abline (coef = c(0, 0.7),lty = 2); text (13.6, 19.2, "-30 percent")
abline (coef = c(0, 1.3),lty = 2); text (28.7, 19.2, "+30 percent")
legend (0, 19.2, c("GR5", "GS3", "SH7"), pch = 19,
col = c("red", "green2", "blue2"), title = "Fuel model")
# Inset Residual plot (not run)
par (fig = c (.57, .98, .07, .55), new = TRUE)
plot (obs, a[[1]] - obs, xlab= "", ylab= "", col = "red",
main= "Residuals", font.main = 1, pch=19, cex=.7)
points (obs, a [[2]] - obs, pch = 19, cex =.7, col = "green2")
points (obs, a [[3]] - obs, pch = 19, cex =.7, col = "blue2")
abline (h = 0)
par (fig = c (0, 1, 0, 1))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.