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

136 CHAPITRE 4. ETUDES DE PROCESSUS OCÉANOGRAPHIQUES 556 Ocean Dynamics (2009) 59:551–563 tel-00545911, version 1 - 13 Dec 2010 note that the thickness of the friction layer is about two times the Ekman layer thickness δ = √ 2ν v /f, confirming the results of Shapiro and Hill (1997). The friction layer typically contains a bottom and an interfacial Ekman layer, which match at its interior. The Ekman layer extends through the entire width of the gravity current, friction layer and vein, and is most conspicuous when looking at the u component of the velocity vector, which is strongly negative (upslope is the positive direction) in the Ekman layer and small in the interior of the vein. Interfacial Ekman layers between the surrounding fluid and the gravity current are also observed, but the gradients involved are smaller than at the ocean floor. In the vein away from the boundaries, the current is close to a geostrophic equilibrium, where the pressure gradient due to the interface of the gravity current is balanced by the Coriolis force due to the velocity of the vein. Near the boundaries, the speed of the gravity current is reduced and the Coriolis force can no longer balance the pressure gradient, which leads to fluid moving down the pressure gradient in the friction layer, that is, on average, down the slope. The variation in the Ekman layer transport creates convergence and divergence and the gravity current water is vertically pumped out and into the vein at the interface of the vein and the Ekman layer (see Fig. 2). These vertical velocities then affect the dynamics in the vein by vortex stretching in the same way as the vertical velocities at the wind induced Ekman layer influence the geostrophic Sverdrup dynamics in the interior ocean (see, e.g. Pedlosky 1998). The influence of friction on the dynamics of the vein is, thus, via the vertical velocities at the base of the vein. Figure 1 also shows that the vein detrains via the friction layer. This detrainment, or Ekman drainage, depends on the (turbulent or eddy) vertical viscosity, which is a key parameter in the gravity current dynamics. In our calculations, the vertical viscosity is constant (see Subsection 2.3); in nature, the vertical eddy viscosity is, however, not constant but varies in space and time. The Ekman drainage has already been found to play an important role in laboratory experiments (Lane-Serf and Baines 2000). In the interior of the vein, inertial oscillations are superposing the slow evolution of the gravity current in all calculations. In the calculations with higher values of the reduced gravity (T > 0.75 K), nonlinear behaviour is observed as can be seen in Fig. 3, where the vertical velocity outside and inside the gravity current shows an involved spatio-temporal variability. The sparse horizontal resolution of our calculations does not allow for the representation of turbulence in the unstratified interior of the gravity current. 3.2 Velocity, angle of descent and broadening Before proceeding with a qualitative analysis of the data, I have to give a precise definition of the gravity current: water parcels having at least half the initial density anomaly are said to be within the gravity current. I verified that the results presented here vary only slightly when this threshold is varied around the value of one-half. The total volume of the gravity current vein friction layer Ekman layer Fig. 2 Schematic of a cross section through a gravity current on an inclined plane as considered in this manuscript. The total gravity current consists of two distinct parts, the vein (the fluid under the blue line) and the friction layer (the fluid under the red line). In both parts, there is an Ekman layer close to the inclined plane (under the thin black line). The green line separates the concave (∂ xx h > 0) from the convex (∂ xx h < 0) part of the vein. The small green arrows show the direction of vertical velocity out of or into the Ekman layer

4.7. ON THE BASIC STRUCTURE OF OCEANIC GRAVITY CURRENTS 137 Ocean Dynamics (2009) 59:551–563 557 Fig. 3 Structure and vertical velocity field at 24 h in exp. G15. The temperature anomaly, with respect to the ambient water, is shown (given in Kelvin). Isolines of the w component of the velocity are given in black for positive values (upward) and red for negative values; the contour lines shown are ±10 −5 , 3 · 10 −5 , 10 −4 , 3 · 10 −4 , 10 −3 , 3 · 10 −3 ; the logarithmic scale was chosen due to the large intermittency of the vertical velocity. Vertical velocities outside the vein are due to wave-like motion of the interface tel-00545911, version 1 - 13 Dec 2010 shows only a small variation during the evolution of the gravity current. This shows that there is no substantial entrainment or detrainment for the total gravity current. The density anomaly and the along-slope velocity (y component) both decrease in the experiments. To leading order, the dynamics is supposed to be given by the geostrophic equilibrium as discussed in Subsection 2.1. To verify this, I calculate the along-slope velocity component in the vein divided by the geostrophic speed of its reduced gravity, given by: v G = (g ′ /f) tan α, for each density class within the vein. When the average for all parcels in the vein which are at least 20 m above ground is analysed, the average along-slope velocity decreases to about 95% the geostrophic velocity after 2 days. This shows that the gravity current progresses along the slope at almost geostrophic speed, in agreement with previous findings of Price and Baringer (1994) and others. Friction leads to a downslope transport of gravity current water. The dynamics of the vein differs completely from the dynamics of the total gravity current, and so does their rate of descent and their spreading, as can be verified in Fig. 4. The spread and descent of the gravity current, which includes the friction layer, is much larger than that of the vein alone. Figure 4 shows that the depth of the upper bound of the vein stays almost constant, a feature that is well documented (see Price and Baringer 1994 and references therein), while the downslope sides of the vein and the friction layer descend. At the upslope side, the y component of the velocity vector is negative, as noted above, the fluid in the Ekman layer should, thus, move upslope. This leads to an arrested Ekman layer, as explained by Garrett et al. (1993), and no friction layer forms at the upslope side of the gravity current. The speed of descent of the downslope front of the friction layer slightly decreases with time due to the decrease of its density anomaly. Indeed, there is a no-slip boundary condition at the ocean floor, the water right above it cannot keep up with the downward speed of the front of the gravity current and gravity current water superposes surrounding water near the downward progressing front. This surrounding water becomes mixed into the friction layer near the front and dilutes the gravity current water near the downward progressing front, a feature common to all gravity currents. The descent of the centre of gravity of the vein is well fitted by a linear law (see Fig. 4). The rates of descent are given in Table 2. The angle of descent of the vein compares well to the theoretically predicted value of 1.2 · 10 −2 , based on a linear force balance Fig. 4 The lower and upper bounds along the slope of the gravity current (black curves) and the vein (red curves) are shown. That is, the gravity current evolves (in time) between the black lines and the vein between the red lines. The vein is defined to be the part of the gravity current more than 20 m above ground. The path of the centre of gravity of the total gravity current (green line) and the vein only (blue line) are also shown. All the results are from exp. G03; other experiments show the same qualitative behaviour

