Structure of the Black Sea Currents Based on the Results of the LADCP Observations in 2004-2014
A.N. Morozov1*, E.M. Lemeshko2**, S.A. Shutov1, V.V. Zima1,
D.V. Deryushkin1
1Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russian Federation
*e-mail:[email protected] 2Black Sea Hydrophysical Polygon, Russian Academy of Sciences, Katsiveli, Russian Federation **e-mail:evgeny. [email protected]
The article presents the core results of lowered acoustic Doppler current profiler (LADCP) applied in the Black Sea within 2004-2014. Average features of the vertical structure of currents in the upper 600 m layer are outlined. The linear relation between vertical distribution of kinetic energy isopycnically averaged over the station ensemble and density is established. This relation is traced till the isopycn occurrence depth with density of 17 kg/m3 (~500 m). The mean current velocity has near constant value about ~ 4 cm/s in lower depth layers (>500 m), which is in a good agreement with autonomic buoys data. The averaged profile of current velocity vertical shears shows two well-separated maxima in the seasonal and main pycnocline layers. LADCP derived velocity shears are more than ten times exceed the geostrophic current shears. Turbulent vertical mixing coefficient, calculated by the G89 model (Gregg 1989) demonstrates the denominated minimum in the permanent pycnocline (4-10-6 m2/s) and increases with depth to achieve value 2-10-5 m2/s at 450 m depth. The profile of current velocity averaged by the ensemble casts in the near bottom boundary layer demonstrates near logarithmic relation.
Keywords: lowered acoustic Doppler current profiler LADCP, vertical structure of currents, vertical shear, vertical mixing, bottom boundary layer, Black Sea.
DOI: 10.22449/1573-160X-2017-1-25-40
© 2017, A.N. Morozov, E.M. Lemeshko, S.A. Shutov, V.V. Zima, D.V. Deryushkin © 2017, Physical Oceanography
Introduction
Acoustic Doppler current profilers (ADCP) have been actively applied in the field studies of the World Ocean dynamics for more than 30 years. These devices assist to solve a great variety of fundamental and applied problems of the modern oceanology. Profilers are produced by many manufacturers, such as E-Teledyne RD Instruments (USA), SonTek (USA), Nortek AS and Aanderaa Data Instruments AS (Norway). The 300 kHz Workhorse Monitor ADCP Profiler manufactured by TRDI (operating frequency - 300 kHz, nominal range - 120 m, resolution - 4 m) has been applied in the field research by Marine Hydrophysical Institute (MHI) since 2004. Overview of the device is shown in Fig. 1.
Application of ADCP in the Lowered ADCP (LADCP) mode is widespread in the study of the distribution of currents in the seas and oceans [1 - 3]. This type of measurements gives a unique opportunity to obtain the pattern of current velocity distribution throughout all the water column mass, repeatedly exceeding the operating range of the instrument [4]. The first LADCP-station was carried out in 1989 [5]. In subsequent years, this method has continued to be developed [6, 7]. Within the framework of the international project WOCE LADCP profiler became a regular oceanographic instrument [4].
During measurements, the ADCP, as a rule, is mounted in the CTD probe frame. In MHI the LADCP is a separate probe consisting of WHM300 profiler, which is mounted in a frame of titanium with a battery pack (Fig. 2, left). Methodical aspects of LADCP measurements in the Black Sea and the problem of data processing issues can be found in [8 - 10].
Fig. 1. WHM300 overview
46-
¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡lip
45
44
43
42
HU.J J i UUiW ........,«1111
■■■Ill
41
28
29
30
31 32 Longitude, °E
33
34
35
Fig. 2. Overview of LADCP (left) and station layout in three expeditions (right): Expedition 1 (green diamonds) - May, 8 - 18, 2004; Expedition 2 (blue triangles) - October, 6 -15, 2005; Expedition 3 (red circles) - July, 18 - 22, 2007
In preparation of this article the data collected in three expeditions was mainly used. These expeditions were carried out on the stations, the location of which is schematically shown in Fig. 2, right: Expedition 1 (green diamonds) - the stations performed in R/V Akademik cruise by the Bulgarian Academy of Sciences on May, 9 - 18, 2004 (GEF BSERP - RER/01/G33 project); Expedition 2 (blue triangles) -R/V Vladimir Parshin, October 6 - 15, 2005 GEF BSERP - RER/01/G33 project); Expedition 3 (red circles) - in R/V Experiment cruise on July, 18 - 22, 2007 (MHI Operational Oceanography project). Total from 2004 to 2014 LADCP had been applied in 17 expeditions to different parts of the Black Sea.
Vertical Structure of Current Velocity Field
Vertical structure of current velocity field in the Black Sea is studied according to the data of instrumental observations in different ways. Dynamic method permits to obtain continuous profiles of one component of current velocity
by CTD profiles at two stations and has the unsolved problem to select a level of the lack of movements [11 - 12]. The autonomous buoy stations (ABS) make it possible to obtain continuous time series of measurements at fixed horizons at the point location of the station [13 - 14]. Akvazond [15] and Akvalog [16 - 18] autonomous probing profilers allow to get regular synchronous profiles of a wide range of aquatic environment parameters (including three components of current velocity) with fine-resolution for several months at the station point location. ARGO Floating buoys provide current velocity estimation at fixed horizons, averaged over the track on the time interval between the buoy releases to the sea surface [19, 20]. As a result of probing, OLT hydrophysical complex permits to obtain continuous profiles of hydrological parameters and two components of current velocity [21, 22].
Fig. 3 vectorially shows the current velocity profiles (colored arrows) and the conditional density oe (gray lines) in the meridional section (31°E), crossing the anticyclonic eddy from north to south (Fig. 3, left). Data was collected in Expedition 1 (May 2004) [23, 24]. CTD data are supplemented by the current velocity profiles. They allow to assess the current state of water dynamics in the measurement area. In particular, according to the LADCP data of Expedition 1 kinematic characteristics of the anticyclonic eddy [9] and the spatial scales of the current velocity field variability in the surrounding area were determined. [24]
Fig. 3. Profiles of the current velocity (colored arrows) and the density (gray lines) in the meridional section, May, 2004
In the continental slope area a pronounced correlation of the vertical current velocity field structure with the vertical distribution of the conditional density is observed. Fig. 4, left, shows a graph of the dependence of the kinetic energy of the
currents (red line) (EK = (u2 + V2)/2) from the conditional density, isopycnally
averaged over the ensemble of the stations in Expedition 1 (EK = ( EK (o0)), where
(...) is an averaging operator); blue line marks the graph of available potential
0
500
energy EP = ( 0,52,2 N2) (where 2 is an isopycnal displacement from the middle
position); green line indicates the frequency graph of buoyancy
N =
' p dz
Fig. 4, right, shows the isopycnally averaged kinetic profiles (red line) and the available potential (blue line) energy, as well as the conditional density (green line). There is an almost linear dependence of the kinetic energy distribution of the conditional density (black dotted line) (EK(oe)) = 0,277 - 0,0164oe (Fig. 4, left) in the upper layer of 20 - 240 m (density of 13.6 - 16.6 kg/m3). In the lower layers the kinetic energy decreases more slowly, and at a depth of -400 m (oe = 16,9 kg/m3), its value goes to an almost stationary level.
Buoyancy frequency (N), rad/s 0 0.006 0.012 0.018 I_1_I_1_I_1_I_1
14 -
Ml M
CD b
15 -
£ J
<U 16 —
17 -1
0 -1
100 -
B 200 -£
<u
Q 300 -
400 -
Density (ctq), kg/m3 14 15 16 17 _J_1_I_1_I_1_I
1 1 I 1 T 0 0.02 0.04 0.06
Energy (EK, EP), J/kg
500 -1
I 1 I 1 I 1 I 0 0.02 0.04 0.06
Energy (EK, EP), J/kg
Fig. 4. Dependence of the kinetic (red line) and the available potential energy (blue line) from the conditional density (left) and depth (right). Green line indicates the buoyancy frequency (left) and the conditional density (right); black dotted line - the dependence of the kinetic energy distribution of the conditional density
Fig. 5 shows the mean current velocity profiles module (\u\ = 42ek )
obtained according to the data of LADCP in 2004 (red line) and autonomous buoy stations in 1955 - 1983 (blue line) [25]. Green solid line indicates the profile of the mean conditional density, dotted lines - profiles at the separate stations.
0 -|
200 -
ft'
<u
Q
400 -
600 -
800 -1
Density (ctq), kg/m3 14 15 16 17 _l_I_I_I_I_I_I
n I 1 I 1 I 0 10 20 30 Velocity (LADCP, Buoy ), cm/s
Fig. 5. Mean profiles of the current velocity according to the of data LADCP in May, 2004 (red line) and buoy stations in 1955 - 1983 (blue line); the conditional density in May, 2004 (green line); profiles on the separate stations (dotted lines)
The trends typical for the World Ocean dynamics are manifested in the vertical distribution of current velocity in the Black Sea. For example, the North Atlantic is characterized by a significant component of barotro-pic currents penetrating to greater depths and reaching an amount of 10 cm/s [26 - 28]. At that, the baroclinic inhomogeneity of the current velocity field is concentrated in the upper 1000 m layer [26]. In the Black Sea the layer of baroclinic/shear currents, according to the LADCP observations is limited by the upper 300 - 500 m layer, or depth of occurrence of isopycnal of oe = 16,9 kg/m3.
In the lower layers sea (750 - 1700 m) the barotropic component of current velocity, according to ARGO floating buoys, is 2 - 3 cm/s [19 - 20]. A similar estimation of the current velocity in the lower sea layer obtained based on the Akvalog profiler data [18].
Vertical shears of currents
Interest to the study of vertical shears of currents is mainly supported by the fact that the shear instability is represented by the principal mechanism of vertical mixing in the seas and oceans [29]. A large number of models/parametrizations
connects the turbulent mixing with the Richardson number
Ri =
N
2 A
Sh
PP81
parameterization [30] is widely applied in numerical simulation. It is similar to previously published summarized dependencies. In the Black Sea, in the words of the authors of the article [22, p. 14], "the main feature of the vertical structure of the currents is sharply-layered structure of the velocity field, at that, the layer boundaries coincide with the boundaries of the thermohaline structures and are located at depths of change of vertical density gradients. The maximum vertical shifts are registered near the core of the cold intermediate layer (CIL), where there is a break in the density gradient. "
Maximal velocity shears, according to [14], are located in layers of 10 - 25 m and 50 - 300 m. The authors of the paper [24] showed that in the Black Sea, the mean profile of the vertical shears has two well-defined maxima in the layers of the seasonal thermocline/pycnocline and main halocline/pycnocline. Maximal vertical shear in the density contrast layer separating the Azov and Black Sea waters is observed in the southern part of the Kerch Strait [31].
Fig. 6 (left) shows the profiles of the latitudinal (U) and meridional (V) current velocity and conditional density components obtained at one of the stations in 2004, where the depth resolution of the ADCP measurements was set to 2 m. In the seasonal thermocline layer, the maximum value of the buoyancy frequency reached 0.044 rad/s, the vertical shear of currents was 0.064 s-1. The corresponding Richardson number was 0.47, which exceeded its critical value of 0.25 [32]. Above the seasonal thermocline, at a depth of less than 10 m from the sea surface, the Ri values become less critical. This phenomenon was also observed in the coastal shelf waters of the Black Sea [15] and in the upper layer of the southern part of the Kerch Strait [33]. The maximum vertical shear of currents in the seasonal thermocline layer is quite often observed. Fig. 6 (right) demonstrates the mean profiles of vertical current shear (red line) and buoyancy frequency (green line). Averaging was performed relative to the depth of occurrence of the seasonal thermocline center. The Ri value at a depth of 15 m is ~1, in the seasonal thermocline —7 and below 35 m —9.
12.7
£
ft
<u
Q
Density (a ), kg/m3
13.1
13.5
13.9
14.3
-30 -20 -10 0
Current velocity (U, V ), cm/s St036 04:31 UTC 16.05.2004
10
Buoyancy frequency (N), rad/s 0 0.006 0.012 0.018 0.024
0 1 0 1
10 - \ i 10 -
V v -
20 - == •c -J. 20 -a . £ 30 CP 30 <u Q " 40 -
/ \
30 - ( V ft
40 - > » ' Í \ -
\ i V 50 -
50 - i s/ 1 60 -
- 2
--
—— ——— == s>-
é
\
\
0.03 J
0.002 0.004 0.006 Shear (Sh), s-1
0.008 0.01
0
Fig. 6. Left - the profiles of the latitudinal (U) and meridional (V) current velocity and conditional density components obtained at one of the stations in 2004. Right - the mean profiles of vertical current shear (red line) and buoyancy frequency (green line), averaged over the station ensemble
Fig. 7 (left) presents an example of the profiles of the components of the relative current velocity (reference horizon is 240 m) and the conditional density obtained at one of the stations in 2012 in the main pycnocline layer. The observed distribution of currents can be explained by the passage of an internal wave (IW) through the main halocline layer. Vertical length of the IW is 50 - 60 m, the amplitude is about 12 cm/s. The rotation of the current velocity vector with depth goes clockwise. For the northern hemisphere it means the downward propagation of the IW [34]. Intensive IW in the main halocline layer, according to LADCP measurements, was observed predominantly in the area of the continental slope during autumn, winter and spring.
14.8
Density (ct ), kg/m3 15.2 15.6 16 16.4
16.8
£
ft
<u
Q
80 -| 100 -120 -140 -160 -180 -200 -220 -240 -
1 )
■—
c.
----
V ——
>
c v
—.
V
\
\
<
\
Buoyancy frequency (N ), rad/s 0 0.006 0.012 0.018 I_,_I_,_I_,_I
I 1 I 1 I 1 I 1 r ■40 -30 -20 -10 0 Current velocity (U, V), cm/s St006 03:03 UTC 14.12.2012
1 10
100-
150-
a 200-
-a -
ct 0) 250-
U -
300-
350-
400-
I 1 I 1 I 0 0.002 0.004
Shear (Sh, Sh ), s-1
0.006
Fig. 7. Left - the profiles of the latitudinal (red line) and meridional (blue line) current velocity and conditional density components (green line) obtained at the separate station. Right - the mean profiles of the buoyancy frequency (green line), the geostrophic (blue line) and the measured (red line) current velocity shears
Fig. 7 (right) shows the mean profiles of the buoyancy frequency (green line), the geostrophic (blue line) and the measured (red line) current velocity shears. The maximum of vertical shears is observed in the main pycnocline layer. For vertical shears of geostrophic currents, Ri is -30 in the main pycnocline layer, then it increases with depth, reaching a value of -300 at the horizon of 400 m. According to the LADCP measured shears in the main pycnocline layer, Ri is -10, in the lower layers the Richardson number decreases to -4 at a depth of 400 m.
The significant difference between the profiles of geostrophic and measured shears is explained by the fact that the value of shears in the lower layers of the Black Sea is mainly determined by ageostrophic processes, such as, for example, internal waves (see Fig. 7, left). When substituting the obtained profiles of geostrophic and measured shears in the PP81 parameterization, a decrease with depth of the vertical turbulent mixing coefficient will be observed in the first case and in the second case - there will be an increase.
Vertical mixing
Vertical turbulent mixing in the seas and oceans has been remained one of the most interesting issues of applied oceanology for a long time. In recent decades, microstructured data collected using free-fall probes are used for the diapicnic exchange estimation. This data is considered to be the most reliable [2, 35 - 38]. Additionally, data of small-scale CTD/ADCP measurements [1, 2, 15, 36, 39 - 41] are also widely used.
Basis of the most models connecting small-scale characteristics of hydrophysical fields with vertical mixing parameters is the assumption that in a statically stationary field of internal waves turbulent kinetic energy is produced at
the velocity equal to the velocity of its transmission along the IW spectrum in the direction of decrease of vertical scales due to the interaction between waves [42]. One of these models, the G89 model [43], is often used to estimate the turbulent kinetic energy dissipation velocity [2, 36 and 40]:
N2 (Sh4)
e = 7-10-10 , /l0 , (1)
N (sh4GM76)
where e is the turbulent kinetic energy dissipation velocity, W/kg;_ N = ,
'p 8z
is
the buoyancy frequency; Sh is the vertical shear of current velocity (further - the
is the mean value of the fourth
shear), Sh = 7(AU/Az)2 + (AV/Az)2 , (Sh4)
degree of shear determined at 10-meter intervals; (ShGM76j is the mean value of the fourth power of the shear, obtained by integrating the spectrum of internal waves GM76 [44, 45]; N0 = 5.24-10"3 rad/s is the buoyancy frequency of GM76. The vertical diffusion coefficient Kp is determined from the following relation [46]
Kp = r—
p N2
(2)
where Kp is the vertical diffusion coefficient, m2/s; r is The the mixing efficiency
coefficient, usually taken as equal to 0.2 [47].
In the Black Sea, the profiles of shears and buoyancy frequencies have a well-pronounced maximum in the main halocline/pycnocline layer [22, 24]. Fig. 8 shows the profiles of the square of the current velocity shear (gray circles, red line) and the square of the buoyancy frequency (gray crosses, green line) averaged over the ensemble of 26 stations of the Expedition 1. The averaging was carried out with respect to the depth of occurrence of the pycnocline center with the subsequent shear of the resulting profile to the value average over the ensemble of stations. Shears and buoyancy frequency were determined at 10-meter intervals.
Fig. 9 (left) demonstrates the average profiles of the velocity of the kinetic turbulent energy dissipation (gray circles mark the initial data, red line - after low-frequency filtration), calculated from the relation (1). The present dependency graph shows the presence of a well-pronounced maximum in the main pycnocline layer, where e reaches 6-10-9 W/kg, in layers below 200 m e decreases with depth and at a horizon of 400 m its value is 9-10-10 W/kg.
Q 300 -
400 -
10-6
-TTTTT1 ........1
10-5 10-4 (Sh2), (N2), c-2
Fig. 8. Mean profiles of the square of the current velocity shear (gray circles, red line) and the square of the buoyancy frequency (gray crosses, green line)
10-10
T
10-' S, W/kg
100-
g 200 -
<u
Q 300 -
400 -
/ ; •
« *] •1
I -N
•• « • 4 X
M • rT i • •
• •«1 ¿ft • • • *• • •• • •
• » &
f **' * • 1 •«4 •
10-8
10-6
10-5 Kp, m2/s
10-4
Fig. 9. Left - the average profile of the velocity of the kinetic turbulent energy dissipation (gray circles mark the initial data, red line - after low-frequency filtration). Right - the profile of the vertical turbulent diffusion coefficient (gray circles are the initial data, the blue line - the data after filtration)
Fig. 9 (right) shows the profile of the vertical turbulent diffusion coefficient (gray circles are the initial data, the blue line - the data after filtration) calculated according to the relation (2). In the main pycnocline layer Kp takes the minimum
values of 4-10-6 m2/s, in layers below 200 m it increases with depth and at the horizon of 300 m reaches a value of 10-5 m2/s. The obtained estimates of the vertical mixing coefficient turned out to be somewhat smaller than in the previously published works by
S.G. Boguslavsky & I.K. Ivashchenko [48] - Kp = (3.1-4.4) • 10-4 m2/s at a depth
of500m;
V.N. Eremeev & V.M. Kushnir [22] - Kp = (3-5) • 10-5 m2/s in the main halocline;
A.S. Samodurov & A.M. Chukharev [49] -Kp = (2-8) • 10-5 m2/s in the main pycnocline;
A G. Zatsepin & co-authors [38] - Kp = 4 • 10-6-5 • 10-5 m2/s in the layer
below the seasonal thermocline to 170 m for summer and autumn.
The values of the vertical mixing coefficient for summer, given in [38], are in good agreement with the estimates obtained according to the data of CTD/LADCP observations.
A considerable dispersion in the estimates of the turbulent diffusion coefficient can be due to both the spatiotemporal variability of vertical mixing processes and the difference in the estimation methods. Taking into consideration the wide application of LADCP together with CTD, the approach can be used to study vertical mixing processes, especially in deep sea layers, where microstructural data is difficult to obtain.
Near-bottom boundary layer
Experimental study of the parameters of the near-bottom boundary layer in the seas and oceans is one of the topical problems of applied oceanography. The results of such studies can be useful for regulating the parameters of numerical models of water circulation and dynamic processes in the sea, in solving problems of soil erosion, sedimentation, etc. Empirical dependences of the near-bottom boundary layer parameters for mean ocean conditions are presently known [50]. In the Black Sea, the empirical characteristics of the near-bottom boundary layer were previously obtained using OLT - a modified hydrological probing complex [51, 52].
In July 2007 and 2009, the specific measurements were carried out in the area of the Black Sea experimental sub-satellite polygon. They were aimed to study dynamics of the near-bottom boundary layer. Fig. 10 schematically shows the location of stations in the expeditions during July 20 - 21, 2007 (squares) and July 29 - 30, 2009 (triangles). ADCP WHM300 was applied as a measuring instrument of current velocity profiles. The discreteness of the measurements over depth was set at 4 m, over time - 1 second, the bottom tracking option was enabled. To obtain current profiles in the near-bottom area, a 10-minute exposure of the instrument was carried out at a distance of 60 to 80 m from the bottom.
44.4 -|
44.35 -
£ 44.3-
3 44.25 -
44.2 -
44.15 -1
33.9 33.95 34 34.05 34.1 34.15 Longitude, °E
Fig. 10. Location map of the stations: July 20 - 21, 2007 (squares) and July 29 - 30, 2009 (triangles)
The results of the measurements showed that in a particular profile of the current velocity, the pronounced patterns of behavior in the near-bottom area are not revealed. This is largely due to the fact that the instrument registers the entire set of dynamic processes, including internal waves, inertial oscillations, etc. [53].
■ m4 l. A
J V B ^ ■ A / \C\ / ^
A L 1 A /
■ A I L ■ A / /
s ■ A| I L As /
A L Ï A / / r
■ A1 ■ L A
A L A, /
■ 1 1 / ■ 2( 307
— - — A2( D09
Fig. 11 presents an example of the profiles of the velocity and density components in the near-bottom layer, presumably during the passage (generation and reflection) of the internal wave.
Density (a0 ), kg/m3
14.4
14.8
15.2
15.6
£ ft
80 -
90 -
f 100 -
110 -
120 -
130 -1
t t tffi / U / V
\ t
1 \ {
\ y --- A
if t-- \
- N
\
-10
10
-5 0 5
Curent velocity (U, V), cm/s St041 13:32 UTC 16.05.2004
Fig. 11. Profiles of the velocity and density components in the near-bottom layer during the passage of the internal wave
In Fig. 12 the black continuous lines mark the current velocity profiles averaged over the ensemble of stations, depending on the distance to the bottom for the expeditions in 2007 and 2009. The direct dotted line is the empirical dependence of the current velocity at the upper boundary of the near-bottom boundary layer [52]. The average velocity profiles of currents can be adequately
approximated by the following logarithmic dependences [29]:
f - \
U = ^ln k
z
(3)
V y
where u* is the friction velocity; k is Karman's constant (0,41); zb is the height of bottom surface roughness; z is the distance to the bottom. The corresponding approximating functions are indicated by gray dotted lines in Fig. 12. The friction velocity was 3.14 cm/s and 2.46 cm/s for 2007 and 2009 respectively. The bottom roughness is the same for the two expeditions and is equal to 0.4 m. The value of the bottom friction coefficient is calculated according to the following ratio
Cb =(uj, (4)
where U is the velocity of current in the upper boundary of the near-bottom boundary layer, that equals -0.007, being more than 4 times higher than given in [51]. The friction velocity estimation, is calculated according to the following ratio [50]
u = 2.5 • H • f, (5)
PHYSICAL OCEANOGRAPHY NO. 1 (2017) 35
where H is the boundary layer thickness [52]; f is the Coriolis parameter that equals 1.5 cm/s for 2007 and 1.14 cm/s for 2009, which is two times less than the values obtained on the basis of ADCP-measurements. The observed discrepancy between near-bottom boundary layer parameter estimations can be explained both by the difference in the methods of instrumental estimation and by the increased dynamic activity in the area of the Black Sea experimental sub-satellite polygon [54].
Current velocity, cm/s
Fig. 12. The mean current velocity profiles in the near-bottom layer (black continuous lines); the empirical dependence of the current velocity at the upper boundary of the near-bottom boundary layer (direct dotted line)
Conclusion
A limited amount of the field material for brief summary of the ten-year LADCP experience is presented within the framework of this article. The results of the profiler application in the expeditionary practice of Marine Hydrophysical Institute broadened the concept of the vertical structure of the currents in the Black Sea. In particular, the characteristics of small-scale processes both in the water column and in the boundary layers were studied in greater detail. Summarizing, it may be concluded that the acoustic Doppler current profilers are undoubtedly the most powerful tool for studying the water dynamics and its variability over a wide range of spatio-temporal scales.
Acknowledgements. The research was carried out within the framework of Marine Hydrophysical Institute Project Operational Oceanography and with a partial financial support of the Ministry of Education and Science of the Russian Federation under the Federal Targeted Programme for Research and Development in Priority Areas of Development of the Russian Scientific and Technological Complex for 2014-2020 (Project ID No. RFMEFI57714X0110).
REFERENCES
1. Polzin, K., Kunze, E., Hummon, J. and Firing. E., 2002. The Finescale Response of Lowered ADCP Velocity Profiles. J. Atmos. Ocean. Techn., [e-journal] 19(2), pp. 205-224. https://doi.org/10.1175/1520-0426(2002)019<0205:TFROLA>2.0.CO;2
2. Fer, I., 2006. Scaling Turbulent Dissipation in an Arctic Fjord. Deep-Sea Res. II., [e-journal] 53(1-2), pp. 77-95. https://doi.org/10.1016Zj.dsr2.2006.01.003
3. Gaardsted, F., Zhou, M., Pavlov, V., Morozov, A. and Tande, K.S., 2010. Mesoscale Distribution and Advection of Overwintering Calanus Finmarchicus off the Shelf of Northern Norway. Deep-Sea Res. I, [e-journal] 57(11), pp. 1465-1473. doi:10.1016/j.dsr.2010.07.003
4. Firing, E., 1998. Lowered ADCP Developments and Use in WOCE. In: International WOCE Newsletter. Southampton, UK: WOCE International Project Office. No. 30, pp. 10-13. Available at: https://www.nodc.noaa.gov/woce/wdiu/wocedocs/newsltr/news30/news30.pdf [Accessed 10 May 2016].
5. Firing, E. and Gordon, R.L., 1990. Deep Ocean Acoustic Doppler Current Profiling. In: Proc. IEEE 4th Working Conf. on Current Measurements. Clinton, MD, USA: IEEE, pp. 192-201. doi: 10.1109/CURM.1990.110905
6. Fisher, J. and Visbeck, M., 1993. Deep Velocity Profiling with Self-Contained ADCPs. J. Atmos. Ocean. Tech.., [e-journal] 10(5), pp. 764-773. https://doi.org/10.1175/1520-0426(1993)010<0764:DVPWSC>2.0. CO;2
7. Visbeck, M., 2002. Deep Velocity Profiling Using Lowered Acoustic Doppler Current Profilers: Bottom Track and Inverse Solutions. J. Atmos. Ocean. Tech.., [e-journal] 19(5), pp. 794-807. https://doi.org/10.1175/1520-0426(2002)019<0794:DVPULA>2.0.CO;2
8. Morozov, A.N. and Lemeshko, E.M., 2005. Opyt Ispol'zovaniya Akusticheskogo Doplerovskogo Izmeritelya Techeniy (ADCP) v Usloviyakh Chernogo Morya [Experience in using of the Acoustic Doppler Current Profiler (ADCP) in the Black Sea Conditions]. In: Ekologicheskaya Bezopasnost' Pribrezhnoy i Shel'fovoy Zon i Kompleksnoe Ispol'zovanie Resursov Shel'fa [Ecological Safety of Coastal and Shelf Zones and Complex Use of Shelf Resources]. Sevastopol: MGI NANU. Iss. 12, pp. 457-476 (in Russian).
9. Morozov, A.N. and Lemeshko, E.M., 2005. Ispol'zovanie Self Contained ADCP dlya Provedeniya Izmereniy s Borta Sudna: Metodicheskie Voprosy i Fizicheskie Rezul'taty [Application Self Contained ADCP to Carry out Measurements from the Ship Board: Methodical Problems and Physical Results]. In: Ekologicheskaya Bezopasnost' Pribrezhnoy i Shel'fovoy Zon i Kompleksnoe Ispol'zovanie Resursov Shel'fa [Ecological Safety of Coastal and Shelf Zones and Complex Use of Shelf Resources]. Sevastopol: MGI NANU. Iss. 13, pp. 425-432 (in Russian).
10. Morozov, A.N. and Lemeshko, E.M., 2006. Methodical Aspects of the Application of Acoustic Doppler Current Profilers (ADCP) in the Black Sea. Physical Oceanography, [e-journal] 16(4), pp. 216-233. doi:10.1007/s11110-006-0027-8
11. Neumann, G., 1942. Die absolute Topographie des physikalischen Meeresniveaus und die Oberflächenströmungen des Schwarzen Meeres. Annalen der Hydrographie und Maritimen Meteorologie, 70(9), pp. 265-282 (in German).
12. Leonov, A.K., 1960. Regional'naya Okeanografiya. Chast' 1. Beringovo, Okhotskoe, Yaponskoe, Kaspiyskoe i Chernoe Morya [Regional Oceanography. Part 1. The Bering, Okhotsk, Japanese, Caspian and Black Seas]. Leningrad: Gidrometeoizdat, 480 p. (in Russian).
13. Boguslavsky, S.G., Bukatov, A.E. and Kazakov, S.I., 2001. Osobennosti Polya Skorosti i Vertikal'nogo Obmena v Chernom More [Pecularities of Velocity Field and Vertical Exchange in the Black Sea]. In: Ekologicheskaya Bezopasnost' Pribrezhnoy i Shel'fovoy Zon i Kompleksnoe Ispol'zovanie Resursov Shel'fa [Ecological Safety of Coastal and Shelf Zones and Complex Use of Shelf Resources]. Sevastopol: MGI NANU. Iss. 3, pp. 62-71 (in Russian).
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
38
Tuzhilkin, V.S., 2008. Thermohaline Structure of the Sea. In: The Black Sea Environment. The Handbook of Environmental Chemistry. Berlin: Springer Berlin Heidelberg. Vol. 5Q (Water Pollution), pp. 217-253. doi:10.1007/698_5_077
Ostrovskii, A.G., Zatsepin, A.G., Derevnin, V.A., Niozov, S.S., Poyarkov, S.G., Tsybulskiy, A.L. and Shvoev, D.A., 2008. Aquazond Moored Automatic Measuring System for vertical profiling of the marine medium. Oceanology, [e-journal] 48(2), pp. 275-283. doi: 10.1134/S000143700802015X
Ostrovskii, A.G., Zatsepin, A.G., Shvoev, D.A. and Soloviev, V.A., 2010. Underwater Anchored Profiler Aqualog for Ocean Environmental Monitoring. In: J.A. Daniels ed., 2010. Advances in Environmental Research. New York, USA: Nova Science Publishers, Inc. Vol. 4, pp. 201-218. Available at:
https://www.researchgate.net/publication/235426516_Underwater_Anchored_Profiler_Aqual og_for_Ocean_Environmental_Monitoring [Accessed 10 July 2016].
Ostrovskii, A. and Zatsepin, A., 2011. Short-Term Hydrophysical and Biological Variability over the Northeastern Black Sea Continental Slope as Inferred from Multiparametric Tethered Profiler Surveys. Ocean Dynamics, [e-journal] 61(6), pp. 797-806. doi:10.1007/s10236-011-0400-0
Ostrovskii, A.G., Zatsepin, A.G., Soloviev, V.A., Tsibulsky, A.L. and Shvoev, D.A., 2013. Autonomous system for vertical probing of the marine environment at a moored station. Oceanology, [e-journal] 53(2), pp. 233-242. doi:10.1134/S0001437013020124
Korotaev, G., Oguz, T. and Riser, S., 2006. Intermediate and Deep Currents of the Black Sea Obtained from Autonomous Profiling Floats. Deep-Sea Res. II, [e-journal] 53(17-19). pp. 1901-1910. https://doi.org/10.1016/j.dsr2.2006.04.017
Gerasimova, S.V. and Lemeshko, E.E., 2011. Otsenka Skorostey Glubokovodnykh Techeniy po Dannym ARGO [Estimation of the Velocities of Deep-Sea Currents according to ARGO Data]. In: V.N. Eremeev ed., 2011. Sistemy Kontrolya Okruzhayushchey Sredy [Monitoring System of Environment]. Sevastopol: MHI. Iss. 15, pp. 187-196 (in Russian).
Drozdov, A.E., Kushnir, V.M., Nikitin, A.V., Timoshenko, N.K. And Shevtsov A.T., 1991. Kompleks Gidrofizicheskiy Zondiruyushchiy dlya Okeanograficheskikh Issledovatel'skikh Sudov [Hydrophysical Probing Complex for Oceanographic Research Vessels]. Zapiski po gidrografii, (226), pp. 49-57 (in Russian).
Eremeev, V.N. and Kushnir, V.M., 1996. Sloistaya Struktura Techeniy i Vertikal'nyy Obmen v Chernom More [The Layered Structure of Currents and the Vertical Exchange in the Black Sea]. Okeanologiya, 36(1), pp. 13-19 (in Russian).
Kondrat'ev, S.I., Romanov, A.S. and Vnukov, Yu.L., 2007. Osobennosti Raspredeleniya Gidrokhimicheskikh Kharakteristik v Rayone Materikovogo Sklona Severo-Zapadnoy Chasti Chernogo Morya [Pecularities of Distribution of Hydrochemical Characteristics in the Region of the Continental Slope of the Northwestern Black Sea]. Morskoy Gidrofizicheskiy Zhurnal, (5), pp. 69-79 (in Russian).
Lemeshko, E.M., Morozov, A.N., Stanichnyi, S.V., Mee, L.D. and Shapiro, G.I., 2008. Vertical Structure of the Field of Current Velocities in the Northwest Part of the Black Sea Based on the LADCP Data for May 2004. Physical Oceanography, [e-journal] 18(6), pp. 319331. doi:10.1007/s11110-009-9029-7
Ivanov, V.A. and Belokopytov, V.N., 2011. Okeanografiya Chernogo Morya [The Black Sea Oceanography]. Sevastopol: ECOSI-Gidrofizika, 212 p. (in Russian).
Bubnov, V.A., 1996. The North Atlantic Current by Atlantex 90 Experiment Data. Oceanology, 34(6), pp. 733-737.
Cisewski, B., Budeus, G. and Krause, G., 2003 Absolute Transport Estimates of Total and Individual Water Masses in the Northern Greenland Sea Derived From Hydrographic and Acoustic Doppler Current Profiler Measurements. J. Geophys. Res., [e-journal] 108(C9), 3298. doi:10.1029/2002JC001530
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
Lherminier, P., Mercier, H., Gourcuff, C., Alvarez, M., Bacon, S. and Kermabon, C., 2007. Transports across the 2002 Greenland-Portugal Ovide Section and Comparison with 1997. J. Geophys. Res., [e-journal] 112(C7), C07003. doi:10.1029/2006JC003716
Bowden, K.F., 1983. Physical Oceanography of Coastal Waters. New York: Ellis Horwood. Ltd., 302 p. (Bowden, K., 1988. Fizicheskaya Okeanografya Pribrezhnykh Vod. Moscow: Mir, 324 p. (in Russian)).
Pacanowski, R.C. and Philander, S.G.H., 1981. Parametrisation of Vertical Mixing in Numerical Models of Tropical Oceans. J. Phys. Oceanogr., [e-journal] 11(11), pp. 1443-1451. https://doi.org/10.1175/1520-0485(1981)011<1443:P0VMIN>2.0.C0;2
Morozov, A.N. and Lemeshko, E.M., 2010. Otsenka Kontsentratsii Vzvesi po Dannym ADCP WHM1200 [Suspension concentration estimation according to ADCP WHM1200]. In: V.N. Eremeev ed., 2011. Sistemy Kontrolya Okruzhayushchey Sredy [Monitoring System of Environment]. Sevastopol: MHI. Iss. 14, pp. 42-46 (in Russian).
Miles, J.W., 1961. On the Stability of Heterogeneous Shear Flows. J. FluidMech., [e-journal] 10(4), pp. 496-508. https://doi.org/10.1017/S0022112061000305
Ivanov, V.A., Morozov, A.N., Kushnir, V.M., Shutov, S.A. and Zima, V.V., 2012. Raspredelenie Techeniy v Kerchenskom Prolive po Dannym ADCP-Nablyudeniy, Sentyabr' 2011 g. [Currents in the Kerch Strait, ADCP-observations, September 2011]. ]. In: Ekologicheskaya Bezopasnost' Pribrezhnoy i Shelfovoy Zon i Kompleksnoe Ispol'zovanie Resursov Shel'fa [Ecological Safety of Coastal and Shelf Zones and Complex Use of Shelf Resources]. Sevastopol: MGI NANU. Iss. 26, Vol. 1, pp. 170-178 (in Russian).
Leaman, K.D. and Sanford, T.B., 1975. Vertical energy propagation of internal waves: a vector spectral analysis of velocity profiles. J. Geophys. Res., [e-journal] 80(15), pp. 19751978, doi:10.1029/JC080i015p01975
Prandke, H. and Stips, A., 1998. Test Measurements with an Operational Micro structure-Turbulence Profiler: Detection Limit of Dissipation Rates. Aquat. Sci., [e-journal] 60(3), pp. 191-209. doi: 10.1007/s000270050036
Cisewski, B., Strass, V.H. and Prandke, H., 2005. Upper-ocean Vertical Mixing in the Antarctic Polar Front Zone. Deep-Sea Res. II, [e-journal] 52(9-10), pp. 1087-1108. https://doi.org/10.1016/j.dsr2.2005.01.010
Sundfjord, A., Fer, I., Kasajima, Y. and Svendsen, H., 2007. Observations of Turbulent Mixing and Hydrography in the Marginal Ice Zone of the Barents Sea. J. Geophys. Res., [e-journal] 112(C5), C05008. doi:10.1029/2006JC003524
Zatsepin, A.G., Golenko, N.N., Korzh, A.O., Kremenetskii, V.V., Paka, V.T., Poyarkov, S.G. and Stunzhas, P.A., 2007. Influence of the Dynamics of Currents on the Hydrophysical Structure of the Water and the Vertical Exchange in the Active Layer of the Black Sea. Oceanology, [e-journal] 47(3), pp. 301-312. doi:10.1134/S0001437007030022
Naveira Garabato, A.C., Oliver, K.I.C., Watson, A.J. and Messias M.-J., 2004. Turbulent Diapycnal Mixing in the Nordic Seas. J. Geophys. Res., [e-journal] 109(C12), C12010. doi: 10.1029/2004JC002411
Morozov, A.N. and Lemeshko, E.M., 2008. Vertikal'noe Peremeshivanie v Chernom More po Dannym CTD/LADCP Nablyudeniy [Vertical Mixing in the Black Sea according to CTD/LADCP Observations]. In: V.N. Eremeev ed., 2011. Sistemy Kontrolya Okruzhayushchey Sredy [Monitoring System of Environment]. Sevastopol: MHI, pp. 266-268 (in Russian).
Morozov, A.N. and Lemeshko, E.M., 2014. Otsenka Vertikal'noy Turbulentnoy Diffuzii po Dannym CTD/LADCP-Izmereniy v Severo-Zapadnoy Chasti Chernogo Morya v Mae 2004 Goda [Estimation of Vertical Turbulent Diffusion Coefficient by CTD/LADCP-Measurements in the Northwestern Part of the Black Sea in May, 2004]. Morskoy Gidrofizicheskiy Zhurnal, (1), pp. 58-67 (in Russian).
Gregg, M.C., Sanford, T.B. and Winkel, D.P., 2003. Reduced Mixing from the Breaking of Internal Waves in Equatorial Waters. Nature, [e-journal] 422, pp. 513-515. doi:10.1038/nature01507
43. Gregg, M.C., 1989. Scaling turbulent dissipation in the thermocline. J. Geophys. Res., [e-journal] 94(C7), pp. 9686-9698. doi:10.1029/JC094iC07p09686
44. Garrett, C. and Munk, W., 1975. Space-Time Scales of Internal Waves: A Progress Report. J. Geophys. Res., [e-journal] 80(3), pp. 291-297. doi:10.1029/JC080i003p00291
45. Cairns, J.L. and Williams, G.O., 1976. Internal Wave Observations from a Midwater Float, 2. J. Geophys. Res., [e-journal] 81(12), pp. 1943-1950. doi:10.1029/JC081i012p01943
46. Osborn, T.R., 1980. Estimates of the Local Rate of Vertical Diffusion from Dissipation Measurements. J. Phys. Oceanogr, [e-journal] 10(1), pp. 83-89, doi:10.1175/1520-0485(1980)010<0083:EOTLRO>2.0.CO;2
47. Moum, J.N., 1996. Efficiency of Mixing in the Main Thermocline. J. Geophys. Res., [e-journal] 101(C5), pp. 12057-12069. doi:10.1029/96JC00508
48. Boguslavsky, S.G., Ivaschenko, I.K., 1990. Vertical Mesostructure of Deep Waters in the Black Sea. Soviet Journal of Physical Oceanography, [e-journal] 1(6), pp. 433-442. doi: 10.1007/BF02196997
49. Samodurov, A.S. and Chukharev, A.M., 2006. Otsenka Intensivnosti Vertikal'nogo Turbulentnogo Obmena v Chernom More po Eksperimental'nym Dannym [Estimation to intensity of vertical turbulent exchange in Upper Layer of the Black Sea according to In Situ Measurements]. In: Ekologicheskaya Bezopasnost' Pribrezhnoy i Shel'fovoy Zon i Kompleksnoe Ispol'zovanie Resursov Shel'fa [Ecological Safety of Coastal and Shelf Zones and Complex Use of Shelf Resources]. Sevastopol: MGI NANU. Iss. 14, pp. 524-529.
50. Weatherly, G.L. and Martin, P.J., 1978. On the Structure and Dynamics of the Oceanic Bottom Boundary Layer. J. Phys. Oceanogr., [e-journal] 8(4), pp. 557-570, doi:10.1175/1520-0485(1978)008<0557:OTSADO>2.0.CO;2
51. Kushnir, V.M., 1999. Turbulent Diffusion in the Near-Bottom Boundary Layer of the Black Sea Shelf Zone. J. Mar. System, [e-journal] 21(1-4), pp. 243-253. https://doi.org/10.1016/S0924-7963(99)00017-2
52. Kushnir, V.M., 2007. Bottom Boundary Layer in the Black Sea: Experimental Data, Turbulent Diffusion, and Fluxes. Oceanology, [e-journal] 47(1), pp. 33-41. doi: 10.1134/S0001437007010067
53. Morozov, A.N., Lemeshko, E.M., Shutov, S.A., and Zima, V.V., 2012. Techeniya v Sevastopol'skoy Bukhte po Dannym ADCP-Nablyudeniy, Iyun' 2008 g. [Currents in the Sevastopol Bay from ADCP-Observations (June 2008)]. Morskoy gidrofizicheskiy zhurnal, (3), pp. 31-43 (in Russian).
54. Morozov, A.N., Lemeshko, E.M., 2009, "Pridonnyy pogranichnyy sloy na shel'fe Yuzhnogo berega Kryma po dannym nablyudeniy [Bottom boundary layer on the shelf of the Southern coast of Crimea according to the observations]", In: V.N. Eremeev ed., 2011. Sistemy Kontrolya Okruzhayushchey Sredy / Sredstva, Informatsionnye Tekhnologii i Monitoring [Monitoring System of Environment / Means, Information Technologies and Monitoring]. Sevastopol: MHI, pp. 270-272 (in Russian).