Journal cover Journal topic
Natural Hazards and Earth System Sciences An interactive open-access journal of the European Geosciences Union
Journal topic
Nat. Hazards Earth Syst. Sci., 18, 2675-2695, 2018
https://doi.org/10.5194/nhess-18-2675-2018
Nat. Hazards Earth Syst. Sci., 18, 2675-2695, 2018
https://doi.org/10.5194/nhess-18-2675-2018

Research article 22 Oct 2018

Research article | 22 Oct 2018

# Implementation and validation of a new operational wave forecasting system of the Mediterranean Monitoring and Forecasting Centre in the framework of the Copernicus Marine Environment Monitoring Service

Implementation and validation of a new operational wave forecasting system
Michalis Ravdas1,*, Anna Zacharioudaki1,*, and Gerasimos Korres1 Michalis Ravdas et al.
• 1Hellenic Centre for Marine Research, P.O. Box 712, 19013 Anavyssos, Hellas, Greece
• *These authors contributed equally to this work.
Abstract

Within the framework of the Copernicus Marine Environment Monitoring Service (CMEMS), an operational wave forecasting system for the Mediterranean Sea has been implemented by the Hellenic Centre for Marine Research (HCMR) and evaluated through a series of preoperational tests and subsequently for 1 full year of simulations (2014). The system is based on the WAM model and it has been developed as a nested sequence of two computational grids to ensure that occasional remote swell propagating from the North Atlantic correctly enters the Mediterranean Sea through the Strait of Gibraltar. The Mediterranean model has a grid spacing of 1∕24. It is driven with 6-hourly analysis and 5-day forecast 10 m ECMWF winds. It accounts for shoaling and refraction due to bathymetry and surface currents, which are provided in offline mode by CMEMS. Extensive statistics on the system performance have been calculated by comparing model results with in situ and satellite observations. Overall, the significant wave height is accurately simulated by the model while less accurate but reasonably good results are obtained for the mean wave period. In both cases, the model performs optimally at offshore wave buoy locations and well-exposed Mediterranean subregions. Within enclosed basins and near the coast, unresolved topography by the wind and wave models and fetch limitations cause the wave model performance to deteriorate. Model performance is better in winter when the wave conditions are well defined. On the whole, the new forecast system provides reliable forecasts. Future improvements include data assimilation and higher-resolution wind forcing.

1 Introduction

In recent years the requirements of the marine industry for real-time wave forecasts have increased substantially. Various sectors such as, for example, maritime transport, shipping, and offshore mineral industry require accurate wave forecasts in order to secure operations at sea, save fleet fuel consumption using more accurate routing, and prevent from potential ship and platform oil spill drift. In addition, a detailed wave information along coastal regions is crucial for coast guards and port authorities as there is the need to anticipate wave conditions that could interfere in ships arriving and leaving the harbours. Furthermore, wave forecasting provides significant advantages for the offshore wind and the wave energy so as to schedule installation and maintenance activities, to define control strategies according to the predominant wave conditions, or to plan for storm events. Scientifically, waves are a bridge between the ocean and the atmosphere, playing a key role in air–sea interaction and are also an important mixing agent with an active role in erosion and resuspension processes.

The dramatic increase in computing power and the enhanced understanding of the physical processes responsible for wave generation, evolution, and dissipation have resulted in third-generation wave models which use first principles in the integration of an action or energy balance equation (Tolman, 1992) based on the sophisticated physics pertaining to wave generation, propagation, and decay mechanisms. Wave models like WAM (WAMDI Group, 1988; Komen et al., 1994), WAVEWATCH III (Tolman, 1991, 1999), and SWAN (Booij et al., 1999) are being used by many meteorological and oceanographic operational centres and have been reasonably successful in operational wave predictions at the global, regional, and coastal scales. One of the pioneers in the implementation and the development of wave analysis and forecast systems is the European Centre for Medium Weather Forecasts (ECMWF), which has provided daily medium-range global wave forecasts up to 10 days ahead since 1992. With time, more centres (e.g. UK Met Office, the Service Hydrographique et Océanographique de la Marine (SHOM), and several others) began to use state-of-the-art third-generation numerical wave models in operational forecasting. In addition, an international wave forecast inter-comparison project was established (Bidlot, 2007), coordinated by the ECMWF, to evaluate forecast quality and performance and identify areas of potential improvement (Breivik et al., 2015).

The EU-funded Copernicus Marine Environmental Monitoring Service (CMEMS) driven by the requirements of a large user panel in need of wave information in all ocean basins has enriched its portfolio since early 2017 with the provision of open, cost-free, and quality-controlled wave products for the global ocean and regional European seas. CMEMS is based on a strong European partnership with more than 50 marine operational and research centres in Europe involved in marine monitoring and forecasting services and their evolution, providing a wide range of marine measurements of social and environmental value such as ocean currents, temperature, salinity, sea level, pelagic biogeochemistry, and waves. The backbone of the CMEMS relies on a Central Information System (CIS) and an architecture of production centres inherited from the MyOcean projects for both observations (Thematic Assembly Centre – TAC) and modelling–assimilation (Monitoring and Forecasting Centre – MFC). The MFCs are distributed according to the marine area covered and generate model-based products including analysis of the current situation, forecasts of the situation a few days in advance, and the provision of retrospective data records (reanalysis). Detailed information on the systems and products is on CMEMS web site: http://marine.copernicus.eu/ (last access: September 2018). The Mediterranean Sea (Med) MFC is composed of INGV (Italy), HCMR (Greece), OGS (Italy), and the Euro-Mediterranean Centre on Climate Change (CMCC, Italy). It is an expert consortium with profound expertise in the Mediterranean phenomenology and dynamics, from waves to currents, and biogeochemistry, and it provides regular and systematic information about the physical state of the ocean and the dynamics of the marine ecosystem of the basin.

In this study we present the wave component of the Med MFC, a high-resolution operational wave forecasting system (hereafter called Med-waves), which has been developed by HCMR and provides daily accurate products – wave simulations and 5-day forecasts – of the wave environment of the Mediterranean Sea to the general public through the CMEMS portal. The system employs the WAM Cycle 4.5.4 (Günther and Behrens, 2012), a modernized and improved version of the WAM model. In this study, a modelling system consisting of a nesting sequence of two computational grids (North Atlantic and Mediterranean) has been developed with the fine-grid model covering the Mediterranean Sea and the coarse-grid model the North Atlantic. A nesting approach of this kind enables us to properly simulate the effect of the remotely generated Atlantic swell into the Mediterranean Sea as it passes through the Strait of Gibraltar. In fact, Cavaleri and Sclavo (2006) pointed out that the narrow Strait of Gibraltar appreciably affects the wave climate in the close-by area of the Alboran Sea and it is often neglected in wave modelling systems of the Mediterranean Sea. Moreover, the system incorporates offline coupling with general circulation models (CMEMS global and Mediterranean analysis and forecast systems) to provide surface currents for wave refraction to both nests of the modelling system. Refraction due to surface currents impacts the wave spectrum to some extent. Its impact on improving the wave forecast is addressed by Osuna and Wolf (2005), Clementi et al. (2013, 2017), and Mao and Xia (2017). Thus, while in the last years few operational centres have already developed and implemented regional wave forecast systems for the entire Mediterranean, none of these take into account the sensitivity of Mediterranean wave dynamics to the nesting with the Atlantic, at the same time incorporating the surface currents effect on wave refraction. In addition, the system offers an extended set of freely available wave products and has a spatial and spectral resolution high enough to describe with sufficient accuracy the wind wave dynamics over the Mediterranean Basin.

