Influence of variable salinity conditions in a tidal creek on riparian groundwater flow and salinity dynamics

While recent studies have revealed that tidal fluctuations in an estuary significantly affect groundwater flows and salt transport in the riparian zone, only seawater salinity in the estuary has been considered. A numerical study is conducted to investigate the influence of estuarine salinity variations on the groundwater flow and salt dynamics in the adjacent aquifer to extend our understanding of these complex and dynamic systems. Tidal salinity fluctuations (synchronous with estuary stage) were found to alter the magnitude and distribution of groundwater discharge to the estuary, which subsequently impacted on groundwater salinity patterns and residence times, especially in the riparian zone. The effects of salinity fluctuations were not fully captured by adopting a constant, time-averaged estuarine salinity. The modelling analysis also included an assessment of the impact of a seasonal freshwater flush in the estuary, similar to that expected in tropical climates (e.g. mean estuary level during flood significantly greater than average), on adjacent groundwater flow and salinity conditions. The three-month freshwater flushing event temporarily disrupted the salt distribution and re-circulation patterns predicted to occur under conditions of constant salinity and tidal water level fluctuations in the estuary. The results indicate that the salinity variations in tidal estuaries impact significantly on estuary–aquifer interaction and need to be accounted for to properly assess salinity and flow dynamics and groundwater residence times of riparian zones.     

Modelling estuarine and coastal flows using an unstructured triangular finite volume algorithm

Details are given of the development and application of a numerical model for predicting free-surface flows in estuarine and coastal basins using the finite volume method. Both second- and third-order accurate and oscillation free explicit numerical schemes have been used to solve the shallow water equations. The model deploys an unstructured triangular mesh and incorporates two types of mesh layouts, namely the ‘cell centred’ and ‘mesh vertex’ layouts, and provides a powerful mesh generator in which a user can adjust the mesh-size distribution interactively to create a desirable mesh. The quality of mesh has been shown to have a major impact on the overall performance of the numerical model.The model has been applied to simulate two-dimensional dam break flows for which transient water level distributions measured within a laboratory flume were available. In total 12 model runs were undertaken to test the model for various flow conditions. These conditions include: (1) different bed slopes (ranging from zero to 0.8%), (2) different upstream and downstream water level conditions, and (3) initially wet and dry bed conditions, downstream of the dam. Detailed comparisons have been made between model predicted and measured water levels and good agreement achieved between both sets of results. The model was then used to predict water level and velocity distributions in a real estuary, i.e. the Ribble Estuary, where the bed level varies rapidly at certain locations. In order to model the whole estuary, a 1-D numerical model has also been used to model the upper part of the estuary and this model was linked dynamically to the 2-D model. Findings from this application are given in detail.     

Analytical and numerical analysis of tides and salinities in estuaries; part II: salinity distributions in prismatic and convergent tidal channels

Estuaries, commonly, are densely populated areas serving the needs of the inhabitants in multiple ways. Often the interests are conflicting and decisions need to be made by the local managers. Intake of fresh water for consumption, agricultural purposes or use by industries may take place within a region not far landward of the limit of salt intrusion. Human interventions (e.g. deepening of the navigation channels) or climate changes (sea level rise, reduction of the river discharge) can bring these intake locations within the reach of saline or brackish water and consequently endanger their function. To support policy and managerial decisions, a profound knowledge of processes associated with the salinity structure in estuaries is required. Although nowadays advanced numerical three-dimensional models are available that are able to cope with the complexity of the physics there is still a need for relatively simple tools for quick-scan actions in a pre-phase of a project or for instructive purposes. The analytical model described in this paper may serve these needs. It computes the maximum salinity distribution using the dispersion coefficient in the mouth as the only model parameter. The model has been calibrated using observational data in a large number of estuaries and experimental data in a tidal flume. The dispersion coefficient was successfully related to geometric and hydrodynamic parameters resulting in an expression that can be used for convergent estuaries as well as prismatic channels, see Eqs. 25a and 25b. Application of the model in a predictive mode showed its promising capabilities. Comparison with three-dimensional numerical models indicates that the channel geometry in the estuary mouth largely influences dispersive processes. The analytical model for salt intrusion may be used in combination with the analytical model for tidal propagation in convergent estuaries and tidal channels by Van Rijn (part I). In this way, input is obtained on the tidal velocity amplitude and the Chézy roughness following calibration of this model on tidal amplitudes along the estuary.     

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.     

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.     

A two-phase numerical model for suspended-sediment transport in estuaries

