\documentclass[epub,pdf]{ebook}
\usepackage[cm]{sfmath}
\usepackage{amsmath}
\usepackage{color}
\usepackage[pdftoolbar=false, pdfpagemode=UseNone, linkbordercolor={0 0 0}]{hyperref}

\ident{ES}{19}{4}{2019}{2019ES000672}
\cpright{}{2019}
\pagestyle{empty}
\definecolor{brown}{RGB}{157,23,4}

\begin{document}

\vspace*{2\baselineskip}

\title{Hindcast of the mesoscale eddy field in the Southeastern Baltic Sea: Model output vs satellite imagery}

\author{Victor Zhurbas$^{1}$, Germo V\"ali$^{2}$, Andrey Kostianoy$^{1,3}$, Olga Lavrova$^{4}$}

\affil{1}{Shirshov Institute of Oceanology RAS, Moscow,  Russia}
\affil{2}{Tallinn University of Technology, Department of Marine Systems, Tallinn,  Estonia}
\affil{3}{S. Yu. Witte Moscow University, Moscow,  Russia}
\affil{4}{Space Research Institute RAS, Moscow,  Russia}

\abstract{An ultra-high resolution circulation model (0.125 nautical mile grid) of the Southeastern Baltic Sea is compiled from the General Estuarine Transport Model (GETM) in order to simulate the mesoscale eddy field in the area. The model results are compared with optical and infrared satellite imagery available for the period of May--August 2015. Mesoscale eddies detected in the satellite images are reasonably well identified in the simulated patterns of the sea surface temperature, currents, and floating Lagrangian \linebreak\newpage}\noindent\centerline{\parbox{100mm}{  particles.}}


\section{1. Introduction}


\begin{sloppypar}
The baroclinic Rossby radius of deformation, $R_{BC}$, is the basic length scale for the size of mesoscale eddies, which determine the weather in both the atmosphere and ocean. A typical value of $R_{BC}$ for the troposphere is $R_{BC} =1000$~km, and the number of sites of quasi-simultaneous (e.g., daily) measurements of meteorological parameters in a $R_{BC} \times R_{BC}$ square, $n$ is sufficiently large $(n \gg 1)$ to construct a weather map. According to recent data, there is one land-based meteorological station for every 1500~km$^2$ on the average
[\itc{Kostianoy,} \reflink{Kost17}{2017}],
i.e., $n=667$. The weather maps/data from yesterday and today can be extrapolated for tomorrow (or, more sophistically, assimilated in atmospheric models), which is a straightforward approach to the weather forecast. A typical value of the baroclinic Rossby radius in the ocean is much smaller than in the atmosphere, namely, $R_{BC}=30$~km in the middle latitudes, so that basically the number of oceanographic measurements in a $R_{BC} \times R_{BC}$  square is much smaller than one, and thousand times more oceanographic measurement sites per unit square are needed to resolve the ocean weather map to the same degree of detail as the weather map in the atmosphere.  Note that the most advanced ocean-observing program Argo consisting currently of 3988 (as of 24 November  2018) free-drifting profiling floats in the World Ocean (see
\href{http://www.argo.ucsd.edu/}{http://www.argo.ucsd.edu/})
provides only one float in a $94 \times R_{BC} \times R_{BC}$ square or $n \approx 0.011$, which is far from being enough to construct the ocean weather map. Nevertheless, further progress towards the ocean weather forecast is believed to be achieved by using numerical modeling and remote sensing technologies, such as, e.g., Argo drifters and satellite altimetry
[\itc{Vignudelli et al.,} \reflink{Vign11}{2011}].
\end{sloppypar}

However, the mesoscale ocean dynamics forecast (i.e., the weather) in relatively shallow inland/coastal seas, where the circulation strongly depends on the wind direction and strength, as well as on the coastline configuration and bottom topography, such as the Baltic Sea, is less challenging than in the open ocean because the process of forming of mesoscale eddies and fronts here becomes predictable to a certain extent. Thus, for example, it is known that in the case of an alongshore wind favorable for coastal upwelling frequently observed in the Baltic Sea
[\itc{Lehmann and Myrberg,} \reflink{Lehm08}{2008};
\itc{Myrberg and Andrejev,} \reflink{Myrb03}{2003}],
cyclonic eddies are formed behind the topographic features of the coastline (headlands or capes) and bottom elevations; these eddies can eventually break away from the formation sites and travel into the open sea
(e.g.
[\itc{Elkin and Zatsepin,} \reflink{Elki13}{2013};
\itc{Ginzburg et al.,} \reflink{Ginz98}{1998};
\itc{Zhurbas et al.,} \reflink{Zhur17}{2017}]).
Moreover, not only the place of probable formation of a cyclonic vortex is known, but also the time itself -- after the culmination of upwelling, when the wind begins to fade
[\itc{V\"{a}li et al.,} \reflink{Vali17}{2017}].
Similarly, in the case of an alongshore wind favorable for coastal downwelling, anticyclonic eddies are formed behind the headlands. In the case of the Southeastern Baltic Sea, prominent features of the coastline behind which the mesoscale vortices can form are the Sambian Rise, Cape Taran, Hel Spit/Peninsula, and Cape Rozewie
(see
 \figref{1}).           %     Figure 1).
%%
\def\topfraction{.9}
\def\textfraction{.1}
%%
\begin{figure}[t]                      % Fig 1
\figurewidth{90mm}
\setimage{}{}{115mm}{}{2019es000672-f01}
\fcaplabel{Left panel: Map of the high resolution model domain (filled colors) with the open boundary locations (black lines). The coarse resolution model domain (blank contours + filled colors) has an open boundary close to the Gothenburg station. Right panel: Close-up of the Southeastern Baltic Sea. 1 -- Sambian Rise, 2 -- Cape Taran, 3 -- Gulf of Gdansk, 4 -- Hel Spit, 5 -- Cape Rozewie, 6 -- Slupsk Furrow, 7 -- Hoburg Channel, 8 -- Gotland Deep.}
\end{figure}
%%
However, even in the case of a straight coastline and no topographic elevations at the bottom, the transient alongshore wind favorable for coastal upwelling (downwelling) drives an alongshore baroclinc front and associated longshore jet-like current of the same direction, which, in turn, can eventually break into randomly located mesoscale cyclonic (anticyclonic) eddies due to baroclinic flow instability
[\itc{Zhurbas et al.,} \reflink{Zhur06}{2006}].
However, eddy formation in this case is much slower than in the case of topographic features/headlands producing finite-amplitude disturbances of the mean flow. Another example of the predictable formation of cyclonic eddies in the Southeastern Baltic Sea is the area of the S{\l}upsk Furrow. During northeasterly winds, saltier and denser water of the North Sea origin rushes to the east in the bottom layer of the furrow
[\itc{Krauss and Br\"{u}gge,} \reflink{Krau91}{1991}],
forming a gravitational current above which cyclonic eddies are generated both in the furrow and east of it
[\itc{Zhurbas et al.,} \reflink{Zhur04a}{2004a}, \reflink{Zhur12}{2012}].
Thus, for the reasons outlined above, it is expected that a numerical model with sufficiently high (submesoscale permitting) horizontal resolution, even without data assimilation, can reproduce the main features of a pattern of real (observed) mesoscale eddies in a relatively shallow coastal sea. In this context, the term ``pattern'' means ``realization'' implying the reproduction of specific mutual location of eddies in space and their characteristics (i.e., size, intensity, and vorticity sign) at given time instants rather than a statistical description of the eddy field, which is inherent to turbulence studies.

A good opportunity of validating the simulated mesoscale eddy patterns can be provided by satellite imagery, since there are a lot of remote sensing images of the Southeastern Baltic Sea displaying mesoscale eddies in radar, infrared, and visible bands
(e.g.,
[\itc{Dabuleviciene et al.,} \reflink{Dabu18}{2018};
\itc{Gurova et al.,} \reflink{Guro13}{2013};
\itc{Ginzburg et al.,} \reflink{Ginz15a}{2015a}, \reflink{Ginz15b}{2015b}, \reflink{Ginz17}{2017};
\itc{Karimova et al.,} \reflink{Kari12}{2012};
\itc{Lavrova et al.,} \reflink{Lavr18}{2018}]).


\begin{sloppypar}
The objective of this study is to implement a circulation model of the Southeastern Baltic Sea with extremely high horizontal resolution (of the order of 100 m) capable of simulating mesoscale/submesoscale eddies, providing simulations for a warm season when the eddies can be better visualized in true color and infrared images due to cyanobacteria blooms and high contrasts of the sea surface temperature (SST), and comparing the simulated eddy patterns against satellite images available.
\end{sloppypar}

\enlargethispage{-.8\baselineskip}

\vspace{-.5\baselineskip}

\section{2. Material and Methods}
\vspace{-.5\baselineskip}

\subsection{2.1. Model Setup}
\vspace{-.5\baselineskip}

\begin{sloppypar}
The General Estuarine Transport Model (GETM)
[\itc{Burchard and Bolding,} \reflink{Burc02}{2002}]
was applied to simulate the meso- and submesoscale variability of temperature, salinity, currents, and overall dynamics in the Southeastern Baltic Sea. GETM is a primitive equation, three-dimensional, free surface, hydrostatic model with an embedded vertically adaptive coordinate scheme
[\itc{Gr\"{a}we et al.,} \reflink{Graw15}{2015};
\itc{Hofmeister et al.,} \reflink{Hofm10}{2010}].
The vertical mixing is parameterized by the two-equation $k-\varepsilon$ turbulence model coupled with an algebraic second-moment closure
[\itc{Burchard and Bolding,} \reflink{Burc01}{2001};
\itc{Canuto et al.,} \reflink{Canu01}{2001}].
The implementation of the turbulence model is performed via the General Ocean Turbulence Model (GOTM)
[\itc{Umlauf and Burchard,} \reflink{Umla05}{2005}].
\end{sloppypar}

\begin{sloppypar}
A uniform horizontal grid of the high-resolution nested model (0.125 nautical miles (n.m.) or approximately 230~m) covers the computational domain, which includes the central Baltic Sea along with the Gulf of Finland and Gulf of Riga
 (\figref{1}),           %     (Figure 1),
while 60 adaptive layers are used in the vertical direction. The digital topography of the Baltic Sea was taken from the Baltic Sea Bathymetry Database
(\href{http://data.bshc.pro/}{http://data.bshc.pro/})
with a 0.5 nautical mile step and interpolated to ensure finer resolution.
\end{sloppypar}

The model simulation run was performed from 1 April  2015 to 9 October 2015. The model domain has the western open boundary in the Arkona Basin and the northern open boundary in the entrance to the Bothnian Sea. The one-way nested approach was used for the open boundary conditions, and results from a coarser resolution model were utilized at the boundaries. The coarse resolution model covered the entire Baltic Sea with the open boundary in the Kattegat and had a horizontal resolution of 0.5~n.m. (approximately 1~km) in the entire model domain. More detailed information about the coarse resolution model is available from
[\itc{Zhurbas et al.,} \reflink{Zhur18}{2018}].\vspace*{.4ex}


\begin{sloppypar}
The atmospheric forcing (the wind stress and surface heat flux components) was calculated from the wind, solar radiation, air temperature, total cloudiness, and relative humidity data generated by the High Resolution Limited Area Model (HIRLAM) maintained operationally by the Estonian Weather Service with the spatial resolution of 11 km and temporal resolution of 1~hour
[\itc{M\"{a}nnik and Merilain,} \reflink{Mann07}{2007}].
The wind velocity components at the 10-m level along with other HIRLAM meteorological parameters were interpolated to the model grids.
\end{sloppypar}


\begin{sloppypar}
The freshwater input from 54 largest Baltic Sea rivers together with their inter-annual variability is taken into account in the coarse resolution model. The original dataset consists of daily climatological values of the discharge for each river, but inter-annual variability is added by adjusting the freshwater input to different basins of the sea to match the values reported annually by HELCOM
[\itc{Johansson,} \reflink{Joha18}{2018}]. The high-resolution model takes into account only those rivers that flow into the sea within the model domain.
\end{sloppypar}


The initial thermohaline field was taken from the coarse resolution model for 1 April  2015 and was interpolated to the high-resolution model grid. The prognostic model runs started from the motionless state and zero sea surface elevation. The spin-up time of the southern Baltic Sea model under atmospheric forcing was expected to be within 10 days
[\itc{Krauss and Br\"{u}gge,} \reflink{Krau91}{1991};
\itc{Lips et al.,} \reflink{Lips16}{2016}],
while the model output was used for comparison with satellite imagery after 45 days of modeling time.

\subsection{2.2. Remote Sensing Data}

The numerical model results were compared with optical and infrared satellite imagery of the Southeastern Baltic Sea acquired in the period between 15 May  and 25 August  2015. The following sensors were used to construct maps in true color, suspended matter, and SST:

\begin{enumerate}
\item
 Visible and Infrared Imager/Radiometer Suite (VIIRS) installed at the Joint Polar Satellite System (JPSS) series of spacecraft, including the Suomi National Polar-orbiting Partnership (S-NPP)  (VIIRS-SNPP radiometer, spatial resolution 750~m);
\item
\begin{sloppypar}
 Landsat Enhanced Thematic Mapper Plus (ETM+) sensor onboard the Landsat 7 satellite (ETM+ Landsat~7, spatial resolution 30~m);
 \end{sloppypar}
\item
 Thermal Infrared Sensor (TIRS) onboard the Landsat 8 satellite (TIRS Landsat~8, spatial resolution 30~m);
\item
 Multispectral Instrument (MSI) onboard the Sentinel-2 satellite (MSI Sentinel-2, spatial resolution 10--60~m);
\item
 Moderate Resolution Imaging Spectroradio\-meter onboard the Terra and Aqua satellites (MODIS TERRA, MODIS AQUA, spatial resolution 250--1000~m).
\end{enumerate}

In total, 44 high resolution and more than 300 moderate resolution optical and infrared satellite images for the period 15 May -- 25 August  2015, were examined, and only eight of them were chosen for further analysis and comparison with the simulations (see Chapter~3).

\section{3. Results}

Hereafter we present a number of satellite images displaying mesoscale eddies in the Southeastern Baltic Sea for the period of May--August 2015 versus simultaneous maps of simulated oceanographic parameters, such as the temperature and current velocity in the surface layer, as well as the patterns formed by the simulated floating Lagrangian particles that were initially uniformly distributed on the sea surface, in one day of advection. The patterns formed by the simulated floating Lagrangian particles were shown to be a powerful tool to visualize mesoscale/submesoscale structures on the sea surface
[\itc{V\"{a}li et al.,} \reflink{Vali18}{2018}].   %%% ??? +
The surface layer was the uppermost adaptive layer of the model
[\itc{Gr\"{a}we et al.,} \reflink{Graw15}{2015}]
whose depth in the Southeastern Baltic Sea varied basically within
0.5--1.5~m. Inertial oscillations that could mask the eddies were removed from the current velocity field by means of a simple procedure described by
[\itc{Zhurbas et al.,} \reflink{Zhur06}{2006}].
In the next figures, the eddies, which are identified both in model and satellite data, are marked as $C_i$ in the case of a cyclone, and $A_i$ in the case of an anticyclone ($i=1, 2, 3, \ldots$ is the serial number of the identified eddy). Tracking the same eddy in a series of model and/or satellite images (note that there is no problem to access the model snapshots with time offset as small as needed, e.g., one day or one hour), we can thoroughly investigate the eddy's fate including the place and time of its origin/disappearance, trajectory, and life time.

\begin{sidewaysfigure}                      % Fig 2
\figurewidth{140mm}
\setimage{}{}{158mm}{}{2019es000672-f02}
\fcaplabel{(a) IR image (SST, \deg C) of the Southeastern Baltic Sea from the VIIRS-SNPP radiometer on 15 May  2015, 10:54 GMT, and (b) modeled current velocity and temperature of the surface layer for the same date and time.}
\end{sidewaysfigure}


The period of May--August 2015 was taken for the analysis because it was well covered by the satellite optical and IR imagery, which made it possible to follow and describe eddies evolution in the Gulf of Gdansk during 1.5 months
[\itc{Ginzburg et al.,} \reflink{Ginz17}{2017}].

The IR image of the Southeastern Baltic Sea from the VIIRS-SNPP radiometer on 15 May  2015
 (\figref{2}a)           %     (Figure 2, a)
displays a pair of a larger anticyclone in the northern part of the Gulf of Gdansk (marked as A1) and a smaller cyclone northwest of it (marked as C1). This vortex pair (dipole) is clearly identified in the modeled current velocity and temperature fields of the surface layer almost at the same location
 (\figref{2}b).           %     (Figure 2, b).
Moreover, two other cyclonic eddies seen in the modeled map can be identified in the IR image (though not as clear as the vortex pair) (marked as C2 and C3 in
 \figref{2}).           %     Figure 2).
Besides, a warm northeast coastal current is well reproduced by the model between the Hel Spit and the Latvian coast, including meanders in the Gulf of Gdansk and filaments offshore the Curonian Spit and northward of Klaipeda
 (\figref{2}).           %     (Figure 2).

\def\dbltopfraction{.95}

\begin{figure}                      % Fig 3
\figurewidth{100mm}
\setimage{}{}{113mm}{}{2019es000672-f03}
\fcaplabel{(a) IR image (SST, \deg C) of the Southeastern Baltic Sea from VIRS-NPP on 31 May  2015 vs simultaneous modeled fields of the surface layer parameters: (b)~current velocity and temperature, (c) temperature, and (d) floating Lagrangian particles.}
\end{figure}


\begin{figure}                      % Fig 4
\figurewidth{90mm}
\setimage{}{}{115mm}{}{2019es000672-f04}
\fcaplabel{(a) IR image (SST, \deg C) of the Southeastern Baltic Sea from the TIRS Landsat~8 on 8 June  2015 vs (b, c, and d) simultaneous modeled patterns of the surface layer parameters (the same as in
\figref{3}).}          % Figure 3).
\end{figure}


The foregoing model images (prior to 15 May  2015) show that the A1--C1 vortex pair was eventually formed on 7 May  2015 from a protuberance that protruded northwards from the tip of the Hel Spit after a southeasterly wind gale. This is a place of frequent generation of dipole structures as observed by the IR, optical, and SAR imagery
[\itc{Ginzburg et al.,} \reflink{Ginz15a}{2015a}, \reflink{Ginz15b}{2015b}, \reflink{Ginz17}{2017}].
On 11 May  soon after generation of the A1--C1 vortex pair, a new cyclone, C2, of the same place of origin, i.e., the tip of the Hel Spit, formed. In the course of time, the C2 cyclone moved northwards for $\sim 70$~km and merged with the C1 cyclone on 20 May (the result of merging of the C1 and C2 cyclonic eddies is referred to as the C2 cyclone, see
 \figref{3}--\figref{8}).           %                %     Figures 3-8).

On 31 May  the A1--C2 vortex pair was clearly identified both in the satellite imagery and model output
 (\figref{3}).           %     (Figure 3).
The C3 cyclone, previously seen near the tip of the Hel Spit
(cf. \figref{2})           %     Figure 2)
moved east to the southeastern slope of the Gulf of Gdansk. Another smaller cyclone, C4, was observed northwest of the large cyclone C2.

 \figref{4}           %     Figure 4
refers to 8 June  showing remarkable correspondence between the vortex field detected by the satellite and the modeled pattern. A large A1--C2 vortex pair (the same as in
 \figref{3})           %     Figure 3)
is observed at the center, and a relatively small cyclone, C4, is located west of it. In the eastern part of the Gulf of Gdansk, the numerical model predicts an intrusion of warm water, which surrounds the Hel Spit from the east
 (\figref{4}b, \figref{4}c).           %     (Figure 4b, c).
In the satellite image, the same area is covered by clouds, which hide the SST pattern including the C3 cyclone, but indirectly confirm the presence of warm water responsible for convective cloudiness right over the warm intrusion
 (\figref{4}a).           %     (Figure 4a).
It is worth noting that the C4 cyclone, being poorly identified in the modeled current velocity, is really well seen in the floating particle simulation
(see
 \figref{3}d
 and
  \figref{4}d).           %                %     Figures 3d and 4d).

\begin{figure}                      % Fig 5
\figurewidth{100mm}
\setimage{}{}{113mm}{}{2019es000672-f05}
\fcaplabel{(a) True color image of the Southeastern Baltic Sea from MSI Sentinel-2 on 28 June  2015 vs (b, c, and d) simultaneous modeled patterns of the surface layer parameters (the same as in
\figref{3}
and
\figref{4}).}          % Figures 3 and 4).
\end{figure}


By 28 June  the C2 cyclone, a former constituent of the A1--C2 vortex pair in
 \figref{3}           %     Figures 3
and
 \figref{4},           %     4,
departed from the A1 anticyclone and went to the north, being replaced by a newborn cyclone, C5, of the same place of origin -- the tip of the Hel Spit
(see
 \figref{5}).           %     Figure 5).
Note that the model on 28 June  displayed another cyclonic eddy, C6, west of the A1--C2 vortex pair, which is not present in the true color satellite image in
 \figref{5}a           %     Figure 5a
as it was located outside the western boundary of the image. The C6 cyclone had a different origin than C1, C2, and C5 -- it came from the S{\l}upsk Furrow. An interesting feature of the true color satellite and the simulated floating particle images of 28 June  is the presence of warm bubble-like structures, which occupy the eastern coastal zone
 (\figref{5}a           %     (Figure 5, a
and
 \figref{5}d).           %     d).
Keeping in mind that the preceding period of 23--26 June  was characterized by the downwelling-favoring southwesterly wind, which promotes a coastal current of the same direction, one may assume that the warm bubbles were formed as a result of coastal downwelling relaxation.

The same A1--C5 vortex pair is most clearly identified in both the satellite imagery and the model output of  1 and 3 July, while the identification of the C2 and C6 cyclones in satellite imagery of 1 July  is less convincing
(see
 \figref{6}           %     Figures 6
 and
  \figref{7}).           %     7).
Nevertheless, the true color image of 3 July  from the OLI Landsat 8 displays three other cyclones in the Eastern Gulf of Gdansk, C7, C8, and C9, that can be also seen in the simulation
 (\figref{7}).           %     (Figure 7).

In the long run, by August 10, according to the model output, the A1 anticyclone disappeared, and a new large anticyclone, A2, formed instead
 (\figref{8}).           %     (Figure 8).
 \figref{8}           %     Figure 8
shows remarkable correspondence between the A2 anticyclone positions detected from the satellite and from the simulation. Several cyclonic eddies of smaller size can be identified around the A2 anticyclone both in satellite and model images, namely, C10, C11, C12, and C13. The most long-living of them, C10, was born near the Hel Spit, made an anticyclonic loop north of the Gulf of Gdansk, and faded not far from the place of origin. The C11 and C13 cyclonic eddies, being most clearly seen in the true color image as spirals
 (\figref{8}a),           %     (Figure 8a),
are weakly identified in the modeled currents
 (\figref{8}b),           %     (Figure 8.b),
but are clearly identified in the simulated floating particles
 (\figref{8}d).           %     (Figure 8d).

\fitopen{.99}
\begin{sloppypar}
Finally, by 23 August  the hydrometeorological conditions in the Southeastern Baltic Sea have drastically changed: a tight northeasterly wind promoted coastal upwelling and a related jet-like coastal current, which, in turn, shed cyclonic eddies behind/downstream topographic features such as the Sambian Rise (C14 and C15 cyclones in
 \figref{9}),           %     Figure 9),
Cape Taran (C16 and C17), and Cape Rozewie (C18).
\end{sloppypar}

\begin{sloppypar}
The information about the evolution of the above-mentioned mesoscale eddies including the date and place of their birth/decay, life path, and life time is summarized in
\tabref{1}.           %    Table 1.
\end{sloppypar}
\fitclose

\subsection{4. Discussion and Conclusions}


\begin{sloppypar}
Satellite monitoring of the Southeastern Baltic Sea is seriously affected by cloudiness, which prevents receiving information on ocean dynamics. According to the authors' estimates based on the data obtained from the online Giovanni system developed and maintained by the Goddard Earth Sciences Data and Information Services Center at NASA
(\href{http://disc.sci.gsfc.nasa.gov/giovanni}{http://disc.sci.gsfc.nasa.gov/giovanni}),
this region is partially or completely covered by clouds approximately 60\%
of the year
[\itc{Ginzburg et al.,} \reflink{Ginz15a}{2015a}].
SAR imagery allows detecting some signatures of ocean dynamics in cloudy conditions, but there are certain restrictions caused by weather conditions and problems in the interpretation of SAR images
[\itc{Lavrova et al.,} \reflink{Lavr11}{2011}, \reflink{Lavr16}{2016}].
In such circumstances, high resolution numerical models seem to be the best tool to reproduce daily mesoscale and sub-mesoscale fields of currents in the cloudy conditions typical for the Baltic Sea.
\end{sloppypar}


The above-presented comparison of the satellite imagery against the model output revealed reasonably good correspondence between the observed and simulated patterns of the mesoscale eddy fields in the surface layer of the Southeastern Baltic Sea for the period of May--August 2015. In our opinion, it has two consequences. First, the ultra-high resolution circulation model we implemented can be used for the hindcast and operational forecast of mesoscale eddy fields in the Southeastern Baltic Sea. Second, as the model has reproduced a handful of satellite images available for the period, one can expect that the simulation provides realistic information on the evolution of mesoscale and submesoscale eddies in the area, including the place and time of birth, size, intensity, track, life time, etc.

\begin{figure}                      % Fig 6
\figurewidth{100mm}
\setimage{}{}{112mm}{}{2019es000672-f06}
\fcaplabel{(a) Satellite image of the Southeastern Baltic Sea from MODIS-AQUA in water-leaving radiance on 1 July  2015 vs (b, c, and d) simultaneous modeled images of the surface layer parameters (the same as in
\figref{3}--\figref{5}).}                  % Figures 3-5).
\end{figure}


\begin{figure}                      % Fig 7
\figurewidth{100mm}
\setimage{}{}{112mm}{}{2019es000672-f07}
\fcaplabel{(a) True color satellite image of the Southeastern Baltic Sea from OLI Landsat~8 on 3 July  2015 vs (b, c, and d) simultaneous modeled patterns of the surface layer parameters (the same as in
\figref{3}--\figref{6}).}                  % Figures 3-6).
\end{figure}



\begin{figure}                      % Fig 8
\figurewidth{100mm}
\setimage{}{}{112mm}{}{2019es000672-f08}
\fcaplabel{(a) True color satellite image of the Southeastern Baltic Sea from MSI Sentinel-2 on 10 August  2015 vs (b, c, and d) simultaneous modeled patterns of the surface layer parameters (the same as in
\figref{3}--\figref{7}).}                  % Figures 3-7).
\end{figure}



\begin{figure}                      % Fig 9
\figurewidth{90mm}
\setimage{}{}{115mm}{}{2019es000672-f09}
\fcaplabel{(a) Suspended matter image of the Southeastern Baltic Sea from MODIS-TERRA (23 August  2015, 09:35 GMT) vs (b, c, and d) simultaneous modeled patterns of the surface layer parameters (the same as in
\figref{3}--\figref{8}).}                  % Figures 3-8).
\end{figure}


According to modeling results presented in
\tabref{1},           %   Table 1,
the tip of the Hel Spit is a special place for the Southeastern Baltic Sea where long-living mesoscale eddies and vortex pairs frequently originate. For example, the A1--C1 vortex pair was generated from a protuberance that protruded northward from the tip of the Hel Spit after a southeasterly wind gale. Then a newborn cyclone, C2, broke away from the tip of the Hel Spit, moved northward to meet and merge with C1, and thereby formed a modified vortex pair, A1--C2. In the course of time, the C2 cyclone detached from the A1 anticyclone and moved further north, while a newborn cyclone, C5, of the same place of origin (the tip of the Hel Spit) approached A1 to form the A1--C5 vortex pair. Note that the formation of a vortex pair north of the Hel Spit was first reported in
[\itc{Zhurbas et al.,} \reflink{Zhur04b}{2004b}]
based on IR satellite images and model results.

\begin{sidewaystable}                                % Table 1
\tablewidth{158mm}
\caption{\scriptsize Evolution of Mesoscale Eddies Identified in the Southeastern Baltic Sea From Satellite Imagery and Model Output for the Period May--August 2015}
\vspace{5pt}
\fitopen{.7}
\tiny
\begin{tabular}
{@{}l@{\hspace{24pt}}
c@{\hspace{24pt}}
c@{\hspace{24pt}}
c@{\hspace{24pt}}
c@{}}
\hline
\\ [-7pt]
Eddy  &   Date and place of birth &     Life path               &       Date and place of decay  &     Life time, \\
ID    &                           &                             &       or position at the end   &     days       \\
         &                           &                             &      of the modeling period    &  \\
[7pt]  \hline   \\ [-4pt]
A1  &     07.05.2015, Tip of the Hel Spit& Within the Gulf               & 20.07.2015, Gulf      &  74 \\
    &                                    & of Gdansk                     &    of Gdansk           & \\
A2  &     01.08.2015, north of the Gulf  & Within the Gulf               & 20.08.2015, Western            &  20 \\
    &     of Gdansk                      & of Gdansk                     & Gulf of Gdansk                 & \\
C1  &     07.05.2015, Tip of  the        & N transition for $\sim 70$ km,& 22.05.2015, North of the       &  15 \\
    &     Hel Spit                       &  merger with C2               &  Gulf of Gdansk                &  \\
C2  &     11.05.2015, Tip of the Hel Spit& NNE transition                &  28.07.2015, Eastern slope     &  78 \\
    &                                    & for $\sim 150$ km             &  of the Hoburg Channel         &  \\
C3  &     13.05.2015, Tip of the Hel Spit& E transition for              &  09.06.2015, Tip of the        &   24  \\
    &                                    & $\sim 50$ km and back         &  Hel Spit                      &  \\
C4  &     21.05.2015, Southern           & S transition for $\sim 80$ km &  12.06.2015, North of the      &   22 \\
    &     Hoburg Channel                 &                               &  Hel Spit                      &  \\
C5  &     14.06.2015, near the           & Anticyclonic semi-loop        &  25.07.2015, near the          &   40 \\
    &     Cape Rozewie                   & North of the Gulf of Gdansk    &  Cape Taran                    &  \\
C6  &     10.06.2015, came from the      & NEE transition for $\sim 100$ km& 09.10.2015, on the beam      &   $>120$ \\
    &     S{\l}upsk Furrow               &                                 & of Klaipeda                  &  \\
C7  &     01.07.2015, near               & Insignificant transition        & 06.07.2015, near             & 5 \\
    &      the Cape Taran                &                                 & the Cape Taran               & \\
C8  &     16.06.2015, near               & SW transition for $\sim 50$ km  & 09.07.2015, Southern Gulf  &22  \\
    &      the Cape Taran                              &                                 &       of Gdansk            & \\
C9  &     30.06.2015, Southern Gulf      & Insignificant transition        & 05.07.2015, Southern Gulf & 5  \\
    &     of Gdansk                      &                                 &   of Gdansk                                 &  \\
C10 &     07.07.2015, near the Hel Spit  & Anticyclonic loop north         & 08.09.2015, Northern Gulf& 61  \\
    &                                    & of the Gulf of Gdansk           &  of Gdansk                                  &  \\
C11 &     03.08.2015, near the Cape Taran& Insignificant transition        & 12.08.2015, near the Cape Taran    & 9  \\
C12 &     06.08.2015, Tip of the Hel Spit& NW transition for $\sim 40$ km  & 15.08.2015, went for West of      & $>8$  \\
    &                                    &                                 & the Cape Rozewie                  &  \\
C13 &    09.08.2015, North of the Hel Spit& Insignificant transition       & 11.08.2015, North of the Hel Spit &   2 \\
C14 &    22.08.2015, near the Sambian Rise& Irregular transition,          & 04.09.2015, near the Sambian Rise &   12  \\
    &                                     & merger with C15                &                                   &  \\
C15 &    19.08.2015, near the Sambian Rise& Irregular transition,          & 09.10.2015, northwest of the      &  $>51$  \\
    &                                     & merger with C14                & Sambian Rise                      &  \\
C16 &    20.08.2015, near the Cape Taran  &SW transition for $\sim 30$ km, & 31.08.2015, Gulf of Gdansk        &  11 \\
    &                                     &merger with C17                 &                                   &  \\
C17 &    20.08.2015, Southeast of the     & N transition for $\sim 110$ km &  09.10.2015, North of the         &   $>50$ \\
    &    Gulf of Gdansk                   &                                &   Gulf of Gdansk                   &  \\
C18 &    21.08.2015, near the Cape Rozewie& Irregular transition           &  09.10.2015, near Cape the Rozewie &   $>49$ \\
[4pt] \hline
\end{tabular}
\end{sidewaystable}

The following question arises: What makes the tip of the Hel Spit a special place where mesoscale eddies are formed? Taking a look at a close-up of the Southeastern Baltic Sea bathymetry
 (\figref{1},           %     (Figure 1,
right panel), one can admit that the tip of the Hel Spit is characterized by an exceptionally steep bottom slope, so that the 70~m depth contour (the level of the upper boundary of the permanent halocline) approaches the shore as close as 1--2 km. In such circumstances, relatively short-term exposure to the upwelling/downwelling favorable wind promotes a perturbation (bending) of the permanent halocline; as a result, mesoscale vortices are more powerful and long-living.

The eddies listed in
\tabref{1}           %   Table 1
are mostly cyclones rather than anticyclones (16 cyclones vs 2 anticyclones). However, this does not mean that the number of cyclonic eddies in the surface sea layer exceeds that of anticyclonic eddies. In our opinion, the only thing that is beyond doubt is that the cyclonic eddies are better visualized on satellite and simulated images, having a larger angular velocity of rotation and a longer life time (e.g.
[\itc{V\"{a}li et al.,} \reflink{Vali17}{2017}]),
and manifesting themselves in the form of easily noticeable spirals
[\itc{Munk et al.,} \reflink{Munk00}{2000}].

Finally, the results of simulation of submesoscale circulation in the Southeastern Baltic Sea are very promising for numerical modeling of oil pollution in this region. Since 2004, this area is under daily satellite monitoring of oil pollution related to oil production at the D-6 offshore oil platform operated by Lukoil-Kaliningradmorneft Ltd
[\itc{Bulycheva et al.,} \reflink{Buly14}{2014}, \reflink{Buly16}{2016};
\itc{Kostianoy,} \reflink{Kost06}{2006}, \reflink{Kost15}{2015};
\itc{Krek et al.,} \reflink{Krek18}{2018}].
In the framework of this satellite monitoring, the Seatrack Web (STW) oil spill model of the Swedish Meteorological and Hydrological Institute (SMHI) is used for modeling the drift and transformation of observed oil spills since 2004. Today, STW is the official HELCOM forecast and hindcast system, which is used by national authorities and research organisations for calculating the fate of oil spills in the Baltic Sea since the 1990s.  Our 15 year-long experience in using STW showed that the model does not provide a correct forecast for oil spill drift in this region in the presence of mesoscale and submesocale eddies, which are shown above to be characteristic of the Gulf of Gdansk
[\itc{Ginzburg et al.,} \reflink{Ginz15b}{2015b}].
Based on the above-reported results, the increase in the horizontal resolution in the SMHI operational circulation model coupled with STW is believed to be capable of significantly improving the quality of the forecast of oil spill drift in mesoscale and submesoscale eddy fields.

\acknowledgements{This research, including numerical modeling and analysis of the results, was performed by Victor Zhurbas in the framework of the Shirshov Institute of Oceanology RAS budgetary financing (Project No. 149-2019-0003) and was supported in part by the Russian Foundation for Basic Research (Grant No. 18-05-80031).
Germo V\"{a}li was supported by institutional research funding IUT 19-6 of the Estonian Ministry of Education and Research and BONUS, the joint Baltic Sea research and development programme (Art 185), funded jointly from the European Union's Seventh Programme for research, technological development and demonstration, and from the Estonian Research Council through grant VEU17107 (BONUS \mbox{INTEGRAL}). An allocation of computing time from the High Performance Computing cluster at the Tallinn University of Technology and University of Tartu is gratefully acknowledged. Andrey Kostianoy and Olga Lavrova (analysis of satellite imagery) were partially supported by budgetary financing of the Shirshov Institute of Ocea\-nology RAS (Project No. 149-2019-0004) and the Space Research Institute RAS (theme ``Monitoring'', state register No. 01.20.0.2.00164), respectively.}

\begin{references}
\reference{Buly16} Bulycheva,  E. V.,  A. V.  Krek,  A. G.  Kostianoy, A. V. Semenov, A. Joksimovich (2016), Oil pollution in the Southeastern Baltic Sea by satellite remote sensing data in 2004--2015, {\itshape Transport and Telecommunication}, {\itshape 17}, no. 22, p. 155--163,  \href{https://doi.org/10.1515/ttj-2016-0015}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Buly14} Bulycheva, E., I. Kuzmenko, V. Sivkov (2014), Annual sea surface oil pollution of the south-eastern part of the Baltic Sea by satellite data for 2006--2013, {\itshape Baltica, Special Issue}, {\itshape 27}, p. 9--14,  \href{https://doi.org/10.5200/baltica.2014.27.10}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Burc01} Burchard, H., K. Bolding (2001), Comparative Analysis of Four Second-Moment Turbulence Closure Models for the Oceanic Mixed Layer, {\itshape J. Phys. Oceanogr.}, {\itshape 31}, p. 1943--1968,  \href{https://doi.org/10.1175/1520-0485(2001)031%3C1943:CAOFSM%3E2.0.CO;2}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Burc02} Burchard, H., K. Bolding (2002), {\itshape GETM -- a general estuarine transport model. Scientific documentation. Technical report EUR 20253en.}, European Commission, Ispra, Italy.

\reference{Canu01} Canuto, V. M., et al. (2001), Ocean Turbulence. Part I: One-Point Closure Model-Momentum and Heat Vertical Diffusivities, {\itshape J. Phys. Oceanogr.}, {\itshape 31}, p. 1413--1426,  \href{https://doi.org/10.1175/1520-0485(2001)031%3C1413:OTPIOP%3E2.0.CO;2}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Dabu18} Dabuleviciene,  T.,  I.   Kozlov,  D.  Vaiciute, I. Dailidiene (2018), Remote sensing of coastal upwelling in the South-Eastern Baltic Sea: Statistical properties and implications for the coastal environment, {\itshape Remote Sensing}, {\itshape 10}, no. 11, p. 1752,  \href{https://doi.org/10.3390/rs10111752}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Elki13} Elkin, D. N., A. G. Zatsepin (2013), Laboratory Investigation of the Mechanism of the Periodic Eddy Formation Behind Capes in a Coastal Sea, {\itshape Oceanology}, {\itshape 53}, no. 1, p. 24--35,  \href{https://doi.org/10.1134/S0001437012050062}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Ginz15a} Ginzburg, A. I., E. V. Bulycheva, A. G. Kostianoy, D. M. Solovyov (2015a), Vortex Dynamics in the Southeastern Baltic Sea from Satellite Radar Data, {\itshape Oceanology}, {\itshape 55}, no. 6, p. 805--813,  \href{https://doi.org/10.1134/S0001437015060065}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Ginz15b} Ginzburg, A. I., E. V. Bulycheva, A. G. Kostianoy, D. M. Solovyov (2015b), On the role of vortices in the transport of oil pollution in the Southeastern Baltic Sea (according to satellite monitoring), {\itshape Current Problems in Remote Sensing of the Earth From Space}, {\itshape 12}, no. 3, p. 149--157 (in Russian).

\reference{Ginz98} Ginzburg, A. I., A. G. Kostianoy, D. M. Soloviev, S. V. Stanichny (1998), Upwelling cyclonic eddies off south-western tip of the Crimea, {\itshape Issledovanie Zemli iz Kosmosa}, {\itshape 3}, p. 83--88 (in Russian).

\reference{Ginz17} Ginzburg,  A. I.,  E. V.  Krek,  A. G.  Kostianoy, D. M. Solovyev (2017), Evolution of mesoscale anticyclonic vortex and vortex dipoles/multipoles on its base in the south-eastern Baltic (satellite information May--July 2015), {\itshape Journal of Oceanological Research}, {\itshape 45}, no. 1, p. 10--22,  \href{https://doi.org/10.29006/1564-2291.JOR-2017.45(1).3}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Graw15} Gr\"{a}we, U., P. Holtermann, K. Klingbeil, H. Burchard (2015), Advantages of vertically adaptive coordinates in numerical models of stratified shelf seas, {\itshape Ocean Modelling}, {\itshape 92}, p. 56--68,  \href{https://doi.org/10.1016/j.ocemod.2015.05.008}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Guro13} Gurova, E., A. Lehmann, A. Ivanov (2013), Upwelling dynamics in the Baltic Sea studied by a combined SAR/infrared satellite data and circulation model analysis, {\itshape Oceanologia}, {\itshape 55}, no. 3, p. 687--707,  \href{https://doi.org/10.5697/oc.55-3.687}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Hofm10} Hofmeister, R., H. Burchard, J.-M. Beckers (2010), Non-uniform adaptive vertical grids for 3D numerical ocean models, {\itshape Ocean Modelling}, {\itshape 33}, no. 1--2, p. 70--86,  \href{https://doi.org/10.1016/j.ocemod.2009.12.003}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Joha18} Johansson, J. (2018), Total and regional runoff to the Baltic Sea, HELCOM Baltic Sea Environment Fact Sheets, Online data, HELCOM, Helsinki, Finland (\href{http://www.helcom.fi/baltic-sea-trends/environment-fact-sheets/}{http://www.helcom.fi/baltic-sea-trends/ environment-fact-sheets/}).

\reference{Kari12} Karimova, S. S., O. Yu. Lavrova, D. M. Solov'ev (2012), Observation of Eddy Structures in the Baltic Sea with the Use of Radiolocation and Radiometric Satellite Data, {\itshape Izvestiya, Atmospheric and Oceanic Physics}, {\itshape 48}, no. 9, p. 1006--1013,  \href{https://doi.org/10.1134/S0001433812090071}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Kost17} Kostianoy, A. G. (2017), Satellite monitoring of the ocean climate parameters. Part 1, {\itshape Fundamental and Applied Climatology}, {\itshape 2}, p. 27--49,  \href{https://doi.org/10.21513/2410-8758-2017-2-63-85}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Kost06} Kostianoy, A. G., et al. (2006), Operational satellite monitoring of oil spill pollution in the southeastern Baltic Sea: 18 months experience, {\itshape Environmental Research, Engineering and Management}, {\itshape 4}, no. 38, p. 70--77.

\reference{Kost15} Kostianoy, A. G., E. V. Bulycheva, A. V. Semenov, A. V. Krainyukov (2015), Satellite monitoring systems for shipping, and offshore oil and gas industry in the Baltic Sea, {\itshape Transport and Telecommunication}, {\itshape 16}, no. 2, p. 117--126,  \href{https://doi.org/10.1515/ttj-2015-0011}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Krau91} Krauss, W., B. Br\"{u}gge (1991), Wind-produced water exchange between the deep basins of the Baltic Sea, {\itshape J. Phys. Oceanogr.}, {\itshape 21}, p. 373--384,  \href{https://doi.org/10.1175/1520-0485(1991)021%3C0373:WPWEBT%3E2.0.CO;2}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Krek18} Krek, E., A. Kostianoy, A. Krek, A. V. Semenov (2018), Spatial distribution of oil spills at the sea surface in the Southeastern Baltic Sea according to satellite SAR data, {\itshape Transport and Telecommunication}, {\itshape 19}, no. 4, p. 294--300,  \href{https://doi.org/10.2478/ttj-2018-0024}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Lavr18} Lavrova, O. Yu., E. V. Krayushkin, K. R. Nazirova, A. Ya. Strochkov (2018), Vortex structures in the Southeastern Baltic Sea: Satellite observations and concurrent measurements, Proc. SPIE 10784, {\itshape Remote Sensing of the Ocean, Sea Ice, Coastal Waters, and Large Water Regions 2018, 1078404 (5 October 2018)},  \href{https://doi.org/10.1117/12.2325463}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Lavr11} Lavrova, O. Yu., A. G. Kostianoy, S. A. Lebedev, M. I. Mityagina, A. I. Ginzburg, N. A. Sheremet (2011), {\itshape Complex Satellite Monitoring of the Russian Seas}, 470 pp., IKI RAN, Moscow (in Russian).

\reference{Lavr16} Lavrova, O. Yu., M. I. Mityagina, A. G. Kostianoy (2016), {\itshape Satellite Methods of Detection and Monitoring of Marine Zones of Ecological Risks}, 336 pp., Space Research Institute, Moscow (in Russian).

\reference{Lehm08} Lehmann, A., K. Myrberg (2008), Upwelling in the Baltic Sea -- A review, {\itshape J. Mar. Syst.}, {\itshape 74}, p. S3--S12,  \href{https://doi.org/10.1016/j.jmarsys.2008.02.010}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Lips16} Lips, U., V. Zhurbas, M. Skudra, G. V\"{a}li (2016), A numerical study of circulation in the Gulf of Riga, Baltic Sea. Part I: Whole-basin gyres and mean currents, {\itshape Cont. Shelf Res.}, {\itshape 112}, p. 1--13,  \href{https://doi.org/10.1016/j.csr.2015.11.008}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Mann07} M\"{a}nnik, A., M. Merilain (2007), Verification of different precipitation forecasts during extended winter-season in Estonia, {\itshape HIRLAM Newsletter }, {\itshape 52}, p. 65--70.

\reference{Munk00} Munk, W. H., L. Armi, K. Fischer, F. Zachariasen (2000), Spirals on the sea, {\itshape Proc. R. Soc. Lond.}, {\itshape A 456}, p. 1217--1280.

\reference{Myrb03} Myrberg, K., O. Andrejev (2003), Main upwelling regions in the Baltic Sea -- A statistical analysis based on three-dimensional modeling, {\itshape Boreal Environment Research}, {\itshape 8}, no. 2, p. 97--112.

\reference{Umla05} Umlauf, L., H. Burchard (2005), Second-order turbulence closure models for geophysical boundary layers. A review of recent work, {\itshape Cont. Shelf Res.}, {\itshape 25}, no. 7--8, p. 795--827,  \href{https://doi.org/10.1016/j.csr.2004.08.004}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Vali17} V\"{a}li, G., V. Zhurbas, U. Lips, J. Laanemets (2017), Submesoscale structures related to upwelling events in the Gulf of Finland, Baltic Sea (numerical experiments), {\itshape J. Mar. Syst.}, {\itshape 171}, no. SI, p. 31--42,  \href{https://doi.org/10.1016/j.jmarsys.2016.06.010}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Vali18} V\"{a}li, G., V. M. Zhurbas, J. Laanemets, U. Lips (2018), Clustering of floating particles due to submesoscale dynanics: a simulation study for the Gulf of Finland, {\itshape Fundamentalnaya i Prikladnaya Gidrofizika}, {\itshape 11}, no. 2, p. 21--35,  \href{https://doi.org/10.7868/S2073667318020028}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Vign11} Vignudelli,  S.,  A.   Kostianoy,  P.  Cipollini, J. Benveniste (eds.) (2011), {\itshape Coastal Altimetry}, 578 pp., Springer-Verlag, Berlin, Heidelberg,  \href{https://doi.org/10.1007/978-3-642-12796-0}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Zhur12} Zhurbas, V., J. Elken, V. Paka, J. Piechura, I. Chubarenko, N. Golenko, S. Shchuka (2012), Structure of unsteady overflow in the Slupsk Furrow of the Baltic Sea, {\itshape J. Geophys. Res. -- Oceans}, {\itshape 117}, no. C04027, p. 1--17,  \href{https://doi.org/10.1029/2011JC007284}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Zhur17} Zhurbas, V. M., N. P. Kuzmina, D. A. Lyzhkov (2017), Eddy formation behind a coastal cape in a flow generated by transient longshore wind (numerical experiments), {\itshape Oceanology}, {\itshape 57}, no. 3, p. 350--359,  \href{https://doi.org/10.1134/S0001437017020229}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Zhur06} Zhurbas, V., I. S. Oh, T. Park (2006), Formation and decay of a longshore baroclinic jet associated with transient coastal upwelling and downwelling: A numerical study with applications to the Baltic Sea, {\itshape J. Geophys. Res.}, {\itshape 111}, no. C04014, p. 1--18,  \href{https://doi.org/10.1029/2005JC003079}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Zhur04a} Zhurbas, V., T. Stipa, P. M\"{a}lkki, I. Hense, V. Sklyarov, et al. (2004a), Generation of subsurface cyclonic eddies in the southeast Baltic Sea: Observations and numerical experiments, {\itshape J. Geophys. Res.}, {\itshape 109}, no. C05033, p. 1--12,  \href{https://doi.org/10.1029/2003JC002074}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\reference{Zhur04b} Zhurbas, V. M., T. Stipa, P. M\"{a}lkki, V. T. Paka, N. P. Kuzmina, V. E. Sklyarov (2004b), Mesoscale Variability of Upwelling in the South-East Baltic: Infrared Images and Numerical Modeling, {\itshape Oceanology}, {\itshape 44}, no. 5, p. 619--628.

\reference{Zhur18} Zhurbas, V., G. V\"{a}li, M. Golenko, V. Paka (2018), Variability of bottom friction velocity along the inflow water pathway in the Baltic Sea, {\itshape J. Mar. Syst.}, {\itshape 184}, p. 50--58,  \href{https://doi.org/10.1016/j.jmarsys.2018.04.008}{\textcolor{brown}{\space\bfseries Cross\-ref}}

\end{references}
\end{document}