Despite the diversity of wave generation models as well as of atmospheric models, obtaining good quality of short-term wave forecasts for the Mediterranean Sea is still a difficult task due to the large spatial and temporal variability in the surface wind field over the basin (e.g. Bidlot, 2017). Wind-wave models are very sensitive to wind field variations, which result in one of the main source of errors in wave predictions. The sensitivity of wave model prediction to variations in wind forcing fields has been studied by several authors (Komen et al., 1994; Teixeira et al., 1995; Holthuijsen et al., 1996; Ponce and Ocampo-Torres, 1998). This is particularly true for the Mediterranean where the limited contribution of swell to the wave spectrum makes the regional wind conditions the most important factor in determining the local wave state. Bearing in mind this context, the complex structure of the Mediterranean Sea due to the presence of large mountainous islands, protruding peninsulas, jagged coastlines, and sharp orography gradients deeply influences the wind and wave dynamics, especially close to the coast, making the local forecast particularly challenging. Many authors have highlighted the fact that forecasted winds in the Mediterranean are not as accurate as in the open ocean (Cavaleri and Bertotti, 2003, 2004; Bolaños et al., 2005; Signell et al., 2005; Ardhuin et al., 2007; Bolaños-Sanchez et al., 2007) and advocate the necessity of further improvement on the wind field quality as well as the increase in its spatial resolution, especially in enclosed basins such as the Mediterranean Sea (Cavaleri et al., 1991; Cavaleri and Bertotti, 2004, 2009a, b; Bentamy et al., 2007; Lionello et al., 2008). Today, the full coupling of different geophysical models (e.g. atmospheric, ocean, wave) is also a favoured approach to reduce the wind forecasting error, especially in dynamically complicated coastal areas, as it accounts for the complex interactions of waves, currents, and the atmosphere (Wahle et al., 2017). For the Mediterranean Sea, the added value of coupled air–sea models has been demonstrated by many research groups (Artale et al., 2009; Renault et al., 2012; Katsafados et al., 2016; Varlas et al., 2018) particularly in the region of the Adriatic Sea (Pullen et al., 2006; Carniel et al., 2016; Licer et al., 2016; Ricchi et al., 2016) where the synoptic wind regimes are difficult to forecast. However, despite the rapid progress and the encouraging results of this field, there are still open scientific issues such as missing physics and poor parameterizations for the coupled atmospheric–ocean–wave processes in the air–sea interaction zone, especially under extreme wind and ocean conditions (Chen et al., 2007; Soloviev et al., 2014). Furthermore, the coupled models are still computationally expensive and even with increased computational potential to date only few major meteorological and oceanographic centres have the capability to incorporate these high computationally demanding systems at fine spatiotemporal resolution and at large spatial scales into their operational weather–ocean forecasting chains. The recent white paper by Cavaleri et al. (2018) analyses in detail the various aspects and problems affecting the performance of wave models in enclosed and inner seas including the Mediterranean Sea. In this framework, the selection of an appropriate wind forcing is a vital step in the operational wave modelling of the Mediterranean Sea. As such, the Med-waves system in particular uses the analysis and forecast wind fields produced by the ECMWF Integrated Forecasting System (IFS). The quality and appropriateness of ECMWF 10 m winds for wave simulations and forecasting in different areas has been demonstrated by many studies. On a global scale, repetitive statistics have shown that the ECMWF products are, and have been for a long while, the best ones in the world (Bertotti et al., 2011). However, for semi-enclosed basins, the quality of ECMWF wind fields decreases and the wave model underestimates the high wave heights because of the underestimation of high wind speeds (Cavaleri and Bertotti, 2004; Signell et al., 2005; Cavaleri and Scalvo, 2006; Saket et al., 2013) and/or overestimates the lower ones because of the overestimation of low wind speeds by ECMWF (Moeini et al., 2010). In our system, this underestimation of wind is compensated for by reducing the energy loss due to whitecapping, performing a fine tuning of the free parameters of the dissipation function.

We present here the first comprehensive documentation of the system and the evaluation of its accuracy over a period of 1 year (2014). The rest of the paper is organized as follows. A detailed description of the Med-waves modelling system is given in Sect. 2. Section 3 outlines the methodology followed in the model validation, and Sect. 4 is devoted to the validation results including both hindcast and forecast skill evaluation against in situ and satellite observations. Finally, a summary and some concluding remarks are given in Sect. 5.

2 The wave forecasting system

As previously stated, Med-waves is based on the WAM Cycle 4.5.4 wave model, a state-of-the-art third-generation wave model which is a modernized and improved version of the well-known and extensively used WAM Cycle 4 wave model (WAMDI Group, 1988; Komen et al., 1994). Cycle 4.5.4 has been released during the MyWave (“A pan-European concerted and integrated approach to operational wave modelling and forecasting – a complement to GMES MyOcean services”) EU FP7 research project and is freely available to the entire research and forecasting community. WAM solves the wave transport equation explicitly without any presumption on the shape of the wave spectrum. Its source–sink terms include the wind input, whitecapping dissipation, non-linear transfer, and bottom friction. The wind input and whitecapping dissipation source terms of the present cycle of the wave model are a further development based on Janssen's quasilinear theory of wind-wave generation (Janssen, 1989, 1991). The non-linear transfer term is a parameterization of the exact non-linear interactions as proposed by Hasselmann and Hasselmann (1985) and Hasselmann et al. (1985). Lastly, the bottom friction term is based on the empirical JONSWAP model of Hasselmann et al. (1973).

The Med-waves set-up includes a coarse-grid domain with a resolution of 1∕6 covering the North Atlantic Ocean (NA) from 75 W to 10 E and from 10 to 70 N and a nested fine-grid domain with a resolution of 1∕24 covering the Mediterranean Sea from 18.125 W to 36.2917 E and from 30.1875 to 45.9792 N. The areas covered by the two grids are shown in Fig. 1, which is a schematic of the Med-waves system. The bathymetric map has been constructed using the GEBCO 30 s bathymetric data set (GEBCO, 2016) for the Mediterranean Sea model and the ETOPO2 2 min bathymetric data set (NGDC, 2006) for the North Atlantic model.

Figure 1Schematic of the Med-waves system.

The Mediterranean Sea model receive the full wave spectrum at hourly intervals at its Atlantic Ocean open boundary from the North Atlantic model. The latter model is considered to have all of its four boundaries closed with no wave energy propagation from the adjacent seas. Because of the wide geographical coverage of the North Atlantic model, the consideration of closed boundaries does not affect the swell propagation towards the open boundary of the Mediterranean model, which is the main interest of this nesting approach.

The wave spectrum is discretized using 32 frequencies, which cover a logarithmically scaled frequency band from 0.04177 to 0.8018 Hz at intervals of $\frac{df}{f}=\mathrm{0.1}$ and 24 equally spaced directions (15 bin size).

The Mediterranean model runs in shallow-water mode considering wave refraction due to depth and currents in addition to depth-induced wave breaking. The North Atlantic model runs in deep-water mode with wave refraction due to currents only. The North Atlantic model additionally considers wave energy damping due to the presence of sea ice.

Following ECMWF (2015), the tunable whitecapping dissipation coefficients Cdis and δ have been altered from their default values. Specifically, the values of Cdis=1.33 (Cdis=2.1 default) and δ=0.5 (δ=0.6 default) have been adopted. The aim of this tuning is to produce results which are in good agreement with data on fetch-limited growth and with data on the dependence of the surface stress on wave age.

The atmospheric forcing of the Med-waves system is the operational 10 m wind analysis and forecast from the ECMWF global meteorological forecasting system (ECMWF, 2017) that are disseminated at a horizontal resolution of 1∕8. The quality of the ECMWF forecasts is regularly evaluated (Haiden et al., 2016). Sea ice coverage fields are also obtained from ECMWF IFS at the same horizontal resolution as the wind fields. An issue here concerns the temporal resolution of the disseminated ECMWF wind fields as these data are provided on a 3-hourly basis for the first 3 days of the forecast while it is expected that higher-temporal-resolution wind forcing can better capture storm events, especially if the latter have a short duration. For example, Alomar et al. (2014) have shown that on the southern Catalan coast (NW Mediterranean Sea), the best wind speed estimations, compared with the observations, corresponded to a grid size of 4 km and a temporal resolution of 1 h.

Surface current forcing is accounted for in Med-waves. The Mediterranean Sea model is forced by surface currents obtained from the physical forecasting system of the CMEMS Med-MFC at 1∕16 horizontal resolution (CMEMS, 2016a) and the North Atlantic model by surface currents obtained from the physical forecasting system of the CMEMS Global-MFC at 1∕12 horizontal resolution (CMEMS, 2016b). Both physical forecasting systems are based on the Nucleus for European Modelling of the Ocean (NEMO) ocean physics model, which is a state-of-the-art free-surface primitive equation model (Madec, 2008). NEMO is free software used by a large community.

Med-waves is run once per day starting at 12:00:00 UTC. It produces 5-day forecast fields initialized by a 1-day hindcast. The wave hindcast is forced by 6-hourly analysis wind fields and daily averaged analysis current fields. The 5-day forecast is forced by 3-hourly forecast wind fields for the first 3 days and 6-hourly forecast wind fields for the rest of the forecast cycle. Daily averaged forecast currents are used over the entire wave forecast. Sea ice coverage fields are updated at daily frequency and remain constant during the forecast cycle.

Figure 2Wave buoys' location and unique ID code.

Med-waves generates hourly wave fields over the Mediterranean Sea at 1∕24 horizontal resolution. These wave fields correspond either to wave parameters computed by integration of the total wave spectrum or to wave parameters computed using wave spectrum partitioning. In the latter case the complex wave spectrum is partitioned into wind sea and primary and secondary swell. Wind sea is defined as those wave components that are subject to wind forcing while the remaining part of the spectrum is termed swell. Wave components are considered to be subject to wind forcing when

$\begin{array}{}\text{(1)}& c\le \mathrm{1.2}×\mathrm{28}{u}_{*}\mathrm{cos}\left(\mathit{\vartheta }-\mathit{\phi }\right),\end{array}$

where c is the phase speed of the wave component, u* is the wind friction velocity, θ is the direction of wave propagation, and φ is the wind direction. As the swell part of the wave spectrum can be made up of different swell systems with quite distinct characteristics, it is further partitioned into the two most energetic wave systems, the so-called primary and secondary swell. Swell partitioning is carried out following the method proposed by Gerling (1992), which finds the lowest energy threshold value at which upper parts of the spectrum get disconnected, with the process repeated until primary and secondary swell is detected.

