Three‐dimensional hydrodynamic and salinity transport modelling of Danshuei River estuarine system and adjacent coastal sea,Taiwan

A three‐dimensional, time‐dependent hydrodynamic and salinity model was applied to the Danshuei River estuarine system and adjacent coastal sea in Taiwan. The model forcing functions consist of tidal elevations along the open boundary and freshwater flows from the main stem and tributaries in the Danshuei River system. The bottom roughness height was calibrated and verified with model simulation of barotropic flow, and the turbulent diffusivities were calibrated through comparison of time‐series of salinity distributions. The overall model verification was achieved with comparisons of residual current and salinity distribution. The model simulation results are in qualitative agreement with the available field data. The model was then used to investigate the tidal current, residual current, and salinity patterns under the low freshwater flow condition in the modelling domain. The results reveal that the extensive intrusion of saline water imposes a significant baroclinic forcing and induces a strong residual circulation in the estuary. The downriver net velocity in the upper layer increases seaward despite the enlargement of the river cross‐section in that direction. Strong residual circulation can be found near the Kuan‐Du station. This may be the result of the deep bathymetric features there. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical determination of residence time and age in a partially mixed estuary using three-dimensional hydrodynamic model

This paper documents a numerical modeling study to calculate the residence time and age of dissolved substances in a partially mixed estuary. A three-dimensional, time-dependent hydrodynamic model was established and applied to the Danshuei River estuarine system and adjacent coastal sea in Taiwan. The model showed good agreement with observations of surface elevation, tidal currents and salinity made in 2002. The model was then applied to calculate the residence time and age distribution response to different freshwater discharges with and without density-induced circulations in the Danshuei River estuarine system. Regression analysis of model results reveals that an exponential equation can be used to correlate the residence time to change of freshwater input. The simulated results show it takes approximately 10, 4.5, and 3 days, respectively, for a water parcel that has entered the headwaters of the estuary to be transported out of the estuary under low, mean, and high flow conditions with density-induced circulation. The calculated age with density-induced circulation is less than that without density-induced circulation. The age of the surface layer is less than that at the bottom layer. Overall the study shows that freshwater discharges are the important factors in controlling the transport of dissolved substances in the Danshuei River estuarine system.     

Effects of channel regulation on salt intrusion and residual circulation of Keelung River

A vertical (laterally integrated) two‐dimensional numerical model was applied to study the hydrodynamic characteristics, salt‐water intrusion and residual circulation in the Danshuei River estuarine system. The cross‐sectional profiles measured in 2001 and 1990 respectively represent the conditions after and before channel regulation in the Keelung River. The model was re‐verified with the available hydrological data measured in 2001. Detailed model re‐verification has been conducted with water surface elevations, tidal current, and salinity distributions measured. The overall performance of the model is in qualitative agreement with the available field data. The model was then used to investigate the change in tidal ranges, salt‐water intrusion, and residual circulation as a result of channel regulation in the Keelung River. The model simulations indicate that more tidal energy propagates into the estuarine system before channel regulation because of the substantial increase in river cross‐sections. The residual circulations before channel regulation are greater than those after channel regulation and result in the limits of the salt intrusion before channel regulation being extended farther inland than those after channel regulation. This may show that channel regulation for flood control in the Keelung River did not contribute to the expansion of the mangrove areas and the disappearance of freshwater marshes at the Kuan‐Du wetlands. Copyright © 2005 John Wiley & Sons, Ltd.     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

The role of salinity in fluvio-deltaic morphodynamics: A long-term modelling study

Salinity difference between terrestrial river discharge and oceanic tidal water plays a role in modifying the local flow field and, as a consequence, estuarine morphodynamics. Although widely recognized, recent numerical studies exploring the long-term morphological evolution of river-influenced estuaries with two-dimensional, depth-averaged models have mostly neglected salinity. Using a three-dimensional morphodynamic model, we aim to gain more insight into the effect of salinity on the morphodynamics of fluvio-deltaic systems. Model results indicate that the resultant estuarine morphology established after 600 years differs remarkably when a salinity gradient is included. A fan-shaped river-mouth delta exhibits less seaward expansion and is cut through by narrower channels when salinity is included. The inclusion of salinity tends to generate estuarine circulation, which favours landward sediment transport and hence limits the growth of the delta while enhancing the development of intertidal areas. The formation of deltaic channel–shoal patterns resulting from morphodynamic evolution tends to strengthen salinity stratification, which is characterized by an increased gradient Richardson number. The direction of the depth-averaged residual sediment transport over a tide may be opposite to the direction of residual velocity, indicating the significant influence of baroclinic effects on the net sediment transport direction (and hence morphological change). The effect of salinity on morphological evolution becomes less profound when the strength of tidal or fluvial forcing is dominant over the other. The effects of sediment type and flocculation, which are particularly important when salinity gradients are present, are also discussed. Overall, this study highlights that neglecting salinity to simulate long-term estuarine morphodynamics requires more careful justification, particularly when the environment is characterized by fine sediment types (favouring suspended transport), and relatively large river discharge and estuarine depth (favouring baroclinic effects). © 2020 John Wiley & Sons, Ltd.     

