Научная статья на тему 'Vertical Fluxes Induced by Weak-Nonlinear Internal Waves in a Baroclinic Flow'

Vertical Fluxes Induced by Weak-Nonlinear Internal Waves in a Baroclinic Flow Текст научной статьи по специальности «Физика»

CC BY
91
12
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Physical Oceanography
WOS
Область наук
Ключевые слова
internal waves / turbulence / Stokes drift

Аннотация научной статьи по физике, автор научной работы — A. A. Slepyshev

In a Bussinesq approximation free internal waves are considered at the account of turbulent viscosity and diffusion in a vertically-non-uniform flow. In linear approximation the dispersive relation and decrement of wave attenuation are found. The equation for amplitude of the vertical velocity, containing small parameter at the higher derivative, proportional to turbulent viscosity, solved by Ljusternik-Vishik asymptotic method. Boundary-layer solution on a vicinity of a bottom and a free surface are found. The non-viscous boundary value problem of the second order is solved numerically according to an Adams's implicit scheme of the third order accuracy. The wave number at the fixed frequency of a wave is found by a shooting method. Stokes drift velocity and vertical wave flux of salt are determined in the second order on wave amplitude. Shift of phases between fluctuations of salinity and vertical velocity with regard to turbulent viscosity and diffusion differs from /2π, therefore the vertical wave flux of salt differs from zero. The dispersive parity, decrement of wave attenuation and wave fluxes are calculated for internal waves observed during the 3rd stage of the 44 cruise of R/V "Mikhail Lomonosov" to the northwest shelf of the Black sea. Critical layers for a current profile were absent at the test site (where the measurements were carried out), i.e. phase rate of internal waves exceeded the current velocity. It is shown that out of a layer with the maximum gradient of salinity, i.e. out of surface layer, the wave flux of salt is comparable in absolute value with the turbulent one. In a surface layer turbulent flux of salt exceeds the wave one. The consideration of current results in slight decrease of the wave flux. Horizontal component of Stokes drift velocity, which is transversal to the wave propagation direction, differs from zero and is one order less than longitudinal one when taking into account the current.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Vertical Fluxes Induced by Weak-Nonlinear Internal Waves in a Baroclinic Flow»

Vertical Fluxes Induced by Weak-Nonlinear Internal Waves in a Baroclinic Flow

A.A. Slepyshev

Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol,

Russian Federation

In a Bussinesq approximation free internal waves are considered at the account of turbulent viscosity and diffusion in a vertically-non-uniform flow. In linear approximation the dispersive relation and decrement of wave attenuation are found. The equation for amplitude of the vertical velocity, containing small parameter at the higher derivative, proportional to turbulent viscosity, solved by Ljusternik-Vishik asymptotic method. Boundary-layer solution on a vicinity of a bottom and a free surface are found. The non-viscous boundary value problem of the second order is solved numerically according to an Adams's implicit scheme of the third order accuracy. The wave number at the fixed frequency of a wave is found by a shooting method. Stokes drift velocity and vertical wave flux of salt are determined in the second order on wave amplitude. Shift of phases between fluctuations of salinity and vertical velocity with regard to turbulent viscosity and diffusion differs from n / 2, therefore the vertical wave flux of salt differs from zero. The dispersive parity, decrement of wave attenuation and wave fluxes are calculated for internal waves observed during the 3rd stage of the 44 cruise of R/V "Mikhail Lomonosov" to the northwest shelf of the Black sea. Critical layers for a current profile were absent at the test site (where the measurements were carried out), i.e. phase rate of internal waves exceeded the current velocity. It is shown that out of a layer with the maximum gradient of salinity, i.e. out of surface layer, the wave flux of salt is comparable in absolute value with the turbulent one. In a surface layer turbulent flux of salt exceeds the wave one. The consideration of current results in slight decrease of the wave flux. Horizontal component of Stokes drift velocity, which is transversal to the wave propagation direction, differs from zero and is one order less than longitudinal one when taking into account the current.

Keywords: internal waves, turbulence, Stokes drift.

DOI: 10.22449/1573-160X-2015-1-59-72

© 2015, A.A. Slepyshev © 2015, Physical Oceanography

