quebec.ldm: Quebec Landscape Dynamic Model
In nuaquilue/QLDM: Quebec Landscape Dynamic Model

Description Usage Arguments Value Examples

Run the landscape dynamic model QLDM that includes the processes of wilfires, clear-cuts, partial-cuts, spruce budworm outbreaks, post-disturbance regeneration, and forest seral succession.

quebec.ldm(
  is.wildfires = FALSE,
  is.sbw = FALSE,
  is.harvesting = FALSE,
  custom.params = NA,
  rcp = NA,
  prec.proj = NA,
  temp.proj = NA,
  timber.supply = "area.based",
  pigni = NA,
  nrun = 1,
  time.step = 5,
  time.horizon = 80,
  save.land = FALSE,
  out.seq = NA,
  out.path = NA,
  ...
)

`is.sbw`	A flag to indicate that spruce budworm outbreaks are a simulated process of the model
`is.harvesting`	A flag to indicate that harvesting by clear and partial cuts is a simulated process of the model
`custom.params`	List with the model paramaters and default and/or user-defined values
`rcp`	Climate projection, either `NA` (default), 'rcp45' or 'rcp85'
`prec.proj`	Data frame with annual precipitation projections for each time step
`temp.proj`	Data frame with mean temperature projections for each time step
`timber.supply`	Approach to do timber supply calculation, either `area.based` or `volume.based`
`pigni`	Optional data frame with probability of fire igntion per grid cell
`nrun`	Number of replicates to run the model
`time.step`	Number of years of each time step
`time.horizon`	Number of years of the model simulation, it has to be a multiple `time.step`
`save.land`	A flag to save as a RDS file the `landscape` data frame at the time step indicated in `out.seq`
`out.seq`	Numeric vector with the time steps the `landscape` is saved
`out.path`	String with the directory path to save the `landscape` data frame at each time step indicated in `out.seq`
`is.widlfire`	A flag to indicate that wildfires are a simualted process of the model

A list with the following items:

SppByAgeClass: A data frame of species abundance per age class, with columns:
- run: Number of replicate.
- year: Year YYYY.
- mgmt.unit: Code of the forest management unit (FMU).
- spp: Code of the species.
- age.class: Code of the age class, C10, C30, C50, C70, C90, and OLD.
- area: Area in km^{2}.
SuitabilityClass: A data frame of suitability of potential species per bioclimatic domain, with columns:
- run: Number of replicate.
- year: Year YYYY.
- bioclim.domain: Code of the bioclimatic domain.
- spp: Code of the species.
- poor: Area in km^{2} of poor environmental suitability.
- med: Area in km^{2} of intermediate environmental suitability.
- good: Area in km^{2} of good environmental suitability.
SppByFireZone: A data frame of species abundance per fire zone, with columns:
- run: Number of replicate.
- year: Year YYYY.
- frz: Code of the fire regime zone.
- spp: Code of the species.
- area: Area in km^{2}.
FuelByFireZone: A data frame of fuel type per fire zone, with columns:
- run: Number of replicate.
- year: Year YYYY.
- frz: Code of the fire regime zone.
- type: Code of the fuel type: low, med or high.
- pct: Relative abundance of the fuel type in the fire regime zone ([0,1]).
Cuts: A data frame of harvestable area and volume per management unit (included if is.harvesting), with columns:
- run: Number of replicate.
- year: Year YYYY.
- mgmt.unit: Code of the forest management unit (FMU).
- a.inc: Area that can be harvested, i.e. non-protected and FMU informed (in km^{2}).
- a.even.age: Area within a.inc of mature stand, i.e.age>age.matu (in km^{2}).
- a.mat.pc: Area to be managed under a partial-cut regime (in km^{2}).
- a.inc.burnt: Area within a.inc burnt in the current time step (in km^{2}).
- a.inc.mat.burnt: Area within a.inc of mature stands burnt in the current time step (in km^{2}).
- a.inc.kill: Area within a.inc killed by SBW in the previous or current time step (in km^{2}).
- a.inc.mat.kill: Area within a.inc of mature stands killed by SBW in the previous or current time step (in km^{2}).
- a.reg.fail.ex: Protected area at defforestation risk (in km^{2}).
- a.reg.fail.in: Non-protected area at defforestation risk (in km^{2}).
- a.salvaged: Salvaged and clear-cut area (in km^{2}).
- a.unaff: Clear-cut area unaffected by disturbances (in km^{2}).
- v.salv: Salvaged and clear-cut volume (in m^{3}).
- v.unaff: Clear-cut volume unaffected by disturbances (in m^{3}).
- a.pcut: Area partial cut (in km^{2}).
- v.pcut: Volume partial cut (in m^{3}).
SppCut: A data frame of area and volum extracted by clear and partial cut per species and management unit (included if is.harvesting), with columns:
- run: Number of replicate.
- year: Year YYYY.
- mgmt.unit: Code of the forest management unit (FMU).
- spp: Code of the species.
- area.ccut: Clear-cut area (in km^{2}).
- vol.ccut: Clear-cut volume (in m^{3}).
- area.pcut: Partial cut area (in km^{2}).
- vol.pcut: Partial cut volume (in m^{3}).
BurntRates: A data frame of target area to be burnt per fire regime zone (included if is.wildfire), with columns:
- run: Number of replicate.
- year: Year YYYY.
- frz: Code of the fire regime zone.
- br: Baseline area to be burnt derived from MFRI (in km^{2}).
- brvar: Baseline area br with inter-period variability added (in km^{2}).
- brfuel: Variable baseline area brvar modified according to zone-level fuel if is.fuel.modifier (in km^{2}).
- brclima: ariable baseline area brvar modified according to zone-level SEP rate if is.clima.modifier (in km^{2}).
- target.area: Target area to be burnt (in km^{2}).
FireRegime: A data frame of number of fires and burnt area per fire regime zone (included if is.wildfire), with columns:
- run: Number of replicate.
- year: Year YYYY.
- frz: Code of the fire regime zone.
- target.area: Target area to be burnt (in km^{2}).
- nfires: Number of fires burnt.
- burnt.area: Area effectively burnt (in km^{2}).
- fire.cycle: Relative fire return interval (in years).
- indx.combust: Fire-zone mean fuel flammability ([0,1]).
- indx.combust.burnt: Mean fuel flammability of actual burnt areas ([0,1]).
Fires: A data frame of wildfires (included if is.wildfire), with columns:
- run: Number of replicate.
- year: Year YYYY.
- frz: Code of the fire regime zone.
- fire.id: Wildfire identificator.
- wind: Main wind direction in degrees.
- target.size: Target area to be burnt (in km^{2}).
- burnt.size: Area effectively burnt (in km^{2}).
- rem: Remanent area to be burnt (in km^{2}).
BurntFuels: A data frame of burnt fuel types per fire regime zone (included if is.wildfire), with columns:
- run: Number of replicate.
- year: Year YYYY.
- frz: Code of the fire regime zone.
- type: Code of the fuel type: low, med or high.
- area: Area burnt of each fuel category (in km^{2}).

## Not run: 
library(QLDM)
data(prec)
data(temp)
# Run one single 80-year replicate with forest management and default RCP 4.5 climate projections 
quebec.ldm(is.harvesting = T, rcp = "rcp45")

## End(Not run)