Predicting Aeolian Sand Transport: Revisiting the White Model

The derivation and history of the frequently cited aeolian transport model of White are considered in light of the continued replication of an error in the original expression. The error may have escaped notice because the expression is still dimensionally correct and it yields predictions that appear reasonable in comparison with both the predictions of other models with field data. The incorrect expression has come to be identified as a distinct model. However, the correct formulation of the ‘White model’ is, in fact, a rearrangement of the Kawamura model with a slightly smaller (c.6%) empirical coefficient. © 1997 John Wiley & Sons, Ltd.     

A continuous model for sand dunes: Review,new developments and application to barchan dunes and barchan dune fields

Basically, sand dunes are patterns resulting from the coupling of hydrodynamic and sediment transport. Once grains move, they modify the surface topography which in turns modifies the flow. This important feedback mechanism lies at the core of continuous dune modelling. Here we present an updated review of such a model for aeolian dunes, including important modifications to improve its predicting power. For instance, we add a more realistic wind model and provide a self‐consistent set of parameters independently validated. As an example, we are able to simulate realistic barchan dunes, which are the basic solution of the model in the condition of unidirectional flow and scarce sediments. From the simulation, we extract new relations describing the morphology and dynamics of barchans that compare very well with existing field data. Next, we revisit the problem of the stability of barchan dunes and argue that they are intrinsically unstable bed‐forms. Finally, we perform more complex simulations: first, a barchan dune under variable wind strength and, second, barchan dune fields under different boundary conditions. The latter has important implications for the problem of the genesis of barchan dunes. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluating overland flow sediment transport capacity

An equation for evaluating the sediment transport capacity of overland flow is a necessary part of a physically based soil erosion model describing sediment detachment and transport as distributed processes. At first, for the hydraulic conditions of small-scale and large-scale roughness, the sediment transport capacity relationship used in the WEPP model is calibrated by Yalin and Govers' equation. The analysis shows that the transport coefficient K_t depends on the Shields parameter, Y, according to a semi-logarithmic (Yalin) or a linear (Govers) equation. The reliability of the semi-logarithmic equation is verified by Smart's, and Aziz and Scott's experimental data. Then the Low's formula, whose applicability is also proved by Smart's, and Aziz and Scott's data, is transformed as a stream power equation in which a stream power coefficient, K_SP, depending on Shields parameter, slope, sediment and water-specific weight, appears. A relationship between transport capacity and effective stream power is also proposed. Finally, the influence of rainfall on sediment transport capacity and the prediction of critical shear stress corresponding to overland flow are examined. © 1998 John Wiley & Sons, Ltd.     

Two-dimensional modeling of contaminant transport in porous media in the presence of colloids

It has long been known that colloids can facilitate the transport of contaminants in groundwater systems by reducing the effective retardation factor. A significant effort has been devoted to study colloid-facilitated contaminant transport during the past decade. Many of the previous studies were restricted to one-dimensional analyses and comparisons with finite-column experiments. In this work, a two-dimensional numerical model is developed and used to study the different interactions between colloids, contaminants, and porous media under homogeneous conditions. The numerical formulation of the model is based on discretizing mass balance equations and reaction equations using finite differences having a third-order, total variance-diminishing scheme for the advection terms. This scheme significantly reduces numerical dispersion and leads to greater accuracy compared to the standard central-differencing scheme. The model is tested against analytical solutions under simplified conditions as well as against experimental data, and the results are favorable. The model is used to investigate the impact of the various reaction rates and parameter values on the movement of contaminant plumes in two dimensions. The model is also used to investigate the hypothesis that colloids may increase the effective retardation factor of contaminant plumes. The analysis shows that assuming kinetic mass exchange between contaminant and colloids with constant reaction rate coefficients that are not related to the concentrations may lead to inaccurate results. These inaccurate results are exemplified in the finding that under the kinetic assumption the ratio of the initial concentration of colloids to the initial concentration of contaminant does not affect the amount of facilitation or retardation that occurs in the system. It is also found that colloids can increase the effective retardation factor for the contaminant under certain combinations of reaction rates and distribution coefficients. A quantitative empirical expression to identify whether colloids retard or facilitate the contaminant movement is presented.     

东亚地区春季黑碳气溶胶源排放及其浓度分布

