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.     

Feedback between residual circulations and sediment distribution in highly turbid estuaries: An analytical model

Motivated by field studies of the Ems estuary which show longitudinal gradients in bottom sediment concentration as high as O(0.01 kg/m⁴), we develop an analytical model for estuarine residual circulation based on currents from salinity gradients, turbidity gradients, and freshwater discharge. Salinity is assumed to be vertically well mixed, while the vertical concentration profile is assumed to result from a balance between a constant settling velocity and turbulent diffusive flux. Width and depth of the model estuary are held constant. Model results show that turbidity gradients enhance tidally averaged circulation upstream of the estuarine turbidity maximum (ETM), but significantly reduce residual circulation downstream, where salinity and turbidity gradients oppose each other. We apply the condition of morphodynamic equilibrium (vanishing sediment transport) and develop an analytical solution for the position of the turbidity maximum and the distribution of suspended sediment concentration (SSC) along a longitudinal axis. A sensitivity study shows great variability in the longitudinal distribution of suspended sediment with the applied salinity gradient and six model parameters: settling velocity, vertical mixing, horizontal dispersion, total sediment supply, fresh water flow, and water depth. Increasing depth and settling velocity move the ETM upstream, while increasing freshwater discharge and vertical mixing move the ETM downstream. Moreover, the longitudinal distribution of SSC is inherently asymmetric around the ETM, and depends on spatial variations in the residual current structure and the vertical profile of SSC.     

Effect of bottom stress formulation on modelled flow and turbidity maxima in cross-sections of tide-dominated estuaries

A three-dimensional numerical model with a prognostic salinity field is used to investigate the effect of a partial slip bottom boundary condition on lateral flow and sediment distribution in a transect of a tidally dominated channel. The transect has a symmetrical Gaussian cross-channel bottom profile. For a deep, well-mixed, tidally dominated channel, partial slip decreases the relative importance of Coriolis deflection on the generation of cross-channel flow patterns. This has profound implications for the lateral distribution of residual salinity that drives the cross-channel residual circulation pattern. Transverse sediment transport, however, is always found to be governed by a balance between advection of residual sediment concentration by residual lateral flow on the one hand and cross-channel diffusion on the other hand. Hence, the changes in the cross-channel distribution of residual salinity modify the lateral sediment distribution. For no slip, a single turbidity maximum occurs. In contrast, partial slip gives a gradual transition to a symmetrical density distribution with a turbidity maximum near each bank. For a more shallow, partially mixed tidal channel that represents the James River, a single turbidity maximum at the left bank is found irrespective of the near-bed slip condition. In this case, semi-diurnal contributions to sediment distribution and lateral flow play an important role in cross-channel sediment transport. As vertical viscosity and diffusivity are increased, a second maximum at the right bank again exists for partial slip.     

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.     

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.     

Transverse structure of tidal and residual flow and sediment concentration in estuaries

An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s^− 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.     

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.     

Suspended sediment fluxes at an intertidal flat: The shifting influence of wave,wind, tidal,and freshwater forcing

Using in situ, continuous, high frequency (8–16 Hz) measurements of velocity, suspended sediment concentration (SSC), and salinity, we investigate the factors affecting near-bed sediment flux during and after a meteorological event (cold front) on an intertidal flat in central San Francisco Bay. Hydrodynamic forcing occurs over many frequency bands including wind wave, ocean swell, seiching (500–1000 s), tidal, and infra-tidal frequencies, and varies greatly over the time scale of hours and days. Sediment fluxes occur primarily due to variations in flow and SSC at three different scales: residual (tidally averaged), tidal, and seiching. During the meteorological event, sediment fluxes are dominated by increases in tidally averaged SSC and flow. Runoff and wind-induced circulation contribute to an order of magnitude increase in tidally averaged offshore flow, while waves and seiching motions from wind forcing cause an order of magnitude increase in tidally averaged SSC. Sediment fluxes during calm periods are dominated by asymmetries in SSC over a tidal cycle. Freshwater forcing produces sharp salinity fronts which trap sediment and sweep by the sensors over short (∼30 min) time scales, and occur primarily during the flood. The resulting flood dominance in SSC is magnified or reversed by variations in wind forcing between the flood and ebb. Long-term records show that more than half of wind events (sustained speeds of greater than 5 m/s) occur for 3 h or less, suggesting that asymmetric wind forcing over a tidal cycle commonly occurs. Seiching associated with wind and its variation produces onshore sediment transport. Overall, the changing hydrodynamic and meteorological forcing influence sediment flux at both short (minutes) and long (days) time scales.     

