Simulation of the 1953 storm surge in the North Sea

The 1953 North Sea floods, the Big Flood, was one of the worst natural disasters in Europe in modern times and is probably one of the most studied severe coastal floods. Several factors led to the devastating storm surge along the southern North Sea coast in combination of strong and sustained northerly winds, invert barometric effect, high spring tide, and an accumulation of the large surge in the Strait of Dover. However, the storm waves and their roles during the 1953 North Sea storm surge are not well investigated. Therefore, the effect of wave setup due to breaking waves in the storm surge processes is investigated through numerical experiments. A coupled process-based tide-wave-surge model was used to investigate and simulate the storm surge in the North Sea during January 31–February 1, 1953 and validated by comparing with historical water level records at tide gauges and wave observations at light vessels in the North Sea. Meteorological forcing inputs for the period, January 27–February 3, 1953 are reproduced from ERA-20C reanalysis data with a constant correction factor for winds. From the simulation results, it is found that, in addition to the high water due to wind setup, wave setup due to breaking waves nearshore play a role of approximately 10% of the storm surge peaks with approximately 0.2 m. The resulting modeling system can be used extensively for the preparedness of the storm surge and wave of extreme condition, and usual barotropic forecast.     

Adjoint free four-dimensional variational data assimilation for a storm surge model of the German North Sea

In this paper, a data assimilation scheme based on the adjoint free Four-Dimensional Variational(4DVar) method is applied to an existing storm surge model of the German North Sea. To avoid the need of an adjoint model, an ensemble-like method to explicitly represent the linear tangent equation is adopted. Results of twin experiments have shown that the method is able to recover the contaminated low dimension model parameters to their true values. The data assimilation scheme was applied to a severe storm surge event which occurred in the North Sea in December 5, 2013. By adjusting wind drag coefficient, the predictive ability of the model increased significantly. Preliminary experiments have shown that an increase in the predictive ability is attained by narrowing the data assimilation time window.     

Climate change and North Sea storm surge extremes: an ensemble study of storm surge extremes expected in a changed climate projected by four different regional climate models

The coastal zones are facing the prospect of changing storm surge statistics due to anthropogenic climate change. In the present study, we examine these prospects for the North Sea based on numerical modelling. The main tool is the barotropic tide-surge model TRIMGEO (Tidal Residual and Intertidal Mudflat Model) to derive storm surge climate and extremes from atmospheric conditions. The analysis is carried out by using an ensemble of four 30-year atmospheric regional simulations under present-day and possible future-enhanced greenhouse gas conditions. The atmospheric regional simulations were prepared within the EU project PRUDENCE (Prediction of Regional scenarios and Uncertainties for Defining EuropeaN Climate change risks and Effects). The research strategy of PRUDENCE is to compare simulations of different regional models driven by the same global control and climate change simulations. These global conditions, representative for 1961–1990 and 2071–2100 were prepared by the Hadley Center based on the IPCC A2 SRES scenario. The results suggest that under future climatic conditions, storm surge extremes may increase along the North Sea coast towards the end of this century. Based on a comparison between the results of the different ensemble members as well as on the variability estimated from a high-resolution storm surge reconstruction of the recent decades it is found that this increase is significantly different from zero at the 95% confidence level for most of the North Sea coast. An exception represents the East coast of the UK which is not affected by this increase of storm surge extremes.     

Results of a 36-hour storm surge prediction of the North Sea for 3 January 1976 on the basis of numerical models

Summary To explore the feasibility of forescasting North Sea storm surges by integrating numerically a combined atmospheric-oceanographic physical model, the severe storm and the resulting water levels occurring on 3 January 1976 were simulated as a first step into this direction. For this purpose, the atmospheric model was run with a resolution of 8 levels in the vertical and a horizontal grid spacing of 1.4° in latitude and 2.8° in longitude on the hemisphere. The initial conditions are based upon observations of 2 January 1976, 12 GMT, i.e. about 24 hours before the storm reached its greatest intensity in the southern parts of the North Sea.The surface wind predicted by the atmospheric model was converted into stress values through a bulk formula which then entered the vertically integrated North Sea model, to yield the desired water elevations in a 22 km-grid. Moreover, also the observed wind, stemming from a careful re-analysis of the storm situation, was handled in the same way. The numerically obtained results were compared with gauge measurements at a number of coastal stations.

