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).     

A wetting and drying algorithm with a combined pressure/free-surface formulation for non-hydrostatic models

A wetting and drying method for free-surface problems for the three-dimensional, non-hydrostatic Navier–Stokes equations is proposed. The key idea is to use a horizontally fixed mesh and to apply different boundary conditions on the free-surface in wet and dry zones. In wet areas a combined pressure/free-surface kinematic boundary condition is applied, while in dry areas a positive water level and a no-normal flow boundary condition are enforced. In addition, vertical mesh movement is performed to accurately represent the free-surface motion. Non-physical flow in the remaining thin layer in dry areas is naturally prevented if a Manning–Strickler bottom drag is used. The treatment of the wetting and drying processes applied through the boundary condition yields great flexibility to the discretisation used. Specifically, a fully unstructured mesh with any finite element choice and implicit time discretisation method can be applied. The resulting method is mass conservative, stable and accurate. It is implemented within Fluidity-ICOM [1] and verified against several idealized test cases and a laboratory experiment of the Okushiri tsunami.     

Historical analysis (2000-2005) of the coastal water quality in San Andrés Island, SeaFlower Biosphere Reserve, Caribbean Colombia

To understand the coastal water quality of San Andrès Island, and provide tools for the management of its marine resources, we present the historical analysis of the island monitoring, which includes ammonia, nitrites, nitrate, phosphates, fecal and total coliforms. The anthropogenic pressure on the coastal system is heavy, with water nutrification, posing at risk seagrass and coral ecosystems. During dry season, biologically available nitrogen is 3-9 times higher than the maximum recommended for coral reefs, while during wet season values are 2-6.4 times the maximum. Biologically available phosphorous is also high, 1-8 times the maximum during dry season, 2-13 times during wet season. In some sites the concentration of pathogenic bacteria is above the limits set by law for primary and secondary contact. It is urgent to improve the management of sewage discharge, the main polluting source of San Andres coastal waters.     

A simple approach to adjust tidal forcing in fjord models

To model currents in a fjord accurate tidal forcing is of extreme importance. Due to complex topography with narrow and shallow straits, the tides in the innermost parts of a fjord are both shifted in phase and altered in amplitude compared to the tides in the open water outside the fjord. Commonly, coastal tide information extracted from global or regional models is used on the boundary of the fjord model. Since tides vary over short distances in shallower waters close to the coast, the global and regional tidal forcings are usually too coarse to achieve sufficiently accurate tides in fjords. We present a straightforward method to remedy this problem by simply adjusting the tides to fit the observed tides at the entrance of the fjord. To evaluate the method, we present results from the Oslofjord, Norway. A model for the fjord is first run using raw tidal forcing on its open boundary. By comparing modelled and observed time series of water level at a tidal gauge station close to the open boundary of the model, a factor for the amplitude and a shift in phase are computed. The amplitude factor and the phase shift are then applied to produce adjusted tidal forcing at the open boundary. Next, we rerun the fjord model using the adjusted tidal forcing. The results from the two runs are then compared to independent observations inside the fjord in terms of amplitude and phases of the various tidal components, the total tidal water level, and the depth integrated tidal currents. The results show improvements in the modelled tides in both the outer, and more importantly, the inner parts of the fjord.     

Estimating open boundary conditions from coastal tidal observations by adjoint approach

A numerical study for estimating the tidal open boundary conditions of a shelf current modrl from tb coastal tidal observations is presented. The method is based on the optimal control/adjoint method. A lrast square fitting of the model state to simulated data is used. Two ideal domains and coastlines are considered. Using the IAP shallow. water model and its adjoint model, some identical twin experiments are carried out to test efficiency and lirnilsd of the method. The results show that the adjoint method can efficiently estimate the open boundary conditions well for gulf/bay like domains. The adjoint method seems to have great potential to improve the accuracy of tide and shelf current modeling in coastal regions. Project supported hy the National Natural Science Fuundation of China (Grant No. 49376256)     

First- and higher-order dynamical controls on water exchanges between tidal basins and the open ocean. A case study for the East Frisian Wadden Sea

