This study addresses the impact of coupling between wave and circulation models on the quality of coastal ocean predicting systems. This is exemplified for the German Bight and its coastal area known as the Wadden Sea. The latter is the area between the barrier islands and the coast. This topic reflects the increased interest in operational oceanography to reduce prediction errors of state estimates at coastal scales, which in many cases are due to unresolved non-linear feedback between strong currents and wind waves. In this study we present analysis of wave and hydrographic observations, as well as results of numerical simulations. A nested-grid modelling system is used to produce reliable nowcasts and short-term forecasts of ocean state variables, including waves and hydrodynamics. The database includes ADCP observations and continuous measurements from data stations. The individual and combined effects of wind, waves and tidal forcing are quantified. The performance of the forecast system is illustrated for the cases of several extreme events. The combined role of wave effects on coastal circulation and sea level are investigated by considering the wave-dependent stress and wave breaking parameterization. Also the response, which the circulation exerts on the waves, is tested for the coastal areas. The improved skill of the coupled forecasts compared to the non-coupled ones, in particular during extreme events, justifies the further enhancements of coastal operational systems by including wave effects in circulation models.
In the last decade, the north European coasts were affected by severe storms which caused serious damages in the North Sea coastal zones. Additionally, different human activities, e.g. offshore wind power industry, oil industry and coastal recreation necessitate information about the sea state in the coastal ocean with high resolution in space and time. There seems to be a consensus that high-quality predictions of extreme events, like storm surges and flooding caused by storms, could substantially contribute to avoid or minimise human and material damages and losses. Therefore, reliable wave forecasts and long-term statistics of extreme wave conditions are of utmost importance for the coastal areas. In many coastal areas the need for reliable risk assessments increases the demand of precise coastal predictions. This cannot be achieved by further neglecting the wave–current interaction in coastal ocean operational forecasting.
Nested grid model domains for the North Sea (top pattern), German Bight (middle pattern) and the east Wadden Sea (bottom pattern). The spatial resolution is 3 nm, 1 km and 200 m, respectively. The geographical location of stations and sections analysed later are shown as well.
Wave–current interaction is recently an important issue in the field of coastal ocean forecasting (Roland and Ardhuin, 2014; Bolaños et al., 2014). Understanding this process is of utmost importance on the road of fully integrating the atmospheric, wave and ocean models and their further coupling with biological, morphological and hydrographical forecasting systems. The uncertainties in most of the presently used models results from the non-linear feedback between the currents, water level variations and wind waves, which can no longer be ignored, in particular in the coastal zone. The joint impact of surges, currents and waves is strongly inter-related (Wolf et al., 2011; Brown et al., 2011) and those cannot be considered separately for coastal ocean predictions.
The ocean waves affect not only the sea level but also the currents and mixing, the latter being of utmost importance for the sediment dynamics (Lettmann et al., 2009). Prandle et al. (2000) demonstrated the need of accounting for surface waves with a significant wave height larger than one metre in the sediment modelling. This is of big importance for sediment dynamics and other ecosystem processes (Wolf and Prandle, 1999). These authors also showed the effects of waves due to surges and tides; on the other side the waves' characteristics are affected by the changes of sea level height due to tides and wind.
The main effects of waves that are commonly considered in the coupled modelling are due to radiation stress and Stokes drift. Babanin et al. (2010) showed that interaction of turbulence and bottom stress is also very important.
Wave–current interaction has been a topic of many studies recently (Ardhuin et al., 2008; Mellor, 2003, 2008, 2011; Kumar et al., 2012; Michaud et al., 2012; Zodiatis et al., 2015). Mellor (2003, 2005, 2008) extended the radiation stress formulation based on the linear wave theory of Longuet-Higgins and Stewart (1964). Bennis and Ardhuin (2011) questioned the method of Mellor and suggested the use of Lagrangian mean framework leading to the so-called vortex force. Vortex force method has been implemented in ROMS-SWAN (Kumar et al., 2012; Lane et al., 2007; McWilliams et al., 2004; Uchiyama et al., 2010). Moghimi et al. (2013) critically compared the two approaches claiming that the radiation stress formulation showed unrealistic offshore directed transport in the wave shoaling regions; on the other hand the results of longshore circulations performed similarly for both methods. Aiki and Greatbatch (2013, 2014) proved that the radiation stress formulation of Mellor is applicable for small bottom slopes. Bolaños et al. (2011, 2014) demonstrated the importance of wave–current interactions in a tidally dominated estuary and showed that the inclusion of wave effects through 3-D radiation stress improves the velocity in the study area. They also compared the different radiation stress methods and concluded that for the tidally dominated area the 3-D version of radiation stress produces better results than the 2-D version. Polton et al. (2005) found that accounting for the Coriolis–Stokes forcing results in encouraging agreement between model and measurements of the mixed layer. Janssen (2012) showed positive impact of wave breaking to the daily cycle of sea surface temperature. Later Breivik et al. (2015) demonstrated reduced bias between modelled and measured water temperature by incorporating the Coriolis–Stokes forcing, turbulence induced by breaking waves and ocean side stress in the NEMO model at global ocean scale. Weber et al. (2006) estimated that the wave induced stress is about 50 % of the total atmospheric stress for moderate to strong wind. Wolff et al. (2011) studied the effects of waves on hydrodynamics; Brown et al. (2013) considered the wave effects on the storm surges; and Roland et al. (2009) studied wave effects on water level for the Adriatic Sea. The importance of ocean depth and velocity variations for the simulated waves in the estuaries is analysed by Pleskachevsky et al. (2011) and Lin and Pierre (2003). However, within the framework of practical coastal ocean forecasting, the interactions between waves and currents are still not yet enough considered.
In this study we will address the coupling between wave and circulation models for coastal ocean prediction systems with the example of the German Bight. We do not plan to analyse the role of different parameterisations used. Rather we will demonstrate the areas of improvements of coastal ocean predictions due to coupling between wave and hydrodynamic models.
The structure of the paper is as follows. The wave and hydrodynamic models and the processes of their interaction are described in Sect. 2. Section 3 addresses the effects of hydrodynamics on wave model performance, while in Sect. 4 we discuss the effects of waves on hydrodynamics and improvement of short-term forecast; followed finally by concluding remarks.
Statistics of the validation. In addition to mean and standard deviation, the coefficients of a linear regression are given.
The General Estuarine Transport Model (GETM, Burchard and Bolding, 2002) was
used in this study to simulate the circulation. This model solves the
primitive equations for momentum, temperature, salinity and water level.
The model setup described here uses the
Elevation amplitude (cm) RMSE and mean errors (model–observations) for the coupled wave–circulation model and GETM-only model for the tide gauge data from British Oceanographic Data Centre (BODC) over the German Bight area.
Time series at the buoy Elbe station (see Fig. 1 for its location) from 1 to 10 July 2011 of
Wave model (WAM) is a third-generation wave model which solves the wave transport
equation explicitly without any presumptions on the shape of the wave
spectrum. The basic physics and numerics of the WAM Cycle 4 wave model,
which is described in Komen et al. (1994) and Guenther et al. (1992) are
kept in the new release WAM 4.5.3. In the coupled model system , the source
function integration scheme of Hersbach and Janssen (1999) and the
reformulated wave model dissipation source function (Bidlot et al., 2005),
later reviewed by Bidlot et al. (2007) and Janssen (2008), are incorporated.
Additionally, depth-induced wave breaking (Battjes and Janssen, 1978) has
been included as source function. Depth and/or current fields can be
non-stationary. It is crucial for strongly tidally forced shallow areas,
like the German Bight, to allow model grid points to fall dry and
account for refraction due to spatially varying current and depth.
These modifications are of utmost importance for the improvement of wave
modelling results in the coastal areas such as the Wadden Sea. The wave
model code is freely available under
Comparison of measured (top) and computed values of the spectral energy density at the buoy Elbe (see Fig. 1 for its location).
Similar to the circulation model, the open boundary conditions for the
German Bight WAM are taken from the regional WAM setup for the North Sea
area (with a spatial resolution of ca. 5 km). The German Bight wave model
has the same horizontal resolution and uses the same topography as the
circulation model GETM. The driving wind fields are the same as the ones
used in the hydrodynamical model. The required boundary information at the
open boundaries of the North Sea model is derived from the regional wave
model EWAM for Europe that is running twice a day in the operational wave
forecast routine of the DWD. Within the framework of Coastal Observing
System for Northern and Arctic seas (COSYNA), a pre-operational wave and
hydrodynamic forecast system has successfully been implemented and is
running continuously since December 2009, providing hindcast and forecast
data freely available on the COSYNA website under
The original version of GETM was modified to account for the depth-dependent radiation stress and Stokes drift. The terms were calculated from the integrated wave parameters according to Mellor (2008, 2011) and Kumar et al. (2012). The gradients of the radiation stresses serve as additional explicit wave forcing in the momentum equations for the horizontal velocity components. Here the Stokes drift components have been subtracted from the wave processes in order to transfer them to the Eulerian framework. Moghimi et al. (2013) studied the effects of the two different approaches utilising the radiation stress (Mellor, 2011) and vortex force (Ardhuin et al., 2008) using WAM–GETM coupled models and showed that the results for the longshore-directed transport are similar for both formulations. Recently Aiki and Greatbatch (2013) showed that the radiation stress parameterization is applicable for small bottom slopes and Grashorn et al. (2015) showed that radiation stress formalism is applicable for shallow areas like the German coastal ones. They also demonstrated that the criterion proposed by Mellor (2013) to test the applicability of the radiation stress method gives reasonable results in this region. This gives us a confidence that, despite the known limitations of the radiation stress formalism it is well applicable for our study area. Additionally, the bottom friction modifications dependent upon bottom roughness and wave properties (Styles and Glenn, 2000) have been implemented. Turbulent kinetic energy due to wave friction (wave breaking/white capping and bottom dissipation) that is wave-enhanced turbulence has also been taken into consideration (Pleskachevsky et al., 2011).
In order to demonstrate the impact of wave–current interaction on coastal model simulations we performed two different experiments. In the first one the wave model WAM and the circulation model GETM have been run separately (we will further refer to it as the non-coupled run). The results have been compared with the WAM–GETM coupled model system, in which all wave-hydrodynamic processes described above are considered. We will further refer to it as the coupled model run. Three case studies have been analysed here, which we consider interesting in terms of both atmospheric conditions/extreme events and observational data availability.
The first analyses period is in July 2011, which was a calm weather period. Two different wind regimes were dominating the atmospheric state in July 2011, which will be addressed separately.
The next two analyses periods are chosen in order to address the effects of two of the most severe storm surges affecting our study region in the last hundred years. The first storm surge is the storm Britta of 31 October–1 November 2006 causing serious damages for the off-shore infrastructures and shipping in the North Sea region. Storm Britta was characterized by a deep low-pressure centre that moved on a trajectory from north of Scotland to western Norway and then eastwards through the Baltic Sea. Severe storm surge damages occurred in the east Wadden Sea. Extreme sea level during this storm surge is considered as a 100-year event (Madsen et al., 2007). In addition to the storm surge, unusually high waves have been measured in the southern North Sea developing on northern North Sea and propagated southward under the influence of strong north winds with a long fetch. The storm Britta has been given particular attention in our analyses for the types of changes that may occur during single event (Bartholomä et al., 2009; Lettmann et al., 2009; Stanev et al., 2011; Grashorn et al., 2015).
The second extreme event that we consider here is the winter storm Xavier on 5–6 December 2013, causing severe flooding and devastation along the German North Sea coast. Besides extreme high water levels along the coasts, extreme sea state conditions have been observed causing serious erosion of dunes and sand displacements on the barrier islands.
At the buoy Elbe, which is located in the open sea (water depth about 21 m, see the middle panel of Fig. 1), two different wind regimes occurred
between 1 and 10 July 2011 (Fig. 2). From 1 to 5 July the
dominating north-western wind did not change its direction (see the red line
in Fig. 2b). However wind speed increased from 7.7 m s
Impact of hydrodynamics on waves: normalized standard deviation (estimated as the difference between the control run and the coupled run relative to the control run values) of significant wave height (Hs, left) and mean period (tm1, right) between coupled wave–circulation model and wave-only model. Averaging is for 1 month (July 2011). The 5 and 10 % isolines are plotted with white lines.
Time series of sea level elevation (SLE) in (m) at four coastal stations of the German Bight (ST1–ST4, see Fig. 1 for the locations). Black line corresponds to tide gauge observations, red line corresponds to the coupled wave–circulation model (WAM–GETM) and green line corresponds to circulation-only model (GETM).
Sea level elevation (SLE) difference (cm) between the coupled wave–circulation model (WAM–GETM) and circulation-only model (GETM) for the German Bight on 3 December 2013 at 01:00 UTC (left) and during the storm Xavier on 6 December 2013 at 01:00 UTC.
The frequency wave spectra from the Elbe buoy and the two runs are shown in Fig. 3 for the first 5 days in July during the strong wind event. Similarly to Fig. 2, the patterns of wave spectra from the measurements and those of the coupled model run are in very good agreement (compare the top and bottom panel of Fig. 3). This is not the case for the non-coupled wave model (the middle panel in Fig. 3). The tidal currents are mainly affecting the tail of the spectra, whereas the energy around the peak is not much different in all three panels. The statistical analysis of the observations and simulations (see Table 1) clearly demonstrates the improvement of the quality of coupled wave–circulation model forecasts for the German Bight in comparison to the non-coupled model.
To quantify the impact of currents, including water depth hydrodynamics on the results of wave model, the SD of Hs and the mean period (tm1), of the coupled run normalized by the mean values of the non-coupled wave model are shown in Fig. 4. The horizontal patterns are given as a 1-month average for July 2011. In the open North Sea area there are no significant differences between the coupled and non-coupled wave modes for both Hs and tm1. However, along the coastal areas, where currents and water level change rapidly under the influence of tides, the impact of coupling seems to be significant. Within the German Bight coastal areas the SD of Hs goes up to 30 %, mainly due to the changes in water depth. The SD of tm1 is about 10–15 % in the coastal area. In particular, in the south-east of the German Bight, where the rivers Elbe and Weser are entering, the impact of coupling on the tm1 period spreads much further off-shore.
Interesting to notice are several relatively small areas, mainly located on
the tidal inlets where the SD of tm1 reaches values of up to 30 %. These
areas are characterized by strong currents, up to 1.5 m s
Top: meridional (left) and zonal (right) velocity time series
(m s
Zonal velocity vertical section (m s
In this section we demonstrate the role of coupling by analysing the impact of waves on hydrodynamics during several extreme events. Sea level variability in four locations (ST1–ST4, see Fig. 1 for their geographical locations) are analysed along the German coast for the period including the extreme event Xavier on 6 December 2013 (see description in Sect. 2). The observations and simulations are shown in Fig. 5 for the tide gauge observations (black line), coupled wave–circulation model simulations (coupled run, red line) and the non-coupled run (circulation-only model, blue line). During normal meteorological conditions, the coupled and non-coupled models fit well with the tide gauge data. However, during the storm Xavier, the sea level predicted by the hydrodynamic-only model is underestimated with more than 40 cm. It appears that the sea level predictions of the coupled model are closer to the measurements (compare the red and black lines). This demonstrates the importance of wave–current interactions also for the hydrodynamics. The root mean square errors (RMSEs) between observations and coupled model have been significantly reduced compared with the RMSE differences between the observations and circulation-only model for all coastal locations (Table 2). Predictions of storm events with coupled models could be of utmost importance for many coastal applications dealing with risk analyses (off-shore wind industry, oil platform operations, etc.) where higher accuracy is needed. This justifies the consideration of waves in operational forecasting.
In order to give an idea of the spatial distribution of the effects resulting from coupling, we show in Figure 6 the differences of sea surface elevation between the coupled and circulation-only models for 3 December 2013 at 01:00 UTC (normal meteorological situation, left panel) and 6 December 2013 at 01:00 UTC (extreme event, right panel). The wave-induced parameterization increases the average water level, which is more pronounced in the coastal area. In the open North Sea the effects of coupling are almost negligible. During normal conditions the difference of the sea level due to the coupling of circulation and wave models reaches a maximum of 10–15 cm in the area of the Elbe estuary. However, during the storm Xavier, the differences of simulated sea level when considering waves are more than 30 cm along the whole German coast. In some of the Wadden Sea areas the increase of water level in the simulations taking into consideration the wave–current interactions was above a half metre. The results shown here are indicative that the uncertainties in most of the presently used non-coupled operational models result from the missing non-linear feedback between strong tidal currents and wind waves. This can no longer be ignored in the operational oceanography, in particular in the coastal zone where the wave–circulation interplay seems to be dominant. The statistical analyses of simulated SLE vs. tide gauge data over the German Bight (Table 2) show that the coupling significantly improves the ocean predictions for the whole German coastal area. The RMSEs during the calm conditions are small in both coupled and circulation-only models. However during the extreme events the RMSE of sea surface elevation are significantly reduced when considering ocean–wave interactions.
In the following we will demonstrate the effect of coupling on the storm
Britta on 1 November 2011. During this storm event (see Fig. 7a),
significant wave height over 10 m has been simulated in the open North Sea
(close to the north-western boundary). The southern North Sea area was
exposed to waves with a magnitude of about 6–7 m. Only 2 days later
significant wave height dropped to 4 m within the German Bight (Fig. 7b). As
an example of the impact of wave effects we show the dissipation of surface
turbulent kinetic energy in the German Bight area at the peak of the storm
at 03:00 UTC on 1 November (Fig. 7c) and under calm meteorological
conditions (Fig. 7d). Along the coast dissipation rates exceed 0.06 m
Predictions of both zonal and meridional velocity have been also improved
due to the coupling between the waves and circulation during storm Britta
(see Fig. 8). The zonal velocity has been under-estimated in the circulation-only
model (green line) and got closer to the ADCP data provided by the BSH for the coupled
wave–circulation model (red line). There is also a very good correlation
between the differences of the predicted velocity and significant wave
height (Fig. 8, bottom patterns). During the storm Britta when the
significant wave height reached almost 8 m in the coastal station the
difference of the zonal velocity between the coupled run and the
hydrodynamic model was more than 40 cm s
The vertical section of the intensification of the longshore currents during the storm Britta is shown on Fig. 9 (the location of the section is plotted in Fig. 1). Not only does the longshore velocity increase but also its vertical structure has been changed through the effects of coupling. Similar behaviour has been also observed by Grashorn et al. (2015).
Wave and hydrodynamic hindcast and forecast for the North Sea and German Bight are of great importance for the management of coastal zones, ship navigation, off-shore wind energy, naval operations, etc. The storms and waves which they generate have direct impact on the coastal and marine environment. The population living in the coastal areas is recently concerned with the impacts of erosion and flooding, and actions aiming at better predictions, impact assessments of minimisation of damages are of greatest importance. Some driving forces that cause serious damages in coastal environments are due to the wave conditions. Their absolute and relative impact can be estimated by using coastal models. In this paper we demonstrated the improvements of coastal ocean predictions due to consideration of wave–current interaction for the North Sea and German Bight regions.
Improved forecast statistics once considering coupling is being demonstrated for both wave and circulation models. The coupled system presented here enables to provide reliable predictions as well as analyse long-term changes of wave and circulation conditions, including extreme events. The performance of the forecasting system was illustrated for the cases of several extreme events along with the effects of ocean waves on coastal circulation. It can be concluded that for our study area the use of radiation stress parameterization produced physically reasonable results. However, the different wave-induced formalisms lead to different limitations and thus no general recommendation should be performed. The improved skill resulting from the recent coupled model developments, in particular during storms, justifies further enhancements of the both forecast applications at operational services and long-term hindcasts and climate analyses for the North Sea and the German Bight.
This work was supported by EU FP7 Project MyOcean 2, Grant agreement 283367 and Horizon2020 Project MyOcean FO, Grant agreement 633085. The authors are thankful to W. Koch for preparing the model forcing and B. Gardeike for assistance with the graphics. We acknowledge the use of the measurement data provided by BSH. The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association. Edited by: H. Wehde