scenarioEmissions: Scenario biomass and emissions
In miquelcaceres/IFNdyn: IFN forest dynamics and scenario analysis

Function scenarioLiveCarbonStock estimates the amount of living biomass per time step. Functions scenarioDeadWood and scenarioProducts estimate the dead and product biomass generated in each time step. Finally, function scenarioEmissions estimates the carbon source and sink effects for each time step of the simulated scenario.

scenarioProducts(scenario, productDestination, variable = "biomass", ...)

scenarioDeadWood(scenario, as.CO2 = FALSE, area = NA)

scenarioLiveCarbonStock(scenario, as.CO2 = FALSE, area = NA)

scenarioEmissions(
  scenario,
  productDestination,
  initialDead = c(roots_dead = 5, coarse_dead = 3, medium_dead = 3, fine_dead = 1),
  initialProducts = c(trituration = 0.5, poles = 0.1, boards = 0.5, construction = 0.5,
    firewood = 0.1, veneer = 0.1),
  deadHalfLives = c(roots_dead = 20, coarse_dead = 10, medium_dead = 5, fine_dead = 2),
  productHalfLives = c(trituration = 2, poles = 25, boards = 15, construction = 35,
    firewood = 2, veneer = 25),
  as.CO2 = FALSE,
  area = NA
)

`scenario`	The result of function `IFNscenario`.
`productDestination`	Data frame of product destination by species and diameter class
`variable`	Either 'biomass' (default) or 'volume'
`...`	Additional parameters to `IFNproducts`.
`as.CO2`	Flag to indicate output as kg of CO2 / ha instead of kg of dry weight / ha. Percentage of carbon per dry weight biomass by species are taken from Montero et al. (2005) (in turn, from Ibáñez et al. 2002).
`area`	Either 'Atlantic' or 'Mediterranean' to specify allometric equations specific to the area (for Pinus pinaster)
`initialDead`	Initial values of dead matter components (in Mg/ha)
`initialProducts`	Initial values of product stocks (in Mg/ha)
`deadHalfLives`	Half-live (in yrs) of dead matter components.
`productHalfLives`	Half-live (in yrs) of products.

Functions scenarioLiveCarbonStock and scenarioDeadWood return a list of two data frames with biomass (dry weight or CO2, depending on as.CO2). The first data data frame contains biomass values in Mg for the target area, whereas the second contains biomass values in Mg/ha.

Function scenarioProducts can estimate both biomass or volumes (by calling function IFNproducts). Hence, and compared to the previous functions, the two data frames can correspond to volumes (m3 and m3/ha, respectively) instead of biomass values.

Function scenarioEmissions returns a list with the following elements:

initialDead: Initial values of dead matter components (in Mg/ha).
initialProducts: Initial values of products (in Mg/ha).
liveWoodChange: Changes in live biomass components (in Mg/ha) per time step.
deadWoodGeneration: Generation of dead biomass components (in Mg/ha) per time step.
productWoodGeneration: Generation of products (in Mg/ha) per time step.
carbonBalance: List of two frames with carbon balance details (including net emissions) in Mg/yr or Mg/ha/yr, respectively.
finalDead: Final values of dead matter components (in Mg/ha).
finalProducts: Final values of products (in Mg/ha).