Spatial variations of flow structure over estuarine hollows

Observations of the flow field over an elongated hollow (bathymetric depression) in the lower Chesapeake Bay showed tidally asymmetric distributions. Current speed increased over the landward side of the hole during flood tides and decreased in the deepest part of the hollow during ebb tides. A simple conceptual analysis indicated that the presence of a horizontal density gradient can generate the asymmetric spatial variations of flow structure depending on the sign of the horizontal density gradient. When water density decreases downstream, the velocity increases over the downstream edge of the hollow. Conversely when water density increases downstream, the flow decreases over the hollow more than a case without a horizontal density gradient. The conceptual analysis is confirmed by numerical experiments of simplified hollows in steady open channel flows and of an idealized tidal estuary. These hollows also alter the local current field of tidally averaged estuarine exchange flows. The residual depth-averaged currents over a hollow show a two-cell circulation when Coriolis forcing is neglected and an asymmetric two-cell circulation, with a stronger cyclonic eddy, when Coriolis forcing is included.     

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.     

MODELING OF THE HIGH CONCENTRATION LAYER OF COHESIVE SEDIMENT UNDER THE ACTION OF WAVES AND CURRENTS

MODELING OF THE HIGH CONCENTRATION LAYEROF COIIESIVE SEDMNT UNDER Tus ACTIONOF WAVES AND CURRENTSQinghe ZHANG', Yongsheng WU', Jiian LIAN1 and Pingxing DING3Abstract:High concenhation layer of cohesive sediment frequenhy occurs in muddy estUaries and coastalzones, and causes raPid siltation of the waterways. A one dimensional vertical coupled modeldescribing the interactions betWeen waves, currentS and suspended cohesive sediment is develoPed inthe pre…     

The influence of asymmetries in overlying stratification on near-bed turbulence and sediment suspension in a partially mixed estuary

Data collected from the York River estuary demonstrate the importance of asymmetries in stratification to the suspension and transport of fine sediment. Observations collected during two 24-h deployments reveal greater concentrations of total suspended solids during the flood phase of the tide despite nearly symmetric near-bed tidal current magnitude. In both cases, tidally averaged net up-estuary sediment transport near the bed was clearly observed despite the fact that tidally averaged residual near-bed currents were near zero. Tidal straining of the along-channel salinity gradient resulted in a stronger pycnocline lower in the water column during the ebb phase of the tide and appeared to limit sediment suspension. Indirect measurements suggest that the lower, more intense, pycnocline on the ebb acted as a barrier, limiting turbulent length scales and reducing eddy diffusivity well below the pycnocline, even though the lower water column was locally well mixed. In order to more conclusively link changes in stratification to properties of near-bed eddy viscosity and diffusivity, longer duration tripod and mooring data from an additional experiment are examined, that included direct measurement of turbulent velocities. These additional data demonstrate how slight increases in stratification can limit vertical mixing near the bed and impact the structure of the eddy viscosity below the pycnocline. We present evidence that the overlying pycnocline can remotely constrain the vertical turbulent length scale of the underlying flow, limiting sediment resuspension. As a result, the relatively small changes in stratification caused by tidal straining of the pycnocline allow sediment to be resuspended higher in the water column during the flood phase of the tide, resulting in preferential up-estuary transport of sediment.Responsible Editor: Iris Grabemann     

Lateral circulation driven by boundary mixing and the associated transport of sediments in idealized partially mixed estuaries

