schmidt.stability: Calculate the Schmidt stability
In rLakeAnalyzer: Lake Physics Tools

Description Usage Arguments Value References See Also Examples

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).

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

`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

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

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.

ts.schmidt.stability lake.number wedderburn.number

	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))