Etudes et évaluation de processus océaniques par des hiérarchies ...

Etudes et évaluation de processus océaniques par des hiérarchies ... Etudes et évaluation de processus océaniques par des hiérarchies ...

tel.archives.ouvertes.fr
from tel.archives.ouvertes.fr More from this publisher
27.12.2013 Views

202 28 CHAPTER 5. DYNAMICS OF THE OCEAN In such situation the last term in equation (5.34) has to be dropped. Exercise 17: suppose (x,y) = (R cos(ωt),Rsin(ωt)) for a fluid particle in the fluid corresponding to fig. 5.4, without exterior forces acting. Calculate ω. Such kind of motion, that is, anti-cyclonic rotation with a period which is half the local rotation period, is indeed often observed in oceanic and atmospheric motion and is called inertial oscillation and their frequency (−f) is called inertial frequency. tel-00545911, version 1 - 13 Dec 2010 ✛ ✲ ✻ Ω Figure 5.4: Cylinder in rotation with a free surface; two fluid particles with centrifugal force (green) and pressure gradient force (red). On earth the same thing happens, the centrifugal force changes the geopotential of the earth, flattens it a little bit, makes it an ellipsoid. 5.7 The Shallow Water Equations in a Rotating Frame ✛ ✲ If we take the results from the previous section we see that we only have to add the Coriolis force in the shallow water equations to obtain the shallow water equations in a rotating frame: ∂ t u + u∂ x u + v∂ y u − fv + g∂ x η = ν∇ 2 u (5.38) ∂ t v + u∂ x v + v∂ y v + fu + g∂ y η = ν∇ 2 v (5.39) ∂ t η + ∂ x [(H + η)u] + ∂ y [(H + η)v] = 0 (5.40) +boundary conditions. The nonlinear terms can be neglected if Rossby number ǫ = u/(fL) is small. The Rossby number compares the distance a fluid particle has traveled in the time f −1 to the length scale of the phenomenon considered. The linear (small Rossby number) version the shallow water equations in a rotating frame is: ∂ t u − fv + g∂ x η = 0 (5.41)

203 5.8. GEOSTROPHIC EQUILIBRIUM 29 ∂ t v + fu + g∂ y η = 0 (5.42) ∂ t η + H(∂ x u + ∂ y v) = 0 (5.43) +boundary conditions. Important: When approaching the equator f tends to zero, so rotation no longer dominates and most of the considerations following are not applicable. Equatorial dynamics is different! In equations (5.41) and (5.42) we have neglected the viscous term which can be safely done as ν(sea water)≈ 10 −6 m 2 /s. Exercise 18: What is the Rossby number of the basin wide circulation in the North Atlantic when u = 10 −1 m/s? What is the Rossby number of a Gulf Stream eddy when u = 1m/s and the radius R = 30km ? tel-00545911, version 1 - 13 Dec 2010 5.8 Geostrophic Equilibrium Large-scale ocean currents usually change on a time scale much larger than f −1 and are thus often well approximated by the stationary versions of eqs. (5.41) – (5.43) which are, fv = g∂ x η (5.44) −fu = g∂ y η (5.45) this is called the geostrophic equilibrium. For a flow in geostrophic equilibrium all variables can be expressed in terms of the free surface, you can easily calculate that the vorticity is given by ζ = (g/f)∇ 2 η. Note that to every function η(x,y) there is a unique flow in geostrophic equilibrium associated to it. In the case with no rotation (f = 0) the stationary solutions of the linearised equations have η = 0 and ∂ x u + ∂ y v = 0. In the case without rotation η(x,y) does not determine the flow. Exercise 19: What happened to the stationary version of equation (5.43) ? Exercise 20: Across the Gulf Stream, which is about 100km wide there is a height difference of approx. 1m. What is the corresponding geostrophic speed of the Gulf Stream. Exercise 21: In a sea surface height (SSH) map, how can you distinguish cyclones from anti-cyclones? What happens on the southern hemisphere? The function Ψ = (gH/f)η is called the geostrophic transport stream-function as ∂ x Ψ = Hv and ∂ y Ψ = −Hu which means that: (i) isolines of Ψ, and also of η, are stream-lines of the geostrophic velocity field, and (ii) Ψ(B) − Ψ(A) is the transport that passes between points A and B. In oceanography transport is usually measured in Sverdrup (1Sv =10 6 m 3 /s), which corresponds to a cube of water of side length 100 meters passing in 1 second. When the stationarity assumption is not made the eqs. (5.41) – (5.43) can be used to derive an equation for η only where u and v can be derived from η. This leads to: Exercise 22: derive eqs. (5.46) – (5.48). ∂ t [ ∂tt η + f 2 η − gH∇ 2 η ] = 0 (5.46) ∂ tt u + f 2 u = −g(∂ tx η + f∂ y η) (5.47) ∂ tt v + f 2 v = −g(∂ ty η − f∂ x η) (5.48) Exercise 23: show that the geostrophic equilibrium is a solution of eqs. (5.46) – (5.48). Exercise 24: show that the only stationary solution of eqs. (5.46) – (5.48) is geostrophic equilibrium.

