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

112 CHAPITRE 4. ETUDES DE PROCESSUS OCÉANOGRAPHIQUES APRIL 2008 W I R T H A N D B A R N I E R 811 FIG. 5. Correlation C y (500 m, U) for experiments (left) E04 and (right) E34. tel-00545911, version 1 - 13 Dec 2010 u w 2 F g fF f 2 F 2 2 f 2 c˜ Ff 2 2 Ff , 1 2 f , 10 11 where c˜ g f 2 /[(f 2 F 2 2 )] is the unperturbed vertical plume speed and is the inverse friction time. We thus see that when friction is negligible (and c˜ is finite) then the velocity vector is directed downward along the axis of rotation. By increasing the friction (while keeping c˜ constant), we see that the rotation vector starts to tilt in the eastward (to first order in /f ) and downward (to second order in /f ) directions. A slight eastward tilt can be seen in Figs. 3, 4, and 5 because the correlation patterns are slightly shifted to the left. The shifts show that /f O(10 1 ), but precise values cannot be obtained from our data. A downward shift is hard to detect because (i) it is only second order and (ii) it is countered by the correlation of isotropic turbulence (cf. to left part of Fig. 3). The higher correlations in the left panel of Fig. 3 for downward shifts are due to the increase of characteristic length scale with depth (see Fig. 9). The southward and eastward deviations have been observed and explained for the case of a single convective plume by Sheremet (2004) and Wirth and Barnier (2006). b. Mean values In this section we focus on the values of horizontal averages, because the problem is statistically homogeneous in both horizontal directions. The averages over a complete horizontal (periodic 8 km 8 km) slice of the ocean are denoted by . h . The typical size of a descending plume in the herein-presented experiments is of the order of 1 km; horizontal averages are as such only averages over a few plumes and the significance of the statistics is limited. To decrease the statistical uncertainty of the averages we also used time averaging, by averaging over the last 30 snapshots of the integration separated by 3 h, that is, the last 90 h of the experiments. We carefully checked that the quantities considered have reached a statistically stationary state before the beginning of the sampling process. These averages are denoted by . h,t , and their spatiotemporal evolution is what a perfect parameterization of the convective process should reproduce in an OGCM. 1) TEMPERATURE The most important variable to consider is either the temperature or buoyancy density [both are linearly related, see Eq. (6)]. These quantities do not reach a statistically stationary state because their mean value varies linearly in time. The dynamics are, however, independent of the mean value (linear equation of state), and the derivatives of temperature reach a statistically stationary state. In Fig. 1 the temporal evolution of the horizontally averaged temperature T h is presented. The qualitative behavior is an initial downwardpropagating front of cold water. Above the front the average temperature is almost constant, except for a boundary layer at the surface. When the front has reached the ocean floor, the horizontally averaged temperature at every depth increases linearly in time at the same rate, leading to a stationary (inverse) stratification. The results shown in Fig. 6 indicate that stratification depends on the strength of forcing but not on the direction of the rotation vector. Figure 6, as well as the isotemperature lines in Fig. 1, shows clearly that the temperature gradient is not constant. The vertical temperature gradient can be well approximated by a linear behavior away from the top and bottom boundary (see Fig. 6). With a dependence

4.5. MEAN CIRCULATION AND STRUCTURES OF TILTED OCEAN DEEP CONVECTION113 812 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 tel-00545911, version 1 - 13 Dec 2010 FIG. 6. Horizontally and temporally averaged temperature gradient z T h, t (see Table 2 for experiments). on the surface buoyancy flux (B 0 ), which is consistent with dimensional analysis, the second derivative of the mean temperature with respect to the vertical is thus constant, zz T h,t 2/3 H 7/3 /(g), where we estimated the dimensionless constant 10. In all of the results on temperature and its vertical gradients no significant difference between the nontilted and tilted case could be detected, except for the differences in the bottom Ekman layer. Also of great importance is the observation that in the lower 500 m the temperature gradient is positive, which means that there is a countergradient flux of heat. This behavior is often named “nonlocal transport,” and it is caused by the buoyancy transport of the convective plumes that extent from the surface to the bottom of the domain. Fitting an affine law to the gradient of the depthaveraged temperature allows the buoyancy flux to be written as Z g z T B 0z H , 12 where denotes the nonlocal part of the flux. Equation (12) is taken from the KPP parameterization (see Large et al. 1994), and the term on the right-hand side is a consequence of the stationarity of the dynamics. Dimensional analysis suggests Z 0 (B 0 H 4 ) 1/3 , where the dimensionless constant 0 1 0.1 best fits our data, leading to Z 40 m 2 s 1 for experiments E01 and E31, and a heat flux of only 250 W m 2 leads to an eddy diffusivity of about Z 5 m 2 s 1 for a convection depth of 1 km. In OGCM calculations a constant eddy diffusivity of Z 10 m 2 s 1 , independent of depth and buoyancy forcing is often employed, which lies in between the values obtained here. If we deduce from Fig. 6 that the temperature gradient is zero at about 500 m from the ground, we can obtain a nonlocal (countergradient) heat flux of about 1/7th of the heat flux at the surface. Dimensional analysis gives 0 (B 0 /H 2 ) 2/3 , with 0 1/7. We would like to emphasize that the total vertical heat flux is always positive (negative temperature perturbations being transported downward), and decreases linearly from the maximum value at the surface to zero at the floor. In our calculations we cannot have an upward heat flux, which is often observed when convection into a stratified medium is considered, because in our experiment the lower boundary is the completely insulating ocean floor. 2) VELOCITY Due to the top and bottom boundaries and the incompressibility of a Boussinesq fluid, the horizontally

