twoStation_dDOdt: Time derivative of DO concentration

In mtalluto/NSmetabolism: N-Station Stream Metabolism Modelling

Description

Time derivative of DO concentration

Usage

twoStation_dDOdt(t, y, parms, data, return_list = FALSE)

Arguments

t	Time step for the ode integration
y	Current state (dissolved oxygen concentration) at both sites
parms	a vector of parameters; see details
data	a list of data for the ode; see details
return_list	logical, should a list (suitable for deSolve) be returned?

Details

params must be a vector. The indices of the vector must map to model parameters as follows; (where parameters with length two are fit per site):

• 1:2 – P1 $(W min g^-1 O_2)$ – inverse of the slope of a photosynthesis–irradiance curve at low light intensity

• 3:4 – P2 $(m^2 min g^-1 O_2)$ – inverse maximum photosynthesis rate; can be zero to assume GPP is linear with light intensity instead of saturating

• 5:6 – k600 – coefficient of gas exchange for a gas with a Schmidt number of 600

• 7:8 – ER24_20 – daily ecosystem respiration rate, standardized at 20 degrees C

data is a list of length 2 (one per site); each item is a list of (constant) data items, including the following:

• PAR a function that takes time as a parameter and returns light intensity $(W/m^2)$

• temp a function that takes time as a parameter and returns water temperature (degrees C)

• P pressure, in atmospheres

• z Depth, in meters

Value

Time derivative of dissolved oxygen

References

Fuß, T. et al. (2017). Land use controls stream ecosystem metabolism by shifting dissolved organic matter and nutrient regimes. Freshw Biol 62:582–599. \ Van de Bogert MC et al. 2007. Assessing pelagic and benthic metabolism using free water measurements. Limnol. Oceanogr.: Methods 5:145–155. \ Uehlinger, U., et al. 2000. Variability of photosynthesis‐irradiance curves and ecosystem respiration in a small river. Freshw Biol 44:493–507.

mtalluto/NSmetabolism documentation built on May 3, 2021, 7:51 p.m.