fbpRaster: Raster-based Fire Behavior Prediction System Calculations

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/fbpRaster.r


fbpRaster calculates the outputs from the Canadian Forest Fire Behavior Prediction (FBP) System (Forestry Canada Fire Danger Group 1992) based on raster format fire weather and fuel moisture conditions (from the Canadian Forest Fire Weather Index (FWI) System (Van Wagner 1987)), fuel type, date, and slope. Fire weather, for the purpose of FBP System calculation, comprises observations of 10 m wind speed and direction at the time of the fire, and two associated outputs from the Fire Weather Index System, the Fine Fuel Moisture Content (FFMC) and Buildup Index (BUI). Raster-based FWI System components can be calculated with the sister function fwiRaster.





The input data, a RasterStack containing fuel types, fire weather component, and slope layers (see below). Each vector of inputs defines a single FBP System prediction for a single fuel type and set of weather conditions. The RasterStack can be used to evaluate the FBP System for a single fuel type and instant in time, or multiple records for a single point (e.g., one weather station, either hourly or daily for instance) or multiple points (multiple weather stations or a gridded surface). All input variables have to be named as listed below, but they are case insensitive, and do not have to be in any particular order. Fuel type is of type character; other arguments are numeric. Missing values in numeric variables could either be assigned as NA or leave as blank.

Required Inputs:

Input Description/Full name Defaults
FuelType FBP System Fuel Type including "C-1", "C-2", "C-3", "C-4", "C-5","C2" "C-6", "C-7", "D-1", "M-1", "M-2", "M-3", "M-4", "NF", S-1", "S-2", "S-3", "O-1a", "O-1b", and "WA", where "WA" and "NF" stand for "water" and"non-fuel", respectively.Because raster format data cannot hold characters, we have to code these fuel types in numeric codes. In sequence, the codes are c(1:19) (see below). FuelType could also be converted as factor and assigned to the raster layer, the function will still work.
Fuel Type code
C-1 1
C-2 2
C-3 3
C-4 4
C-5 5
C-6 6
C-7 7
D-1 8
M-1 9
M-2 10
M-3 11
M-4 12
NF 13
O-1a 14
O-1b 15
S-1 16
S-2 17
S-3 18
WA 19
LAT Latitude [decimal degrees] 55
LONG Longitude [decimal degrees] -120
FFMC Fine fuel moisture code [FWI System component] 90
BUI Buildup index [FWI System component] 60
WS Wind speed [km/h] 10
GS Ground Slope [percent] 0
Dj Julian day 180
Aspect Aspect of the slope [decimal degrees] 0

Optional Inputs (1): Variables associated with certain fuel types. These could be skipped
if relevant fuel types do not appear in the input data.

Input Full names of inputs Defaults
PC Percent Conifer for M1/M2 [percent] 50
PDF Percent Dead Fir for M3/M4 [percent] 35
cc Percent Cured for O1a/O1b [percent] 80
GFL Grass Fuel Load [kg/m^2] 0.35

Optional Inputs (2): Variables that could be ignored without causing major impacts to the
primary outputs

Input Full names of inputs Defaults
CBH Crown to Base Height [m] 3
WD Wind direction [decimal degrees] 0
Accel Acceleration: 1 = point, 0 = line 0
ELV* Elevation [meters above sea level] NA
BUIEff Buildup Index effect: 1=yes, 0=no 1
D0 Julian day of minimum Foliar Moisture Content 0
hr Hours since ignition 1
ISI Initial spread index 0
CFL Crown Fuel Load [kg/m^2] 1.0
FMC Foliar Moisture Content if known [percent] 0
SH C-6 Fuel Type Stand Height [m] 0
SD C-6 Fuel Type Stand Density [stems/ha] 0
theta Elliptical direction of calculation [degrees] 0

FBP output offers 3 options (see details in Values section):

Outputs Number of outputs
Primary (default) 8
Secondary 34
All 42

Selected outputs


m Optimal number of pixels at each iteration of computation when ncell(input) >= 1000. Default m = NULL, where the function will assign m = 1000 when ncell(input) is between 1000 and 500,000, and m=3000 otherwise. By including this option, the function is able to process large dataset more efficiently. The optimal value may vary with different computers.


Number of CPU cores (integer) used in the computation, default is 1. By signing cores > 1, the function will apply parallel computation technique provided by the foreach package, which significantly reduces the computation time for large input data (over a million grid points). For small dataset, cores=1 is actually faster.

* Elevation is only used in the calculation of Foliar Moisture Content (FMC). However, FMC can also be calculated without elevation input. The default is to not use elevation in the calculation of FMC.


