THREE-DIMENSIONAL NUMERICAL MODELING OF COHESIVE SEDIMENT TRANSPORT BY TIDAL CURRENT IN PEARL RIVER ESTUARV

IINTRODUCTIONEstUariesareprominentcoastalfeatUres.Estuariesareofgreateconomicssignificancetomankind.Attheseareas,manyharborsandwaterchannelshavetobebuiltforeconomicpurposes.ThedesignandconstrUctionofcoastalstrUctUresinestUariesrequireknowledgeofhydrodynamicsaswellassedimenttransportinsuchregions.ThenatUreofestuariesiscontrolledbyvariouscoastalhydrodynamicprocesses.Undertheactionofhydrodynamics,sedimentdepositionsorerosionswilloccurinestuariesornearcoastalstrUCtures.Tomaintainnavigati…     

A three-dimensional numerical model of particulate transport for coastal waters

A Lagrangian marker particle in Eulerian finite difference cell solution to the three-dimensional incompressible mass transport equation was developed for predicting particulate transport in coastal and estuarine waters. Special features of the solution procedure include a finite difference grid network which translates horizontally and vertically with the mean particle motion and expands with the dispersive growth of the marker particle cloud. The cartesian vertical coordinate of the three-dimensional mass transport equation has been transformed, using instantaneous water column depth to allow adaptation to flow situations with a temporally and spatially varying bottom topography and free surface. Results from this model for turbulent diffusion and advection of a uniform plug flow of sediment in an unbounded uniform flow field with various sediment settling velocities were in excellent agreement with the corresponding analytic solutions. Using current information from a two-dimensional vertically averaged hydrodynamic's model, the model was utilized to predict the long term diffusion and advection of dilute neutrally and negatively buoyant suspended sediment clouds resulting from a hypothetical instantaneous release of dredge spoil waste at Brown's Ledge in Rhode Island Sound.     

A mixed-dimensional finite volume method for two-phase flow in fractured porous media

We present a vertex-centered finite volume method for the fully coupled, fully implicit discretization of two-phase flow in fractured porous media. Fractures are discretely modeled as lower dimensional elements. The method works on unstructured, locally refined grids and on parallel computers with distributed memory. An implicit time discretization is employed and the nonlinear systems of equations are solved with a parallel Newton-multigrid method. Results from two-dimensional and three-dimensional simulations are presented.     

An improved method for evaluating the seasonal variability of total suspended sediment flux field in the Yellow and East China Seas

The suspended sediment flux field in the Yellow and East China Seas(YECS) displays its seasonal variability.A new method is introduced in this paper to obtain the flux field via retrieval of ocean color remote sensing data,statistical analysis of historical suspended sediment concentration data,and numerical simulation of three-dimensional(3D) flow velocity.The components of the sediment flux field include(i) surface suspended sediment concentration inverted from ocean color remote sensing data;(ii) vertical distribution of suspended sediment concentration obtained by statistical analysis of historical observation data;and(iii) 3D flow field modeled by a numerical simulation.With the improved method,the 3D suspended sediment flux field in the YECS has been illustrated.By comparison with the suspended sediment flux field solely based on the numerical simulation of a suspended sediment transport model,the suspended sediment flux field obtained by the improved method is found to be more reliable.The 3D suspended sediment flux field from ocean colour remote sensing and in situ observation are more closer to the reality.Furthermore,by quantitatively analyzing the newly obtained suspended sediment flux field,the quantity of sediment erosion and deposition within the different regions can be evaluated.The sediment exchange between the Yellow Sea and the East China Sea can be evident.The mechanism of suspended sediment transport in the YECS can be better understood.In particular,it is suggested that the long-term transport of suspended sediment is controlled mainly by the circulation pattern,especially the current in winter.     

THREE-DIMENSIONAL GEO-MORPHOLOGICAL MODELLING OF ESTUARINE WATERS

1 INTRODUCTION Estuaries and coastal zones have been used as means of navigation, disposal of waste material, fishing and many commercial and economic activities over the centuries. One of the most important phenomena in these regions is the suspended sediment transport, which may cause erosion and deposition, and hence changes in the estuarys morphology. In turn, such changes may lead to problems relating to navigation and estuarine management. When the bed boundary of an estuary change…     

Fluvial suspended sediment transport from cold and warm-based glaciers in Svalbard

