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.     

Improving hydrodynamic modeling of river networks by incorporating data assimilation using a particle filter

Numerical modeling is a well-recognized method for studying the hydrodynamic processes in river networks. Multi-source measurements also offer abundant information on the patterns and mechanisms within the processes. Therefore, improving hydrodynamic modeling of river networks through the use of data assimilation techniques has become a hot research topic in recent years. The particle filter (PF) is a commonly used data assimilation method and has been proven to be applicable to various nonlinear and non-Gaussian models. In the current study, an improved numerical hydrodynamic model for large-scale river networks is established by incorporating the advanced PF algorithm. Furthermore, the PF method based on the Gaussian likelihood function (GLF) and the method based on the Cauchy likelihood function (CLF) are compared for a complex river network scenario. The feasibility of the PF-based methods was evaluated through application to the Yangtze-Dongting River-lake Network (YDRN) by assimilating water stage data collected at six hydrometric stations during the entire hydrodynamic process in 2003. Additionally, the parameters used in the likelihood function, which affect the assimilation performance, also were explored in the current study. The study results found that the accuracy of the model-derived water stage data was improved when the PF-based methods are utilized, with improvement not only at the data assimilation (calibration) sites but also at three hydrometric stations not used in the data assimilation (i.e., verification sites). The highest average Nash-Sutcliffe Efficiency result for the six assimilation sites were 0.98 while the lowest summed root-mean-square-error result was 1.801 m. The comparison results also indicated that the CLF-based PF outperformed the GLF-based PF when high-accuracy observed data are available. Specifically, the CLF can effectively resolve the filtering failure problem and the dispersion problem of PFs, and further improve the accuracy of the filtering results for a river network scenario. In summary, the CLF-based PF method along with high-accuracy observation data shows promise to provide reliable reference and technical support for hydrodynamic modeling of large-scale river networks.     

Fine grid tidal modeling of the Yellow and East China Seas

A fine grid tidal modeling experiment is carried out in order to investigate the tidal regimes for major five tidal constituents, the nonlinear tidal phenomena in terms of M₄ and MS₄ generation, and the independent tide by the tide generating force in the Yellow and East China Seas (YECS). In this study a two-dimensional numerical model based upon a subgrid-scale (SGS) stress modeling technique is used with the tide generating force included. The model was validated with recently observed tide and current data. The calculated tidal charts for tidal elevation show a generally good agreement with existing ones, with more accurate feature of the M₂ cotidal chart in comparison with both the observed data and the existing tidal charts. A careful comparison of the computed diurnal amplitude with observations suggests that the diurnal constituents seem to be overdamped especially in the Kyunggi Bay region, for the case when quadratic bottom friction law is used.Propagation features of the M₄(MS₄) tides are discussed in the YECS, based upon the analyses of the observed and calculated results. The amphidromic system of the M₄ is quite complicated and one noticeable characteristic is that the propagation direction of the M₄ tidal wave along the west coast of Korean peninsula is opposite to that of the M₂ tidal wave. This result coincides with observations. The propagation feature of the MS₄ is almost similar to that of the M₄, but with lesser amplitude. The responses of the M₄ tidal features to momentum diffusion term and depth-dependent form of the friction coefficient are also discussed.It is also shown that when the independent tide (Defant, 1960) arising from tide generating force (TGF) coexists with tidal waves (co-oscillating tide) arising from external boundary forcing, the TGF tidal waves are dissipated and their amphidromes tend to move westward. This may be interpreted as a process whereby the incident and reflected TGF tidal waves are damped by co-oscillating tide arising from external force at open boundaries. The TGF amplitude is found to be up to 10 cm and 4 cm in the Kyunggi Bay for the M₂ and S₂ constituents while those for all diurnal constituents are less than 1 cm over the entire model domain.     

Effect of tidal forcing on a subterranean estuary