The Canadian Forest Fire Behavior Prediction (FBP) System (Forestry Canada Fire Danger Group 1992) is a subsystem of the Canadian Forest Fire Danger Rating System, which also includes the Canadian Forest Fire Weather Index (FWI) System. The FBP System provides quantitative estimates of head fire spread rate, fuel consumption, fire intensity, and a basic fire description (e.g., surface, crown) for 16 different important forest and rangeland types across Canada. Using a simple conceptual model of the growth of a point ignition as an ellipse through uniform fuels and under uniform weather conditions, the system gives, as a set of secondary outputs, estimates of flank and back fire behavior and consequently fire area perimeter length and growth rate.

The FBP System evolved since the mid-1970s from a series of regionally developed burning indexes to an interim edition of the nationally develop FBP system issued in 1984. Fire behavior models for spread rate and fuel consumption were derived from a database of over 400 experimental, wild and prescribed fire observations. The FBP System, while providing quantitative predictions of expected fire behavior is intended to supplement the experience and judgment of operational fire managers (Hirsch 1996).

The FBP System was updated with some minor corrections and revisions in 2009 (Wotton et al. 2009) with several additional equations that were initially not included in the system. This fbp function included these updates and corrections to the original equations and provides a complete suite of fire behavior prediction variables. Default values of optional input variables provide a reasonable mid-range setting. Latitude, longitude, elevation, and the date are used to calculate foliar moisture content, using a set of models defined in the FBP System; note that this latitude/longitude-based function is only valid for Canada. If the Foliar Moisture Content (FMC) is specified directly as an input, the fbp function will use this value directly rather than calculate it. This is also true of other input variables.

Note that Wind Direction (WD) is the compass direction from which wind is coming. Wind azimuth (not an input) is the direction the wind is blowing to and is 180 degrees from wind direction; in the absence of slope, the wind azimuth is coincident with the direction the head fire will travel (the spread direction azimuth, RAZ). Slope aspect is the main compass direction the slope is facing. Slope azimuth (not an input) is the direction a head fire will spread up slope (in the absence of wind effects) and is 180 degrees from slope aspect (Aspect). Wind direction and slope aspect are the commonly used directional identifiers when specifying wind and slope orientation respectively. The input theta specifies an angle (given as a compass bearing) at which a user is interested in fire behavior predictions; it is typically some angle off of the final spread rate direction since if for instance theta=RAZ (the final spread azimuth of the fire) then the rate of spread at angle theta (TROS) will be equivalent to ROS.


fbpRaster returns a RasterStack with primary, secondary, or all output variables, a combination of the primary and secondary outputs. Primary FBP output includes the following 8 raster layers:


Crown Fraction Burned by the head fire


Crown Fuel Consumption [kg/m^2]


Fire description (1=Surface, 2=Intermittent, 3=Crown)


Head Fire Intensity [kW/m]


Spread direction azimuth [degrees]


Equilibrium Head Fire Rate of Spread [m/min]


Surface Fuel Consumption [kg/m^2]


Total Fuel Consumption [kg/m^2]

Secondary FBP System outputs include the following 34 raster layers. In order to calculate the reliable secondary outputs, depending on the outputs, optional inputs may have to be provided.


BUI effect on spread rate


Slope Factor (multiplier for ROS increase upslope)


Initial Spread Index


Fine fuel moisture code [FWI System component]


Foliar Moisture Content [%]


Julian Date of minimum FMC


Critical spread rate for crowning [m/min]


Critical Surface Intensity for crowning [kW/m]


Equilibrium Flank Fire Rate of Spread [m/min]


Equilibrium Back Fire Rate of Spread [m/min]


Head Fire Rate of Spread at time hr [m/min]


Flank Fire Rate of Spread at time hr [m/min]


Back Fire Rate of Spread at time hr [m/min]


Flank Fire Crown Fraction Burned


Back Fire Crown Fraction Burned


Equilibrium Spread Flank Fire Intensity [kW/m]


Equilibrium Spread Back Fire Intensity [kW/m]


Flank Fire Total Fuel Consumption [kg/m^2]


Back Fire Total Fuel Consumption [kg/m^2]


Head Fire Spread Distance after time hr [m]


Back Fire Spread Distance after time hr [m]


Flank Fire Spread Distance after time hr [m]


Time to Crown Fire Initiation [hrs since ignition]


Time to Flank Fire Crown initiation [hrs since ignition]


Time to Back Fire Crown initiation [hrs since ignition]


Length to Breadth ratio


Length to Breadth ratio after elapsed time hr


Net vectored wind speed [km/hr]


Equilibrium Rate of Spread at bearing theta [m/min]


Rate of Spread at bearing theta at time t [m/min]


Crown Fraction Burned at bearing theta


