An idealized model study of flocculation on sediment trapping in an estuarine turbidity maximum

Flocculation has an important impact on particle trapping in estuarine turbidity maximum (ETM) through associated increases in particle settling velocity. To quantify the importance of the flocculation processes, a size-resolved flocculation model is implemented into an ocean circulation model to simulate fine-grained particle trapping in an ETM. The model resolves the particle size from robust small flocs, about 30 μm, to very large flocs, over 1000 μm. An idealized two-dimensional model study is performed to simulate along-channel variations of suspended sediment concentrations driven by gravitational circulation and tidal currents. The results indicate that the flocculation processes play a key role in generating strong tidal asymmetrical variations in suspended sediment concentration and particle trapping. Comparison with observations suggests that the flocculation model produces realistic characteristics of an ETM.     

Three-dimensional numerical modelling of wind-driven circulation in a homogeneous lake

A three-dimensional numerical model has been developed to study wind-induced circulation patterns in a shallow homogeneous lake with a complex bathymetry. The governing equations are the unsteady Reynolds-Averaged Navier–Stokes equations in which the non-hydrostatic pressure distribution has been included. The model was tested against analytical solutions and laboratory data for wind-induced currents and then applied to Esthwaite Water, a small lake in Cumbria, UK. The model was used to study the main model parameters and to generate typical circulation patterns for a variety of conditions in the lake. Simulations showed that a non-hydrostatic pressure distribution did not have any noticeable influence on the overall circulation pattern in the lake. However, comparisons with field data at some measurement stations in the near-shore region with sharply varying bottom topography showed that the hydrodynamic pressure component had some influence on the vertical velocity profile.     

The influence of asymmetries in overlying stratification on near-bed turbulence and sediment suspension in a partially mixed estuary

Data collected from the York River estuary demonstrate the importance of asymmetries in stratification to the suspension and transport of fine sediment. Observations collected during two 24-h deployments reveal greater concentrations of total suspended solids during the flood phase of the tide despite nearly symmetric near-bed tidal current magnitude. In both cases, tidally averaged net up-estuary sediment transport near the bed was clearly observed despite the fact that tidally averaged residual near-bed currents were near zero. Tidal straining of the along-channel salinity gradient resulted in a stronger pycnocline lower in the water column during the ebb phase of the tide and appeared to limit sediment suspension. Indirect measurements suggest that the lower, more intense, pycnocline on the ebb acted as a barrier, limiting turbulent length scales and reducing eddy diffusivity well below the pycnocline, even though the lower water column was locally well mixed. In order to more conclusively link changes in stratification to properties of near-bed eddy viscosity and diffusivity, longer duration tripod and mooring data from an additional experiment are examined, that included direct measurement of turbulent velocities. These additional data demonstrate how slight increases in stratification can limit vertical mixing near the bed and impact the structure of the eddy viscosity below the pycnocline. We present evidence that the overlying pycnocline can remotely constrain the vertical turbulent length scale of the underlying flow, limiting sediment resuspension. As a result, the relatively small changes in stratification caused by tidal straining of the pycnocline allow sediment to be resuspended higher in the water column during the flood phase of the tide, resulting in preferential up-estuary transport of sediment.Responsible Editor: Iris Grabemann     

Near-bed sediment transport in a heavily modified coastal plain estuary

Numerous estuaries of the world have been strongly modified by human activities.These interferences can make great adjustments of not only sediment transport processes,but also the collective behavior of the estuary.This paper provides a typical case of a heavily modified coastal plain estuary of Sheyang on the China coast,where a sluice barrage was built in 1956 to stop the intrusions of storm surges and saline water.Four sets of instrumented tripods were simultaneously deployed along a cross-shore transect to continuously observe near-bed flow currents and sediment transport.The in-situ surveys lasted over a spring and neap tide cycle when a strong wind event occurred in the neap tide.Comparisons of flows and sediment transport between tide-dominated and wind-dominated conditions demonstrated the important role of episodic wind events in flows and sediment transport.The wind-induced currents,bottom stresses,and sediment transport rates were significantly greater when wind was present than corresponding quantities induced by the tides.The long-shore sediment transport induced by winds exceeds the cross-shore component,especially near the river mouth bar.These results indicate the noticeable importance of wave-dominated coastal processes in shaping topographic features.A regime shift of estuarine evolution under highly intense human forcing occurs from fluvial to marine processes.This finding suggests that the management strategy of the estuarine system should focus on the restoration of estuarine processes,rather than the present focus on inhibition of marine dynamics.     

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.     