Total spectrum and partitioned wave parameters produced by Med-waves and disseminated though CMEMS include spectral significant wave height (Hm0), spectral moments (−1, 0) wave period (Tm−10), spectral moments (0, 2) wave period (Tm02), wave period at the spectral peak or peak period (Tp), mean wave direction (Mdir), wave principal direction at spectral peak, surface Stokes drift U, surface Stokes drift V, spectral significant wind wave height, spectral moments (0, 1) wind wave period, mean wind wave direction, spectral significant primary swell wave height, spectral moments (0, 1) primary swell wave period, mean primary swell wave direction, spectral significant secondary swell wave height, spectral moments (0, 1) secondary swell wave period, and mean secondary swell wave direction.

3 Validation framework

Med-waves has been validated against in situ and satellite observations, focusing on its performance in the Mediterranean Sea. Model output and observations corresponding to the year 2014 have been compared, focusing on the fundamental wave parameters of significant wave height, Hs, and mean wave period, Tm.

In situ measurements of Hs and Tm for 2014 were extracted from the Copernicus In Situ Thematic Assemble Centre (INS-TAC), a component of CMEMS which aims at providing a research and operational framework to develop and deliver in situ observations and derived products based on such observations. Hs measurements from 32 wave buoys within the Mediterranean Sea were available in the examined year. Figure 2 depicts their location and unique ID code. Tm measurements were available from a subset of the depicted buoys which excludes all the buoys offshore from the Italian coastline. To collocate model output and buoy measurements, in space model output was taken at the grid point nearest to the buoy location. In time, buoy measurements within a time window of ±1 h from model output times at 3 h intervals (0, 3, 6 …, etc.) were averaged. Prior to model–buoy collocation, the in situ observations were filtered so as to remove those values accompanied by a bad quality flag (quality flags included in the data files provided by the INS-TAC). After collocation, visual inspection of the data was carried out, which led to some further filtering of spurious data points. In addition, Tm data below 2 s were omitted from the statistical analysis since 0.5 Hz (T=2 s) is a typical cut-off frequency for wave buoys. It is noted that WAM, in contrast to the wave buoys, does not implement a high-frequency cut-off in the computation of Tm but instead includes a high frequency tail extending the calculation to infinity. As a result, model Tm is anticipated to be somewhat biased towards lower values compared to measured Tm.

Figure 3Mediterranean Sea subregions for qualification metrics.

Satellite observations of significant wave height, Hs, and wind speed, U10 (used to obtain HsU10 quality associations), for the year 2014 were obtained from a merged altimeter wave height database set-up at CERSAT – IFREMER (France). This database contains altimeter measurements that have been filtered and corrected (Queffeulou and Croizé-Fillon, 2013). Here, measurements from three satellite missions, Jason-2, Cryosat-2, and SARAL, were used. To collocate model output and satellite observations, the first two were interpolated in time and space to the individual satellite tracks. For each track, corresponding to one satellite pass, along-track pairs of satellite measurements and interpolated model output were averaged over ∼50 km (0.5) grid cells, centred at grid points of the forcing wind model (0.125× 0.125). This averaging is intended to break any spatial correlation present in successive 1 Hz (∼7 km) observations and/or in neighbouring model grid output (Pierre Queffeulou, personal communication, 2015).

Metrics that are commonly applied to assess numerical model skill and are in alignment with the recommendations of the EU FP7 project MyWave (Saulter, 2014) have been used to qualify the Med-waves system within the Mediterranean Sea. These include the RMSE, bias, scatter index (SI), Pearson correlation coefficient (CORR), and best-fit slope (slope). The SI, defined here as the standard deviation of errors (model – observations) relative to the observed mean, as well as being dimensionless, is more appropriate to evaluate the relative closeness of the model output to the observations at different locations compared with the RMSE, which is representative of the size of a “typical” error. The slope corresponds to a best-fit line forced through the origins (zero intercept). In addition to the aforementioned core metrics, merged density scatter and quantile–quantile (QQ) plots are provided. Metrics are computed for the Mediterranean Sea as a whole, for the individual wave buoy locations shown in Fig. 2, and for 17 subregions of which one is in the Atlantic Ocean and 16 are in the Mediterranean Sea (Fig. 3): (atl) Atlantic, (alb) Alboran Sea, (swm1) West South-West Med, (swm2) East South-West Med, (nwm) North West Med, (tyr1) North Tyrrhenian Sea, (tyr2) South Tyrrhenian Sea, (adr1) North Adriatic Sea, (adr2) South Adriatic Sea, (ion1) South-West Ionian Sea, (ion2) South-East Ionian Sea, (ion3) North Ionian 3, (aeg) Aegean Sea, (lev1) West Levantine, (lev2) North-Central Levantine, (lev3) South-Central Levantine, and (lev4) East Levantine. All metrics are evaluated over a period of 1 year (2014). In addition, metrics associated with the full Mediterranean Sea are evaluated seasonally.

4 Validation results

## 4.1 Hindcast significant wave height

### 4.1.1 Comparison with in situ observations

Table 1 shows results of the comparison between hindcast Hs (model data) and in situ observations (reference data), for the Mediterranean Sea as a whole, for the entire year of 2014, and seasonally. In the table, “Entries” refers to the number of model–buoy collocation pairs; i.e. to the sample size available for the computation of the relevant statistics, $\overline{R}$ is the mean reference value, $\stackrel{\mathrm{‾}}{M}$ is the mean model value, and SD R and SD M are the standard deviations of the reference and model data respectively. The remaining quantities are the qualification metrics defined in the previous section. Figure 4 is the respective merged QQ–scatter plot (left), together with the QQ plot alone (right) for visual clarity, for the full 1-year period. In the figure, the QQ–plot is depicted with black crosses. Also shown are the best-fit line forced through the origin (red solid line, left plot) and the 45 reference line (red dashed line).

Table 1Med-waves Hs evaluation against wave buoy Hs, for the full Mediterranean Sea, for a 1-year period (2014) and seasonally.

Figure 4QQ–scatter plots (a) and QQ plot alone (b) of Med-waves output Hs versus wave buoy observations, for the full Mediterranean Sea, for a 1-year period (2014): QQ plot (black crosses), 45 reference line (dashed red line), least-squares best-fit line (red line, a).

Table 1 shows that the typical error (RMSE) varies from 0.17 m in summer to 0.25 m in winter. However, the scatter in summer (0.26) is about 2 % higher than the scatter in winter (0.24) whilst a lower correlation coefficient is associated with the former season. This suggests that the model follows the observations in “stormy” conditions better, with well-defined patterns and higher waves. A similar conclusion has been derived by other studies (Cavaleri and Sclavo, 2006; Ardhuin et al., 2007; Bertotti et al., 2013) with respect to wind and wave modelling performance in the Mediterranean Sea. Like summer, autumn is characterized by lower mean wave height, higher scatter, and a lower correlation coefficient compared to winter and spring. Relative bias (BIAS/$\stackrel{\mathrm{‾}}{R}$) has an approximate variation of 2.5 % (spring) to 5 % (autumn) and is always negative. Accordingly, slopes are also below unity with small variation among seasons. These values are indicative of an underestimation of the wave height in the Mediterranean Sea by the model, a result which is in agreement with the results of a number of operational or preoperational models for the Mediterranean Sea (e.g. Bidlot, 2015; Donatini et al., 2015) and is linked to an underestimation of the wind speed by the ECMWF forcing wind model (see Fig. 7). Overall, the spring statistics are the ones closest to the year-long statistics for the Mediterranean Sea.

Figure 4 depicts the pattern of the agreement between hindcast and observed Hs for different Hs value ranges. The figure reveals that the Hs underestimation by the model mainly occurs for wave heights below 4 m and is rather small. It also shows a QQ plot that is really close to the reference line over most of the Hs range observed, which means that waves of a specific wave height have a very similar probability of occurrence in the hindcast and in the observations. The “outliers” present in the scatter plot, i.e. a number of measured waves of 2–4.5 m in height which are not simulated by the model (not enough evidence was found to remove the depicted outliers from the calculation of the statistics as faulty), correspond to buoy location 61218 in the Adriatic Sea (Fig. 2) and mostly belong to a single storm.

Table 2Med-waves Hs evaluation against wave buoy Hs, for each individual buoy location, for a 1-year period (2014).

Figure 5QQ–scatter plots of Med-waves output Hs versus wave buoy observations at specific wave buoy locations, for a 1-year period (2014): QQ plot (black crosses), 45 reference line (dashed red line), and least-squares best-fit line (red line).

Table 3Med-waves Hs evaluation against satellite Hs, for the full Mediterranean Sea, for a 1-year period (2014).

Figure 6Med-waves Hs evaluation against satellite Hs, for each Mediterranean Sea subregion shown in Fig. 3, for a 1-year period (2014).