This paper presents a full 2-D X/Z numerical model for sediment transport in open channels and estuaries using a two-phase (fluid–solid particle) approach. The physical concept and the mathematical background of the model are given and test-cases have been carried out to validate the proposed model. In order to illustrate its feasibility for a real estuary, the model has been applied to simulate the suspended-sediment transport and the formation of turbidity maximum in the Seine estuary. The numerical results show that the main characteristics of estuarine hydro-sediment dynamics in the Seine estuary are in fact reproduced by the proposed model. A qualitative agreement between the numerical results and the actual observations has been obtained and is presented in this paper.     

DEPTH-AVERAGED 2-D CALCULATION OF TIDAL FLOW,SALINITY AND COHESIVE SEDIMENT TRANSPORT IN ESTUARIES

A depth-averaged 2-D numerical model for unsteady flow, salinity and cohesive sediment transport in estuaries is established using the finite volume method on the non-staggered, curvilinear grid. The convection terms are discretized by upwind schemes, the diffusion terms are by the central difference scheme, and the time derivative terms are by the three-time-level implicit scheme. The coupling of flow velocity and water level in the 2-D shallow water equations is achieved by the SIMPLEC algorithm with the Rhie and Chow's momentum interpolation method. The sediment model calculates the settling, deposition, erosion and transport of cohesive sediment, taking into account the influence of sediment size, sediment concentration, salinity and turbulence intensity on the flocculation of cohesive sediment. The flow model is first tested against the measurement data in the Tokyo Bay and San Francisco Bay, showing good agreements. And then, the entire model of flow, salinity and sediment transport is verified in the Gironde Estuary. The water elevation, flow velocity, salinity and sediment concentration are well predicted.     

A model study of tide- and wind-induced mixing in the Columbia River Estuary and plume

A numerical simulation of circulation in the Columbia River estuary and plume during the summer of 2004 is used to explore the mixing involved as river water is transformed into shelf water. The model is forced with realistic river flow, tides, wind stress, surface heat flux, and ocean boundary conditions. Simulated currents and water properties on the shelf near the mouth are compared with records from three moorings (all in 72 m of water) and five CTD sections. The model is found to have reasonable skill; statistically significant correlations between observed and modeled surface currents, temperature, and salinity are all 0.42–0.72 for the mooring records. Equations for the tidally averaged, volume-integrated mechanical energy budget (kinetic and potential) are derived, with attention to the effects of: (i) Reynolds averaging, (ii) a time varying volume due to the free surface, and (iii) dissipation very close to the bottom. It is found that convergence of tidal pressure work is the most important forcing term in the estuary. In the far field plume (which has a volume 15 times greater than that of the estuary), the net forcing is weaker than that in the estuary, and may be due to either tidal currents or wind stress depending on the time period considered. These forcings lead to irreversible mixing of the stratification (buoyancy flux) that turns river water into shelf water. This occurs in both the plume and estuary, but appears to be more efficient (17% vs. 5%), and somewhat greater (4.2 MW vs. 3.3 MW), in plume vs. estuary. This demonstrates the importance of both wind and tidal forcing to watermass transformation, and the need to consider the estuary and plume as part of a single system.     

Mixing,hypersalinity and gradients in Hervey Bay,Australia

Hervey Bay, a large coastal embayment situated off the central eastern coast of Australia, is a shallow tidal area (average depth = 15 m), close to the continental shelf. It shows features of an inverse estuary, due to the high evaporation rate (approx. 2 m/year), low precipitation (less than 1 m/year) and on average almost no freshwater input from rivers that drain into the bay. The hydro- and thermodynamical structures of Hervey Bay and their variability are presented here for the first time, using a combination of four-dimensional modelling and observations from field studies. The numerical studies are performed with the Coupled Hydrodynamical Ecological Model for Regional Shelf Seas (COHERENS). Due to the high tidal range (>3.5 m), the bay is considered as a vertically well-mixed system, and therefore, only horizontal fronts are likely. Recent field measurements, but also the numerical simulations, indicate characteristic features of an inverse/hypersaline estuary with low salinity (35.5 psu) in the open ocean and peak values (>39.0 psu) in the head water of the bay. The model further predicts a nearly persistent mean salinity gradient of 0.5 psu across the bay (with higher salinities close to the shore). The investigation further shows that air temperature, wind direction and tidal regime are mainly responsible for the stability of the inverse circulation and the strength of the salinity gradient across the bay. Due to an ongoing drying trend, the occurrence of severe droughts at the central east coast of Australia and, therefore, a reduction in freshwater supply, the salinity flux out of the bay has increased, and the inverse circulation has also strengthened.     

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.     