A 3-D numerical study of salinity variations in the Bohai Sea during the recent years

Salinity is an important component of the marine system. Previous studies indicated that the mean salinity in the Bohai Sea had increased by 2.0 psu in the second half of the 20th century, mainly due to a sharp decrease in the Yellow River runoff, and also the effects of large-scale climatic variations and the intrusions of the North Yellow Sea Water (NYSW). Since 2002, the Yellow River Conservancy Commission has carried out the flow regulation at the beginning of every flood season, resulting in more discharge of the Yellow River freshwater into the Bohai Sea. In this study, the variations of salinity in the Bohai Sea during the recent years are investigated using a well-established three-dimensional baroclinic model, HAMburg Shelf Ocean Model (HAMSOM). The simulation results show that the Yellow River diluted water was mainly discharged into the Laizhou Bay, so the remarkable increase in the Yellow River runoff after 2002 led to a regime shift of salinity in the Laizhou Bay. However, in other parts of the Bohai Sea, salinity variation was influenced by the surrounding rivers or the intrusions of NYSW, and has little relation with the Yellow River runoff. As a whole, advection is more important than diffusion in the salinity distribution, and seasonal oscillation is the main feature of salinity variation. Via several case studies, evaporation and precipitation rates are found to be important in the long-term simulation of salinity.     

Numerical investigation of the factors influencing the vertical profiles of current,salinity, and SSC within a turbidity maximum zone

The flow-sediment interaction plays a considerable role on the vertical (internal) profiles of current,salinity and suspended sediment concentration (SSC) within a turbidity maximum zone (TMZ).Numerical modeling provides valuable insights into the complex estuarine physical processes.By combining numerical modeling with field observations,the influencing factors of fine sediment dynamics within the TMZ of Yangtze Estuary have been explored in this study.Firstly,during the neap tide,the measured data present that the current is too weak to break the density stratification,and the vertical flow structure is effectively altered.Secondly,a three-dimensional numerical model based on the Delft3D has been developed and a range of numerical sensitivity analyses were carried out to distinguish the dominant mechanisms and physical processes responsible for the phenomena observed from the measurement data.The numerical investigation highlights the following findings.(1) The vertical profile of currents within the TMZ is largely affected by saltwater intrusion,especially during lower currents when the baroclinic pressure gradient can significantly reshape the local vertical profiles of velocity.(2) The baroclinic effects are primarily determined by the stratification of salinity.(3) In addition to salinity,SSC also influences the local density stratification when its contribution to fluid density is comparable to that of salinity.(4) The settling velocity determines the overall sediment distribution and vertical profiles of the SSC in the water column.The SSC-dependent settling velocity (including the flocculation-induced acceleration and hindered settling deceleration phases) affects the longitudinal movement of the sediments.(5) The vertical profiles of current,salinity and SSC within the TMZ are highly associated with the turbulence determined by the model.The approach to modulate the vertical eddy viscosity in the model,based on the empirical dependency between Rig and Prt,may lead to a numerical instability in the stratified flow.In order to improve the stratification of SSC,additional turbulence damping effect is suggested to be implemented in the model.     

The history of water salinity in the Pearl River estuary,China, during the Late Quaternary