Groundwater flow and chemical transport in subterranean estuaries are poorly understood despite their potentially important implications for chemical fluxes from aquifers to coastal waters. Here, a numerical study of the dynamics in a subterranean estuary subject to tidal forcing is presented. Simulations show that salt transport associated with tidally driven seawater recirculation leads to the formation of an upper saline plume in the intertidal region. Computed transit times and flow velocities indicate that this plume represents a more active zone for mixing and reaction than the dispersion zone of the lower, classical salt wedge. Proper conceptualisation of this surficial mixing zone extends our understanding of processes within the subterranean estuary. Numerical tracer simulations reveal that tidal forcing may reduce the threat of a land-derived contaminant discharging to the marine environment by modifying the subsurface transport pathway and local geochemical conditions. Mixing and stratification in the subterranean estuary are strongly affected by both inland and tidal forcing. Based on the estuarine analogy we present a systematic classification of subterranean estuaries.     

Three-dimensional modelling of estuarine turbidity maxima in a tidal estuary

A new numerical model for simulating estuarine dynamics is introduced here. This model, called General Estuarine Transport Model (GETM), has been specifically designed for reproducing baroclinic, bathymetry-guided flows where the tidal range may exceed the mean water depth in large parts of the domain such that drying and flooding processes are relevant. Several physical and numerical features of the model support exact and stable results for such domains. For the physics, high-order turbulence closure schemes guarantee proper reproduction of vertical exchange processes. Among the specific numerical features, generalised vertical coordinates, orthogonal curvilinear horizontal coordinates, high-order TVD advection schemes and stable drying and flooding algorithms have been implemented into GETM. The model is applied here to simulate the dynamics of estuarine turbidity maxima (ETMs), a complex feature present in most tidal estuaries. First, idealised simulations for a two-dimensional domain in the xz space will be shown to reproduce the basic generation mechanisms for ETMs. Then, a realistic three-dimensional simulation of the Elbe estuary in Northern Germany will be carried out. It is demonstrated that for a given forcing situation the model reproduces a stable ETM at the correct location.Responsible Editor: Phil Dyke     

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary’s entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429–35) for two-dimensional non-interacting tidal waves.     

A Lagrangian-trajectory study of a gradually mixed estuary

When modelling is used for investigating estuarine systems, a choice generally has to be made between applying simple mass-balance considerations or using a process-resolving three-dimensional (3-D) numerical circulation model. In the present investigation of the Gulf of Finland, a gradually mixed estuary in the Baltic Sea, it is demonstrated how Lagrangian-trajectory analysis applied to the output from a 3-D model minimizes the disadvantages associated with both of the modelling techniques referred to above. This formalism made it possible to demonstrate that the main part of the Gulf is dominated by water originating from the Baltic proper, and that the most pronounced mixing with fresh water from the river Neva takes place over a limited zone in the inner part of the Gulf. Dynamical insights were furthermore obtained by using the Lagrangian formalism to construct overturning stream-functions for the two source waters.     

Interactions of river discharge and tidal modulation in a tropical estuary, NE Brazil

Observations of thermohaline properties and currents were undertaken in the Curimataú River estuary (6°18′S), Rio Grande do Norte state (RN), Brazil, during consecutive neap–spring tidal cycles in the austral autumn rainy season. Highly asymmetric neap tide along channel velocities (−0.4 to 0.9 m s⁻¹) and highly stratified conditions were generated by an increase of the buoyancy energy from the freshwater input (R _iE≈5.6). During the spring-tidal cycle the river discharge decreased and the longitudinal velocity components were higher, less asymmetrical (−0.8 to 1.1 m s⁻¹) and semidiurnal, associated with moderately stratified conditions (R _iE≈0.1) due to the increase of the kinetic tidal energy forcing mechanism. The overall salinity variation from surface to bottom during two tidal cycles was from 20.5 to 36.3 and 29 to 36.7 in the neap and spring tide experiments, respectively; in the last experiment, the tropical water (TW) mass intrusion was enhanced. The net salt transport reversed from down to up estuary during the neap and spring tide experiments, respectively, varied from 6.0 to –2.0 kg m⁻¹ s⁻¹, an indication of changes in the main forcing of the estuary dynamics. Evaluation of a classical steady analytical model, in comparison with nearly steady experimental vertical profiles of velocity, shows an agreement classifiable as reasonably fair.     