Observational data, high-resolution numerical modelling results and a simple analytical theory are combined in this paper to demonstrate the dependence of the volume transports through tidal inlets on topographical or morphological parameters of a Wadden Sea system. The area of interest covers the East Frisian Wadden Sea and consists of seven weakly connected tidal basins. The observations include time series of tidal gauge data and surface currents measured at a pile station in the backbarrier basin of the island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The simple theoretical concepts presented illustrate the effect of topography (hypsometry) in the tidal basins on the temporal variability of the exchange of water. This topographic control is effectuated through the bottom slope in the areas prone to drying and flooding. For our study area it takes about twice as long from slack water to maximum flood current than from slack water to maximum ebb current. The underlying physics of this signal modulation from a more or less harmonic forcing at the open-sea boundary and the quantification of the contributing physical processes are the major results of this paper. Estimates based simply on volume conservation are consistent with observations and results from numerical modelling, but they do not completely capture the actual non-linear tidal response. Our analysis shows that at least during part of the tidal cycle characteristic topographic parameters of the inlet/bay system have a major impact on the rate of exchange of waters between the Wadden Sea and the open ocean. This impact is especially strong during the transition between flood and ebb conditions. The possible morphodynamic responses are also addressed focusing on some common (universal) topographic features in seven tidal basins.Responsible Editor: Hans Burchard     

The tidal and wind induced hydrodynamics of the composite system Adriatic Sea/Lagoon of Venice

A shallow water hydrostatic 2D hydrodynamic numerical model, based on the boundary conforming coordinate system, was used to simulate aspects of both general and small scale oceanic features occurring in the composite system constituted by the Adriatic Sea and the Lagoon of Venice (Italy), under the influence of tide and realistic atmospheric forcing. Due to a specific technique for the treatment of movable lateral boundaries, the model is able to simulate efficiently dry up and flooding processes within the lagoon. Firstly, a model calibration was performed by comparing the results of the model, forced using tides and ECMWF atmospheric pressure and wind fields, with observations collected for a set of 33 mareographic stations uniformly distributed in the Adriatic Sea and in the Lagoon of Venice. A second numerical experiment was then carried out by considering only the tidal forcing. Through a comparison between the results obtained in the two experiments it was possible to assess the reliability of the estimated parameter through the composite forcing. Model results were then verified by comparing simulated amplitude and phase of each tidal constituent as well as tidal velocities simulated at the inlets of the lagoon and in the Northern Adriatic Sea with the corresponding observed values. The model accurately reproduces the observed harmonics: mean amplitude differences rarely exceed 1 cm, while phase errors are commonly confined below 15°. Semidiurnal and diurnal currents were correctly reproduced in the northern basin and a good agreement was obtained with measurements carried out at the lagoon inlets. On this basis, the outcomes of the hydrodynamic model were analyzed in order to investigate: (i) small-scale coastal circulation features observed at the interface between the adjoining basins, which consist often of vortical dipoles connected with the tidal flow of Adriatic water entering and leaving the Lagoon of Venice and with along-shore current fields connected with specific wind patterns; (ii) residual oscillations, which are often connected to meteorological forcing over the basin. In particular, it emerges that small-scale vortical features generated near the lagoon inlet can be efficiently transported toward the open sea, thus contributing to the water exchange between the two marine regions, and a realistic representation of observed residual oscillations in the area would require a very detailed knowledge of atmospheric as well as remote oceanic forcing.     

Tidal straining effect on the suspended sediment transport in the Huanghe (Yellow River) Estuary, China

Tidal straining effect on sediment transport dynamics in the Huanghe (Yellow River) estuary was studied by field observations and numerical simulations. The measurement of salinity, suspended sediment concentration, and current velocity was conducted during a flood season in 1995 at the Huanghe river mouth with six fishing boats moored at six stations for 25-h hourly time series observations. Based on the measurements, the intra-tidal variations of sediment transport in the highly turbid river mouth was observed and the tidal straining effect occurred. Our study showed that tidal straining of longitudinal sediment concentration gradients can contribute to intra-tidal variability in sediment stratification and to asymmetries in sediment distribution within a tidal cycle. In particular, the tidal straining effect in the Huanghe River estuary strengthened the sediment-induced stratification at the flood tide, thus producing a higher bottom sediment concentration than that during the ebb. A sediment transport model that is capable of simulating sediment-induced stratification effect on the hydrodynamics in the bottom boundary layers and associated density currents was applied to an idealized estuary to demonstrate the processes and to discuss the mechanism. The model-predicted sediment processes resembled the observed characteristics in the Huanghe River estuary. We concluded that tidal straining effect is an important but poorly understood mechanism in the transport dynamics of cohesive sediments in turbid estuaries and coastal seas.     

