A very important aspect in the planning, design, and maintenance of a harbor is to determine the response of the harbor basin to incident waves. The Saras harbor in South Sardinia (Italy) has been experiencing significant wave disturbance that affects the safety of mooring and operational activities. In the framework of a large research, this paper summarises a comprehensive description of new energy absorbing structures that can be seen as an efficient alternative to more traditional structures when limited by economic or technical constraints. Specifically, the paper presents the results of a graphical preliminary approach and a numerical modelling that solves the enhanced Boussinesq equations in two horizontal dimensions.

1. Introduction

Harbor tranquillity is an essential condition for safe anchorage and efficient port operation. Measuring the distribution of wave height in a given harbour layout should be related to the ship motion [1]. Here, we focus on the motion of service boats (tug, rescue boat, and coast guard) moored in water of finite depth and the viewpoint of maritime operation is not considered. Service boats remain within the harbor during all weather conditions and navigate at the harbour entrance particularly during severe storm conditions. Planning and design of harbor facilities should be based on an comprehensive understanding of the offshore wave modelled from an extreme wave analysis [2] and the estimation of the waves in the harbor for different layout configurations of breakwaters and other protective facilities. Facilities include energy dissipaters in the form of dentated and serrated blocks distributed on the vertical face of quay structures inside the harbor [3]. The amount of wave energy dissipation depends greatly on the slope provided. However, the milder the slope is the more expensive the structure is. In addition, mild slope may not be a technically feasible solution particularly in the case of upgrading of existing port-related waterfront structures [4]. Serrated and dentated vertical quay structures are expected to hydrodynamically replace the mild slope structures from dissipating the incident wave energy and hence expected to reduce the wave reflection [5].

This paper presents a comprehensive application of graphical and numerical methods to the estimation of harbor tranquillity at the Saras harbor located on the south-west coast of Sardinia, Italy, where the Saras Group owns and operates integrated refinery and gasification plants. A refinery jetty extending 1.6 km into the Mediterranean was constructed to accept tankers of all sizes up to 300′000 DWT. The 10 berths of the jetty, which are about 10 miles south of the port of Cagliari, can take tankers with drafts of 20–39.6 ft. All crude for the refinery and petroleum products is transported by tanker. Service boats (tug, rescue boat, and coast guard) in the Saras harbor provide direct tanker assist service at the jetty. The main aim of the paper is to assess the improvement of harbor tranquillity at the Saras harbor when vertical dentated wall is designed in a large portion of the waterfront where a vertical plane wall exists. The manuscript is organized as follows. Section 2 describes the vertical dentated wall and its design in the Saras harbor. An extreme wave analysis from a 20-year wave dataset is presented in Section 3. The problem of excessive wave disturbance at the Saras harbor was explored in both graphical and numerical models. In Section 4, a graphical method is applied for the estimation of wave height distribution in the harbor. Section 5 shows the application of a Boussinesq numerical model and a comparison of results with the graphical method is presented. The final section gives some suggestions for future developments.

2. Wave Absorbing Structures

The hydrodynamic performance of rough structures has been investigated in previous physical and numerical model studies (e.g., [69]), considering a wide range of waves, both regular and irregular, and level conditions. It is found that the investigation of [3, 5] on wave reflection from a dentated and serrated sloped seawall is a unicum in the literature. The physical modelling results have confirmed that the introduction of a number of rectangular blocks in the form of dents or serrations on a plane vertical seawall hydrodynamically replaces the mild sloped seawall from dissipating the incident wave energy and hence reduces the wave reflection. Specifically, they investigated the effect of wave height (in terms of wave steepness, Hs/Lp, being Hs the significant wave and Lp the wavelength corresponding to the peak period), wave period (in terms of relative water depth, d/Lp), and surf similarity parameter (ξ) on the reflection coefficient (Kr). The coefficient of wave reflection Kr is defined as the ration of reflected wave height Hr to incident wave height Hi. Here, both the incident and reflected wave heights refer to the significant heights. Due to the higher energy dissipation, Kr values decrease significantly with decrease of the seawall slope for steeper waves either in the case of plane, dentated, or serrated walls. Here, we focus on the results on random waves acting on vertical seawalls. The effect of seawall type on the reflection coefficient Kr for different peak wave periods is provided in Figure 1. For d/Lp = 0.09 (long wave), the reflection characteristics of all three types of seawall are almost the same. For d/Lp > 0.09, compared to plane and serrated seawall, dentated seawall assures a significant reduction for a large range of d/Lp > 0.09, while the performance of dentated and serrated is almost similar for d/Lp = 0.45. In addition, Neelamani and Sandhya [5] proposed predictive equations for the reflection coefficient Kr for regular and irregular waves, vertical and sloped walls, and plane, dentated, and serrated walls. Here, the predictive equation for vertical dentated wall with random waves is shown below:

