3.1. Wind Fields
The National Oceanic and Atmospheric Administration (NOAA) and United States Navy government website provides several datasets, that have been widely be used as model wind data inputs in the several oil spill test-cases. Firstly, at global scale, the NOAA and the National Center for Environmental Prediction (NCEP) supports the Climate Forecast System Reanalysis (CFSR) model, which was created and implemented as a worldwide, high-resolution, linked atmosphere-ocean-land surface-sea-ice system to properly predict the conditions of these coupled domains during a 32-year period (January 1979—March 2011) [
75,
76]. The CFSR has a spatial horizontal resolution of 0.5° (~56 km) with hourly time step (
https://www.ncei.noaa.gov/products/weather-climate-models/climate-forecast-system). The CFSR has been applied for oil spill modeling studies by French-McCay et al. [
66] and Meza-Padill et al. [
77]. Moreover, NOAA and the US Navy provide meteorological model outputs of the Navy Operational Global Atmospheric Prediction System (NOGAPS) with horizontal resolution of 0.5
o (~56 km) and temporal resolution of 6 hours, globally (ttps://
www.ncdc.noaa.gov/data-access/model-data/model-datasets/navyoperational-global-atmospheric-prediction-system). The NOGAPS has been integrated in several oil spill test cases, as in King et al. [
78], Brushett et al. [
79], Le Hénaff et al. [
80], Vaz et al. [
81] and French-McCay et al. [
66]. Similarly, the NCEP (National Centers for Environmental Prediction) provides atmospheric data from the GFS (Global Forecasting System) for dozens of atmospheric and land-soil variables, including water temperature, winds, precipitation, soil moisture, and atmospheric ozone concentration [
82]. NCEP-GRF produces forecasts at three spatial resolutions of 0.25°, 0.5°, and 1° (
https://www.nco.ncep.noaa.gov/pmb/products/gfs/) covering the whole globe [
83]. Most oil spill models have been forced with the 0.25° horizontal resolution [
44,
50,
70,
83]. The temporal resolution of GFS is 3 hours and the NCEP contains wind velocity 10 meters above sea level, for the entire planet from May 6, 2011 to the present [
84].
At regional scales, the National Center for Atmospheric Research (NCAR) and the National Centers for Environmental Prediction (NCEP) maintains the Weather Research and Forecasting (WRF) Model, which is a cutting-edge mesoscale numerical weather prediction system, intended for both atmospheric research and operational forecasting. In oil spill modeling, the WRF has been implemented by Zacharias et al. [
44] in a region from 20° S, 50° W to 10° N, 20° W, with 1-hour time interval and 0.15° horizontal resolution [
85]. Moreover, the NOAA NCEP system (
https://www.ncei.noaa.gov/products/weather-climate-models/north-american-mesoscale) provides atmospheric forecasts for the North America through the North American Mesoscale Forecast System (NAM). NAM creates different grids (or domains) of weather forecasts with varying horizontal resolutions [
86]. Temperature, precipitation, light intensity, and turbulent kinetic energy are just a few of the weather elements estimated at each grid cell. NAM employs additional numerical weather models to create high-resolution predictions over fixed regions and, on occasion, to track major weather events, such as hurricanes. In oil spill modeling NAM has been applied at 12 km horizontal resolution, with 1-hour time step, by French-McCay et al. [
66]. In addition, the NOAA NCEP provides the North American Regional Reanalysis (NARR) system for weather reanalysis, (
http://www.emc.ncep.noaa.gov/mmb/rreanl/) having 3-hr time step and 0.3° (~ 32 km) spatial resolution over the North America [
87,
88]. NARR project has been applied in oil spill modeling by French-McCay et al. [
66,
89]. Moreover, NARR system is based on a version of the Eta Model and its 3D variational data assimilation system (EDAS) has been operational since April 2003 [
90].
In addition, the short-term model results produced by the European Centre for Medium-Range Weather Forecasts [
92] have been widely used as forcing in oil spill modelling. ECMWF provides reliable daily global atmospheric forcing with three-hourly winds and spatial resolution of 0.125° (approx. 27 km). Specifically, ERA5 contains wind forcing reanalysis data (
https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5) and is derived from fifth-generation ECMWF atmospheric reanalysis of the global climate, which integrates multi-source measurements with numerical simulations, using an assimilation model. This dataset has been produced using the 4D-Var data assimilation scheme and model forecasts in CY41R2 of the ECMWF Integrated Forecast System (IFS) [
48]. It has a high temporal-spatial resolution (1 hour - 0.25°) and a long-time span from January 1, 1979 to the present [
93]. Data can be obtained by visiting
https://cds.climate.copernicus.eu/. Recent model upgrades have improved the overall performance of the forecasting system throughout the medium range. Further details on model description and verification can be found in the works of Ehard et al. [
94], Haiden et al. [
95] and Hersbach et al. [
96]. ERA5 has been used as wind boundary forcing in various oil spill scenarios, as in Zhang et al. [
83], Abdallah & Chantsev [
69,
71], Davis Morales et al. [
51] and Liu et al. [
42]. In parallel, in case of simulating retrospective oil spills, the ERA-Interim dataset could be used [
97], a global reanalysis data product covering the data-rich period since 1979. Originally, ERA-Interim was run from 1989, but the 10-year extension for 1979-1988 was produced in 2011, providing data every 6 hours, with a 1/8° spatial resolution [
57,
98] [
99]. ERA-Interim has been applied in oil spill simulations in several test cases all over the world [
46,
47,
48,
49,
53,
54,
56,
57]. Moreover, Kampouris et al. [
55] has used the ECMWF ensemble prediction system at ∼ 9 km and ~18 km horizontal resolution for wind forcing.
The Eta/NCEP model [
100,
101] has been in operational use at the Hellenic National Meteorological Service and at the University of Athens (
http://forecast.uoa.gr). Moreover, the high frequency winds from the SKIRON non-hydrostatic forecasting model [
102], with 5 and 10 km spatial resolution has been utilized during real oil spill pollution [
58,
62] and in numerous oil spill test-cases, such as in Zafirakou-Koulouris et al. [
103], Ribotti et al. [
104], Zodiatis et al. [
60,
63], Goldman et al. [
105], De Dominicis et al. [
19] and Sepp-Neves et al. [
53], both in the Mediterranean and the Black Sea. In parallel, the HCMR (Hellenic Centre for Marine Research) provides meteorological forecasts via the POSEIDON weather forecasting system [
106], also based on SKIRON/Eta model [
102], which covers an area broader than the Mediterranean basin, with a horizontal resolution of ~5 km. POSEIDON has been coupled with oil spill models, as in Annika et al. [
107], and Zodiatis et al. [
108] . Finally, The University of Malta (UoM) provides meteorological forecasts through the MALTA Maria ETA Model (
http://www.capemalta.net/maria/regional/results.html) with a horizontal resolution of ~4 km, covering the Central Mediterranean Sea and the Maltese Islands [
109]. The model has temporal resolution of 3 hours, providing forecasts for only 1 day in advance. MARIA/Eta High Resolution Atmospheric Forecasting System is based on the atmospheric limited area NCEP/Eta model [
101,
110] and it has been applied in oil spill case studies, like in Drago et al. [
109].
The Spanish Met Office, AEMET (Spanish State Meteorological Agency) (
https://www.aemet.es/es/portada), produces meteorological forcing forecasts using the HIRLAM (High Resolution Limited Area Model) [
111,
112]. This forecasting system runs with a horizontal resolution at 1/7° (~15 km) over the whole Western Mediterranean, providing hourly data every 6 hours, for the next 72-hours [
113]. HIRLAM has been coupled to the oil spill model TESEO, and has been applied in the Prestige oil spill accidental release in the Bay of Biscay [
113,
114]. Météo-France contributes with atmospheric data through the ARPEGE model (Action de Recherche Petite Echelle Grande Echelle) for the entire Mediterranean basin (
http://www.umr-cnrm.fr/spip.php?article121&lang=en). ARPEGE is a numerical model for global general circulation. Météo-France developed it in collaboration with ECMWF (Reading, UK) for operational numerical weather forecasting [
115,
116,
117]. The ARPEGE model has incorporated the four-dimensional variational assimilation (4D-Var). The spatial resolution of the ARPEGE model is ~10 km in the Mediterranean Sea and the temporal resolution is 3 hours [
118]. Recently, the model was upgraded in its vertical grid, composing of 105 levels, with horizontal grid of ~5 km over Europe and 24 km elsewhere [
119]. This fine resolution 5 km edition has not yet been applied to oil spill modeling. Oil spill models have been coupled only to the 10 km resolution version.
Figure 1.
General concept of Operational oil spill modeling.
Figure 1.
General concept of Operational oil spill modeling.
Table 1.
Atmospheric models with the corresponding domains and horizontal resolutions used in oil spill modeling.
Table 1.
Atmospheric models with the corresponding domains and horizontal resolutions used in oil spill modeling.
Wind |
Provider |
Geographical area |
Spatial Resolution |
Data Type |
Reference |
Poseidon |
HCMR |
Mediterranean |
~5 km |
Forecast |
[106] |
HIRLAM |
AEMET |
Western Mediterranean |
~ 15 km |
Forecast |
[108,112,113,114] |
ARPEGE |
Meteo-France |
Mediterranean |
~10 km |
Forecast |
[108,118] |
SKIRON |
UOA |
Mediterranean and Black Sea |
~5 and 10 km |
Forecast |
[53,102,120] [19,58,61,103,104,121] |
MALTA/Maria ETA model |
UOM |
Central Mediterranean |
~4 km |
Forecast |
[108,109] |
NAM |
NOAA/NCEP |
North America |
12km |
Forecast |
[66,86] |
NARR |
NOAA/NCEP |
North America |
0.3° (32 km) |
Reanalysis |
[66,89,90] |
NCOM AMSEAS |
NOAA/FNMOC |
Gulf of Mexico and Caribbean |
1/36° (~ 3km) |
Hindcast |
[74] |
WRF |
NCAR/NCEP |
Regional |
0.15° (~16 km) |
Forecast |
[44] |
NOGAPS |
NOAA/United States Navy |
Global |
0.5° (~56 km) |
Forecast |
[66,81] |
CFSR |
NOAA/NCEP |
Global |
0.5° (~56 km) |
Reanalysis |
[66,77] |
GFS |
NOAA/NCEP |
Global |
0.25° (~27 km) |
Forecast |
[50,69,70,71,83] |
ERA5 |
ECMWF |
Global |
0.25° (~27 km) |
Reanalysis |
[42,51] |
Era-Interim |
ECMWF |
Global |
0.125° (~12.5 km) |
Reanalysis |
[46,47,48,49,53,54,56,57] |
3.2. Hydrodynamics
The Copernicus Marine Environmental Monitoring Service (CMEMS) provides several hydrodynamic datasets at global and six regional seas. In the present study, only the data products relevant to oil spill modeling will be discussed. Firstly, the Global Ocean 1/12° Physics Analysis and Forecast model, provides daily and monthly-mean data for sea temperature, salinity, currents, sea level, mixed layer depth and ice parameters, from the top to the bottom of the global ocean [
122]. In addition, it produces the hourly-mean surface fields for sea level height, temperature and currents. The global ocean output files have a horizontal resolution of 1/12° (~ 9 km) and a regular longitude/latitude equirectangular projection. This dataset has been widely applied in oil spill simulations, like in the studies of Devis Morales et al. [
51], Sepp Neves et al. [
53] and Siqueira et al. [
57]. Moreover, CMEMS provides the dataset Mediterranean Sea Physics Analysis and Forecast (MEDSEA_ANALYSISFORECAST_PHY_006_013) [
123] which is being produced from a coupled hydrodynamic-wave model, implemented over the entire Mediterranean Basin. It consists the physical component of the Mediterranean Forecasting System (Med-Currents), with a horizontal grid resolution of 1/24° (approximately 4 km) and it has 141 unevenly spaced vertical levels. This dataset has been implemented in oil spill boundary forcing, e.g., in Liubartseva et al. [
54], Kampouris et al. [
55] and Keramea et al. [
50]. The hydrodynamics are provided by the Nucleus for European Modelling of the Ocean (NEMO v3.6). The model solutions are corrected by a variational data assimilation scheme (3DVAR) of temperature and salinity vertical profiles, as well as along track satellite Sea Level Anomaly observations [
124]. Finally, CMEMS supports the GLO-CPL dataset (Global Ocean 1/4° physics analysis and prediction) which is a data assimilation and forecast system that provides 10 days of 3D global ocean forecasts at ~27 km spatial resolution. The system employs the Met Office Unified Model v10.6 atmosphere configuration (at 40 km resolution) that is hourly coupled to NEMO v3.4 [
125] ocean configuration and the CICE v4.1 multi-thickness category sea ice model (both on the ORCA025 grid) [
126]. The GLO-CPL dataset has been used as forcing input in the GNOME oil spill simulations of Abdallah & Chantsev [
71] for the Red Sea.
The NOAA National Ocean Service (NOAA/NOS) Coast Survey Development Laboratory (CSDL) runs the NOS GOM Nowcast/Forecast Model (NGOM) [
127], which is an application of the POM model [
128] in the Gulf of Mexico. Moreover, the spatial resolution of NGOM is 5–6 km in the northeastern and central GoM, with 37 levels in the vertical (
https://www.bco-dmo.org/dataset/831523). Furthermore, the forecasts are obtained every 3 hours. NGOM has been used as forcing data in oil spill simulations, as in the case of the Deepwater Horizon buoyant plume simulation in combination with the OILMAP DEEP model [
66,
67]. In parallel, the NOAA and FNMOC (Fleet Numerical Meteorology and Oceanography Center) provides operational ocean predictions using the Navy Coastal Ocean Model (NCOM), with horizontal resolution of 1/36° (~3km) and 40 levels in the vertical. The model is capable to produce 4-day forecasts at 3-hour time steps. French-McCay et al. [
67] and Grubesic et al. [
74] have applied NCOM in oil spill simulations of OIL-MAPDEEP and BLOSOM, respectively. The AMSEAS ocean prediction system assimilates all quality-controlled observations, including satellite sea surface temperature and altimetry, as well as surface and profile temperature and salinity data, using the NRL-developed Navy Coupled Ocean Data Assimilation (NCODA) system [
129].
The hydrodynamic model, Hybrid Coordinate Ocean Model (HYCOM; hycom.org), uses as outer model the operational GLoBal HYCOM (GLB-HYCOM) with horizontal resolution 1/12° (approximately 9 km) and 32 vertical layers (
https://www.nrl.navy.mil/) [
130]. The GLB-HYCOM model has been used in oil spill simulations, as in the case of a Brazilian oil spill model implementation, using the Spill, Transport, and Fate Model (STFM) [
44], and in offshore India, coupled with GNOME [
70]. In the Gulf of Mexico, the Gom-HYCOM model has 1/25° horizontal resolution, vertical resolution of 20 hybrid layers, and current predictions every 3 hours [
131]. Retrospective model results are included in the reanalysis dataset of GoM-HYCOM, i.e., the HYCOM-NRL Reanalysis product (GOMu0.04/expt_50.1) produced by the US Naval Research Laboratory’s (NRL). The product has 1/25°spatial resolution (~ 3.5 km) at mid-latitudes, 36 vertical layers, and contains current predictions for the Gulf of Mexico every 3 hours. This dataset can be downloaded from these links:
http://tds.hycom.org/thredds/catalog/datasets/GOMu0.04/expt_50.1/data/netcdf/catalog.html, http://hycom.org/data/gomu0pt04/expt-50pt1. Similarly, the GoM-HYCOM includes the Real-time dataset, the HYCOM-NRL Real-time [
130] that uses the product HYCOM + NCODA GOM 1/25° with spatial horizontal resolution of 1/25° and 36 vertical layers, producing hourly 3D outputs in netCDF format (
https://www.hycom.org/data/goml0pt04/expt-31pt0). These two datasets, HYCOM-NRL Reanalysis and HYCOM-NRL Real-time, have been used as forcing inputs in several oil spill case-studies, as in French-McCay et al. [
66,
67]. On the other hand, the GLB-HYCOM 1/12° is used to provide boundary conditions to a regional implementation for the GoM, having higher horizontal resolution (1/50°) with 32 hybrid vertical layers (GoM-HYCOM 1/50°) for the Atlantic Ocean areas, over the Southeastern US Continental Shelf. The GoM-HYCOM model has been implemented in a near-real time mode, by the Coastal and Shelf Modeling Group at the Rosenstiel School of Marine and Atmospheric Science (RSMAS), University of Miami (
https://coastalmodeling.rsmas.miami.edu/), together with the Ocean Modeling and OSSE Center (OMOC) between RSMAS and the NOAA Atlantic Oceanographic and Meteorological Laboratory (AOML). The model covers the entire Gulf of Mexico, as well as a portion of the Caribbean Sea, the Florida Straits, and a portion of the Atlantic Ocean, along Florida, Georgia, and the Bahamas Islands. Le Hénaff and Kourafalou [
22] and Androulidakis et al. [
132] conducted detailed descriptions of the technical characteristics of the GoM-HYCOM 1/50° simulation (parameterizations, initial, boundary, and atmospheric forcing) and extended evaluations against non-assimilated in situ and satellite observations. The model user's manual contains additional information about the HYCOM model (
www.hycom.org). The GoM-HYCOM with 1/50° horizontal resolution has been used in oil spill simulations, like in the studies of Hole et al. [
47], Androulidakis et al. [
46], and Kourafalou et al. [
49]. Moreover, the FKeys-HYCOM model, based on HYCOM, is a high-resolution forcing model for oil spill simulations, covering the Southern Florida coastal and shelf areas and the Straits of Florida, with horizontal resolution 1/100° (~ 1 km) and 26 vertical layers. Besides, it has enabled new findings in eddy variability, with Kourafalou et al. [
133,
134] presenting more detailed information and data-based evaluation of its simulations. FKeys-HYCOM has been integrated in the oil spill simulations of Hole et al. [
48] and Androulidakis et al. [
46].
The North Carolina State University (NCSU) developed the South Atlantic Bight and Gulf of Mexico (SABGOM) hydrodynamic model, based on the Regional Ocean Modeling System (ROMS). A model implementation for the GoM exists [
135,
136] with horizontal resolution ~5 km and 36 vertical layers. French-McCay, et al. [
66] used SABGOM in their oil spill model. Moreover, SABGOM has now been replaced by the Coupled Northwest Atlantic Prediction System (CNAPS) (
http://omgsrv1.meas.ncsu.edu:8080/CNAPS/), covering a larger area than SABGOM [
137]. CNAPS is a three-dimensional marine environmental nowcast and forecast model, developed by OOMG (Ocean Observing and Modeling Group). The model computes the daily fields of ocean circulation, wave, and atmospheric variables. In addition, the SABGOM developed the Intra-Americas Sea Regional Ocean Modeling System (IAS ROMS) with a horizontal resolution of ~ 6 km and 30 vertical levels. Chao et al. [
138] developed an IAS ROMS simulation (version "4C") for year 2010, including a 2-km nested grid within the coarser and extended IAS ROMS domain, as part of the trustees' NRDA program.
SANIFS (Southern Adriatic Northern Ionian coastal Forecasting System) is an operational coastal-ocean model, developed by the CMCC-OPA (Euro-Mediterranean Centre for Climate Change), producing short-term forecasts. The operational chain is based on a downscaling approach that begins with a large-scale system for the entire Mediterranean basin (MFS, Mediterranean Forecasting system, e.g., Oddo et al. [
139]; Tonani et al. [
140]) for the derivation of the open-sea fields. SANIFS is based on the finite-element three-dimensional hydrodynamic SHYFEM model, using an unstructured grid [
141,
142]. The horizontal resolution ranges from 3 km in open-sea to 500-50 m in coastal areas. The model configuration has been outlined to provide reliable hydrodynamics and active tracer forecasts in the mesoscale-shelf-coastal waters of South Eastern Italy (Apulia, Basilicata, and Calabria regions). The model is forced in two ways: (a) at the lateral open boundaries, using a full nesting strategy, directly imposed by the MFS (temperature, salinity, sea surface height, and currents) and the OTPS (tidal forcing) fields; and (b) at the sea surface using two alternative atmospheric forcing datasets (ECMWF-12km and COSMOME-6km) through the MFS-bulk-formulae [
143,
144]. SANIFS open-sea features were validated by comparing model results to observed data, such as Argo floats, CTDs, XBTs, and satellite SSTs, as well as MFS operational products. The model's large-scale oceanographic dynamics are completely consistent with MFS [
145]. SANIFS model results have been imported as sea surface boundary condition in MEDSLIK-II model [
56].
Moreover, the NorShelf model, developed by the Norwegian Meteorological Institute, provides forecasted ocean currents for the Norwegian Shelf Sea, based on the Regional Ocean Modeling System (ROMS), with a 4D-variational (4D-Var) DA assimilation scheme (MET Norway). To accommodate the scale of the available observations and to compromise the need to resolve high resolution eddy dynamics, while confining non-linearities that limit the 4D-Var DA capabilities, a horizontal model resolution of 2.4 km was chosen. The model is intended to be used for the forecasting of ocean circulation and hydrography beyond the coastal area, including the entire shelf sea and the dynamics of the North Atlantic current at the shelf slope [
146]. Röhrs et al., [
23] used the NorShelf model outputs in OpenOil simulations.
In parallel, the IRD (French Institute of Research for the Development) provides the oceanic modeling system, CROCO (Coastal and Regional Ocean COmmunity model), which is based on the ROMS-UCLA model [
147] and ROMS AGRIF model [
148]. CROCO is a free surface hydrostatic C-grid model with a terrain-following coordinate system and an efficient split-explicit approach for distinguishing between barotropic and baroclinic terms. It is the oceanic component of a complex coupled system that includes the atmosphere, surface waves, marine sediments, biogeochemistry, and ecosystems, among others [
68,
149]. CROCO also offers MATLAB-based pre-processing tools [
150] for creating a model grid, interpolating atmospheric and oceanic data as boundary and forcing input, and setting tides and rivers in the grid model. MPI and OPENMP computations are supported by CROCO. In GNOME oil spill simulations, CROCO has been implemented with 1 km horizontal resolution and 25 vertical layers in the study of Nugroho et al. [
68].
In addition, the DHI (Danish Hydraulic Institute) developed the MIKE 21 model (
https://www.mikepoweredbydhi.com/areas-of-application), a comprehensive simulation system for hydraulics and hydraulic-related phenomena in estuaries, coastal waters, and seas [
151]. It simulates unsteady two-dimensional flows in a single layer of fluid (vertically homogeneous) [
152]. In oil spill modeling MIKE 21 has been applied by Liu et al. [
42]. Hydrodynamics, advection-dispersion, short waves, sediment transport, water quality, eutrophication, and heavy metals are all simulated by the system's various modules. The system has many engineering and environmental applications, including coastal hydraulics, oceanography, wave dynamics, harbors, rivers, environmental hydraulics, and sediment processes. A comprehensive pre- and post-processing module allows for the analysis and graphical presentation of data and model results [
153]. The versatile, interactive menu system simplifies data entry, model input, and program execution.
SHYFEM (Shallow water HYdrodynamic Finite Element Model) is 3D hydrodynamic model, developed by ISMAR-CNR (Institute of Marine Sciences - National Research Council). The model solves the system of primitive equations by vertically integrating them across each vertical layer, using the Boussinesq approximation horizontally and the hydrostatic approximation vertically. It employs the generic ocean turbulence model [
154] to calculate vertical diffusivity and viscosity. It was integrated with a transport simulation module and has a spatial resolution ranging from 25 m in extremely coastal areas or shallow waters, to a few kilometers offshore [
142,
155]. In oil spill modeling SHYFEM has been used by Cucco and Daniel [
156], integrating the Lagrangian trajectory and weathering module (FEMOIL) into the operational forecasting system (BOOM), and Ribotti et al. [
104], coupling it with the MEDSLIK-II model.
The POSEIDON System was established by the Hellenic Center for Marine Research (HCMR) and provides its hydrodynamic forecasting products [
157] also in the MEDESS4MS format to suit the needs of oil spill models for the entire Mediterranean and the Aegean Sea. The POSEIDON hydrodynamical forecasts are released through the implementation of the Princeton Ocean model (POM) [
128] with a spatial resolution of 10 km and 24 vertical sigma layers [
19,
108]. The POSEIDON Mediterranean model provides boundary conditions at the POSEIDON Aegean model’s western and eastern open boundaries [
107,
157]. Every week, the POSEIDON Aegean Sea model is re-initialized using the HCMR Mediterranean model analysis at 3.5 km horizontal resolution [
108].
The CYCOFOS (Cyprus Coastal Ocean Forecasting System) provides hydrodynamic data [
158] for the Eastern Mediterranean, covering the Aegean Sea and the Levantine basin. CYCOFOS hydrodynamic model is based on a modified version of POM [
128,
159] with a spatial resolution of 2 km and 30 vertical sigma layers, while for oil spill models produce 15 vertical z-layers in MEDESS4MS format [
108]. The CYCOFOS hydrodynamical model is nested to the Copernicus marine service, while the surface forcing is provided by the ECMWF. CYCOFOS generates daily, 6-hourly mean forecasts for the following four and a half days in a dedicated NetCDF files, designed specifically to cover the needs of oil spill models, known as MEDESS4MS format [
108]. The CYCOFOS ocean forecasts have been extensively validated and compared to the parent models, as well as satellite remote SST and in-situ observations [
60,
158]. The CYCOFOS hydrodynamical forecasts have been used in real oil spill pollution incidents [
58,
62] and in several oil spill test-cases [
59,
108,
160].
The INGV (Istituto Nazionale di Geofisica e Vulcanologia) provides hydrodynamic data products utilizing the MFS (Mediterranean Forecasting System) model [
144,
161] and the high-resolution Adriatic Forecasting System (AFS) model [
162]. The MFS is based on the NEMO Ocean General Circulation Model (OGCM), which is applied on a model domain with a spatial resolution of 6.5 km. The model spans the whole Mediterranean Sea, and provides daily 10-day forecasts. It has been implemented in many oil spill scenarios, such as the works of Coppini et al. [
163] and De Dominicis et al. [
19]. The OGCM is linked to a Wave Watch III implementation [
164] for the entire Mediterranean Sea, at the same resolution with the hydrodynamics. Furthermore, the Adriatic Forecasting System (AFS) model receives the initial and lateral boundary conditions for temperature, salinity, and velocity from the MFS to produce high resolution (~2.2 km) forecasting outputs for the Adriatic Sea with ECMWF forcing. AFS has been coupled with MEDSLIK-II in oil spill simulations via the De Dominicis et al. [
29].
The CNR, Institute for the Marine and Coastal Environment (Naples) (CNR - IAMC) computes hydrodynamics using the Western Mediterranean Model (WMED) [
155]. This forecasting system of the marine circulation at sub-regional scale (about 3.5 km) covers the Western Mediterranean, around Sardinia. The model is based on the three- dimensional primitive equation, finite difference hydrodynamic model, named POM [
128]. It has been coupled with MEDSLIK-II oil spill model [
19,
108] and MOTHY model [
156].
The IASA provides forecasting ocean products for the Eastern Mediterranean for the next 5 days via the Aegean Levantine Eddy Resolving Model (ALERMO) [
165]. This model is a high-resolution implementation of the POM used in the Aegean-Levantine basins, with horizontal resolution of 3.5 km and 25 logarithmic sigma levels in the vertical. The Copernicus Med-MFC (
https://marine.copernicus.eu/about/producers/med-mfc) is used to define the one-way nested open boundary conditions [
166]. ALERMO has been applied in oil spill modeling as forcing data by Zafirakou-Koulouris et al. [
103] and Zodiatis et al. [
108,
158]. Moreover, the IFREMER (French Research Institute for Exploitation of the Sea) provides oceanographic forecasts, produced by the PREVIMER-MENOR model, which cover the northern part of the Western Mediterranean Sea with 1.2 km horizontal resolution and 60 sigma levels, refined near the surface. The boundary conditions of PREVIMER-MENOR are provided via the Copernicus Med-MFC model. This configuration, which is based on a primitive equation model devoted to regional and coastal modeling, is utilized for both operational and academic reasons [
167,
168]. The PREVIMER-MENOR model has been coupled to oil spill simulations by De Dominicis et al. [
19] and Zodiatis et al. [
108].
Table 2.
Hydrodynamics models with the corresponding domains and horizontal resolutions.
Table 2.
Hydrodynamics models with the corresponding domains and horizontal resolutions.
Hydrodynamics |
Provider |
Geographic Coverage |
Spatial Resolution |
Reference |
Poseidon Med Model |
HCMR |
Mediterranean |
~ 10 km |
[19,108] |
Poseidon Aegean Model |
HCMR |
Aegean Sea |
~ 3.5 km |
[108] |
CYPOM |
CYCOFOS |
Aegean-Levantine |
~ 2 km |
[108,160,163] |
|
|
|
|
|
WMED |
CNR IAMC |
Western Mediterranean |
~ 3.5 km |
[19,108,156] |
ALERMO |
IASA |
Eastern Mediterranean |
~ 3.5 km |
[103,108,158] |
MFS |
INGV |
Mediterranean |
~ 6.5 km |
[19,108,163] |
AFS |
INGV |
Adriatic Sea |
~ 2.2 km |
[29,108] |
PREVIMER MENOR |
IFREMER |
North Western Mediterranean |
~ 1.2 km |
[19,108] |
ΜΙΚΕ21 |
DHI |
Regional |
- |
[42] |
CROCO |
IRD |
Regional |
1 km |
[68] |
NorShelf |
Norwegian Meteorological Institute |
Norwegian Shelf Sea |
2.4 km |
[23] |
SANIFS |
CMCC-OPA |
Mediterranean basin |
3 km |
[56] |
SHYFEM |
ISMAR-CNR |
Regional |
4 km, 1km |
[104,155] |
NEMO |
CMEMS |
Mediterranean Global Global |
(1/24°) ~ 4 km (1/12°) ~ 9 km (1/4°)~27 km |
[50,54,55] [51,53,57] [69] |
NGOM |
NOAA- CSDL |
North-eastern and Central GOM |
5-6 km |
[66] |
NCOM |
NOAA FNMOC |
American Seas and Alaska |
3 km |
[66,74] |
SABGOM |
NCSU |
GOM |
~ 5 km |
[66] |
IASROMS |
NCSU |
GOM |
~ 2 km |
[66] |
GLB-HYCOM |
NOAA NRL |
Global |
1/12° (~ 9 km) |
[44,70] |
GoM-HYCOM |
NOAA NRL |
GOM |
(~ 4 km) |
[66,67] |
GoM-HYCOM |
NOAA NRL |
GOM |
1/50° (~ 2 km) |
[46,47,49] |
Fkeys-HYCOM |
NOAA NRL |
South Florida coastal, shelf areas and Straits of Florida |
1/100° (~ 1 km) |
[46,48] |
3.3. Waves
The ECMWF provides global wave forecasts with a spatial resolution 1/8° using the third generation spectral WAve Model (WAM), [
95,
169,
170]. With 25 frequencies and 24 directions, the WAM model computes the two-dimensional wave distribution. In addition, WAM model has a 1/8° horizontal resolution with outputs every 12-hours. It is the first wave model to solve the complete action density equation, which includes non-linear wave-wave interactions. The WAM model is used operationally in global and regional applications to forecast the sea state. The model may be used for a variety of applications, including ship routing and offshore activities, as well as the validation and interpretation of satellite observations. The output of this wave data product is provided in netCDF format (
https://www.ecmwf.int/en/forecasts/datasets/set-ii).
Moreover, the CMEMS, Copernicus Marine System [
171] performs operational wave simulations and provides several wave datasets. In this review two datasets are considered: (a) The Global Reanalysis, which is a global wave reanalysis product, describing historic sea states, since 1993. This dataset is based on the WAVERYS model (WAVeReanalYSis) [
172], subject to the MFWAM (Météo-France WAve Model) model [
173], a third-generation wave model that calculates the directional wave spectrum (i.e., the distribution of sea state energy in frequency and direction) on a 1/5° irregular grid. Average wave quantities derived from this wave spectrum are the significant wave height, the average wave period and the sea surface wave Stokes drift (u and v velocities). These are provided on a regular 1/5° grid, with 3-hour time step. Besides, WAVERYS incorporates oceanic currents from the GLORYS12 physical ocean reanalysis [
174] as well as the significant wave height, observed from historical altimetry missions and directional wave spectra from Sentinel 1 SAR, since the beginning of 2017. Detailed information of this dataset is presented by Law-Chune et al. [
172]. (b) The main wave product of the Mediterranean Sea Forecasting system (MEDSEA_ANALYSISFORECAST_WAV_006_017) [
175], which is made up of hourly wave parameters, with a horizontal resolution of 1/24° that cover the Mediterranean Sea and extend up to 18.125W into the Atlantic Ocean [
176]. The wave forecast component (Med-Waves system) is based on the upgraded WAM Cycle 4.6.2. (
https://github.com/mywave/WAM). With 24 directional and 32 logarithmically distributed frequency bins, the Med-Waves modeling system resolves the prognostic part of the wave spectrum, and the model solutions are corrected by an optimal interpolation data assimilation scheme of all available along track satellite significant wave height observations. The model employs wave spectra from the GLOBAL ANALYSIS FORECAST WAV 001 027 product (
https://data.marine.copernicus.eu/product/GLOBAL_ANALYSISFORECAST_WAV_001_027/description) for Open Boundary Conditions [
177,
178]. The wave system includes two forecast cycles that provide a Mediterranean wave analysis twice per day and wave forecasts for the next ten days.
The POSEIDON System, established by Hellenic Center for Marine Research (HCMR) [
157], provides the sea state forecasting outputs for the Mediterranean and Aegean/Ionian Seas via WAM Cycle 4 model with spatial resolutions of 10 km and 3.5 km, respectively [
179]. It is a third-generation wave model capable to compute the spectra of randomly produced short-crested wind waves. The wave forecasting system generates wave forecasts for the next five days, based on hourly analysis and forecasted winds, generated by the POSEIDON weather prediction system.
The Cyprus Coastal Ocean Forecasting and Observing System products (CYCOFOS) [
180,
181] provides hourly wave forecasting data using the WAM4 model for the Mediterranean and the Black Sea, at 5 km horizontal resolution [
19,
181], in MEDESS4MS netCDF format to suit the oil spill models.
The UoM provides the wave forecasting data produced by the MALTA Maria WAM model for the central Mediterranean region, with 12.5 km horizontal resolution, using the 3rd generation spectral wave model WAM Cycle 4 [
109,
182].
The IFREMER provides wave data for the whole Mediterranean, as distributed to the MEDESS-4MS service, through the PREVIMER-MENOR-WW3 model, with 10 km spatial resolution. This wave model is based on the WaveWatch III configuration [
183]. Also, it is forced by the atmospheric model of the French Metoffice, the ARPEGE. In parallel, the Puertos del Estado (PdE) (
https://www.mitma.gob.es/empresas-fomento/puertos-del-estado ) distributes wave data through the PdE WAM wave forecasting model, with 8 km horizontal resolution, covering the western Mediterranean domain [
108].
The Météo-France supports a global forecasting wave system, named MFWAM [
184], which is based on the wave model WAM [
177]. It is a global model that has 1/12° horizontal resolution and 3-hour instantaneous temporal resolution [
83].
Table 3.
Wave models and wave products with the corresponding domains and horizontal resolutions.
Table 3.
Wave models and wave products with the corresponding domains and horizontal resolutions.
Wave system |
Provider |
Geographical area |
Spatial Resolution |
Data Type |
Reference |
Poseidon WAM Cycle 4 Med |
HCMR |
Mediterranean |
~ 10 km |
Forecast |
[108,179] |
Poseidon WAM Cycle 4 Aegean |
HCMR |
Aegean |
~ 3.5 km |
Forecast |
[108,179] |
WAM4 |
CYCOFOS |
Mediterranean and Black Sea |
~ 5 km |
Forecast |
[19,108,181] |
PdE-WAM |
PdE |
Western Mediterranean |
~ 8 km |
Forecast |
[108] |
PREVIMER-MENOR-WW3 |
IFREMER |
Mediterranean |
~ 10 km |
Forecast |
[19,108] |
MALTA/Maria WAMI |
UOM |
Central Mediterranean |
~ 12.5 km |
Forecast |
[108,109] |
WAM Cycle 6, WAM 4.6.2 |
CMEMS |
Mediterranean |
1/24° (~ 4.5 km) |
Forecast/Reanalysis |
[50,176] |
MFWAM |
Meteo-France |
Global |
1/12° (~9 km) |
Forecast |
[83,177] |
WAVERYS |
CMEMS |
Global |
1/5° (~ 22 km) |
Reanalysis |
[51,172] |
WAM |
ECMWF |
Global |
0.125° (~ 13 km) |
Forecast |
[23,46,47,48,49,185] |