Influence of thermal stratification on the surfacing and clustering of floaters in free surface turbulence

The dispersion of floaters, small organic particles lighter than water, on the free surface of an open turbulent channel flow subject to thermal stratification is studied by Direct Numerical Simulation (DNS) of turbulence and Lagrangian Particle Tracking (LPT). Constant heat flux is maintained at the free surface of the channel, the bottom wall is adiabatic and the turbulent flow is driven by a pressure gradient. This archetypal flow setup mimics an environmentally plausible situation which can be found in terrestrial water bodies. The free surface turbulence characteristic of such flows has a strong influence on the distribution of the floaters: the objective of this work is to study the effect of different regimes of stable stratification on the surface distribution of floaters. The distribution of the floaters can possibly influence the transfer of chemical species across the water/atmosphere interface. Our results show that the modification of turbulence due to the thermal stratification strongly influences the settling velocity of floaters in the bulk of the flow. At the surface, stratification effects are also observed on the clustering of the floaters: the filamentary patterns of floaters observed in unstratified turbulence are progressively lost as thermal stratification increases, and the distribution of the floaters remains roughly two-dimensional.     

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

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

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.     

Bowhead whales feed on plankton concentrated by estuarine and tidal currents in Academy Bay,Sea of Okhotsk

Academy Bay in the Sea of Okhotsk is an important summertime feeding ground for pelagic-feeding Bowhead whales (Balaena mysticetus) in the western subarctic North Pacific. The present work combines satellite observations with physical (CTD, currents, tides) and biological (zooplankton sampling) measurements. Data obtained aboard the RV Lugovoe in August–September 2003 and July–August 2004 show dense populations of zooplankton (such as copepods Calanus glacialis, Pseudocalanus sp., pteropods Limacina helicina, and chaetognaths Sagitta sp.) that are concentrated by physical mechanisms within this critically important ecosystem. We show that near-bottom, cold water intrusions from the northern Sea of Okhotsk associated with the estuarine circulation advect arctic herbivorous calanoids (Calanus glacialis) and mollusks Limacina helicina into the region.     

Long term trends of Hg uptake in resident fish from a polluted estuary

Mercury contamination of fish is dependent upon a system’s ability to transform inorganic Hg into biologically available forms; however, fish biometrics also play an important role. To assess long term trends in Hg concentrations in sand flathead (Platycephalus bassensis) a polynomial model, corrected for fish length, was used to evaluate temporal trends and spatial variability, while growth rates were estimated using the Von Bertalanffy length-at-age model. Hg concentrations showed no decrease over time, and generally remained near recommended consumption levels (0.5 mg kg⁻¹). Previously reported spatial differences in Hg concentrations were not supported by the data once the models were corrected for fish length. Growth rate variation accounted for a large part of the previously published spatial differences. These results suggest that inclusion of fish biometrics is necessary to facilitate an accurate interpretation of spatial and temporal trends of contaminant concentrations in long term estuarine and marine monitoring programs.     

Numerical modelling of channel migration with application to laboratory rivers

The paper presents the development of a morphological model and its application to experimental model rivers.The model takes into account the key processes of channel migration,including bed deformation,bank failure and wetting and drying.Secondary flows in bends play an important role in lateral sediment transport,which further affects channel migration.A new formula has been derived to predict the near-bed secondary flow speed,in which the magnitude of the speed is linked to the lateral water level gradient.Since only non-cohesive sediment is considered in the current study,the bank failure is modelled based on the concept of submerged angle of repose.The wetting and drying process is modelled using an existing method.Comparisons between the numerical model predictions and experimental observations for various discharges have been made.It is found that the model predicted channel planform and cross-sectional shapes agree generally well with the laboratory observations.A scenario analysis is also carried out to investigate the impact of secondary flow on the channel migration process.It shows that if the effect of secondary flow is ignored,the channel size in the lateral direction will be seriously underestimated.     