Benchmark on the numerical simulations of the hydrodynamic and morphodynamic evolution of tidal channels and tidal inlets

A correct understanding of the hydrodynamics and morphodynamics of tidal basins is of fundamental importance for the fate of the Venice Lagoon, Italy. If on one hand, the development of sophisticated numerical models is called for in order to reproduce the complexity of the mechanisms governing the morphodynamic evolution of many natural environments, including lagoons, on the other hand, a clear knowledge of the reliability and limits of the results provided by these models is crucial in order to establish the condition under which they can be safely applied. To this aim, researchers involved in numerical modeling in the framework of the recent Corila research programmes, agreed to perform an accurate comparison of results provided by three different numerical models, applying them to the test case offered by the experimental investigations performed under controlled conditions by Tambroni et al. (2005a). Here, we consider the following numerical models: (i) a 2D finite element hydrodynamic model coupled with a 2D finite volume morphodynamic model (5 and 3); (ii) a 2D finite element morphodynamic model (Ferrarin et al., 2008); (iii) a 2D depth-averaged model for the inlet region, coupled with a 1D model for the channel (Tambroni et al., 2005b). A first set of simulations concerns the fixed bed case and shows that all the models provide similar results: in particular, they are able to predict the observed free surface oscillations satisfactorily, while comparison with the measured velocity field is less satisfactory. Moreover, as far as the flow field at the inlet is concerned, the models describe accurately the potential flow into the channel during the flood phase, while they are not able to adequately reproduce the occurrence of the fine structure of the shear layers shed by the inlet edges during the ebb phase. This limit is related to the shallow water character of the models. As for the morphodynamics, the long term equilibrium configurations of the bottom of the channel and of the near inlet region show qualitative agreement with the experimental observations, although in this case the differences between the results provided by the distinct numerical approaches are more marked.     

The attenuation of tidal and subtidal oscillations in the Patos Lagoon estuary

The single, long and narrow channel that usually connects choked coastal lagoons to the ocean can serve as a natural hydraulic low-pass filter that reduces or eliminates tidal and subtidal effects inside the lagoon. This study proposes an alternative method of estimating the attenuation of the tidal and subtidal oscillations throughout the Patos Lagoon estuary. The attenuation is estimated for conditions of contrasting river runoff and weather (summer and winter). A high-pass/low-pass filter (fast fourier transformation technique – FFT) is applied to time series of sea-surface elevation (SSE) measured at the mouth of the Patos Lagoon. The resulting high-frequency (tidal) and low-frequency (subtidal) signals are used in independent simulations to force the TELEMAC-2D model. Attenuation of the tidal and subtidal signals throughout the estuary is estimated by applying cross-spectral analysis between the model-generated SSE time series at different locations throughout the estuary and the filtered SSE time series measured at the mouth. Results from the proposed method suggest that: (1) the low-frequency (subtidal) oscillations are less attenuated and propagate further than the high-frequency (tidal) oscillations in the Patos Lagoon estuary; (2) the filtering capability of the Patos Lagoon estuary is expected to follow a seasonal pattern, although further investigations on an interannual time scale are recommended in order to confirm this hypothesis; (3) the influence of the oceanic boundary processes on the SSE dynamics of the lagoon is restricted to the lower estuary. Further inland, the local forcing generated by the wind and freshwater input is likely to be the main forcing effect controlling the dynamics of the system. The proposed method proved to be an efficient and alternative way of estimating the attenuation of energy in the tidal and subtidal bands throughout the access channel of a choked coastal lagoon located in an area of reduced tidal influence.Responsible Editor: Iris Grabemann     

Trapped internal waves over undular topography in a partially mixed estuary

