Effective Longitudinal Dispersion Research Articles

Calculation of Chlorine and Fluoride Diffusion/Dispersion Coefficients in Water Supply Pipes

The present study was designed to report experimental results for effective diffusion and longitudinal dispersion coefficients for both chlorine and fluoride added to water supply pipes. This is the first extensive study to find effective diffusion and dispersion coefficients for chlorine and fluoride added to water supply pipes. A pipe of field diameter and length was used to find dispersion coefficients for both chlorine and fluoride. The study covers the Reynolds number (R) from laminar and transitional flows (R=200–20,000), for which more correct information about the dispersion coefficients might be needed in the literature. The paper first presents the solution of a two-dimensional (2D) advection-diffusion-reaction equation of concentration movement in a pipe. Second, it presents the analytical solution of a one-dimensional (1D) advection-dispersion-reaction equation. Third, by using the numerical and analytical solutions, effective diffusion and longitudinal dispersion coefficients of chlorine and fluoride are calculated. At the end, simulation of part of a real network with the results reported in this study is presented.

Upstream swimming and Taylor dispersion of active Brownian particles

Locomotion of self-propelled particles such as motile bacteria or phoretic swimmers often takes place in the presence of applied flows and confining boundaries. Interactions of these active swimmers with the flow environment are important for the understanding of many biological processes, including infection by motile bacteria and the formation of biofilms. Recent experimental and theoretical works have shown that active particles in a Poiseuille flow exhibit interesting dynamics including accumulation at the wall and upstream swimming. Compared to the well-studied Taylor dispersion of passive Brownian particles, a theoretical understanding of the transport of active Brownian particles (ABPs) in a pressure-driven flow is relatively less developed. In this paper, employing a small wave-number expansion of the Smoluchowski equation describing the particle distribution, we explicitly derive an effective advection-diffusion equation for the cross-sectional average of the particle number density in Fourier space. We characterize the average drift (specifically upstream swimming) and effective longitudinal dispersion coefficient of active particles in relation to the flow speed, the intrinsic swimming speed of the active particles, their Brownian diffusion, and the degree of confinement. In contrast to passive Brownian particles, both the average drift and the longitudinal dispersivity of ABPs exhibit a nonmonotonic variation as a function of the flow speed. In particular, the dispersion of ABPs includes the classical shear-enhanced (Taylor) dispersion and an active contribution called the swim diffusivity. In the absence of translational diffusion, the classical Taylor dispersion is absent and we observe a giant longitudinal dispersion in the strong flow limit. Our continuum theory is corroborated by a direct Brownian dynamics simulation of the Langevin equations governing the motion of each ABP.

Random walk particle tracking simulation on scalar diffusion with irreversible first-order absorption boundaries.

Scalar transport in an open channel with irreversible first-order absorption boundaries is investigated through random walk particle tracking method (RWPT). We provide the pre-asymptotic behavior of scalar transport as well as a comparison with existing asymptotic formulations. The RWPT method is based on the development of the probability that a particle is absorbed at the boundary. The random walk model is applied to simulate three main parameters of the scalar transport in the laminar flow with sorption boundaries. The three parameters are (1) the attenuation coefficient reflecting the depletion rate of the total mass; (2) the effective velocity coefficient representing the effective advection of the total mass; and (3) the effective longitudinal dispersion coefficient expressing the effective spreading rate of the solute plume. When the Peclet number is large, numerical results agree well with previous theoretical solutions. Results show that (1) the three coefficients reach their asymptotic values when the dimensionless time τ (= Dt/H2) is about 0.5, where t is the time, H is the water depth, and D is the molecular diffusion coefficient and (2) when the Damkohler number reaches a value around 100, the boundary can be treated as a totally absorptive boundary.

Regional-scale analysis of karst underground flow deduced from tracing experiments: examples from carbonate aquifers in Malaga province, southern Spain

