schmidt.stability {rLakeAnalyzer}

Calculate the Schmidt stability

Description

Schmidt stability, or the resistance to mechanical mixing due to the potential energy inherent in the stratification of the water column.

Schmidt stability was first defined by Schmidt (1928) and later modified by Hutchinson (1957). This stability index was formalized by Idso (1973) to reduce the effects of lake volume on the calculation (resulting in a mixing energy requirement per unit area).

Usage

schmidt.stability(wtr, depths, bthA, bthD, sal = 0)

Arguments

wtr

a numeric vector of water temperature in degrees C

depths

a numeric vector corresponding to the depths (in m) of the wtr measurements

bthA

a numeric vector of cross sectional areas (m^2) corresponding to bthD depths

bthD

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

sal

a numeric vector of salinity in Practical Salinity Scale units

Value

a numeric vector of Schmidt stability (J/m^2)

References

Schmidt, W., 1928. Ueber Temperatur and Stabilitaetsverhaltnisse von Seen. Geo- graphiska Annaler 10, 145-177.

Hutchinson, G.E., 1957. A Treatise on Limnology, vol. 1. John Wiley & Sons, Inc., New York.

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

See Also

ts.schmidt.stability lake.number wedderburn.number

Examples

	bthA	<-	c(1000,900,864,820,200,10)
	bthD	<-	c(0,2.3,2.5,4.2,5.8,7)
	
	wtr	<-	c(28,27,26.4,26,25.4,24,23.3)
	depths	<-	c(0,1,2,3,4,5,6)
	
	cat('Schmidt stability for input is: ')
	cat(schmidt.stability(wtr, depths, bthA, bthD))