Tidal and Fluvial Influence on Shallow Groundwater Fluctuation in Coastal Wetlands in Yellow River Delta,China

In coastal wetland, groundwater is influenced by both tidal processes and land hydrological processes. To study the influences of tidal processes and river runoff on the shallow groundwater dynamic in coastal wetland of the Yellow River Delta, surface and shallow groundwater depth were monitored. It was found that in the east part of the study area where close to the shoreline, surface is characteristic of obvious periodicity, the fluctuation of surface water was mainly controlled by tidal processes. As to the shallow groundwater, fluctuation of water depth is also influenced by the tidal process and has obvious periodicity similar to the surface water. In the north part of the study area, the shallow groundwater is mainly influenced by the Yellow River runoff. The influence distance of tidal processes and river runoff on shallow groundwater system is 14 700 and 11 600 m, respectively.     

Seasonal variation in cascade‐driven hyporheic exchange,northern Honduras

A characterization of hyporheic exchange for dry and wet season baseflow, as well as partially dewatered discharge, was done in Prieta Creek, a first‐order cascade in northern Honduras. The cascade had discharges from 1 to 15 l s^?1, had average slopes of 12%, pool spacing of 3 m, and shallow substrate of sand and gravel. Tracer tests were conducted in a 15‐m sub‐reach, a length considered to be adequate for the experiment based on the DaI test, a ratio of exchange and transport processes. In the three tests, between 9 and 18% of tracer was not recovered, possibly due to entrainment in flowpaths passing beneath the downstream monitoring location. Tracer data were analysed by the one‐dimensional transport with inflow and storage (OTIS) transient storage model (TSM) to derive standard exchange parameters, and by the solute transport in rivers (STIR) model to examine hyporheic residence time distributions (RTDs). The best fit of the observed tracer breakthrough curves was obtained by using the STIR model with a combination of two exponential RTDs to represent hyporheic retention. With increasing discharge, the OTIS model predicted increasing storage exchange fluxes and exchange coefficients and decreasing storage zone areas and transient storage times, which are trends supported by riparian and streambed piezometric head data. Riparian water levels rose during the transition from the dry to wet season, which could constrict the hyporheic storage zone. Thirteen of the 19 streambed piezometers recorded seasonal changes in hydraulic gradients and flux direction, with fewer yet stronger upwelling zones during higher discharges. The MODFLOW model missed the observed seasonal changes, possibly due to subtle changes in the seasonal change in water surface profiles. We conclude that partially dewatered dry season exchange, compared to wet season exchange, was initiated and terminated with smaller pressure gradients and, in different streambed locations, was smaller in volume, had longer residence times, and may connect with deeper and longer flow paths. Copyright © 2010 John Wiley & Sons, Ltd.     

Surface water–groundwater exchange dynamics in a tidal freshwater zone

In coastal rivers, tides can propagate for tens to hundreds of kilometres inland beyond the saltwater line. Yet the influence of tides on river–aquifer connectivity and solute transport in tidal freshwater zones (TFZs) is largely unknown. We estimate that along the TFZ of White Clay Creek (Delaware, USA), 11% of river water exchanges through tidal bank storage zones. Additional hyporheic processes such as flow through bedforms likely contribute even more exchange. The turnover length associated with tidal bank storage is 150 km, on the order of turnover lengths for all hyporheic exchange processes in non‐tidal rivers of similar size. Based on measurements at a transect of piezometers located 17 km from the coast, tides exchange 0.36 m³ of water across the banks and 0.86 m³ across the bed per unit river length. Exchange fluxes range from ?1.66 to 2.26 m day^?1 across the bank and ?0.84 to 1.88 m day^?1 across the bed. During rising tide, river water infiltrates into the riparian aquifer, and the downstream transport rate in the channel is low. During falling tide, stored groundwater is released to the river, and the downstream transport rate in the channel increases. Tidal bank storage zones may remove nutrients or other contaminants from river water and attenuate nutrient loads to coasts. Alternating expansion and contraction of aerobic zones in the riparian aquifer likely influence contaminant removal along flow paths. A clear need exists to understand contaminant removal and other ecosystem services in TFZs and adopt best management practices to promote these ecosystem services. Copyright © 2015 John Wiley & Sons, Ltd.     

Beach water table fluctuations due to spring–neap tides: moving boundary effects

