lake.number: Calculate Lake Number

Description Usage Arguments Value References See Also Examples

View source: R/lake.number.R

Description

The Lake Number, defined by Imberger and Patterson (1990), has been used to describe processes relevant to the internal mixing of lakes induced by wind forcings. Lower values of Lake Number represent a higher potential for increased diapycnal mixing, which increases the vertical flux of mass and energy across the metalimnion through the action of non-linear internal waves. Lake Number is a dimensionless index.

Lake Number has been used, for example, to estimate the flux of oxygen across the thermocline in a small lake (Robertson and Imberger, 1994), and to explain the magnitude of the vertical flux of ammonium in a lake (Romero et al., 1998). In Imberger and Patterson (1990), Lake Number was defined as Ln = (g * St * (zm - zT)) / (rho_0 * u*^2 * A0^3/2 * (zm - zg)) with all values referenced from the lake bottom, e.g., zm being the height of the water level, zT the height of metalimnion and zg the height of center volume. Our calculations assume that the reference is at the lake surface, therefore: height of metalimnion becomes metalimnion depth (average of meta top and bot): (zm - zT) –> (metaT + metaB)/2 height of center of volume depth becomes center of volume depth Zcv: (zm - zg) –> Zcv Further, we note that in that original work St was defined as St = int (z - zg) A(z) rho(z) dz and rLakeAnalyzer defines St as St = g/A0 int (z - zg) rho(z) dz Therefore, we calculate St_uC = St * Ao / g

Usage

1
lake.number(bthA, bthD, uStar, St, metaT, metaB, averageHypoDense)

Arguments

bthA

a numeric vector of cross sectional areas (m2) corresponding to bthD depths, hypsographic areas

bthD

a numeric vector of depths (m) which correspond to areal measures in bthA, hypsographic depths

uStar

a numeric array of u* (m/s), water friction velocity due to wind stress

St

a numeric array of Schmidt stability (J/m2), as defined by Idso 1973

metaT

a numeric array of the top of the metalimnion depth (m from the surface)

metaB

a numeric array of the bottom of the metalimnion depth (m from the surface)

averageHypoDense

a numeric array of the average density of the hypolimnion (kg/m3)

Value

A numeric vector of Lake Number [dimensionless]

References

Imberger, J., Patterson, J.C., 1990. Physical limnology. Advances in Applied Mechanics 27, 303-475.

Idso, S.B., 1973. On the concept of lake stability. Limnology and Oceanography 18, 681-683.

See Also

ts.lake.number wedderburn.number

Examples

1
2
3
4
5
6
7
8
9
	bthA	<-	c(1000,900,864,820,200,10)
	bthD	<-	c(0,2.3,2.5,4.2,5.8,7)
	uStar	<-	c(0.0032,0.0024)
	St	<-	c(140,153)
	metaT	<-	c(1.34,1.54)
	metaB	<-	c(4.32,4.33)
	averageHypoDense	<-	c(999.3,999.32)
	cat('Lake Number for input vector is: ')
	cat(lake.number( bthA, bthD, uStar, St, metaT, metaB, averageHypoDense) )

GLEON/rLakeAnalyzer documentation built on May 3, 2021, 1:32 p.m.