---

Ergebnisse einer 36stündigen Vorhersage der Nordsee-Sturmflut am 3. Januar 1976 aufgrund von numerischen Modellen

Zusammenfassung Zur Untersuchung der Möglichkeiten, Sturmfluten mit physikalischen Modellen der Atmosphäre und der Nordsee vorherzusagen, wurde als erster Schritt in dieser Richtung die schwere Sturmflut an der deutschen Nordseeküste vom 3. Januar 1976 nachgerechnet. Hierzu wurde das im Sonderforschungsbereich 94 entwickelte Atmosphärenmodell auf der Grundlage der hydro- und thermodynamischen Gleichungen mit einer Auflösung von 8 Flächen in der Vertikalen sowie 1,4° in der geographischen Breite und 2,8° in der geographischen Länge auf der gesamten Nordhalbkugel über 36 Stunden numerisch integriert. Als Eingangsdaten dienten die Beobachtungen vom 2. Januar 1976, 12 Uhr. dieser Termin liegt ca. 24 Stunden vor Eintritt des stärksten Sturmes im Bereich der Deutschen Bucht.Die auf diese Weise prognostizierte Windverteilung wurde über der Nordsee in Schubspannungswerte überführt, woraus das Nordseemodell mit den hydrodynamischen Gleichungen die Wasserstandswerte in einem 22-km-Netz berechnete. In entsprechender Weise wurde mit einem aus Beobachtungen während dieses Zeitraums gewonnenen Windfeld verfahren; diesem lag eine besonders sorgfältige, nachträglich angefertigte Analyse des Seewetteramts Hamburg zugrunde. Die so errechneten Wasserstände sind dann mit Pegelmessungen verglichen worden.

---

Résultats d'une prédiction d'onde de tempête de 36 heures de la mer du Nord pour le 3 janvier 1976 sur la base de modèles numériques

Résumé Afin d'étudier la possibilité de prévior les ondes de tempête en mer du Nord en intégrant numériquement un modèle physique combiné atmosphère-océan, la forte tempête du 3 janvier 1976 et des niveaux d'eau qui en ont résulté ont été simulés comme premier pas dans cette voie. Dans ce but, on a fait fonctionner le modèle atmosphérique avec une résolution de 8 niveaux dans le plan vertical et un réseau horizontal aux pas de 1,4° en latitude et de 2,8° en longitude sur l'hémisphère. Les conditions initiales sont basées sur les observations du 2 janvier 1976 à 12h GMT, c'est-à-dire environ 24 heures avant que la tempête n'atteignît sa plus forte intensité dans les parties méridionales de la mer du Nord.Le vent en surface prédit par le modèle atmosphérique était converti par l'intermédiaire d'une formule globale en valeurs de tensions qui étaient alors introduites dans le modèle de la mer du Nord intégré verticalement, pour obtenir les hauteurs d'eau désirées, dans un réseau de 22 km de côté. En outre, le vent observé, issu d'une nouvelle analyse soignée de la situation de tempête, était aussi traité de la même manière. Les résultats obtenus numériquement furent comparés aux mesures fournies par les marégraphes en un certain nombre de stations côtières.

     

Source term balance in a severe storm in the Southern North Sea

