The Flow Generator Relations for Water Renewal through the Flushing Culverts in Marinas

: An understanding of the relationship between the ﬂow generators and harbor design could result in an acceptable seawater renewal into an enclosed basin. The ﬂushing culverts located in the breakwater are the most cost-effective technique for improving water renewal in port engineering. The aim of the paper was to identify and quantify the relation between ﬂow generators and culvert ﬂow to propose an optimal position for the ﬂushing culverts. Mentioned relations are important for designing and positioning the ﬂushing culverts in marinas to achieve optimal sea exchange in the marina basin. Field measurements in the marina Opatija were performed and the data served for calibrating a 3D circulation model and a wave generation/propagation model. The absolute volume of the ﬂow through the culverts caused by wind waves has the more important contribution to the water renewal in the marina basin than the coastal springs have. The optimal ﬂow rates through the culverts are performed by positioning the ﬂushing culvert axis at the mean sea level in accordance with the main incident direction of the propagating waves. The model the phenomena of wave generation due to wind action, nonlinear wave interaction, refraction as well as the dissipative process caused by surface white-capping. The multi sequential Euler explicit used for the discretization of the temporal term and the ‘upwind’ numerical scheme for the convective flux term.


Introduction
The fraction of people inhabiting the seashore is projected to rise from the present 60% to 75% by the year 2025, and the population as a whole is also anticipated to grow to 11.1 billion people by the year 2100 [1,2]. With the growth of the inhabitants and metropolitan areas near the shore, an anthropogenic pressure to the marine environment is also on the rise. Coastal infrastructure scattered along the shoreline are often pivotal points where the water quality is under biggest threat from anthropogenic activities [3]. Communities use ports in the vicinity of cities to trade goods, hence during cargo manipulation unanticipated accidents may well occur [4]. In the sea, near densely populated areas, further possibly harmful pollutants could accumulate as a consequence of human activity. Harmful bacteria like E. coli, heavy metals like copper and pesticides are likely to be deposited in the seabed sediment which poses a threat to both human and animal life at sea, especially if dredging or seabed erosion resuspends them again [5][6][7][8]. Stagnant seawater regularly is characterized by a reduction in dissolved oxygen, which favors the growth of harmful algal blooms and afterward could further compromise the marine ecosystem [9,10]. Pollution concentrations should be continuously monitored at sensitive coastal locations, to prevent the undesirable effects of various pollutants. Harmful effects of high pollutant concentrations within a semi-enclosed basin can be avoided if a sufficient water renewal rate is upheld or the sources of pollution are prevented [11]. Water renewal in semi-enclosed basins has proved to be a critical aspect influencing the pollutant concentrations and therefore life in a marine ecosystem [12]. The water renewal rate in human-made or naturally enclosed basins is conditioned by seawater circulation generators (e.g., wind, waves, tidal oscillations, inner freshwater

Field Measurements
Field measurements were performed inside and in the vicinity of ACI Marina Opatija (45°19′/14°17.7′) in northern Croatia, near the city of Rijeka ( Figure 2). The marina is linked to the Adriatic Sea with a 30 m wide inlet (marina entrance). In the primary breakwater, eight parallel flushing culverts (circular cross sections with a 1 m diameter) are embedded to improve water renewal inside the basin. The mean marina depth is approximately 5 m, with the depth decreasing in the vicinity of the breakwater. Detailed bathymetry surveillance was conducted in the marine area, to establish a sea circulation numerical model. Oceanographic parameters (current velocity and direction, salinity, sea temperature, significant wave height, sea level) were monitored during the winter (15 February 2017-

Field Measurements
Field measurements were performed inside and in the vicinity of ACI Marina Opatija (45 • 19 /14 • 17.7 ) in northern Croatia, near the city of Rijeka ( Figure 2). The marina is linked to the Adriatic Sea with a 30 m wide inlet (marina entrance). In the primary breakwater, eight parallel flushing culverts (circular cross sections with a 1 m diameter) are embedded to improve water renewal inside the basin. The mean marina depth is approximately 5 m, with the depth decreasing in the vicinity of the breakwater.

Field Measurements
Field measurements were performed inside and in the vicinity of ACI Marina Opatija (45°19′/14°17.7′) in northern Croatia, near the city of Rijeka ( Figure 2). The marina is linked to the Adriatic Sea with a 30 m wide inlet (marina entrance). In the primary breakwater, eight parallel flushing culverts (circular cross sections with a 1 m diameter) are embedded to improve water renewal inside the basin. The mean marina depth is approximately 5 m, with the depth decreasing in the vicinity of the breakwater. Detailed bathymetry surveillance was conducted in the marine area, to establish a sea circulation numerical model. Oceanographic parameters (current velocity and direction, salinity, sea temperature, significant wave height, sea level) were monitored during the winter (15 February 2017- Detailed bathymetry surveillance was conducted in the marine area, to establish a sea circulation numerical model. Oceanographic parameters (current velocity and direction, salinity, sea temperature, significant wave height, sea level) were monitored during the winter ( marina Opatija basin ( Figure 3). Obtained data sets were utilized in order to define the initial conditions, boundary conditions and to verify the numerical model results. Currents and the sea surface elevation (Figure 4) were monitored continuously during the winter and summer measurement periods with the use of five Acoustic Doppler Current Profilers (ADCP at positions 1-5) (Figure 3). For the duration of the summer measurement period, all the devices were positioned at the same locations as during the winter except ADCP 3, which was repositioned next to ADCP 4 at the marina entrance in an attempt to estimate the flow rate through the entrance inlet. The sea temperature (T) and salinity (S) were measured throughout the vertical sea column at sites 1-4 nearby the marina basin ( Figure 3), during measuring device deployment and recovery. Wave conditions (significant wave heights H S , peak period T P and incident wave direction, Figure 5) were obtained using ADCPs at station 1 and 2 ( Figure 3). The PCM device (Pro Flow Meter for portable ultrasonic flow measurements) was located at the entrance of one flushing culvert (station 2, Figure 3) and used for culvert-flow rate measurements. A positive flow rate through the flushing culvert was assigned to inflows into the marina and negative to outflows from the marina. A CTD device (Conductivity/Temperature/Depth, Teledyne Marine group, Poway, CA, USA) was fixed at site 2, next to the flushing culvert, in order to measure temperature and salinity during measurement periods. Data regarding wind velocities were obtained from an anemometer (positioned on the lighthouse at the head of the breakwater) 10 m above the ground (Figures 3 and 5). Sampling resolutions of deployed instruments are shown in Table 1. The data on precipitation ( Figure 4) was obtained from a local meteorological station located near the marina (Rijeka station, φ = 45 • 20 , λ = 14 • 27 ).  20 March 2017) and summer (4 July 2017-29 August 2017) measurement periods at several locations at the ACI marina Opatija basin ( Figure 3). Obtained data sets were utilized in order to define the initial conditions, boundary conditions and to verify the numerical model results. Currents and the sea surface elevation ( Figure 4) were monitored continuously during the winter and summer measurement periods with the use of five Acoustic Doppler Current Profilers (ADCP at positions 1-5) (Figure 3). For the duration of the summer measurement period, all the devices were positioned at the same locations as during the winter except ADCP 3, which was repositioned next to ADCP 4 at the marina entrance in an attempt to estimate the flow rate through the entrance inlet. The sea temperature (T) and salinity (S) were measured throughout the vertical sea column at sites 1-4 nearby the marina basin ( Figure 3), during measuring device deployment and recovery. Wave conditions (significant wave heights HS, peak period TP and incident wave direction, Figure 5) were obtained using ADCPs at station 1 and 2 ( Figure 3). The PCM device (Pro Flow Meter for portable ultrasonic flow measurements) was located at the entrance of one flushing culvert (station 2, Figure 3) and used for culvert-flow rate measurements. A positive flow rate through the flushing culvert was assigned to inflows into the marina and negative to outflows from the marina. A CTD device (Conductivity/Temperature/Depth, Teledyne Marine group, Poway, CA, USA) was fixed at site 2, next to the flushing culvert, in order to measure temperature and salinity during measurement periods. Data regarding wind velocities were obtained from an anemometer (positioned on the lighthouse at the head of the breakwater) 10 m above the ground (Figures 3 and 5). Sampling resolutions of deployed instruments are shown in Table 1. The data on precipitation ( Figure 4) was obtained from a local meteorological station located near the marina (Rijeka station, ϕ = 45°20′, λ = 14°27′).    The incident wind waves approaching the marina in a band from 90° (E) to 180° (S) provide the most intense culvert inflow (marked by red line in Figure 6). During the winter and summer measurement periods, noteworthy wave situations were documented (Tables 2 and 3) assuming the criteria: HS ≥ 0.3 m, TP ≥ 2.8 s and incident wave direction in the band from 90° (E) to 180° (S), measured at the ADCP station 1. A total of seven wave situations were acknowledged during the winter period (Table 2) and 10 in the course of the summer measurement periods ( Table 3). The starting date of each proper situation, its duration, corresponding average significant wave height, peak period and maximum significant wave height are listed in Tables 2 and 3. The relevant wave conditions are also indicated in Figure 6. In the course of the winter, waves ranging from 100° to 174° caused an increased culvert inflow +QPCM (+QPCM denotes inflow-seawater entering into the marina basin Figure 3).  The incident wind waves approaching the marina in a band from 90° (E) to 180° (S) provide the most intense culvert inflow (marked by red line in Figure 6). During the winter and summer measurement periods, noteworthy wave situations were documented (Tables 2 and 3) assuming the criteria: HS ≥ 0.3 m, TP ≥ 2.8 s and incident wave direction in the band from 90° (E) to 180° (S), measured at the ADCP station 1. A total of seven wave situations were acknowledged during the winter period (Table 2) and 10 in the course of the summer measurement periods ( Table 3). The starting date of each proper situation, its duration, corresponding average significant wave height, peak period and maximum significant wave height are listed in Tables 2 and 3. The relevant wave conditions are also indicated in Figure 6. In the course of the winter, waves ranging from 100° to 174° caused an increased culvert inflow +QPCM (+QPCM denotes inflow-seawater entering into the marina basin Figure 3). The incident wind waves approaching the marina in a band from 90 • (E) to 180 • (S) provide the most intense culvert inflow (marked by red line in Figure 6). During the winter and summer measurement periods, noteworthy wave situations were documented (Tables 2 and 3) assuming the criteria: H S ≥ 0.3 m, T P ≥ 2.8 s and incident wave direction in the band from 90 • (E) to 180 • (S), measured at the ADCP station 1. A total of seven wave situations were acknowledged during the winter period (Table 2) and 10 in the course of the summer measurement periods ( Table 3). The starting date of each proper situation, its duration, corresponding average significant wave height, peak period and maximum significant wave height are listed in Tables 2 and 3. The relevant wave conditions are also indicated in Figure 6. In the course of the winter, waves ranging from 100 • to 174 • caused an increased culvert inflow +Q PCM (+Q PCM denotes inflow-seawater entering into the marina basin Figure 3).    Table 3. Start date of each relevant wave situation (adopted criteria HS ≥ 0.3 m; TP ≥ 2.8 s; incident waves ranging from 90° to 180°), its duration, corresponding average significant wave height, peak period, maximum significant wave height and as wave direction recorded at ADCP 1 for the summer period.    Northeastern wind generates wave fields parallel to the main breakwater with lower significant wave heights and periods ( Figure 5) (situation: 7 March and 10 March, Figure 5) wherein a culvert inflow is absent. The winds could cause currents parallel to the main breakwater and in consequence raise the sea surface elevation in the SW part of the marina, where culverts are constructed, so greater outflow could occur (situation on 10 March, marked by magenta line- Figure 6). Contrary to expectations, during a period of low wind velocities in the summer period (from 2 August to 4 August marked by yellow circle- Figure 6), a dominant culvert inflow was observed. It is probably a result of the specific flow pattern in the bay of Rijeka that influenced the unforeseen inflow.

Summer Measurement Campaign Deepwater Conditions Duration HS-AV TP-AV HS-MAX dirAV
In the Opatija area, i.e., the Ičići area, a large number of coastal springs have been documented ( Figure 7). This phenomenon is a consequence of the characteristic karstic geological condition with the occurrence of surface and subsurface springs in the observed area [27]. According to the hydrological observations at one of the most dominant coastal springs Kristal, during the dry summer periods, a discharge of 500 L/s was documented.  Figure 5) wherein a culvert inflow is absent. The winds could cause currents parallel to the main breakwater and in consequence raise the sea surface elevation in the SW part of the marina, where culverts are constructed, so greater outflow could occur (situation on 10 March, marked by magenta line- Figure 6). Contrary to expectations, during a period of low wind velocities in the summer period (from 2 August to 4 August marked by yellow circle- Figure 6), a dominant culvert inflow was observed. It is probably a result of the specific flow pattern in the bay of Rijeka that influenced the unforeseen inflow.
In the Opatija area, i.e., the Ičići area, a large number of coastal springs have been documented ( Figure 7). This phenomenon is a consequence of the characteristic karstic geological condition with the occurrence of surface and subsurface springs in the observed area [27]. According to the hydrological observations at one of the most dominant coastal springs Kristal, during the dry summer periods, a discharge of 500 L/s was documented. During the winter measurement period, coastal spring activities were observed along the marina coastline ( Figure 3). The coastal spring intensity is unknown, but visual inspection has determined its activation in days after precipitation occurrence. Further examinations showed that the spring intensity correlates with precipitation intensity. During the winter period, PCM device recorded a dominant outflow through the flushing culvert in the days after precipitation at 10 mm/day (marked by green arrows- Figure 6). The most intense precipitation intensity during field measurements was registered on the 25 February (38.6 mm/day, Figure 4) when the PCM device recorded a maximum culvert outflow of −0.112 m 3 /s. In the summer measurement period, the effect of spring activity on the culvert outflow is not as dominant as during the winter. It is evident that after the rain on 27 July (27 mm/day), a culvert outflow is recorded, but in a reduced magnitude in comparison to the winter observations. The summer period is characterized by higher temperatures and higher evapotranspiration rates, so much more precipitation is needed for coastal spring full activation. Measured precipitation at the meteorological station Rijeka does not necessarily have to cause the coastal spring outflow around the marina Opatija, since the movement of the air masses, which bring precipitation, is not always in the same direction. During the winter measurement period, coastal spring activities were observed along the marina coastline ( Figure 3). The coastal spring intensity is unknown, but visual inspection has determined its activation in days after precipitation occurrence. Further examinations showed that the spring intensity correlates with precipitation intensity. During the winter period, PCM device recorded a dominant outflow through the flushing culvert in the days after precipitation at 10 mm/day (marked by green arrows- Figure 6). The most intense precipitation intensity during field measurements was registered on the 25 February (38.6 mm/day, Figure 4) when the PCM device recorded a maximum culvert outflow of −0.112 m 3 /s. In the summer measurement period, the effect of spring activity on the culvert outflow is not as dominant as during the winter. It is evident that after the rain on 27 July (27 mm/day), a culvert outflow is recorded, but in a reduced magnitude in comparison to the winter observations. The summer period is characterized by higher temperatures and higher evapotranspiration rates, so much more precipitation is needed for coastal spring full activation.

Sea Circulation and Wave Generation Numerical Model
Measured precipitation at the meteorological station Rijeka does not necessarily have to cause the coastal spring outflow around the marina Opatija, since the movement of the air masses, which bring precipitation, is not always in the same direction.

Sea Circulation Model Setup
A three-dimensional Mike 3 FM model [28] was used in order to study the area of the marina Opatija. Its hydrodynamics are solved by utilizing the 3D RANS equations using the Boussinesq and hydrostatic approximations. The numerical model divides the domain in the vertical direction utilizing the standard sigma coordinate approach [29]. Governing equations are resolved within a finite volume grid, based on a discretization of the continuum with non-overlapping elements [30]. An approximate Riemann solver [31] is employed to compute the convective terms in the governing equations, ensuring numerical stability even in cases with steep gradients. The Smagorinsky scheme [32] and k-ε models [33] are the turbulence models employed in the horizontal and vertical directions, respectively.
Simulations with the Mike 3 FM model were run using the following parameter values: minimum time step of external mode ∆t = 0.1 s, maximum time step of internal mode ∆t = 30 s, with a critical threshold CFL of 0.8. dispersion coefficients for the temperature and salinity scalar fields were specified with a proportionality factor of 0.85 in the horizontal 0.09 in the vertical in accordance to the scaled eddy viscosity. Roughness and Smagorinsky coefficients were set as spatially and temporally constant at 0.01 and 0.2, respectively. The wind stress coefficient which models the influence of the wind on the sea surface was defined at 0.0013 [34]. In Dalton's law, 0.5 was set for the wind coefficient and 0.9 for the evaporation coefficient. The heat flux absorption profile in the short-wave radiation is described by a modified version of Beer's law. In the surface layer, the energy absorption constant equals 0.15, and light decay coefficient corresponds to 0.4 in the vertical direction. The hourly time series of air temperature, air humidity and cloud clearness from the nearby meteorological station Rijeka were used.
Numerical At the open boundaries of the numerical model domain (B-1 and B-2) the sea temperature and vertical salinity profiles were set. These datasets were obtained with the mobile CTD device at station 1 (Figures 3 and 4) during the deployment and retrieval of measurement devices. The initial conditions were defined with current velocity values set to 0 m/s in all three directions throughout the whole domain. For the initial conditions, a homogeneous scalar field of sea temperature and salinity in the horizontal direction has been set in agreement to the measured vertical profiles with the CTD at station 1 (Figures 3 and 4).
Coastal springs were implemented into the model utilizing 13 concentrated water sources along the marina coastline ( Figure 3). The source inflow was defined with a constant temperature value of 10 • C and salinity of 0 PSU. The source inflow dynamics depend on precipitation occurrence in the area (Figure 4). A flow rate of one spring source can be calculated using Q SOURCE = k·P, where P is a value of measured precipitation and k is a proportionality coefficient which depends on precipitation and other geophysical factors like soil type and temperature. For the current examination, only precipitation was taken into consideration to define k. Coefficient k is determined through model calibration ( Figure 9) and it has a different value for the winter and summer period because of mentioned climatological conditions. Exponential correlation was determined between precipitation P and coefficient k with the R 2 = 0.7 for winter and R 2 = 0.85 for the summer period ( Figure 9). Each day 25% of the total water volume created by precipitation was pushed into the domain through the sources at the seawall. In the course of 4 days, the source provided freshwater into the basin after which the freshwater sources subsided. m in total length. The culvert top edge was positioned at the mean sea level (Figure 3). The breakwater is included as a non-permeable structure (no flow land boundary) having zero wave reflection. At the open boundary B-1 ( Figure 8) velocities measured at ADCP station 1 are imposed to get resemblance between measured and modelled values. The measured surface elevation at ADCP 1 is used as forcing at the open boundary B-2 (Figure 8). At the boundary B-3 the normal velocity is set to zero.
At the open boundaries of the numerical model domain (B-1 and B-2) the sea temperature and vertical salinity profiles were set. These datasets were obtained with the mobile CTD device at station 1 (Figures 3 and 4) during the deployment and retrieval of measurement devices. The initial conditions were defined with current velocity values set to 0 m/s in all three directions throughout the whole domain. For the initial conditions, a homogeneous scalar field of sea temperature and salinity in the horizontal direction has been set in agreement to the measured vertical profiles with the CTD at station 1 (Figures 3 and 4).   Coastal springs were implemented into the model utilizing 13 concentrated water sources along the marina coastline ( Figure 3). The source inflow was defined with a constant temperature value of 10 °C and salinity of 0 PSU. The source inflow dynamics depend on precipitation occurrence in the area (Figure 4). A flow rate of one spring source can be calculated using QSOURCE = k·P, where P is a value of measured precipitation and k is a proportionality coefficient which depends on precipitation and other geophysical factors like soil type and temperature. For the current examination, only precipitation was taken into consideration to define k. Coefficient k is determined through model calibration ( Figure 9) and it has a different value for the winter and summer period because of mentioned climatological conditions. Exponential correlation was determined between precipitation P and coefficient k with the R 2 = 0.7 for winter and R 2 = 0.85 for the summer period ( Figure 9). Each day 25% of the total water volume created by precipitation was pushed into the domain through the sources at the seawall. In the course of 4 days, the source provided freshwater into the basin after which the freshwater sources subsided. The pipe culvert is implemented in the model as an explicit external routine based upon the difference of sea surface levels between the leeward and seaward side of the breakwater at the culvert position. The culvert is defined in the model as a connected sink-source function inside the domain. Depending on the sea surface difference, the function employs the free surface flow or pressurized flow through the culvert, utilizing the Manning equation in both cases (Manning coefficient n = 0.013 s/m 1/3 ). The head loss factors determine the energy (Head) loss over the culvert. The total head loss ΔH trough culvert is given by: where Q is discharge, ζ1 is the entrance or contraction loss, ζ2 is the outlet or expansion loss, ζf is the friction loss calculated using Manning equation, ζb is the bend loss coefficient. The suffix '1' and '2' represent the area on the inflow and outflow side of the structure respectively, and 's' the area in the structure itself [28]. Each head loss factor is defined by a suffix where '+' indicates positive flow direction and the suffix '−' indicates negative flow direction [28]. The inflow and outflow (positive and negative) head loss factors are constant during a simulation (positive inflow coefficient: ζIN + = 0.9, Figure 9. Correlation between the proportionality coefficient k and precipitation P (mm) used for source inflow rate calculation.
The pipe culvert is implemented in the model as an explicit external routine based upon the difference of sea surface levels between the leeward and seaward side of the breakwater at the culvert position. The culvert is defined in the model as a connected sink-source function inside the domain. Depending on the sea surface difference, the function employs the free surface flow or pressurized flow through the culvert, utilizing the Manning equation in both cases (Manning coefficient n = 0.013 s/m 1/3 ). The head loss factors determine the energy (Head) loss over the culvert. The total head loss ∆H trough culvert is given by: (1) where Q is discharge, ζ 1 is the entrance or contraction loss, ζ 2 is the outlet or expansion loss, ζ f is the friction loss calculated using Manning equation, ζ b is the bend loss coefficient. The suffix '1' and '2' represent the area on the inflow and outflow side of the structure respectively, and 's' the area in the structure itself [28]. Each head loss factor is defined by a suffix where '+' indicates positive flow direction and the suffix '−' indicates negative flow direction [28].  (Tables 2 and 3) the head loss factor was variable. The head loss factor values were adopted according to the integrated culvert discharge measured by the PCM device during each wave situation. The head loss factors were varied until the similarity between integrated measured and modelled discharges through the culvert were obtained for each wave situation.

Sea Circulation Model Results
Numerical analyses of sea currents and water mass exchange between the marina basin and the surrounding sea were carried out. A comparison of the measured and modeled hourly averaged current velocities at the position ADCP 1 (Figure 3) is shown in Figures 10 and 11. To evaluate the model performances normalized root-mean-square error (NRMSE) [35] and Nash-Sutcliffe efficiency index (NSE) [36] of the measured and modelled velocity values were used ( Table 4). The value of perfect fit for NRMSE is 0 and for NSE is 1 [35]. Obtained values for NRMSE and NSE indicate better resemblance between measured and modelled velocities for the summer than for the winter period. In the summer months, there were less intense transitional conditions of the atmosphere, especially the wind, during which the greater discrepancy between measured and modeled velocities is noticed.
Baroclinic flow is recognized in the marina entrance area. NW, NE and E winds caused the surface layer inflow, while winds of S, SW and W direction caused surface layer outflow. Figure 12 shows the hourly averaged current fields in characteristic situations during the occurrence of baroclinic compensation currents. N and SE winds coincide with the rainfall occurrence and the coastal spring activation. In the winter period two similar situations were detected during which N wind, contrary to expected, caused outflow current in the marina entrance area (on 7 March 15:00-21:00 h with average speed of 3.97 m/s, and on 8 March 5:00 h-8:00 h with average speed of 5 m/s). In the mentioned situations coastal springs were active and had more influence on surface layer flow than wind had.

Sea Circulation Model Results
Numerical analyses of sea currents and water mass exchange between the marina basin and the surrounding sea were carried out. A comparison of the measured and modeled hourly averaged current velocities at the position ADCP 1 (Figure 3) is shown in Figures 10 and 11. To evaluate the model performances normalized root-mean-square error (NRMSE) [35] and Nash-Sutcliffe efficiency index (NSE) [36] of the measured and modelled velocity values were used ( Table 4). The value of perfect fit for NRMSE is 0 and for NSE is 1 [35]. Obtained values for NRMSE and NSE indicate better resemblance between measured and modelled velocities for the summer than for the winter period. In the summer months, there were less intense transitional conditions of the atmosphere, especially the wind, during which the greater discrepancy between measured and modeled velocities is noticed.
Baroclinic flow is recognized in the marina entrance area. NW, NE and E winds caused the surface layer inflow, while winds of S, SW and W direction caused surface layer outflow. Figure 12 shows the hourly averaged current fields in characteristic situations during the occurrence of baroclinic compensation currents. N and SE winds coincide with the rainfall occurrence and the coastal spring activation. In the winter period two similar situations were detected during which N wind, contrary to expected, caused outflow current in the marina entrance area (on 7 March 15:00-21:00 h with average speed of 3.97 m/s, and on 8 March 5:00 h-8:00 h with average speed of 5 m/s). In the mentioned situations coastal springs were active and had more influence on surface layer flow than wind had.      The surface mixed-layer in a stratified water body is defined as that part of a water column which is directly influenced by momentum and turbulence introduced by a wind surface stress [37]. The spatial distribution of turbulent kinetic energy (TKE) (Figure 13) indicates the location of its production and diffuse propagation. The winds of N, SW and NE direction had the most intensive impact on the diffusion of TKE down through the water column ( Figure 13). The mixed-layer deepens, and a homogenous density distribution is established.  The surface mixed-layer in a stratified water body is defined as that part of a water column which is directly influenced by momentum and turbulence introduced by a wind surface stress [37]. The spatial distribution of turbulent kinetic energy (TKE) (Figure 13) indicates the location of its production and diffuse propagation. The winds of N, SW and NE direction had the most intensive impact on the diffusion of TKE down through the water column ( Figure 13). The mixed-layer deepens, and a homogenous density distribution is established. The surface mixed-layer in a stratified water body is defined as that part of a water column which is directly influenced by momentum and turbulence introduced by a wind surface stress [37]. The spatial distribution of turbulent kinetic energy (TKE) (Figure 13) indicates the location of its production and diffuse propagation. The winds of N, SW and NE direction had the most intensive impact on the diffusion of TKE down through the water column ( Figure 13). The mixed-layer deepens, and a homogenous density distribution is established. Oakey and Elliott [38] reported a relationship between the dissipation rate of mechanical energy in the surface mixed-layer and the wind energy. The correlation between the energy dissipation rate and wind forcing showed a linear relationship between εINT and (U10) 3 . Here, εINT is an average value of the dissipation rate integrated over the depth, and U10 is an average wind speed measured at 10 m above the water surface. Figure 14 shows the mentioned linear relationship for the winds of N, NE and SW directions with the R 2 = 0.99 for all three directions. The goal of the model calibration was to obtain the best similarity between the integrated discharge measured by PCM and modelled through one culvert pipe. Figure 15 shows the results of the measured and modelled integrated discharge through one culvert pipe. Oakey and Elliott [38] reported a relationship between the dissipation rate of mechanical energy in the surface mixed-layer and the wind energy. The correlation between the energy dissipation rate and wind forcing showed a linear relationship between ε INT and (U 10 ) 3 . Here, ε INT is an average value of the dissipation rate integrated over the depth, and U 10 is an average wind speed measured at 10 m above the water surface. Figure 14 shows the mentioned linear relationship for the winds of N, NE and SW directions with the R 2 = 0.99 for all three directions. Oakey and Elliott [38] reported a relationship between the dissipation rate of mechanical energy in the surface mixed-layer and the wind energy. The correlation between the energy dissipation rate and wind forcing showed a linear relationship between εINT and (U10) 3 . Here, εINT is an average value of the dissipation rate integrated over the depth, and U10 is an average wind speed measured at 10 m above the water surface. Figure 14 shows the mentioned linear relationship for the winds of N, NE and SW directions with the R 2 = 0.99 for all three directions. The goal of the model calibration was to obtain the best similarity between the integrated discharge measured by PCM and modelled through one culvert pipe. Figure 15 shows the results of the measured and modelled integrated discharge through one culvert pipe. The goal of the model calibration was to obtain the best similarity between the integrated discharge measured by PCM and modelled through one culvert pipe. Figure 15 shows the results of the measured and modelled integrated discharge through one culvert pipe. In the winter measurement campaign, the integrated measured discharge through the one culvert and obtained from the model shows better agreement than in the summer period. The results of the model integrated discharges include the contribution of waves and coastal springs along the marina coastline (as explained in Section 2). The coastal springs flow rate was calculated according to curves in Figure 9. In the winter period, all measured precipitation was taken into account for source inflow calculation. In the summer period, it was assumed that the rain on 7 August (yellow bar in Figure 15) did not activate springs in the marina and there was no noticeable culvert outflow measured by the PCM device. As mentioned before, in the summer period during the low wind speeds of NW direction, significant culvert inflow due to an unknown flow generator is recorded and cannot be obtained by model.

Wave Generation Model Setup
A regional model encompassing the entire Adriatic Sea was established in previous research in order to validate the numerical wave generation model. The Figure 16 shows the domain included within the numerical model and how the domain was meshed using finite volumes. The distance between the numerical nodes is ranging from 10 km in the deep-water area to 250 m in the shallow water.
The spectral numerical model Mike 21/SW [39,40] was used for the numerical analysis with which it is possible to simulate the generation, deformation and attenuation of gravitational wind waves in open-sea and coastal areas. The model computes the evolution of wave action density by solving the action balance equation [41,42]. To discretize the wave spectra in the frequency domain, a logarithmic scale was employed with a minimum frequency of 0.08 Hz (wave period 12.5 s) up to a maximum frequency of 0.95 Hz (wave period 1.05 s), with 28 discrete steps. The model applies the phenomena of wave generation due to wind action, nonlinear wave interaction, refraction as well as the dissipative process caused by surface white-capping. The multi sequential Euler explicit method was used for the discretization of the temporal term and the 'upwind' numerical scheme for the convective flux term.
Initial condition (15 February 2017) is defined with zero wave spectra, meaning the absence of initial wave movement in the modeled area. During the whole analyzed period (15 February 2017 to In the winter measurement campaign, the integrated measured discharge through the one culvert and obtained from the model shows better agreement than in the summer period. The results of the model integrated discharges include the contribution of waves and coastal springs along the marina coastline (as explained in Section 2). The coastal springs flow rate was calculated according to curves in Figure 9. In the winter period, all measured precipitation was taken into account for source inflow calculation. In the summer period, it was assumed that the rain on 7 August (yellow bar in Figure 15) did not activate springs in the marina and there was no noticeable culvert outflow measured by the PCM device. As mentioned before, in the summer period during the low wind speeds of NW direction, significant culvert inflow due to an unknown flow generator is recorded and cannot be obtained by model.

Wave Generation Model Setup
A regional model encompassing the entire Adriatic Sea was established in previous research in order to validate the numerical wave generation model. The Figure 16 shows the domain included within the numerical model and how the domain was meshed using finite volumes. The distance between the numerical nodes is ranging from 10 km in the deep-water area to 250 m in the shallow water. The spectral numerical model Mike 21/SW [39,40] was used for the numerical analysis with which it is possible to simulate the generation, deformation and attenuation of gravitational wind waves in open-sea and coastal areas. The model computes the evolution of wave action density by solving the action balance equation [41,42]. To discretize the wave spectra in the frequency domain, a logarithmic scale was employed with a minimum frequency of 0.08 Hz (wave period 12.5 s) up to a maximum frequency of 0.95 Hz (wave period 1.05 s), with 28 discrete steps. The model applies the phenomena of wave generation due to wind action, nonlinear wave interaction, refraction as well as the dissipative process caused by surface white-capping. The multi sequential Euler explicit method was used for the discretization of the temporal term and the 'upwind' numerical scheme for the convective flux term.
Initial condition (15 February 2017) is defined with zero wave spectra, meaning the absence of initial wave movement in the modeled area. During the whole analyzed period (15 February 2017 to 23 March 2017) the regional model was forced using wind fields from the atmospheric model Aladin-HR with spatial resolution from 4 km and time resolution from 3 h [43,44].
Water 2018, 10, x FOR PEER REVIEW 14 of 22 23 March 2017) the regional model was forced using wind fields from the atmospheric model Aladin-HR with spatial resolution from 4 km and time resolution from 3 h [43,44]. Two sources were used for calibration of the wave generation model. The first was measurements from the wave rider station V1 situated in the northern Adriatic open-sea area (see Figure 16 (ϕ = 44°44.57′ N, λ = 13°17.98′ E). Figure 17 shows the comparison of measured and modelled significant wave heights at the wave recording station during the period from 15 November 2007 to 15 November 2008. Statistical errors of modelled values for half-hourly average wave heights at the wave recording station in relation to measured values for the simulation period have the following feature(s): average error AE = 0.064; root mean square error RMSE = 0.0038. The second source on wind-wave data represents measurements at wave rider site V2 in front of the main breakwater in Rijeka Harbor (see Figure 16, station V2). Figure 18 shows the relationship between measured significant wave heights HS and the peak period TP (wave rider site V2) against measured wind speed Vwind, for wind activity from SSE direction.  Two sources were used for calibration of the wave generation model. The first was measurements from the wave rider station V1 situated in the northern Adriatic open-sea area (see Figure 16 (φ = 44 • 44.57 N, λ = 13 • 17.98 E). Figure 17 shows the comparison of measured and modelled significant wave heights at the wave recording station during the period from 15 November 2007 to 15 November 2008. Statistical errors of modelled values for half-hourly average wave heights at the wave recording station in relation to measured values for the simulation period have the following feature(s): average error AE = 0.064; root mean square error RMSE = 0.0038. The second source on wind-wave data represents measurements at wave rider site V2 in front of the main breakwater in Rijeka Harbor (see Figure 16, station V2). Figure 18 shows the relationship between measured significant wave heights H S and the peak period T P (wave rider site V2) against measured wind speed V wind , for wind activity from SSE direction.
After the calibration process of the wave generation model was done, it was used for the wave generation during the winter and summer period along which the culvert discharge was measured in the marina Opatija. Figure 19 shows the values of the significant wave heights at station 1 (Figure 3).
The differences of measured (station 1) and modelled wave height have been analyzed with a normalized root mean square error (NRMSE) approach [42], gaining the NRMSE = 0.63 for winter measurement period and NRMSE = 0.51 for summer measurement period.
2007 to 15 November 2008. Statistical errors of modelled values for half-hourly average wave heights at the wave recording station in relation to measured values for the simulation period have the following feature(s): average error AE = 0.064; root mean square error RMSE = 0.0038. The second source on wind-wave data represents measurements at wave rider site V2 in front of the main breakwater in Rijeka Harbor (see Figure 16, station V2). Figure 18 shows the relationship between measured significant wave heights HS and the peak period TP (wave rider site V2) against measured wind speed Vwind, for wind activity from SSE direction.  After the calibration process of the wave generation model was done, it was used for the wave generation during the winter and summer period along which the culvert discharge was measured in the marina Opatija. Figure 19 shows the values of the significant wave heights at station 1 (Figure 3).  The differences of measured (station 1) and modelled wave height have been analyzed with a normalized root mean square error (NRMSE) approach [42], gaining the NRMSE = 0.63 for winter measurement period and NRMSE = 0.51 for summer measurement period. Figure 19 and the data presented in Tables 2 and 3   After the calibration process of the wave generation model was done, it was used for the wave generation during the winter and summer period along which the culvert discharge was measured in the marina Opatija. Figure 19 shows the values of the significant wave heights at station 1 (Figure 3).  The differences of measured (station 1) and modelled wave height have been analyzed with a normalized root mean square error (NRMSE) approach [42], gaining the NRMSE = 0.63 for winter measurement period and NRMSE = 0.51 for summer measurement period. Figure 19 and the data presented in Tables 2 and 3 Table 5 shows an average number of hours per year (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001) in which wind waves are characterized by: HS > 0.3 m and TP > 2.8 s and incident wave direction in the band from 90° (E) to 180° (S). Also, the average number of hours with recognised adopted wave situations during the winter and summer period are shown, if the months from May to October are assumed as the summer period and the rest as the winter period. An average amount of hours with significant wave situations in the course of one year is 801 h, which is 9% of the year (approximately one month-32 days). According to the Tables 3 and 4, the total volume of the sea exchange (through one culvert pipe) caused by wave action is 13% of the total seawater exchange in 10 years (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001). The mentioned percentages show that even small waves have a significant impact on the sea exchange through the flushing culverts.

Field Measurements and Results from the Referenced Physical Model
In previous research [25,26], the optimum vertical placement of a culvert in regard to the mean sea level was studied in order to achieve the maximum flushing using the influence of waves. A series of tests on a physical model of breakwater in a wave flume were conducted with the goal of investigating the effect of irregular wave fields on velocities inside them.  Table 5 shows an average number of hours per year (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001) in which wind waves are characterized by: H S > 0.3 m and T P > 2.8 s and incident wave direction in the band from 90 • (E) to 180 • (S). Also, the average number of hours with recognised adopted wave situations during the winter and summer period are shown, if the months from May to October are assumed as the summer period and the rest as the winter period. An average amount of hours with significant wave situations in the course of one year is 801 h, which is 9% of the year (approximately one month-32 days). According to the Tables 3 and 4, the total volume of the sea exchange (through one culvert pipe) caused by wave action is 13% of the total seawater exchange in 10 years (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001). The mentioned percentages show that even small waves have a significant impact on the sea exchange through the flushing culverts.

Field Measurements and Results from the Referenced Physical Model
In previous research [25,26], the optimum vertical placement of a culvert in regard to the mean sea level was studied in order to achieve the maximum flushing using the influence of waves. A series of tests on a physical model of breakwater in a wave flume were conducted with the goal of investigating the effect of irregular wave fields on velocities inside them.
Concerning the fully submerged culvert, the velocity generally fluctuated around 0 m/s throughout the testing period. The absolute values of the mean velocities were mainly below 0.1 m/s in more than 90% of the test cases, which suggests a wave damping effect of a fully submerged culvert regarding velocity values. The stated hindering effect was only detected if the water level was at the top of the culvert, irrespective of the incident wave field in approaching the breakwater. On the other hand, if the culvert was partially submerged in water, the mean velocity had a substantially higher value that extends from 0 m/s to 1 m/s. The velocity of the water oscillates predictably in accordance with the sinusoidal function. In conclusion, the velocity magnitudes recorded at the end of the culvert are substantially lower in the case of a fully submerged flushing culvert than in the case of a partially submerged.
The influence of several wave parameters and geometrical culvert features to the velocity at the culvert exit was studied. The authors provided the following empirical expression to evaluate the culvert velocity [25]: where: v p is the velocity at the exit out of the culvert, g is the gravitational constant, d is the culvert diameter, H and L correspond to significant wave height and peak wave length respectively, L c is the culvert length, and w is the submergence coefficient which is defined as the ratio of water level measured from the bottom of the culvert and culvert diameter, consequently it is a dimensionless value. The proposed Equation (3) is now further utilized for the prediction of the average flow through the culvert assuming the following relations: Q M = v p ·A, where a uniform velocity distribution is assumed across the whole wetted cross-section. In accordance to the Froude similarity criterion and length ratio of α = 10, discharge could be defined in the prototype scale Q P using the following ratio factor: Q R = Q P /Q M = α 5/2 = 316.2 (Table 6). Calculated values of Q P (Table 6) indicate the important contribution of waves on seawater exchange through the flushing culvert. Discharges Q P are evaluated primarily for the incident waves that are perpendicular to the breakwater.  (Table 6). For easier comparison, situations with a similar wave field were considered to evaluate the influence of the culvert submergence. The corresponding significant wave height matches 0.6 m, and wave length equals 15 m. The surface elevation (SE) changed between the observed situations.
As observed by the field measurements, the higher the water elevation is above the top of the culvert the smaller the flow rates are. Sometimes flow rate values could also be negative, indicating that water is exiting the marina through the culvert. As the sea level decreases, the flow rate becomes more positive which is observed in both field measurements and laboratory experiments. This could be viewed on the 28 February where the largest inflow is observed (0.08 m 3 /s), and also the largest inflow is evaluated using Equation (3) (0.04 m 3 /s). On the 4 March, the outflow was observed (−0.02 m 3 /s) while the culvert was submerged (SE = 0.17 m) and the outflow was evaluated (−0.07 m 3 /s) using Equation (3). The difference between the observed and evaluated flow (using Equation (3)) is probably due to incident wave angle, which is always perpendicular to the breakwater during the laboratory tests, which is not the case during the field measurements. An additional influence that could cause a flow value discrepancy is the fact that the measured flow rates in the marina Opatija were affected by wind, tidal oscillations, density gradients and other atmospheric conditions. Due to discrepancies, the proposed equation (obtained by laboratory measurements) may serve as the first approximation in the engineering assessment of the sea exchange through the culvert pipes.
The results displayed in Figure 21 also confirm the influence of the surface elevation on the flow rate through the flushing culvert. With the decline of the surface elevation below the top of the culvert, the flow rate tends to the inward direction. On the other hand, the rising of the surface elevation above the top of the culvert produces outflow from the marina. were affected by wind, tidal oscillations, density gradients and other atmospheric conditions. Due to discrepancies, the proposed equation (obtained by laboratory measurements) may serve as the first approximation in the engineering assessment of the sea exchange through the culvert pipes. The results displayed in Figure 21 also confirm the influence of the surface elevation on the flow rate through the flushing culvert. With the decline of the surface elevation below the top of the culvert, the flow rate tends to the inward direction. On the other hand, the rising of the surface elevation above the top of the culvert produces outflow from the marina.

Flushing Culverts Position
The change in the flushing efficiency by varying the culvert position, when the waves were not taken into consideration, is analyzed in the previous research [19]. The culvert's position, opposite to the main harbor entrance, was recognized as optimal for the marina flushing. In the marina Opatija, where culverts are positioned following the main direction of the incident wave propagation, coastal springs and wind waves are recognized as flow generators affecting the marina flushing. Statistical values regarding the flow measured inside the flushing culvert are shown in Table 7. The entire sea exchange is calculated as the sum of absolute volume passing through the culvert and the marina entrance cross section (Vabs = |V1-culvert| + |Ventrance|). The absolute volume passing through the marina entrance cross section is obtained by the model (winter: |Ventrance| = 1,567,178 m 3 , summer: |Ventrance| = 2,576,703 m 3 ). The absolute amount of the sea exchange through the one flushing culvert is 3% of total sea exchange for the winter period and 2% for the summer period ( Table 7). The sea exchange caused by the wave action (Tables 1 and 3) over the total absolute volume obtained through the one culvert pipe is 15% in the winter and 4% in the summer period. During the winter period, the absolute sum of water passing through the one culvert pipe is 5% lower if not taking into account the impact of inflow caused by freshwater sources. Mentioned percentages indicate that flows caused by wind waves have more significant impact on the seawater exchange than freshwater sources. The absolute volume of sea exchange through the eight culverts obtained by model (|V8 culvert| = 341,603 m 3 ) is seven times higher compared to the volume through one culvert pipe. The absolute volume of sea exchange through the eight pipes is 18% of the total sea exchange, which has an important role in the marina water dynamics.

Flushing Culverts Position
The change in the flushing efficiency by varying the culvert position, when the waves were not taken into consideration, is analyzed in the previous research [19]. The culvert's position, opposite to the main harbor entrance, was recognized as optimal for the marina flushing. In the marina Opatija, where culverts are positioned following the main direction of the incident wave propagation, coastal springs and wind waves are recognized as flow generators affecting the marina flushing. Statistical values regarding the flow measured inside the flushing culvert are shown in Table 7. The entire sea exchange is calculated as the sum of absolute volume passing through the culvert and the marina entrance cross section (V abs = |V 1-culvert | + |V entrance |). The absolute volume passing through the marina entrance cross section is obtained by the model (winter: |V entrance | = 1,567,178 m 3 , summer: |V entrance | = 2,576,703 m 3 ). The absolute amount of the sea exchange through the one flushing culvert is 3% of total sea exchange for the winter period and 2% for the summer period ( Table 7). The sea exchange caused by the wave action (Tables 1 and 3) over the total absolute volume obtained through the one culvert pipe is 15% in the winter and 4% in the summer period. During the winter period, the absolute sum of water passing through the one culvert pipe is 5% lower if not taking into account the impact of inflow caused by freshwater sources. Mentioned percentages indicate that flows caused by wind waves have more significant impact on the seawater exchange than freshwater sources. The absolute volume of sea exchange through the eight culverts obtained by model (|V 8 culvert | = 341,603 m 3 ) is seven times higher compared to the volume through one culvert pipe. The absolute volume of sea exchange through the eight pipes is 18% of the total sea exchange, which has an important role in the marina water dynamics.

Estimation of the Marina Flushing Characteristics
Adequate flushing time of the marina basin is one of most important measures for the water quality in marinas. The tidal prism analytical model [38,45] is used for predicting the marina flushing time (T f ). Schwartz and Imberger [11] reported the comparison of the tidal prism model with the measured dilution rate. The tidal prism model describes the general water exchange characteristics of a semi-enclosed basin using a simplified analytical model (Equation (4)). The present formulation considers only the tide-induced and freshwater flushing effects: where: P indicates the tidal prism (in this case: 19,292 m 3 ), V is the mean basin volume (192,925 m 3 ), T is the tidal period (12 h), Q R designates the freshwater discharge (0.5 m 3 /h), and b is the return flow factor (0.7 according to [16]). Taking into account the mentioned values appropriate for ACI marina Opatija, the flushing time of the marina according to Equation (4) equals 439 h. Furthermore, this value was compared with the flushing time calculated by the use of the established numerical model. The flushing time was evaluated utilizing the CSTR model which considers the time needed for the concentration to fall to e −1 (37%) as the flushing time of the semi-enclosed basin [45]. A flushing time of 65 h is calculated by the numerical model utilizing the transport module for the tracer scalar field in conjunction with the calibrated hydrodynamic model, which is significantly lower than the time calculated by the analytical model. For the studied marina, the calculated T f using the analytical model was seven times higher than that evaluated using the numerical model. Due to this significant discrepancy between the results, it could be concluded that the analytical model does not offer a satisfactory value of marina flushing time. Some other researches also reported discrepancies where the mean flushing times were being overestimated by the analytical model [11,45].

Conclusions
Water quality of semi-enclosed basins depends on the flushing efficiency. Flushing culverts are the most cost-effective and relatively simple solution to improve the water renewal within an enclosed maritime area. Field measurements were performed inside and in the vicinity of ACI Marina Opatija where eight parallel flushing culverts are embedded in the primary breakwater. The 3D sea circulation and wave generation model were calibrated for the winter and summer measurement period. Flow generators contribution to the flow through the flushing culverts were analyzed.
The incident wind waves approaching the marina in a band from 90 • (E) to 180 • (S) provide the most intense culvert inflow. Even small waves (H S ≥ 0.3 m; T P ≥ 2.8 s) have a significant impact on the flow through the culverts. The sea exchange caused by the wave action over the total absolute volume obtained through one culvert pipe is 15% in the winter and 4% in the summer period. When the observed period is expanded to 10 years, it is concluded that wind waves could produce a considerable volume of the water exchange through the culverts. In relative relation to other flow generators, waves caused the 13% of the total seawater exchange. Coastal springs flow is recognised along the marina coastline after rainfall events and is more dominant in the winter than in the summer period. Their activation effects the culvert outflow and surface layer current. Its absolute volume accounts 5% of the total absolute volume obtained through the one culvert pipe in the winter period. The results from the physical model indicate that the flow through the culvert pipe depends on the culvert submergence. Field measurement confirmed the influence of the surface elevation on the flow rate through the culvert pipe. The higher the water elevation is above the top of the culvert, the smaller the flow rates are.
In conclusion, to achieve more efficient seawater exchange in the marina basin, it is suggested to position flushing culvert axis at the mean sea level following the main direction of the incident wave propagation.

Conflicts of Interest:
The authors declare no conflict of interest.