本研究收集整理了东亚地区一个黑碳气溶胶的排放源,并利用一个三维化学输送模式系统模拟研究了东亚地区春季黑碳气溶胶的输送和演变过程.黑碳气溶胶浓度的模拟值与地面观测站和观测船获取的监测值相比,模拟值与监测值具有良好的一致性,这表明本研究收集整理的黑碳气溶胶的排放量及其分布是基本合理的,模式系统较好地反映了黑碳气溶浓度的时空分布和变化趋势,再现了许多观测到的重要特征.模拟结果表明,东亚地区春季间隙性的长距离输送对于西北太平洋海域的黑碳气溶胶浓度分布具有非常重要的影响.     

Dispersion dependence on retardation in a real fracture geometry using lattice-gas cellular automaton

Fractures have been recently identified in potential host rock for high level nuclear waste disposal, like indurated argilite formations. These fractures appear as potential rapid pathways for radionuclides transport and hydrodynamic properties of the transport inside these systems must thus be characterized. Miscible non-sorbing and sorbing tracers displacements were performed on a 2-D model derived from a real fracture geometry observed in the Tournemire argilite formation with a lattice-gas cellular automaton (LGA). LGA was shown to easily handle the complex geometry of such a fracture. The numerical breakthrough curves obtained were inverted with the 1-D CDE and MIM transport models. Two main conclusions were drawn: (i) at the length scale of the study, the non-sorbing tracer transport in our fracture geometry was more accurately interpreted in terms of the MIM model rather than in terms of the classical CDE model; (ii) in order to correctly model the sorbing tracers migration, the hydrodynamic dispersion coefficient value was found to increase with the increase of the retardation factor. A semi-empirical relation based on the Taylor–Aris theory was then used to describe this dependency.     

Interbasin freshwater,heat and salt transport through the boundaries of the East and South China Seas from a variable-grid global ocean circulation model

It has long been recognized that the circulation in the East China Sea (ECS) and Japan/East Sea (JES) is closely related with that in Pacific, especially with the Kuroshio (e.g., Nitani[1], Hi-daka[2]). Based on current measurements in the Taiwan Strait a…     

Modeling Sediment Transport Dynamics in Thompson Island Pool,Upper Hudson River

Two-dimensional, vertically-averaged hydrodynamic and sediment transport models were developed and applied as part of a PCB fate and transport modeling study of Thompson Island Pool (TIP), Upper Hudson River. Mechanistic formulations were used to simulate cohesive and non-cohesive suspended load transport; site-specific data were extensively used to determine model inputs. This modeling approach is compared and contrasted to non-mechanistic solids transport sub-models used in other contaminant fate studies. A minimum number of model parameters were adjusted to calibrate the sediment transport model using data collected during the 1994 spring flood. The model was validated during the 1997 spring flood and for a 22-year (1977–1998) period. Successful calibration and validation of the model showed that: (1) deposition and resuspension processes were realistically and accurately formulated in the model; (2) the model is an effective diagnostic tool for quantitatively evaluating net deposition and erosion from various areas of TIP; and (3) sediment transport results can be coupled with a PCB fate model with a high degree of confidence.     

A two-dimensional horizontal wave propagation and mud mass transport model

The present study offers a two-dimensional horizontal wave propagation and morphodynamic model for muddy coasts. The model can be applied on a general three-dimensional bathymetry of a soft muddy coast to calculate wave damping, fluid mud mass transport and resulting bathymetry change under wave actions. The wave propagation model is based on time-dependent mild slope equations including the wave energy dissipation due to the wave-mud interaction of bottom mud layers as well as the combined effects of the wave refraction, diffraction and breaking. The constitutive equations of the visco-elastic–plastic model are adopted for the rheological behavior of fluid mud. The mass transport velocity within the fluid mud layer is calculated combining the Stokes’ drift, the mean Eulerian velocity and the gravity-driven mud flow. The results of the numerical model are compared against a series of conducted wave basin experiments, wave flume experiments and field observations. Comparisons between the computed results with both the field and laboratory data reveal the capability of the proposed model to predict the wave transformation and mud mass transport.     

Sediment transport at the network scale and its link to channel morphology in the braided Vjosa River system