The flow of a stratified fluid over small-scale topographic features in an estuary may generate significant internal wave activity. Lee waves and upstream influence generated at isolated topographic features have received considerable attention during the past few decades. Field surveys of a partially mixed estuary, the Rotterdam Waterway, in 1987, also showed a plethora of internal wave activity generated by isolated topography, banks and groynes. Additionally it revealed a spectacular series of resonant internal waves trapped above low-amplitude bed waves. The internal waves reached amplitudes of 3–4 m in an estuary with a mean depth of 16 m. The waves were observed during the decreasing flood tide and are thought to make a significant contribution to turbulence production and mixing. However, while stationary linear and finite amplitude theories can be used to explain the presence of these waves, it is important to further investigate their time-dependent and non-linear behaviour. With the development of advanced non-hydrostatic models it now becomes possible to further investigate these waves through numerical experimentation. This is the focus of the work presented here. The non-hydrostatic finite element numerical model FINEL3D developed by Labeur was used in the experiments presented here. The model has been shown to work well in a number of stratified flow investigations. Here, we first show that the model reproduces the field data and for idealised stationary flow scenarios that the results are in agreement with the resonant response predicted by linear theory. Then we explore the effects of non-linearity and time dependence and consider the importance of resonant internal waves for turbulence production in stratified coastal environments.Responsible Editior: Hans Burchard     

Turbulence and mixing in a freshwater-influenced tidal bay: Observations and numerical modeling

In situ observations and numerical simulations of turbulence are essential to understanding vertical mixing processes and their dynamical controls on both physical and biogeochemical processes in coastal embayments. Using in situ data collected by bottom-mounted acoustic Doppler current profilers(ADCPs) and a free-falling microstructure profiler, as well as numerical simulations with a second-moment turbulence closure model, we studied turbulence and mixing in the Xiamen Bay, a freshwater-influenced tidal bay located at the west coast of the Taiwan Strait. Dynamically, the bay is driven predominantly by the M2 tide, and it is under a significant influence of the freshwater discharged from the Jiulong River. It is found that turbulence quantities such as the production and dissipation rates of the turbulent kinetic energy(TKE) were all subject to significant tidal variations, with a pronounced ebb-flood asymmetry. Turbulence was stronger during flood than ebb. During the flooding period, the whole water column was nearly well mixed with the depth-averaged TKE production rate and vertical eddy viscosity being up to 5?10?6 W kg?1 and 2?10?2 m2 s?1, respectively. In contrast, during the ebb strong turbulence was confined only to a 5?8 m thick bottom boundary layer, where turbulence intensity generally decreases with distance from the seafloor. Diagnosis of the potential energy anomaly showed that the ebb-flood asymmetry in turbulent dissipation and mixing was due mainly to tidal straining process as a result of the interaction between vertically shared tidal currents and horizontal density gradients. The role of vertical mixing in generating the asymmetry was secondary. A direct comparison of the modeled and observed turbulence quantities confirmed the applicability of the second-moment turbulence closure scheme in modeling turbulent processes in this weakly stratified tidally energetic environment, but also pointed out the necessity of further refinements of the model.     

Observation and numerical modeling of tidal dune dynamics

Tidal sand dune dynamics is observed for two tidal cycles in the Arcachon tidal inlet, southwest France. An array of instruments is deployed to measure bathymetric and current variations along dune profiles. Based on the measurements, dune crest horizontal and vertical displacements are quantified and show important dynamics in phase with tidal currents. We observed superimposed ripples on the dune stoss side and front, migrating and changing polarity as tidal currents reverse. A 2D RANS numerical model is used to simulate the morphodynamic evolution of a flat non-cohesive sand bed submitted to a tidal current. The model reproduces the bed evolution until a field of sand bedforms is obtained that are comparable with observed superimposed ripples in terms of geometrical dimensions and dynamics. The model is then applied to simulate the dynamics of a field of large sand dunes of similar size as the dunes observed in situ. In both cases, simulation results compare well with measurements qualitatively and quantitatively. This research allows for a better understanding of tidal sand dune and superimposed ripple morphodynamics and opens new perspectives for the use of numerical models to predict their evolution.     