Up to now, the overall performance of the Med-waves modelling system at the different wave buoy locations has been analysed independently of the severity of the conditions. Figure 5, similar to Fig. 4, shows the QQ–scatter plots of hindcast Hs versus measured Hs at three buoy locations, exhibiting model performance over the different wave height ranges. The results at these three locations are reasonably representative of the different behaviours of the wave model at the different wave buoy locations in the Mediterranean Sea shown in Fig. 2. Thus, the top left plot shows the behaviour of the model at location 61188, offshore from the border between France and Spain, backed on the west by the Pyrenees Mountains. It is seen that model underestimation occurs throughout the measured Hs range except from the highest percentiles of Hs at which model overestimation is observed. This distribution, with a smaller or larger model underestimate and with a more or less pronounced convergence or overestimate towards the highest waves, is observed at the majority of the wave buoy locations. The bottom left plot corresponds to location 61221, south of the island of Sardinia. There, the model overestimates the observed Hs over the entire Hs range, even more so in the upper end of this range. Considerable model overestimation, mostly over the middle and higher Hs ranges, is observed at all wave buoys associated with unresolved bathymetric features in their surroundings (e.g. 61021, SARON). Over the lower Hs range, convergence or underestimation is also observed in these conditions. At SARON (not shown), the surrounding topography is highly complex, including both orographic and bathymetric effects, resulting in highly scattered data around the reference line. The right plot shows results at buoy location 61197 east of the Balearic Islands. This is a well-exposed offshore wave buoy; consequently, the behaviour of the model at this location is expected to be representative of its performance at well-exposed offshore sites. Relatively small scatter of the data points is shown in the plot with QQ crosses and the best-fit line laying close to the reference line. More specifically, the model converges to the observations for wave heights below about 2 m, somewhat underestimates the observations for wave heights between 2 and 4 m, and tends to overestimate Hs for higher waves. A very similar distribution is found for location 61430 west of the Balearic Islands. Other well-exposed offshore locations present QQ–scatter patterns that are not far from the one shown for location 61197.

### 4.1.2 Comparison with satellite observations

This subsection starts with the comparison of Med-waves hindcast Hs with satellite observations of Hs separately for each satellite. This is carried out for a 1-year period (2014) for the full Mediterranean Sea and for the different subregions defined in Fig. 3. Respective results are shown in Table 3 and Fig. 6.

Table 3 shows that even though the model–satellite comparison behaves similarly for the three different satellites in terms of SI and CORR, a substantially more (>10 %) negative model bias associated with a considerably lower slope is found for Cryosat-2. RMSE is also higher for this satellite. Figure 6 shows that these results are largely consistent among the different Mediterranean subregions although they are more pronounced in the western Mediterranean and the Adriatic Sea. A lower model underestimate of the Cryosat-2 measurements is observed in the Ionian Sea and the eastern Mediterranean. The statistics of model–Jason-2 and model–SARAL comparisons are comparable, with the model exhibiting its best performance when compared to the observations of SARAL.

Table 4Med-waves Hs evaluation against satellite Hs (Jason-2 and SARAL), for the full Mediterranean Sea, for a 1-year period (2014) and seasonally.

Figure 7QQ–scatter plots of (a) ECMWF forcing wind speed U10 versus satellite U10 (Jason-2) and (b) Med-waves Hs versus satellite Hs (Jason-2 and SARAL), for the full Mediterranean Sea, for a 1-year period (2014).

It was decided to exclude the observations of Cryosat-2 from the analysis. Apart from the aforementioned discrepancies, there are other results in the literature to support this decision. Specifically, satellite–buoy comparisons performed by Sepulveda et al. (2015) have shown that SARAL Hs is of better quality than Jason-2 and Cryosat-2 Hs at both open-ocean and coastal buoy sites. In fact, SARAL data are of very high quality with no need of corrections whilst corrections are applied to Jason-2 and Cryosat-2 Hs observations (corrected data are used herein). After corrections, Jason-2 Hs has been found to approximate SARAL Hs well whilst less accurate results have been obtained for Cryosat-2, particularly for wave heights below 1.5 m. For these reasons, in what follows, the comparison of Med-waves Hs is performed against merged satellite observations from the SARAL and Jason-2 satellites, which are of similar accuracy.

Table 4 shows statistics from the comparison of the Med-waves hindcast Hs and satellite observations of Hs, for the full Mediterranean Sea, for a 1-year period and seasonally. Figure 7b shows the corresponding QQ–scatter plot for a 1-year period, for the full Mediterranean Sea. Figure 7a shows an equivalent QQ–scatter plot resulting from the comparison of the ECMWF forcing wind speeds, U10, and Jason-2 measurements of U10 (no U10 available from SARAL or Cryosat-2).

Figure 7a shows that the ECMWF forcing wind model mostly underestimates observed U10, even more so at high wind speeds. An overall model underestimation of 8 % associated with a slope of 0.9 has been computed. Figure 7b also shows an overall Med-waves model underestimation of observed Hs by about 5 % associated with a slope of 0.96. Nevertheless, in this case, the model somewhat underestimates observed Hs over the lower Hs range (<2 m), converges to the observed Hs over the middle Hs range (2–3.5 m), and, generally, somewhat overestimates the larger waves in the data records. This apparent discrepancy between wind and wave scatter distributions is a consequence of the modification of the default values of the whitecapping dissipation coefficients in WAM as described in Sect. 3. A QQ–scatter plot obtained before this modification (not shown) is indeed very similar to the one of the ECMWF wind speeds in Fig. 7. On the whole, Fig. 7 shows that the performance of Med-waves at offshore locations in the Mediterranean Sea (satellite records near the coast are mostly filtered out as unreliable) is very good. Comparing to the equivalent results obtained from the model–buoy comparison (Fig. 4), a very similar pattern of scatter distribution is observed in the two plots, also evident from the orientation of the best-fit lines and the curvature of the QQ plots. A smaller scatter (by about 6 %) with a larger overall bias (by about 2 %) is associated with the model–satellite comparison. SI values compare well at the more exposed wave buoys in the Mediterranean Sea.

Table 5Med-waves Hs evaluation against satellite Hs (Jason-2 and SARAL), for each individual Mediterranean Sea subregion shown in Fig. 3, for a 1-year period (2014).

Table 4 shows the seasonal variation in the Med-waves model performance. RMSE varies from 0.17 m in summer to 0.24 m in winter. SI is highest in summer (0.2) and lowest in winter (0.17). Correlation coefficient varies accordingly. In general, as explained in the previous subsection, a lower scatter with a higher correlation is expected the more well-defined the weather conditions are. Similar to in the model–buoy comparison, bias is negative in all seasons. Its highest relative value (BIAS/$\stackrel{\mathrm{‾}}{R}$) of 7.7 % is computed for autumn and its lowest of 3.5 % for summer. Slope varies from 0.95 in spring and autumn to 0.97 in winter. Overall, Table 4, like Table 1, reveals that the statistics of spring are the most representative of the year-long statistics for the Mediterranean Sea.

Table 5 shows the statistics of the comparison of the Med-waves hindcast Hs and satellite observations of Hs for the different subregions of the Mediterranean Sea defined in Fig. 3. For visualization purposes, Fig. 8b maps the statistics shown in Table 5. In addition, equivalent statistics are mapped for the ECMWF – satellite comparison of wind speeds (Fig. 8a). It is noted that the relative bias (BIAS/$\stackrel{\mathrm{‾}}{R}$) is displayed in the figure. This quantity allows for a more straightforward comparison among the different sub-basins in terms of percentage deviations from the observed mean value. It is also highlighted that the spatial coverage of the model–satellite wind collocations (measurements only from Jason-2) is much more limited than the spatial coverage of the model–satellite wave collocations (measurements from both SARAL and Jason-2). As a consequence, the wave statistics are expected to be more representative of the subregions under consideration compared to the wind statistics. This is particularly true for the Adriatic, the Ligurian, and the Alboran Seas. In addition, the wave statistics have been computed using a sample size of at least 400 data points whilst the wind statistics have been obtained with a minimum sample requirement of 200 data points. Thus, the confidence associated with the wave statistics is higher than the confidence associated with the wind statistics. For the above reasons, the wind metrics presented in Fig. 8 are interpreted with caution.

Figure 8ECMWF U10 (a) and Med-waves Hs (b) evaluation against satellite U10 (Jason-2) and satellite Hs (Jason-2 and SARAL) respectively: maps of metric values over the Mediterranean Sea subregions shown in Fig. 3, for a 1-year period (2014).

