Altimetry-derived tide model for improved tide and water level forecasting along the European continental shelf

Ocean Dynamics - With the continued rise in global mean sea level, operational predictions of tidal height and total water levels have become crucial for accurate estimations and understanding of...     

Effects of self-attraction and loading at a regional scale: a test case for the Northwest European Shelf

The impact of the self-attraction and loading effect (SAL) in a regional 2D barotropic tidal model has been assessed, a term with acknowledged and well-understood importance for global models but omitted for boundary-forced, regional models, for which the implementation of SAL is non-trivial due to its non-local nature. In order to understand the impact of the lack of SAL effects in a regional scale, we have forced a regional model of the Northwest European Continental Shelf and the North Sea (continental shelf model (CSM)) with the SAL potential field derived from a global model (GTSM), in the form of a pressure field. Impacts have been studied in an uncalibrated setup and with only tidal forcing activated, in order to isolate effects. Additionally, the usually adopted simple SAL parameterization, in which the SAL contribution to the total tide is parameterized as a percentage of the barotropic pressure gradient (typically chosen 10%), is also implemented and compared to the results obtained with a full SAL computation. A significant impact on M2 representation is observed in the English Channel, Irish Sea and the west (UK East coast) and south (Belgian and Dutch Coast) of the North Sea, with an impact of up to 20 cm in vector difference terms. The impact of SAL translates into a consistent M2 amplitude and propagation speeds reduction throughout the domain. Results using the beta approximation, with an optimal domain-wide constant value of 1.5%, show a somewhat comparable impact in phase but opposite direction of the impact in amplitude, increasing amplitudes everywhere. In relative terms, both implementations lead to a reduction of the tidal representation error in comparison with the reference run without SAL, with the full SAL approach showing further impacted, improved results. Although the overprediction of tidal amplitudes and propagation speeds in the reference run might have additional sources like the lack of additional dissipative processes and non-considered bottom friction settings, results show an overall significant impact, most remarkable in tidal phases. After showing evidence of the SAL impact in regional models, the question of how to include this physical process in them in an efficient way arises, since SAL is a non-local effect and depends on the instantaneous water levels over the whole ocean, which is non-trivial to implement.     

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.     

Coupled 1D–3D hydrodynamic modelling,with application to the Pearl River Delta

Within the hydrodynamic modelling community, it is common practice to apply different modelling systems for coastal waters and river systems. Whereas for coastal waters 3D finite difference or finite element grids are commonly used, river systems are generally modelled using 1D networks. Each of these systems is tailored towards specific applications. Three-dimensional coastal water models are designed to model the horizontal and vertical variability in coastal waters and are less well suited for representing the complex geometry and cross-sectional areas of river networks. On the other hand, 1D river network models are designed to accurately represent complex river network geometries and complex structures like weirs, barrages and dams. A disadvantage, however, is that they are unable to resolve complex spatial flow variability. In real life, however, coastal oceans and rivers interact. In deltaic estuaries, both tidal intrusion of seawater into the upstream river network and river discharge into open waters play a role. This is frequently approached by modelling the systems independently, with off-line coupling of the lateral boundary forcing. This implies that the river and the coastal model run sequentially, providing lateral discharge (1D) and water level (3D) forcing to each other without the possibility of direct feedback or interaction between these processes. An additional disadvantage is that due to the time aggregation usually applied to exchanged quantities, mass conservation is difficult to ensure. In this paper, we propose an approach that couples a 3D hydrodynamic modelling system for coastal waters (Delft3D) with a 1D modelling system for river hydraulics (SOBEK) online. This implies that contrary to off-line coupling, the hydrodynamic quantities are exchanged between the 1D and 3D domains during runtime to resolve the real-time exchange and interaction between the coastal waters and river network. This allows for accurate and mass conserving modelling of complex coastal waters and river network systems, whilst the advantages of both systems are maintained and used in an optimal and computationally efficient way. The coupled 1D–3D system is used to model the flows in the Pearl River Delta (Guangdong, China), which are determined by the interaction of the upstream network of the Pearl River and the open waters of the South China Sea. The highly complex upstream river network is modelled in 1D, simulating river discharges for the dry and wet monsoon periods. The 3D coastal model simulates the flow due to the external (ocean) periodic tidal forcing, the salinity distribution for both dry and wet seasons, as well as residual water levels (sea level anomalies) originating from the South China Sea. The model is calibrated and its performance extensively assessed against field measurements, resulting in a mean root mean square (RMS) error of below 6% for water levels over the entire Pearl River Delta. The model also represents both the discharge distribution over the river network and salinity transport processes with good accuracy, resolving the discharge distribution over the main branches of the river network within 5% of reported annual mean values and RMS errors for salinity in the range of 2 ppt (dry season) to 5 ppt (wet season).