Empirical observations and numerical modelling of tides,channel morphology,and vegetative effects on accretion in a restored tidal marsh

Tidal marshes form at the confluence between estuarine and marine environments where tidal movement regulates their developmental processes. Here, we investigate how the interplay between tides, channel morphology, and vegetation affect sediment dynamics in a low energy tidal marsh at the Paul S. Sarbanes Ecosystem Restoration Project at Poplar Island. Poplar Island is an active restoration site where fine-grained material dredged from navigation channels in the upper Chesapeake Bay are being used to restore remote tidal marsh habitat toward the middle bay (Maryland, USA). Tidal currents were measured over multiple tidal cycles in the inlets and tidal creeks of one marsh at Poplar Island, Cell 1B, using Acoustic Doppler Current Profilers (ADCP) to estimate water fluxes throughout the marsh complex. Sediment fluxes were estimated using acoustic backscatter recorded by ADCPs and validated against total suspended solid measurements taken on site. A high-resolution geomorphic survey was conducted to capture channel cross sections and tidal marsh morphology. We integrated simple numerical models built in Delft3d with empirical observations to identify which eco-geomorphological factors influence sediment distribution in various channel configurations with differing vegetative characteristics. Channel morphology influences flood-ebb dominance in marshes, where deep, narrow channels promote high tidal velocities and incision, increasing sediment suspension and reducing resilience in marshes at Poplar Island. Our numerical models suggest that accurately modelling plant phenology is vital for estimating sediment accretion rates. In-situ observations indicate that Poplar Island marshes are experiencing erosion typical for many Chesapeake Bay islands. Peak periods of sediment suspension frequently coincide with the largest outflows of water during ebb tides resulting in large sediment deficits. Ebb dominance (net sediment export) in tidal marshes is likely amplified by sea-level rise and may lower marsh resilience. We couple field observations with numerical models to understand how tidal marsh morphodynamics contribute to marsh resilience. © 2019 John Wiley & Sons, Ltd.     

Non-linear channel–shoal dynamics in long tidal embayments

The dynamics of finite-amplitude bed forms in a tidal channel is studied with the use of an idealized morphodynamic model. The latter is based on depth-averaged equations for the tidal flow over a sandy bottom. The model considers phenomena on spatial scales of the order of the tidal excursion length. Transport of sediment mainly takes place as suspended load. The reference state of this model is characterized by a spatially uniform M₂ tidal current over a fixed horizontal bed. The temporal evolution of deviations from this reference state is governed by amplitude equations: these are a set of non-linear equations that describe the temporal evolution of bed forms. These equations are used to obtain new morphodynamic equilibria which may be either static or time-periodic. Several of these bottom profiles show strong similarity with the tidal bars that are observed in natural estuaries. The dependence of the equilibrium solutions on the value of bottom friction and channel width is investigated systematically. For narrow channels (width small compared to the tidal excursion length) stable static equilibria exist if bottom friction is slightly larger than r_cr. For channel widths more comparable to the tidal excursion length, multiple stable steady states may exist for bottom friction parameter values below r_cr. Regardless of channel width, stable time-periodic equilibria seem to emerge as the bottom friction is increased.Responsible Editor: Jens Kappenberg     

Short-term dispersion and turbulence in a complex-shaped estuarine embayment

Initial dispersion of material in complex-shaped embayments is examined using observations and scaling based in Crail Bay, Pelorus Sound, New Zealand in autumn of 2005. These observations show the highly variable nature of dispersive transport in an embayment with multiple headlands. Acoustic current profiler-derived typical flow speeds were around 0.05 m s⁻¹ which resulted in drifter-derived short-term (<6 h) horizontal eddy diffusivities of the order of 1 m² s⁻¹ which is somewhat larger than the empirical paradigm. Microstructure estimates of the turbulent kinetic energy dissipation rate were in the range 10⁻⁹–10⁻⁷ m² s⁻³, with some evidence that sidewalls influence the variability and that headlands increase the dissipation rate by at least and order of magnitude. A new parameter relating horizontal diffusivity to circulation and the tide is proposed. This and other scaling comparisons indicate that the headlands in Crail Bay create similar effects to those studied in other systems. However, the long decay times estimated for eddies here implies that they likely interact with other headlands, unlike some previously studied examples.     

Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction

