At the Hydro-meteo-climate service of the Regional environmental agency of Emilia-Romagna, Italy (Arpae-SIMC), the oceanographic numerical model AdriaROMS is used in the operational forecasting suite to compute sea level, temperature, salinity and 3-D current fields of the Adriatic Sea (northern Mediterranean Sea). In order to evaluate the performance of the sea-level forecast and to study different configurations of the ROMS model, two marine storms occurred on the Emilia Romagna coast during the winter 2015–2016 are investigated.
The main focus of this study is to analyse the sensitivity of the model to the horizontal resolution and to the meteorological forcing. To this end, the model is run with two different configurations and with two horizontal grids at 1 and 2 km resolution. To study the influence of the meteorological forcing, the two storms have been reproduced by running ROMS in ensemble mode, forced by the 16-members of the meteorological ensemble COSMO-LEPS system. Possible optimizations of the model set-up are deduced by the comparison of the different run outputs.
Arpae-SIMC is the weather service of the Emilia-Romagna region and provides weather warnings for the Region and for the Italian Civil Protection Department. For these activities, Arpae-SIMC runs a chain of operational meteorological, oceanographic and wave numerical models at different horizontal resolution. Additionally, it provides tools, products and indicators for the evaluation of the impacts of storms on the coast that are discussed every day by the regional task force, together with the Civil Protection Agency and the regional Geological Seismic and Soil Survey. Weather warnings are based, for coastal events, on the combined information of sea level and significant wave height: the correct forecasting of waves and sea level, which is accomplished by the operational wave and ocean models, is therefore of crucial importance for accurate marine early warnings (Harley et al., 2016).
AdriaROMS (Chiggiato and Oddo, 2008) is the oceanographic model that runs operationally at Arpae-SIMC and it is an implementation of the ROMS model (Regional Ocean Modeling System, Shchepetkin and McWilliams, 2005; Haidvogel et al., 2008) for the Adriatic Sea, in the Mediterranean (see Fig. 1). ROMS is a free-surface, terrain-following ocean model that solves the primitive equations of the ocean: it evolves the 3-D equations for momentum and tracers together with the equations for free surface and for the vertically integrated momentum.
Computational grid and bathymetry of the Adriatic Sea
With the aim to update the oceanographic model suite, a study has been carried out to evaluate and compare different ROMS implementations. Also the Coupled-Ocean-Atmosphere-Wave-Sediment Transport Modeling System (COAWST, Warner et al., 2010), which is a model that couples an ocean-sediment model with waves and atmospheric models, is here used as a pure oceanographic model, i.e. as an alternative implementation of ROMS. This study aims to investigate the performance of the models to properly reproduce sea level in case of extreme conditions, especially during storm surge events. In this study we focus in particular on the sensitivity of the models to the horizontal resolution and to the meteorological forcing.
Two moderate storms that occurred in the winter 2015–2016 have been analysed in this study.
The first event is a one-day event, which started in the late afternoon of
the 21 November 2015 (see the event observations in Fig. 2). It is triggered
by a Bora wind event, with cold winds from NE that reached about 20 m s
21 November 2015 storm.
The second storm, of moderate intensity, occurred in February 2016 following
Sirocco winds from SE. It lasted two days, from 28 to 29 February
2016. The winds in Ravenna and in Ferrara area reached the velocity of 15 m s
Simulation runs for the two event cases with the models Ah and Aop, for the grids G1 (1 km horizontal resolution) and G2 (2 km horizontal resolution).
AdriaROMS has been running operationally at Arpae-SIMC since 2005 and it
underwent a major update in 2010 (Russo et al., 2013). AdriaROMS has a
regular grid with 2 km horizontal resolution and 20 vertical
An alternative implementation of ROMS for the Adriatic Sea has been
developed within the HAZADR project (HAZADR final publication, 2015) using
the COAWST system (similar to the work of Benetazzo et al., 2013), using the
ROMS model only. In order to distinguish between the two different
implementations of ROMS, we refer to the operational implementation of
AdriaROMS as “Aop” and to the implementation using COAWST for HAZADR as
“Ah”. The main differences of the two implementations are the following:
the the generic length scale vertical mixing implementation: the model Aop uses
a the tidal forcing, which includes diurnal and semi-diurnal tidal
constituents: five of Aop (O1, K1, N2, M2, S2) and eight of Ah (O1, K1, K2,
N2, M2, S2, P1, Q1).
To infer the influence of the horizontal resolution, the Ah setup was run on
two different computational grids: the grid G2 (which is currently used for
the operational Aop), characterized by 2 km horizontal resolution and a
minimum bathymetry of 4 m, and the grid G1, with 1 km horizontal resolution
and a minimum bathymetry of 1 m. The events have been simulated also using
the Aop model on the G2 grid. It is thus possible to compare the differences
between the effects on the grid resolution only, and the effects on the two
different model configurations separately.
Additionally, in order to evaluate the meteorological forcing effect, the
Aop and Ah models have been run in ensemble mode, forced by the Arpae
16-members COSMO-LEPS ensemble system (Montani et al., 2011). All test runs
have a forecast length of
The outputs of the model Ah runs on the two different grids G1 and G2 are
shown in Fig. 3 for the November 2015 event. All runs start at 20 November
2015, 00:00 UTC with a forecast length of
21 November 2015 storm. Measured sea level in Ancona and Vieste (grey line) and forecasted sea level with models Ah (with grids G1 and G2) and Aop with grid G2.
Comparison of the ECMWF analysis (black) and the deterministic run
of COSMO-LEPS (red) after
The outputs are also compared with sea-level data retrieved from the Italian ISPRA-RMN tide gauge network for Ancona and Vieste tide-gauge stations (see also map in Fig. 1): in this case, Ancona is strongly affected by a storm surge during the event, while Vieste is not affected by the storm and presents a dominant tidal cycle. It is clear that, for both grids G1 and G2, the surge in the north-central part of the Adriatic is not correctly reproduced by the models, which underestimates the sea level in Ancona of about 0.3 m at the surge peak. In Vieste, the tidal oscillation is in good agreement with observed data. A correct astronomical tide modelling, but a missing storm surge component in the sea level forecast is confirmed comparing simulations against the data recorded by the other available tide gauges on the Italian Adriatic side (not shown here).
In Fig. 3, the sea level forecasted by the model Ah run on the grid G2 can be compared with the sea level forecasted by the model Aop on the same grid. It is evident that the tidal trend of the model Aop is shifted respect to the Ah model and to the data, while the trend of the model Ah is correct. The comparison between the two models show that 5 astronomical components do not seems sufficient to correctly reproduce tides.
In order to verify the influence of the meteorological forcing in the sea level forecasting, the models Aop and Ah have been run in ensemble mode with respectively grid G2 and G1, forced by the output fields of the 16 members COSMO-LEPS suite. Since the atmospheric forcing plays a significant role for the sea level modelling, it is important to take into account that the meteorological forecast presents some discrepancies with the real situation in both cases. By comparing the fields of COSMO-LEPS with the ECMWF analysis (see Fig. 4) for the event of November 2015, the minimum of pressure located on the Adriatic sea is shifted towards north-east in the forecast, which leads to different wind speeds and pattern along the central-northern Adriatic coasts. The same comparison for the 2016 event shows that the forecasts of the ensemble runs start drifting from the analysis starting at 29 February, 00:00 and that the different location of the minimum of pressure on the Mediterranean near Corsica leads to a different rotation of winds in the central Mediterranean and along the Adriatic basin.
The ensemble mean on the whole Adriatic basin and the ensemble spread (defined as the standard deviation of a variable with respect to the ensemble mean) of the wind speed and mean sea level pressure is in general agreement with the spread of the COSMO-LEPS system and it is shown in Fig. 5 for the two storm events, together with the mean of the same fields of the deterministic run of COSMO-LEPS and of the ECMWF analysis.
Mean values mean sea level pressure
In both cases, the storm surge component is not reproduced correctly. The November 2015 event interested mainly the northern part of the basin, as it can be seen by comparing the tide gauge data (in Fig. 6) of Ancona and Vieste. The comparison with sea-level data at the Ancona station for the 2015 event shows that the storm surge component of sea level is highly underestimated and that the spread in the forecasted sea level, generated by the meteorological ensemble, is not sufficient to include the measurements. However, it is important to take into account that this event was not correctly reproduced by the meteorological forecast and that this discrepancy contributes to the sea level error. The maximum forecasted sea level spread is computed for selected stations for the November 2015 event in Table 2: it increases from south to north to reach about 0.1 m, but it remains quite low respect to the storm surge, that can be estimated of about 0.25–0.3 m in Ancona from tide-gauge data (Fig. 6, upper panels).
Comparison of the sea level data (grey line) recorded in Ancona
and Vieste with the ensemble output for models Aop (in blue) and Ah (in red)
for the 2015 event
The event of February 2016 generated a small storm surge in the northern Adriatic basin (see Fig. 6, lower panels). Despite the moderate intensity of the event, the ensemble predictions for this case give a higher sea level spread than for the November 2015 event: the sea level ensemble spread has higher values not only for selected locations (see Table 3), but also considering the mean on the whole basin and coasts (see Fig. 7). Also for this case, the sea level trend is not correctly reproduced by the forecast, and it seems overestimated for both Ancona and Vieste (see Fig. 6, lower panels).
Maximum sea level spread for selected tide-gauge stations for the November 2015 event.
Ensemble spread for sea level for the whole basin (thick lines) and
for coastal grid cells (thin lines) for the 2015 event
The sea level spread of the two models Ah and Aop does not present significant differences, as it can be seen also in Fig. 7, where the sea level spread for the whole Adriatic basin, and for the coastal grid cells (defined as the cells connected to a land cell) is shown for the two events. The maximum spread for the whole surface basin is about 2 cm for both events and it occurs simultaneously with the storm peaks. The spread for the coastal region is higher than the one for the whole Adriatic basin, and it is of about 3 cm for the 2015 event case and of about 4 cm for the 2016 event case. Both models have a similar variability to the meteorological ensemble forcing (Fig. 7), but a different performance, especially for the February 2016 event (Fig. 6): the different model parametrization might contribute to the missing storm surge modelling, but its influence must be carefully investigated.
Maximum sea level spread for selected tide-gauge stations for the February 2016 event.
In order to investigate the possible causes of the storm surge underestimation, a study was carried out testing two different ROMS model parametrization, Ah and Aop, to reproduce two marine storms that occurred during the winter 2015–2016 in the central-northern part of the Adriatic Sea. The tests carried out included the comparison of two runs of the model with two different horizontal resolution (1 and 2 km) and the comparison of two runs with two different model parametrization on the same grid. An additional test was carried out taking into account the uncertainty of the meteorological forcing by forcing the ocean model with meteorological ensemble fields.
The first test that compares the runs with two different grids highlight that the model horizontal resolution does not seem to be an important factor for these type of events. From the test that compares the two model parametrizations, it is possible to conclude that the configuration differences do not account for the missing storm surge, but that the tidal forcing is not correctly reproduced with only 5 tidal components.
The meteorological forcing is a critical factor and, at least for the 2015 marine event, the meteorological forecast present some discrepancies respect to the observed conditions. As direct consequence, the oceanographic models hardly reproduce the sea level and the corresponding sea level spread of the oceanic ensemble seems to be inadequate for both cases. Despite the meteorological ensemble runs show a generally adequate meteorological spread in terms of the wind and pressure when averaged over the entire basin, the difference in the position of the minimum of the mean sea level pressure between the ensemble and the analysis may be responsible for the discrepancies in the sea level prediction at the considered locations.
It is important to be aware that the meteorological forcing is not the only cause for the missing surge in the northern Adriatic sea. In order to correctly forecast storm surges, other possible sources of error that should be investigated carefully are the forcing imposed at the southern boundary, the horizontal resolution of the meteorological forcing and the model parametrization.
Wind and atmospheric data, the tide gauge data of Porto Garibaldi
and the wave data of Arpae Nausicaa wave buoy are public and are stored in the Arpae database,
accessible with the web interface
The authors declare that they have no conflict of interest.
We thank Davide Cesari and the collegues of Arpae-SIMC for the useful discussions. Edited by: M. M. Miglietta Reviewed by: A. Benetazzo and one anonymous referee