Numerical modeling of the penetration area of reversed submerged jet by the methods of finite elements and finite volumes

Water jet flow has many usages in the field of management and water resource operation that can be applied in mixing, dilution and aerification. The current study has calculated the flow velocity, length and height of jet penetration area (jet and main flow are of opposite directions) by the use of methods of finite element (FEM) and finite volume (FVM) and k–ε model. In order to evaluate and verify this turbulent model, the results of the numerical model have been compared with the experimental results. This model has been studied for consideration of various jet flow velocities and thicknesses. The conclusions have indicated that the length and the height of the penetration area have linear relationship with jet flow velocity; therefore, as the jet flow velocity increases, the length of jet penetration increases as well. The comparison of the results of numerical method with the experimental data have demonstrated that the FVM holds less convergence time and better results compared with FEM.     

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.     

THREE DIMENSIONAL MATHEMATICAL MODEL OF SEDIMENT TRANSPORT IN ESTUARINE REGIONS -A CASE STUDY OF THE HAIHE RIVER MOUTH

1 INTRODUCTIONUnderstanding the flow characters in open channel or other water area, whether natufal water bodies(such as rivers and estuaries) or man-made strUctures (such as navigational channels, harbors), isimportant for addressing numerous hydraulic engineering problems. These include the selechons ofsuitable waste disposal sites, contndnant transport, sediment transport and other ecological problems.The major difficulty for solving these problems is that it often requires an optimal…     

A simple approach to estimating freshwater discharge in branched estuarine systems

In tidal estuaries, quantifying freshwater discharge is still a difficult problem that has not yet been overcome due to the inherent difficulty in measuring and analysing the tidal discharge, especially during periods of low river flow. Because observations are often made in the stations further upstream, where the ratio of river to tidal discharge is large, it remains difficult to determine the discharge rate in the saline region. Freshwater discharge estimation is even more difficult in a branched estuary system having multiple diversion channels that connect with each other at a junction. To date, several methods have been developed for estimating freshwater discharge in estuaries. The most widely used are analytical and conceptual models that employ salinity as the principal trace and numerical simulations. However, these methods are very time consuming and costly as they require large sets of observations before the computations can take place. This paper presents a simple approach to investigating the discharge distribution over branched channels by considering the energy loss due to friction. We develop an analytical model that can obtain the discharge rate quantitatively at a junction where the main flow bifurcates into two branches. The model uses the bed roughness, tidal water level, and cross‐sectional profile under tidally averaged conditions as input data. Two selected estuarine systems in the Hiroshima delta in Japan and the Mekong delta in Vietnam have been investigated. Computations of the newly developed model show good agreement with earlier published results computed by sophisticated analytical and numerical models.     

Interaction between suspended sediment and tidal amplification in the Guadalquivir Estuary

Water level records at two stations in the Guadalquivir Estuary (Spain), one near the estuary mouth (Bonanza) and one about 77 km upstream (Sevilla), have been analysed to study the amplification of the tide in the estuary. The tidal amplification factor shows interesting temporal variation, including a spring-neap variation, some extreme low values, and especially the anomalous behaviour that the amplification factor is larger during a number of periods. These variations are explained by data analysis combined with numerical and analytical modelling. The spring-neap variation is due to the quadratic relation between the bottom friction and the tidal flow velocity. The river flood events are the direct causes of the extreme low values of the amplification factor, and they trigger the non-linear interaction between the tidal flow and suspended sediment transport. The fluvial sediment input during a river flood causes high sediment concentration in the estuary, up to more than 10 g/l. This causes a reduction of the effective hydraulic drag, resulting in stronger tidal amplification in the estuary for a period after a river flood. After such an event the tidal amplification in the estuary does not always fall back to the same level as before the event, indicating that river flood events have significant influence on the long-term development of this estuary.     

Implementation of a wetting-and-drying model in simulating the Kennebec–Androscoggin plume and the circulation in Casco Bay

