180 likes | 191 Views
This article provides an introduction to surface ocean modeling, focusing on the physical processes and dynamic equations involved. It explores the Reynolds decomposition method and discusses turbulent fluxes, TKE equation, spectral considerations, and various turbulence closure models. The article also highlights the importance of initial and boundary conditions in modeling surface ocean properties.
E N D
SOPRAN GOTM School Warnemünde: 10.-11.09.07 Hans Burchard Baltic Sea Research Institute Warnemünde, Germany Introduction to surface ocean modelling
How to model all this ? Surface ocean physical processes (Thorpe, 1995)
Basic physical principles: • conservation of volume (incompressibility) • conservation of mass (water, salt, …) • conservation of momentum (velocity) • conservation of angular momentum • conservation of total energy (mechanical & thermodynamic) • plus • material laws for water (viscosity, …) • gives • Dynamic equations for momentum, heat, salt, …
These dynamic equation are valid on all scales, and all these scale are relevant. Problem: in numerical models, we cannot resolve from millimeter to kilometer. Therefore, equations are statistically separated into mean (expected) and fluctuating part (Reynolds decomposition):
The Reynolds decomposition allows to derive dynamic equations for mean-flow quantaties, but … … new (unknown) terms are introduced, the turbulent fluxes: Eddy viscosity Reynolds stresses Eddy diffusivity Turbulent salt flux Turbulent heat flux
The TKE equation Reynolds stresses cause loss of kinetic energy from mean flow, which is a source of turbulent kinetic energy (TKE, k). TKE may is produced as large eddy sizes and is dissipated into heat by small eddy sizes at rate e (dissipation rate). At stable stratification, TKE is converted into potential energy (vertical mixing, depening of mixed layer). Unstable stratification converts potential energy into TKE (convective mixing).
Spectral Considerations From Schatzmann
How to calculate the eddy viscosity / eddy diffusivity ? Turbulent macro length scale The well-known k-e model uses dynamic equations for the TKE and its dissipation rate. There are however many other models in use …
Bulk models Differential models Kraus-Turner type models Empirical models Statistical models KPP Ri number depending models Flow depending models Algebraic stress models Full Reynolds stress models Treatment of algebraic stresses Treatment of TKE and length scale Non-equilibrium models Quasi-equilibrium models Zero-equation models One-equation models Two-equation models Generic length scale Mixing length formulations Blackadar-type length scale Gaspar et al. (1990) type models k-w k-e MY Tree of turbulence closure models (extention of Haidvogel & Beckmann, 1999)
Convenient approximations: • Hydrostatic approximation (vertical velocity • dynamically irrelevant) • One-dimensional approximation (horizontal homogeneity, • far away from coasts and fronts)
The dynamic equations for momentum, temperature, etc., are PDEs (Partial Differential Equations), and therefore need initial and boundary conditions. Initial conditions are either from observations, idealised, or simply set to dummy values (because they may be forgotten after a while). Surface boundary conditions for physical properties come from atmospheric conditions:
Wind stress vector Latent heat flux Sensitive heat flux
Simulated short-wave radiation profile in water, I(z) Surface radiation Attenuation lengths Weighting Bio-shading Short-wave radiation in water The local heating depends on the vertical gradient of I(z).
Example: Station P in Northern Pacific Time series of SST (observed and simulated)