This research reconstructed the Late Quaternary salinity history of the Pearl River estuary, China, from diatom records of four sedimentary cores. The reconstruction was produced through the application of a diatom–salinity transfer function developed based on 77 modern surface sediment samples collected across the estuary from shallow marine environment to deltaic distributaries. The statistical analysis indicates that the majority of sediment samples from the cores has good modern analogues, thus the reconstructions are reliable. The reconstructed salinity history shows the older estuarine sequence formed during the last interglacial was deposited under similar salinity conditions to the younger estuarine sequence, which was formed during the present interglacial. Further analysis into the younger estuarine sequence reveals the interplays between sea level, monsoon‐driven freshwater discharge, and deltaic shoreline movement, key factors that have influenced water salinity in the estuary. In particular, a core from the delta plain shows the effects of sea‐level change and deltaic progradation, while cores from the mouth region of the estuary reveal changes of monsoon‐driven freshwater discharge. This study demonstrates the advantages of quantitative salinity reconstructions to improve the quality of reconstruction and allow direct comparison with other quantitative records and the instrumentally observed values of salinity. Copyright © 2010 John Wiley & Sons, Ltd.     

Controls on monthly estuarine residuals: Eulerian circulation and elevation

The Dee Estuary, at the NW English–Welsh border, is a major asset, supporting: one of the largest wildlife habitats in Europe, industrial importance along the Welsh coastline and residential and recreational usage along the English coast. Understanding of the residual elevation is important to determine the total water levels that inundate intertidal banks, especially during storms. Whereas, improved knowledge of the 3D residual circulation is important in determining particle transport pathways to manage water quality and morphological change. Using mooring data obtained in February–March 2008, a 3D modelling system has been previously validated against in situ salinity, velocity, elevation and wave observations, to investigate the barotropic–baroclinic wave interaction within this estuary under full realistic forcing. The system consists of a coupled circulation-wave-turbulence model (POLCOMS-WAM-GOTM). Using this modelling system the contribution of different processes and their interactions to the monthly residuals in both elevation and circulation is now assessed. By studying a tidally dominated estuary under wave influence, it is found that baroclinicity induced by a weak river flow has greater importance in generating a residual circulation than the waves, even at the estuary mouth. Although the monthly residual circulation is dominated by tidal and baroclinic processes, the residual estuarine surface elevation is primarily influenced by the seasonal external forcing to the region, with secondary influence from the local wind conditions. During storm conditions, 3D radiation stress becomes important for both elevation and circulation at the event scale but is found here to have little impact over monthly time scales.     

HYDRODYNAMIC MODELLING OF ESTUARINE CIRCULATION USING A SPLINE METHOD OF FRACTIONAL STEPS

I. INTRODUCTIONNumerical models of estuaries incorporating simulation of the density induced circulation are consideredimPOrtant tools in the assessment of their dispersive characteristics. Most models to date have used finite difference approximations to the governing equations (see for example thericke and Hogan, 1977, Blumberg,1 977, Perrels and Karelse, 1978). A numeriCal model of estuarine circulation using a SADI (Spline Alternating Direction implicit) procedure together with a…     

Nitrogen over enrichment in subtropical Pearl River estuarine coastal waters: Possible causes and consequences

Hong Kong is surrounded by estuarine, coastal and oceanic waters. In this study, monthly averages over a 10 year time series of salinity, temperature, chlorophyll a (chl a), dissolved oxygen (DO), dissolved inorganic nitrogen (DIN), silicate (SiO₄) and orthophosphate (PO₄) at three representative stations around Hong Kong were used to examine if excess nitrogen in estuarine influenced waters is due to P limitation. The monthly distribution clearly shows the dominant influence of the seasonal change in river discharge in the Pearl River estuary and adjacent coastal waters. In winter, the river discharge is small and more oceanic waters are dominant and as a result, salinity is high, and chlorophyll and nutrients are low. In summer, when the river discharge is high, salinity decreases and nutrients increase. DIN is very high, reaching 100 μM in the estuary. This indicates over enrichment of nitrogen relative to P and consequently there is an excess of N in coastal waters of Hong Kong. P remains low (∼1 μM) and can potentially limit both phytoplankton biomass and N utilization which was demonstrated in field incubation experiments. P limitation would result in excess N being left in the estuarine influenced waters south of Hong Kong. Phosphate concentration is lower in the Pearl River estuary than in many other eutrophied estuaries. Therefore, this relatively low PO₄ concentration should be a significant factor limiting a further increase in the magnitude of algal biomass and in the degree of eutrophication in the Pearl River estuary. The export of the excess N offshore into the northern South China Sea may result in an increase in the size of the region that is P limited in summer.     

Phytoplankton distribution and productivity in a highly turbid, tropical coastal system (Bach Dang Estuary, Vietnam)