A three-dimensional, hydrostatic, primitive equation numerical model with modern turbulence closures is used to explore lateral circulation and the associated transport of sediments in idealized, moderately to highly stratified estuaries. The model results suggest that boundary mixing on a sloping bottom can drive a significant amount of lateral circulation. This mechanism has received little attention to date in the estuarine literature. Good agreement with an analytical solution and similar vertical structures of lateral flows to observations from the Hudson River estuary support the importance of the boundary mixing mechanism. Boundary mixing is at least as important as differential advection for the modeled scenarios, when the two mechanisms are evaluated using the salt balance equation for model runs without rotation. Linearly superposing analytical solutions for lagged boundary mixing lateral flow and Ekman-forced lateral flow yields a good representation of the near-bottom lateral flow from the model with rotation. The 2 h lag required for the boundary mixing solution is roughly equal to the vertical diffusion time scale, indicating that lateral flow adjustment depends on development of a bottom mixed layer. Sediment dynamics at cross sections seaward and landward of the salt intrusion are very different. Seaward of the salt intrusion, sediments are eroded in the channel and preferentially deposited on the right slope (looking seaward), mainly due to the combination of high sediment concentration in the channel during flood with strong up-slope transport on that side (tidal pumping). Lateral sediment re-distribution landward of the salt intrusion is negligible due to weak residual lateral circulation.     

Sediment dynamics in an offshore tidal channel in the southern Yellow Sea

The geomorphology of the southern Yellow Sea（SYS） is characterized by offshore radial sand ridges（RSR）.An offshore tidal channel（KSY Channel） is located perpendicular to the coast,comprised of a main and a tributary channel separated by a submarine sand ridge（KSY Sand Ridge） extending seaward.In order to investigate the interactions among water flow,sediment transport,and topography,current velocity and suspended sediment concentration（SSC） were observed at 11 anchor stations along KSY Channel in RSR during a spring tide cycle.High resolution bottom topography was also surveyed.Residual currents and tidally averaged suspended sediment fluxes were calculated and analyzed by using the decomposition method.Results suggested that the water currents became stronger landward but with asymmetrical current speed and temporal duration of flood and ebb tides.Residual currents showed landward water transport in the nearshore channel and a clockwise circulation around the KSY Sand Ridge.Tidally-averaged SSC also increased landward along the channel.The main mechanisms controlling SSC variations were resuspension and horizontal advection,with spatial and temporal variations in the channel,which also contributed to sediment redistribution between channels and sand ridges.Residual flow transport and the tidal pumping effect dominated the suspended sediment flux in the KSY Channel.The KSY Sand Ridge had a potential southward migration due to the interaction between water flow,sediment transport,and topography.     

A coupled analytical model for salt intrusion and tides in convergent estuaries

ABSTRACT

In this paper we develop a coupled analytical model for salinity and tidal propagation in estuaries where the cross-sectional area varies exponentially. A simple analytical model for tidal dynamics has been used to estimate the tidal excursion, which has an important influence on the salt intrusion process since it determines the extreme salinities (i.e. salinity distribution for high water slack and low water slack). The objective of the coupling is to reduce the number of calibration parameters, which subsequently strengthens the reliability of the salt intrusion model. Moreover, the coupling enables us to assess the potential impacts of external changes, both human-induced interventions (e.g. dredging) and natural changes (e.g. global sea level rise), on the salt intrusion process. In addition, the fully analytical solution for hydrodynamics allows immediate estimation of the tidally averaged depth and friction coefficient for given water level recordings and salinity measurements. This is particularly useful when a geometric survey is not available. The coupled model has been applied to six previously unsurveyed estuaries in Malaysia and the results show that the correspondence between analytical estimations and observations is very good. Thus, the coupled model proves to be a useful tool to obtain estimates of salt intrusion in estuaries based on a minimum amount of information required and for assessing the effect of human-induced or natural changes.

EDITOR D. Koutsoyiannis ASSOCIATE EDITOR B. Dewals     

Flood–ebb and spring–neap variations of mixing,stratification and circulation in Chesapeake Bay

A three-dimensional hydrodynamic model is used to investigate intra-tidal and spring–neap variations of turbulent mixing, stratification and residual circulation in the Chesapeake Bay estuary. Vertical profiles of salinity, velocity and eddy diffusivity show a marked asymmetry between the flood and ebb tides. Tidal mixing in the bottom boundary layer is stronger and penetrates higher on flood than on ebb. This flood–ebb asymmetry results in a north–south asymmetry in turbulent mixing because tidal currents vary out of phase between the lower and upper regions of Chesapeake Bay. The asymmetric tidal mixing causes significant variation of salinity distribution over the flood–ebb tidal cycle but insignificant changes in the residual circulation. Due to the modulation of tidal currents over the spring–neap cycle, turbulent mixing and vertical stratification show large fortnightly and monthly fluctuations. The stratification is not a linear function of the tidal-current amplitude. Strong stratification is only established during those neap tides when low turbulence intensity persists for several days. Residual circulation also shows large variations over the spring–neap cycle. The tidally averaged residual currents are about 50% stronger during the neap tides than during the spring tides.     