A high-resolution coastal ocean model was developed to simulate the temporal/spatial variability of the Kennebec–Androscoggin (K–A) river plume and the circulation in Casco Bay. The model results agree favorably with the moored and shipboard observations of velocity, temperature, and salinity. The surface salinity gradient was used to distinguish the plume from the ambient coastal water. The calculated plume thickness suggests that the K–A plume is surface trapped. Its horizontal scales correlate well with Q ^0.25, where Q is the volume discharge of the rivers. Directional spreading is affected by the wind with the upwelling favorable wind transporting the plume water offshore. Both the wind and the tide also enhance mixing in the plume. The inclusion of a wetting-and-drying (WAD) scheme appears to enhance the mixing and entrainment processes near the estuary. The plume becomes thicker near the mouth of the estuary, the outflow velocity of the plume is weaker, and the radius of the river plume shrinks. The flow field in the model run with the WAD is noisier, not only in shallow areas of Casco Bay but also in the plume and even on the shelf. We speculate that the WAD processes can affect much larger areas than the intertidal zones, especially via a river plume that feeds into a coastal current.     

METHOD FOR DETERMINATION OF VELOCITY AND DEPTH FROM SEISMIC REFLECTION DATA1

The estimation of velocity and depth is an important stage in seismic data processing and interpretation. We present a method for velocity-depth model estimation from unstacked data. This method is formulated as an iterative algorithm producing a model which maximizes some measure of coherency computed along traveltimes generated by tracing rays through the model. In the model the interfaces are represented as cubic splines and it is assumed that the velocity in each layer is constant. The inversion includes the determination of the velocities in all the layers and the location of the spline knots. The process input consists of unstacked seismic data and an initial velocity-depth model. This model is often based on nearby well information and an interpretation of the stacked section. Inversion is performed iteratively layer after layer; during each iteration synthetic travel-time curves are calculated for the interface under consideration. A functional characterizing the main correlation properties of the wavefield is then formed along the synthetic arrival times. It is assumed that the functional reaches a maximum value when the synthetic arrival time curves match the arrival times of the events on the field gathers. The maximum value of the functional is obtained by an effective algorithm of non-linear programming. The present inversion algorithm has the advantages that event picking on the unstacked data is not required and is not based on curve fitting of hyperbolic approximations of the arrival times. The method has been successfully applied to both synthetic and field data.     

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.     

Factors affecting longitudinal dispersion in estuaries of different scale

Traditionally, the overall diluting capacity of an estuary is characterized using a coefficient of longitudinal dispersion, K_xe, which is given by the ratio of the dispersive flux of a dissolved substance to its tidally averaged longitudinal gradient. A steady-state model, which assumes a balance between the dispersive and non-dispersive fluxes and an exponential increase in estuary cross-sectional area towards the sea, has been used to derive expressions for the axial salinity distribution and the dispersive flux of salt. The model was set up assuming either a constant dispersion coefficient along the estuary or one that increased with distance towards the sea. By comparing salinity predictions with data from five UK estuaries, estimates were made of the maximum dispersive salt flux and the corresponding maximum salinity gradient for each system. The results indicated that there was an approximately linear relation between the fluxes and gradients, and the slope of a line plotted through the origin provided an estimate of a common K_xe for all five estuaries. The magnitude of K_xe was found to be about 90 m² s^–1 with a standard deviation of approximately ±32 m² s^–1. It is concluded that a representative value of 100 m² s^–1 for K_xe is a reasonable first choice when setting up a cross-sectionally averaged estuary model. The results also showed that larger systems, such as the Thames, had lower salinity gradients and lower dispersive salt fluxes, whilst smaller estuaries displayed the opposite characteristics. The model was used to predict the variation in the non-advective flux of salt along an estuary. The distribution was found to be similar to the corresponding flux distribution estimated from observations at the seaward end of the Tees estuary, despite appreciable spatial variations in the individual flux components. Allowing for a small decrease in freshwater flow, the model indicated that there was a decrease in the maximum dispersive flux between neap and spring tides. It is argued that such a reduction in flux can result in a seaward shift in the salinity distribution to a region of greater cross-section, where the freshwater transport per unit area again balances the reduced upstream dispersive flux, as found in the neap to spring response in the Tees estuary.Responsible Editor: Hans Burchard     

2D多尺度非线性地震速度成像

将遗传算法和单纯形算法相结合,得到了一种高效、健全的2D混合地震走时反演方法.把速度场划分为不同的空间尺度,定义网格节点上的速度作为待反演参数,采用双三次样条函数模型参数化,正问题采用有限差分走时计算方法,反问题采用多尺度混合反演方法.首先在较大的空间尺度内反演,然后减小空间尺度,将大尺度的反演结果作为次一级尺度反问题的初始模型,再进行混合反演,如此类推逐次逼近全局最优解.一个低速度异常体的数值模拟试验和抗走时扰动试验表明该方法是有效和健全的.我们将该方法应用到青藏高原东北缘阿尼玛卿缝合带东段上部地壳速度结构研究中,并与前人的成果进行了对比.