Tidal water table fluctuations in a coastal aquifer are driven by tides on a moving boundary that varies with the beach slope. One-dimensional models based on the Boussinesq equation are often used to analyse tidal signals in coastal aquifers. The moving boundary condition hinders analytical solutions to even the linearised Boussinesq equation. This paper presents a new perturbation approach to the problem that maintains the simplicity of the linearised one-dimensional Boussinesq model. Our method involves transforming the Boussinesq equation to an ADE (advection–diffusion equation) with an oscillating velocity. The perturbation method is applied to the propagation of spring–neap tides (a bichromatic tidal system with the fundamental frequencies ω₁andω₂) in the aquifer. The results demonstrate analytically, for the first time, that the moving boundary induces interactions between the two primary tidal oscillations, generating a slowly damped water table fluctuation of frequency ω₁−ω₂, i.e., the spring–neap tidal water table fluctuation. The analytical predictions are found to be consistent with recently published field observations.     

Verification of mid-ocean ballast water exchange using naturally occurring coastal tracers

We examined methods for verifying whether or not ships have performed mid-ocean ballast water exchange (BWE) on four commercial vessels operating in the Pacific and Atlantic Oceans. During BWE, a ship replaces the coastal water in its ballast tanks with water drawn from the open ocean, which is considered to harbor fewer organisms capable of establishing in coastal environments. We measured concentrations of several naturally occurring chemical tracers (salinity, six trace elements, colored dissolved organic matter fluorescence and radium isotopes) along ocean transects and in ballast tanks subjected to varying degrees of BWE (0-99%). Many coastal tracers showed significant concentration changes due to BWE, and our ability to detect differences between exchanged and unexchanged ballast tanks was greatest under multivariate analysis. An expanded dataset, which includes additional geographic regions, is now needed to test the generality of our results.     

Tidal and wind-induced circulation within the Southeastern limit of the Bay of Biscay: Pasaia Bay,Basque Coast

The Basque coastal area, in the southeastern Bay of Biscay, can be characterised as being more influenced by land climate and inputs, than other typically ‘open sea’ areas. The influence of coastal processes, together with the presence of irregular and steep topography, complicate greatly the water circulation patterns. Water movement along the Basque coastal area is not well understood; observations are scarce and long-term current records are lacking. The knowledge available is confined to the surface currents: the surface water circulation is controlled mainly by wind forcing, with tidal and density currents being weak. However, there is a lack of knowledge available on currents within the lower levels of the water column; likewise, on the main time-scales involved in the water circulation. This study quantifies the contribution of the tidal and wind-induced currents, to the overall water circulation; it identifies the main time-scales involved within the tidal and wind-induced flows, investigating difference in such currents, throughout the water column, within Pasaia Bay (Basque coast). Within this context, extensive oceanographic and meteorological data have been obtained, in order to describe the circulation. The present investigation reveals that the circulation, within the surface and the sub-surface waters, is controlled mainly by wind forcing fluctuations, over a wide range of meteorological frequencies: third-diurnal, semidiurnal and diurnal land–sea breezes; synoptic variability; frequencies, near fortnightly periods; and seasonal. At the lower levels of the water column, the main contribution to the water circulation arises from residual currents, followed by wind-induced currents on synoptic time-scales. In contrast, tidal currents contribute minimally to the overall circulation throughout the water column.     

Seasonal,tidal, and geomorphic controls on sediment export to Amazon River tidal floodplains

Mainstem–floodplain material exchange in the tidal freshwater reach of major rivers may lead to significant sequestration of riverine sediment, but this zone remains understudied compared to adjacent fluvial and marine environments. This knowledge gap prompts investigation of floodplain-incising tidal channels found along the banks of tidal rivers and their role in facilitating water and suspended-sediment fluxes between mainstem and floodplain. To evaluate this role, and how it evolves along the tidal river and with time, we measured water level, flow velocity, temperature, and suspended-sediment concentration (SSC) in four tidal channels along the tidal Amazon River, Brazil. Eleven deployments were made during low, rising, high, and falling seasonal Amazon discharge. Generally, channels export high-SSC water from the mainstem to the tidal floodplain on flood tides and transfer low-SSC water back to the mainstem on ebbs. Along the length of the tidal river, the interaction between tidal and seasonal water-level variations and channel–floodplain morphology is a primary control on tidal-channel sediment dynamics. Close to the river mouth, where tides are large, this interaction produces transient flow features and current-induced sediment resuspension, but the importance of these processes decreases with distance upstream. Although the magnitude of the exchange of water and sediment between mainstem and floodplain via tidal channels is a small percentage of the total mainstem discharge in this large tidal-river system, tidal channels are important conduits for material flux between these two environments. This flux is critical to resisting floodplain submergence during times of rising sea level. © 2019 John Wiley & Sons, Ltd.     

