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

110 CHAPITRE 4. ETUDES DE PROCESSUS OCÉANOGRAPHIQUES APRIL 2008 W I R T H A N D B A R N I E R 809 FIG. 2. Isosurfaces of vertical velocity w 0.05 m s 1 ( red, blue) looking upward from the ocean floor (x to the right, y downward) at the end of the experiments at t 168 h: (left) E04 and (right) E34. The elongation of structures in the y direction is conspicuous in experiment E34. tel-00545911, version 1 - 13 Dec 2010 front descends into the unstratified ocean until it reaches the ocean floor after about 18–30 h, depending on the strength of the forcing and rotation. The subject of this paper is the subsequent statistically stationary stage of convective dynamics in the entire water column. During this stage descending plumes surrounded by rising water mix the water column as can be seen in Fig. 2. A shallow boundary layer develops at the ocean surface. As in the case of a single plume, the influence of the (no slip) bottom slows down the front propagation 500 m (the typical horizontal plume scale) above the bottom, emphasizing once more the importance of the Ekman layer dynamics at the ocean floor. Figure 2 exposes that the differences in the horizontally averaged temperature structure between the tilted and nontilted dynamics are small. a. The signature of the Taylor–Proudman–Poincaré theorem The analysis presented in this section is based on the discussion of the derivation, generality, and consequences of the Taylor–Proudman–Poincaré theorem by Colin de Verdière (2002). When friction and nonlinearity are negligible the dynamics of a stratified Boussinesq fluid, subject to rotation, is governed by the TPP theorem: 2 y y z z u y w g 0 x . 0 7 The TPP theorem states that in the direction of the axis of rotation the velocity component aligned with gravity (w) does not change. While the constraint in the direction of the axis of rotation on the other two components of the velocity vector is influenced by changes in the density structure, the constraint of the TPP theorem on the horizontal components of the velocity vector is in fact the classical thermal wind relation in the case of a vertical rotation vector, and we will call it the generalized thermal wind relation for arbitrary directions of the rotation vector. There are numerous discussions of the TPP theorem when applied to both large- and mesoscale ocean dynamics, that is, when the Rossby number associated with the dynamics is (very) small. In our case, however, the Rossby number Ro w rms /(fL) is an order of unity or larger, and the dynamics are clearly in a three-dimensional regime. The TPP theorem nevertheless leaves an imprint in the turbulent dynamics. To investigate this imprint we calculate the correlation of a component of the velocity vector in two parallel planes that are separated in the y direction by a distance y. The southward plane, which spans the entire (periodic) domain in the x direction and 1 km in the z direction (from a depth of 2250 to 1250 m), is kept fixed, while the northward plane, which is of equal size, is shifted in the x and z direction by the amount (x, z). The correlation between a velocity component in the two planes is then calculated. In mathematical terms (applied to the U component), we calculate the following correlations:

4.5. MEAN CIRCULATION AND STRUCTURES OF TILTED OCEAN DEEP CONVECTION111 810 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 38 FIG. 3. Correlation C y (500 m, W) for experiments (left) E04 and (right) E34. C˜ yy, U, y 0 , t Ux, y 0 , z, tUx x, y 0 y, z z, t dx U P dzP 2 x, y 0 , z, t dx dz P U 2 x x, y 0 y, z z, t dx dzC y y, U C˜ yy, U, y 0 , t y0 ,t, 8 tel-00545911, version 1 - 13 Dec 2010 where . y0 ,t is an averaging over values y 0 separated by 1 km and 30 consecutive time instances separated by 3 h during the last 90 h of the experiment. The results of this analysis are presented in Figs. 3, 4, and 5. The predictions of Eq. (7) are nicely confirmed as follows: (i) The maximum correlation of the w component is displaced about 500 m and 1 km upward for a y distance y of the two planes of 500 m and 1 km, respectively, which demonstrates the high correlation along the axis of rotation; (ii) the same is approximately true for the u component, but the correlations are only about half, due to the decorrelation by the generalized thermal wind. The cases with stronger forcing look qualitatively the same. Quantitative changes are due to the increased horizontal scales in the cases with stronger forcing [rotation is less efficient in stopping the horizontal growth of plume structures; see Wirth and Barnier (2005)]. If we suppose that the convective motion is composed of descending plumes surrounded by upwardmoving fluid and an isotropic turbulent component, we can build an analytical model of the convective process, based on the balance of rotation (f and F), turbulent friction (), and reduced gravity (g) acting on a convective parcel of fluid moving with the velocity (u, , w). More precisely, f F f 0 F 0 u w 0 0 g . Solving the linear system leads to 9 FIG. 4. Correlation C y (1 km, W) for experiments (left) E04 and (right) E34.