Figure 8b shows that the typical error (RMSE) varies from 0.18 m in the South-Central Levantine Basin (lev3 in Fig. 3) to 0.24 m in the Alboran Sea, the North West Mediterranean (nwm), and the North Adriatic subregions up to 0.29 m in the Atlantic subregion. In terms of SI, the highest value (0.28) is obtained in the North Adriatic Sea followed by the Aegean Sea (0.24). The South Adriatic, Alboran, Ligurian (tyr1), and East Levantine (lev4) seas also have relatively high SI values (0.21–0.23). The lowest values are found over the south-eastern Mediterranean Sea (0.15–0.16) and in the Atlantic subregion (0.15). SI and CORR have a similar pattern of variation, a notable difference being that the correlation coefficient obtains its worst value in the East Levantine and, in general, it has relatively lower values in the well-exposed regions of the Levantine Basin compared to the well-exposed regions to its west. In accordance with the above results, Ratsimandresy et al. (2008), examining model–satellite agreement over coastal locations of the western Mediterranean Sea, found the worst correlations in the Alboran Sea and east of Corsica. Bertotti et al. (2013), in a comparison of high-resolution wind and wave model output with satellite data over different subregions of the Mediterranean Sea, also found the largest scatter and lowest correlations in the Adriatic and the Aegean seas. In agreement, Zacharioudaki et al. (2015), focusing on the Greek seas, have shown a considerably larger scatter in the Aegean Sea than in the surrounding seas when model output was compared to satellite observations. As explained in the previous subsection (model–buoy comparison), it is difficult for wind models to reproduce orographic effects and/or local sea breezes well and difficult for wave models to resolve complicated bathymetry that introduces errors in these fetch-limited, enclosed regions, often characterized by a complex topography, well. Indeed, comparison with the equivalent results for the ECMWF wind speeds confirms these difficulties. For example, the pattern of SI and CORR variation for U10 largely resembles that for Hs, corroborating the conclusion of many studies that errors in wave height simulations by sophisticated wave models are mainly caused by errors in the generating wind fields (e.g. Komen et al., 1994; Ardhuin et al., 2007). Nevertheless, some differences do exist. For instance, the Hs SI in the Aegean Sea is relatively higher than the corresponding U10 SI. This is most probably because in this region of highly complicated bathymetry with many little islands the error of the wave model increases in relation to the error of the wind model. Similarly, in the East Levantine, Hs SI is lower than that implied by U10 SI. In this case, the wind model may not simulate local wind patterns, characterized by local sea breezes and easterly directions, well (Galil et al., 2006); however, the wave regime which is dominated by waves from the west sector (Galil et al., 2006) is better reproduced by the wave model. Negative bias and slope below unity are the case in all subregions except for the Atlantic Ocean and Alboran Sea.

Table 6Med-waves Tm evaluation against wave buoys' Hs, for the full Mediterranean Sea, for a 1-year period (2014) and seasonally.

## 4.2 Hindcast mean wave period

Table 6 presents the statistics of the comparison between the Med-waves hindcast mean wave period, Tm, and in situ observations of mean wave period, for the full Mediterranean Sea, for a 1-year period (2014) and seasonally. Figure 9 shows the corresponding QQ–scatter plot for the year-long statistics. It is shown that the model exhibits greater variability than the observations (SD in Table 6). RMSE varies from 0.8 s in summer to 1.07 s in winter. In relation to the mean of the observations, the error is about 17 %–19 %, with winter and spring being at the low end of this range and autumn at the high end. SI varies from 0.12 in winter and spring to 0.14 in summer and autumn. The non-trivial deviation of SI from relative RMSE (RMSE/$\stackrel{\mathrm{‾}}{R}$) indicates that a substantial part of the error is caused by bias. CORR has its minimum value (0.78) in summer and its maximum (0.87) in winter and spring. As before, these results indicate that the model wave period, like the model wave height, better follows the observations in well-defined wave conditions of higher waves and larger periods. Bias is negative with values that correspond to a model underestimate of about 11.5 %–13 %. Correspondingly, slope has a small variation of 0.87–0.89. Like for wave height, spring statistics are the most representative of the year-long statistics. Figure 9 clearly shows that the wave model underestimates the observations throughout the observed Tm range. Measurements of Tm<4.5 s are especially underestimated while those of relatively high Tm are better approximated by the model. As mentioned in Sect. 3, part of the model underestimation of observed Tm, especially over the lower Tm range, may be attributed to the absence of a high frequency cut-off in the model in contrast to the observations.

Figure 9QQ–scatter plots of Med-waves output versus wave buoy observations, for the full Mediterranean Sea, for a 1-year period (2014).

Table 7Med-waves Tm evaluation against wave buoy Tm, for each individual buoy location, for a 1-year period (2014).

Table 7 gives the statistics of the model–buoy comparison at the individual wave buoy locations. The typical error relative to the mean of the observations (RMSE/$\stackrel{\mathrm{‾}}{R}$) has its lowest values of 12 %–16 % over the western part of the Mediterranean Sea, west and south of France, with the two locations nearest to the Gibraltar Straight being at the low end of this range. Otherwise, this error is 17 %–23 %, reaching up to 29 % at location 61187 near the French–Italian border. At this location, all qualification metrics obtain their worst value. This is because wave buoy 61187 is located at a distance less than 2 km from coast and is affected by winds blowing from land. As already explained in the model–buoy wave height comparison, in this situation, the simulated fetch may differ substantially from the actual fetch because of the wave model grid size approximation; moreover, wind speed and wave height are considerably underestimated. RMSE is mainly caused by bias, which is negative at all locations. Thus, according to the RMSE, the relative bias is below 10 % over the western Mediterranean and is 14 %–20 % otherwise reaching up to 23 % at location 61187. It is only at location 61197, offshore from the eastern Balearic Islands, that the scatter of the data appears to contribute more to the typical error than the bias. This is a well-exposed offshore location where bias (2 %) and slope (0.97) have their best values and where model performance has been found to be optimal for wave height. The relatively high SI (0.15) and moderate correlation (0.86) at this location could be associated with the appearance of two density peaks in the density scatter plot (not shown), indicative of a double-peaked frequency spectrum. Density scatter plots with two peaks, although less distinct, have also been obtained for locations 61289 and 61021, offshore from France. In general, a close examination of the QQ–scatter plots (not shown) corresponding to the different locations has revealed that the model largely underestimates the observed Tm over the lower wave period range at all locations. Over the higher range, the model converges or overestimates the observed Tm in the western Mediterranean Sea, west and south of France. Otherwise, the model underestimates all observed Tm values with some convergence towards higher values. Slope mostly follows the pattern of variation in bias with values between 0.76 and 0.97. SI is relatively small with values between 0.09 (ATHOS) and 0.18 (61187) while CORR varies from 0.65 (61187) to 0.9 (61430, 68422). Generally, similar to the wave height results, the lowest correlations are found at coastal locations affected by fetch differences between model and reality due to a complex surrounding topography. Conversely, the highest correlations are obtained at the most exposed wave buoy locations.

## 4.3 Forecast skill

In the previous section, the performance of the Med-waves system has been characterized through the comparison of hindcast wave parameters with observations. In this section, the forecast skill of the Med-waves system is explored by comparing forecast wave parameters with observations at different forecast lead times. Hence, Fig. 10 shows Med-waves forecast skill for Hs (Fig. 10b, c) together with ECMWF forecast skill for U10 (Fig. 10a). The latter is evaluated against satellite observations, the former is evaluated against satellite (Fig. 10b) and buoy (Fig. 10c) observations. It is noted that in the model–buoy comparisons, each lead time represents a single point in time, whist, in the model–satellite comparisons, each lead time represents a full forecast day; i.e. +0 h represents forecast day 1 containing data from 0 to 24 h forecast. This approach is dictated by the scarcity of satellite observations in time.

Figure 10ECMWF U10 forecast skill evaluated against satellite observations (a) and Med-waves Hs forecast skill evaluated against satellite (b) and buoy (c) observations, for the full Mediterranean Sea, for a 1-year period (2014).

Figure 11Med-waves Tm forecast skill evaluated against buoy observations, for the full Mediterranean Sea, for a 1-year period (2014).

Figure 10 shows that Hs SI grows with forecast lead time at a constantly increasing growth rate. At the same time, CORR decreases with forecast lead time, with the decrease being more notable after the third day of forecast (+48 h). These patterns, which are consistent between model–buoy and model–satellite observations, mostly agree with the equivalent U10 patterns and manifest the deterioration of the forecast in time. A small difference between U10 and Hs forecast skill is the somewhat more linear increase in SI with forecast lead time in the first case, which results in a smaller overall U10 deterioration over the length of the forecast (14 %) compared to the respective Hs deterioration (19 % for model–satellite comparisons). This is indicative of the sensitivity of wave height to even limited variations in the input wind intensity. Conversely, waves seem to be less sensitive to wind misfits in time and space, which is manifested by the higher and more persistent Hs CORR over the forecast range compared to the respective U10 CORR. Contrary to SI and CORR but also to U10 bias, the evolution of Hs bias with forecast lead time in not monotonic. This apparent discrepancy between wind and wave bias evolution is attributed to the modification of the default values of the whitecapping dissipation coefficients in WAM, which, as shown in Sect. 4.1.2, have an impact on the bias of the wave model output relative to the observations. In any case, for both winds and waves, the variation in bias with forecast range is small and does not exceed 2 %.

Figure 11, like Fig. 10c, shows Med-waves forecast skill for Tm evaluated against wave buoy observations. Similar to Hs, SI increases with forecast range, CORR decreases, and bias exhibits a non-monotonic variation analogous to the one of Hs. In this case however, the variation in SI over the forecast period is small (5 %) compared to the respective Hs variation (25 %). This agrees with the finding that Tm errors are mainly caused by bias (Sect. 4.2).

5 Conclusions