Tracer concentration data from field experiments conducted in several carbonate aquifers (Malaga province, southern Spain) were analyzed following a dual approach based on the graphical evaluation method (GEM) and solute transport modeling to decipher flow mechanisms in karst systems at regional scale. The results show that conduit system geometry and flow conditions are the principal factors influencing tracer migration through the examined karst flow routes. Solute transport is mainly controlled by longitudinal advection and dispersion throughout the conduit length, but also by flow partitioning between mobile and immobile fluid phases, while the matrix diffusion process appears to be less relevant. The simulation of tracer breakthrough curves (BTCs) suggests that diffuse and concentrated flow through the unsaturated zone can have equivalent transport properties under extreme recharge, with high flow velocities and efficient mixing due to the high hydraulic gradients generated. Tracer mobilization within the saturated zone under low flow conditions mainly depends on the hydrodynamics (rather than on the karst conduit development), which promote a lower longitudinal advection and retardation in the tracer migration, resulting in a marked tailing effect of BTCs. The analytical advection-dispersion equation better approximates the effective flow velocity and longitudinal dispersion estimations provided by the GEM, while the non-equilibrium transport model achieves a better adjustment of most asymmetric and long-tailed BTCs. The assessment of karst underground flow properties from tracing tests at regional scale can aid design of groundwater management and protection strategies, particularly in large hydrogeological systems (i.e. transboundary carbonate aquifers) and/or in poorly investigated ones.

Pore-Network Simulation of Unstable Miscible Displacements in Porous Media

We develop a comprehensive pore-network model for simulation of miscible displacements in laboratory-scale porous media. The goal is to study the effect of pore-scale heterogeneity, and in particular pore-size distribution and pore connectivity, as well as mobility ratio M on miscible displacements, and compute the effective longitudinal dispersion coefficient as a function of M. At the pore level the phenomenon is described by the convective-diffusion equation with a Taylor–Aris dispersion coefficient. Extensive simulations indicate the strong effect of the connectivity and broadness of the pore-size distribution on the efficiency of the displacement process. Moreover, the model enables us to compute the effective dispersion coefficient as a function of the mobility ratio of the two fluids that, to our knowledge, has not been calculated before. Simple generalization of the model will make it possible to study more complex related phenomena, and in particular injection of $$\mathrm{CO}_2$$ into a pore space that contains saline water.

A Subtidal Model of Temperature for a Well-Mixed Narrow Estuary: the Guadalquivir River Estuary (SW Spain)

This paper analyzes thermal energy transport in the narrow and tidally energetic Guadalquivir River Estuary (SW Spain). Measurements from a comprehensive monitoring campaign (2008–2011) reveal the forcing factors of the temperature field and its spatio-temporal variability. The along-channel thermal energy gradient reaches magnitudes of ∼375 J/m4 near the mouth during the summer and winter. The water temperature is primarily controlled by shortwave radiation, latent heat transfer through the free surface, and tidal advection, whereas it depends less on freshwater discharge and longitudinal dispersion. The tidally averaged effective longitudinal thermal dispersion coefficient was evaluated at several stretches for each tidal cycle. The mean values of the coefficient tend to increase landward and are on the order of ∼103, larger than (but of the same order of magnitude as) the salinity coefficient values. Based on these analyses, a deterministic operational model for thermal energy transport was developed. The model solves the tidally and cross-sectionally averaged advection–dispersion equation for the thermal energy balance and obtains accurate fits of the subtidal temperature field at any location within the estuary. The modeled water temperatures agreed well with the observations at all the stations (coefficients of determination, R 2 greater than 0.98), even after the seasonal oscillation in radiation was removed (R 2 > 0.77).

A field-scale long-term study on radionuclide transport through weathered granites at a site in southern China