This paper presents the results of a wave hindcast of a severe storm in the Southern North Sea to verify recently developed deep and shallow water source terms. The work was carried out in the framework of the ONR funded NOPP project (Tolman et al. 2013) in which deep and shallow water source terms were developed for use in third-generation wave prediction models. These deep water source terms for whitecapping, wind input and nonlinear interactions were developed, implemented and tested primarily in the WAVEWATCH III model, whereas shallow water source terms for depth-limited wave breaking and triad interactions were developed, implemented and tested primarily in the SWAN wave model. So far, the new deep-water source terms for whitecapping were not fully tested in shallow environments. Similarly, the shallow water source terms were not yet tested in large inter-mediate depth areas like the North Sea. As a first step in assessing the performance of these newly developed source terms, the source term balance and the effect of different physical settings on the prediction of wave heights and wave periods in the relatively shallow North Sea was analysed. The December 2013 storm was hindcast with a SWAN model implementation for the North Sea. Spectral wave boundary conditions were obtained from an Atlantic Ocean WAVEWATCH III model implementation and the model was driven by hourly CFSR wind fields. In the southern part of the North Sea, current and water level effects were included. The hindcast was performed with five different settings for whitecapping, viz. three Komen type whitecapping formulations, the saturation-based whitecapping by Van der Westhuysen et al. (2007) and the recently developed ST6 whitecapping as described by Zieger et al. (2015). Results of the wave hindcast were compared with buoy measurements at location K13 collected by the Dutch Ministry of Transport and Public Works. An analysis was made of the source term balance at three locations, the deep water location North Cormorant, the inter-mediate depth location K13 and at location Wielingen, a shallow water location close to the Dutch coast. The results indicate that at deep water the source terms for wind input, whitecapping and nonlinear four-wave interactions are of the same magnitude. At the inter-mediate depth location K13, bottom friction plays a significant role, whereas at the shallow water location Wielingen also depth-limited wave breaking becomes important.     

Surge modelling in the eastern Irish Sea: present and future storm impact

It is believed that, in the future, the intensity and frequency of extreme coastal flooding events may increase as a result of climate change. The Natural Environment Research Council (NERC) Flood Risk from Extreme Events (FREE) project, Coastal Flooding by Extreme Events and EU FP7 Morphological Impacts and Coastal Risks Induced by Extreme Storm Events project are investigating the flood risks in the eastern Irish Sea, an area that includes most of England’s coastal types. Using a previously modelled and validated historical extreme surge event, in November 1977, we now investigate the changes in peak surge as a result of possible future climate conditions. In order to simulate the surge, we have set up a one-way nested approach, using the Proudman Oceanographic Laboratory Coastal Ocean Modelling System 3D baroclinic model, from a domain covering the whole NW European continental shelf, through to a 1.85 km Irish Sea model; both areas are forced by tides, atmospheric pressure and winds. We use this modelling system to investigate the impact of enhanced wind velocities and increased sea levels on the peak surge elevation and residual current pattern. The results show that sea level rise has greater potential to increase surge levels than increased wind speeds.     

Improved water-level forecasting for the Northwest European Shelf and North Sea through direct modelling of tide, surge and non-linear interaction

In real-time operational coastal forecasting systems for the northwest European shelf, the representation accuracy of tide–surge models commonly suffers from insufficiently accurate tidal representation, especially in shallow near-shore areas with complex bathymetry and geometry. Therefore, in conventional operational systems, the surge component from numerical model simulations is used, while the harmonically predicted tide, accurately known from harmonic analysis of tide gauge measurements, is added to forecast the full water-level signal at tide gauge locations. Although there are errors associated with this so-called astronomical correction (e.g. because of the assumption of linearity of tide and surge), for current operational models, astronomical correction has nevertheless been shown to increase the representation accuracy of the full water-level signal. The simulated modulation of the surge through non-linear tide–surge interaction is affected by the poor representation of the tide signal in the tide–surge model, which astronomical correction does not improve. Furthermore, astronomical correction can only be applied to locations where the astronomic tide is known through a harmonic analysis of in situ measurements at tide gauge stations. This provides a strong motivation to improve both tide and surge representation of numerical models used in forecasting. In the present paper, we propose a new generation tide–surge model for the northwest European Shelf (DCSMv6). This is the first application on this scale in which the tidal representation is such that astronomical correction no longer improves the accuracy of the total water-level representation and where, consequently, the straightforward direct model forecasting of total water levels is better. The methodology applied to improve both tide and surge representation of the model is discussed, with emphasis on the use of satellite altimeter data and data assimilation techniques for reducing parameter uncertainty. Historic DCSMv6 model simulations are compared against shelf wide observations for a full calendar year. For a selection of stations, these results are compared to those with astronomical correction, which confirms that the tide representation in coastal regions has sufficient accuracy, and that forecasting total water levels directly yields superior results.     

Application of an unstructured mesh model to storm surge propagation in the Mersey estuary region of the Irish Sea