An analysis of temporal variability in proglacial suspended sediment concentration is undertaken using time series data collected from three Svalbard basins which include one largely cold-based glacier (Austre Brøggerbreen), one largely warm-based glacier (Finsterwalderbreen) and one intermediate polythermal glacier (Erdmannbreen). The temporal variability in proglacial suspended sediment concentration is analysed using multiple regression techniques in which discharge is supplemented by other predictors acting as surrogates for variability in sediment supply at diurnal, medium-term and seasonal timescales. These multiple regression models improve upon the statistical explanation of suspended sediment concentration produced by simple sediment rating curves but need to account for additional stochastic elements within the time series before they may be considered successful. An interpretation of the physical processes which are responsible for the regression model characteristics is offered as a basis for comparing the different arctic glaciofluvial suspended sediment transport systems with that of their better known temperate glaciofluvial counterparts. It is inferred that the largely warm-based glacier is dominated by sediment supply from subglacial reservoirs which evolve in a similar manner to temperate glaciers and which cause a pronounced seasonal exhaustion of suspended sediment supply. The largely cold-based glacier, however, is dominated by sediment supply from marginal sources which generate a responsive system at short time scales but no significant seasonal pattern. The intermediate polythermal glacier basin, which was anticipated to be similar to the warm-based glacier, instead shows a highly significant seasonal increase in suspended sediment supply from an unusual subglacial reservoir emerging under pressure in the glacier foreland. The temperate model of glaciofluvial suspended sediment transport is therefore found to be of limited use in an arctic context. Copyright © 1999 John Wiley & Sons, Ltd.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

COMPUTATION OF THE ELECTRIC POTENTIAL IN THREE-DIMENSIONAL STRUCTURES*

A method has been developed for computation of the electrical DC potential of an arbitrary three-dimensional resistivity structure using a finite difference method. The threedimensionality is necessary for interpretation of geoelectrical soundings with controlled point sources over a laterally inhomogeneous medium. Lateral inhomogeneities should be considered in resistivity soundings with large layouts. The results obtained with the described method permit a more realistic representation of geological features. The resolution of the method is determined by the number of elements in the resistivity network. The problem of core memory space has been resolved by using random access disc files. The results computed using a Fortran program are in good agreement with analytically obtained solutions.     

模拟震源动力学破裂过程的数值方法研究

为模拟震源的动力学破裂过程,本文讨论一种求解剪切裂纹动力学扩展问题的有限差分法。在研究二维反平面破裂的基础上,我们把所用方法推广到三维问题。研究了许多简单震源模型,通过对比同一问题的解析解和其它数值解,对方法的正确性和计算精度进行了检验,结果表明我们的数值方法是可行的。最后,作为初步应用的例子,我们研究了非均匀断层的自发破裂问题。     

CASE STUDIES OF THREE-DIMENSIONAL NUMERICAL SIMULATION FOR TOTAL SEDIMENT TRANSPORT

1 INTRODUCTION With advances in computer technology and numerical methods, three-dimensional (3D) mathematical models for sediment transport are gradually applied more often and for more practical projects of hydraulic engineering. Three-dimensional mathematical river models can describe not only the secondary flow, but also the transport, deposition, and erosion of sediment in the river channel and the flood plain. However, at present the theories of sediment transport are not as well de…     

A MATHEMATICAL MODEL FOR UNSTEADY SEDIMENT TRANSPORT IN THE LOWER YELLOW RIVER

1 INTRODUCTIONMany mathematical models for sediment transport have been developed for solving practical Problemsin hydraulic engineering. HoweveT, most of them are not able to simulate the hyper-concentratd flows inthe Lower Yellow Xiver because of the extremely high load concentration of the flows. This paper isdevoted tO the simulation of unsteady sediment trallsport in the Lower Yellow mveLIn this paPer, the riverbed deformation equation is modified and the new expressions for sedi…     

模拟震源动力学破裂过程的数值方法研究

为模拟震源的动力学破裂过程,本文讨论一种求解剪切裂纹动力学扩展问题的有限差分法。在研究二维反平面破裂的基础上,我们把所用方法推广到三维问题。研究了许多简单震源模型,通过对比同一问题的解析解和其它数值解,对方法的正确性和计算精度进行了检验,结果表明我们的数值方法是可行的。最后,作为初步应用的例子,我们研究了非均匀断层的自发破裂问题。     

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.     

