Function growth is a forest growth model that calculates changes in leaf area, sapwood area and structural variables for all plant cohorts in a given forest stand during a period specified in the input climatic data.
growth(x, soil, meteo, latitude = NA, elevation = NA, slope = NA, aspect = NA)An object of class growthInput.
A list containing the description of the soil (see soil).
A data frame with daily meteorological data series. Row names of the data frame should correspond to date strings with format "yyyy-mm-dd" (see Date). When x$TranspirationMode = "Granier" the following columns are required:
MeanTemperature: Mean temperature (in degrees Celsius).
Precipitation: Precipitation (in mm).
Radiation: Solar radiation (in MJ/m2/day), required only if snowpack = TRUE.
PET: Potential evapotranspiration (in mm).
WindSpeed: Wind speed (in m/s). If not available, this column can be left with NA values.
When x$TranspirationMode = "Sperry" the following columns are required:
MeanTemperature: Mean temperature(in degrees Celsius).
MinTemperature: Minimum temperature (in degrees Celsius).
MaxTemperature: Maximum temperature (in degrees Celsius).
MinRelativeHumidity: Minimum relative humidity (in percent).
MaxRelativeHumidity: Maximum relative humidity (in percent).
Precipitation: Precipitation (in mm).
Radiation: Solar radiation (in MJ/m2/day).
WindSpeed: Wind speed (in m/s). If not available, this column can be left with NA values.
Latitude (in degrees). Required when x$TranspirationMode = "Sperry".
Elevation above sea level (in m), slope (in degrees) and aspect (in degrees from North). Required when x$TranspirationMode = "Sperry". Elevation is also required for 'Granier' if snowpack dynamics are simulated.
A list of class 'growth' with the following elements:
"latitude": Latitude (in degrees) given as input.
"topography": Vector with elevation, slope and aspect given as input.
"spwbInput": A copy of the object x of class spwbInput given as input.
"soilInput": A copy of the object soil of class soil given as input.
"WaterBalance": A data frame where different variables (in columns) are given for each simulated day (in rows):
"LAIcell": The LAI of the stand (accounting for leaf phenology) (in m2/m2).
"Cm": The water retention capacity of the stand (in mm) (accounting for leaf phenology).
"Lground": The proportion of PAR that reaches the ground (accounting for leaf phenology).
"PET": Potential evapotranspiration (in mm).
"Rainfall": Input precipitation (in mm).
"NetPrec": Net precipitation, after accounting for interception (in mm).
"Infiltration": The amount of water infiltrating into the soil (in mm).
"Runoff": The amount of water exported via surface runoff (in mm).
"DeepDrainage": The amount of water exported via deep drainage (in mm).
"Etot": Evapotranspiration (in mm).
"Esoil": Bare soil evaporation (in mm).
"Eplanttot": Plant transpiration (considering all soil layers) (in mm).
"Eplant.1", ..., "Eplant.k": Plant transpiration from each soil layer (in mm).
"Soil": A data frame where different variables (in columns) are given for each simulated day (in rows):
"W.1", ..., "W.k": Relative soil moisture content (relative to field capacity) in each soil layer.
"ML.1", ..., "ML.k": Soil water volume in each soil layer (in L/m2).
"MLTot": Total soil water volume (in L/m2).
"psi.1", ..., "psi.k": Soil water potential in each soil layer (in MPa).
"PlantTranspiration": A data frame with the amount of daily transpiration (in mm) for each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantPhotosynthesis": A data frame with the amount of daily photosynthesis (in g C <U+00B7> m-2) for each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantRespiration": A data frame with the amount of daily maintenance respiration (in g C <U+00B7> m-2) for each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantCStorageFast": A data frame with the daily amount of fast-dynamics carbon reserves (as proportion of total storage capacity) for an average individual of each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantCStorageSlow": A data frame with the daily amount of slow-dynamics carbon reserves (as proportion of total storage capacity) for an average individual of each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantSAgrowth": A data frame with the daily amount of newly created sapwood area (in cm2) for an average individual of each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantSA": A data frame with the daily amount of sapwood area (in cm2) for an average individual of each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantStress": A data frame with the amount of daily stress suffered by each plant cohort (relative whole-plant conductance). Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantPsi": A data frame with the average daily water potential of each plant (in MPa). Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantLAIdead": A data frame with the dead leaf area index (in m2 <U+00B7> m-2) for each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
"PlantLAIlive": A data frame with the live leaf area index (in m2 <U+00B7> m-2) for each plant cohort. Days are in rows and plant cohorts are in columns. Columns in this data frame correspond to the elements in 'SP'.
Detailed model description is available in the vignettes section. Simulations using the 'Sperry' transpiration mode are computationally much more expensive than those using the simple transpiration mode.
# NOT RUN {
#Load example daily meteorological data
data(examplemeteo)
#Load example plot plant data
data(exampleforest)
#Default species parameterization
data(SpParamsMED)
#Initialize soil with default soil params (2 layers)
examplesoil = soil(defaultSoilParams(2))
#Initialize control parameters
control = defaultControl()
#Initialize input
x = forest2growthInput(exampleforest,examplesoil, SpParamsMED, control)
#Call simulation function
G1<-growth(x, examplesoil, examplemeteo, elevation = 100)
# }
Run the code above in your browser using DataLab