Effects of bottom slope, flocculation and hindered settling on the coupled dynamics of currents and suspended sediment in highly turbid estuaries, a simple model

This study aims at gaining basic understanding about two specific phenomena that are observed in the highly turbid estuaries tidal Ouse, Yangtze and Ems, i.e. (1) the accumulation of suspended matter in the deeper parts of the estuaries and (2) the relatively high values of turbidity near the surface in the area of the turbidity maximum. A semi-analytical model is analysed to verify the hypothesis that these phenomena result from bottom slope-induced turbidity currents and from hindered settling, respectively. The model governs the dynamics of residual flow, driven by fresh water discharge, salinity gradients and turbidity gradients. It further uses the condition of morphodynamic equilibrium (no divergence of net sediment transport) to compute the residual sediment concentration. New aspects are that depth variations on flow and mixing processes, as well as flocculation and hindered settling of sediment, are explicitly accounted for. Tides act as a source of mixing and erosion of sediment only, thus processes like tidal pumping are not considered. Model results show that the estuarine turbidity maximum (ETM) shifts in the down-slope direction, compared to the case of a constant depth. Slope-induced turbidity currents, which are directed down-slope near the bottom and up-slope near the surface, are responsible for this shift, thereby confirming the first part of the hypothesis above. The down-slope shift of the ETM is reduced by currents resulting from gradients in depth-dependent mixing, which counteract turbidity currents, but which are always weaker. Including flocculation and hindered settling yields increased surface sediment concentrations in the area of the turbidity maximum, compared to the situation of a constant settling velocity, thereby supporting the second part of the hypothesis. Sensitivity experiments reveal that the conclusions are not sensitive to the values of the model parameters.     

Interactions among waves,current, and mud: Numerical and laboratory studies

Interactions between waves, current, mud and turbulence are very complicated in the coastal and estuarine turbid waters. It is still necessary to improve our understanding of the fundamental physical processes governing the cohesive sediment transport in the coastal and estuarine waters. A numerical model is developed to study the interactions among waves, current, and mud. An eddy viscosity model for wave and current is proposed in order to close the equations of wave motion or of current motion in a combined flow, respectively. The equations of mud transport are derived based on the visco-elastic properties of mud. Coupling the equations of wave motion or of current motion for water layer with those of mud layer can give (1) wave height; (2) distributions of current velocities in the water layer; (3) distributions of transport velocities at the water–mud interface; and (4) distributions of mass transport velocities within the mud layer. These modeled results are in a reasonable agreement with experimental results. Results suggest that (1) the rate of wave attenuation increases in the opposing currents (currents against in the direction in which the waves propagate) and decreases in the following currents (currents in the same direction as that in which the waves propagate); (2) the opposing currents would have more significant effects on the rate of wave height attenuation than the following currents; (3) the effect of current on the rate of wave attenuation on the muddy bottom is larger than that on the rigid bottom; (4) mud transport rate increased in the following currents but decreased in the opposing currents; and (5) the rate of wave height attenuation on the mud bottom is one order of magnitude larger than that on the rigid bottom.     

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.     

THREE DIMENSIONAL NUMERICAL MODELLING OF FLOW AND SEDIMENT TRANSPORT IN RIVERS

The 3D numerical model, ECOMSED (open source code), was used to simulate flow and sediment transport in rivers. The model has a long history of successful applications to oceanic, coastal and estuarine waters. Improvements in the advection scheme, treatment of river roughness parameterization and shear stress partitioning were necessary to reproduce realistic and comparable results in a river application. To account for the dynamics of the mobile bed boundary, a model for the bed load transport was included in the code. The model reproduced observed secondary currents, bed shear stress distribution and erosion-deposition patterns on a curved channel. The model also successfully predicted the general flow patterns and sediment transport characteristics of a 1-km long reach of the River Klar?lven, located in the north of the county of V?rmland, Sweden.