Introduction to surface ocean modelling

1. SOPRAN GOTM School Warnemünde: 10.-11.09.07 Hans Burchard Baltic Sea Research Institute Warnemünde, Germany Introduction to surface ocean modelling

2. Surface ocean processes

3. How to model all this ? Surface ocean physical processes (Thorpe, 1995)

4. 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, …

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

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

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

8. Spectral Considerations From Schatzmann

9. 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 …

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

11. Convenient approximations: • Hydrostatic approximation (vertical velocity • dynamically irrelevant) • One-dimensional approximation (horizontal homogeneity, • far away from coasts and fronts)

12. 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:

13. Wind stress vector Latent heat flux Sensitive heat flux

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

15. Have we modelled all this ?

16. Example: Station P in Northern Pacific

17. Example: Station P in Northern Pacific

18. Example: Station P in Northern Pacific Time series of SST (observed and simulated)