4.7. ON THE BASIC STRUCTURE OF OCEANIC GRAVITY CURRENTS 137<br />

Ocean Dynamics (2009) 59:551–563 557<br />

Fig. 3 Structure and vertical velocity field at 24 h in exp. G15.<br />

The temperature anomaly, with respect to the ambient water,<br />

is shown (given in Kelvin). Isolines of the w component of the<br />

velocity are given in black for positive values (upward) and<br />

red for negative values; the contour lines shown are ±10 −5 , 3 ·<br />

10 −5 , 10 −4 , 3 · 10 −4 , 10 −3 , 3 · 10 −3 ; the logarithmic scale was chosen<br />

due to the large intermittency of the vertical velocity. Vertical<br />

velocities outsi<strong>de</strong> the vein are due to wave-like motion of the<br />

interface<br />

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

shows only a small variation during the evolution of<br />

the gravity current. This shows that there is no substantial<br />

entrainment or d<strong>et</strong>rainment for the total gravity<br />

current.<br />

The <strong>de</strong>nsity anomaly and the along-slope velocity (y<br />

component) both <strong>de</strong>crease in the experiments. To leading<br />

or<strong>de</strong>r, the dynamics is supposed to be given by the<br />

geostrophic equilibrium as discussed in Subsection 2.1.<br />

To verify this, I calculate the along-slope velocity component<br />

in the vein divi<strong>de</strong>d by the geostrophic speed of<br />

its reduced gravity, given by: v G = (g ′ /f) tan α, for each<br />

<strong>de</strong>nsity class within the vein. When the average for all<br />

<strong>par</strong>cels in the vein which are at least 20 m above ground<br />

is analysed, the average along-slope velocity <strong>de</strong>creases<br />

to about 95% the geostrophic velocity after 2 days.<br />

This shows that the gravity current progresses along<br />

the slope at almost geostrophic speed, in agreement<br />

with previous findings of Price and Baringer (1994) and<br />

others.<br />

Friction leads to a downslope transport of gravity<br />

current water. The dynamics of the vein differs compl<strong>et</strong>ely<br />

from the dynamics of the total gravity current,<br />

and so does their rate of <strong>de</strong>scent and their spreading,<br />

as can be verified in Fig. 4. The spread and <strong>de</strong>scent of<br />

the gravity current, which inclu<strong>de</strong>s the friction layer, is<br />

much larger than that of the vein alone.<br />

Figure 4 shows that the <strong>de</strong>pth of the upper bound<br />

of the vein stays almost constant, a feature that is<br />

well documented (see Price and Baringer 1994 and<br />

references therein), while the downslope si<strong>de</strong>s of the<br />

vein and the friction layer <strong>de</strong>scend. At the upslope<br />

si<strong>de</strong>, the y component of the velocity vector is negative,<br />

as noted above, the fluid in the Ekman layer should,<br />

thus, move upslope. This leads to an arrested Ekman<br />

layer, as explained by Garr<strong>et</strong>t <strong>et</strong> al. (1993), and no<br />

friction layer forms at the upslope si<strong>de</strong> of the gravity<br />

current. The speed of <strong>de</strong>scent of the downslope front<br />

of the friction layer slightly <strong>de</strong>creases with time due<br />

to the <strong>de</strong>crease of its <strong>de</strong>nsity anomaly. In<strong>de</strong>ed, there<br />

is a no-slip boundary condition at the ocean floor, the<br />

water right above it cannot keep up with the downward<br />

speed of the front of the gravity current and gravity<br />

current water superposes surrounding water near the<br />

downward progressing front. This surrounding water<br />

becomes mixed into the friction layer near the front<br />

and dilutes the gravity current water near the downward<br />

progressing front, a feature common to all gravity<br />

currents.<br />

The <strong>de</strong>scent of the centre of gravity of the vein is<br />

well fitted by a linear law (see Fig. 4). The rates of<br />

<strong>de</strong>scent are given in Table 2. The angle of <strong>de</strong>scent of<br />

the vein com<strong>par</strong>es well to the theor<strong>et</strong>ically predicted<br />

value of 1.2 · 10 −2 , based on a linear force balance<br />

Fig. 4 The lower and upper bounds along the slope of the gravity<br />

current (black curves) and the vein (red curves) are shown. That<br />

is, the gravity current evolves (in time) b<strong>et</strong>ween the black lines<br />

and the vein b<strong>et</strong>ween the red lines. The vein is <strong>de</strong>fined to be the<br />

<strong>par</strong>t of the gravity current more than 20 m above ground. The<br />

path of the centre of gravity of the total gravity current (green<br />

line) and the vein only (blue line) are also shown. All the results<br />

are from exp. G03; other experiments show the same qualitative<br />

behaviour

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

Saved successfully!

Ooh no, something went wrong!