The storm surge period of 13–16 November 1977 when there was a major positive surge followed by a negative surge in the Irish Sea is investigated using a two-dimensional unstructured mesh model of the west coast of Britain. The model accounts for tidal and external surge forcing across its open boundaries which are situated in the Celtic Sea and off the west coast of Scotland. Although this period has been examined previously using a uniform finite-difference model, and a finite element model, neither of these could resolve the Mersey estuary which is the focus of the present study. By using a finite element model with very high mesh resolution within the Mersey, the spatial variability of surge elevations and currents within the Mersey to rapidly changing surge dynamics can be examined. The mesh in the model varies from about 7 km in deep water, to the order of 100 m in the Mersey, with the largest mesh length reaching 17 km in deep offshore regions, and smallest of order 26 m occurring in shallow coastal regions of the Mersey estuary. The model accounts for wetting/drying which occurs in shallow water coastal areas. Calculations showed that during the positive surge period, the amplitude and speed of propagation of the surge was largest in the deep water channels. This gave rise to significant spatial variability of surge elevations and currents within the estuary. As wind stresses decreased over the Irish Sea, a negative surge occurred over Liverpool Bay and at the entrance to the Mersey. However, within the Mersey there was a local positive surge which continued to propagate down the estuary. This clearly showed that although the large scale response of the Irish Sea to changing wind fields occurred rapidly, the response in the Mersey was much slower. These calculations with a west coast variable mesh model that included a high-resolution representation of the Mersey revealed for the first time how elevations and currents within the Mersey responded to Irish Sea surges that rapidly changed from positive to negative.     

Coupled wave and surge modelling for the eastern Irish Sea and implications for model wind-stress

We revisit the surge of November 1977, a storm event which caused damage on the Sefton coast in NW England. A hindcast has been made with a coupled surge-tide-wave model, to investigate whether a wave-dependent surface drag is necessary for accurate surge prediction, and also if this can be represented by an optimised Charnock parameter. The Proudman Oceanographic Laboratory Coastal Modelling System-Wave Model (POLCOMS-WAM) has been used to model combined tides, surges, waves and wave-current interaction in the Irish Sea on a 1.85 km grid. This period has been previously thoroughly studied, e.g. Jones and Davies [Jones, J.E., Davies, A.M., 1998. Storm surge computations for the Irish Sea using a three-dimensional numerical model including wave-current interaction. Continental Shelf Research 18(2), 201–251] and we build upon this previous work to validate the POLCOMS-WAM model to test the accuracy of surge elevation predictions in the study area. A one-way nested approach has been set up from larger scale models to the Irish Sea model. It was demonstrated that (as expected) swell from the North Atlantic does not have a significant impact in the eastern Irish Sea. To capture the external surge generated outside of the Irish Sea a (1/9° by 1/6°) model extending beyond the continental shelf edge was run using the POLCOMS model for tide and surge.     

Identification of extreme storm surges with high-impact potential along the German North Sea coastline

Planning and design of coastal protection rely on information about the probabilities of very severe storm tides and the possible changes that may occur in the course of climate change. So far, this information is mostly provided in the form of high percentiles obtained from frequency distributions or return values. More detailed information and assessments of events that may cause extreme damages or have extreme consequences at the coast are so far still unavailable. We describe and compare two different approaches that may be used to identify highly unlikely but still physically possible and plausible events from model simulations. Firstly, in the case when consistent wind and tide-surge data are available, different metrics such as the height of the storm surge can be derived directly from the simulated water levels. Secondly, in cases where only atmospheric data are available, the so called effective wind may be used. The latter is the projection of the horizontal wind vector on that direction which is most effective in producing surges at the coast. Comparison of events identified by both methods show that they can identify extreme events but that knowledge of the effective wind alone does not provide sufficient information to identify the highest storm surges. Tracks of the low-pressure systems over the North Sea need to be investigated to find those cases, where the duration of the high wind is too short to induce extreme storm tides. On the other hand, factors such as external surges or variability in mean sea level may enhance surge heights and are not accounted for in estimates based on effective winds only. Results from the analysis of an extended data set suggest that unprecedented storm surges at the German North Sea coast are possible even without taking effects from rising mean sea level into account. The work presented is part of the ongoing project “Extreme North Sea Storm Surges and Their Consequences” (EXTREMENESS) and represents the first step towards an impact assessment for very severe storm surges which will serve as a basis for development of adaptation options and evaluation criteria.     

