massLiveRoots: Mass of living roots

In tilbud/MEMs: R-Cohort Marsh Equilibrium Model

Mass of living roots

Description

This function returns the mass of the living roots between two layers. Given a mass per area of roots and a maximum rooting depth, this calculates the total mass between two layers with a specified distribution. Currently only a linear distribution is implemented.
The linear algorithm is as follows:

• #mass_per_depth = slope * depth + intercept

• slope <- -2 * totalRootMass / (rootDepthMax^2)

• intercept <- 2 * totalRootMass / rootDepthMax

• #mass = integral(mass_per_depth, depth)

• rootMass <- intercept * (layerBottom-layerTop) + slope/2 * (layerBottom ^2-layerTop^2)

Usage

massLiveRoots(
  layerBottom,
  layerTop,
  totalRootMassPerArea,
  rootDepthMax,
  soilLength = 1,
  soilWidth = 1,
  shape,
  expDecayRatePerMaxDepth = log(0.05),
  ...
)

Arguments

layerBottom	an array of depths from soil top to the bottom of the layer, generally in cm
layerTop	an array of depths from soil top to the top of the layer, generally in cm
totalRootMassPerArea	an integer that is the total mass per area of the roots, generally in g cm-3
rootDepthMax	an integer that is the maximum root depth, generally in cm
soilLength	unit length of soil area, generally 1
soilWidth	unit width of soil area, generally 1.
shape	flag of root shape, only linear is implimented.
expDecayRatePerMaxDepth	a numeric, value for exponential decay rate per maximum root depth
...

Value

an array of the associated live root mass for each layer specified.

Examples

massLiveRoots(layerBottom = 1:10, layerTop = 0:9, 
              totalRootMassPerArea = 0.3, rootDepthMax=5,
              soilLength=1, soilWidth=1, shape='linear')

tilbud/MEMs documentation built on April 3, 2024, 2:12 p.m.