Optimal implicit staggered‐grid finite‐difference schemes based on the sampling approximation method for seismic modelling

We propose new implicit staggered‐grid finite‐difference schemes with optimal coefficients based on the sampling approximation method to improve the numerical solution accuracy for seismic modelling. We first derive the optimized implicit staggered‐grid finite‐difference coefficients of arbitrary even‐order accuracy for the first‐order spatial derivatives using the plane‐wave theory and the direct sampling approximation method. Then, the implicit staggered‐grid finite‐difference coefficients based on sampling approximation, which can widen the range of wavenumber with great accuracy, are used to solve the first‐order spatial derivatives. By comparing the numerical dispersion of the implicit staggered‐grid finite‐difference schemes based on sampling approximation, Taylor series expansion, and least squares, we find that the optimal implicit staggered‐grid finite‐difference scheme based on sampling approximation achieves greater precision than that based on Taylor series expansion over a wider range of wavenumbers, although it has similar accuracy to that based on least squares. Finally, we apply the implicit staggered‐grid finite difference based on sampling approximation to numerical modelling. The modelling results demonstrate that the new optimal method can efficiently suppress numerical dispersion and lead to greater accuracy compared with the implicit staggered‐grid finite difference based on Taylor series expansion. In addition, the results also indicate the computational cost of the implicit staggered‐grid finite difference based on sampling approximation is almost the same as the implicit staggered‐grid finite difference based on Taylor series expansion.     

Estimation of vertical water fluxes from temperature time series by the inverse numerical computer program FLUX‐BOT

The application of heat as a hydrological tracer has become a standard method for quantifying water fluxes between groundwater and surface water. The typical application is to estimate vertical water fluxes in the shallow subsurface beneath streams or lakes. For this purpose, time series of temperatures in the surface water and in the sediment are measured and evaluated by a vertical 1D representation of heat transport by advection and conduction. Several analytical solutions exist to calculate the vertical water flux from the measured temperatures. Although analytical solutions can be easily implemented, they are restricted to specific boundary conditions such as a sinusoidal upper temperature boundary. Numerical solutions offer higher flexibility in the selection of the boundary conditions. This, in turn, reduces the effort of data preprocessing, such as the extraction of the diurnal temperature variation from the raw data. Here, we present software to estimate water fluxes based on temperatures—FLUX‐BOT. FLUX‐BOT is a numerical code written in MATLAB that calculates vertical water fluxes in saturated sediments based on the inversion of measured temperature time series observed at multiple depths. FLUX‐BOT applies a centred Crank–Nicolson implicit finite difference scheme to solve the one‐dimensional heat advection–conduction equation. FLUX‐BOT includes functions for the inverse numerical routines, functions for visualizing the results, and a function for performing uncertainty analysis. We present applications of FLUX‐BOT to synthetic and to real temperature data to demonstrate its performance.     

A comparative study of three different learning algorithms applied to ANFIS for predicting daily suspended sediment concentration

The modeling and prediction of suspended sediment in a river are key elements in global water recourses and environment policy and management. In the present study, an Adaptive Neuro-Fuzzy Inference System model trained with the Levenberg-Marquardt learning algorithm is considered for time series modeling of suspended sediment concentration in a river. The model is trained and validated using daily river discharge and suspended sediment concentration data from the Schuylkill River in the United States. The results of the proposed method are evaluated and compared with similar networks trained with the common Hybrid and Back-Propagation algorithms, which are widely used in the literature for prediction of suspended sediment concentration. Obtained results demonstrate that models trained with the Hybrid and Levenberg-Marquardt algorithms are comparable in terms of prediction accuracy. However, the networks trained with the Levenberg-Marquardt algorithm perform better than those trained with the Hybrid approach.     

Modelling-based assessment of suspended sediment dynamics in a hypertidal estuarine channel

We investigate the dynamics of suspended sediment transport in a hypertidal estuarine channel which displays a vertically sheared exchange flow. We apply a three-dimensional process-based model coupling hydrodynamics, turbulence and sediment transport to the Dee Estuary, in the north-west region of the UK. The numerical model is used to reproduce observations of suspended sediment and to assess physical processes responsible for the observed suspended sediment concentration patterns. The study period focuses on a calm period during which wave-current interactions can reasonably be neglected. Good agreement between model and observations has been obtained. A series of numerical experiments aim to isolate specific processes and confirm that the suspended sediment dynamics result primarily from advection of a longitudinal gradient in concentration during our study period, combined with resuspension and vertical exchange processes. Horizontal advection of sediment presents a strong semi-diurnal variability, while vertical exchange processes (including time-varying settling as a proxy for flocculation) exhibit a quarter-diurnal variability. Sediment input from the river is found to have very little importance, and spatial gradients in suspended concentration are generated by spatial heterogeneity in bed sediment characteristics and spatial variations in turbulence and bed shear stress.     