Goda [1] reports approximate values of reflection coefficients from various sources (Table 1). The range in the coefficient for a vertical wall depends on the degree of wave overtopping.

Along the problem of energy dissipation, harbour tranquillity is related to the estimation of the propagation of the reflected waves. Based on the theory of geometrical optics, each component of a random sea is assumed to be reflected at an angle equal to the angle of incidence and to continue to propagate in that direction. An exception is due to the phenomenon called Mach-stem reflection [10, 11]. Having reflected waves a finite length along their crest lines, reflected waves disperse during propagation away from the source of reflection in a manner similar to the phenomenon of wave diffraction. Goda [1] proposed an engineering approximation of the spatial distribution of reflected wave height by making use of diffraction diagrams of random sea waves for an opening between breakwaters. This allows properly taking into account directional spectral characteristics in 2D propagation of sea waves. Then the effect of reflected waves on the harbor agitation can be estimated in terms of wave heights by the principle of summation of energy component as in the following:where HT denotes the significant height of the superposed waves and (Hr)1, (Hr)2, … represent the significant heights of the waves reflected at various structures. Specifically, Equation (2) is applicable at a distance from the reflective structure more than one wavelength.

2.1. Saras Harbor Case Study

The Saras harbor is located in the Gulf of Cagliari on the south-west coast of Sardinia at a latitude 39° 4′57″N and longitude 9° 1′21″E (Figure 2). The harbor protection is achieved by the construction of two main shore-connected rubble-mound breakwaters (moles A1 and A2 in Figure 2) with one layer of armor quarry rocks. Waterfront structures have a complex layout and are composed of seawall (S) and quaywall (Q) (with a berthing face for boats) both of gravity, vertical and impermeable structure types. Unfortunately, the basin confined by the structures does not provide safe and efficient port operation due to wave disturbance. As suitable berthing conditions are not provided, minor structural modifications of waterfront structures should be considered. At this time, major modifications as a new harbor layout cannot be proposed for economic and authorization constraints. The pattern of wave propagation into the harbour basin is mainly governed by diffraction and reflection.

The proposed design for the improvement of the Saras harbor tranquillity was to provide with wave-absorbing revetments. Specifically in the portion of the waterfront from where the outer sea can be viewed, the vertical plane seawall will be covered by a dentated structure of marine-grade stainless steel (Figure 3). Being the entrance too wide, as can be seen in Figure 2, the outer sea can be viewed from every basin location through the harbor entrance. Then, the design incoming waves are assumed those that would have arrived in the same locations in the absence of the harbor. Specifically, incoming waves with mean wave direction at the harbour entrance between 0°N and 120°N are simulated in the following.

3. Wave Analysis: Extreme and Ordinary Climate

Wave data from this analysis come from the buoy data provided by the Italian buoy network RON (Rete Ondametrica Nazionale, Italy). Along the Italian coasts, RON includes 15 directional wave buoys collecting 20-minute series of wave and meteorological parameters every 30 minutes. Significant wave height, peak period, and mean direction were used in this study from December 15, 2009 to December 31, 2014 collected from the Watchkeeper buoy with particle-following sensor technology deployed at a latitude 39°06′54″N and longitude 9°24′18″E. In Figure 4, a rose chart of significant wave heights illustrates the percentage of wave measurements coming from different directions. The rose shows that two main directions are consistently predominant in magnitude and frequency: 170–180°N (%) and 260–270°N (%), corresponding to the Scirocco (southeast wind) and Libeccio (southwest wind) winds in the Mediterranean Sea, respectively.

The extreme wave analysis (EWA) was performed on sea store values. Following Boccotti [12], a sea storm as “a sequence of sea states in which the significant wave height Hs exceeds the threshold Hcrit and does not fall below this threshold for a continuous time interval greater than 12 hours” is defined in this study. Regarding Hcrit, it depends upon the given location. For enclosed seas in the Mediterranean Sea, Sulis et al. [2] considered an Hcrit equal to 2.0 m. Each storm was represented by its peak wave height during the storm. The wave dataset was composed of 53 wave heights (Figure 5). For maritime design, Goda’s peak over threshold method [1] is one of the main stable models applied to extreme wave analysis [13, 14] and the techniques for EWA were applied for the use of POT values. POT takes in account for the adoption of a second censoring threshold (censoring parameter ν) in order to achieve the highest stochastic independence and homogeneity possible between the samples. As ν value selection is critical [15], a sensitive analysis was performed, and ν = 1 was selected giving a large data.