Fire Intensity at bearing theta [kW/m]


Total Fuel Consumption at bearing theta [kg/m^2]


Time to Crown Fire initiation at bearing theta [hrs since ignition]

*These outputs represent fire behaviour at a point on the perimeter of an elliptical fire defined by a user input angle theta. theta represents the bearing of a line running between the fire ignition point and a point on the perimeter of the fire. It is important to note that in this formulation the theta is a bearing and does not represent the angle from the semi-major axis (spread direction) of the ellipse. This formulation is similar but not identical to methods presented in Wotton et al (2009) and Tymstra et al (2009).


Xianli Wang, Alan Cantin, Marc-André Parisien, Mike Wotton, Kerry Anderson, and Mike Flannigan


1. Hirsch K.G. 1996. Canadian Forest Fire Behavior Prediction (FBP) System: user's guide. Nat. Resour. Can., Can. For. Serv., Northwest Reg., North. For. Cent., Edmonton, Alberta. Spec. Rep. 7. 122p.

2. Forestry Canada Fire Danger Group. 1992. Development and structure of the Canadian Forest Fire Behavior Prediction System. Forestry Canada, Ottawa, Ontario Information Report ST-X-3. 63 p. http://cfs.nrcan.gc.ca/pubwarehouse/pdfs/10068.pdf

3. Wotton, B.M., Alexander, M.E., Taylor, S.W. 2009. Updates and revisions to the 1992 Canadian forest fire behavior prediction system. Nat. Resour. Can., Can. For. Serv., Great Lakes For. Cent., Sault Ste. Marie, Ontario, Canada. Information Report GLC-X-10, 45p. http://publications.gc.ca/collections/collection_2010/nrcan/Fo123-2-10-2009-eng.pdf

4. Tymstra, C., Bryce, R.W., Wotton, B.M., Armitage, O.B. 2009. Development and structure of Prometheus: the Canadian wildland fire growth simulation Model. Nat. Resour. Can., Can. For. Serv., North. For. Cent., Edmonton, AB. Inf. Rep. NOR-X-417.

See Also

fbp, fwiRaster, hffmcRaster


# The dataset is the standard test data for FBP system
# provided by Wotton et al (2009), and randomly assigned
# to a stack of raster layers
test_fbpRaster <- stack(system.file("extdata", "test_fbpRaster.tif", package="cffdrs"))
# Stack doesn't hold the raster layer names, we have to assign
# them:
names(input)<-c("FuelType","LAT","LONG","ELV","FFMC","BUI", "WS","WD","GS","Dj","D0","hr","PC",
# Primary outputs:
system.time(foo<-fbpRaster(input = input))
# Using the "select" option:
system.time(foo<-fbpRaster(input = input,select=c("HFI","TFC", "ROS")))
# Secondary outputs:
system.time(foo<-fbpRaster(input = input,output="S"))
# All outputs:
#system.time(foo<-fbpRaster(input = input,output="A"))

### Additional, longer running examples  ###
# Keep only the required input layers, the other layers would be
# assigned with default values:
# keep only the required inputs:
dat0<-input[[c("FuelType","LAT","LONG","FFMC","BUI","WS","GS", "Dj","Aspect")]]
system.time(foo<-fbpRaster(input = dat0,output="A"))

Example output

Loading required package: rgdal
Loading required package: sp
rgdal: version: 1.5-18, (SVN revision 1082)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 3.0.4, released 2020/01/28
Path to GDAL shared files: /usr/share/gdal
GDAL binary built with GEOS: TRUE 
Loaded PROJ runtime: Rel. 6.3.1, February 10th, 2020, [PJ_VERSION: 631]
Path to PROJ shared files: /usr/share/proj
Linking to sp version:1.4-4
To mute warnings of possible GDAL/OSR exportToProj4() degradation,
use options("rgdal_show_exportToProj4_warnings"="none") before loading rgdal.
Loading required package: raster
Loading required package: foreach
FD = 1,2,3 representing Surface (S),Intermittent (I), and Crown (C) fire
   user  system elapsed 
  0.351   0.004   0.374 
Warning messages:
1: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced
2: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced
   user  system elapsed 
  0.223   0.000   0.229 
Warning messages:
1: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced
2: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced
   user  system elapsed 
  1.101   0.000   1.115 
Warning messages:
1: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced
2: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced
FD = 1,2,3 representing Surface (S),Intermittent (I), and Crown (C) fire
   user  system elapsed 
  1.625   0.004   1.649 
Warning messages:
1: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced
2: In log(1 - (RSF/a[FUELTYPE])^(1/c0[FUELTYPE])) : NaNs produced

cffdrs documentation built on July 1, 2020, 6:04 p.m.