EcosystemProduction.20 {StreamMetabolism} | R Documentation |
Calculates Ecosytem Production with respiration corrected for temperature.
EcosystemProduction.20(velocity, meandepth, temperature, DO, day, startday, endday, sunrise.time, sunset.time, num.readings)
velocity |
velocity in meters/s |
meandepth |
mean stream channel depth in meters |
temperature |
temperature in degrees Celcius |
DO |
Dissolved Oxygen in mg/L |
day |
date of the day of interest must be in quotes |
startday |
time of the start of the "day" usually 00:00:00 must be in quotes |
endday |
time of the end of the "day" usually 23:45:00 must be in quotes |
sunrise.time |
time of sunrise in the form "04:22:00" must be in quotes |
sunset.time |
time of sunset in the form "19:23:00" must be in quotes |
num.readings |
number of readings if data is collected in 15min intervals then there are 96 readings per day |
**The only difference between this function and EcosystemProduction is the correction for temperature for respiration. I think EcosystemProduction.20 (The function documented in this page) should be favoured over just EcosystemProduction because of it should better estimate Metabolism.**
O'Conner Dobbins eauation is used to calculate Rearation coefficient. This function has only been tested with zoo objects read in with the read.production function. Should work for others, but untested. The values that are returned are only as good as the assumptions that went into collecting the data and should be viewed with scrutiny. Also the Reaeration coefficient is the most important number used in the calculation of Production and therefore needs to be as correct as possible. In the future (as the need arises- or enough requests come in) I plan on incorperating different ways of calculating this, but for now you are stuck with my assumptions for the task at hand.
CR24 |
Community Repiration 24hours (sum of pre-dawn and post-dusk change in Oxygen curve corrected for Reaeration) |
NDM |
Net Daily Metabolism (sum of total day change in Oxygen Curve corrected for Reaeration) |
abs.CR24 |
absolute value of Community Repiration 24hours (because CR24 is always negative) |
GPP |
Gross Primary Productivity=NDM+abs.CR24 |
P.R |
P/R ratio |
Units |
areal units (gOxygen/squared meter/day) of everything except P.R which is unitless |
Stephen A Sefick Jr.
Odum, H. T. (1956). "Primary production in flowing waters." Limnology and Oceanography, 1: 102-117.
Thyssen, N., Erlandsen, M., Jeppesen, E., Holm T. F.,1983. Modelling the reaeration capacity of low-land streams
M.R. Grace and S.J. Imberger. 2006. "Stream Metabolism: Performing & Interpreting Measurements". Water Studies Centre Monash University, Murray Darling Basin Commission and New South Wales Department of Environment and Climate Change. 204 pp. Accessed at http://www.sci.monash.edu.au/wsc/docs/tech-manual-v3.pdf
#zoo real data #velocity 0.6 #depth 0.4572 #sunrise 6:00AM #sunset 8:15PM #96 is the number of intervals in a day for fifteen minute data data(DOTemp) prod.20 <- EcosystemProduction.20(0.6, 0.4572, DOTemp[,1], DOTemp[,2], "8/18/70", "00:00:00", "23:45:00", "06:00:00", "20:15:00", 96) prod.20 #compare to EcosystemProduction prod <- EcosystemProduction(0.6, 0.4572, DOTemp[,1], DOTemp[,2], "8/18/70", "00:00:00", "23:45:00", "06:00:00", "20:15:00", 96) #notice that NDM stays the same...# prod prod.20