4.5. MEAN CIRCULATION AND STRUCTURES OF TILTED OCEAN DEEP CONVECTION113<br />

812 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<br />

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

FIG. 6. Horizontally and temporally averaged temperature gradient z T h, t (see Table 2 for<br />

experiments).<br />

on the surface buoyancy flux (B 0 ), which is consistent<br />

with dimensional analysis, the second <strong>de</strong>rivative of the<br />

mean temperature with respect to the vertical is thus<br />

constant, zz T h,t 2/3 H 7/3 /(g), where we estimated<br />

the dimensionless constant 10.<br />

In all of the results on temperature and its vertical<br />

gradients no significant difference b<strong>et</strong>ween the nontilted<br />

and tilted case could be d<strong>et</strong>ected, except for the<br />

differences in the bottom Ekman layer. Also of great<br />

importance is the observation that in the lower 500 m<br />

the temperature gradient is positive, which means that<br />

there is a countergradient flux of heat. This behavior is<br />

often named “nonlocal transport,” and it is caused by<br />

the buoyancy transport of the convective plumes that<br />

extent from the surface to the bottom of the domain.<br />

Fitting an affine law to the gradient of the <strong>de</strong>pthaveraged<br />

temperature allows the buoyancy flux to be<br />

written as<br />

Z g z T B 0z<br />

H ,<br />

12<br />

where <strong>de</strong>notes the nonlocal <strong>par</strong>t of the flux. Equation<br />

(12) is taken from the KPP <strong>par</strong>am<strong>et</strong>erization (see Large<br />

<strong>et</strong> al. 1994), and the term on the right-hand si<strong>de</strong> is a<br />

consequence of the stationarity of the dynamics. Dimensional<br />

analysis suggests Z 0 (B 0 H 4 ) 1/3 , where<br />

the dimensionless constant 0 1 0.1 best fits our<br />

data, leading to Z 40 m 2 s 1 for experiments E01<br />

and E31, and a heat flux of only 250 W m 2 leads to an<br />

eddy diffusivity of about Z 5 m 2 s 1 for a convection<br />

<strong>de</strong>pth of 1 km. In OGCM calculations a constant eddy<br />

diffusivity of Z 10 m 2 s 1 , in<strong>de</strong>pen<strong>de</strong>nt of <strong>de</strong>pth and<br />

buoyancy forcing is often employed, which lies in b<strong>et</strong>ween<br />

the values obtained here.<br />

If we <strong>de</strong>duce from Fig. 6 that the temperature gradient<br />

is zero at about 500 m from the ground, we can<br />

obtain a nonlocal (countergradient) heat flux of about<br />

1/7th of the heat flux at the surface. Dimensional analysis<br />

gives 0 (B 0 /H 2 ) 2/3 , with 0 1/7.<br />

We would like to emphasize that the total vertical<br />

heat flux is always positive (negative temperature perturbations<br />

being transported downward), and <strong>de</strong>creases<br />

linearly from the maximum value at the surface to zero<br />

at the floor. In our calculations we cannot have an upward<br />

heat flux, which is often observed when convection<br />

into a stratified medium is consi<strong>de</strong>red, because in<br />

our experiment the lower boundary is the compl<strong>et</strong>ely<br />

insulating ocean floor.<br />

2) VELOCITY<br />

Due to the top and bottom boundaries and the incompressibility<br />

of a Boussinesq fluid, the horizontally

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

Saved successfully!

Ooh no, something went wrong!