Phytoplankton diversity, primary and bacterial production, nutrients and metallic contaminants were measured during the wet season (July) and dry season (March) in the Bach Dang Estuary, a sub-estuary of the Red River system, Northern Vietnam. Using canonical correspondence analysis we show that phytoplankton community structure is potentially influenced by both organometallic species (Hg and Sn) and inorganic metal (Hg) concentrations. During March, dissolved methylmercury and inorganic mercury were important factors for determining phytoplankton community composition at most of the stations. In contrast, during July, low salinity phytoplankton community composition was associated with particulate methylmercury concentrations, whereas phytoplankton community composition in the higher salinity stations was more related to dissolved inorganic mercury and dissolved mono and tributyltin concentrations. These results highlight the importance of taking into account factors other than light and nutrients, such as eco-toxic heavy metals, in understanding phytoplankton diversity and activity in estuarine ecosystems.     

Multi-scale modeling of Puget Sound using an unstructured-grid coastal ocean model: from tide flats to estuaries and coastal waters

Water circulation in Puget Sound, a large complex estuary system in the Pacific Northwest coastal ocean of the United States, is governed by multiple spatially and temporally varying forcings from tides, atmosphere (wind, heating/cooling, precipitation/evaporation, pressure), and river inflows. In addition, the hydrodynamic response is affected strongly by geomorphic features, such as fjord-like bathymetry and complex shoreline features, resulting in many distinguishing characteristics in its main and sub-basins. To better understand the details of circulation features in Puget Sound and to assist with proposed nearshore restoration actions for improving water quality and the ecological health of Puget Sound, a high-resolution (around 50 m in estuaries and tide flats) hydrodynamic model for the entire Puget Sound was needed. Here, a three-dimensional circulation model of Puget Sound using an unstructured-grid finite volume coastal ocean model is presented. The model was constructed with sufficient resolution in the nearshore region to address the complex coastline, multi-tidal channels, and tide flats. Model open boundaries were extended to the entrance of the Strait of Juan de Fuca and the northern end of the Strait of Georgia to account for the influences of ocean water intrusion from the Strait of Juan de Fuca and the Fraser River plume from the Strait of Georgia, respectively. Comparisons of model results, observed data, and associated error statistics for tidal elevation, velocity, temperature, and salinity indicate that the model is capable of simulating the general circulation patterns on the scale of a large estuarine system as well as detailed hydrodynamics in the nearshore tide flats. Tidal characteristics, temperature/salinity stratification, mean circulation, and river plumes in estuaries with tide flats are discussed.     

IndoTROPICS studies on the plume of the Mamberamo river into the Bismarck Sea, West Papua, Indonesia

The tropical river ocean processes in coastal settings (TROPICS) program in Indonesia (Indotropics) was carried out in the Mamberamo estuary on May–June 1999 and August 2000. The Mamberamo River flows northward from the high mountains of West Papua to the narrow and steep continental slope of the Pacific Ocean. The data for the 1999 Mamberamo estuary cruises show the dispersal of fresh water coming out of the Mamberamo River into the Bismarck Sea, and the plume is clearly defined by the pattern of salinity, turbidity, and nutrient distribution. The Mamberamo River and other nearby rivers supply high concentrations of phosphate and nitrate to the New Guinea Coastal Current, so that nutrient concentrations are higher in the surface estuarine plume, compared to offshore areas. The distribution of fishes and benthic animals were also sampled over this area.     

Modeling lateral circulation and its influence on the along-channel flow in a branched estuary

A numerical modeling study of the influence of the lateral flow on the estuarine exchange flow was conducted in the north passage of the Changjiang estuary. The lateral flows show substantial variabilities within a flood-ebb tidal cycle. The strong lateral flow occurring during flood tide is caused primarily by the unique cross-shoal flow that induces a strong northward (looking upstream) barotropic force near the surface and advects saltier water toward the northern part of the channel, resulting in a southward baroclinic force caused by the lateral density gradient. Thus, a two-layer structure of lateral flows is produced during the flood tide. The lateral flows are vigorous near the flood slack and the magnitude can exceed that of the along-channel tidal flow during that period. The strong vertical shear of the lateral flows and the salinity gradient in lateral direction generate lateral tidal straining, which are out of phase with the along-channel tidal straining. Consequently, stratification is enhanced at the early stage of the ebb tide. In contrast, strong along-channel straining is apparent during the late ebb tide. The vertical mixing disrupts the vertical density gradient, thus suppressing stratification. The impact of lateral straining on stratification during spring tide is more pronounced than that of along-channel straining during late flood and early ebb tides. The momentum balance along the estuary suggests that lateral flow can augment the residual exchange flow. The advection of lateral flows brings low-energy water from the shoal to the deep channel during the flood tide, whereas the energetic water is moved to the shoal via lateral advection during the ebb tide. The impact of lateral flow on estuarine circulation of this multiple-channel estuary is different from single-channel estuary. A model simulation by blocking the cross-shoal flow shows that the magnitudes of lateral flows and tidal straining are reduced. Moreover, the reduced lateral tidal straining results in a decrease in vertical stratification from the late flood to early ebb tides during the spring tide. By contrast, the along-channel tidal straining becomes dominant. The model results illustrate the important dynamic linkage between lateral flows and estuarine dynamics in the Changjiang estuary.     