MODELLING FLOOD HYDRAULICS AND OVERBANK DEPOSITION ON RIVER FLOODPLAINS

This paper outlines a numerical model for the prediction of floodplain inundation sequences, overbank deposition rates and deposit grain size distributions. The model has two main components: first, a simplified hydraulic scheme which predicts floodwater flow depths and velocities, and second, a sediment transport element which employs a mass balance relation describing suspended sediment dispersion by convective and diffusive processes and sediment deposition as a function of particle settling rates. These relationships are solved numerically on a finite difference grid that accurately replicates the complex topographic features typical of natural river floodplains. The model is applied to a 600 m reach of the River Culm, Devon, U.K. using data derived from a range of field and laboratory techniques. Continuous records of river stage and suspended sediment concentration provide the model's upstream boundary input requirements. These are supplemented by measurements of the in situ settling characteristics of the suspended sediment load. The model's sediment transport component is calibrated with the aid of a dataset of measured overbank deposition amounts derived from flood events over a 16 month period. The model is shown to predict complicated floodwater inundation sequences and patterns of suspended sediment dispersion and deposition, which are largely a product of the complex topography of the floodplain. These results compare favourably with observations of overbank processes and are an improvement over those of previous models which have employed relatively simple representations of floodplain geometry. © 1997 by John Wiley & Sons, Ltd.     

Effects of hydrological processes on the chemical composition of riverine suspended sediment in the Zhujiang River,China

The chemical composition of riverine suspended sediment is the integration of the weathering crust minerals, soil organic matter and erosion agency within a specific drainage basin, which has been largely disturbed by the human activities. Selected metal elements of the riverine suspended sediment in the Zhujiang River were analysed using inductively coupled plasma–atomic emission spectrometry (ICP–AES) in three different hydrological phases from 1997 to 1998 at Makou and Sanshui hydrographic gauge stations, located at the lower reaches of the two main tributaries of the Zhujiang River, i.e. the Xijiang and the Beijing Rivers respectively. Organic carbon and nitrogen were also analysed using a conventional element analyser. The results demonstrate that the chemical composition of the riverine suspended sediment show obvious variability in different hydrological phases, which closely correlate to the organic matter content in suspended sediment. Intensified erosion in the flood phase results in lower concentration of the organic matter than that in the lower water level phase. The riverine suspended sediment with rich organic matter in the lower water level phase adsorbs some metal elements from the river water. Copyright © 2003 John Wiley & Sons, Ltd.     

Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems

GeoFEST (Geophysical Finite Element Simulation Tool) is a two- and three-dimensional finite element software package for the modeling of solid stress and strain in geophysical and other continuum domain applications. It is one of the featured high-performance applications of the NASA QuakeSim project. The program is targeted to be compiled and run on UNIX systems, and is running on diverse systems including sequential and message-passing parallel systems. Solution to the elliptical partial differential equations is obtained by finite element basis sampling, resulting in a sparse linear system primarily solved by conjugate gradient iteration to a tolerance level; on sequential systems a Crout factorization for the direct inversion of the linear system is also supported. The physics models supported include isotropic linear elasticity and both Newtonian and power-law viscoelasticity, via implicit quasi-static time stepping. In addition to triangular, quadrilateral, tetrahedral and hexahedral continuum elements, GeoFEST supports split-node faulting, body forces, and surface tractions. This software and related mesh refinement strategies have been validated on a variety of test cases with rigorous comparison to analytical solutions. These include a box-shaped domain with imposed motion on one surface, a pair of strike slip faults in stepover arrangement, and two community-agreed benchmark cases: a strike slip fault in an enclosing box, and a quarter-domain circular fault problem. Scientific applications of the code include the modeling of static and transient co- and post-seismic earth deformation, Earth response to glacial, atmospheric and hydrological loading, and other scenarios involving the bulk deformation of geologic media.