calc_SOCeq: Equilibrium Soil Organic Carbon (SOCeq)

In ruan-de-wet/snapr: SNAP model in R

Description

To calculate the change in SOC for year t (deltaSOCt), we need the combination of PDSOCt and DDSOCt, but also the maximum rate of microbial respiration for year t (MRESPt). A key input to MRESPt is WETDAYS, which is based on RAIN, the Sand % and SOCt.

Usage

calc_SOCeq(PDSOCt, DDSOCt, SAND, RAIN)

Arguments

PDSOCt	Output of calc_PDSOCt()
DDSOCt	Output of calc_DDSOCt()
SAND	Sand % in top 30 cm soil
RAIN	MAP for year t (mm/year)

Details

SOCt is unknown to us though, so it's not helpful to have to rely on it as an input... We're interested in the SOCt at equilibrium (SOCeq from Eq 12).

By setting deltaSOCt = 0, we can solve for the SOCt term in MRESPt to yield SOCeq.

Not generally called directly but incorporated in wrapper function.

See Also

SNAP which wraps this function.

calc_PDSOCt

calc_DDSOCt

ruan-de-wet/snapr documentation built on Dec. 31, 2020, 3:14 a.m.