Introduction. Internal waves are presented everywhere in the ocean due to effect of sources which generate them. Internal waves can exist in stratified medium when water density increases with the depth. Below the upper mixed layer such situation is typical for the World Ocean. Atmospheric pressure perturbations, wind stresses on the sea surface, interaction of tides and currents with bottom topography inhomogeneities [1], eddy currents may be referred to internal wave energy sources.

Actuality of the problems is related to the fact that internal waves can contribute to vertical transfer in the ocean. Usually, vertical transfer in marine environment is connected with small-scaled turbulence that has an intermittent character, i. e. the turbulence is presented in the form of "patches", generated by hydrodynamic instability of currents and breaking of internal waves. Vertical transfer plays an important role in admixture transport, oxygen diffusion to the deep layers of the sea and hydrogen sulphide diffusion from deep layers of the Black Sea.

The internal waves with regard for eddy viscosity and diffusion were considered in a number of works [1 - 3], where a decrement of wave attenuation on

PHYSICAL OCEANOGRAPHY NO. 1 (2015) 59

turbulence was determined. Non-linear effects at wave propagation with no regard for eddy viscosity and diffusion were considered in [4, 5]. In these works a mean wave-induced current and non-oscillatory correction to the average density were determined.

Internal waves with regard for nonlinearity and viscosity were considered in [6], but vertical wave fluxes of heat and salt were not taken into account.

Vertical transfer in stratified seawater column is associated with breaking of small-scale internal wave [7], and turbulent energy dissipation rate and a coefficient of vertical eddy diffusivity are estimated. It is shown that in the area of continental slope at the edge of the Black Sea shelf the intensification of vertical transfer, related with internal wave amplitude increase at their propagation to the shallow water area [8], takes place.

Vertical flows, determined by weak-nonlinear internal waves, were considered in [9]. Such flows exist due to vertical velocity and temperature (salinity) oscillation phase shift, which differs from n/2 when diffusion and eddy viscosity are taken into account.

In the given work the vertical flows caused by internal waves on the baroclinic flow are considered. It is of interesting to compare wave fluxes with the corresponding turbulent fluxes both in the presence of current and in the absence of it. Moreover, the Stokes drift speed in both these cases is found.

Formulation of the problem. Free internal waves on the baroclinic flow with account for eddy viscosity and diffusion in the Boussinesq approximation are considered. Internal wave amplitude distribution, dispersive relation and wave attenuation decrement are calculated in linear approximation. The Stokes drift speed and wave fluxes of salt are found in the second order of the wave amplitude.

We introduce dimensionless variables by the following formulas (dimensional physical quantities are denoted by wavy lines above the symbols):

t ~ k

t = —, Ik = —, cc = ú)„Ú) , U~1 = u1Hc, u2 = u2Hc u~3 = u3Ha„,

C H

~ H n H n ~

P = p0H VP, p = p0V-, p0 = p0V-0, xi = Hxi, Kt = ,

g g

Mi = Mm (i = 1,2,3), K = K2, Ml = M2, M =^2,

where /u1, ju3 - are the characteristic values of horizontal and vertical eddy viscosity; g - is an acceleration of gravity; x1, x2, x3 - are two horizontal and one vertical coordinates, vertical axis x3 is directed vertically upwards; p and P - are the wave perturbations of the density and pressure; p0 - is an unperturbed average water density; p0 - is a depth-averaged density; u1, u2, u3 - are two horizontal and one vertical components of velocity wave disturbances; K1, K3, M1, M3 - are horizontal and vertical coefficients of eddy viscosity and

diffusion, respectively; H - is a sea depth; v - is a characteristic wave

60

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

frequency. The coefficients of vertical eddy diffusivity and viscosity, and two components of the mean current velocity U0, V0 are assumed to be dependent on

the vertical coordinate. The coefficients of horizontal eddy diffusivity and viscosity are assumed to be constant. The system of hydrodynamic equations for the wave disturbances in the Boussinesq approximation has the following form:

Cki, Cki,

_1__L i i _1

dt

u

' * U0

ox.

du, ^ du,

____L T/ __

dx_ 0 cX2

dUn

u

dP dx1

s2 Kii+d

dx_ cX2

d

^ dx3 ^ 3 dx3 ,

du1

(la)

du-, du,, - u

dt

' dx.

-U

du2 dx.

du2 dVr,

- V+ u

dx.

dx.

dx.