Previous research [2] performed a detailed analysis of extreme wave height estimates by making use of different cumulative distribution functions (Gumbel, Fisher–Tippett type 2, Weibull, and Generalized Pareto) and fitting methods (least squares method, the method of L-moments, the maximum likelihood method, and the method of moments). The correlation coefficient was used as degree of goodness to fit in order to assign the “optimal” distribution parameter values that best fit the dataset for each distribution function and best-fitting method; then the probability plot correlation coefficient values were used to compare goodness-of-fit test among distribution functions that passed the Chi-square test. Sulis et al. [2] demonstrated that the least squares method is the best for the Gulf of Cagliari in terms of unbiasedness and efficiency, and the GPD and Weibull are generally better than Gumbel and FT-II. In the following, the application of Weibull distribution fitted with the least squares method is presented:with the shape parameter k = 1, the scale parameter A = 0.77 and the location parameter B = 1.33. The fitted Weibull distribution in Equation (3) was used to estimate wave height quantiles for return periods of 2, 5, 10, 20, 50, and 100 years (Figure 6). Specifically, the wave heights Hst estimated for 2- and 5-year return periods are equal to 3.85 m and 4.56 m. To consider the effects of the selected structure over different wave periods, two wave scenarios were constructed with mean periods T0 = 4 and 9 seconds.

4. Application of a Graphical Method

The selected wave conditions are summarised in Table 2. Specifically, three incoming waves were considered with mean wave direction MWD = 0°N, 90°N, and 120°N at the harbor entrance in order to be representative of the all possible directions (0–120°N). Considering a water depth h = 6.5 m at mean higher high water level, wavelengths at the harbor entrance are Lp ≈ 24 and 68 m and the apparent opening width of the entrance from the wave direction of MWD = 0°N is B' ≈ 70 m, thus the relative width B'/Lp ≈ 2.9 and 1.0. In the following, the paper focuses on the MWD = 0°N that produces the worst case of wave disturbance. A preliminary analysis was carried out with a graphical method called the mirror-imaged diffraction method [1, 16]. The method uses diffraction diagrams for random sea waves [1], transfers the geometry of harbor layout in the plane by a mirror image along the boundary of wave reflection, and examines the reflected waves as waves travelling in the mirror-image plane. Goda [1] presented several examples of the method application and the reader is referred to this textbook for an exhaustive presentation. Points Aj (j = 1, 2, … , 11) are the locations at which the wave height is to be estimated (Figure 7). Being B'/Lp ≈ 2.9, for smax = 75, the random wave diffraction diagram is given by [1] at page 85. It is to be noticed that the application slightly underestimates the disturbance phenomenon as the secondary wave reflections are omitted. The resultant wave height at each of the eleven points is estimated using the principle of summation of wave energy of Equation (2) where the reflected wave heights (Hr)1, (Hr)2, … are those calculated in the plane of the two mirror images (Zone B and Zone C) with the vertical seawall and quaywall as reflective structures. When the absorbing structures are considered, the diffraction coefficients in Zone B and Zone C are reduced by a rate equal to the reflection coefficient. Table 3 shows three different cases in the estimation of superimposed waves: Case A is the current configuration with vertical reflective structures and Case B and Case C are related to the design of dentated vertical seawalls in place of the existing seawall and quaywall for incident waves with T0 = 4 and 9 seconds, respectively. Using Figure 1, it can be seen that the reflection coefficient Kr for T0 = 4 and 9 seconds is equal to 0.7 and 0.8. Comparing Case A with Case B and Case C gives useful preliminary information to harbor managers and designers on the improvement in wave disturbance that the construction of these energy dissipating vertical structures can produce. Considering the wider harbor area with eleven points, it can be noticed that the improvement can assume values in a wide range, from 8% to 30%. Specifically, considering the mooring points (A7, A8, A9, A10, and A11), the improvement is higher in the shadow zone (A10 and A11), while points A7, A8, and A9 are mainly interested by incident waves coming from the wide harbor entrance, and the reflected components are less significant. In the latter, a lower reflection coefficient in the designed absorbing structures will not significantly improve the wave disturbance. As expected, dentated structures have a higher performance in reducing the wave disturbance in the case of incoming random sea waves with higher d/Lp values. In particular, Case B and Case C are related to d/Lp = 0.27 and 0.10. In the mooring points (A7, A8, A9, A10, and A11), the energy dissipation for shorter waves is always higher than 15%. From this preliminary graphical approach, the results confirmed that designed energy absorbing structures are efficient in shadow area of the harbor hit by short storm waves of T0 = 4. For longer waves (T0 > 9 seconds), the improvement cannot be very significant.