A regional strategy for safety disposal of low- and intermediate-level radioactive wastes (LILW) has been implemented in China to protect humans and the environment. A joint onsite and laboratory investigation was conducted for a field site in southern China to assess the probability for safe disposal of LILW, which requires the understanding of long-term radionuclide transport behaviors under field conditions. This study presents the field-scale modeling of radionuclide transport through weathered granites for a conservative, a weakly sorbing and a strongly sorbing tracer by incorporating laboratory and field data. The field-scale radionuclide transport model was developed on the basis of a validated long-term groundwater flow model and field-measured dispersion coefficient, as well as laboratory-characterized strontium and cesium distribution coefficients in the weathered granites. The model was then used to perform the long-term transport prediction and risk assessment of radionuclide pollution for both the natural site setting and the graded site setting. Model simulation reveals that the numerical sensitivities of calculated concentrations are tracer dependent and changing with time. The conservative radionuclide is most sensitive to changes in hydraulic conductivity (K) while slightly sensitive to changes in effective porosity (φ), specific yield (μ) and longitudinal dispersion coefficient (D L), indicating advection is the main transport process of conservative radionuclide. The weakly and strongly sorbing tracers, on the other hand, are most sensitive to changes in the distribution coefficient (K d) and less sensitive to changes in the rest of model parameters, revealing that sorption is the main process for controlling the transport of sorbing tracers. A conservative radionuclide plume moves at an average velocity of about 54 m/a, which is too fast to be considered as safe disposal under the natural site setting. However, the plume of the conservative radionuclide could be slowed down to a velocity around 5.3 m/a due to the reduction of the hydraulic gradient under the graded site setting. Therefore, the conservative radioactive wastes could be disposed at the mid-eastern part of the site under the graded site setting because the transport path has been prolonged and thus no conservative radionuclides could migrate out of the site in a reasonable timeframe. For the sorbing tracers, however, results of the computed transport distance are 40 and 2 m at 500 years, respectively. Therefore, they can be disposed safely at the site under both natural and graded site settings. This study provides an insight to the field-scale long-term behaviors of radionuclide transport. The integrated modeling method presented in this study is most useful for the environmental impact assessment of the site conditions relevant to the safe disposal of hazardous wastes.

Flow and axial dispersion in a sinusoidal-walled tube: Effects of inertial and unsteady flows

In this work, we consider a sinusoidal-walled tube (a three-dimensional tube with sinusoidally-varying diameter) as a simplified conceptualization of flow in porous media. Direct numerical simulation using computational fluid dynamics (CFD) methods was used to compute velocity fields by solving the Navier–Stokes equations, and also to numerically solve the volume averaging closure problem, for a range of Reynolds numbers (Re) spanning the low-Re to inertial flow regimes, including one simulation at Re=449 for which unsteady flow was observed. The longitudinal dispersion observed for the flow was computed using a random walk particle tracking method, and this was compared to the longitudinal dispersion predicted from a volume-averaged macroscopic mass balance using the method of volume averaging; the results of the two methods were consistent. Our results are compared to experimental measurements of dispersion in porous media and to previous theoretical results for both the low-Re, Stokes flow regime and for values of Re representing the steady inertial regime. In the steady inertial regime, a power-law increase in the effective longitudinal dispersion (DL) with Re was found, and this is consistent with previous results. This rapid rate of increase is caused by trapping of solute in expansions due to flow separation (eddies). One unsteady (but non-turbulent) flow case (Re=449) was also examined. For this case, the rate of increase of DL with Re was smaller than that observed at lower Re. Velocity fluctuations in this regime lead to increased rates of solute mass transfer between the core flow and separated flow regions, thus diminishing the amount of tailing caused by solute trapping in eddies and thereby reducing longitudinal dispersion. The observed tailing was further explored through analysis of concentration skewness (third moment) and its assymptotic convergence to conventional advection–dispersion behavior (skewness=0). The method of volume averaging was applied to develop a skewness model, and demonstrated that the skewness decreases as a function of inverse square root of time. Our particle tracking simulation results were shown to conform to this theoretical result in most of the cases considered.

Spatio‐temporal distribution, along‐channel transport, and post‐riverflood recovery of salinity in the Guadalquivir estuary (SW Spain)

This paper presents an experimental analysis of the salinity distribution, the salt balance, and the variation of the saline intrusion in comparison to the freshwater discharge in the Guadalquivir estuary, which is a mesotidal system regulated and normally subjected to extremely low river flows. In such low‐flow conditions, it is positive, well‐mixed, and tidally dominated. The estuary is also characterized by a nonstationary, effective longitudinal dispersion coefficient, whose probability density becomes increasingly narrower and whose mean value is higher further upstream. The tidal‐averaged salt flux is controlled by the following mechanisms (in order of importance): the nontidal transport, the Stokes transport, and the tidal pumping induced by the covariance between the current and salinity. These three factors account for more than 98% of the flux variation. In high river‐flow conditions, the subtidal response and recovery of the estuary to changes in the river flow is analyzed. The increase in the tidal‐averaged salinity during the first 2 weeks of the post‐riverflood recovery in the middle and upper sections of the estuary is found to be linear in time. During that time, the celerity of the salt intrusion front was 4cm/s. The 2 psu isohaline salt intrusion X2 exhibits a complex dependence on the river flow Qd, including the effects of human interventions in the estuary. Three regimes are identified for the intrusion: X2=57.0± 2.1km for discharges of less than 20m 3/s, X2 proportional to between 20 and 1000m 3/s, and X2proportional to for larger discharges.

