StandAttributes: Calculate stand attributes

In ptompalski/MISTR: Implementation of Ontario's forest growth and yield models.

Calculate stand attributes

Description

A wrapper function to calculate all available stand attributes using other MISTR functions.

Usage

StandAttributes(
  SI,
  age,
  species = NA,
  SFU = NA,
  stocking = 1,
  Origin = "N",
  PlantedSpp = NA,
  bhage = NA
)

Arguments

SI	site index
age	stand total age
species	Id of dominant species. May be NA (default) if SFU is provided.
SFU	Id of Standard Forest Unit (SFU). May be NA (default) if species is provided. It is recommended to enter both SFU and species if attributes are calculated for SFU.
stocking	Stocking. Default=1.
Origin	Origin of the stand. Use "N" for natural, "P" for plantation. Default="N".
PlantedSpp	Opional information on planted species required for some plantations.
bhage	Optional. Breast-height age. If not provided calculated automatically based on species and age.

Details

This function calculates a set of stand attributes including top height (H), basal area (BA), quadratic mean diameter (QMD), gross total volume (GTV), gross merchantable volume (GMV), net merchantable volume (NMV), volume-to-basal area ratio (VBAR), and cull.

Attributes may be calculated for:

• dominant species (argument species is provided, argument SFU is NA)

• standard forest unit (both arguments SFU and species are provided (recommended)). If species is missing, the function will try to identify the species automatically, which is possible for some SFUs. If the species is successfully identified the function will inform the user and return all stand attributes. If the species cannot be identified automatically, only a limited set of attributes will be returned (BA, QMD, N).

Value

A data frame containing the calculated stand attributes.

Examples

#-------------
#Calculate attributes for a 30-year-old White Ash stand, with site index of 20:
StandAttributes(SI = 20, age = 30, species = "AW")

#Calculate attributes for a Red pine stand with site index of 18, for a sequence of ages between 10 and 50:
StandAttributes(SI = 18, age = 10:50, species = "PR")

#Calculate attributes for a stand classified as PJ1 forest unit, dominated by Jack Pine (PJ), with site index of 18, and age of 85
StandAttributes(SI = 18, age = 85, species = "PJ", SFU = "PJ1")

#In some cases the species can be omitted.
StandAttributes(SI = 20, age = 50, SFU = "PJ1")

#In some cases species is required to calculate all attributes. Without species only QMD, BA, and N are calculated.
StandAttributes(SI = 20, age = 50, SFU = "SF1")

#-------------
#Calculating stand attributes for planted stands:
StandAttributes(SI = 20, age = 50, SFU = "PJ1", species = "PJ", Origin = "P")

#For some planted stands PlantedSpp argument is requried
StandAttributes(SI = 20, age = 50, SFU = "SP1", species = "SW", Origin = "P", PlantedSpp = "SW")

ptompalski/MISTR documentation built on March 18, 2024, 10:13 p.m.