5. Application to a Numerical Method

The deepwater waves from the EWA with the statistics of Table 2 were transformed to shallow water using nearshore spectral wind-wave (NSW) model of DHI [17]. NWS is a stationary, directionally decoupled, parametric model which describes the propagation and decay of short period and short crested waves in nearshore areas. Sulis and Annis [18] provide details of the nearshore wave transformation. A detailed bathymetric model at the entrance and in the inner harbor was generated using a single-beam echo-sounder (Figure 8). For examining the inner harbor wave propagation in a good trade-off between computational times, and accuracy and stability, a bathymetry with a grid resolution of 4 m × 4 m was chosen. A model area covered approximately 120 m × 120 m and was rotated of 35° in such a way that the offshore boundary is parallel to the deepwater contours (Figure 9). This wave and bathymetric configuration was used to assess wave disturbance into the Saras harbor using Boussinesq wave (BW) model of DHI [19]. The BW model solves the enhanced Boussinesq expressed in two horizontal dimensions [20] in terms of the free surface elevation, ξ and the depth-integrated velocity-components, along x-axis (P) and y-axis (Q). Boussinesq equations are solved by implicit finite difference techniques and the convective terms are discretised using central differences and the normal ADI algorithm for numerical integration. Partial reflection from revetment, vertical seawall, and quaywall were included in the simulations by specifying porosity values corresponding to reflection characteristics of the specified structures. Wave height and ratio of superimposed waves were simulated for different combinations of incident wave heights, periods, and directions as reported in Table 2. In the following, the paper summaries on the worst condition of a MWD = 0°N that produces the most significant wave penetration into the harbor. Considering the Saras harbor orientation with respect to the Cagliari Gulf, a wave of T0 = 4 and Hst = 3.85 m was chosen for MWD = 0°N. Two structural configurations were simulated: the current configuration with vertical impermeable seawall (Conf 1) and quaywall and the configuration with the dentated structures (Conf 2). The maps of wave disturbance coefficient for both simulations are reported in Figure 10 for an immediate visual comparison. It is evident, also by a qualitative point of view, that the introduction of energy dissipating structures in forms of dentated vertical structures significantly reduces the disturbance of superimposed waves specifically in the mooring area where Conf 1 does not allow service boats for an essential condition for safe anchorage and efficient port operation. In the mooring area, the incoming waves are significantly amplified through reflection at the quaywall reaching coefficient values higher than 0.7 (orange and red values in the legend). Conf 2 will assure lower wave disturbance with the coefficient in the mooring area always lower than 0.6. Considering the A2 point where the maximum superimposition was simulated by the BW model, the maximum wave height Hmax in the current configuration Conf 1 was equal to 95.5 cm, while in the designed configuration Conf 2 results show Hmax = 65.1 cm with an improvement in wave disturbance higher than 31%.

6. Conclusions

This paper summarizes a large research in the field of wave disturbance applied to the Saras harbor located in South Sardinia (Italy) where the current harbor configuration does not assure safe operational conditions. Due to technical restrictions, traditional wave absorber cells cannot be realized and new, experimental approach should be considered. Among these approaches, the dentated and serrated vertical seawalls appear to give the possibility of reducing the high wave disturbance particularly for short wave. A graphical method and a Boussinesq numerical method were applied to simulate the degree of the improvement when the designed new structures are constructed. While the graphical approach gives an approximated solution, it can be seen as a preliminary step in the structure design when no detailed information on bathymetry or additional forces (e.g., wind and tides) are available. When these additional forces become available, more accurate numerical methods (e.g., [21]) can be used also provided by well-known commercial software. In this paper, the application of a Boussinesq wave model is presented. The results suggest that the dentated structural solution can assure safer condition for service boats in the Saras harbor. In addition, the comparison between the graphical and numerical models shows that the latter simulates a significantly higher improvement in the all mooring area.

Data Availability

Data are public (e.g., wave data).

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.