In this article we apply the CASCADE network-scale sediment connectivity model to the Vjosa River in Albania. The Vjosa is one of the last unimpaired braided rivers in Europe and, at the same time, a data scarce environment, which limits our ability to model how this pristine river might respond to future human disturbance. To initialize the model, we use remotely sensed data and modeled hydrology from a regional model. We perform a reach-by-reach optimization of surface grain size distribution (GSD) and bedload transport capacity to ensure equilibrium conditions throughout the network. In order to account for the various sources of uncertainty in the calculation of transport capacity, we performed a global sensitivity analysis. The modeled GSD distributions generated by the sensitivity analysis generally match the six GSDs measured at different locations within the network. The modeled bedload sediment fluxes increase systematically downstream, and annual fluxes at the outlet of the Vjosa are well within an order of magnitude of fluxes derived from previous estimates of the annual suspended sediment load. We then use the modeled sediment fluxes as input to a set of theoretically derived functions that successfully discriminate between multi-thread and single-thread channel patterns. This finding provides additional validation of the model results by showing a clear connection between modeled sediment concentrations and observed river morphology. Finally, we observe that a reduction in sediment flux of about 50% (e.g., due to dams) would likely cause existing braided reaches to shift toward single thread morphology. The proposed method is widely applicable and opens a new avenue for application of network-scale sediment models that aid in the exploration of river stability to changes in water and sediment fluxes.     

Incremental distributed modelling investigation in a small agricultural catchment: 2. Erosion and phosphorus transport

Distributed physically based erosion and phosphorus (P) transport models, run by the overland flow model described in Taskinen and Bruen (2006. Hydrological Processes 20 : this issue), are described. In the erosion model, the additional components to the basic model were the outflow of the particles by infiltration and a new model component, i.e. deposition when rainfall stops. Two ways of calculating the shielding factor due to the flow depth were compared. The P transport model had both dissolved P (DP) and particulate P (PP) components. The processes included in the DP model were desorption from the soil surface, advection, storage in the overland flow and infiltration. The PP model accounted for advection, storage in the flow, infiltration, detachment from the soil surface by flow and rainfall and deposition both when transport capacity of suspended solids (SS) is exceeded and when rainfall ceases. When the models were developed and validated in small agricultural fields of cohesive soil types in southern Finland, comparisons were made between corresponding processes and the significance of added components were estimated in order to find out whether increased model complexity improves the model performance. The sedigraphs were found to follow the dynamics of rainfall, emphasizing the importance of the rainfall splash component. The basic model was too slow to react to changes in rainfall and flow rates, but infiltration and deposition that acts during the cessation in rainfall improved the model significantly by enabling the modelled SS to fall sharply enough. The shielding effect of flow depth from the splash detachment was found to play a significant role. Transport capacity should also be included in erosion models when they are applied to cohesive soils. In this study, the Yalin method worked well. A strong correlation was obvious between the measured SS and total P concentrations, indicating that the main form of P in runoff is PP. This emphasizes the importance of a good sediment transport model in P transport modelling. The submodel used for DP desorption from the soil surface produced plausible results without any calibration. Copyright © 2006 John Wiley & Sons, Ltd.     

One‐dimensional transitional behaviour in saltation

One‐dimensional simulations of the unsteady saltation process show that the transport rate's response depends on the amplitude and frequency of the wind fluctuations. At frequencies higher than f ≈ 0·5 Hz the transport rate was found not to respond to the wind changes. The initial overshoot reported by previous investigators was found not to appear for simulation heights smaller than 50 to 60 cm. This is due to the fast propagation of the grains' influence upward in the flow and the immediate deceleration of the wind. Confirmation of these findings comes from reports of experiments conducted in wind tunnels of different sizes. Further test calculations show that the discretization time step size Δt has an influence on the model's temporal behaviour. The reason for this is the better coupling of the wind–sand system when a smaller Δt is used. The implications of bed area choice on the statistical accuracy of predicted transport rate is also demonstrated. In the one‐dimensional case the grain cloud's total forward momentum equals transport rate, which is independent of model geometry. Copyright © 2000 John Wiley & Sons, Ltd.     

On the temporal resolution of mass balance models for soluble chemicals in soils