The CMEMS Mediterranean wave forecasting system, Med-waves, has been operational since April 2016, providing short-range forecasts over the Mediterranean Sea at hourly intervals and at a horizontal resolution of 1∕24. The development and the evaluation of the performance of the system has been presented in detail in this paper. In the framework of this evaluation, the wave parameters of significant wave height and mean wave period have been evaluated against in situ and satellite observations over a period of 1 year (2014). Both hindcast quality and forecast skill have been assessed. In the former case, evaluation statistics have been provided for the Mediterranean Sea as a whole, at individual buoy locations and over predefined Mediterranean subregions. In the latter case, evaluation statistics have been provided only for the entire Mediterranean Sea. The main findings of this evaluation assessment are summarized below.

Overall, the significant wave height is accurately simulated by the model. Considering the Mediterranean Sea as a whole, the RMSE is 0.21 m and the bias is −0.03 m (3.7 %) when the model is compared to in situ observations and −0.06 m (5.5 %) when it is compared to satellite observations. In general, the model somewhat underestimates the observations for wave heights below 4 m whilst it mostly converges to the observations for higher waves. The scatter index, indicative of the scatter of the data around their regression line, is 25 % for the model–in situ comparison and 19 % for the model–satellite comparison, demonstrating a reduced scatter off the shore (where satellite measurements are mostly located) compared to near the shore (where in situ measurements are mostly obtained). The correlation coefficient is 0.95–0.96 and so is the best-fit slope. Model performance is better in winter when the wave conditions are well defined. Spatially, the model performs optimally at offshore wave buoy locations and well-exposed Mediterranean subregions. Within enclosed basins and near the coast, unresolved topography by the wind and wave models and fetch limitations cause the wave model performance to deteriorate. In particular, the model has an optimal performance along most of the southern Mediterranean Sea. Its performance is less good in the Alboran, Ligurian, Adriatic, Aegean, and Eastern Levantine seas, with the worst evaluation statistics obtained in the Adriatic. In terms of bias, the model overall underestimates the measurements in the Mediterranean Sea. The smallest underestimate is observed in the Aegean Sea while overestimate is observed in the Alboran Sea and in the Atlantic.

Poor wave model statistics were strongly linked to poor wind forcing statistics. Naturally, it is expected that an improved orographic representation of the ECMWF forecasting system will improve the quality of the surface wind fields near the coast while a higher temporal resolution of wind forcing would be beneficial to the model to resolve the high wind and wave variability in the Mediterranean Sea, providing more accurate wave fields. Moreover, the atmospheric model of ECMWF, which was not coupled with an ocean model (coupling was performed only with the wave component) by the time of this study, did not consider some vigorous air–sea interaction processes (large heat fluxes or strong ocean mixing processes) that occur in regions which are usually affected by extreme weather events such as the northern part of the Adriatic during strong wind events (Bora, Sirocco). As a result, it failed to properly reproduce the spatial structure of the wind fields. To overcome this limitation, many studies (Carniel et al., 2016; Ricchi et al., 2016, 2017) show that the use of a fully coupled atmosphere–ocean–wave model can be considered appropriate for these regions for properly representing the air–sea interactions and for providing a more realistic and consistent evolution of the atmospheric and oceanic fields. It is noted at this point that substantial progress has been made since the year 2014 – the year the results of this study are obtained for – by ECMWF regarding spatial and temporal resolutions and coupling. Specifically, the spatial resolution of the ECMWF winds has increased since spring 2016 from 16 to 9 km. Also, in June 2018, it was decided that ECMWF can provide hourly forecasts up to a forecast step of 90 h (Jean Bidlot, personal communication, 2018). Finally, in June 2018, the ECMWF forecasting system is a fully coupled atmosphere–waves–ocean–sea ice system (ECMWF, 2018). As a consequence of these improvements, a future validation of the Med-waves system is expected to yield better validation results.

The mean wave period is reasonably well simulated by the model. The RMSE is 0.7 s and is mainly caused by model bias, which has a value of −0.48 s (12 %). In general, the model underestimates the observed mean wave period and exhibits greater variability than the observations. A relatively larger model underestimate is found for mean wave periods below 4.5 s. The scatter index is 13 %, the correlation coefficient is 0.85 and the best-fit slope is 0.88. Model performance is a little better in winter when wave conditions are well defined. Spatially, the model somewhat overestimates the highest mean wave period values in the western Mediterranean Sea, west and south of France. Otherwise, model underestimate is widespread. Similar to the wave height, the model performance is best at well-exposed offshore locations and deteriorates near the shore mainly due to fetch limitations.

The forecast skill of the model over the Mediterranean Sea deteriorates with forecast range. The growth of error in the wave forecast is mainly due to the growth of error in the forcing wind fields. The scatter index of the significant wave height deteriorates by 19 and 25 % over the 5-day forecast for model–satellite and model–buoy comparisons respectively. The equivalent deterioration for mean wave period is only 5 % (model–buoy comparison). A monotonic decrease in correlation is also observed. On the contrary, the evolution of bias with forecast range shows some variability with no clear trend. Nevertheless, this variability does not exceed 3 % over the forecast period.

In the next version of the system an optimal interpolation data assimilation scheme is added to the Med-waves system in order to blend satellite along-track significant wave height measurements with model background forecasts. Although wave data assimilation is known not to be particularly beneficial in areas where wind sea conditions are dominant, we expect that wave forecasts in certain sub-areas of the Mediterranean Sea where swell propagation is quite frequent, will be improved at +24 h and perhaps +48 h lead time. The enhanced Med-waves system with the data assimilation system module is going to produce 3-hourly wave analyses on a daily basis for the Mediterranean Sea by assimilating Sentinel-3 and Jason-3 altimeter measured significant wave heights and surface winds. The assimilation is based on the inherent data assimilation scheme of WAM Cycle 4.5.4 model, which generates an updated wave field by distributing the information from the observed significant wave height and surface wind data within a given time window over the entire model grid. The Med-waves data assimilation component is integrated into the Med-waves system since April 2018.

Lastly, more work will be devoted to improve the offline coupling between waves and currents by including as a next step the variations in the sea level as predicted by the physical component of the Med MFC system. In parallel the full online two-way coupling between Mediterranean waves and currents will be developed and implemented into the Med MFC system, with a target to enter into the operational chain in future versions of the Copernicus Marine Service, improving the forecasting skill of the models in various coastal areas of the basin.

Data availability
Data availability.

The in situ wave buoy observations used in this study have been obtained from the Copernicus Marine Environment Monitoring Service (CMEMS) IN-SITU Thematic Assembly Centre (INS TAC) archive and are available from http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=INSITU_MED_NRT_OBSERVATIONS_013_035 (CMEMS, 2018a). The satellite observations have been obtained from a merged altimeter wave height database set-up at CERSAT – IFREMER, France, and are available from ftp://ftp.ifremer.fr/ifremer/cersat/products/swath/altimeters/waves/data/ (CERSAT-IFREMER, 2017). The model outputs for the year 2014 are available upon request from the authors. Model outputs since 2016 are available through CMEMS from http://marine.copernicus.eu/services-portfolio/access-to-products/ (CMEMS, 2018b).

Author contributions
Author contributions.

