IFNscenario: IFN scenario simulation engine

In miquelcaceres/IFNdyn: IFN forest dynamics and scenario analysis

Description

Performs simulations of forest dynamics including management prescriptions (on a demand-basis) for a set of forest plots

Usage

IFNscenario(
  treeData,
  plotData,
  prescription,
  strata,
  demand,
  nb = NULL,
  numSteps = 1,
  plotDataDyn = NULL,
  demandOffset = NULL,
  stepOffset = 0,
  as.CO2 = TRUE,
  sequence = FALSE,
  dead.table = FALSE,
  cut.table = FALSE,
  verbose = FALSE,
  ...
)

Arguments

treeData	A data frame with tree records in rows and columns 'ID', 'Species', 'DBH' (in cm), 'H (in m)' and 'N' (ha-1)
plotData	A data frame with plot 'ID' in rows and columns 'SWHC' (mm), 'Rad' (MJ), 'Temp' (degrees C), 'Prec' (mm), 'PET' (mm), 'slope' (degrees), 'Managed' (boolean) and 'Estrato' (stratum code). It can also include a column 'Area' (ha) with the surface corresponding to each forest plot. In this case, 'strata' will not be used.
prescription	A data frame with plots in rows (row names should be equal to 'ID' in tree data) and prescription options in columns (see IFNmanagement)
strata	A data frame with surface areas per stratum (see estratosIFN3). Used to determine the surface corresponding to each plot.
demand	A data frame with annual volume to be extracted by species (see defaultDemand)
nb	An object of class listw (see functions IFNknn and nb2listw).
numSteps	Number of steps to simulate (1 step = 10 years)
plotDataDyn	A data frame with dynamic plot data, with columns 'ID', 'Step', 'Temp' (degrees C), 'Prec' (mm) and 'PET' (mm)
demandOffset	Initial demand offset (for simulations continuing previous ones).
stepOffset	Offset of time steps to read dynamic plot data (for simulations continuing previous ones).
as.CO2	Boolean flag to indicate output as kg of CO2 instead of kg of dry weight (see IFNbiomass).
sequence	Boolean flag to indicate that the sequence of tree lists by step should be returned. If sequence = TRUE, an additional data frame is included in the output, called treeDataSequence. This table includes a column 'Step' to indicate the simulation step.
dead.table	Boolean flag to indicate that a table of trees dead during the simulation is also desired. This table includes a column 'Step' to indicate at which step the trees died.
cut.table	Boolean flag to indicate that a table of trees cut during the simulation is also desired. This table includes a column 'Step' to indicate at which step the trees were cut.
verbose	Boolean flag to indicate extra console output (reporting management interventions in each plot)
...	Additional options to be passed to function IFNdynamics

Value

A list with the following items:

• plots: number of IFN plots included in the simulation.

• totalArea: total area (in ha) that the set of plots represent.

• plotDataSuppl: A data frame with supplementary plot data, with columns:

  • Area: Area (in ha) represented by each plot.

  • DominantSpecies: Code corresponding to the dominant species of each plot.

  • DemandRow: Demand species row index corresponding to each plot.

  • PrescriptionRow: Prescription row index corresponding to each plot.

  • finalPreviousStage: Previous stage of the final cuts (for plots with regular management).

• steps: number of 10-yr steps simulated.

• extractedDemand: a data frame with the volume (m3) extracted for each demand species and time step. The three last columns include: AnnualDemand the mean annual extraction (simulated), Demand the annual demand (input), and Coverage the coverage of the demand (in percent), for each species. The last row ('Total') reports statistics across species.

• managedArea: a data frame with the area (ha) managed for each demand species and time step. The last row ('Total') reports statistics across species.

• volume, density, rootBiomass, stemBiomass, branchBiomass: simulation results by variable. Volume is in m3, density in individuals and biomass in Mg. Each of these elements is in turn a list with the following elements:

  • summary: data frame with a summary of results for each time step (standing, variation, dead and extracted). Here values are given per hectare.

  • standingSpp: volume, biomass or density of living trees, per species and time step.

  • standingPlot: volume, biomass or density of living trees, per plot and time step.

  • standingDomSpp: volume, biomass or density of living trees, per dominant species (i.e. pooling all plots with the same dominant species) and time step.

  • variationSpp: variation (i.e. difference between the current time step and the previous one) in volume, biomass or density of living trees, per species and time step.

  • variationPlot: variation (i.e. difference between the current time step and the previous one) in volume, biomass or density of living trees, per dominant species (i.e. pooling all plots with the same dominant species) and time step.

  • variationDomSpp: variation (i.e. difference between the current time step and the previous one) in volume, biomass or density of living trees, per plot and time step.

  • standgrowthSpp: stand growth in volume, biomass or density of living trees, per species and time step.

  • standgrowthPlot: stand growth in volume, biomass or density of living trees, per plot and time step.

  • standgrowthDomSpp: stand growth in volume, biomass or density of living trees, per dominant species (i.e. pooling all plots with the same dominant species) and time step.

  • deadSpp: volume, biomass or density of dead trees (i.e. dead during the corresponding time step), per species and time step.

  • deadPlot: volume, biomass or density of dead trees (i.e. dead during the corresponding time step), per plot and time step.

  • deadDomSpp: volume, biomass or density of dead trees (i.e. dead during the corresponding time step), per dominant species (i.e. pooling all plots with the same dominant species) and time step.

  • extractedSpp: volume, biomass or density of trees extracted (during the corresponding time step), per species and time step.

  • extractedPlot: volume, biomass or density of trees extracted (during the corresponding time step), per plot and time step.

  • extractedDomSpp: volume, biomass or density of trees extracted (during the corresponding time step), per dominant species (i.e. pooling all plots with the same dominant species) and time step.

• treeDataNew: tree data frame at the end of the simulated period.

• treeDataSequence: tree data frame for the whole sequence of steps (if sequence = TRUE).

• treeDataDead: tree data frame with dead trees for every step (if dead.table = TRUE).

• treeDataCut: tree data frame with cut trees for every step (if cut.table = TRUE).

Examples

## Not run: 
library(IFNread)
data(defaultDemand)
data(defaultPrescription)
data(estratosIFN3)
data(plotDataIFN)

accessfilesIFN3 = list.files("D:/Recerca/Datasets/IFN/Sources/IFN3/MSAccess/",
pattern = "*[.]accdb", full.names = T)

prov = 27 # Lugo

treeData = readPiesMayoresIFN3(accessfilesIFN3[prov])

IFNscenario(treeData, plotDataIFN, defaultPrescription,
            estratosIFN3[[prov]], defaultDemand[[prov]], numSteps = 3)

## End(Not run)

miquelcaceres/IFNdyn documentation built on Feb. 1, 2021, 10:55 a.m.