dP 2 ( d2u^ d2u,, - + s K.

dxx

+ s2p2 Mk3 ——2

dx, I dx,

d

(_b)

du3 ~dt

du, ^ du, ^ du, u. —3 + U0 —3 + V0 —-

d

dx1

s2 p2 -L

dx.

dx,.

K du3

dP dx.

■+s2 Ki

i ~i2 ~,2 \ du, du,

--H--3

v dx_ dx\ j

dx.

3 j

s2p2 àK— du3

dx3 dx3

P,

(1c)

dp dp dp dp 2-Kt

+ u —— + Un—+ = s M,

dt ' dx. 0 dx 0 dxn ' 1 2

fd2p d2p ^ y Sx_2 5x2 j

2 o2 d (dp "\ dp,

+ s p --- u 3 ——

dx3 ^ dx3 j dx3

du,

dx.

= 0.

, (1d)

(_e)

Here s2 = M - is a small parameter proportional to the value of the horizontal

H 2a,

eddy viscosity coefficient p2 = M-, while << s2.

Mi

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Boundary conditions on the surface are the solid top condition and the absence of tangential stresses [2]:

u-(0) = 0, p2K- du_ + Kl du- = 0

ox3 ox1

p2 K3 ^u2+K_ ^u- = 0

ôx dx.

(2)

x3 =0

x3 =0 3 2

Boundary conditions on the bottom are the adhesion conditions:

u_( -1) = 0, u2(-1) = 0, u3(-1) = 0. (3)

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

61

Linear approximation. We seek the solutions of linear approximation in the following form:

U = u10(x3)Ae'e + c.c., u0 = u20(x3)Ae'e + c.c., ul = u30(x3)Ae'e + c.c.,

p = Po(X3)Aew + c.c., p = pio(X3)Aew + c.c., (4)

where c.c. - are complex conjugate terms; A - amplitude factor; d - wave phase;

= k; = —v (k - wave number, V - frequency). It is assumed that the

dx1 dt

wave propagates along the x1 axis.

Substituting (4) into the system of equations (1), we obtain the relation of amplitude functions u10, Pw, u30and equation for u20(x3), p10(x3):

i du

u10 =--—, Q = v- kU0, (5a)

k dx3

P10 = (s2p2K3 ^ + s2p2 dK3— + i'Q-s2K1k2)^^ + LdU.u30, (5b) dx^ d^3 d^3 k d^3 k d^3

(s202M3 ddxY + s2p2 Mdr s2Mk2 + i^)A0 = u30 dp, (5c)

CIJ\'3 dJ\'3 CUx'3 dJ\'3

s2p2K3 ^ + s2P2 ^ ^ + (i^s2Klk2)u20 =u30 dpi. (5d)

CIJ\'3 dJ\'3 dJ\'3 dJ\'3

(5d) equation should be supplemented by boundary conditions arising from (2), (3):

du.

20

dx,

3 x3=0

The function u30 satisfies the equation

-N2u30 = (s2^2M3 ^ + s1p1 -e1Mlk2 + iQ) x

= 0; U20 X=_j = 0. (6)

x <j iQu30 —— I dx3

dx> 3 dx3 dX3

K3 d2- + s2?2 - s2 K1k2 + iQ)^- du30 + u3{

dx 3 dx3 dx3 k dx3 k dx3

(7)

-s2K1k2u30 + s2$2K3 ^ + 2s2^2 dK3 du30

dX3 dX3 A3

on the surface at x3 = 0

From the boundary conditions (2), (3) we obtain the conditions for u30 :

d 2u3

dx32

62 PHYSICAL OCEANOGRAPHY NO. 1 (2015)

u30 = 0, —f- = 0, (8a)

on the bottom at x3 = -1

u3o = ^ = 0. (8b)

dx3

Equation (7) has a small parameter at the higher derivative ~(s3)4 . We solve this equation, following [10, 11], by the Lusternik - Vishik asymptotical method, expanding u30, < in series by s? small parameter:

U30 (X3) = X (spy ui (X3, s) + s?Z (s//) v + (s/?)2X (s//)'vi, (9a)

< = <01(s) + s/<11(s) + (s/)2<21(s) + ..., (9b)

where v1 ((1 + x3) / sp) - boundary layer solution in the neighborhood of the bottom, v°(x3/sp) - in the neighborhood of free surface. Boundary layer

corrections are the functions, rapidly decreasing with the distance from the boundary, which provide the execution of the boundary conditions. In the general case the functions u3i(x3,s) depend on the parameter s , because it is contained

in the equations for these functions.

Substituting the expansions (9) in (7), we obtain the boundary value problem for the function u3°0 in the first order on s3 parameter:

(s2Mik2 -/Qoi)J/CV0o

I dx,

(/Q01 -^Kt)+~dJT"0o -^KU k dx k dx

3 3

:#2Ml,(io)

where N2 = - - square of the Brent-Vaissala frequency; Qoi = co1 - kUo -dx3

wave frequency with Doppler shift.

Boundary conditions for u°o are the following:

uox=o = o, = o. (11)

The equation (1o) has a small parameters and following the method, described in [1o], we seek the solution and frequency co1 in a form of asymptotic series by this parameter:

u3o(x3,s) = wo(x3) + ew1( x3) + s2w2(x3) +..., (12a)

co1 = co +sc1 + s2c2 +.... (12b)

PHYSICAL OCEANOGRAPHY NO. 1 (2o15)

63

In a zero order on parameter s function w0 satisfy boundary value problem

(N2 -«2).

d2w0 2(N2 -n2) k d2Ut ' + k

dx,

a

w0 = 0.

0 1x3 =0

2 W0 +

a0 dx3

w = 0.

2 "0

w0 = 0,

0 lx3 =-1

(13a)

(13b)

Where Q0 = <0 - kU0. In the absence of singularity Q0 = 0 the boundary value problem (13) has a numerable set of eigenfunctions - a set of modes and for any value of frequency <0 < max(N) corresponds a certain wave number k for given mode.

The next term in the expansion (12a) is determined from the following equation

d2w ,2(N2 -a2) k d2U0 m

-2- + k --w1 +--20 w1 = —

dx3 a2 a0 ^^0

Boundary conditions for w1 function

2kw0 - 2 ^.2

d w kwn d U,

2

= f ( X3 )(14)

W-

1L =0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

= 0,

W-

1L =-1

= 0.

The solvability condition of the boundary value problem (14), (15)

u

J f1W0dx3 = 0.

(15)

(16)

For < ^ 0 this condition is not generally fulfilled and the boundary value problem (14), (15) has no solutions.

The next approximation w2 in the parameter s satisfies the following equation

a2

dxl

w2 + N 2k2 w2 +ka0 d-u20 w2 = (®2 +iMk2 )|k 2a0w0 -a0 -kd-U0 w„

dx, V 'I dx

-a

®2ik 2W0 \iKt2 dx0+Kk 4W0

= o.

(17)

Boundary conditions for w2 have the following form

w.

2 lx3 =0

= 0,

w.

2 Ix, =-1

= 0.

Solvability condition of the boundary value problem (17), (18)

u

J Ow0 dx3 = 0 .

(18)

(19)

64

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

Hence we find the expression for®2 :

a)2 = -i

2

0 f Mk\r2 ^ , 2 d2w0 ,4 ^

/[m^ - K'k *t+kKw ,

a

dx,

N2 2 k d U0 ^

— 2k2 + —-^

a3 a dx2

w0 dx3

3 J

Let us find a boundary layer solution v0 in the expansion (9a) to satisfy the boundary conditions (8a), (8b) on the surface. Substituting the expansion (9a) in the equation (7), and expanding K3 ,M3 and U0, V0 in Taylor series in the

neighborhood x3 = 0 we obtain the equation for v0(n), in the zero order on parameter sf :

K3 (0)M3 (0) d6v0 + /Qo(0)(K3(0) + M3(0)) ^ = Q0(0) d2v0, n = ^ (21)

dn dn dn sf

The solution of the equation (21), attenuating with the distance from the surface, has the following form:

V00 = C0 exp(^°n) + F,0 exp(^2n), (22)

where

,0 _ d2Wo(0)_1 _

dx23 ®2 - (40)2'

Cu = w- T0V»y_^_ Fu = _ c0 (93)

Here A,0 , 4 are calculated according to the formulas

40 = tofa_,), 4 = (WLf0_0. (24)

" 1 2M:,(0) I V ' 2 I 2^3(0). V ' K '

'J

The equation for the boundary layer solution v0((1 + x3)/s2) has the following form:

dV+ffl0( _i) Mi±KL dV dV=0, (25)

dtf M3K: dtf M3K: dtf

where n0 = ~-). The solution of the given equation, attenuating with the

sj3

distance from the bottom and satisfying the boundary conditions (8b), is determined by formula:

V0 = D° exp(_4n) + G exp(_4ni), (26a)

where

D0 = ^ -_-, G0 =_ D0. (26b)

caj^: /to

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

65

The values \ , 12 are calculated according to the formulas (24), only the eddy

transfer coefficients and wave frequency with the Doppler shift are calculated at the lower boundary.

The equation for salinity wave perturbations 5 has the form

ds

ds

ds

ds

dSn

d

f

— + (u1 + U0)— + (u2 + V0)-+ u3-+ u3—0 = s—

dt

dx,

M,

Sx,

ds

1 dx,

dx.

Sx

dx,

dx,

ds

\

M1 — . dx1 y

+ (sfifé

2 y

dx,

M

ds

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

dx.

(27)

3 y

where S0(x3) - is an average salinity profile. The solutions of linear approximation will be sought in the following form:

= j10(x3)Ae'e + c.c., (28)

where 510 satisfies the equation

dSn

d _ , ds,,

-1^10 + = s2 (-k2M^o +P2^- (M3 . (29)

dx,

d^3 (^3

Equation (29) has a small parameter at the higher derivative, and we are to seek a solution at the form similar to (9a):

s10 = s/o + (s/)2 sjj +

x + 1

(s/3)2 sfo +s/32 s^^/S +(s/)2 2 . (30)

s3

s3

From (29) it follows that where

j 2

s10 = s0 + S s2 + ... :

iw0 dS0

Q0

s2 =-(®2i + k 2M1 )

w dSn dSn

££2 <^3 ££0 d^3

(31) (32a) (32b)

0

here a>2i = — . Boundary layer solutions s n 0 and s n0 in the neighborhood of the

bottom and free surface satisfy the following equations, respectively

d (M (0)d )/„„ = dS0(0)-dn

(0)sn0 +d (M3(0)d)s°n0 = «S^1 v0(n)-.

dn

d

d

i£0(-1)sn0 + (M3(-1)—yn0 =

dS,(-1) ,

v0(n).

Boundary conditions have the form

s10(0) = s10(-1) = 0.

(33 a)

(33b)

(34)

i=0

i=0

66

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

Solution of the equation (33a) in the neighborhood of the sea surface is the following:

•C (n) = C0S exp(^0n) + Q0 exp(^0n) + nP0 exp(A), (35)

where

P = C dSM, (36a)

S 2A!0M3 (0) dx3

Q0 =_C_iS(0>, (36b)

Q 0(0)[1 - M3(0)/ K3(0)] dx3

C0 =-Q0. (36c)

Boundary layer solution of the equation (33b) in the neighborhood of the bottom has the following form:

S n 0 (ni ) = C12 S

exp(-A, ni) + Qsi exp(-A2n,) + n, PSi exp(-A, n, ), (37)

where

Ps,=__dS0£=i), (38a)

Si 2^M3( -i) dx3

Qs, =_D_dSai-i> . (38b)

-Sl Q,(-1)[1 -M,(-1)/Kj(-1)] dx,

From the boundary condition (34) it follows

CL = -Qsi. (39)

Nonlinear effects. The velocity of the Stokes drift of the fluid particles is determined by the formula [i2]:

u s = J udrVu, (40)

0

where u - is a field of Eulerian velocities, a horizontal bar means the averaging over the wave period. Horizontal velocity component of the Stokes drift (directed along the wave vector) in the second order of amplitude, has the following form:

„ = AiAi i 200 du30 du30 + «30 d U30 + «30 d U30\ (4i)

uis =~1 ( 1 + + * A 2 ) ' ^ ^

k oo dx3 dx3 6) dx3 6) dx3

Horizontal velocity component of the Stokes drift (transversal to the direction of wave propagation) is calculated according to the formula

u2s = d(U30U30) + c.c., A, = Ae6 . (42)

s o dx3

At the presence of the mean current, whose velocity component V0 (transversal

to the direction of wave propagation) depends on vertical coordinate, the value u2s differs from zero and in the non-viscous case is calculated as follows:

PHYSICAL OCEANOGRAPHY NO. i (20,5) 67

2 A A* d ,

U2s =--("

wL dVf)

dx^ ^^ q dx^

(43)

Vertical wave flux of salt u3S taking into account the expansions (9a), (3Q), (31) determined by formula:

U3s /| Ai2| = Wq(s2 s 2 +SJ3S n Q + (sJ3)2 O* + s Q)(sJ3v10 + (sj3)2 vQ) + s2 W2 sQ + c.c. (44)

Calculation results. Let us calculate salt wave fluxes for the internal waves, which had been observed southwestward of Yevpatoria during the experiment, performed at the third stage of the 44th R/V Mikhail Lomonosov cruise.

In Fig. 1 we have plotted four realizations of the temperature contour elevations calculated according to the data of the GRAD instruments, that is to say, the gradient-distributed temperature gauges [13]. The first device was situated in 5 - 15 m layer, the second - in 15 - 25 m layer, the third - in 25 - 35 m layer and the fourth - in 35 - 6Q m layer. It is obvious that strong oscillations with 15 m period in 25 - 6Q m layer are in against phase with the oscillations in 25 - 6Q m layer, which indicates the oscillations of the second mode. The maximum amplitude by swells made up Q.5 m. In Fig. 2, a, b the vertical profiles of the Brent-Vaissala frequency and two components of the current velocity are represented. The boundary value problem (13) for internal waves is solved numerically by the implicit Adams scheme of the third order accuracy. Wave number is found by the shooting method from the necessity of satisfying the boundary conditions (13b). The eigenfunction of 15 min internal waves of the second mode is shown in Fig. 2, b.

Fig. 1. Time variation of the vertical displacements temperature contours

This implies

max^3

A =

1 2 max |wQ / Qq i

68

PHYSICAL OCEANOGRAPHY NO. 1 (2Q15)

Thus, the amplitude of vertical displacements is proportional to w0. According to the experiment data (Fig. 1; 2, b), the maximum of w0 function correspond to the maximum vertical displacements, i.e. in the experiment the second mode was observed. We calculate a vertical eddy diffusion coefficient M3 by the empirical formula, which is valid in the area of continental slope at the Black sea north-west shelf [7]: M3 = 8.4•10-4Ni-1 m2/s, Nc corresponds to Brent-Vaissala frequency in cph. Wavelength of the second mode 15-min internal waves is 195.6 m, typical value of horizontal eddy diffusion coefficient M1 is 1 m2/s. The dispersion curves of the first two modes in the presence and absence of current are shown in Fig. 3.

Fig. 2. Vertical profile of Brent-Vaissala frequency - a, an eigenfunction of 15 min internal waves -b, vertical profiles of current velocity components jj (■■■), y (-)-c

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

69

л о.

к, rad/m

Fig. 3. Dispersion curves: the first mode without current (-) and with current (---); the second

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

mode without current (.........) and with current (-----)

Рис. 4. Stokes drift velocity in the presence of current (-) and in the absence of current (.........)

Boundary value problem (17), (18) on the determination of w2 function is solved numerically by the implicit Adams scheme of the third accuracy order at K3 = 2M3, Kj = 2Mj. We find the only solution, orthogonal wQ and wave

attenuation decrement 8(0 = (i . In case of 15 min internal waves of the second

mode 8( = -1.46 -1Q-3 rad/s. Vertical profiles of Stokes drift velocity horizontal component in the presence and absence (dashed curve) of current are shown in Fig. 4. In the presence of current the Stokes drift velocity in modulus is lesser. Stokes drift velocity component, transversal to the direction of wave propagation, differs from zero only with regard to the current (Fig. 5).

Fig. 5. Stokes drift velocity component, transversal to the direction of wave propagation Fig. 6. Vertical salinity profile

70

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

Fig. 7. Profiles of wave and turbulent salt fluxes: wave flux in the presence of current (-) and in

the absence of current (.........); turbulent salt flux (---)

Wave salt flux is calculated at the vertical salinity profile, shown in Fig. 6.

dS

Turbulent salt flux is calculated according to the formula sf = -M3 —0. Vertical

f 3 dz

wave salt flux u3S for 15 min internal waves of the second mode and a turbulent

flux are represented in Fig. 7. In the presence of current the wave salt flux is somewhat lesser than in its absence. In the absence of current a wave flux (dashed line) is comparable in modulus with the turbulent flux at a depth more than 25 m, i.e. outside the surface layer of the maximum salinity gradients.

Conclusions.

1. With account of eddy viscosity and diffusion the wave salt flux u3S

differs from zero, and in the presence of current it differs somewhat smaller than in its absence.

2. In the absence of current the wave salt flux is comparable in absolute value with the turbulent flux out of the subsurface layer (where the salinity gradients are maximal).

3. The component of Stokes drift velocity (transversal to the wave propagation direction) differs from zero when the current is taken into account.

REFERENCES

1. LeBlond, P., Mysak, L., 1981, "Volny v okeane", Moscow, Mir, part. 2, 363 p. (in Russian).

2. LeBlond, P.H., 1966, "On damping of internal gravity waves in a continuously stratified ocean", J. FluidMech., vol. 25, part. 1, pp. 121-142.

3. Ostrovsky, L.A., Soustava, I.A., 1979, "Verkhniy peremeshannyy sloy kak stok energii vnutrennikh voln [The upper mixed layer of the ocean as an energy sink of internal waves]", Okeanologiya, vol. 19, no. 6, pp. 973-981 (in Russian).

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

71

4. Borisenko, Yu.D., Voronovich, A.G. & Leonov, A.I. [et al.], 1976, "K teorii nestatsionarnykh slabonelineynykh vnutrennikh voln v stratifitsirovannoy zhidkosti [On the theory of non-stationary weak nonlinear internal waves in stratified fluid]", Izv. AN SSSR, Fizika atmosfery i okeana, vol. 12, no. 3, pp. 293-301 (in Russian).

5. Grimshaw, R., 1977, "The modulation of an internal gravity wave packet and the resonance with the mean motion", Stud. In Appl. Math., vol. 56, pp. 241-266.

6. Grimshaw, R., 1982, "The effect of a dissipative processes on mean flows induced by internal gravity-waves packets", J. FluidMech. vol. 115, pp. 347-378.

7. Ivanov, V.A., Samodurov, A.S. & Chuharev, A.M. [et al.], 2008, "Intensifikatsiya vertikal'nogo turbulentnogo obmena v rayonakh sopryazheniya shel'fa i kontinental'nogo sklona v Chernom more [Intensification of the vertical turbulent transfer near the interface between the shelf and the continental slope in the Black Sea]", Dop. NAN Ukrainy, no. 6, pp. 108-112 (in Russian).

8. Samodurov, A.S., Nosova, A.V. & Slepyshev, A.A., 2011, "Fizicheskie mekhanizmy intensifikatsii vertikal'nogo obmena v zone sopryazheniya shel'fa i svala glubin [Physical mechanisms of an intensification of a vertical transfer in a zone of interface of a shelf and the continental slope]", Ekologicheskaya bezopasnost'pribrezhnoy i shel'fovoy zon i kompleksnoe issledovanie resursov shel'fa, iss. 25. vol. 2, pp. 190-203 (in Russian).

9. Slepyshev, A.A., 1997, "Protsessy perenosa, obuslovlennye slabonelineynymi vnutrennimi volnami pri nalichii turbulentnosti [Transport Processes Induced by Weakly Nonlinear Internal Waves in the Presence of Turbulence]", Izv. RAN. Fizika atmosfery i okeana, vol. 33, no. 4, pp. 536-548 (in Russian).

10. Vishik, M.I., Lyusternik, L.A., 1957, "Regulyarnoe vyrozhdenie i pogranichnyy sloy dlya lineynykh differentsial'nykh uravneniy s malym parametrom [Regular Degeneration and Boundary Layer for Linear Differential Equations with a Small Parameter]", Uspekhi matematicheskikh nauk, vol. 12, no. 5 (77), pp. 3-122 (in Russian).

11. Cherkesov, L.V., 1980, "Gidrodinamika voln", Kiev, Naukova Dumka, 259 p. (in Russian).

12. Longuet-Higgins, M.S., 1969, "On the transport of mass by time varying ocean current", Deep-SeaRes., vol. 16, no. 5, pp. 431-447.

13. The report on works in 44th cruise (3rd stage) of research vessel "Michael Lomonosov" August, 7th - September, 15th 1985, Sevastopol, MGI NAN Ukrainy, 1985, vol. 1, 135 p. (in Russian).

72

PHYSICAL OCEANOGRAPHY NO. 1 (2015)

i Надоели баннеры? Вы всегда можете отключить рекламу.