Some environmental impacts of a storm surge barrier

In 1976 the Dutch government decided to close the Eastern Scheldt estuary with a storm surge barrier instead of a dam as was previously envisaged. In this way safety is bought and the excellent Eastern Scheldt ecological system is preserved. The cost of this compromise is estimated at 8.7 billion Dutch guilders. However what does this compromise mean for the Eastern Scheldt ecosystem?     

The statistics of wave height and crest elevation during the December 2012 storm in the North Sea

This paper concerns new field measurements of wave height and crest elevation probability distributions as measured in the North Sea during a storm in December 2012. The water surface elevation was recorded by Saab WaveRadar REX instruments mounted on eight fixed-jacket platforms in addition to a Datawell Directional Waverider buoy. The storm generated an easterly sea state which peaked well in excess of the 100-year wave height for that direction in the region. Furthermore, 19 freak waves occurred during the storm according to the definition as reported by Haver (2000). The present study demonstrates that the significant steepness and spectral bandwidth during the storm remain almost constant. Consequently, there is little change in the commonly applied design wave height and crest elevation probability distributions throughout the storm. Whilst the bulk of the recorded data was in good agreement with the theoretical distributions, it was demonstrated that when the wind speed was larger than 25 m/s, the measured crest elevation lies above the second-order Forristall distribution.     

Discontinuous Galerkin methods for modeling Hurricane storm surge

Storm surge due to hurricanes and tropical storms can result in significant loss of life, property damage, and long-term damage to coastal ecosystems and landscapes. Computer modeling of storm surge can be used for two primary purposes: forecasting of surge as storms approach land for emergency planning and evacuation of coastal populations, and hindcasting of storms for determining risk, development of mitigation strategies, coastal restoration and sustainability.Storm surge is modeled using the shallow water equations, coupled with wind forcing and in some events, models of wave energy. In this paper, we will describe a depth-averaged (2D) model of circulation in spherical coordinates. Tides, riverine forcing, atmospheric pressure, bottom friction, the Coriolis effect and wind stress are all important for characterizing the inundation due to surge. The problem is inherently multi-scale, both in space and time. To model these problems accurately requires significant investments in acquiring high-fidelity input (bathymetry, bottom friction characteristics, land cover data, river flow rates, levees, raised roads and railways, etc.), accurate discretization of the computational domain using unstructured finite element meshes, and numerical methods capable of capturing highly advective flows, wetting and drying, and multi-scale features of the solution.The discontinuous Galerkin (DG) method appears to allow for many of the features necessary to accurately capture storm surge physics. The DG method was developed for modeling shocks and advection-dominated flows on unstructured finite element meshes. It easily allows for adaptivity in both mesh (h) and polynomial order (p) for capturing multi-scale spatial events. Mass conservative wetting and drying algorithms can be formulated within the DG method.In this paper, we will describe the application of the DG method to hurricane storm surge. We discuss the general formulation, and new features which have been added to the model to better capture surge in complex coastal environments. These features include modifications to the method to handle spherical coordinates and maintain still flows, improvements in the stability post-processing (i.e. slope-limiting), and the modeling of internal barriers for capturing overtopping of levees and other structures. We will focus on applications of the model to recent events in the Gulf of Mexico, including Hurricane Ike.     

Spectral analysis of storm surge in Hong Kong Victoria Harbour

Based on a linear model the dynamic characteristics of Victoria Harbour (Hong Kong) is obtained by means of spectral analysis of the storm surge hydrographs. The results show that the harbour is an ideal one which has a small gain factor and a flat response in the frequency range from 0 to 6 × 10⁻⁵ Hz. The results also show that the power spectra possess the narrow band features which indicates that the periodic components associated with tidal motions are predominant over the random components. The power spectrum corresponding to a frequency of 2.3 × 10⁻⁵ Hz is likely to be associated with the astronomical tides. The peaks in the power spectra at zero frequency suggest that the pumping mode of oscillations is dominant in a storm surge. This mode of oscillations represents the temporal variations in mean sea level. To demonstrate the full potential of the present model, more case studies should be conducted when surge as well as non-surge data are available.