%% Copernicus Publications Manuscript Preparation Template for LaTeX Submissions
%% ---------------------------------
%% This template should be used for copernicus.cls
%% The class file and some style files are bundled in the Copernicus Latex Package which can be downloaded from the different journal webpages.
%% For further assistance please contact the Copernicus Publications at: publications@copernicus.org
%% http://publications.copernicus.org


%% Please use the following documentclass and Journal Abbreviations for Discussion Papers and Final Revised Papers.


%% 2-Column Papers and Discussion Papers
\documentclass[os, manuscript]{copernicus}



%% Journal Abbreviations (Please use the same)

% Atmospheric Chemistry and Physics (acp)
% Advances in Geosciences (adgeo)
% Advances in Statistical Climatology, Meteorology and Oceanography (ascmo)
% Annales Geophysicae (angeo)
% ASTRA Proceedings (ap)
% Atmospheric Measurement Techniques (amt)
% Advances in Radio Science (ars)
% Advances in Science and Research (asr)
% Biogeosciences (bg)
% Climate of the Past (cp)
% Drinking Water Engineering and Science (dwes)
% Earth System Dynamics (esd)
% Earth Surface Dynamics (esurf)
% Earth System Science Data (essd)
% Fossil Record (fr)
% Geographica Helvetica (gh)
% Geoscientific Instrumentation, Methods and Data Systems (gi)
% Geoscientific Model Development (gmd)
% Geothermal Energy Science (gtes)
% Hydrology and Earth System Sciences (hess)
% History of Geo- and Space Sciences (hgss)
% Journal of Sensors and Sensor Systems (jsss)
% Mechanical Sciences (ms)
% Natural Hazards and Earth System Sciences (nhess)
% Nonlinear Processes in Geophysics (npg)
% Ocean Science (os)
% Primate Biology (pb)
% Scientific Drilling (sd)
% SOIL (soil)
% Solid Earth (se)
% The Cryosphere (tc)
% Web Ecology (we)


\begin{document}

\linenumbers

\title{Modelling origin and transport fate of waste materials on the south-eastern Adriatic coast (Croatia)}


% \Author[affil]{given_name}{surname}