A finite-element method for computing the electric field in a 3-D conductivity model of the Earth for plane wave sources, thus enabling magnetotelluric responses to be calculated, is presented. The method incorporates in the iterative solution of the electric-field system of equations the divergence correction technique introduced for finite-difference solutions by Smith (1996). The correction technique accelerates the development of the discontinuity of the normal component of the approximate electric field across conductivity discontinuities. The convergence rate of the iterative solution is improved significantly, especially for low frequencies. The correction technique involves computing the divergence of the current density for the approximate electric field, computing the static potential whose source is this divergence of the current density, and ‘correcting’ the approximate electric field by subtracting from it the gradient of the potential. This is repeated at regular intervals during the iterative solution of the electric-field system of equations. For the method presented here, the Earth model is discretised using a rectilinear mesh comprising uniform cells. Edge-element basis functions are used to approximate the electric field and nodal basis functions are used to approximate the correction potential. The Galerkin method is used to derive the systems of equations for the approximate electric field and correction potential from the respective differential equations. A bi-conjugate gradient solver was found to be adequate for the system of equations for the correction potential; a generalised minimum residual solver was found to be better for the electric-field system of equations. The method is illustrated using the COMMEMI 3D-1A and 3D-2A models.     

On the effects of wind and tides on the hydrodynamics of a shallow mediterranean estuary

A study is made of the effect of wind and tides on the hydrodynamics of the shallow inner basins of mediterranean estuaries. The paper includes a case study of Harvey Estuary in southwestern Australia where salinity and temperature data exist for 11 years during the 1980s and 1990s when that estuary experienced massive annual blue-green algal blooms. An analysis is made of salt exchange through the channels that join estuarine basins of this class to either the ocean or, as in the case of Harvey Estuary, to another shallow estuarine basin. A detailed three-dimensional numerical model is also implemented for the basin of Harvey Estuary. It is concluded that exchange through the channel is dominated by the (mainly diurnal) tides, despite the general micro-tidal nature of this class of estuary, although the efficiency of this process is found to be controlled by the length of the channel. Wind set-up in the basin also produces channel exchange and for Harvey Estuary this is about 20% of the exchange due to tides. Baroclinic flow through the channel is also capable of producing significant exchange but this is suppressed by the tidal currents in the channel except immediately after riverflow. Salt transport along the basins of this class of estuary is mainly driven by the longitudinal density gradient and the strength of this process is controlled by vertical mixing from the wind. However, there is also significant salt transport from wind-induced advection, the effect of which changes seasonally with the direction of the salt gradient.     

Lateral dynamics and associated transport of sediment in the upper reaches of a partially mixed estuary,Chesapeake Bay,USA

Data from time series of transects made over a tidal period across a section of the upper Chesapeake Bay, USA, reveal the influence of lateral dynamics on sediment transport in an area with a deep channel and broad extents of shallower flanks. Contributions to lateral momentum by rotation (Coriolis plus channel curvature), cross channel density gradients and cross channel surface slope were estimated, and the friction and acceleration terms needed to complete the balance were compared to patterns of observed lateral circulation. During ebb, net rotation effects were larger because of river velocity and reinforcement of Coriolis by curvature. During flood, stratification was greater because of landward advection of strong vertical density gradients. Together, the ebb intensified lateral circulation and flood intensified stratification focused sediment and sediment transport along the left side of the estuary (looking seaward). The tendency for greater stratification on flood and net sediment flux toward the left-hand shoal are contrary to more common models which, in the northern hemisphere, predict greater resuspension on flood and move sediment toward the right-hand shoal. These tidal asymmetries interact with the lateral circulation to focus net sediment flux on the left side of the estuary, and to produce net ebb directed sediment transport at the surface of the same order of magnitude as net flood directed sediment transport at the bottom.