Unraveling Interactions between Asymmetric Tidal Turbulence, Residual Circulation, and Salinity Dynamics in Short, Periodically Weakly Stratified Estuaries

Asymmetric tidal turbulence (ATT) strongly influences estuarine health and functioning. However, its impact on the three-dimensional estuarine dynamics and the feedback of water motion and salinity distribution on ATT remain poorly understood, especially for short estuaries (estuarine length ≪ tidal wavelength). This study systematically investigates the abovementioned interactions in a short estuary for the first time, considering periodically weakly stratified conditions. This is done by developing a three-dimensional semi-analytical model (combining perturbation method with finite element method) that allows a dissection of the contributions of different processes to ATT, estuarine circulation, and salt transport. The generation of ATT is dominated by (i) strain-induced periodic stratification and (ii) asymmetric bottom-shear-generated turbulence, and their contributions to ATT are different both in amplitude and phase. The magnitude of the residual circulation related to ATT and the eddy viscosity–shear covariance (ESCO) is about half of that of the gravitational circulation (GC) and shows a “reversed” pattern as compared to GC. ATT generated by strain-induced periodic stratification contributes to an ESCO circulation with a spatial structure similar to GC. This circulation reduces the longitudinal salinity gradients and thus weakens GC. Contrastingly, the ESCO circulation due to asymmetric bottom-shear-generated turbulence shows patterns opposite to GC and acts to enhance GC. Concerning the salinity dynamics at steady state, GC and tidal pumping are equally important to salt import, whereas ESCO circulation yields a significant seaward salt transport. These findings highlight the importance of identifying the sources of ATT to understand its impact on estuarine circulation and salt distribution.

IberWQ: A GPU Accelerated Tool for 2D Water Quality Modeling in Rivers and Estuaries

Numerical models are useful tools to analyze water quality by computing the concentration of physical, chemical and biological parameters. The present work introduces a two-dimensional depth-averaged model that computes the most relevant and frequent parameters used to evaluate water quality. High performance computing (HPC) techniques based on graphic processing unit (GPU) parallelization have been applied to improve the efficiency of the package, providing speed-ups of two orders of magnitude in a standard PC. Several test cases were analyzed to show the capabilities and efficiency of the model to evaluate the environmental status of rivers and non-stratified estuaries. IberWQ will be freely available through the package Iber.

IberWQ: new simulation tool for 2D water quality modelling in rivers and shallow estuaries

This paper presents a new freeware simulation tool (IberWQ) for 2D water quality modelling in rivers and non-stratified estuaries. The model computes the spatial and temporal evolution of several species and variables which are relevant for the evaluation of the environmental status of rivers and estuaries, including: Escherichia coli, dissolved oxygen, carbonaceous biochemical oxygen demand, organic nitrogen, ammoniacal nitrogen, nitrate–nitrite nitrogen, water temperature and salinity. A depth-averaged transport equation is solved for each variable with a mass conservative unstructured finite volume solver. IberWQ is fully coupled to the hydrodynamic module of the software Iber, a freeware simulation tool for solving the 2D shallow water equations. Both models are integrated in the same windows graphical environment, including several tools which allow the user to generate unstructured meshes adapted to the site topography, define spatially variable input parameters and visualize model outputs. We present four application examples to illustrate the possibilities of the software for the dynamic simulation of water quality in rivers and estuaries.

Decomposition of Residual Circulation in Estuaries

The residual currents in estuaries are produced by a variety of physical mechanisms. To understand the contribution of each individual mechanism to the creation of residual circulation, it is necessary to separate the effect of one particular mechanism from the others. In this study, a method based on dynamics is developed to decompose the residual circulation into individual components corresponding to different forcing mechanisms. Specifically, residual flows are partitioned based on the separate contributions by river discharge, horizontal density gradient, internal tidal asymmetry, advection, semi–Stokes transport, and wind. The method includes the effects of the earth’s rotation and can be applied for general conditions. Under the precondition that the ratio between width and length of the estuary is small, the continuity equation can be simplified such that the method only requires the data at a cross-estuary section to decompose residual currents. This makes the method practicable for real estuaries. Results from a generic numerical model are used to illustrate the decomposition method and to demonstrate its validity for weakly stratified estuaries.

Residual Sediment Fluxes in Weakly-to-Periodically Stratified Estuaries and Tidal Inlets

In this idealized numerical modeling study, the composition of residual sediment fluxes in energetic (e.g., weakly or periodically stratified) tidal estuaries is investigated by means of one-dimensional water column models, with some focus on the sediment availability. Scaling of the underlying dynamic equations shows dependence of the results on the Simpson number (relative strength of horizontal density gradient) and the Rouse number (relative settling velocity) as well as impacts of the Unsteadiness number (relative tidal frequency). Here, the parameter space given by the Simpson and Rouse numbers is mainly investigated. A simple analytical model based on the assumption of stationarity shows that for small Simpson and Rouse numbers sediment flux is down estuary and vice versa for large Simpson and Rouse numbers. A fully dynamic water column model coupled to a second-moment turbulence closure model allows to decompose the sediment flux profiles into contributions from the transport flux (product of subtidal velocity and sediment concentration profiles) and the fluctuation flux profiles (tidal covariance between current velocity and sediment concentration). Three different types of bottom sediment pools are distinguished to vary the sediment availability, by defining a time scale for complete sediment erosion. For short erosion times scales, the transport sediment flux may dominate, but for larger erosion time scales the fluctuation sediment flux largely dominates the tidal sediment flux. When quarter-diurnal components are added to the tidal forcing, up-estuary sediment fluxes are strongly increased for stronger and shorter flood tides and vice versa. The theoretical results are compared to field observations in a tidally energetic inlet.

Estuarine geomorphology and low salinity requirement for fertilisation influence spawning site location in the diadromous fish, Galaxias maculatus

In habitats such as estuaries, which are characterised by large and fluctuating gradients in abiotic variables, finding appropriate habitat for successful spawning and egg development can be critical to a species’ survival. We explored how salinity requirements for successful fertilisation may govern the distribution of estuarine spawning habitat for the diadromous fish, Galaxias maculatus, which spawns in inundated vegetation on estuary banks during spring tides. Artificial fertilisation experiments confirmed that successful fertilisation only occurs at low salinities (<20). Thus, we predicted that egg distributions would depend upon the extent of low-salinity surface waters in an estuary. Using estuary geomorphology classification schemes, which classify estuaries by physical and chemical characteristics such as their salinity dynamics, we hypothesised that stratified estuaries would provide a greater extent of low salinity surface water than well-mixed estuaries. This prediction was supported by surveys of egg distributions in five estuaries in Victoria, Australia. Eggs were distributed over a greater proportion of ‘stratified’ v. ‘mixed’ estuary types. We suggest that combining knowledge of the spawning requirements of a species and physical properties of the habitat, such as those encapsulated in estuary geomorphic classification schemes, can greatly facilitate efforts to identify critical habitats and thus aid in species management and conservation.