Spatial variation of the longitudinal dispersion coefficient in an estuary

Abstract. The effective longitudinal dispersion is a primary tool for determining property distributions in estuaries. Most previous studies have examined the longitudinal dispersion coefficient for the average tidal condition. However, information on spatial and temporal variations of this coefficient at low and high tides is scarce. Three years of hydrographic data taken at low and high tide along the main axis of the Sumjin River Estuary (SRE), Korea are used to estimate the spatial and temporal variation of the effective longitudinal dispersion coefficient. The range of the dispersion coefficient is rather broad at high water slack (HWS) and narrower at low water slack (LWS) because of the different tidal amplitudes. The spatially varying dispersion coefficient has maximal values (>300 m2 s−1) near the mouth at high water and decreases gradually upstream, with fluctuations. The temporally varying dispersion coefficient appears to be positively correlated with river discharges at both low and high tide. The dispersion varies with the square root of river discharges at HWS and LWS. The dispersive salt fluxes increases with increasing river discharges and decreases with decreasing river discharges at HWS and LWS. Estimation of the numerical values of the effective longitudinal dispersion coefficient in the SRE can be useful for better understanding of the distributions of other tracers in the SRE as well as for developing and testing hypotheses about various mixing mechanisms.

Hydrodynamic dispersion in steady buoyancy‐driven geological flows

An analytical model is developed to evaluate mixing induced by natural convection in a fluid‐saturated porous medium. First, the velocity and concentration fields are decoupled to generate a steady state velocity field and initiate a naturally convective system. In order to decouple the velocity and concentration fields, a steady thermal natural convection is established by imposing a destabilizing vertical temperature gradient across a porous layer and then introducing a passive tracer into the system. Based on the steady velocity field, effective longitudinal and transverse dispersion coefficients are evaluated using the shear flow dispersion theory, and convective mixing of the passive tracer is obtained using the developed analytical mixing model. The estimated dispersion coefficients and convective mixing are then characterized by the system Rayleigh and Sherwood numbers. The mixing obtained by the analytical model is then compared with high‐resolution numerical simulations. The results reveal that the simple analytical solution represents the nonlinear mixing involved in such a system and agrees with the numerical results. The developed model has potential applications in geophysical and geothermal buoyancy‐driven flows.

Stochastic modeling of reactive transport in wetlands

This study describes the development of a general model for reaction in and performance of spatially heterogeneous bioreactors such as treatment wetlands. The modeled domain possesses local-scale velocities, reaction rates and transverse dispersion coefficients that are functions of an underlying heterogeneity variate representing one or more controlling biophysical attributes, for example, reactive surface area (submerged plant) density. Reaction rate coefficients are treated as related to local velocities in an inverse square fashion via their mutual dependence upon the variate. The study focuses on the solution for the steady-state case with constant inlet concentration. Results compare well with exact solutions developed for laterally-bounded systems in which the heterogeneity is represented explicitly. Employing the bicontinuum analogue of a second-order model, an expression for an effective longitudinal dispersion coefficient as a function of travel distance is developed using the method of moments. The result provides insights into the behavior of concentration as transverse mixing drives the system asymptotically toward Fickian longitudinal dispersion. The model may represent an improvement over other approaches for characterizing treatment wetland performance because it accounts for evolving shear flow dispersion, and because parameters are few in number, physically based, and invariant with mean velocity.

Asymptotic dispersion in 2D heterogeneous porous media determined by parallel numerical simulations

We determine the asymptotic dispersion coefficients in 2D exponentially correlated lognormally distributed permeability fields by using parallel computing. Fluid flow is computed by solving the flow equation discretized on a regular grid and transport triggered by advection and diffusion is simulated by a particle tracker. To obtain a well‐defined asymptotic regime under ergodic conditions (initial plume size much larger than the correlation length of the permeability field), the characteristic dimension of the simulated computational domains was of the order of 103 correlation lengths with a resolution of ten cells by correlation length. We determine numerically the asymptotic effective longitudinal and transverse dispersion coefficients over 100 simulations for a broad range of heterogeneities σ2 ∈ [0, 9], where σ2 is the lognormal permeability variance. For purely advective transport, the asymptotic longitudinal dispersion coefficient depends linearly on σ2 for σ2 < 1 and quadratically on σ2 for σ2 > 1 and the asymptotic transverse dispersion coefficient is zero. Addition of homogeneous isotropic diffusion induces an increase of transverse dispersion and a decrease of longitudinal dispersion.