GK and MR decided on the set-up and initial configuration of the wave forecasting system and designed a number of sensitivity runs to help reach a conclusion on the final model configuration. MR collected the necessary model inputs and performed all simulations. AZ performed all validation so as to decide together with the co-authors on the final version of the wave model. She also prepared the paper with contributions from all co-authors. GK coordinated all activities.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This work has been supported by HCMR (https://www.hcmr.gr/en/, last access: September 2018) and Copernicus Marine Environment Monitoring Service (CMEMS) (http://marine.copernicus.eu/, last access: September 2018). CMEMS is implemented by Mercator Ocean through a delegation agreement with the European Union. The authors want to acknowledge the CMEMS for providing ocean current data and in situ observations, the Italian Meteorological Service for providing ECMWF wind data, and CERSAT – IFREMER for the satellite observations for validating the system.

Edited by: Piero Lionello
Reviewed by: Jean-Raymond Bidlot and one anonymous referee

References

Alomar, M., Sánchez-Arcilla, A., Bolaños, R., and Sairouni, A.: Wave growth and forecasting in variable, semi-inclosed domains, Cont. Shelf Res., 87, 1–6, 2014.

Ardhuin, F., Bertotti, L., Bidlot, J. R., Cavaleri, L., Filipetto, V., Lefevre, J. M., and Wittmann, P.: Comparison of wind and wave measurements and models in the Western Mediterranean Sea, Ocean Eng., 34, 526–541, https://doi.org/10.1016/j.oceaneng.2006.02.008, 2007.

Artale, V., Calmanti, S., Carillo, A., Dell'Aquila, A., Herrmann, M., Pisacane, G., Ruti, P., Sannino, G., and Struglia, M.: An atmosphere–ocean regional climate model for the Mediterranean area: assessment of a present climate simulation, Ocean Model., 30, 56–72, 2009.

Bentamy, A., Ayina, H.-L., Queffeulou, P., Croize-Fillon, D., and Kerbaol, V.: Improved near real time surface wind resolution over the Mediterranean Sea, Ocean Sci., 3, 259–271, https://doi.org/10.5194/os-3-259-2007, 2007.

Bertotti, L. and Cavaleri, L.: Wind and Wave Predictions in the Adriatic Sea, J. Mar. Syst., 78, S227–S234, 2009a.

Bertotti, L. and Cavaleri, L.: Large and small scale wave forecast in the Mediterranean Sea, Nat. Hazards Earth Syst. Sci., 9, 779–788, https://doi.org/10.5194/nhess-9-779-2009, 2009b.

Bertotti, L., Canestrelli, P., Cavaleri, L., Pastore, F., and Zampato, L.: The Henetus wave forecast system in the Adriatic, Nat. Hazards Earth Syst. Sci., 11, 2965–2979, https://doi.org/10.5194/nhess-11-2965-2011, 2011.

Bertotti, L., Cavaleri, L., Loffredo, L., and Torrisi, L.: Nettuno: Analysis of a Wind and Wave Forecast System for the Mediterranean Sea, Mon. Weather Rev., 141, 3130–3141, https://doi.org/10.1175/MWR-D-12-00361.1, 2013.

Bidlot, J. R.: Inter-comparison of operational wave forecasting systems, in: Proceedings of the 10th International Workshop on Wave Hindcasting and Forecasting and Coastal Hazard Symposium, available at: http://www.waveworkshop.org (last access: June 2017), 2007.

Bidlot, J. R.: Intercomparison of operational wave forecasting systems against buoys: data from ECMWF, MetOffice, FNMOC, MSC, NCEP, MeteoFrance, DWD, BoM, SHOM, JMA, KMA, Puerto del Estado, DMI, CNR-AM, METNO, SHN-SM, JCOMM Technical Report, SPA_ETWS_verification201401_201412, available at: https://www.jcomm.info/index.php?option=com_oe&task=viewDocumentRecord&docID=14526 (last access: September 2018), 2015.

Bidlot, J. R.: Twenty-one years of wave forecast verification, ECMWF Newsletter, 150, 31–36, 2017.

Bolaños, R., Sanchez-Arcilla, A., Gómez, J., and Cateura, J.: Limits of operational wave prediction in the North-Western Mediterranean, in: Proceedings of the 29th International Conference Coastal Engineering 2004, Vol. 4, April 2005, Lisbon, Portugal, 818–829, 2005.

Bolaños-Sanchez, R., Sanchez-Arcilla, A., and Cateura, J.: Evaluation of two atmospheric models for wind–wave modelling in the NW Mediterranean, J. Mar. Syst., 65, 336–353, 2007.

Booij, N., Ris, R. C., and Holthuijsen, L. H.: A third-generation wave model for coastal regions, 1, model description and validation, J. Geophys. Res., 104, 7649–7666, 1999.

Breivik, O., Swail, V., Babanin, A., and Horsburgh, K.: The International Workshop on Wave Hindcasting and Forecasting and the Coastal Hazards Symposium, arXiv:1503.00847v1, 2015.

Carniel, S., Benetazzo, A., Bonaldo, D., Falcieri, M. F., Miglietta, M. M., Ricchi, A., and Sclavo, M.: Scratching beneath the surface while coupling atmosphere, ocean and waves: Analysis of a dense water formation event, Ocean Model., 101, 101–112, https://doi.org/10.1016/j.ocemod.2016.03.007, 2016.

Cavaleri, L. and Bertotti, L.: The characteristics of wind and wave fields modelled with different resolutions, Q. J. Roy. Meteorol. Soc., 129, 1647–1662, 2003.

Cavaleri, L. and Bertotti, L.: Accuracy of the modelled wind and wave fields in enclosed seas, Tellus A, 56, 167–175, 2004.

Cavaleri, L. and Sclavo, M.: The calibration of wind and wave model data in the Mediterranean Sea, Coast. Eng., 53, 613–627, 2006.

Cavaleri, L., Bertotti, L., and Lionello, P.: Wind wave cast in the Mediterranean Sea, J. Geophys. Res., 96, 10739–10764, https://doi.org/10.1029/91JC00322, 1991.

Cavaleri, L., Abdalla, S., Benetazzo, A., Bertotti, L., Bidlot, J.-R., Breivik, Ø., Carniel, S., Jensen, R. E., Portilla-Yandun, J., Rogers, W. E., Roland, A., Sanchez-Arcilla, A., Smith, J.-M., Staneva, J., Toledo, Y., van Vledder, G. P., and van der Westhuysen, A. J.: Wave modelling in coastal and inner seas, Prog. Oceanogr., 167, 164–233, https://doi.org/10.1016/j.pocean.2018.03.010, 2018.

CERSAT-IFREMER: Dataset, available at: ftp://ftp.ifremer.fr/ifremer/cersat/products/swath/altimeters/waves/data/, last access: June 2017.

Chen, S. S., Zhao, W., Donelan, M. A., Price, J. F., and Walsh, E. J.: The CBLAST-hurricane program and the next-generation fully coupled atmosphere–wave–ocean models for hurricane research and prediction, B. Am. Meteorol. Soc., 88, 311–317, https://doi.org/10.1175/BAMS-88-3-311, 2007.

Clementi, E., Oddo, P., Korres, G., Drudi, M., and Pinardi, N.: Coupled wave–ocean modelling system in the Mediterranean Sea, in: 13th International Workshop on Wave Hindcasting and Forecasting and 4th Coastal Hazards Symposium, At Banff, Canada, 2013.

Clementi, E., Oddo, P., Drudi, M., Pinardi, N., Korres, G., and Grandi, A.: Coupling hydrodynamic and wave models: first step and sensitivity experiments in the Mediterranean Sea, Ocean Dynam., 67, 1293–1312, https://doi.org/10.1007/s10236-017-1087-7, 2017.

CMEMS: QUality Information Document For Global Sea Physical Ananalysis and Forecasting Product, Ref: CMEMS-GLO-QUID-001-024, available at: http://marine.copernicus.eu/documents/QUID/CMEMS-GLO-QUID-001-024.pdf (last access: April 2018), 2016a.

CMEMS: QUality Information Document for Med Currents analysis and forecast product, Ref: CMEMS-Med-QUID-006-001-V2, available at: http://marine.copernicus.eu/documents/QUID/CMEMS-MED-QUID-006-001.pdf (last access: December 2017), 2016b.

CMEMS: Dataset: INSITU_MED_NRT_OBSERVATIONS_013_035, available at: http://marine.copernicus.eu/services-portfolio/30access-to-products/?option (last access: September 2018), 2018a.

CMEMS: Dataset: MEDSEA_ANALYSIS_FORECAST_WAV_006_017, available at: http://marine.copernicus.eu/services-portfolio/access-to-products/?option (last access: September 2018), 2018b.

Donatini, L., Lupieri, G., Contento, G., Feudale, L., Pedroncini, A., Cusati, L., and Crosta, A.: A high resolution wind and wave forecast model chain for the Mediterranean – Adriatic Sea, in: Towards Green Marine Technology and Transport, CRC Press, London, UK, 859–866, 2015.

ECMWF: ifs.documentation CY41R1. Part VII: ECMWF Wave – Model documentation 2015, available at: https://www.ecmwf.int/en/forecasts/documentation-and-support/changes-ecmwf-model/ifs-documentation (last access: March 2018), 2015.

ECMWF: ifs.documentation CY43R3, available at: https://www.ecmwf.int/en/forecasts/documentation-and-support/changes-ecmwf-model/ifs-documentation (last access: March 2018), 2017.

Galil, B., Herut, B., Rosen, D. S., and Rosentroup, Z.: A Consise Physical, Chemical and Biological Characterizasion of Eastern Mediterranean with Emphasis on the Israeli coast, IOLR Report H07, Israel Oceanographic and Limnological Research LTD, Israel, 2006.

GEBCO: GEBCO 30 arc-second grid, available at: http://www.gebco.net/data_and_products/gridded_bathymetry_data/gebco_30_second_grid/ (last access: December 2017), 2016.

Gerling, T. W.: Partitioning Sequences and Arrays of Directional Ocean Wave Spectra into Component Wave Systems, J. Atmos. Ocean. Tech., 9, 444–458, https://doi.org/10.1175/1520-0426(1992)009<0444:PSAAOD>2.0.CO;2, 1992.

Günther, H. and Behrens, A.: The WAM model. Validation document Version 4.5.4, Intitute of Coastal Research Helmholtz-Zentrum Geesthacht (HZG), Geesthacht, Germany, 2012.

Haiden, T., Janousek, M., Bidlot, J., Ferranti, L., Prates, F., Vitart, F., Bauer, P., and Richardson, D. S.: Evaluation of ECMWF forecasts, including the 2016 resolution upgrade, ECMWF Technical Memoranda, ECMWF, England, 2016.

Hasselmann, K., Barnett, T. P., Bouws, E., Carlson, H., Cartwright, D. E., Enke, K., Ewing, J. A., Gienapp, H., Hasselmann, D., Kruseman, P., Meerburg, A., Müller, P., Olbers, D. J., Richter, K., Sell, W., and Walden, H.: Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP), Deutches Hydrographisches Institut, Hamburg, Germany, 1973.

