Since 2010, the P.P. Shirshov Institute of Oceanology RAS (SIO RAS) in Gelendzhik maintains the research (observational) site for year round multi-disciplinary studies and monitoring of the marine environment in the coastal zone. Analysis of the data obtained at the observational site revealed the existence of well pronounced short-term variability of coastal zone hydrodynamics at time scales from 1-3 days to 1-2 weeks. The paper examines the role of external forcing (including the impact of adjoined open sea dynamics and wind stress) in the short-term variability of hydrodynamics and upper mixed layer evolution.
Black Sea, coastal zone, short-term variability, hydrodynamics, wind forcing.
Studies of coastal zone hydrodynamics in the semi-enclosed and marignal seas are primarily motivated by conventional wisdom that the coastal currents is main driver for shelf water ventilation, self-cleaning and pollution-removing processes. Coastal water pollution is often serious problem for major ports, coastal cities and recreational facilities. In the Russian part of the Black Sea, one of the regions under strong anthropogenic stress is the shelf zone in vicinity of cities Gelendzhik and Novorossiysk. Since 2010, the Southern Branch of the P.P. Shirshov Institute of Oceanology RAS (SIO RAS) in Gelendzhik maintains the observational site for multi-disciplinary studies and monitoring of the marine environment in the coastal zone.
The site occupies an area of 10*15 km2 (Fig. 1, ). The continuous year-round monitoring is carried out using autonomous measurement systems both bottom mounted and moored buoys. Three types of autonomous platforms are used: 1) acoustic Doppler velocity profiler (ADCP) at the bottom mount; 2) thermo-chain at the mooring line; 3) moored profiler Aqualog . By these platforms the time series of hydrophysical and bio-optical data (vertical profiles of temperature, salinity, density, current velocity, acoustic backscatter, water transmissivity, chlorophyll-a fluorescence, etc.) of high spatial and temporal resolution are obtained. The data from some of the autonomous platforms is transferred via a communication line to the coastal center for real-time operability. Also the ship-born multidisciplinary surveys are carried out on a regular basis.
The observational data are used for research on marine environment and biota, cross-shelf exchange processes and etc. In particular, it has been shown that there is a significant short-term variability of the coastal zone hydrodynamics and upper mixed layer thickness at the scale from 1-3 days to 1-2 weeks [3-5]. It was found that this variability is to a certain extent due to the influence of external forcing including the open ocean dynamics and the local wind forcing.
Fig. 1. Scheme of the observational site at the Black Sea coastal zone near Gelendzhik. Triangles are coastal meteorological stations, black stars – joint stations with bottom mounted ADCP and moored thermo-chains, small black circles – CTD and multidisciplinary stations of ship-borne surveys, large black circle is moored automatic profiler Aqualog. Blue star is a lost station. Dotted lines with arrows show the R/V tracks for measurements of the current velocity using towed ADCP. Red stars – autonomous stations to be deployed in 2016. Red arrows - HF Doppler radar. Green diamond – proposed location of the moored surface meteorological buoy.
- Influence of open ocean dynamics
A deep part of the Black Sea is characterized by cyclonic circulation, encircled by the Rim Current (RC) . The basin scale cyclonic gyre and the RC are formed under the impact of cyclonic vorticity of the wind over the Black Sea region . In autumn and winter, when wind cyclonic vorticity is enhanced, RC is strong and pressed against the continental slope. In spring and summer when cyclonic vorticity of the wind decreases, RC becomes weak; owing to baroclinic instability, its flow forms meanders and often breaks into mesoscale eddies [8-10]. In autumn and winter, the growth of the RC meanders to large amplitude is restrained by the Ekman pumping . However, the RC meandering can be developed in any season of the year.
Figure 2 shows the time-depth diagram of specific density and alongshore current velocity measured by the profiler Aqualog on moored buoy station at a depth of 270 m (the upper part of the continental slope) from November 5 till December 27, 2012). The alongshore velocity changes its direction in 5-15 days. In its turn the isopycnals change their position coherently with the velocity. They deepen with the development of north-western current, and rise with the development of the south-eastern current. The amplitude of vertical displacement of isopycnals including the contours of specific density 16.2 kg/m3, which coincides with the upper boundary of hydrogen sulfide layer, is up to 50 m.
Fig.2. Variability of vertical distributions of the specific density and alongshore current velocity in the upper 200 m from November 5 to December 27, 2012 according to profiler Aqualog data. Profiling frequency was 4 times per day. Black lines are contours of specific density (isopicnals), orange line is isopycnal 15.8 kg/m3, lower limit of the oxygen-containing layer, red line is isopycnal 16.2 kg/m3, the upper boundary of hydrogen sulfide layer. Color marks the direction of longshore velocity (color bar at the right vertical axis), its negative values are for the SE current direction, positive are for the NW current direction.
The nature of the described above quasi-periodic oscillations is the RC meandering (Fig. 3). Due to displacement of the RC meanders along the shore in the NW direction, moored station with Aqualog profiler periodically appears in the area of cyclonic and anticyclonic meanders. Being in cyclonic meander profiler detects the NW flow and lowering of isopycnals due to geostrophic adjustment. After several days or a week, the moored station is in the anticyclonic meander and the Aqualog profiler measures the SE flow and the rise of isopycnals .
Quasi-periodic oscillations both of direction and magnitude of the alongshore current velocity at the upper part of the continental slope have strong impact on the hydrodynamics of the shelf. Analysis of satellite and towed ADCP velocity survey data, revealed that in case of a cyclonic RC meander anticyclonic submesoscale eddies are formed at the shelf, while in case of the anticyclonic meander the cyclonic submesoscale eddies are generated [3, 4].
Fig. 3. Scheme of quasi-periodic (5-15 days) oscillations of the alongshore current velocity. The black arrows show position of the RC axis. The mark C denotes a cyclonic meander, the mark A - anticyclonic. Thick red arrow indicates the direction of meanders alongshore displacement. Red cross indicates the position of the Aqualog profiler mooring.
The typical diameter of the submesoscale eddy is 2-10 km, orbital velocityis 20-50 cm/s, alongshore displacement (translation) velocity is 10-20 cm/s. The lifetime of submesoscale eddy does not exceed a few days . Unlike mesoscale eddy, which is characterized by small Rossby number, submesoscale eddy is characterized by the Rossby number close to 1 or even more thereby they are ageostrophic eddies.
Submesoscale anticyclones often accumulate passive tracers and ensure rapid transport of pollutants from the coast to the outer shelf (Fig. 4), whereas submesoscale cyclones produce not only horizontal, but also intense vertical circulation and contribute to the rise of the rich chlorophyll thermocline waters to the sea surface (Fig. 5, ). These eddies play an important role in the ventilation of the coastal zone especially in areas with a concave shoreline where circulation is weaker in comparison with that in the shelf regions with a flat or convex shore line.
Physical mechanisms of submesoscale eddies formation were studied by means of laboratory modelling [12-14]. Several experiments were conducted in the cylindrical tank on top of rotating platform filled by homogeneous or stratified fluid. First, the mechanism of eddies formation behind the capes was studied . It was shown that periodic eddy formation behind the cape and separation of eddies from the coast occur only in the case of decelerating cyclonic alongshore current in rotating fluid. The explanation for the nature of this effect based on the account of bottom Ekman friction was suggested. The dependences of eddy basic characteristics (diameter, orbital velocity, etc.) on time and external parameters were established.
The submesoscale eddy formation process due to horizontal velocity shear instability in the rotating fluid was investigated in cases of cyclonic and anticyclonic shear between the alongshore current axis and round tank wall. It was shown that in presence of a current with cyclonic shear, chain of eddies was formed in wide range of the non-dimensional shear value. In case of a current with the anticyclonic shear, the chain of eddies was observed only at rather small values of the shear value and by contrast, the current had chaotic turbulent structure at the large non-dimensional shear value. The asymmetry of the eddy-like structure formation in the current with cyclonic and anticyclonic shear was explained by simple physical model .
Results of the laboratory experiments were used successfully to interpret the observations of the submesoscale eddies at the narrow shelf of the NE Black Sea.
As a characteristic period of the oscillations of the alongshore current direction due to meandering of the RC is about 10 days, during one half of this period the conditions for formation of submesoscale cyclones are favorable, and during other half of the period they are favorable for submesoscale anticyclones formation.
It should be noted that in the periods of restructuring of the circulation upon the continental slope (transition from a cyclonic RC meander to anticyclonic one or vise verse) shelf circulation could be rather weak. In such situations, contaminants may accumulate at the source for several days, resulting in their concentration grouth.
- coastal upwelling and downwelling and their relation to the wind forcing and the sea currents
According to Ekman's theory, the events of the coastal upwelling and downwelling should occur in NE Black Sea if the wind has significant alongshore component: NW or SE, respectively. It is known that Ekman coastal upwelling plays important role in maintaining a high level of bio-productivity.
In the NE part of the Black Sea, strong and continuous northwesterly winds are rare. Therefore complete Ekman upwelling events occur no often than a few times during the warm season of the year . Here, according to , we mean that upwelling is complete if the thermocline water rises to the sea surface. When thermocline rises, but does not reach the sea surface, the upwelling is incomplete. Unlike complete upwelling, incomplete upwelling is observed quite often in the coastal area of the NE Black Sea: several times per month.
One of the complete Ekman upwelling events, which was followed by well developed downwelling, was observed at the SIO RAS research site on September 24-28, 2013.
Figure 6 illustrates the interrelations between wind forcing, water temperature and alongshore current velocity during the period of the upwelling and downwelling events mentioned above. The time based diagrams show wind velocity vectors measured by the automatic weather station (Fig. 6a) and sea temperature vertical distribution measured by the moored thermochain, as well as absolute value and direction of the current velocity vertical distributions measured by bottom mount ADCP (Fig. 6b–6d, respectively). During September 24-26, prevailing northwesterly wind led to the temperature decrease by more than 10°C in the entire water column. The complete upwelling event was accompanied by development of the coastal current of southeastern direction. On September 26 (second half of the day), the wind direction changed to the southerly then to the southeasterly and the wind velocity increased. The wind change resulted in relaxation of the upwelling and quick development of the downwelling manifested itself in temperature increase up to 19.5°C in the entire water column. The downwelling phase was accompanied by the development of the coastal current of the NW direction, which velocity peaked at 70 cm/s.
A criterion of the complete Ekman upwelling in two-layer fluid was obtained in :
txt/frwHRd > 1 (1)
Here tx is the wind stress alongshore component , t is time from the beginning of the upwelling favorable wind, f is the Coriolis parameter, and rw is the seawater density, H - initial upper mixed layer (UML) depth, Rd = (DrgH/rw)0.5/f is local baroclinic deformation radius and Dr is the density jump between the upper and lower layers. This criterion should be verified.
Fig. 6. Observational data of September 18–27, 2013: (a) wind velocity measured by automatic weather station; (b) vertical distribution of water temperature based on thermistor chain data; (c) absolute value of the current velocity measured by the bottom-mount ADCP; (d) direction of the current measured by on the bottom-mount ADCP.
Figure 6b demonstrates also an incomplete upwelling event on September, 18-22. Cold thermocline water did not rich the surface instead it was observed at the depth of about 10 m. The cause of this upwelling event was not a wind forcing. As could be seen from Fig. 6a, wind was weak and variable during those days. However, the existence of the alternating alongshore and eddy currents, whose nature is not directly related to local wind forcing, should be a widespread cause of coastal incomplete upwellings and downwellings. Indeed, if a NW baroclinic alongshore jet exists in the upper due to geostrophic adjustment a downwelling should be observed near the coast and the upwelling would occur in the case of a SE alongshore current.
This upwelling of dynamical origin could enhance the wind induced turbulent entrainment of thermocline water into UML. Such events could increase the primary production and chlorophyll-a concentration of UML because the thermocline water is rather rich in nutrients. An example of the impact of wind-induced turbulent entrainment on the chlorophyll-a concentration is given in the next section of this paper.
Concluding this section we note that the analysis of variations of sea temperature and current velocity along with wind forcing allows us to understand the origin of the observed coastal upwelling and downwelling. This approach will be used in future investigations of upwelling and downwelling processes in the Black Sea.
IV. Wind-induced turbulent entrainment of thermocline water into upper mixed layer as a source of maintaining primary production
Wind-induced turbulent entrainment is an important process that impacts the thickness and temperature of UML in a warm period of a year. The velocity We of turbulent entrainment is strongly nonlinear function of wind forcing. According to the parametrization 
We = 2.5U*3/(gHΔρ/ρw), (2).
Here U* = (t/ρw)0.5 = (Сdρa/ρw)0.5V is dynamical frictional velocity in water, t is wind stress, Сd is wind drag coefficient, ρa is air density, V is absolute value of wind velocity at 10 m height, Δρ – density jump across the thermocline.
Fig.7. Two sequential thermal images (NOAA-15) of the NE Black Sea in July 2006 (GMT, hh:min): (a) July 01, 19:03; (b) July 03, 14:32. The thick arrow denotes the warm alongshore current.
It follows from (2) that the rate of turbulent entrainment is proportional to the cube of wind velocity and inversely proportional to the thickness of UML. If everything else being equal, an intensive turbulent entrainment occurs when the wind is strong, and the initial thickness of the UML is small. This situation was observed at the beginning of July in 2006 when after a long period of weak wind, strong burst of the northeasterlies of the speed about 15 m/s was observed. Initially thickness of the UML was about 10 m. After two days of the strong wind, the UML thickness increased to 25 meters while its temperature dropped from 22-25°C to 15-18°C (Fig.7), see also . The application of (2) to evaluate the temporal evolution of the UML under the real wind forcing yielded results similar to those observed (see Fig. 8).
Fig.8. Time evolution of absolute value of the wind velocity V based on the NCEP data, the UML thickness H and its temperature T during July 1–7, 2006, in the area of strong winds between the Kerch straight and city of Tuapse.
One important phenomenon related to strong northeasterly wind need to be mentioned also is a formation of warm alongshore cyclonic current (Fig. 7). This fast current originates from Tuapse-Sochi region just after the wind decay and reaches Gelendzhik – Novorossiyisk region in 2-3 days. The physical origin and characteristics of the current are described in . Here we would like to notice that quasi-periodic formation of such currents after the northeasterly wind events could be an efficient mechanism of alongshore transport of nutrients and pelagic species from the SE sea towards the NW sea.
Fig. 9. Chlorophyll-a concentration at two sequential satellite images (MODIS-Aqua) in the NE Black Sea in July 2006 (GMT, hh:min): (a) July 01, 19:03; (b) July 05, 14:32.
An intense turbulent entrainment of thermocline water into UML led to increase of satellite measured chlorophyll-a concentration (CC) two days after the wind decay (Fig. 9). Using the empirical relation between the CC values derived from the satellite data and in situ measurements was calculated. The plot of these values is shown in Fig. 10 along with that of the wind velocity. As follows from Fig. 10 the increase in CC lasted for a few days. In a week after the decay of strong wind, CC recovered to its initial level. An impact of wind forcing on CC and primary production should be studied in more details in the future.
Fig. 10. Time evolution of chlorophyll “a” concentration (Cl_a) together with wind velocity (V).
The development of a moored profiler was supported by the Program №43 of the RAS, observational part of the work - by the RNF grant №14-50-00095, the data processing and analysis – by the RFBR grants №№14-05-00159 and 14-05-00792.
1. Zatsepin, A. G., Ostrovskii, A. G., Kremenetskiy, V. V., Nizov, S. S., Piotukh, V. B., Soloviev, V. A., et al. Subsatellite polygon for studying hydrophysical processes in the Black Sea shelf-slope zone // Izvestiya RAS, Atmospheric and Oceanic Physics. 2014. V.50 №1. P. 13-25.
2. Ostrovskii A. G., Zatsepin A. G., Solovyev V. A., Tsibul’skii A. L., and Shvoev D. A. Autonomous system for vertical profiling of the marine environment at a moored station // Oceanology. 2013. V. 53 № 2. P. 233–242.
3. Zatsepin A. G., Korzh A. O., Kremenetskii V. V., Ostrovskii A. G., Poyarkov S. G., and Solov’ev D. M. Studies of the hydrophysical processes over the shelf and upper part of the continental slope of the Black Sea with the use of traditional and new observation techniques // Oceanology. 2008. V.48. № 4. P. 446–475.
4. Zatsepin A. G., Kondrashov A. A., Korzh A. O., Kremenetskii V. V., Ostrovskii A. G., and Solov’ev D. M. Submesoscale eddies at the Caucasus Black Sea shelf and the mechanisms of their generation // Oceanology. 2011. V. 51. №4. P. 554–567.
5. Zatsepin A. G., Piotukh V. B., Korzh A. O., Kukleva O. N., and Solov’ev D. M. Variability of currents in the coastal zone of the Black Sea from long-term measurements with a bottom mounted ADCP // Oceanology. 2012. V. 52. №5. P. 579–592.
6. Bogatko, O. N., S. G. Boguslavskii, Yu. M. Belyakov, and R. I. Ivanov. Surface currents in the Black Sea – In: Kompleksnye Issledovaniya Chernogo Morya, Marine Hydrophysical Institute. 1979. Sevastopole. Ukraine. P. 25–33. (in Russian).
7. Stanev E.V. On the mechanisms of the Black Sea circulation. Earth-Sci. Rev. 1990 V. 28. P. 285–319.
8. Zatsepin A. G., Denisov E. S., Emel'yanov S. V., Kremenetskiy V. V., Poyarkov S. G., Stroganov O. Yu., et al. // Effect of bottom slope and wind on the near-shore current in a rotating stratified fluid: laboratory modeling for the Black Sea // Oceanology. 2005. V. 45. №1, P. S13-S26.
9. Ginzburg, A. I., A. G. Kostianoy, D. M. Soloviev, and S. V. Stanichny. Remotely sensed coastal/deep-basin water exchange processes in the Black Sea surface layer – In: Satellites, Oceanography and Society, Ed. D. Halpern, Elsevier Oceanography Series, 63, 273–287, 2000.
10. Zatsepin, A. G., Ginzburg, A. I., Kostianoy, A. G., Kremenetskiy, V. V., Krivosheya, V. G., Stanichny, S. V., & Poulain, P. M. Observations of Black Sea mesoscale eddies and associated horizontal mixing // Journal of Geophysical Research: Oceans. 2003., V. 108. №. C8. P. 1978–2012.
11. Zatsepin A.G., Ostrovskii A.G., KIremenetskiy V.V., Piotoukh V.B., Kuklev S.B., Kukleva O.N., et al. On the nature of short-period oscillations of the main Black Sea pycnocline, submesoscale eddies and response of the marine environment to the catastrophic shower of 2012 // Izvestiya RAS, Atmospheric and Oceanic Physics. 2013. V. 50 №1. P. 717-732.
12. Elkin D.N., Zatsepin A.G. Laboratory Investigation of the Mechanism of the Periodic Eddy Formation Behind Capes in a Coastal Sea // Oceanology. 2013. V. 53. No.1 .P. 24-35.
13. Elkin D.N., Zatsepin A.G. Laboratory study of a shear instability of an alongshore sea current // Oceanology. 2014. V.54. No.5. P. 576-582.
14. Elkin D.N., Zatsepin A.G. Laboratory studies of the of eddy formation in rotating and non-rotating fluid due to spatially non-uniform wind forcing // Selected papers of EMECS11 – Sea Coasts XXVI Joint Conference. 2016 (this volume).
15. Zatsepin, A. G., Silvestrova, K. P., Kuklev, S. B., Piotoukh V. B., and Podymov O. I. Observations of a Cycle of Intense Coastal Upwelling and Downwelling at the Research Site of the Shirshov Institute of Oceanology in the Black Sea // Oceanology. 2016. V.56. No.2. P. 203-214.
16. Zhurbas, V. M. Oh, I. S. and Park, T. Formation and decay of a longshore baroclinic jet associated with transient coastal upwelling and downwelling: a numerical study with applications to the Baltic Sea // J. Geophys. Res. 2006. V. 111. P. 1--18 doi: 10.1029/2005JC003079.
17. Kato, H. and Phillips, O. M. On the Penetration of a Turbulent Layer into a Stratified Fluid // J. Fluid Mech. 1969. V. 37. P. 643–655.
18. Zatsepin A.G., Kremenetskiy V.V., Piotoukh V.B., Poyarkov S.G., Ratner Yu.B., Soloviev D.M., et al. Formation of the Coastal Current in the Black Sea Caused by Spatially Inhomogeneous Wind Forcing upon the Upper Quasi-Homogeneous Layer // Oceanology. 2008. V.48. No.2. P. 176-192.