Decadal evolution of tidal flats and channels in the Outer Weser estuary,Germany

This study focuses on the medium scale morphodynamics of the tidal flat and channel system Fedderwarder Priel, located in the Outer Weser estuary (Wadden Sea, Germany). Tidal channels and adjacent flats are highly dynamic systems whose morphologic evolution are driven by tidal, wind, and wave forcings. These coastal environments are an important ecosystem and react to changes in hydrodynamic conditions in various spatial and temporal scales. Based on annual medium-resolution digital elevation models from 1998 to 2016, we describe changes in the surface area over depth with hypsometries and use vertical dynamic trends in order to analyze and visualize the morphologic evolution of the Fedderwarder Priel and adjacent tidal channels. It is shown that several intertidal flats rise in the order of 1.3 to 5.6 cm/year. The findings indicate that the Outer Weser estuary was not in an equilibrium state for the investigated period, and tidal flats accreted with a rate exceeding mean sea level rise.     

The effect of asymmetric dune roughness on tidal asymmetry in the Weser estuary

The bed of estuaries is often characterized by ripples and dunes of varying size. Whereas smaller bedforms adapt their morphological shape to the oscillating tidal currents, large compound dunes (here: asymmetric tidal dunes) remain stable for periods longer than a tidal cycle. Bedforms constitute a form roughness, that is, hydraulic flow resistance, which has a large-scale effect on tidal asymmetry and, hence, on hydrodynamics, sediment transport, and morphodynamics of estuaries and coastal seas. Flow separation behind the dune crest and recirculation on the steep downstream side result in turbulence and energy loss. Since the energy dissipation can be related to the dune lee slope angle, asymmetric dune shapes induce variable flow resistance during ebb and flood phases. Here, a noncalibrated numerical model has been applied to analyze the large-scale effect of symmetric and asymmetric dune shapes on estuarine tidal asymmetry evaluated by residual bed load sediment transport at the Weser estuary, Germany. Scenario simulations were performed with parameterized bed roughness of symmetric and asymmetric dune shapes and without dune roughness. The spatiotemporal interaction of distinct dune shapes with the main drivers of estuarine sediment and morphodynamics, that is, river discharge and tidal energy, is shown to be complex but substantial. The contrasting effects of flood- and ebb-oriented asymmetric dunes on residual bed load transport rates and directions are estimated to be of a similar importance as the controls of seasonal changes of discharge on these net sediment fluxes at the Lower Weser estuary. This corroborates the need to consider dune-induced directional bed roughness in numerical models of estuarine and tidal environments.     

Fluctuations in the tidal limit of the Yangtze River estuary in the last decade

The tidal limit is the key interface indicating whether water levels will be affected by tidal waves, which is of great significance to navigation safety and regional flood control. Due to limitations in research methods, recent changes in the Yangtze River tidal limit, caused by sea level rise and large-scale engineering projects, urgently need to be studied. In this study,spectrum analysis was undertaken on measured water level data from downstream Yangtze River hydrological stations from2007 to 2016. The bounds of the tidal limit were identified through comparisons between the spectra and red noise curves, and the fluctuation range and characteristics were summarized. The results showed that:(1) During the extremely dry period, when the flow rate at Jiujiang station was about 8440 m3 s-1, the tidal limit was near Jiujiang; whereas during the flood season, when the flow rate at Jiujiang station was about 66700 m3 s-1, the tidal limit was between Zongyang Sluice and Chikou station.(2)From the upper to lower reach, the effect of the Jiujiang flow rate on the tidal limit weakens, while the effect of the Nanjing tidal range increases. The tidal limit fluctuates under similar flow rates and tidal ranges, and the fluctuation range increases with increasing flow rate and decreasing tidal range.(3) With the continued influence of rising sea levels and construction in river basin estuaries, the tidal limit may move further upstream.     