Multi-tracer assessment of seasonal water source changes in coastal water systems along the southeastern coast of Ivory Coast (West Africa)

ABSTRACT

Proper management of coastal freshwater resources depends on an understanding of processes controlling their chemistry and seasonal flowpaths. A quantitative approach involving the coupling of major solutes and isotopes (δ¹⁸O, δ²H) of 180 samples in end-member mixing analysis (EMMA) was adopted to elucidate seasonal patterns of hydraulic exchanges amongst coastal waters along the Ebrié Lagoon catchment, Ivory Coast. The results show that the Ebrié Lagoon is a hydrologically dynamic system. In the dry season, evaporation and seawater inflow are the dominating processes, while in the wet season, river discharge is the main water source in the lagoon. Regional geology plays a significant role in aquifer recharge patterns. The Quaternary aquifer responds faster to precipitation, while the Mio-Pliocene aquifer is recharged indirectly via floodplain seepages. Salinization of over 90% of wells arises from hydrological exchanges with the Ebrié Lagoon. A diluted seawater effect was recorded in wells during the wet season owing to the relative increase in freshwater inflow.     

Large seasonal changes in the recharge of seawater in a subterranean estuary revealed by a radon tracer

Activities of radon (²²²Rn) in groundwater were continuously monitored in a saline aquifer from September 2010 to July 2011. The activities of ²²²Rn ranged from 200 to 4300 Bq m^?3, with a large seasonal variation. Because the activity of ²²²Rn in seawater is low, ²²²Rn in saline groundwater must be produced in the aquifer from radium (²²⁶Ra) in rocks and sediments. The ²²²Rn activities were higher in the warm‐dry seasons (September–November 2010 and April–May 2011) when the saline aquifer was stable. In contrast, the lowest ²²²Rn activities were observed in the cold‐dry season (December 2010 and January–March 2011), because of the effective exchange between groundwater and seawater. In addition, sudden decreases of ²²²Rn activities coincided with episodic drops in groundwater temperatures. These results reveal that lower seawater temperature in winter may result in density‐driven seawater intrusion. During the wet season (June–July 2011), ²²²Rn activities were more clearly affected by semi‐monthly and diurnal tidal pumping, showing higher ²²²Rn activities during low and spring tides. Such a tidal effect was not clearly observed during the warm‐dry and cold‐dry seasons. This result implies that the residence time of SGD in coastal zones is significantly affected by seasonal changes in driving forces such as tidal pumping and seawater intrusion. Copyright © 2016 John Wiley & Sons, Ltd.     

Tidal residual eddies and their effect on water exchange in Puget Sound

Tidal residual eddies are one of the important hydrodynamic features in tidally dominant estuaries and coastal bays, and they could have significant effects on water exchange in a tidal system. This paper presents a modeling study of tides and tidal residual eddies in Puget Sound, a tidally dominant fjord-like estuary in the Pacific Northwest coast, using a three-dimensional finite-volume coastal ocean model. Mechanisms of vorticity generation and asymmetric distribution patterns around an island/headland were analyzed using the dynamic vorticity transfer approach and numerical experiments. Model results of Puget Sound show that a number of large twin tidal residual eddies exist in the Admiralty Inlet because of the presence of major headlands in the inlet. Simulated residual vorticities near the major headlands indicate that the clockwise tidal residual eddy (negative vorticity) is generally stronger than the anticlockwise eddy (positive vorticity) because of the effect of Coriolis force. The effect of tidal residual eddies on water exchange in Puget Sound and its subbasins was evaluated by simulations of dye transport. It was found that the strong transverse variability of residual currents in the Admiralty Inlet results in a dominant seaward transport along the eastern shore and a dominant landward transport along the western shore of the inlet. A similar transport pattern in Hood Canal is caused by the presence of tidal residual eddies near the entrance of the canal. Model results show that tidal residual currents in Whidbey Basin are small in comparison to other subbasins. A large clockwise residual circulation is formed around Vashon Island near entrance of South Sound, which can potentially constrain the water exchange between the Central Basin and South Sound.