Kinetic partitioning of Co,Mn, Cs,Fe, Ag,Zn and Cd in fresh waters (Loire) mixed with brackish waters (Loire estuary): experimental and modelling approaches

To simulate the behavior of radionuclides along a salinity gradient, in vitro sorption and desorption kinetics of Co, Mn, Cs, Fe, Ag, Zn and Cd were studied in Loire river water and the macrotidal Loire estuarine water over two different seasons. Partitioning between the dissolved phase and suspended solids were followed up over 100 h after adding radioactive tracers to freshly collected freshwater (sorption stage); this stage was followed by desorption in fresh and estuarine waters. A kinetic model describing the interactions between trace metals and particles under a salinity gradient was developed and calibrated. Among parameters and/or processes that control the fate and behavior of contaminated particles during their transfer in estuarine systems, this study shows that the speciation of trace metals is controlled by: (i) the chemical water composition: for all the elements except for Fe, desorption increased with salinity; however, the amplitude of such an effect strongly depended on the element and/or on the composition of the particulate phase (and consequently on the season); (ii) the possibility for a given element to form (or not) stable surface particle moieties such as oxides or inner-sphere complexes; (iii) the distribution of a given element among different types of sites characterised by different binding forces that can lead (or not) to re-adsorption processes after mixing of contaminated particles with uncontaminated water.Our model enabled the quantification of the contribution and the characteristic time of reactions that took place over short and long periods on the global partitioning between particulate and dissolved phases during sorption and desorption and to determine the extent to which these reactions were modified by the salinity.     

An empirical model for salinity intrusion in alluvial estuaries

The main parameters that affect the salinity intrusion in estuaries are their geometric, hydrologic and hydrodynamic characteristics. The recognition of effective parameters and understanding their roles in the salinity intrusion are required for estuarine water management. In this study, the governing equations of the salinity intrusion processes were scaled to derive the effective dimensionless parameters. Then, a previously verified model, CE-QUAL-W2, was utilized as a virtual laboratory to investigate the effects of different governing parameters on the salinity intrusion. Analysis of the results showed that logarithmic functions can be used to describe the effect of dimensionless parameters obtained by scaling of governing equations. Finally, a formula was suggested to predict the salinity intrusion length based on geometrical and hydrodynamic characteristics of alluvial estuaries.     

Assessing the geochemical reactivity of inorganic phosphorus along estuaries by means of laboratory simulation experiments

Phosphate behaviour in natural estuarine systems can be studied by performing field measurements and by undertaking laboratory simulation experiments. Thus, in this paper we describe the use of a dynamic automated estuarine simulator to characterize the geochemical reactivity of phosphate in varying salinity gradients in order to study possible mechanisms of phosphate removal from the dissolved phase (e.g. formation of some kind of apatite) and how changes in pH and salinity values influence this removal. Six laboratory assays, representing various salinity and pH gradients (average pH values between 7 and 8), were carried out. The geochemical equilibrium model MINTEQA2 was employed to characterize removal of phosphate. Among the minerals from which dissolved phosphate can originate, it seems that hydroxyapatite is by far the mineral that shows the greatest saturation indexes in the experiments. Thus, there is evidence that a type of calcium phosphate (hydroxyapatite) is involved in phosphate removal in the assays. Phosphate removal by Ca²⁺ occurs sharply at salinity values of 1–2, whereas by Fe³⁺ it is relatively gradual, at least until a salinity value of 7. Copyright © 2006 John Wiley & Sons, Ltd.