CC BY
17Display of the ^-effect in the Black Sea Two-Layer Model
A.A. Pavlushin, N.B. Shapiro, E.N. Mikhailova
Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russian Federation
e-mail: pavlushin@mhi-ras.ru
The research is a continuation of a series of numerical experiments on modeling formation of wind currents and eddies in the Black Sea within the framework of a two-layer eddy-resolving model. The main attention is focused on studying the ^-effect role. The stationary cyclonic wind is used as an external forcing and the bottom topography is not considered. It is shown that at the ^-effect being taken into account, the Rossby waves propagating from east to west are observed both during the currents' formation and at the statistical equilibrium mode when the mesoscale eddies are formed. In the integral flows' field the waves are visually manifested in a form of the alternate large-scale cyclonic gyres and zones in which the meso-scale anti-cyclones are formed. This spatial pattern constantly propagates to the west that differs from the results of calculations using the constant Coriolis parameter when the spatially alternate cyclonic and anti-cyclonic vortices are formed, but hold a quasi-stationary position. The waves with the parameters of the Rossby wave first barotropic mode for the closed basin are most clearly pronounced. Interaction of the Rossby waves with large-scale circulation results in intensification of the of the currents' hydrodynamic instability and in formation of the mesoscale eddies. Significant decrease of kinetic and available potential energy as compared to the values obtained at the constant Coriolis parameter is also a consequence of the eddy formation intensification.
Keywords: the Black Sea, eddy-resolving model, numerical experiment, ^-effect, Rossby waves.
DOI: 10.22449/1573-160X-2016-5-3-23
© 2016, A.A. Pavlushin, N.B. Shapiro, E.N. Mikhailova © 2016, Physical Oceanography
Introduction. To investigate the effect of various factors (internal and external) on the formation and evolution of the hydrophysical fields in the Black
Sea, a series of targeted numerical experiments with eddy-resolving two-layer isopycnic model was carried out [1, 2]. The model is based on a system of primitive equations in the Boussinesq, hydrostatic and yS-plane approximations. The sea comprising of two layers with the upper layer density p = const and the lower layer one p2 = const. The layers are not mixed with each other. The friction between the layers can be taken into account in the model. The tangential wind stress t is set on the sea surface. The energy sink is carried out at the expense of the bottom friction and horizontal turbulent viscosity.
The equations of motion and continuity, vertically integrated within each layer, have the following form
(U )t + (uU )x + (v1U1 )y - V = ghZx +T - RX + A, V(h1Vu1),
(V )t + (u1 V )x + (vV1 )y + fU1 = (1)
= gh^Zy +Ty - Ry + a,lV(h1 VV1),
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
3
(U2) +{u2U2)x +(v2U2)y -= gh2Z + g\(h1 )x + RX -RX + Atv(h2Vu2), (V2 )t +UV2 )x + )y + fU 2 = = gh2^y + g \ (hi )y + Ry - Rb + A,v(h2vv2), (1)
(hi ) +(Ui )x +(Vi )y = 0,
(h2 ) +(U 2 )x +(V2 )y = 0,
where indices 1 and 2 denote the upper and lower layer correspondingly; the lower indices x, y and t designate differentiation; ui, vi are the /-layer horizontal
components of the current velocity; hi, h2 are the layer thickness; U = uihi, Vi = Vjhj are the flow components; Z is the sea level; R^i = ra (ui - U2), Ry = ra(v - V2) are the components of the friction force between the layers; Rb = rbU2, Rb = rbV2 are the bottom friction force components; ra, rb are the constant coefficients; f = f0 + f3y is the Coriolis parameter, f0 = 10-4 1/s, P = 2 -10-13 1/(cm^s); g = 980 g-cm/s2 is the free fall acceleration; g ' = g (p2 - pi)/p2 ; T, Ty are the tangential wind stress components; Ai is the horizontal eddy viscosity coefficient.
The integral continuity equation in the rigid lid approximation terminates the equations. It permits to introduce the stream function y for the total flows:
Ui + U 2 =-¥y, Vi + V2 =Yx .
At the side basin boundaries the no-slip conditions are set. River runoff into the sea and water exchange through the straits are not taken into account. Initially, the water is at rest, the interface layer and the sea surface are horizontal.
In the finite-difference model representation a time two-layer numerical scheme is applied. It is based on the 5-grid box method (in the terminology of Arakawa), the implicit approximation of the Coriolis force and friction forces on the section and the bottom surface. Advective members in the equations of continuity are approximated by the first accuracy order (directed differences) scheme and in the equations of motion - by the second accuracy order scheme of the (Lax - Wendroff).
The first stage of the works involved the study of the wind circulation formation mechanisms, due to the fact that the wind is a major factor in determining the Black Sea seasonal and interannual variability. Taking into account the characteristics and specifics of the model applied, it is possible to identify the important factors that influence the dynamics of currents, such as the wind spatio-temporal variability, the horizontal turbulent viscosity, the bottom friction force, the friction force between the layers, the basin shape, especially the coastline, yS-effect, the bottom topography.
4
PHYSICAL OCEANOGRAPHY NO. 5 (20i6)
The paper [1] presents the analysis of the results of experiments with different coefficients of horizontal turbulent viscosity Ai and bottom friction n at the constant Coriolis parameter f = 10-4 1/s, corresponding to the Black Sea latitude. It is shown that satisfactory results in the model are obtained using the coefficients Ai ~ 105 cm2/s and n in the range of 0.01 - 0.1 cm/s. Under these parameters in the basin alongside with the stable quasi-stationary cyclonic gyres the mesoscale anticyclonic eddies associated with the hydrodynamic instabilities of currents sporadically appear.
In the paper [2] the effect of the basin shape on circulation was studied. The results of calculations of current fields in basins with different configurations in the absence of yS-effect were analyzed. In the stretched basins, provided the bottom friction is quite weak, the large-scale circulation is divided into individual cyclonic eddies, even without the protruding coastline elements available. The number of these eddies depends on the ratio of the length of the basin to its width.
In the papers [1, 2] the energy inflow in the sea under the stationary wind t (t) = const is also shown to be regulated by the oscillatory process [3] associated with changes in the work of the tangential wind stress WT = |t| • \u1\ • cos(«). This is due to the transformation of the surface current field u1 by mesoscale eddies generated by hydrodynamic (barotropic and/or baroclinic) large-scale circulation instability. The period of oscillation in the yS-effect absence is mainly determined by the lifetime of mesoscale eddies.
In the present article the research of the yS-effect role on the formation of wind-driven circulation in the Black Sea is carried out. According to the theory [4 - 6] latitude change of the Coriolis parameter leads to appearance of the planetary Rossby waves, which can have a very wide range of manifestations depending on various conditions. The Black Sea specificity is that it is almost a closed basin. Its size is less than the barotropic (external) Rossby deformation radius, but more than baroclinic (internal) deformation radius.
There are a number of works that are directly or indirectly related to the modeling and description of processes in the Black Sea, connected with the manifestation of planetary waves [7 - 14]. First of all, the works of E. Stanev, N. Rachev [7, 8] should be paid attention to. There the analysis of the numerical simulation of wind circulation in the Black Sea applying the level eddy-resolving Bryan - Cox model is described. The movement in the sea as well as in our case was excited by stationary wind. As a result of computations performed in the Black Sea long-term fluctuations were obtained. They were described as barotropic Rossby waves generated in a closed basin, whose dimensions are smaller than the outer radius of deformation. According to the authors of [8], the Rossby waves are the dominant form of wave motion in the Black Sea and the variability of the circulation associated with planetary modes can be compared in some cases with the variability caused by baroclinic instability.
Similar results were obtained in the paper [9], dedicated to the modeling of wind circulation in the Black Sea applying quasi-geostrophic model. During the numerical integration under the stationary wind influence a problem solution came PHYSICAL OCEANOGRAPHY NO. 5 (2016) 5
out on a quasi-periodic mode with the oscillations in the total flow fields with periods of 43 and 83 days. According to opinion of the authors, the first of these periods is associated with the manifestation of baroclinic instability of the main current, and the second one can be caused by manifestation of barotropic Rossby waves.
In [i0] based on the processing of satellite observations of the Black Sea level changes the long-period oscillations propagating in the western direction were monitored. The authors explain this process by radiation of the baroclinic Rossby waves, caused by seasonal fluctuations of the wind vorticity from, the eastern coast of the sea.
In [ii] based on the analysis of hydrophysical and morphometric characteristics of the Black Sea, the conclusion about the possible existence of the Rossby waves in the basin with periods of 80 - 200 days and phase velocities of 2 - 8 cm/s is drawn.
As the observational data indicating the presence of planetary waves in the Black Sea basin, the works [i2, i3] can also be mentioned. They describe the anticyclonic mesoscale eddies generated in the Black Sea during summer. Monitoring of the eddies showed that during the observation period, they had a velocity component of ~ 2 km/day, directed to the west, which may indicate their relation to the Rossby waves.
Generally, the features of the appearance and propagation of planetary waves in the Black Sea have not been studied enough.
Description of the experiments. As it was previously mentioned, the experiments considered in the present paper, are devoted to /-effect and its effect on the formation of wind circulation in the Black Sea. The experiments do not take into account the bottom topography and friction between the layers. t wind power
is set stationary at a constant cyclonic vorticity 0,5 • i0-7 N/m3. The sea depth is set equal to H = 2200 m. The upper layer initial thickness is h = i75 m, the Coriolis parameter on the southern boundary is /3 = 2-10_i3 s_icm_i, the Rossby parameter is /3 = 2-10_i3 s_icm_i, horizontal eddy viscosity coefficient is Ai = i05 cm2s_i.
In previous studies [i, 2] it was shown that the results of solving the problem under the constant Coriolis parameter are very sensitive to the value of bottom friction, while rb ratio variation range is quite wide. Therefore, computations with different values of the rb coefficient were carried out.
For the convenience of description, in the present paper we introduced the following alphanumeric designations of the experiments. In the numerical experiments, indicated by the letter B, /3 > 0, experiments with /3 = 0 are indicated by the letter A. The paper discusses the results of six experiments: Bi, B2, B3 and Ai, A2, A3. In addition to /-parameter, the difference between them is in the use of different values of the rb bottom friction coefficient: in Bi, Ai, B3, A3 experiments rb = 0.0i cm/s; B2, A2 ones - rb = 0.icm/s.
Furthermore, in B3 and A3 experiments the model basin shape ("sausage"-type), the same as in Fi experiment in [2], is used except the actual Black Sea
6
PHYSICAL OCEANOGRAPHY NO. 5 (20i6)
configuration. This makes it possible to demonstrate the influence of the effect of the basin configuration and coastline features on the formation of currents, eddies and Rossbi waves.
The computations were performed on a 3 * 3 km square grid with a time step of At = 3 min. The use of a high spatial resolution compared to the previous works [i, 2], which applied a 4 * 3.5 km grid and time step of At = 6 min., permitted to obtain the smaller eddies with an increase in their quantity in the study area
Duration of the computations in each experiment was not less than i0 years. The experiment was stopped after the decision came out on a statistically equilibrium mode [2], in which the averaged time characteristics of the model has changed little over time.
Results of the numerical experiments. The performed computations resulted in obtaining instantaneous and averaged over time fields of the upper layer
thickness hi, the sea level Z, currents in the upper and lower layer ui, u 2 and integrated stream function y.
In addition, graphs of the time changes of area averaged values were constructed. These values were the kinetic energy of the upper and lower layer KEi KE2, available potential energy DPE and tangential wind stress work Wt calculated according to the following formulas
KEi = (pAu + vi2) /2), KE2 = ph2(u2 + v22)/2), DPE = (pg'(hi -h0)2/2), WT=(pi(u1zx + viTy)),
where hi, h2 are the upper and lower layer thickness; ui, vi, u2 and v2 horizontal components of the current velocities in the respective layers; g' = g (p2 - pi)/p2 = 3,2 cm/s2; Tx, Ty are the tangential wind stress horizontal components; the angle brackets denote averaging over the area.
Results of the computations obtained under /3 > 0 were compared with results of the experiments where /-parameter was assumed to be zero.
Fig. i shows the instantaneous spatial distributions of the features in Bi and Ai experiments. Ai experiment (/ = 0) is characterized by the presence of all the fields in the two large-scale cyclonic gyres (Fig. i, a - e). One cycle is located in the western part of the sea, the other - in the east one. These sub-basin cyclones define large-scale circulation in the basin. Between the cyclones and the coast and between mesoscale cyclones anticyclonic eddies occasionally appear. They are
clearly visible in the hi, Zand ui fields (Fig. i, f - h). Their lifetime is about 50 days.
PHYSICAL OCEANOGRAPHY NO. 5 (20i6)
7
There are clearly distinguished two cyclonic eddies averaged over several in the fields M and y (Fig. 2, a, b). They were also observed in the instantaneous distributions. This testifies their stationary nature. With regard to mesoscale anticyclonic eddies, due to their non-stationarity in time and space, they are weakly manifested in the averaged field y and do not appear at all in the field M.
0 200 400 600 800 1000 A", km 0 200 400 600 800 I000.\',km
Fig. 1. The instantaneous spatial distributions of the features in the experiment A1 (left) u B1 (right) 8 PHYSICAL OCEANOGRAPHY NO. 5 (2016)
AI experiment B1 experiment
0 200 400 600 800 1000 X, km 0 200 400 600 800 1000 A', km
Fig. 2. Average distributions hi and y in the experiments A1 (left) and B1 (right)
In the experiment (B1) with /-effect taken into account the large-scale cyclones in the sea central part are also found in the instantaneous distributions hi,
Zand ui (Fig. 1, a - c) (B1), but, unlike A1 experiment, they are weaker, and their number is more than two. Therefore, we can rather speak of the presence of a large-scale cyclonic gyre with several peaks in the central part of the given basin. Mesoscale anticyclonic eddies are observed along the coast.
The single cyclonic and anticyclonic eddies are more distinguished in the
fields y and u 2 in B1 experiment than in the upper layer.
Eddy formations, available in instantaneous fields, are constantly moving. It leads appearance of a vast cyclonic area t in the medium fields M and y in the central part of the basin (Fig. 2. c, d). As a consequence of the intensification of western currents, there is a slight thickening of the isolines h1 and y near the western coast.
In both experiments considered the qualitative coincidence of current fields in the upper and lower layers, indicating that the process barotropization is noteworthy. The direction of currents in the upper and lower layers is substantially the same. The velocity of currents in the lower layer is about 40 % of the one in the upper layer.
Integrated stream function y, which characterizes the barotropic currents,
correlates well with the field u 2, and the isolines h are close to the current lines
U1and the isolines Z This peculiarity in the spatial distributions of the features permits to apply the field hi to visualize the circulation in the upper layer and y -in the lower one.
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
9
0 200 400 600 800 1000A', km 0 200 400 600 800 1000A', km
Fig. 3. Consequent fields h1 (left) and y (right) with 15 day period in B1 experiment. The dotted line represents the section Y = 270 km
As it has already been mentioned, the eddy formations occurring under / > 0, do not have a stationary position. If the sequence of instantaneous fields h1 and y, created for B1 experiment, is considered (Fig. 3), the movement of eddy formation has mainly the western direction, which is not observed in A1 experiment.
In graphs M and y (Fig. 4), plotted for the latitudinal section, going through the central sea part along the line Y = 270 km, the waves moving at 8 cm/s velocity with 80 - 100 day period in the western direction, can bee seen. In field y these
10
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
waves are better manifested than in Mone, as in the upper layer the currents and mesoscale eddies associated with baroclinic instability are superposed on them.
у/, *106 nvVs h\, m
100 200 300 400 500 600 700 800 900 A",km 100 200 300 400 500 600 700 800 900 X,km
Fig. 4. Time graphs y (left) and hi (right) on the section Y = 270 km in B1 experiment
The yS-effect role is also manifested in the behavior of energy. Fig. 5, a, b, c, d show the graphs ofthe area averaged KE\, KE2, DPE, W, calculated for B1 and A1 experiments.
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
11
Fig. 5. Energy graphs of KEi, KE2, DPE in B1 (a) and A1 (d) experiments; the tangential wind stress work Wt in B1 (b) and A1 (e) experiments; mean vorticity B1 (c) and A1 f
Fig. 5 a, d, show the solution reaches a statistically equilibrium mode, the characteristics change within the relatively limited certain mean values in the both experiments.
In B1 experiment under / > 0 the KE1, KE2 and DPE energy values resulted significantly smaller than in A1 experiment. The work Wt, which provides the energy inflow in the sea, also appeared smaller. As Wt is determined by oscillatory process, depending on the stability of currents in the upper layer, it can be assumed that the /-effect influences the currents, making them more instable.
It can be seen if we calculate separately the average area values of cyclonic and anticyclonic relative vorticity in the different layers of the model. It is necessary to sum over areas cyclonic and anticyclonic vorticity values in each layer separately and divide by the corresponding areas. The resulting mean values for
cyclonic and anticyclonic vorticity will be denoted as 44, £2 44 and 42 :
£ =-C iff^-— \xdy, =±r iff^-^W , a = 1, 2),
Sf fJ{dx dy, y' Sf JJ^cX dF J
12
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
where Ci, Ai are areas occupied by cyclones and anticyclones respectively; SC, Sf are areas under cyclones and anticyclones; i is number of the layer.
Due to the adhesion conditions on the lateral borders of the total cyclonic and anticyclonic vorticity in the basin are equal in absolute value: jjfSf + jjASA = 0 . But the
area under the anticyclones and cyclones may be different, so the mean values of the cyclonic and anticyclonic vorticity in absolute value are not equal. The values jjf and ¿A can be applied as an integral vorticity field characteristics in the sea.
Fig. 5, c, f show the graphs of jjC , j j and j, calculated in accordance with the instantaneous current velocity values in B1 and A1 experiments. As it can be seen from the aforementioned figures, in moth cases the mean anticyclonic vorticity in the upper layer is 1.5 times more than the cyclonic one in absolute value. Such pattern is obviously caused by cyclonic wind effect. For the lower layer the same vorticity characteristics are almost equal in absolute value. It is noteworthy that both in experiments respective quantities j, j jj1 and j take close values unlike energy and tangential wind stress work. This means that in B1 experiment under в> 0 the mean velocity field vorticity, similar by the magnitude in A1 experiment (в = 0) is achieved at a lower energy level.
In B2 and A2 experiments in order to relieve barotropic component of currents was set (гь = 0.1 cm/s) more by an order of the bottom friction. As a result, in A2 experiment (Fig 6, a - e) under the constant Coriolis parameter the large-scale circulation separation into individual cyclonic eddies, as in A1 Experiment, did not take place. There was a large cyclonic gyre observed in the sea. It occupied the central part of the basin. On its periphery mesoscale anticyclones periodically arose. They moved along the coast in the direction of the main stream. The areas of anticyclonic vorticity (hydrodynamic instability of currents) were mainly tied to the protruding features of the coastline.
Taking into consideration в-effect in В2 experiment (Fig. 6, f - j), in comparison with А2 one, the separate eddy formations began to be more clearly identified in the instantaneous fields, especially in у and U2 fields (Fig. 6, i, j). Although in general circulation changed little.
In the mean fields у and M in В2 experiment (Fig. 7) the west intensification of currents is visible and the area of anticyclonic vorticity at the eastern boundary of the basin appeared. In addition, in the averaged field M the cyclonic circulation in the basin center has two centers ("Knipovich glasses"), which was not observed in B1 experiment.
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
13
Fig. 6. Instantaneous distributions of the characteristics in A2 (left) and B2 (right) experiments
Like in B1 experiment, the eddy formations represented in the instantaneous fields moved in the western direction. It is well observed on the time diagrams y and M, built for the section Y = 270 km (Fig. 8). In diagram y the waves moving from the east to the west with the phase velocity of 7 - 8 cm/s are clearly identified. In the field M (Fig. 8, b) the phase velocity of the observed waves is less and equals ~6.5 cm/s.
14
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
0 200 400 600 800 1000 X, km 0 200 400 600 800 1000 Ji,km
Fig. 7. Mean distributions hi and y in A2 (left) and B2 (right) experiments
Fig. 9 shows graphs of KEi, KE2, DPE, W, tf , 42c # and , plotted for A2 and B2 experiments. It is clear that as compared with A1 and B1 experiments (Fig. 5, a, d) the more intense bottom friction leads to a significant decrease of KE2 and increase of DPE, with practically unchanged KEi.
/-effect in B2 experiment (Fig. 9, a) caused the energy decrease (especially DPE) comparing with A2 experiment (Fig. 9, d), at that the vorticity features (Fig. 9, c, f) did not changed significantly.
Summarizing the results of B1 and B2 experiments, we can assume that in the Black Sea under the cyclonic wind influence the role of /-effect leads to the formation of long-period waves, which are manifested in the form of alternating eddies moving in the western direction. These waves can be most clearly traced in the field y.
Similar results have already been obtained and are described in the papers [7 -9]. In [8] to determine the nature of the observed processes, the authors compared the model results with the normal modes of own vibrations of a rectangular basin with a homogeneous liquid, with the length, width and depth roughly corresponding to the Black Sea parameters.
The solution of the system of equations of natural vibrations in /-plane for the rectangular two-layer ocean geostrophic approximation is given in [10]. In this case the system of equations has the following form
|Afa + Pd-fa- "Ei|fa -fa) = 0,
dt dx dt (2)
d A ¥2 + P^T-- E2 -YWi) = a
dt dx dt
where yi, y2 are stream functions, for vertically averaged velocities in the upper
and lower layer; Ei = fV{ghz); e2 = fV(sh2s); s = (p2-Pi)/p>; y = PiIPI
PHYSICAL OCEANOGRAPHY NO. 5 (20i6)
i5
y/, *106m3/s hi, m
100 200 300 400 500 600 700 800 900 A"km 100 200 300 400 500 600 700 800 900X,km Fig. 8. Time diagrams y (left) and h1 (right) on the section Y = 270 mm in B2 experiment
16
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
Fig. 9. Graphs of the KEi, KE2 and DPE energy in B2 - a and A2 - d experiments; tangential wind stress work Wt in B2 - b and A2 - e experiments; mean vorticity in B2 - c and A2 - f experiments
Multiplying the second equation of the system (2) by an arbitrary constant a, and adding it to the first one, we obtain the following
d_
dt
A[ + a[) + p-d(wi +a\ff2)-Ei-^[[1(1 -aa{)y) + V2(aao -1)]= 0
d
dx
dt'
(3)
where a0 = h1jh2 .
Then we are to define a new dependent variable F and an arbitrary constant G: F = [1 +a[2 =[[ (l - aa0y) + [2 (aa0 -1)] G , (4)
as a result we obtain the following equation of the normal vibrations
— AF + - E dF = 0. dt dx G dt
(5)
From (3) we can find the values of the arbitrary constants a and G:
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
17
G = 1 , . (6)
1 -aa0y 1 -aa0Y
To determine a we obtain a quadratic equation as follows
2 a0 -1 1
a + —-a--= 0,
2aoY aoY
from which we find
a0 -1 _
a, 2 =—0—+.
2aoY
(ao -1)2 + 4aoY (7)
22 4ao Y
Taking into account F = y/l + aiy2 (i = 1, 2) the equation (5) is rewritten in the following way:
d_
~dt
( E ^
G
i J
F = 0. (8)
dx
The expression (8) is a well-known equation of planetary waves [15]. It defines two normal fluctuations (barotropic and baroclinic modes) with frequencies a in a rectangular basin with the a x b dimensions at Fi = 0 on the boundary expressed by the formula below
a.= — 1 2
ß( m V n V E^-1X 2 + —^ + —
. a2 b2 G . v t J
(9)
where i = 1, 2; m and n are whole numbers.
Taking into account (6) and (7), two values for EjGt can be found:
El = f2 EL = f2 (h + ^2)
G1 g(hl + h2 y G2 g£hA
The first expression defines a value for the barotropic oscillation mode with an equivalent depth of H = h + fe, the second one - for the baroclinic modes with
h h
equivalent depth of H = s-^—.
(h + h2)
Solution of the equation (8) for a rectangular basin with the a x b dimensions under Fi = 0 on a solid boundary was first obtained by Longuet-Higgins [15] and has the following form
F = sin\ ^^x I sin\ ^^y I cos(y x + at), (i = 1, 2), (10)
I a J \ b J where x and y are coordinates; m and n are whole numbers; yt .
The expression (2) describes the Rossby waves, which are the normal modes of natural oscillations in the basin [6, 7]. Under i = 1 the barotropic mode is
18 PHYSICAL OCEANOGRAPHY NO. 5 (2016)
obtained, under i = 2 - the baroclinic one. The phase velocity of the wave ci =-aji _1, the distance between the moving nodes (half wavelength) Ai = 2ncrip~l and the wave period Ti = 2n/ci.
For the rectangular basin with the dimensions of 1125 * 285 km and the depth H = 2200 m, with two water layers having the parameters, analogical to the experiments considered, for the principal (m = 1, n =1) barotropic normal mode we obtain: C1 = -7.7 cm/s, M = 276 km, T1 = 82.8 days, for the principal baroclinic mode: C1 = -0.5 cm/s, M = 69 km, T1 = 331 days. The wave, related to the main barotropic normal mode has maximum phase velocity. Baroclinic normal modes under the same values of m and n are have smaller phase velocities and wavelengths, than the barotropic modes.
In B1 and B2 experiments, phase velocities of the observed waves (Fig. 4, 8) were obtained close enough to the above calculated velocity of the main barotropic modes of natural oscillations in a rectangular basin. But such a comparison may not look correct due to the strong differences in the shapes of the compared basins. In addition, the calculations using the same formula (10) with different sizes of a rectangular basin showed that the basin width b (the size according to the 7-axis) has the greatest effect on the results, and in case of a real configuration it is very difficult to determine the equivalent width of the sea. Therefore additional B3 experiment was carried out in an elongated basin with constant width and boundaries rounded off to the east and west (Fig. 10). As shown in [2], this form allows you to minimize the nonlinearity effect associated with the basin asymmetry and the coastline roughness. The length, width and depth matched the dimensions of the rectangular basin, the normal modes were calculated for. In addition, to compare and control A3 experiment was carried out, which differs from B3 one by zero value of /-parameter.
Instantaneous distributions M and y, obtained after the problem entrance in the statistical equilibrium mode in B3 and A3 experiments and giving an idea of the circulation in the upper and lower layer, are shown in Fig. 10, a, b, d and e. Distributions of the other characteristics are not presented, as there is a good cross-
correlation between fields M, Z, U1 and y, u 2, as well as in the previously considered pairs of the experiments. As it can be seen in Fig. 10, a - d, in A3 experiment the large-scale circulation was divided into three rounded shape cyclones with a diameter equal to the basin width. Two zones are located between the cyclones. There are mesoscale anticyclonic eddies periodically encountered (Fig. 10 a, b). In the averaged field M (Fig. 10, c, d) cyclonic circulations are only observed. Some anticyclonic eddies are not expressed because of their unsteadiness. In y field (Fig. 10, d) two anticyclonic vorticity areas are observed between the three stationary cyclones.
A3 experiment peculiar feature is the presence of the central symmetry in the spatial distribution of all the characteristics - the instantaneous and averaged ones. In B3 experiment taking into account /-effect the circulation separation in the individual eddies takes place (Fig. 10 d, e), but there is no symmetry of the fields. In the medium fields M and y (Fig. 10, f, h) the separate eddy structures observed
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
19
in the instantaneous fields, are not manifested. A slight thickening of the contour as a result of the intensification of western currents is observed near the western coast.
Fig. 10. Instantaneous and mean fields in A3 experiment: hi - a; hi - b; y - c; ^ - d, Instantaneous and mean fields in B3 experiment: hi - e; hi -f; y - g; ^ - h
In B3 experiment all the eddy formations move in the western direction. This is clearly shown in Fig. ii (left). There the successive fields y, built in the period of 15 days, are demonstrated. If they are compared with the main barotropic normal mode distribution, calculated with the same frequency (Fig. ii, right), great similarity will be found. Also, time diagrams y and hi, built for the latitudinal section passing through the central axis of the basin (Fig. i2), are in a very good agreement with the barotropic normal mode. All the diagrams clearly show wave with a period of 90 days, moving in the western direction with the phase velocity of 7.5 cm/s. Thus, it is possible to make presumptive conclusion that by accounting for the influence of S-effect in Bi, B2 and B3 numerical experiments the Rossby waves were obtained. They correspond to the normal barotropic modes of natural oscillations in a closed basin. These waves can be best seen in the field y, which characterizes the barotropic component of currents. In diagram hi (Fig. i2, c) you can see traces of the movement of mesoscale anticyclonic eddies on the barotropic
20
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
Rossby wave background in the areas of positive values of y, which correspond to anticyclones.
Fig. 11. Consistent spatial distributions y with the period of 15 days in B3 experiment (left) and normal mode in a rectangular basin with m =n = 1 (right). The dotted line indicates the latitudinal sections along the central axis of the basins
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
21
Fig. 12. Consistent spatial distributions at the latitudinal sections along the central axis of the basin for y, corresponding to the normal mode (left) in B3 experiment (center), hi in B3 experiment (right)
As for the Rossby waves caused by baroclinic modes of the normal vibrations, which, according to theory, should also be present in the solution, they were failed to be identify explicitly in the considered experiments. Most likely, this is due to the fact that the self-oscillating process, changing the sink of energy from the wind and resulting in the oscillation of the observed waves, has the frequency the barotropic Rossby waves predominantly occur at. Theoretical ocean response research on the different types of external exposure is considered in [14]. It is shown that the external effect at the intervals of several weeks - months mostly results in the barotropic Rossby waves. But for the baroclinic waves to appear the oscillating force frequency should be more than a year.
Conclusions. A good match of features of the waves, obtained in the experiment, with the analytic solution of the planetary wave equation gives reason to believe that in the study area there are Rossby waves corresponding to the main normal barotropic mode of the natural basin oscillations. Under steady external effect the initial velocity field oscillations are due to hydrodynamic instability of currents and are further supported by the self-oscillating process that regulates energy sink into the sea by the wind. The amplitude of Rossby wave baroclinic mode was small compared with the amplitude of the barotropic mode, and was difficult to distinguish in the obtained fields.
22 PHYSICAL OCEANOGRAPHY NO. 5 (2016)
Analysis of the behavior of energy characteristics showed that the /-effect reduces the tangential wind stress work, which in its turn determines the lower level of the energy components calculated in the model. Rossby waves have a stimulating effect on the processes of hydrodynamic instabilities that lead to the appearance of opposite wind vorticity mesoscale eddies.
On the basis of the performed experiments we can agree with the conclusion of the authors of the article [8] that the /-effect and associated Rossby waves in the Black Sea can have a significant impact on the formation of large-scale circulation, comparable to the baroclinic instability effect.
REFERENCES
1. Pavlushin, A.A., Shapiro, N.B., Mikhaylova, E.N. & Korotaev, G.K., 2015, "Two-layer eddy-resolving model of wind currents in the Black Sea", Physical oceanography, no. 5, pp. 3-12.
2. Pavlushin, A.A., Shapiro, N.B., Mikhaylova, E.N, 2016, "Effect of the basin shape on formation of the Black Sea circulation", Physical oceanography, no. 2, pp. 3-15.
3. Seidov, D.G., 1989, "Sinergetika okeanskikh protsessov [Synergetics of the ocean processes]", Leningrad, Gidrometeoizdat, 288 p. (in Russian).
4. Gill, A., 1986, "Dinamika atmosfery i okeana. T. 1 [Atmosphere and ocean dynamics. Vol. 1]", Moscow, Mir, 396 p. (in Russian).
5. Kamenkovich, V.M., Monin, A.S., 1978, "Fizika okeana. T. 2-Gidrodinamika okeana [Ocean physics. T2-Ocean hydrodynamics]", Moscow, Nauka, 435 p. (in Russian).
6. Kamenkovich, V.M, Koshlyakov, M.N. & Monin, A.S., 1987, "Sinopticheskie vikhri v okeane [Synoptical eddies in the ocean]", Leningrad, Gidrometeoizdat, 510 p. (in Russian).
7. Rachev, N.H., Stanev, E.V., 1997, "Eddy processes in semi-enclosed seas: A case study for the Black Sea", J. Phys. Oceanogr., no. 27, pp. 1581-1601.
8. Stanev, E.V., Rachev, N.H., 1999, "Numerical study on the planetary Rossby waves in the Black Sea", J. Mar. Sys., no. 21, pp. 283-306.
9. Demyshev, S.G., Korotaev, G.K. & Kuftarkov, A.Yu., 1996, "Analiz nachal'nogo perioda prisposobleniya geofizicheskikh poley Chernogo morya k novoy chislennoy modeli [Analysis of the initial adaptation period of the Black Sea geophysical fields to the new numerical model]", Izv. RAN. FAO, vol. 32, no. 5, pp. 635-644 (in Russian).
10. Korotaev, G.K., Saenko, O.A. & Koblinsky, C.J., 2001, "Satellite altimetry observation of the Black Sea level", J. Geophys. Res., iss. 106, no. C1, pp. 911-933.
11. Blatov, A.S., Bulgakov, N.P. & Ivanov, V.F. [et al.], 1984, "Izmenchivost' gidrofizicheskikh poley Chernogo morya [Variability of hydrophysical fields of the Black Sea]", Leningrad, Gidrometeoizdat, 239 p. (in Russian).
12. Latun, V.S., 1989, "Antitsiklonicheskie vikhri v Chernom more letom 1984 goda", Morskoy gidrofizicheskiy zhurnal, no. 3, pp. 27-34 (in Russian).
13. Latun, V.S., 1990, "Energosnabzhenie glubokikh antitsiklonicheskikh vikhrey Chernogo morya [Power supply of the Black Sea deep anticyclonic eddies. Complex oceanographic research of the Black Sea: Proceedings]", Kompleksnye okeanograficheskie issledovaniya Chernogo morya: Sb. nauch. tr., pp. 10-21 (in Russian).
14. Safronov, G.F., 1985, "Vozbuzhdenie dlinnykh voln v okeane krupnomasshtabnymi izmeneniyami v pole kasatel'nogo napryazheniya vetra [Oscillation of long waves in the ocean by large-scale changes in the tangential wind stress]", Moscow, Gidrometeoizdat, 108 p. (in Russian).
15. Longuet-Higgins, M.S., 1965, "Planetary waves on a rotating sphere II", Proc. Roy. Soc. A., vol. 284, iss. 1396, pp. 40-68.
PHYSICAL OCEANOGRAPHY NO. 5 (2016)
23
