Equation of State 1. Physical Approach. 2. Polynomial fits. • Tumlirz equation • Eckart’s method (Eckart 1958) • Wright’s equation of state (Wright 1996) *nonlinear least square fit to IES80 over the full parameter range. *use potential temperature rather than temperature. • IES80 (UNESCO 1981, Gill 1982) 0 0 0 0 0 p p p p p p • Range of Parameters • Ref. Formula • Approximation Polynomial Example (Friedrich and Levitus 1972) : Knudsen-Ekman formula (Bradshaw and Schleicher, 1970) % 38 % 30 , 30 2 S C T C o o ST C T C S C T C C S T 5 2 4 3 2 1 ) , (
16
Embed
Equation of State 1. Physical Approach. 2. Polynomial fits. Tumlirz equation Eckart’s method (Eckart 1958) Wright’s equation of state (Wright 1996) *nonlinear.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
• Wright’s equation of state (Wright 1996) *nonlinear least square fit to IES80 over the full parameter range. *use potential temperature rather than temperature.
• IES80 (UNESCO 1981, Gill 1982)
00
000 pp
pppp
• Range of Parameters
• Ref. Formula
• Approximation Polynomial
Example (Friedrich and Levitus 1972) :
Knudsen-Ekman formula (Bradshaw and Schleicher, 1970)
%38%30,302 SCTC oo
STCTCSCTCCST 52
4321),(
Effective Calculation of Sea Water Density 1. Sanderson, Dietrich, Stilgoe 2001.
2. Accuracy and Computational Advantage over UES80 and Wright’s E.O.S based on the choice of ocean model and problem being solved.
In a 3D ocean model solving a deep convection problem: 7% drop and 15% drop of computational cost for nonhydro. and hydro. cases relative to UES80. 3. Computational cost is high because of calculation for coefficients. *updating only many time steps. *In deep ocean, where temperature and Salinity are nearly constant. Taylor expansion about the time mean state.
Simple Wright’s UNESCO Linear Nonlinear1 Nonlinear2
2FP 21FP 85FP+1SB 4FP+1SB 6FP+1SB 8FP+1SB
Single Double Single Single Single Single
1 2 3 4 5 6
Table 1. Brief comparison of different E.O.S. employed in Ocean Model
• Moreover, they found κ change significantly with depth.
• By fitting surfaces to the observed density field,
1. In the North Atlantic, between 800 m and 2000 m, κ= 1.0x10-5 m2/s (Olbers etal., 1985)
2. In the Southern Ocean, between 100 m and 2500 m, κ= 1.0x10-4 m2/s (with some regions reach much larger values, 1.0x10-3 m2/s )
3. Direct Estimates, by bulk fluid properties.
[Wunsch & Ferrari, 2004]
Interior Vertical Eddy Diffusivity, “κ ” • Microstructure Measurements. (Osborn & Cox, 1972; Osborn, 1980) “Shears”. • Mixing in the Open ocean area V.S. Boundaries.
1. In the open ocean upper 1000 m, below the mixed layer, κ can’t exceed 1.0x10-5 m2/s. (Gregg, 1987)
2. Use tracer-releasing method. (Ledwell et al., 1998,2000).3. For supporting the observed circulation, one tenth of the value might be adopted.4. More mixing happens at boundaries. (Ledwell & Hickey, 1995)
Mixing of the Upper Ocean
Abyssal Mixing
1. High-accuracy measurements in 1990’s.2. κ< O(10-5 m2/s) over abyssal plains and other simple structures. (Toole er
al., 1994; Polzin et al., 1997) Contradiction to dependence on N(z).3. κ increases with depth (Ledwell et al., 2000); over continental slopes
(Polzin, 2002); and in some abyssal passages.
[Sanderson et al., 2001]
Deep Ocean Diapycnal Diffusivity ( 950-1450 m )1. Hibiya, Nagasawa, and Niwa, 2006, Hibiya and Nagasawa, 2004.
2. Expendable Current Profiler (XCP) -> An empirical relationship between estimated diapycnal diffusivity and the numerically predicted, available energy density of the semidiurnal internal tide. -> use numerically predicted energy density incorporating into the relationship to obtain the global distribution of diapycnal diffusivity in the thermohaline.
3. The estimated diapycnal diffusivities are dependent on the latitude.
Values of the estimated diapycnal diffusivities:
• Near the Hawaiian Ridge(~25oN) & Izu-Ogasawara Ridge (~28oN) : 1.5x10-4 m2s-1
• As the latitude exceeds 30oN : 0.2x10-4 m2s-1
• Near the Aleutian Ridge (~52oN) & the Emperor Seamount (~38oN) : 0.1x10-4 m2s-
1
Vertical Mixing Parameterization (Only in Mixed layer)Bulk Formula
1. Kraus and Turner, 1967; Price, et al.,1986; Chen, et al.,1994
2. Considering the solar radiation, energy input from winds, and dissipation within one layer.
3. Compute layer depth and temperature.
Vertical Mixing Formula
1. Pacanowski and Philander, 1981 (PP); Mellor and Yamada, 1982; Large et al., 1994 (KPP); Canoto and Dubovikov, 1996.
2. Mixing length theory + turbulence closure. Partly based on the Richardson’s concept.
3. Employed in our ocean model.
• Meso-scale eddies, Gent and McWilliams, 1990 (GM90)• Parameterization of Eddy Fluxes near Ocean Boundaries, Ferrari, 2007• Suppression of Eddy Diffusivity across Jets in the Southern Ocean, Ferrari, 2010
KPP. (Large, 1994) Closure Model. (Canuto, 2000)
KPPPresent ModelLevitus data
Parameterizations Employed in ModelsVertical Horizontal Comments