A daily step model of chemical mass balance in the topsoil is presented and validated at the three experimental sites in Europe, and subsequently applied to perform two distinct numerical experiments. First, an experiment was run using hypothetic soluble chemicals with half‐lives ranging from 10^?1 to 10⁴, with a range of representative European climate and soil properties, assuming uniform constant emissions of the chemicals throughout the year. Chemical mass in soil from the daily step model calculations can be surrogated by the monthly step model consistently parameterized in terms of absolute values, patterns and inter‐monthly variability with decreasing accuracy at higher chemical half‐lives. Leaching fluxes can be also surrogated by the monthly step calculation, although with higher errors. Runoff is correct in the order of magnitude, but it shows only a weak correlation with the monthly mean of the daily model output. For leaching and runoff, the accuracy depends mainly on soil properties. Variability is well reproduced for both leaching and runoff. The second experiment represented a pulse emission of chemicals discharged on a single day in a 12‐month period. Results from the annual average mass of chemicals in the soil, annual runoff and leaching fluxes from the daily step model were compared with the results obtained from the experiment assuming constant‐removal rates for the year. The two values are within a factor of 10 for half‐lives longer than 10 days; therefore, it is possible to emulate the daily step model with a simple constant‐removal rate model for screening‐level assessment. The experiments suggest that simpler schemes may be a practical screening‐level approximation of detailed daily step models for both continuous and pulse emissions, two cases providing extreme bounds of variation to real world emissions. Copyright © 2010 John Wiley & Sons, Ltd.     

A process-based model for sediment transport under various wave and current conditions

The purpose of this study is to investigate the capability of a newly developed process-based model for sediment transport under a wide variety of wave and current conditions.The model is based on the first-order boundary layer equation and the sediment advection-diffusion equation.In particular,a modified low Reynolds number k-e model is coupled to provide the turbulence closure.Detailed model verifications have been performed by simulating a number of laboratory experiments,covering a considerable range of hydrodynamic conditions such as sinusoidal waves,asymmetric waves and wave-current interactions.The model provides satisfactory numerical results which agree well with the measured results,including the time-averaged/dependent sediment concentration profiles and sediment flux profiles,as well as the time series of concentration at given elevations.The observed influences of wave orbital velocity amplitude,wave period and sediment grain size are correctly reproduced,indicating that the fundamental physical mechanisms of those processes are properly represented in the model.It is revealed that the present model is capable of predicting sediment transport under a wide range of wave and current conditions,and can be used to further study the morphodynamic processes in real coastal regions.     

An ELLAM approximation for advective–dispersive transport with nonlinear sorption

We consider an Eulerian–Lagrangian localized adjoint method (ELLAM) applied to nonlinear model equations governing solute transport and sorption in porous media. Solute transport in the aqueous phase is modeled by standard advection and hydrodynamic dispersion processes, while sorption is modeled with a nonlinear local-equilibrium model. We present our implementation of finite volume ELLAM (FV-ELLAM) and finite element (FE-ELLAM) discretizations to the reactive transport model and evaluate their performance for several test problems containing self-sharpening fronts.     

Numerical analysis of bacterial transport in saturated porous media

A mathematical model to describe bacterial transport in saturated porous media is presented. Reversible/irreversible attachment and growth/decay terms were incorporated into the transport model. Additionally, the changes of porosity and permeability due to bacterial deposition and/or growth were accounted for in the model. The predictive model was used to fit the column experimental data from the literature, and the fitting result showed a good match with the data. Based on the parameter values determined from the literature experimental data, numerical experiments were performed to examine bacterial sorption and/or growth during bacterial transport through saturated porous media. In addition, sensitivity analysis was performed to investigate the impact of key model parameters for bacterial transport on the permeability and porosity of porous media. The model results show that the permeability and porosity of porous media could be altered due to bacterial deposition and growth on the solid matrix. However, variation of permeability due to bacterial growth was trivial compared with natural permeability variation. Copyright © 2005 John Wiley & Sons, Ltd.     

APPLICATION OF THE MULTIGRID METHOD IN NUMERICAL SIMULATION OF SEDIMENT TRANSPORT

IINTRODUCTIONThemultigridmethod(Brandt,1988)iswidelyusedinmanydisciplinesandengineeringsubjects,especiallyintheareaofComputationalFluidMechanicsbecauseofitshighefficiencyincomputation.Themultigridmethodhasbeenprovenintheorytobeoneofbestnumericalmethodsforsolvinglinearellipticproblemsbecausetheloadofcomplltationisindirectproportiontothenumberofgridnodes;theconvergentspeedofmultigridmethodisirrelevanttothesizeofgird.Soitissuitableforthenumericalsimulationoffullscaleengineeringprojects.The…     

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…