Hasselmann, S. and Hasselmann, K.: Computations and Parameterizations of the Nonlinear Energy Transfer in a Gravity-Wave Spectrum. Part I: A New Method for Efficient Computations of the Exact Nonlinear Transfer Integral, J. Phys. Oceanogr., 15, 1369–1977, 1985.

Hasselmann, S., Hasselmann, K., Allender, J., and Barnett, T.: Computations and Parameterizations of the Nonlinear Energy Transfer in a Gravity-Wave Specturm. Part II: Parameterizations of the Nonlinear Energy Transfer for Application in Wave Models, J. Phys. Oceanogr., 15, 1378–1391, 1985.

Holthuijsen, L. H., Booji, N., and Bertotti, L.: The propagation of wind errors through ocean wave hindcasts, J. Offshore Mech. Arct. Eng., 118, 184–189, 1996.

Janssen, P.: Wave-Induced Stress and the Drag of Air Flow over Sea Waves, J. Phys. Oceanogr., 19, 745–754, 1989.

Janssen, P.: Quasi-linear Theory of Wind-Wave Generation Applied to Wave Forecasting, J. Phys. Oceanogr., 21, 1631–1642, 1991.

Katsafados, P., Papadopoulos, A., Korres, G., and Varlas, G.: A fully coupled atmosphere–ocean wave modeling system for the Mediterranean Sea: Interactions and sensitivity to the resolved scales and mechanisms, Geosci. Model Dev., 9, 161–173, https://doi.org/10.5194/gmd-9-161-2016, 2016.

Komen, G. J., Cavaleri, L., Donelan, M., Hasselmann, K., Hasselmann, S., and Janssen, P.: Dynamics and modelling of ocean waves, Cambridge University Press, Cambridge, 1994.

Licer, M., Smerkol, P., Fettich, A., Ravdas, M., Papapostolou, A., Mantziafou, A., Strajnar, B., Cedilnik, J., Jeromel, M., Jerman, J., Petan, S., Malacic, V., and Sofianos, S.: Modeling the ocean and atmosphere during an extreme bora event in northern Adriatic using one-way and two-way atmosphere–ocean coupling, Ocean Sci., 12, 71–86, https://doi.org/10.5194/os-12-71-2016, 2016.

Lionello, P., Cogo, S., Galati, M. B., and Sanna, A.: The Mediterranean surface wave climate inferred from future scerario simulations, Global Planet. Change, 63, 152–162, 2008.

Madec, G.: NEMO ocean engine, Note du Pole de modelisation, Institut Pierre-Simon Laplace (IPSL), France, 2008.

Mao, M. and Xia, M.: Dynamics of wave-current-surge interactions in Lake Michigan: A model comparison, Ocean Model., 110, 1–20, https://doi.org/10.1016/j.ocemod.2016.12.007, 2017.

Moeini, M. H., Etemad-Shahidi, A., and Chegini, V.: Wave modeling and extreme value analysis off the northern coast of the Persian Gulf, Appl. Ocean Res., 32, 209–218, 2010.

NGDC: 2-minute Gridded Global Relief Data (ETOPO2) v2, National Geophysical Data Centre, NOAA, available at: https://www.ngdc.noaa.gov/mgg/global/etopo2.html (last access: November 2016), 2006.

Osuna, P. and Wolf, J.: A numerical study of the effect of wave-current interaction process in the hydrodynamics of the Irish Sea Ocean Waves Measurements and Analysis, in: Fifth International Symposium WAVES 2005, 3–7 July 2005, Madrid, Spain, 2005.

Ponce, S. and Ocampo-Torres, F. J.: Sensitivity of a wave model to wind variability, J. Geophys. Res., 103, 3179–3201, 1998.

Pullen, J., Doyle, J. D., and Signell, R. P.: Two-Way Air–Sea Coupling: A Study of the Adriatic, Mon. Weather Rev., 134, 1465–1483, https://doi.org/10.1175/MWR3137.1, 2006.

Queffeulou, P. and Croizé-Fillon, D.: Glogal altimeter SWH data set – May 2013, Technical Report, available at: ftp://ftp.ifremer.fr/ifremer/cersat/products/swath/altimeters/waves/documentation/ (last access: December 2016), 2013.

Ratsimandresy, A. W., Sotillo, M. G., Carretero Albiach, J. C., Álvarez Fanjul, E., and Hajji, H.: A 44-year high-resolution ocean and atmospheric hindcast for the Mediterranean Basin developed within the HIPOCAS Project, Coast. Eng., 55, 827–842, https://doi.org/10.1016/j.coastaleng.2008.02.025, 2008.

Renault, L., Chiggiato, J., Warner, J. C., Gomez, M., Vizoso, G., and Tintoré, J.: Coupled atmosphere–ocean–wave simulations of a storm event over the Gulf of Lion and Balearic Sea, J. Geophys. Res.-Oceans, 117, C09019, https://doi.org/10.1029/2012JC007924, 2012.

Ricchi, A., Miglietta, M. M., Falco, P. P., Benetazzo, A., Bonaldo, D., Bergamasco, A., Sclavo, M., and Carniel, S.: On the use of a coupled ocean–atmosphere–wave model during an extreme Cold Air Outbreak over the Adriatic Sea, Atmos. Res., 172–173, 48–65, https://doi.org/10.1016/j.atmosres.2015.12.023, 2016.

Ricchi, A., Miglietta, M. M., Barbariol, F., Benetazzo, A., Bergamasco, A., Bonaldo, D., Cassardo, C., Falcieri, F. M., Modugno, G., Russo, A., Sclavo, M., and Carniel, S.: Sensitivity of a Mediterranean tropical-like cyclone to different model configurations and coupling strategies, Atmosphere, 8, 92, https://doi.org/10.3390/atmos8050092, 2017.

Saket, A, Etemad-Shahidi, A., and Moeini, M. H.: Evaluation of ECMWF wind data for wave hindcast in Chabahar zone, J. Coast. Res., 1, 380–385, 2013.

Saulter, A.: Proposal for wave verification within a Marine Core Service, MyWave Report D4.4, October 2014.

Sepulveda, H. H., Queffeulou, P., and Ardhuin, F.: Assessment of Saral/AltiKa wave height measurements relative to buoy, Jason-2, and Cryosat-2 Data, Mar. Geod., 38, 449–465, https://doi.org/10.1080/01490419.2014.1000470, 2015.

Signell, R. P. S., Carniel, S., Cavaleri, L., Chiggiato, J., Doyle, J. D., Pullen, J., and Scavo, M.: Assessment of wind quality for oceanographic modelling in semi-closed basins, J. Mar. Syst., 53, 217–233, 2005.

Soloviev, A. V., Lukas, R., Donelan, M. A., Haus, B. K., and Ginis, I.: The air–sea interface and surface stress under tropical cyclones, Nat. Scient. Rep., 4, 5306, https://doi.org/10.1038/srep05306, 2014.

Teixeira, J. C., Abreu, M. P., and Guedes Soares, C.: Uncertainty of ocean wave hindcasts due to wind modelling, J. Offshore Mech. Arct. Eng., 117, 294–297, 1995.

Tolman, H. L.: A third-generation model for wind waves on slowly varying, unsteadyand inhomogeneous depth and currents, J. Phys. Oceanogr., 21, 782–797, 1991.

Tolman, H. L.: Effects of numerics on the physics in a third generation wind-wave model, J. Phys. Oceanogr., 22, 1095–1111, 1992.

Tolman, H. L.: User manual and system documentation of WAVEWATCH-III version 1.18, NOAA/NWS/NCEP/OMB Technical Note 166, p. 110, available at: http://polar.ncep.noaa.gov/mmab/papers/tn166/OMB_166.pdf (last access: October 2018), 1999.

Varlas, G., Katsafados, P., Papadopoulos, A., and Korres, G.: Implementation of a two-way coupled atmosphere-ocean wave modeling system for assessing air–sea interaction over the Mediterranean Sea, Atmos. Res., 208, 201–217, https://doi.org/10.1016/j.atmosres.2017.08.019, 2018.

Wahle, K., Staneva, J., Koch, W., Fenoglio-Marc, L., Ho-Hagemann, H. T., and Stanev, E.: An atmosphere-wave regional coupled model: improving predictions of wave heights and surface waves in the Southern North Sea, Ocean Sci., 13, 289–301, https://doi.org/10.5194/os-13-289-2017, 2017.

WAMDI Group: The WAM Model – A Third Generation Ocean Wave Prediction Model, J. Phys. Oceanogr., 18, 1775–1810, https://doi.org/10.1175/1520-0485(1988)018<1775:TWMTGO>2.0.CO;2, 1988.

Zacharioudaki, A., Korres, G., and Perivoliotis, L.: Wave climate of the Hellenic Seas obtained from a wave hindcast for the period 1960–2001, Ocean Dynam., 65, 795–816, https://doi.org/10.1007/s10236-015-0840-z, 2015.