Estimating effective longitudinal dispersion in the Chesapeake Bay

An analysis of Environmental Protection Agency's Chesapeake Bay Program hydrographic dataset shows that the bay responds coherently to variability in freshwater flux. Mean salinity and salinity stratification both respond to variability in freshwater flux on time scales of roughly 90 days. Stratification is also influenced by local wind forcing but on much shorter (4–5 day) time scales. The volume of available data allows the effective longitudinal dispersion coefficient to be estimated as a function of either time or space. Values for this dispersion coefficient vary between 200 and 1000 m 2 s −1, with mean values around 650 m 2 s −1. The spatially dependent structure has a maximum roughly 75 km from the head of the estuary, and decreases gradually towards the mouth. The temporally varying effective dispersion varies spatially as the inverse of the estuarine cross-section, and temporally as the cube root of the freshwater flux, and is at least qualitatively consistent with models of estuarine circulation and results of previous field studies. Estimates of the numerical values of the dispersion are useful for better understanding distributions of other tracers within the bay, as well as providing another metric against which numerical models should be measured.

Enhanced dispersion in cylindrical packed beds

The effective longitudinal dispersion constant, D(L)(eff), in cylindrical packed beds is larger than in the bulk due to the existence of radial inhomogeneities induced by the cylinder walls. For dense random packed beds, D(L)(eff) can be several times larger than the bulk value, even for arbitrarily large cylinder radius, R. The time-scale for attaining asymptotic dispersion rates in a cylindrical geometry is neither the convective nor the diffusive time-scale, but rather D(T)/R(2), where D(T) is the bulk transverse dispersion rate. Similar effects are predicted for packed beds confined in ducts of any cross-sectional geometry. The case of a rectangular duct, compared with an infinite slit, provides an intuitive model for the influence of walls in the limit as R goes to infinity.

Solute transport in variable‐aperture fractures: An investigation of the relative importance of Taylor dispersion and macrodispersion

Dispersion of solutes in a variable aperture fracture results from a combination of molecular diffusion and velocity variations in both the plane of the fracture (macrodispersion) and across the fracture aperture (Taylor dispersion). We use a combination of physical experiments and computational simulations to test a theoretical model in which the effective longitudinal dispersion coefficient DL is expressed as a sum of the contributions of these three dispersive mechanisms. The combined influence of Taylor dispersion and macrodispersion results in a nonlinear dependence of DL on the Peclet number (Pe = V〈b〉/Dm, where V is the mean solute velocity, 〈b〉 is the mean aperture, and Dm is the molecular diffusion coefficient). Three distinct dispersion regimes become evident: For small Pe (Pe ≪1), molecular diffusion dominates resulting in DL ∝ Pe0; for intermediate Pe, macrodispersion dominates (DL ∝ Pe); and for large Pe, Taylor dispersion dominates (DL ∝ Pe2). The Pe range corresponding to these different regimes is controlled by the statistics of the aperture field. In particular, the upper limit of Pe corresponding to the macrodispersion regime increases as the macrodispersivity increases. Physical experiments in an analog, rough‐walled fracture confirm the nonlinear Pe dependence of DL predicted by the theoretical model. However, the theoretical model underestimates the magnitude of DL. Computational simulations, using a particle‐tracking algorithm that incorporates all three dispersive mechanisms, agree very closely with the theoretical model predictions. The close agreement between the theoretical model and computational simulations is largely because, in both cases, the Reynolds equation describes the flow field in the fracture. The discrepancy between theoretical model predictions and DL estimated from the physical experiments appears to be largely due to deviations from the local cubic law assumed by the Reynolds equation.

The influence of surface sorption on dispersion in parallel plate fractures

