View source: R/massLiveRoots.R
massLiveRoots | R Documentation |
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)
massLiveRoots(
layerBottom,
layerTop,
totalRootMassPerArea,
rootDepthMax,
soilLength = 1,
soilWidth = 1,
shape,
expDecayRatePerMaxDepth = log(0.05),
...
)
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 |
expDecayRatePerMaxDepth |
a numeric, value for exponential decay rate per maximum root depth |
... |
an array of the associated live root mass for each layer specified.
massLiveRoots(layerBottom = 1:10, layerTop = 0:9,
totalRootMassPerArea = 0.3, rootDepthMax=5,
soilLength=1, soilWidth=1, shape='linear')
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.