\Author[1]{Martina}{Tudor}
\Author[2]{Ivica}{Janekovi\'{c}}

\affil[1]{Croatian Meteorological and Hydrological Service, Gri\v{c} 3, Zagreb, Croatia}
\affil[2]{Rudjer Bo\v{s}kovi\'{c} Institute, Bijeni\v{c}ka 54, Zagreb, Croatia}

%% The [] brackets identify the author with the corresponding affiliation. 1, 2, 3, etc. should be inserted.



\runningtitle{Modelling waste transport}

\runningauthor{Tudor and Janekovi\'{c}}

\correspondence{Martina Tudor (tudor@cirus.dhz.hr)}



\received{}
\pubdiscuss{} %% only important for two-stage journals
\revised{}
\accepted{}
\published{}

%% These dates will be inserted by Copernicus Publications during the typesetting process.


\firstpage{1}

\maketitle



\begin{abstract}
The south-eastern parts of the Adriatic Sea coastline were severely polluted by large amounts of accumulated waste material in the second half of November 2010.  The waste,  reported by major news agencies, accumulated dominantly during 21 of November 2010 by favourable wind - ocean current transport systemi (the East Adriatic Current from southeast).
In the study we analysed meteorological and oceanographic conditions that lead to the waste deposition using available in situ measurements, remote sensing data as well as numerical models of the ocean and the atmosphere.
The measured data reveal that an intensive rainfall event from 7 till 10 of November 2010, over the parts of Montenegro and Albania, was followed by a substantial increase of the river water levels indicating a possibility of flash floods that could have splashed the waste material into a river and after to the Adriatic Sea (or to the sea directly).
The currents that can bring this waste to Croatia could have been intensified by the strong wind from south-east.
In order to test these two hypotheses we set a number of numerical drifter experiments with trajectories initiated off the coast of Albania during the intensive rainfall events following their faith in space and time.
The numerical drifter trajectory experiments that resulted with drifters that reached the right position (south--eastern Adriatic coast) and time (exactly by the time the waste was observed) were initiated on 00 and 12 UTC of 10 November 2010, during the mentioned flash flood event.
\end{abstract}



\introduction  %% \introduction[modified heading if necessary]
On the 21 of November 2010, a dramatic waste accumulation has been widely reported by public media (web news agencies, television, radio, daily papers) at the south--eastern coast of Croatia, particularly area of Pelje\v{s}ac Peninsula; islands Mljet, Kor\v{c}ula and Lastovo as well as in numerous inlets and beaches northwest of Dubrovnik (see map of the area in Fig~\ref{terrain}).
The heaps of waste were composed mostly of plastic packages, glass bottles, clothes and other typical floating municipal garbage while labels suggested that some part of the waste arrived from Albania. It is a country on the south-eastern Adriatic coast 100 km south-east of the area where waste accumulated. 
The labels are no proof that the waste originated from Albania.
One can easily imagine a heap of waste accumulated on a sea shore anywhere on Earth with labels "Made in China" since it is well known for its massive production and export.
This study describes a possible chain of events that lead to waste accumulation on beaches in south-east Croatia.
It is not unusual that a few pieces of waste reach Croatian coast in a late autumn, however the event was several orders of magnitude larger than any other in previous years (according to the local officials, there were no reports in the media).

The meteorological conditions during October and November in southern Adriatic included several episodes of intensive precipitation that initiated flash floods in Montenegro and Croatia (there were no reports avalable for Albania).
A flash flood event could have washed the waste to the sea (or first to a lake or a river that would eventually take it to the sea).
There was no rainfall data from Albania available through standard international data exchange so remote sensing data and NWP model data were used to estimate which intensive precipitation events could have initiated a flash flood there.
The high precipitation events that could have initiated flash flood in the area are identified by combination of in-situ and remote sensing data.

Numerical modeling studies that examine how a floating entity reached a certain position by means of atmosphere and sea driven currents have been done before \citep{beg,doos,liu}. The subjects range from explanation of how floating sweet potato reached Polynesia from South America \citep{monte2008}, spread of oil spills such as the one following the Deepwater Horizon disaster has received more attention ( for a collection of articles see http://deep waterhorizon.nooa.gov/) as well as the floating debris that was washed to the sea by tsunami following the Tohoku 9Mw eathquake on 11 March 2011 (see http://www.marinedebris.noaa.gov/).

In our case we use meteorological and ocean models to explore how waste items, once washed to the sea, floated to the area affected by the accumulated waste.
During this study we assume that strong south-easterly wind intensified the sea current system that was favourable to bring the floating waste materials. 
Previously, we also assumed that the waste was splashed into the sea by strong flash floods as a consequence of the severe torrential rain. These hypotheses are further investigated using all available meteorological, oceanographic and
hydrological observational data as well advanced meteorological and oceanographic numerical models.

The next section describes the geographical characteristics of the area, meteorological and oceanographic conditions, measured data and models used in this study. Results of model simulations are presented in Section 3 and summarized in conclusions in Section 4.

\section{Region, Data, and Models}
\subsection{Region}
The Adriatic Sea is a narrow sea, connected to the Mediterranean only by the Otranto Strait at the southern part. Bathymetry varies over the basin, with the northern part as the shallowest, with mean depth of 35 m, the central and southern Adriatic are significantly deeper and divided by Palagru\v{z}a sill (Fig~\ref{terrain}).
The central region reaches up to 280 m depth in Jabuka Pit.
The southern region is the deepest, up to 1200 m in the South Adriatic Pit (SAP).

The Adriatic Sea surface flow is predominately cyclonic orientation \citep{cushman2001} with
distinct current regime of East Adriatic Current (EAC) flowing northwest along the eastern coast characterized with salty and warm water from the Ionian Sea.
During the rain seasons EAC is further intensified with the outflow of the Albanian rivers creating region of fresh water (ROFI) dynamics \citep{burrage2008,burrage2009}. In the central region the sea surface flow typically bifurcates east of the Palagruza Sill (eg. \citet{wolf1887}) enhancing the cyclonic circulation in the southern Adriatic \citep{artegiani,horton1997}.
On the other side of the Adriatic Sea there is a Western Adriatic Current (WAC) holding fresher and colder water along the western coast. It carries a signature of Po river outflow, the most important source of fresh water in the whole Adriatic Sea.

On the land, the area is surrounded by Apennines in the west, Dinaric Alps and high mountains of Montenegro and Albania along eastern coast while on the northern coast reaches low and flat Po valley. Mountains are much closer to the shore on the eastern side of the Adriatic Sea, with several peaks higher
than 1.5 km located less than 10 km from the coast (Fig~\ref{terrain}).
Those mountains have a strong effect on the air flow and atmospheric dynamics \citep{mesing1981} and consequently define the sea current response as well.

Mediterranean cyclones often traverse the area \citep{horvath2009} from south.
However cyclones often form in the Genoa Bay, at the northwest \citep{mesing1981} traverse the Tyrhennian Sea and continue to the east possibly supporting cyclone development and intensification in the Adriatic Sea at the east and Ionnian Sea at the south \citep{alpert}.
These cyclones usually cross the Adriatic Sea but in a certain synoptic conditions can support development of a separate Adriatic cyclone \citep{horvath2008} and other mesoscale weather activity.  In that case two of these cyclones can coexist forming a system of twin cyclones in which moist air converges
and generates large quantities of available precipitable water \citep{lionello2006}.

The intensive dynamics found in the area also supports strong wind \citep{horvath2011,bajic2007,brankovic2008} development with the most severe
and gusty wind from northeast named bura (see \cite{grisog2009} for a review), as well the local wind from southeast refered as jugo \citep{jurcec1996}.
Strong bura or strong jugo can last for several days inducing strong response in the Adriatic Sea \citep{kuzmic, dorman}.
Onset, duration and spatial distribution of wind strength in a bura or jugo episode is controlled by
an interaction of the synoptic and/or mesoscale forcing with local topography \citep{ivateks2004,pasaric2007,tudor2010}.
The bura strength varies significantly in space and time forming usually several stripes of strong wind across the Adriatic Sea separated by areas of milder wind \citep{grubi2004} related to mountain gaps and ridges upstream.
On the other hand, jugo blows along shore, it is more steady and relatively warm wind related to a Genoa cyclone \citep{jurcec1996} or mesoscale cyclone above northern Adriatic \citep{brzomah1999,brzovic1999}.
Jugo can be associated with sirocco. However, the terms jugo and sirocco are not synonyms. The latter is a southern wind blowing from Sahara in advance of low pressure moving eastward across southern Mediterranean Sea. Jugo is SE wind over Adriatic, while in sirocco episodes, wind that brings the Sahara desert dust over Adriatic can be from south or southwest.

Climatologicaly, Southern Adriatic region is characterized with warm and dry summers and mild and wet winters \citep{zaninovic2008}.
The area receives abundant precipitation amounts as Crkvice in Montenegro holds the maximum measured
on the European continent \citep{magas2002}.
Parts of the northern Albania is rich in precipitation but unevenly distributed in space and time.
In those regions precipitation can be further intensified by increased aerosol concentration \citep{koren2012} usually advected to the Adriatic Sea area from the Sahara desert by the sirocco wind.

It is worth to say that wind forcing, when pronounced, dominate over all other forcing contributions and dynamically shape the sea surface currents system found in the Adriatic Sea.
The surface wind jets and wakes of the bura wind have a profound effect on the surface currents
\citep{orlic1994,pullen2003}, while jugo wind is well known to influence WAC flow reversals
\citep{orlic2007,poulain2004}.
It is therefore important to drive the ocean model with a high resolution wind feild that resolves high resolution wind features that develop due to interaction of large scale dynamics with local mountains surrounding the sea.

\subsection{Data}

In order to test our hypothesis and numerical model results we used available remote sensing data and in situ measurements.
For the meteorological part we used SYNOP, climatological and rain-gauge measurements from Croatia, Montenegro, Italy, Greece and Macedonia.
At the time of the event (November 2010) there were no in-situ measured data available from Albania thorugh standard meteorological network, the data from airport did not contain measurements of precipitation. 
The hydrological analysis was based on the water level measurements on relevant major rivers in Montenegro and Macedonia used to confirm intensive precipitation as possible cause of the flush flood events.
Since these rivers flow through Albania or along its border on ther way to the sea, this is a good indicator of the state of the rivers in Albania.

Remote sensing data, used in this study, originate from Meteosat Second Generation (MSG), specifically  from The EUMETSAT Network of Satellite Application Facillities (NWC SAF)\footnote{products available on the  http://www.eumetrain.org/}.
Satellite derived precipitation data are used as provided from the Tropical Rainfall Measuring Mission (TRMM, \citet{huffman2007}), in particular we used the diurnal accumulated precipitation data from the 3B42 product and 3 hourly precipitation intensity data from 3B40RT product.
Precipitation can be enhanced by the presence of aerosols.
The two sets of aerosol data presented in this study are the aerosol optical thickness (AOT) from Moderate Resolution Imaging Spectroradiometer (MODIS, \citet{remer2008}) aboard Aqua satellite and Ozone Monitoring Instrument
aboard NASA's Earth Observing System (EOS)  Aura satellite (OMI, \citet{torres2002},\citet{veihel2007})\footnote{The OMI and TRMM data are available from Giovanni web server interface \citep{acker2007} on http://disc.sci.gsfc.nasa.gov}.
The wind over the sea surface derived from MetOp ASCAT \citep{bentamy2012,bentamy2012a} was used to evaluate 10 meter wind field from the meteorological model.

\subsection{Models}
\subsubsection{Atmospheric model -- ALADIN}
The numerical weather prediction (NWP) model data used in this study originate from the operational 8 km resolution forecast runs using ALADIN limited area model (Aire Limit\'{e}e Adaptation Dynamique d\'{e}veloppement InterNational, \citet{aladin}) with a specific local 3D--var data assimilation \citep{stanesic2011}.
In autumn 2010, operational forecast run twice per day up to 72 hours in advance starting from 00 and 12 UTC analyses.
The model forecast in 8 km resolution used initial and boundary conditions from global model ARPEGE (Action de Recherche Petite Echelle Grande Echelle, \citet{cassou2001}) run operationally in Meteo France.
The operational high-resolution dynamical adaptation \citep{ivateks2004}
provides forecast of 10 m wind adapted to local and upstream topography \citep{horvath2011}.
Unfortunately, this method provides only wind field at high resolution, but not the other meteorological variables needed to force the ocean model.
The meteorological model 10 m wind field is obtained by vertical interpolation from the lowest model level (17m above sea), see \cite{geleyn1988} for more details. 

%The model wind field quality was evaluated using the wind over the sea surface derived from MetOp ASCAT \citep{bentamy2012,bentamy2012a}.
In order to simulate the mesoscale characteristics and development of the low pressure field, a 2km resolution forecast using the non-hydrostatic set-up of the ALADIN model  and the full parametrization set, including radiation, microphysics and convection schemes \citep{tudor2010} was used to model the state of the atmosphere.  The high-resolution forecast uses scale selective digital filter initialization \citet{termonia2008} and no data assimilation to initialize the model fields. It is coupled to the ALADIN 8 km resolution with 3 hour interval.
This might be insufficient to prevent the fastest of the meteorological features to enter the domain unnoticed by the lateral boundary coupling procedure \citep{tudor2010b} and possibly miss or undersample a storm rapidly entering the domain through the lateral boundaries.
Since the southern Adriatic is not very far from the southern lateral boundary of the high resolution domain, model could have underestimated a storm arriving from south through Otranto strait, however this would be a short duration event related to a flash flood but too short to affect the sea currents substantially.

\subsubsection{Ocean model}

The quality of simulated currents on the ocean surface depends on the wind field.
Wind field over Adriatic is variable in both space and time, depends on surrounding topography and events with strong and severe wind are better forecast in high resolution NWP models \citep{ivateks2004,brankovic2008,tudor2010}.
The ocean dynamics as a response to the atmospheric forcing was computed using REgional Ocean modelling System (ROMS, \cite{schep2005}) numerical model. ROMS model belongs to free surface, Boussinesq and hydrostatic approximation models that solves primitive equations using curvilinear finite difference grids. Model was
forced with ALADIN meteorological model data (10m wind, 2m temp and relative humidity, sea level pressure, rainfall rate, short wave radiation and cloud fraction), climatological values for the Adriatic river run-offs and open boundary values with daily temperature, salinity, currents and sea level information from AREG (INGV) Mediterranean model. The advection scheme for tracers (temperature and salinity) is based on multidimensional positive definite advection transport algorithm - MPDATA \citep{smolarkiewitz} while for momentum on 3rd order upwind scheme. More details of model implementation for the Adriatic Sea are described in \cite{janekovic}.

The sea surface currents, responsible for waste transport, are computed using 2 km resolution ROMS ocean model and were used for virtual drifter trajectory simulations. Drifters are set to the surface layer, without vertical dynamics, ensuring representation of floating waste material.
For computing numerical drifer tajectories we used Roms 
OFfline Floats (ROFF) package \citep{carr}, available for ROMS 
comunity (http://web.atmos.ucla.edu/~capet/Myresearch/my\_research\_floats.html)

\section{Results and discussion}

\subsection{Atmospheric model results}

To estimate the convective rainfall rate and precipitating clouds we used derived fields from the NWC SAF products focused on studied area and time, rain-gauge measurements and TRMM rainfall data.
The NWC SAF precipitating clouds (PC, \citet{pc}) field provides precipitation probabilities and the convective rainfall rate in mm/hour (CRR \citet{crr}) is computed assuming that clouds being both high and with a large vertical extent are more likely to induce rain (see http://www.nwcsaf.org/ for more details).
CRR gives estimate of intensive rainfall from convective clouds, but PC is usefull estimate of rainfall from other types of clouds (eg. nimbostratus).
According to the available rain-gauge measurements, 6 hourly PC and CRR fields and TRMM rainfall data, there were several heavy rainfall events in the month preceding 21 November 2010 that could have caused flash floods in the area of southeast Adriatic coast and inland.

We identify those events as four episodes: 23-25 October (E1), 2-3 (E2), 8-10 (E3) and 17-18 of November 2010 (E4).

The large scale synoptic conditions responsible for meteorological setup are described using ERA Interim \citep{dee2011} re-analysis fields. It turns out that
on 24 of October 2010, low pressure system entered western Mediterranean from Atlantic, deepened and formed a cyclone, centred over Genoa bay.
The next day the pressure decreased further and the associated southern wind
strengthened from northern Africa to Adriatic causing intensive rainfall over the eastern Adriatic coast (E1). The cyclone moved southeast on 26 and initiated severe bura wind first on northern Adriatic and later spread over the whole Adriatic Sea by 27 of October 2010.
An ensemble of trajectories initiated over southeast Adriatic on 12 UTC, 25 of October 2010 were used to test if this severe rainfall event was the one that flushed the waste to the sea. The results of these trajectory computations are described later in the text (as experiment 2 in Section 3.3.1)

Another cyclone from the 1 to the 4 of November 2010 (E2) moved from the Genoa bay southeastward, causing strong jugo wind over the Adriatic Sea (Fig~\ref{wind}).  The rainfall was the most intensive over the northern Italy and central Adriatic region with most of the rainfall above the Adriatic sea.
Northern Adriatic received more than 100 mm of precipitation within 24 hours, while the rain was weak in the southeastern region of our interest (Fig~\ref{raingauge}). Consequently, E2 case was omitted from further analysis as was too weak to initiate a flash flood in southeast Adriatic. In the following days, meteorological situation was stable with weak pressure gradient, low wind as well high pressure over western Mediterranean inducing moderate winds from northwest.

The weather changed again in the period from 7 till 10 of November 2010 (E3), dominated by a large scale cyclone (Fig~\ref{erai850}) that
arrived from the Northern Atlantic causing sirocco wind over Mediterranean (the colour of the wind vectors in Fig~\ref{erai850} indicates wind speed in m/s) that brought warm air and Sahara dust from the northern Africa (often found in rain gauges after intensive rainfall events), aerosols are also shown in Fig~\ref{erai850}. 
Over Adriatic, the wind was strong to severe from southwest and south direction (Fig~\ref{wind}).
The wind direction was well forecast by the model, but at Palagru\v{z}a, Dubrovnik, Prevlaka wind speed was underestimated, while at Mljet and Biokovo the observed wind speed was correctly modelled.
Pressure measurements reveal that during this event the Adriatic Sea was subject to a deep cyclone that last for several days (Fig~\ref{wind}) with a strong pressure gradient over the Adriatic sea.

The precipitation intensity was estimated using the PC and CRR fields that showed strong convection and rainfall in the afternoon and evening having periods with weak to moderate rain intensity in the night and early morning.
The precipitating clouds covered much of the area, while the convective rainfall rate is far more localized and very intensive.
It is important to note that PC and CRR fields were available on 6 hour interval, while heavy rainfall could have occurred outside the sampling interval and easily could have been missed.
The 24 hourly precipitation exceeds 100 mm over northern Albania in TRMM precipitation estimates (shown as squares in Fig~\ref{rain1}) and measurements at several rain-gauges in Montenegro for two consecutive days (raingauge measurements are shown as circles in Fig~\ref{rain1}) as well as in the model forecast (shaded background in Fig~\ref{rain1}). 
Measurements from the rain-gauges showed that rainfall during the E3
was the most intensive on stations in Montenegro (larger circles and stars on Fig~\ref{raingauge}) hence on southeast Adriatic coast and significantly more intensive than in other episodes in November 2010.
For example, on 10 November 2010, there was 188.1 mm of rain measured at Cetinje and 143 mm measured at Golubovci station, both in Montenegro.
 Accumulated precipitation data are shown on maps for 9 and 10 November 2010 (for all available stations, circles in Fig~\ref{rain1}). There are 3 intensive precipitation events before 21 November 2012, and measured precipitation exceeded 100mm/24 hours in Montenegro only in the event 7 to 10 November 2012 (Fig~\ref{raingauge}. 
The forecast of the accumulated 24 hourly rainfall corresponds to the values measured on rain-gauges, although the model exaggerated slightly the rainfall on the coastline and underestimated the rainfall on several locations further inland (Fig~\ref{rain1}).
Wind measurements (Fig~\ref{wind}) show that wind in E3 episode was from south direction, more energetic and lasted longer than for other strong wind episodes during November 2010.

After E3, in the period from 11 till 15 of November 2010, the weather was mostly dry with weak to moderate wind and direction typical of the sea breeze diurnal cycle (Fig~\ref{wind}). In the next days a cyclone formed in the Genoa bay (15 of November 2010)
supporting again strong jugo wind over the whole Adriatic Sea (Fig~\ref{eraiuvnp}). Moreover, the wind strengthened (Fig~\ref{wind}) with prevailing direction from southeast as measured in Dubrovnik, Mljet and Prevlaka (Fig~\ref{wind16}). Precipitation was intensive with peaks above 100 mm within 24 hours on 17 November (E4), but the maxima were localized on northeast Adriatic. The ALADIN model forecast had similar rainfall distribution, as a result, E4 was omitted from detailed analysis as a period favourable with respect to the flash flood. However, it is important as the wind field has driven the currents on the sea surface. During E4 wind was stronger on the coast (Dubrovnik and Prevlaka) than in the off-shore region (Mljet), as a consequence of channelling effect the coastal mountains (Fig~\ref{wind16}). Model yields stronger wind above open sea (thiner arrows on regular grid in Fig~\ref{wind16}) over southern Adriatic than MetOp ASCAT wind data (thicker arrows in Fig~\ref{wind16}) for 16 November 2010, but the wind strength and direction are correct for 17 November 2010 (Fig~\ref{wind16}). The global pressure gradient over the Mediterranean and Central Europe supported the wind regime from south and southeast over the whole southern Adriatic (Fig~\ref{eraiuvnp}).
Later, by 19 of November 2010, the wind changed direction to southwest.
A cyclone moved from Atlantic southeast, to the western Mediterranean.
The wind changed to strong and severe jugo wind on the 21 of November 2010.

Based on the analysis of available precipitaion fields, it turns out that E3 episode was the one when intensive rainfall occurred over Albania and surrounding countires and was the most likely event that triggered a flash flood.


\subsection{Hydrology}

Annual river run-off distribution for the Albanian rivers usually varies for an order of magnitude during the year
(embedded diagram in Fig~\ref{rivers}) with one pronounced peak in November and another in January.
Bojana river collects the water flowing from Drim river and Skadar lake and flows into Adriatic along the border between Albania and Montenegro.

The largest lake in the region, the Skadar lake, is filled by river Mora\v{c}a and Crnojevica in Montenegro and drained into Bojana river (the name is Bojana in Montenegro and Buna in Albania).
Bojana river also receives Drim river as a major tributary on the way to the Adriatic Sea.
Drim (Drim in Montenegro, Drin in Albania) river
powers 3 hydroelectric power plants in Albania. Downstream it splits into two flows,
the smaller one reaches the Adriatic sea directly, and the larger part flows into Bojana river.

The river water level measurements on  the rivers in Montenegro that belong to the Adriatic Sea catchment (Fig~rivers \ref{rivers}) increase substantialy for the E3 episode.
The water level surge is most intense for the rivers that fill the Skadar lake.
The level of Bojana river raised as well during the same event.
This is followed by a rise of 1.5 m in the water level of the Skadar lake.
Bojana river level rose before the level of Lake Skadar, this could have happened due to an increase in contribution from the Drim river tributary. The water levels of Skadar lake and Bojana river stayed high until the end of November 2010.

Those measurements suggest that the event (E3) from 8 till 10 of November 2010 was capable of flushing the waste material into the Adriatic sea or any of the rivers in the area that flow into it.

\subsection{Ocean model results}

Ocean model results show (Fig~\ref{roms}) consistent development of strong surface northwest currents after strong jugo wind episodes and small eddies close to the eastern Adriatic coast in the periods of weak wind forcing. As stated before in the text, during the November 2010, we can find three periods of different wind conditions over the Adriatic Sea.  The first one from 7 till 11, when strong south-east wind generated strong north-west current system in south-east Adriatic (Fig~\ref{roms}a) followed by a weak wind period when the sea-current transport was weaker (Fig~\ref{roms}b) and the ocean model formed a pool of colder water in southeast Adriatic in the area where Bojana river enters the Adriatic Sea. 
Finally, the period with moderate to strong south-east wind (Fig~\ref{roms}c) strengthened the north-west current. This was most likely responsible for waste transport and deposition.

\subsubsection{Numerical drifter trajectories}

%% Previously, we assumed that a severe rainfall event initiated flash-flood that washed the waste to the sea.
The waste consisted of floating items, so its movement was computed as virtual floating drifters released in southeast Adriatic. 
Trajectories of virtual drifters were computed using ROFF package and surface currents from 2km resolution ROMS run.
The computation of the trajectories of the drifters stopped when they reached coastline for the first time.
A number of numerical drifter experiments was set in which trajectories were initiated in southeast Adriatic on the 12 UTC, 19 (experiment 1) and 25 of October 2010 (experiment 2), and then sequentially at 00 and 12 UTC on each day starting from 8 till 12 of November 2010 (experiments 3 -- 11).
All virtual drifters were released within a polygon covering an area over southern Adriatic.
The initial points of virtual drifter trajectories are separated by 0.01 degree ($\approx$ 1 km) along longitude and latitude (total 3071 drifters) and fill a polygon with longitude and latitude coordinates of south-west corner (19.1,41.0) and north-east corner (19.4,41.9) that cover a portion of south-eastern Adriatic Sea in the vicinity of coast of Albania and Montenegro (Fig~\ref{trajs}).
Furthermore we divided the polygon into 9 areas (A1,...A9) to better cluster track different regions, hence possible source origin.
The drifters starting from different areas are plotted in different colours, as marked on the Fig~\ref{trajs}.
The longitude and latitude coordinates of sotuthwest and northeast points of the polygons are
\begin{itemize}
\item A1: SW (19.1,41.0), NE (19.2,41.3) shown in red,
\item A2: SW (19.2,41.0), NE (19.3,41.3) shown in green,
\item A3: SW (19.3,41.0), NE (19.4,41.3) shown in blue,
\item A4: SW (19.1,41.3), NE (19.2,41.6) shown in magenta,
\item A5: SW (19.2,41.3), NE (19.3,41.6) shown in cyan,
\item A6: SW (19.3,41.3), NE (19.4,41.6) shown in yellow,
\item A7: SW (19.1,41.6), NE (19.2,41.9) shown in dark green,
\item A8: SW (19.2,41.6), NE (19.3,41.9) shown in orange,
\item A9: SW (19.3,41.6), NE (19.4,41.9) shown in gray,
\end{itemize} 
A plot of drifter positions was done with 6 hourly interval for each experiment (not shown), and the summary trajectories are shown in Fig~\ref{trajs}.

\begin{itemize}
\item In the experiment 1 the drifters were released at 12 UTC, 19 of October 2010 and were first pushed offshore into EAC. It turns out that a considerable number of drifters originated from regions A7, A8 and A9 reached Croatian coast and Mljet island already on 27 of October 2010.
Drifters from regions A4, A5 and A6 reached Mljet channel by 3 of November, but were pushed back south--east in the following days. Those drifters continued further to the north--west and finally accumulated on the islands much further northwest than observed (Fig~\ref{trajs}a). There were no reports of significant accumulation of waste on the Croatian coast that would be a consequence of this event. In that sense, we can reject the hypothesis that this rainfall event was the one that caused the flash flood that got the waste material to the sea.

\item For experiment 2 the drifters were released  at 12 UTC on 26 of October 2010.
Soon the dirfters were advected in the westward direction. When drifters entered EAC, they moved more to the north--west and were deposed on Mljet Island already on 9 of November 2010. However, several drifters starting from A6 region entered Mljet channel on 18 and later deposed on Pelje\v{s}ac on 21 of November (Fig~\ref{trajs}b). Based on those results we can assume that it is possible but unlikely that the rainfall event on 26 October 2010 has initiated
the chain of events that led to severe waste disposal in the region.

\item Drifters initiated on 8 November (both 00 and 12 UTC -- experiments 3 and 4) mostly arrived to south--east Adriatic Sea coast, the northern Albania and Montenegro as soon as on the 11 of November 2010 (Figures \ref{trajs}c and \ref{trajs}d) as a consequence of sea current system supported by strong southern and SSW wind blowing on 8 and 9 of November 2010.
Furthermore, strong wind changed direction into NW on 12 and 13th of November 2010 generating currents that transported numerical drifters off the coast, resulting only with a small number of them (initiated from A4 region), to reach Mljet island and coastline north--west of Dubrovnik by 18 of November 2010.
The rest of the drifters dominantly stayed in the south--east region, while only a small number of them moved north--eastward not entering Mljet Channel, but instead floated much closer to coast, at the end finally accumulated in the Kolo\v{c}ep Channel and Ston bay area  on 25 of November 2010.

\item Quickest drifters initiated from A7, A4 and A8 regions at 00 and 12 UTC on 9 November 2010 (experiments 5 and 6) reached Mljet Island and entered Mljet Channel already on 17 and 18, while a majority of drifters
from other areas accumulated in the Ston bay after 21 of November 2010 (Figures \ref{trajs}e and \ref{trajs}f).

\item Small number of drifters from A7, A8 and A9 regions, released at 00 and 12 UTC on 10 November 2010 (experiments 7 and 8) reached Mljet, Dubrovnik and Kolo\v{c}ep channel by 18, while other drifters initiated from the same area accumulate on the Croatian shores on 21 of November 2010 (Figures \ref{trajs}g and \ref{trajs}h). The drifters initiated further south lagged behind former ones and approached the affected area on 22 of November 2010. Those dates were the most commonly reported in the media as the onset of severe pollution at the Croatian coast.

\item The last three sets of drifters, released at 00 and 12 UTC on 11 and 00 UTC on 12 November 2010 (experiments 9--11) were first pushed westward off the Albanian coast and stayed in the area off shore of Albania and
Montenegro for a few days. Later they were transported into southward direction on 14 to 16 November 2012 (Fig~\ref{trajs}i). Apparently, EAC was detached, at that time, from the shore and it's typical path. As a result, drifters from A1 and A2 regions arrived further north--westward than drifters initiated more to the north or closer to the coastline.

\end{itemize}


One should bear in mind that the drifter trajectories do not allow us to assign a single event in space and time as the moment when the waste was disposed to the sea. There is ambiguity in the position as a result of unresolved physics, unperfect meteorological model and initial conditions used to force ocean counterpart, missing dynamics in the ocean model introduced with a lack of wave--current interaction, spatial model resolutions in the narrow channels etc.
However, the results do show that there is a possibility and is the most probably that the heavy rain on 9 and 10 of November 2010 washed the waste into the sea (or first to a river that carried it to the sea by that date).
The computations further show that not all the numerical drifters initiated in south-east Adriatic inevitably ended on the coastline of south--east of Croatia. Surface sea currents enhanced by the wind forcing can carry the waste back to the shore, or to the closer coastline of Montenegro.
Otherwise, different meteo--ocean conditions can push the waste off shore, and EAC can carry the waste to central or even north Adriatic, or in some cases back to the southern regions of the Adriatic Sea.
However, none of the trajectories initiated in our experiments crossed the Adriatic Sea and approached to Italy, which is probably due to an absence of intensive bura events during the studied period.

\conclusions  %% \conclusions[modified heading if necessary]

The oceanographic and meteorological conditions that lead to a severe deposition of waste material on the south--eastern Adriatic Sea coast on 21 of November 2010, are studied using ALADIN -- metorological and ROMS -- ocean numerical models along with available measurements. We tried to answer what, where and when was the cause for the event, the initial points presented in this study were limited since the labels that indicated some of the items originate from Albania and the local current system is north-west.

Based on the meteorological simulations and satellite derived precipitation reveal several intensive rainfall events that could have initiated flash floods in Albania and presumably flush the waste material to the rivers and later to the Adriatic Sea. Moreover, measured and NWP model rainfall data shows that the rain was more intensive over the Albania in the event from 8 till 10 of November 2010 (E3) than in the other intensive rainfall events that occurred in the studied area during the 4 weeks before the reported waste accumulation.

Measured wind speed during the episode E3 was strong to severe from southern direction, however slightly underestimated by the operational ALADIN model forecast at several land locations. Improvements in the atmospheric model resolution could resolve those issues as noted in \cite{signell}.
Since strong wind influences the surface currents that advect the drifters, this could have an impact on the computed trajectories of virtual drifters.
It is interesting to note that based on the ASCAT estimated wind data, for the 16 of November 2010, ALADIN model wind speed was larger than the measured one (Fig~\ref{wind16}).
During the last studied period (E4), wind from observations as well from the model, was from south--east with weaker magnitudes than in the E3 period. In the E4 case the strongest wind was found over the open sea, in the south--east region of Mljet Island as well south of Dubrovnik (Fig~\ref{wind16}). This event is a typical jugo wind episode which further enhanced the sea surface current system -- responsible for transport of the waste material to the Croatian waters and finally to the shore on 21 of November 2010.



\begin{acknowledgements}
The authors would like to thank the Meteorological and Hydrological Service of Montenegro for providing the meteorological and water level measurements as well AREG -- INGV, Bologna, Italy.
The authors are grateful to NASA for providing valuable satellite derived products throught the GIOVANNI web interface as well as TRMM, OMI and MODIS scientists and developers. The NOAA NESDIS CoastWatch and NOAA SWFSC ERD are acknowledged for providing Metop ASCAT wind data. Dr. Janekovi\'{c} was sponsored through project 098-0982705-2707 of Croatian Ministry of Science. Martina Tudor thanks the Croatian Ministry of Science for support through the grant 004-1193086-3036.
\end{acknowledgements}


%% REFERENCES

%% The reference list is compiled as follows:

\begin{thebibliography}{}

\bibitem[{{Acker and Leptoukh}(2007)}]{acker2007}
Acker, J.G. and G. Leptoukh (2007) Online Analysis Enhances Use of NASA Earth Science Data, {\it Eos, Trans. AGU,} {Vol. 88}, No. 2.
\bibitem[{{ALADIN International Team}(1997)}]{aladin}
ALADIN International Team (1997), The ALADIN project: Mesoscale modelling seen as a basic tool for weather forecasting and atmospheric research, {\it WMO Bull.,} {46}, 317-324.
\bibitem[{{Alpert et al.}(1990)}]{alpert}
Alpert, P., B. U. Neeman, and Y. Shay-El (1990), Intermonthly variability of cyclone tracks in the Mediterranean, {\it J. Climate,} {3}, 1474--1478.
\bibitem[{{Artegiani et al.}(1997)}]{artegiani}
Artegiani, A., D. Bregant, E. Paschini, N. Pinardi, F. Raicich, and A. Russo (1997), The Adriatic Sea general circulation. Part II: baroclinic circulation structure, {\it Journal of Physical Oceanography,} {27}, 1515--1532.
\bibitem[{{Baji\'{c} et al.}(2007)}]{bajic2007}
Baji\'{c}, A., S. Ivatek-\v{S}ahdan, and K. Horvath (2007), Spatial distribution of wind speed in Croatia obtained using the ALADIN model, {\it Cro. Met. J.,} {42}, 67--77.
\bibitem[{{Beg Paklar et al.}(2008)}]{beg}
Beg Paklar, G., \v{Z}agar, N., \v{Z}agar, M., Vellore, R., Kora\v{c}in, D., Poulain, P.-M., Orli\'{c}, M., Vilibi\'{c}, I., and Dadi\'{c}, V., Modeling the trajectories of satellite-tracked drifters in the Adriatic Sea during a summertime bora event, J. Geophys. Res., 113, C11S04, doi:10.1029/2007JC004536, 2008.
\bibitem[{{Bentamy and Croiz\'{e}-Fillon}(2012)}]{bentamy2012a}
Bentamy, A. and D. Croiz\'{e}-Fillon (2012), Gridded surface wind fields from Metop/ASCAT measurements. {\it Int. Jour. of Rem. Sens.}, {33}, 6, 1729--1754.
\bibitem[{{Bentamy et al.}(2012)}]{bentamy2012}
Bentamy, A., S.A. Grodsky, J.A. Carton, D. Croiz\'{e}-Fillon and B. Chapron (2012), Matching ASCAT and QuikSCAT winds, {\it Journal of Geophysical Research} {117}, C02011, 15 pp. doi:10.1029/2011JC007479
%\bibitem[{{Brankovi\'{c} et al.}(2007)}]{brankovic2007}
%Brankovi\'{c}, \v{C}., B. Matja\v{c}i\'{c}, S. Ivatek-\v{S}ahdan, and R. Buizza (2007), Dynamical downscaling of ECMWF EPS forecasts applied to cases of severe weather in Croatia, {\it ECMWF RD Technical Memorandum,} {507}, 38 pp.
\bibitem[{{Brankovi\'{c} et al.}(2008)}]{brankovic2008}
Brankovi\'{c}, \v{C}., B. Matja\v{c}i\'{c}, S. Ivatek-\v{S}ahdan, and R. Buizza (2008), Downscaling of ECMWF Ensemble Forecasts for Cases of Severe Weather: Ensemble Statistics and Cluster Analysis, {\it Mon. Wea. Rev.,} {136}, 3323--3342.
\bibitem[{{Brzovi\'{c}}(1999)}]{brzovic1999}
Brzovi\'{c}, N. (1999), Factors affecting the Adriatic cyclone and associated windstorms, {\it Contr. Atmos. Phys.,} {72}, 51--65.
\bibitem[{{Brzovi\'{c} and Strelec Mahovi\'{c}}(1999)}]{brzomah1999}
Brzovi\'{c}, N., and N. Strelec Mahovi\'{c} (1999), Cyclonic activity and severe jugo in the Adriatic, {\it Phys. Chem. Earth,} {24}, 653--657.
\bibitem[{{Burrage at al.}(2008)}]{burrage2008}
Burrage, D., J. Wesson, C. Martinez, C. P\'{e}rez, O. M\:{o}ller, and A. Piola (2008), Patos Lagoon outflow within the R\'{i}o de la Plata plume using an airborne salinity mapper: observing an embedded plume, {\it Continental Shelf Research,} {28}, 1625--1638.
\bibitem[{{Burrage at al.}(2009)}]{burrage2009}
Burrage, D.M., J.W. Book, and P.J. Martin (2009), Eddies and filaments of the Western Adriatic Current near Cape Gargano: Analysis and prediction, {\it Journal of Marine Systems,} {78}, 205--226.
\bibitem[{{Carr et al.}(2008)}]{carr}
Carr, S.D., X.J. Capet, J.C. McWilliams, J.T. Pennington, and F.P. 
Chavez, 2008: The influence of diel vertical migration on zooplankton 
transport and recruitment in an upwelling region: Estimates from a 
coupled behavioral-physical model. Fisheries Ocean., 17, 1-15.
\bibitem[{{Cassou and Terray}(2001)}]{cassou2001}
Cassou, C., and L. Terray (2001), Oceanic forcing of the wintertime low-frequency atmospheric variability in the north atlantic european sector: a study with the arpege model, {\it J. Climate,} {14}, 4266--4291, doi: http://dx.doi.org/10.1175/1520-0442(2001)014<4266:OFOTWL>2.0.CO;2
\bibitem[{{Cushman-Roisin at al.}(2001)}]{cushman2001}
Cushman-Roisin, B., M. Ga\v{c}i\'{c}, P.-M. Poulain, and A. Artegiani (Eds.) (2001), Physical Oceanography of the Adriatic Sea, {\it Kluwer Academic Publishers,} Dordrecht.
%\bibitem[{{Campinis at al.}(2000)}]{campini2000}
%Campinis, J., A. Genoves, A. Jansa, J. A. Guijarro, and C. Ramis (2000), A catalogue and a classification of surface cyclones for the Western Mediterranean, {\it Int. J. Climatol.,} {20}, 969--984.
%\bibitem[{{Catry at al.}(2007)}]{catry2007}
%Catry B., J.-F. Geleyn, M. Tudor, P. B\'{e}nard, and A. Trojakova (2007), Flux conservative thermodynamic equations in a mass-weighted framework, {\it Tellus, Ser. A,} {59}, pp 71--79
\bibitem[{{Dee et al.}(2011)}]{dee2011}
Dee, D. P., S.M. Uppala, A.J. Simmons, P. Berrisford, P. Poli, S. Kobayashi, U. Andrae, M.A. Balmaseda, G. Balsamo, P. Bauer, P. Bechtold, A.C.M. Beljaars, L. van de Berg, J. Bidlot, N. Bormann, C. Delsol, R. Dragani, M. Fuentes, A.J. Geer, L. Haimberger, S.B. Healy, H. Hersbach, E.V. H\'{o}lm, L. Isaksen, P. K\`{a}llberg, M. K\"{o}hler, M. Matricardi, A.P. McNally, B.M. Monge-Sanz, J.-J. Morcrette, B.-K. Park, C. Peubey, P. de Rosnay, C. Tavolato, J.-N. Th\'{e}paut, and F. Vitart (2011), The ERA-Interim reanalysis: configuration and performance of the data assimilation system. {\it Quart. J. Roy. Meteor. Soc.}, {137}, 553--597. doi: 10.1002/qj.828
\bibitem[{{D\"{o}\"{o}s et al.}(2011)}]{doos}
D\"{o}\"{o}s, K., Rupolo, V., Brodeau, L., Dispersion of surface drifters and model-simulated trajectories, Ocean Modelling, 39, 3-4, 301-310, 2011.
\bibitem[{{Dorman, et al.}(2006)}]{dorman}
Dorman, C., Book, J., Carniel, S., Cavaleri, L., Chiggiato, J.,  Doyle, J., Grbec, B.,  Haack, T.,  Janekovi\'{c}, I.,  Lee, C., Mala\v{c}i\v{c},  V., Orli\'{c}, M.,  Pullen, J., Russo, N., Paschini, E., Sclavo, M., Vilibi\'{c}, I., (2006),  Winter 2003 Marine Atmospheric Conditions and the Bora over the Northern Adriatic , {\it Journal of Geophysical Research}, {11 1} , C03S03, doi:10.1029/2005JC003134.
%\bibitem[{{Doswell at al.}(1996)}]{doswell1996}
%Doswell, C. A., H. E. Brooks, and R. A. Maddox (1996), Flash flood forecasting: An ingredients-based methodology, {\it Wea. Forecasting,} {11}, 560--581.
\bibitem[{{Geleyn}(1988)}]{geleyn1988}
Geleyn, J.-F. (1988), Interpolation of wind, temperature and humidity values from model levels to the height of measurement. {\it Tellus, Ser. A,} {40}, pp.347--351.
%\bibitem[{{Grbec et al.}(2007)}]{grbec2007}
%Grbec, B., I. Vilibi\'{c}, A. Baji\'{c}, M. Morovi\'{c}, G. Beg Paklar, F. Mati\'{c}, and V. Dadi\'{c} (2007), Response of the Adriatic Sea to the atmospheric anomaly in 2003, {\it Annales Geophysicae,} {25}, 835--846.
\bibitem[{{Grisogono and Belu\v{s}i\'{c}}(2009)}]{grisog2009}
Grisogono, B., and D. Belu\v{s}i\'{c} (2009), A review of recent advances in understanding the meso- and microscale properties of the severe Bora wind, {\it Tellus, Vol A,} {61}, 1--16.
\bibitem[{{Grubi\v{s}i\'{c}}(2004)}]{grubi2004}
Grubi\v{s}i\'{c}, V. (2004), Bora-driven potential vorticity banners over the Adriatic, {\it Quart. J. Roy. Meteor. Soc.,} {130}, 2571--2603.
%\bibitem[{{Holton}(1992)}]{holton1992}
%Holton, J.R. (1992), An Introduction to Dynamic Meteorology, Third Edition,{\it Academic Press,} 511 pp.
\bibitem[{{Horton et al.}(1997)}]{horton1997}
Horton, C., M. Clifford, J. Schmitz, and L.H. Kantha (1997), A real-time oceanographic nowcast/forecast systemfor the Mediterranean Sea, {\it Journal of Geophysical Research,} {102 (C11)}, 25123--25156.
\bibitem[{{Horvath et al.}(2011)}]{horvath2011}
Horvath, K., A. Baji\'{c}, and S. Ivatek-\v{S}ahdan (2011), Dynamical downscaling of wind speed in complex terrain prone to bora-type flows, {\it J. Appl. Meteor. Climatol.,} {50}, 1676--1691.
\bibitem[{{Horvath et al.}(2009)}]{horvath2009}
Horvath, K., S. Ivatek-\v{S}ahdan, B. Ivan\v{c}an-Picek, and V. Grubi\v{s}i\'{c} (2009), Evolution and structure of two severe cyclonic bora events: contrast between the northern and southern Adriatic, {\it Weather and forecasting,} {24}, 946--964.
\bibitem[{{Horvath et al.}(2008)}]{horvath2008}
Horvath, K., Y.-H. Lin, and B. Ivancan-Picek (2008), Classification of cyclone tracks over the Apennines and the Adriatic Sea, {\it Mon. Wea. Rev.,} {136}, 2210--2227.
\bibitem[{{Huffman et al.}(2007)}]{huffman2007}
Huffman, G.J., D.T. Bolvin, E.J. Nelkin, D.B. Wolff, R.F. Adler, G.Gu, Y.Hong, K.P. Bowman, E.F. Stocker (2007), The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-Global, Multiyear, Combined-Sensor Precipitation Estimates at Fine Scales {\it J. Hydrometeorol.}, {8}, 38--55.
%\bibitem[{{Ivatek-\v{S}ahdan and Ivan\v{c}an-Picek}(2006)}]{ivateks2006}
%Ivatek-\v{S}ahdan, S. and B. Ivan\v{c}an-Picek (2006), Effects of different initial and boundary conditions in ALADIN/HR simulations during MAP IOPs, {\it Meteorol. Z.} {15}, 187--197.
\bibitem[{{Ivatek-\v{S}ahdan and Tudor}(2004)}]{ivateks2004}
Ivatek-\v{S}ahdan, S., and M. Tudor (2004), Use of high-resolution dynamical adaptation in operational suite and research impact studies, {\it Meteorol Z} {13(2)}, 1--10.
\bibitem[{{Janekovi\'{c} et al.}(2010)}]{janekovic}
Janekovi\'{c}, I., M. A. Dutour Sikiri\'{c}, I. Toma\v{z}i\'{c}, and M. Kuzmi\'{c} (2010), Hindcasting the Adriatic Sea surface temperature and
salinity: A recent modeling experience. {\it Geofzika}, {27(2)}, 85--
100.
\bibitem[{{Jur\v{c}ec et al.}(1996)}]{jurcec1996}
Jur\v{c}ec, V., B. Ivan\v{c}an-Picek, V. Tuti\v{s}, and V. Vuki\v{c}evi\'{c} (1996), Severe Adriatic jugo wind, {\it Meteorol. Z.,} {N.F. 5}, 67--75.
\bibitem[{{Koren et al.}(2012)}]{koren2012}
Koren, I., O. Altaratz, L.A. Remer, G. Feingold, J. V. Martins and R.H. Heiblum (2012) Aerosol-induced intensification of rain from the tropics to the midlatitudes. {\it Nat. Geosci.,} {5(2)}, 118-122.
\bibitem[{{Kuzmi\'{c} et al.}(2006)}]{kuzmic}
Kuzmi\'{c}, M., I. Janekovi\'{c}, J.W. Book, P.J. Martin, and J.D. Doyle (2006), Modeling the northern Adriatic double-gyre response to intense bora wind: A revisit. {\it Journal of Geophysical Research - Oceans} {111}, C03S12. doi:10.1029/ 2005JC003377 (printed 112(C3), 2007).
\bibitem[{{Lionello et al.}(2006)}]{lionello2006}
Lionello, P., and Coauthors, 2006: Cyclones in the Mediterranean region: Climatology and effects on the environment, in Mediterranean Climate Variability, edited by P. Lionello, P. Malanotte-Rizzoli, and R. Boscolo, pp. 325--372, Elsevier, .
\bibitem[{{Liu and Weisberg}(2011)}]{liu}
Liu, Y., and Weisberg R. H., Evaluation of trajectory modeling in different dynamic regions using normalized cumulative Lagrangian separation, J. Geophys. Res., 116, C09013, doi:10.1029/2010JC006837, 2011.
\bibitem[{{Maga\v{s}}(2002)}]{magas2002}
Maga\v{s}, D. (2002) Nautral-geographic characteristics of the Bok Kotorska area as the basis of development. {\it Geoadria,} {7/1}, 51--81.
%\bibitem[{{Martin et al.}(2009)}]{martin2009}
%Martin, P.J., J.W. Book, D.M. Burrage, C.D. Rowley, and M. Tudor (2009), Comparison of model simulated and observed currents in the central Adriatic during DART, {\it Journal of Geophysical Research - Oceans} {114}, C01S05. doi:10.1029/2008JC004842.
\bibitem[{{Mesinger and Strickler}(1981)}]{mesing1981}
Mesinger, F., and R. F. Strickler (1981), Effect of mountains on Genoa cyclogeneses, {\it Journal of Meteorological Society of Japan,} {60}, 326--337.
\bibitem[{{Montenegro et al.}(2008)}]{monte2008}
Montenegro, A., C. Avis and A.J. Weaver (2008), Modelling the pre-historic arrival of the sweet potato in Polynesia. {\it Journal of Archaeological Science}, {35}, 355--367.
\bibitem[{{Orli\'{c} et al.}(1994)}]{orlic1994}
Orli\'{c}, M., M. Kuzmi\'{c}, and Z. Pasari\'{c} (1994), Response of the Adriatic Sea to the bora and sirocco forcing, {\it Cont. Shelf Res.,} {14}, 91--116.
\bibitem[{{Orli\'{c} et al.}(2007)}]{orlic2007}
Orli\'{c}, M., V. Dadi\'{c}, B. Grbec, N. Leder, A. Marki, F. Mati\'{c}, H. Mihanovi\'{c}, G. Beg Paklar, M. Pasari\'{c}, Z. Pasari\'{c}, and I. Vilibi\'{c} (2007), Wintertime buoyancy forcing, changing seawater properties and two different circulation systems produced in the Adriatic, {\it Journal of Geophysical Research - Oceans,} {112}, C03. doi:10.1029/2005JC003271.
\bibitem[{{Pasari\'{c} et al.}(2007)}]{pasaric2007}
Pasari\'{c}, Z., D. Belu\v{s}i\'{c}, and Z. Benceti\'{c} Klai\'{c} (2007), Orographic influences on the Adriatic sirocco wind, {\it Annales Geophysicae} {25}, 1263-1267.
%\bibitem[{{Platnick et al.}(2003)}]{platnick2003}
%Platnick, S. (2003) The MODIS cloud product
\bibitem[{{Poulain et al.}(2004)}]{poulain2004}
Poulain, P.-M., E. Mauri, and L. Ursella (2004), Unusual upwelling event and current reversal off the Italian Adriatic coast in summer 2003, {\it Geophys. Res. Lett.}, {31}, L05303, doi:10.1029/2003GL019121.
\bibitem[{{Pullen et al.}(2003)}]{pullen2003}
Pullen, J., J.D. Doyle, R. Hodur, A. Ogston, J.W. Book, H. Perkins and R. Signell (2003). Coupled ocean-atmosphere nested modelling of the Adriatic Sea during winter and spring 2001. {\it J. Geophys. Res.}, {108}C10, 3320, doi:10.1029/2003JC001780.
\bibitem[{{Remer et al.}(2008)}]{remer2008}
Remer, L.A., et al. (2008) Global aerosol climatology from MODIS satellite sensors. {\it J. Geophys. Res.}, {113}, D14S07.
\bibitem[{{Rodiriguez and Marcos}(2012)}]{crr}
Rodriguez, A. and C. Marcos (2012), Product user manual for the "convective rainfall rate" (CRR - PGE05 v3.1.1), {\it NWC SAF and AEMET}, {available from www.nwcsaf.org}, pp. 32.
\bibitem[{{Signell et al.} (2005)}]{signell}
Signell, R. P., S. Carniel, L. Cavaleri, J. Chiggiato, J. D. Doyle, J. Pullen,
and M. Sclavo (2005), Assessment of wind quality for oceanographic
modelling in semi-enclosed seas, {\it J. Mar. Syst.}, {53}, 217--233.
\bibitem[{{Shchepetkin and McWilliams}(2005)}]{schep2005}
Shchepetkin, A.F., and J.C. McWilliams (2005), The regional ocean modelling system: a splitexplicit, free-surface, topography-following-coordinate oceanic model, {\it Ocean Model.}, {9}, 347--404.
\bibitem[{{Smolarkiewicz and Margolin}(1998)}]{smolarkiewitz}
Smolarkiewicz, P. K. and Margolin, L. G. (1998), MPDATA: A finite-difference solver for geophysical flows. {\it J. Comput. Phys.}, {140}, 459--480.
\bibitem[{{Stane\v{s}i\'{c}}(2011)}]{stanesic2011}
Stane\v{s}i\'{c}, A. (2011), Assimilation system at DHMZ: development and first verification results. {\it Cro. Met. Jour.,} {44/45}, 3--17.
\bibitem[{{Termonia}(2008)}]{termonia2008}
Termonia, P. (2008), Scale-selective digital filter initialization. {\it Mon. Wea. Rev.}, {136}, pp. 5246--5255.
\bibitem[{{Thoss}(2012)}]{pc}
Thoss, A. (2012), Product user manual for SAFNWC/MSG "precipitating cloud" (PC - PGE04 v1.5), {\it NWC SAF and SMHI}, {available from www.nwcsaf.org}, pp. 18.
%\bibitem[{{Tompkins et al.}(2007)}]{tompk2007}
%Tompkins, A. M., K. Gierens and G. R\:{a}del (2007), Ice supersaturation in the ECMWF integrated forecast system. {\it Quart. J. Roy. Meteor. Soc.,} {133}, 53--63. doi: 10.1002/qj.14
\bibitem[{{Torres et al.}(2002)}]{torres2002}
Torres, O., R. Decae, J.P. Veefkind and G. de Leeuw (2002) OMI Aerosol Retrieval Algorythm, {\it Clouds, Aerosols and Surface UV Irradiance}, P. Stammes (ed.), {Vol III}, Version2.0.
\bibitem[{{Tudor and Ivatek-\v{S}ahdan}(2010)}]{tudor2010}
Tudor, M., and S. Ivatek-\v{S}ahdan (2010), The case study of bura of 1 and 3 February 2007, {\it Meteorol. Z.,} {19 (5)} 453--466.
\bibitem[{{Tudor and Termonia}(2010)}]{tudor2010b}
Tudor, M. and P. Termonia (2010), Alternative formulations for incorporating
lateral boundary data into limited area models, {\it Mon. Wea. Rev.}, {138}, pp. 2867--2882.
\bibitem[{{Veihelmann et al.}(2007)}]{veihel2007}
Veihelmann, B., P.F. Levelt, P. Stammes and J.P. Veefkind (2007), Simulation study of the aerosol information content in OMI spectral reflectance measurements, {\it Atmos. Chem. Phys.}, {7}, 3115--3127.
\bibitem[{{Wolf and Luksch}(1887)}]{wolf1887}
Wolf, J., and J. Luksch (1887), Bericht an die Seebeh\:{o}rde in Fiume \:{u}ber die Vorexpedition Nautilus im Sommer 1874 zum Zwecke physikalischer Untersuchungen des adriatischen Meeres, {\it I, Bericht. Fiume.} 55 pp.
\bibitem[{{Zaninovi\'{c} et al.}(2008)}]{zaninovic2008}
Zaninovi\'{c}, K., Gaji\'{c}-\v{C}apka, M., Per\v{c}ec Tadi\'{c}, M., Vu\v{c}eti\'{c}, M., Milkovi\'{c}, J., Baji\'{c}, A., Cindri\'{c}, K., Cvitan, L., Katu\v{s}in, Z., Kau\v{c}i\'{c}, D., Likso, T., Lon\v{c}ar, E., Lon\v{c}ar, \v{Z}., Mihajlovi\'{c}, D., Pand\v{z}i\'{c}, K., Patar\v{c}i\'{c}, M., Srnec, L., Vu\v{c}eti\'{c}, V., (2008) Climate atlas of Croatia: 1961. - 1990. : 1971. - 2000. . Zagreb : CMHS Monograph, pp. 200.
%\bibitem[{{Zavatarelli et al.}(2002)}]{zavat2002}
%Zavatarelli, M., N. Pinardi, V.H. Kourafalou, and A. Maggiore (2002), Diagnostic and prognostic model studies of the Adriatic Sea general circulation: seasonal variability, {\it Journal of Geophysical Research,} {107} (C1), 3004. doi:10.1029/2000JC000210.
%\bibitem[{{Zore}(1956)}]{zore1956}
%Zore, M. (1956), On the gradient currents in the Adriatic Sea, {\it Acta Adriatica,} {8 (6)}, 1--38.
%\bibitem[{{Zore-Armanda and Bone}(1987)}]{zore1987}
%Zore-Armanda, M., and M. Bone (1987), The effect of bottom topography on the current system of the open Adriatic Sea,{\it  Bolletino di Oceanologia Teorica ed Applicata,} {vol. 5}, pp. 3--18.

%\bibitem[AUTHOR(YEAR)]{LABEL}
%REFERENCE 1

%\bibitem[AUTHOR(YEAR)]{LABEL}
%REFERENCE 2

\end{thebibliography}

%% Since the Copernicus LaTeX package includes the BibTeX style file copernicus.bst,
%% authors experienced with BibTeX only have to include the following two lines:
%%
%% \bibliographystyle{copernicus}
%% \bibliography{example.bib}
%%
%% URLs and DOIs can be entered in your BibTeX file as:
%%
%% URL = {http://www.xyz.org/~jones/idx_g.htm}
%% DOI = {10.5194/xyz}


%% LITERATURE CITATIONS
%%
%% command                        & example result
%% \citet{jones90}|               & Jones et al. (1990)
%% \citep{jones90}|               & (Jones et al., 1990)
%% \citep{jones90,jones93}|       & (Jones et al., 1990, 1993)
%% \citep[p.~32]{jones90}|        & (Jones et al., 1990, p.~32)
%% \citep[e.g.,][]{jones90}|      & (e.g., Jones et al., 1990)
%% \citep[e.g.,][p.~32]{jones90}| & (e.g., Jones et al., 1990, p.~32)
%% \citeauthor{jones90}|          & Jones et al.
%% \citeyear{jones90}|            & 1990


%% FIGURES

%% ONE-COLUMN FIGURES

%f
%\begin{figure}[t]
%\includegraphics[width=8.3cm]{FILE NAME}
%\caption{TEXT}
%\end{figure}

%% TWO-COLUMN FIGURES

%f
%\begin{figure*}[t]
%\includegraphics[width=12cm]{FILE NAME}
%\caption{TEXT}
%\end{figure*}

 \begin{figure}
 \includegraphics[width=12cm]{orogpeljzoom.eps}
 \caption{The South Adriatic region with locations of meteorological stations used in this study. The shaded background represents terrain height and sea depth in kilometers according to the scale on the right side. The "Bojana" mark shows position of the point where Bojana river enters the Adriatic Sea (this is also the border between Montenegro and Albania on the coast).}
 \label{terrain}
 \end{figure}
%
 \begin{figure}
 \includegraphics[width=12cm]{prec2.eps}
 \caption{Measured accumulated 24 hourly precipitation on rain guagues in southeastern Croatia (smaller sybols) and Montenegro (larger symbols) during November 2012. Precipitation shown for a certain day is measured at 6 UTC accumulated from the previous 24 hours. Only stations close to Adriatic are shown.}
 \label{raingauge}
 \end{figure}
%
 \begin{figure}
 \includegraphics[angle=270,width=12cm]{erai_850uvgt20101107_12.eps}
 \caption{ERA Interim 850 hPa wind (colour of the vectors shows wind speed in m/s as on the colour bar below), geopotentiel (blue isolines) and temperature (red isolines) with measured aerosol optical thickness at 12 UTC 7 November 2010 from MODIS (circles) and OMI (squares), AOT is shown in colour according to the colour bar on the right.}
 \label{erai850}
 \end{figure}
%
 \begin{figure}
 \includegraphics[angle=270,width=12cm]{prog_meas24h09NOV2010.eps}
 \includegraphics[angle=270,width=12cm]{prog_meas24h10NOV2010.eps}
 \caption{Measured accumulated 24 hourly precipitation on rain guagues in Croatia, Montenegro and Macedonia (circles), TRMM rainfall data (squares) and 8km ALADIN forecast data (shaded), the precipitation is accumulated for the period from 06 UTC on 8 until 06 UTC on 9 (top) and from 06 UTC on 9 until 06 UTC on 10 (bottom) November 2010. The correcponding rainguage data is shown in Fig~\ref{raingauge} as measured on 9 and 10 November 2010.}
 \label{rain1}
 \end{figure}
%
 \begin{figure}
 \includegraphics[angle=270,width=12cm]{erai_uvnp20101116_12.eps}
 \caption{ERA Interim 10m wind (colour of the vectors shows wind speed in m/s as on the colour bar below), mean sea level pressure (blue) and cloudiness (shades of gray) for 12 UTC on 16 November 2010.}
 \label{eraiuvnp}
 \end{figure}
%
 \begin{figure}
 \includegraphics[angle=270,width=12cm]{vjea_hrp.eps} 
 \caption{Measured wind speed (dark blue, scale on the left) and mean sea level pressure (red) reduced by 980 hPa (to fit the scale on the left) and wind direction (light blue in degrees from north clockwise, scale on the right) for November 2010 (Sv. Jure is marked as Biokovo in Fig~\ref{terrain}).}
 \label{wind}
 \end{figure}
%
 \begin{figure}
 \includegraphics[width=6.0cm]{wind20101116.eps}
 \includegraphics[width=6.0cm]{wind20101117.eps}
 \includegraphics[width=6.0cm]{wind20101116a.eps}
 \includegraphics[width=6.0cm]{wind20101117a.eps}
 \caption{Forecast 10m wind in 8km (top) and 2km (bottom) resolution and measured wind speed and direction (arrows) from MetOp ASCAT data (above the sea surface), SYNOP and automatic stations for 12 UTC on 16 (left) and 12 UTC on 17 (right) November 2010. Colour of the vectors shows wind speed in m/s as on the colour bars. Model data are shown as thin vectors on a denser grid in the background, measured wind is shown as thicker vectors on the location of measurement. }
 \label{wind16}
 \end{figure}
%
 \begin{figure}
 \includegraphics[width=12cm]{cgrivers.eps}
 \caption{Measured water levels of rivers in Montenegro and Skadar lake in October and November 2010 and climatological river run-off of Albanian rivers (embeded figure). Moraca and Crnojevica rivers are tributaries to the Skadar lake, Bojana river takes the outflow from the Skadar lake and flows into the Adriatic sea at the position marked as Bojana in Fig~\ref{terrain}.}
 \label{rivers}
 \end{figure}
%
 \begin{figure}
 \includegraphics[width=8.3cm]{roms_fig.eps}
 \caption{Surface currents (vectors) and sea surface temperature (shaded background) from ROMS for 00 UTC on 11 (a), 15 (b) and 19 (c) November 2010.}
 \label{roms}
 \end{figure}
%
 \begin{figure}
 \includegraphics[width=12cm]{trajs.eps}
 \caption{Trajectories of drifters released at 12 UTC 19 (a) and 12 UTC 26 October (b), 00 (c) and 12 (d) UTC 8, 00 (e) and 12 (f) UTC 9, 00 (g) and 12 (h) UTC 10 and 00 UTC 11 (i) November 2010. Trajectories initiated off different parts of a rectangle in southeast Adriatic are plotted in different colours.}
 \label{trajs}
 \end{figure}
%


%% TABLES
%%
%% The different columns must be seperated with a & command and should
%% end with \\ to identify the column brake.

%% ONE-COLUMN TABLE

%t
%\begin{table}[t]
%\caption{TEXT}
%\begin{tabular}{column = lcr}
%\tophline

%\middlehline

%\bottomhline
%\end{tabular}
%\belowtable{} % Table Footnotes
%\end{table}

%% TWO-COLUMN TABLE

%t
%\begin{table*}[t]
%\caption{TEXT}
%\begin{tabular}{column = lcr}
%\tophline

%\middlehline

%\bottomhline
%\end{tabular}
%\belowtable{} % Table Footnotes
%\end{table*}


%% NUMBERING OF FIGURES AND TABLES
%%
%% If figures and tables must be numbered 1a, 1b, etc. the following command
%% should be inserted before the begin{} command.

\addtocounter{figure}{-1}\renewcommand{\thefigure}{\arabic{figure}a}


%% MATHEMATICAL EXPRESSIONS

%% All papers typeset by Copernicus Publications follow the math typesetting regulations
%% given by the IUPAC Green Book (IUPAC: Quantities, Units and Symbols in Physical Chemistry,
%% 2nd Edn., Blackwell Science, available at: http://old.iupac.org/publications/books/gbook/green_book_2ed.pdf, 1993).
%%
%% Physical quantities/variables are typeset in italic font (t for time, T for Temperature)
%% Indices which are not defined are typeset in italic font (x, y, z, a, b, c)
%% Items/objects which are defined are typeset in roman font (Car A, Car B)
%% Descriptions/specifications which are defined by itself are typeset in roman font (abs, rel, ref, tot, net, ice)
%% Abbreviations from 2 letters are typeset in roman font (RH, LAI)
%% Vectors are identified in bold italic font using \vec{x}
%% Matrices are identified in bold roman font
%% Multiplication signs are typeset using the LaTeX commands \times (for vector products, grids, and exponential notations) or \cdot
%% The character * should not be applied as mutliplication sign


%% EQUATIONS

%% Single-row equation

%\begin{equation}

%\end{equation}

%% Multiline equation

%\begin{align}
%& 3 + 5 = 8\\
%& 3 + 5 = 8\\
%& 3 + 5 = 8
%\end{align}


%% MATRICES

%\begin{matrix}
%x & y & z\\
%x & y & z\\
%x & y & z\\
%\end{matrix}


%% ALGORITHM

%\begin{algorithm}
%\caption{�}
%\label{a1}
%\begin{algorithmic}
%�
%\end{algorithmic}
%\end{algorithm}


%% CHEMICAL FORMULAS AND REACTIONS

%% For formulas embedded in the text, please use \chem{}

%% The reaction environment creates labels including the letter R, i.e. (R1), (R2), etc.

%\begin{reaction}
%% \rightarrow should be used for normal (one-way) chemical reactions
%% \rightleftharpoons should be used for equilibria
%% \leftrightarrow should be used for resonance structures
%\end{reaction}


%% PHYSICAL UNITS
%%
%% Please use \unit{} and apply the exponential notation


\end{document}
