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 ...
206 32 CHAPTER 5. DYNAMICS OF THE OCEAN 5.11 Potential Vorticity (non-linear) tel-00545911, version 1 - 13 Dec 2010 Similar calculations for the non-linear equations (5.38) – (5.40) lead to ( ) d ζ + f = 0. (5.57) dt H + η This means that every fluid parcel, or in this case every fluid column, conserves its potential vorticity Q sw = (ζ+f)/(H+η), that is, potential vorticity is transported by the two dimensional flow. The part f/(H + η) which does not depend explicitly on the velocity is called planetary potential vorticity, while ζ/(H + η) is called the dynamical part. Example: Eddy over sea mount. Exercise 30: derive eq. (5.57). Exercise 31: If you make the assumption of linearity, can you obtain Q lin sw from Q sw ? Exercise 32: The moment of inertia of a cylinder of mass m, radius r and height H is given by I = mr 2 /2 the angular momentum is given by L = Iω. If a cylinder stretches or flatness without any forces acting from the outside its angular momentum is conserved: H Show, that during this process ω/H is conserved. The previous exercise demonstrates, that the conservation of potential vorticity is nothing else than the conservation of angular momentum applied to a continuum in a rotating frame. 5.12 The Beta-plane H So far we supposed the earth to be flat! The dominant difference, induced by the spherical shape of the earth, for the large scale ocean dynamics, at mid- and low latitudes, is the change of the (locally) vertical component of the rotation vector. For the large scale circulation a major source of departure from geostrophy is the variation of f = f 0 with latitude. So far we have considered f to be constant we will now approximate it by f = f 0 + βy, where f 0 = 2|Ω| sin(θ 0 ) and β = 2(|Ω|/R) cos(θ 0 ) are constant where R is the radius of the earth, it takes its maximum value β max = 2.3 × 10 −11 m −1 s −1 at the equator. The geometry with a linearly changing Coriolis parameter is called the β-plane. The change of f with latitude, the so called β-effect, compares to the effect of constant rotation for phenomena with horizontal extension L ≈ f/β = R tan(θ) or larger. Replacing f by f 0 + βy in equations (5.41), (5.42) and (5.43): ∂ t (∂ x v − ∂ y u) + f(∂ x u + ∂ y v) + βv = 0 (5.58)
207 5.13. A FEW WORDS ABOUT WAVES 33 leading to: ∂ t ζ − f∂ z w + βv = 0. (5.59) which states, that the vorticity ζ is changed by the vertical gradient of the vertical velocity (vortex stretching) and the planetary change, due to β and the latitudinal velocity). Exercise 33: what is the sign of f and β on the northern and southern hemisphere, respectively? Exercise 34: what is the value of f and β on the equator, north and south pole? Exercise 35: discuss the importance of f and β for equatorial dynamics. 5.13 A few Words About Waves tel-00545911, version 1 - 13 Dec 2010 As mentioned in the preface we do not explicitly consider wave dynamics in this introductory text. I like to make, nevertheless, some “hand waving” arguments about the role of waves in the ocean. The ocean and atmosphere dynamics at large scales are always close to a geostrophic balance. There are, however, different sources of perturbations of the geostrophically balanced state: • variation of the Coriolis parameter f • non-linearity • topography • instability • forcing (boundary conditions) • friction • other physical processes (convection, ..) As the geostrophic adjustment process happens on a much faster time scale than the geostrophic dynamics, these perturbations lead not so much to a departure from the geostrophic state but more to its slow evolution. In this adjustment process, discussed in section 5.10, (gravity) waves play an important role. It is an important part of research in geophysical fluid dynamics (GFD) (DFG, en français) to find equations that reflect the slow evolution of the geostrophic state, without explicitly resolving the geostrophic adjustment process. Such equations are called balanced equations, and are based on the evolution of PV. The best known system of balanced equations are the quasi-geostrophic equations. The problem in constructing such equations is how to calculate the velocity field from PV, a process usually referred to as inversion. The fast surface gravity waves influenced by rotation, Poincaré waves have no PV signature and thus do not appear in the balanced equations, which leads to a large simplification for analytical and numerical calculations. Balanced equations such rely on the assumption that the ocean dynamics can be separated into fast wave motion and slow vortical motion with no or negligible interactions between the two. They describe the dynamics on time-scales longer than the period of gravity waves, typically several inertial periods f −1 . The balanced equations are not valid when approaching the equator, as f −1 → ∞. The dynamics described by the balanced equations is said to represent the slow dynamics or to evolve on the slow manifold .
- Page 161 and 162: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 163 and 164: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 165 and 166: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 167 and 168: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 169 and 170: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 171 and 172: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 173 and 174: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 175 and 176: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 177 and 178: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 179 and 180: Quatrième partie tel-00545911, ver
- Page 181 and 182: 175 tel-00545911, version 1 - 13 De
- Page 183 and 184: 177 Contents 1 Preface 5 tel-005459
- Page 185 and 186: 179 Chapter 1 Preface tel-00545911,
- Page 187 and 188: 181 Chapter 2 Observing the Ocean t
- Page 189 and 190: 183 Chapter 3 Physical properties o
- Page 191 and 192: 185 3.3. θ-S DIAGRAMS 11 3.3 θ-S
- Page 193 and 194: 187 3.6. HEAT CAPACITY 13 tel-00545
- Page 195 and 196: 189 3.7. CONSERVATIVE PROPERTIES 15
- Page 197 and 198: 191 Chapter 4 Surface fluxes, the f
- Page 199 and 200: 193 4.2. FRESH WATER FLUX 19 water.
- Page 201 and 202: 195 Chapter 5 Dynamics of the Ocean
- Page 203 and 204: 197 5.2. THE LINEARIZED ONE DIMENSI
- Page 205 and 206: 199 5.4. TWO DIMENSIONAL STATIONARY
- Page 207 and 208: 201 5.6. THE CORIOLIS FORCE 27 Whic
- Page 209 and 210: 203 5.8. GEOSTROPHIC EQUILIBRIUM 29
- Page 211: 205 5.10. LINEAR POTENTIAL VORTICIT
- Page 215 and 216: 209 Chapter 6 Gyre Circulation tel-
- Page 217 and 218: 211 6.1. SVERDRUP DYNAMICS IN THE S
- Page 219 and 220: 213 6.2. THE EKMAN LAYER 39 In the
- Page 221 and 222: 215 6.3. SVERDRUP DYNAMICS IN THE S
- Page 223 and 224: 217 Chapter 7 Multi-Layer Ocean dyn
- Page 225 and 226: 219 7.3. GEOSTROPHY IN A MULTI-LAYE
- Page 227 and 228: 221 7.5. EDDIES, BAROCLINIC INSTABI
- Page 229 and 230: 223 Chapter 8 Equatorial Dynamics t
- Page 231 and 232: 225 Chapter 9 Abyssal and Overturni
- Page 233 and 234: 227 9.2. MULTIPLE EQUILIBRIA OF THE
- Page 235 and 236: 229 9.3. WHAT DRIVES THE THERMOHALI
- Page 237 and 238: 231 Chapter 10 Penetration of Surfa
- Page 239 and 240: 233 10.2. TURBULENT TRANSPORT 59 If
- Page 241 and 242: 235 10.5. ENTRAINMENT 61 instabilit
- Page 243 and 244: 237 Chapter 11 Solution of Exercise
- Page 245 and 246: 239 65 Exercise 32: The moment of i
- Page 247 and 248: 241 INDEX 67 Transport stream-funct
- Page 249 and 250: Annexe A Attestation de reussite au
- Page 251 and 252: Annexe B Rapports du jury et des ra
- Page 253 and 254: tel-00545911, version 1 - 13 Dec 20
- Page 255 and 256: tel-00545911, version 1 - 13 Dec 20
- Page 257 and 258: 251 utilisés avec pertinence. Sur
- Page 259 and 260: 253
- Page 261 and 262: tel-00545911, version 1 - 13 Dec 20
206<br />
32 CHAPTER 5. DYNAMICS OF THE OCEAN<br />
5.11 Potential Vorticity (non-linear)<br />
tel-00545911, version 1 - 13 Dec 2010<br />
Similar calculations for the non-linear equations (5.38) – (5.40) lead to<br />
( )<br />
d ζ + f<br />
= 0. (5.57)<br />
dt H + η<br />
This means that every fluid <strong>par</strong>cel, or in this case every fluid column, conserves its potential<br />
vorticity Q sw = (ζ+f)/(H+η), that is, potential vorticity is transported by the two dimensional<br />
flow. The <strong>par</strong>t f/(H + η) which does not <strong>de</strong>pend explicitly on the velocity is called plan<strong>et</strong>ary<br />
potential vorticity, while ζ/(H + η) is called the dynamical <strong>par</strong>t.<br />
Example: Eddy over sea mount.<br />
Exercise 30: <strong>de</strong>rive eq. (5.57).<br />
Exercise 31: If you make the assumption of linearity, can you obtain Q lin<br />
sw from Q sw ?<br />
Exercise 32: The moment of inertia of a cylin<strong>de</strong>r of mass m, radius r and height H is given<br />
by I = mr 2 /2 the angular momentum is given by L = Iω. If a cylin<strong>de</strong>r str<strong>et</strong>ches or flatness<br />
without any forces acting from the outsi<strong>de</strong> its angular momentum is conserved:<br />
H<br />
Show, that during this process ω/H is conserved.<br />
The previous exercise <strong>de</strong>monstrates, that the conservation of potential vorticity is nothing<br />
else than the conservation of angular momentum applied to a continuum in a rotating frame.<br />
5.12 The B<strong>et</strong>a-plane<br />
H<br />
So far we supposed the earth to be flat! The dominant difference, induced by the spherical<br />
shape of the earth, for the large scale ocean dynamics, at mid- and low latitu<strong>de</strong>s, is the change<br />
of the (locally) vertical component of the rotation vector.<br />
For the large scale circulation a major source of <strong>de</strong><strong>par</strong>ture from geostrophy is the variation<br />
of f = f 0 with latitu<strong>de</strong>. So far we have consi<strong>de</strong>red f to be constant we will now approximate it<br />
by f = f 0 + βy, where f 0 = 2|Ω| sin(θ 0 ) and β = 2(|Ω|/R) cos(θ 0 ) are constant where R is the<br />
radius of the earth, it takes its maximum value β max = 2.3 × 10 −11 m −1 s −1 at the equator. The<br />
geom<strong>et</strong>ry with a linearly changing Coriolis <strong>par</strong>am<strong>et</strong>er is called the β-plane. The change of f<br />
with latitu<strong>de</strong>, the so called β-effect, com<strong>par</strong>es to the effect of constant rotation for phenomena<br />
with horizontal extension L ≈ f/β = R tan(θ) or larger.<br />
Replacing f by f 0 + βy in equations (5.41), (5.42) and (5.43):<br />
∂ t (∂ x v − ∂ y u) + f(∂ x u + ∂ y v) + βv = 0 (5.58)