110 CHAPITRE 4. ETUDES DE PROCESSUS OCÉANOGRAPHIQUES<br />

APRIL 2008 W I R T H A N D B A R N I E R 809<br />

FIG. 2. Isosurfaces of vertical velocity w 0.05 m s 1 ( red, blue) looking upward from the ocean floor (x<br />

to the right, y downward) at the end of the experiments at t 168 h: (left) E04 and (right) E34. The elongation<br />

of structures in the y direction is conspicuous in experiment E34.<br />

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

front <strong>de</strong>scends into the unstratified ocean until it<br />

reaches the ocean floor after about 18–30 h, <strong>de</strong>pending<br />

on the strength of the forcing and rotation.<br />

The subject of this paper is the subsequent statistically<br />

stationary stage of convective dynamics in the entire<br />

water column. During this stage <strong>de</strong>scending plumes<br />

surroun<strong>de</strong>d by rising water mix the water column as can<br />

be seen in Fig. 2. A shallow boundary layer <strong>de</strong>velops at<br />

the ocean surface. As in the case of a single plume, the<br />

influence of the (no slip) bottom slows down the front<br />

propagation 500 m (the typical horizontal plume scale)<br />

above the bottom, emphasizing once more the importance<br />

of the Ekman layer dynamics at the ocean floor.<br />

Figure 2 exposes that the differences in the horizontally<br />

averaged temperature structure b<strong>et</strong>ween the tilted and<br />

nontilted dynamics are small.<br />

a. The signature of the Taylor–Proudman–Poincaré<br />

theorem<br />

The analysis presented in this section is based on the<br />

discussion of the <strong>de</strong>rivation, generality, and consequences<br />

of the Taylor–Proudman–Poincaré theorem by<br />

Colin <strong>de</strong> Verdière (2002).<br />

When friction and nonlinearity are negligible the dynamics<br />

of a stratified Boussinesq fluid, subject to rotation,<br />

is governed by the TPP theorem:<br />

2 y y z z u y <br />

w g 0 x<br />

. 0 7 The TPP theorem states that in the direction of the<br />

axis of rotation the velocity component aligned with<br />

gravity (w) does not change. While the constraint in the<br />

direction of the axis of rotation on the other two components<br />

of the velocity vector is influenced by changes<br />

in the <strong>de</strong>nsity structure, the constraint of the TPP theorem<br />

on the horizontal components of the velocity vector<br />

is in fact the classical thermal wind relation in the<br />

case of a vertical rotation vector, and we will call it the<br />

generalized thermal wind relation for arbitrary directions<br />

of the rotation vector.<br />

There are numerous discussions of the TPP theorem<br />

when applied to both large- and mesoscale ocean dynamics,<br />

that is, when the Rossby number associated<br />

with the dynamics is (very) small. In our case, however,<br />

the Rossby number Ro w rms /(fL) is an or<strong>de</strong>r of unity<br />

or larger, and the dynamics are clearly in a three-dimensional<br />

regime. The TPP theorem nevertheless<br />

leaves an imprint in the turbulent dynamics. To investigate<br />

this imprint we calculate the correlation of a component<br />

of the velocity vector in two <strong>par</strong>allel planes that<br />

are se<strong>par</strong>ated in the y direction by a distance y. The<br />

southward plane, which spans the entire (periodic) domain<br />

in the x direction and 1 km in the z direction<br />

(from a <strong>de</strong>pth of 2250 to 1250 m), is kept fixed, while<br />

the northward plane, which is of equal size, is shifted in<br />

the x and z direction by the amount (x, z). The correlation<br />

b<strong>et</strong>ween a velocity component in the two<br />

planes is then calculated. In mathematical terms (applied<br />

to the U component), we calculate the following<br />

correlations:

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

Saved successfully!

Ooh no, something went wrong!