Solute transport in a parallel plate fracture is simulated using a random walk model which accounts explicitly for sorption onto the fracture walls. With the hypothesis that solute must move into the vicinity of the fracture wall in order to participate in the sorption process, three implications follow. (1) In comparison to a nonreactive solute, a sorbing solute requires a greater entrance length along the fracture before transverse homogenization is established. The increase in this entrance length is proportional to the surface retardation factor. (2) At short transport distances, surface sorption leads to a high degree of non-uniform retardation. The transport distance required to establish uniform retardation of the entire solute mass varies greatly with geometric conditions and sorption strength. (3) Surface sorption results in enhanced longitudinal spreading of the solute mass in transport regimes which favor advective transport along the fracture relative to transverse diffusion across the fracture aperture. At distances greater than that required for transverse homogenization, an effective longitudinal dispersion coefficient can be defined that describes this enhanced dispersion for a wide range of fluid velocities. The magnitude of enhanced dispersion increases with sorption strength.

Stochastic Theory for Irregular Stream Modeling. II: Solute Transport

A stochastic theory is developed for longitudinal dispersion in natural streams. Irregular variations in river width and bed elevation are conveniently represented as one-dimensional random fields. Longitudinal solute migration is described by a one-dimensional stochastic solute transport equation. When boundary variations are small and statistically homogeneous, the stochastic transport equation is solved in closed-form using a stochastic spectral technique. The results show that large scale longitudinal transport can be represented as a gradient dispersion process described by an effective longitudinal dispersion coefficient. The effective coefficient reflects longitudinal mixing due to flow variation both within the river cross section and along the flow and can be considerably greater than that of corresponding uniform channels. The discrepancy between uniform channels and natural rivers increases as the variances of river width and bed elevation increase, especially when the mean flow Froude number is high.

Dispersion of a passive solute in non-ergodic transport by steady velocity fields in heterogeneous formations

An inert solute is convected by a steady random velocity field, which is associated with flow through a heterogeneous porous formation. The log conductivity and the velocity are stationary random space functions. The log conductivity Y is assumed to be normal, with an isotropic two-point correlation of variance σY2 and of finite integral scale I. The solute cloud is of a finite input zone of lengthscale l. The transport is characterized with the aid of the spatial moments of the solute body. The effective dispersion coefficient is defined as half of the rate of change with time of the second spatial moment with respect to the centroid. Under the ergodic hypothesis, which is bound to be satisfied for l/I [Gt ] 1, the centroid moves with the mean velocity U and the longitudinal dispersion coefficient [dscr ]L tends to its constant, Fickian, limit. Under a Lagrangian first-order analysis in σY2 it has been found that [dscr ]L = σY2UI.This study addresses the computation of the effective longitudinal dispersion coefficient for a finite input zone, for which ergodic conditions may not be satisfied. In this case the centroid trajectory and the second spatial moments are random variables. In line with a previous work (Dagan 1990) the effective dispersion coefficient DL is defined as half the rate of change of the expected value of the second spatial moment for large transport time. The aim of the study is to derive DL and its dependence upon l/I and in particular to determine the conditions under which it tends to the ergodic limit [dscr ]L. The computation is carried out separately for a thin body aligned with the mean flow and one transverse to it. In the first case it is found that DL is equal to zero, i.e. the streamlined body does not disperse in the mean. This result is explained by the correlation between the trajectories of the leading and trailing edges, respectively, once the latter reaches the position of the first. The relatively modest increase of the mean second spatial moment is effectively computed. In the case of a thin body initially transverse to the mean flow, DL may reach the ergodic limit [dscr ]L for a ratio l/I of the order 102. For smaller values, DL is found to be bounded from above, and its maximum depends on l but not on I. The uncertainty caused by the randomness of the velocity field is manifested in the trajectory of the centroid rather than in the effective dispersion.

Characterization of a Malaysian Mangrove Estuary

The Sungai Merbok estuary, in wet tropical Peninsular Malaysia, borders the Straits of Malacca. Tide, current, and salinity data are used to describe the salient hydrographic features of the mangrove-fringed system. The Sungai Merbok estuary is characterized by a 1.7 m semidiurnal tide with a 0.16 form number, peak currents of 1.3 m s−1, and mean freshwater discharge of 20 m3 s−1. The system is classified as 2a/2b estuary (Hansen and Rattray 1966) or 1a/1b during periods of low runoff. Gravitational circulation is highly variable (but coincides with the neap stratification) and vertical stratification varies from 10−2 to 1. The estuary displays a pronounced fortnightly neap-spring stratification-destratification cycle. The effective longitudinal dispersion coefficient is approximately 100 m2 s−1.