Statistical analyses and modeling approaches to hydrodynamic characteristics in alluvial aquifer

This study is aimed at quantifying the difference in aquifer's response to recharge between some different locations in a fan aquifer and a delta aquifer for a preliminary study of revealing mechanisms of water transport in alluvial aquifer. The aquifer's response to recharge is statistically quantified with the two viewpoints: (1) timing and volume of recharge and (2) time length of aquifer's holding water. For the first point, a statistical model that links precipitation and groundwater level is introduced, and its parameters are identified using correlation analysis. Our results show that the recharge rate at the toe is higher than that at the apex and at the delta. For the second point, the concept of ‘memory effect’ of aquifer is adopted and quantified using the autocorrelation and spectral analyses. Our results show that the memory effect is longer at the toe of fan than at the apex, and thus, a temporary increase of water level has about five times as long‐term influence on subsequent water levels at the toe of the fan as at the apex. This study demonstrates that the statistical analyses and modeling of hydrological data are useful for characterizing aquifer's hydrodynamics. Copyright © 2012 John Wiley & Sons, Ltd.     

Two-phase hydrodynamic and sediment transport modeling of wave-generated sheet flow

This numerical investigation was carried out to advance mechanistic understanding of sediment transport under sheet flow conditions. An Euler–Euler coupled two-phase flow model was developed to simulate fluid–sediment oscillatory sheet flow. Since the concentration of sediment particles is high in such flows, the kinematics of the fluid and sediment phases are strongly coupled. This model includes interaction forces, intergranular stresses and turbulent stress closure. Each phase was modeled via the Reynolds-Averaged Navier–Stokes equations, with interphase momentum conservation accounting for the interaction between the phases. The generation and transformation of turbulence was modeled using the two-equation k–ε$k''–''\epsilon$ turbulence model. Concentration and sediment flux profiles were compared with experimental data for sheet flow conditions considering both symmetric and asymmetric oscillatory flows. Sediment and fluid velocity variations, concentration profiles, sediment flux and turbulence parameters of wave-generated sheet flow were studied numerically with a focus on sediment transport characteristics. In all applications, the model predictions compared well with the experimental data. Unlike previous investigations in which the flow is driven by a horizontal pressure gradient, the present model solves the Navier–Stokes equations under propagating waves. The model’s ability to predict sediment transport under oscillatory sheet flow conditions underscores its potential for understanding the evolution of beach morphology.     

Fine sediment transport by tidal asymmetry in the high-concentrated Ems River: indications for a regime shift in response to channel deepening

This paper describes an analysis of the observed up-river transport of fine sediments in the Ems River, Germany/Netherlands, using a 1DV POINT MODEL, accounting for turbulence-induced flocculation and sediment-induced buoyancy destruction. From this analysis, it is inferred that the net up-river transport is mainly due to an asymmetry in vertical mixing, often referred to as internal tidal asymmetry. It is argued that the large stratification observed during ebb should be attributed to a profound interaction between turbulence-induced flocculation and sediment-induced buoyancy destruction, as a result of which the river became an efficient trap for fine suspended sediment. Moreover, an asymmetry in flocculation processes was found, such that during flood relative large flocs are transported at relative large flow velocity high in the water column, whereas during ebb, the larger flocs are transported at smaller velocities close to the bed??this asymmetry contributes to the large trapping mentioned above. The internal tidal asymmetry and asymmetry in flocculation processes are both driven by the pronounced asymmetry in flow velocities, with flood velocities almost twice the ebb values. It is further argued that this efficient trapping is the result of a continuous deepening of the river, and occurs when concentrations in the river become typically a few hundred mg/l; this was the case during the 1990 survey analyzed in this paper. We also speculate that a second regime shift did occur in the river when fluid mud layers become so thick that net transport rates are directly related to the asymmetry in flow velocity itself, probably still in conjunction with internal asymmetry as well. This would yield an efficient mechanism to transport large amounts of fine sediment far up-river, as currently observed.