203<br />

5.8. GEOSTROPHIC EQUILIBRIUM 29<br />

∂ t v + fu + g∂ y η = 0 (5.42)<br />

∂ t η + H(∂ x u + ∂ y v) = 0 (5.43)<br />

+boundary conditions.<br />

Important: When approaching the equator f tends to zero, so rotation no longer dominates<br />

and most of the consi<strong>de</strong>rations following are not applicable. Equatorial dynamics is different!<br />

In equations (5.41) and (5.42) we have neglected the viscous term which can be safely done<br />

as ν(sea water)≈ 10 −6 m 2 /s.<br />

Exercise 18: What is the Rossby number of the basin wi<strong>de</strong> circulation in the North Atlantic<br />

when u = 10 −1 m/s? What is the Rossby number of a Gulf Stream eddy when u = 1m/s and<br />

the radius R = 30km ?<br />

tel-00545911, version 1 - 13 Dec 2010<br />

5.8 Geostrophic Equilibrium<br />

Large-scale ocean currents usually change on a time scale much larger than f −1 and are thus<br />

often well approximated by the stationary versions of eqs. (5.41) – (5.43) which are,<br />

fv = g∂ x η (5.44)<br />

−fu = g∂ y η (5.45)<br />

this is called the geostrophic equilibrium. For a flow in geostrophic equilibrium all variables can<br />

be expressed in terms of the free surface, you can easily calculate that the vorticity is given<br />

by ζ = (g/f)∇ 2 η. Note that to every function η(x,y) there is a unique flow in geostrophic<br />

equilibrium associated to it. In the case with no rotation (f = 0) the stationary solutions of<br />

the linearised equations have η = 0 and ∂ x u + ∂ y v = 0. In the case without rotation η(x,y)<br />

does not d<strong>et</strong>ermine the flow.<br />

Exercise 19: What happened to the stationary version of equation (5.43) ?<br />

Exercise 20: Across the Gulf Stream, which is about 100km wi<strong>de</strong> there is a height difference<br />

of approx. 1m. What is the corresponding geostrophic speed of the Gulf Stream.<br />

Exercise 21: In a sea surface height (SSH) map, how can you distinguish cyclones from<br />

anti-cyclones? What happens on the southern hemisphere?<br />

The function Ψ = (gH/f)η is called the geostrophic transport stream-function as ∂ x Ψ = Hv<br />

and ∂ y Ψ = −Hu which means that: (i) isolines of Ψ, and also of η, are stream-lines of the<br />

geostrophic velocity field, and (ii) Ψ(B) − Ψ(A) is the transport that passes b<strong>et</strong>ween points<br />

A and B. In oceanography transport is usually measured in Sverdrup (1Sv =10 6 m 3 /s), which<br />

corresponds to a cube of water of si<strong>de</strong> length 100 m<strong>et</strong>ers passing in 1 second.<br />

When the stationarity assumption is not ma<strong>de</strong> the eqs. (5.41) – (5.43) can be used to <strong>de</strong>rive<br />

an equation for η only where u and v can be <strong>de</strong>rived from η. This leads to:<br />

Exercise 22: <strong>de</strong>rive eqs. (5.46) – (5.48).<br />

∂ t<br />

[<br />

∂tt η + f 2 η − gH∇ 2 η ] = 0 (5.46)<br />

∂ tt u + f 2 u = −g(∂ tx η + f∂ y η) (5.47)<br />

∂ tt v + f 2 v = −g(∂ ty η − f∂ x η) (5.48)<br />

Exercise 23: show that the geostrophic equilibrium is a solution of eqs. (5.46) – (5.48).<br />

Exercise 24: show that the only stationary solution of eqs. (5.46) – (5.48) is geostrophic<br />

equilibrium.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!