How tides and river flows determine estuarine bathymetries

For strongly tidal, funnel-shaped estuaries, we examine how tides and river flows determine size and shape. We also consider how long it takes for bathymetric adjustment, both to determine whether present-day bathymetry reflects prevailing forcing and how rapidly changes might occur under future forcing scenarios.Starting with the assumption of a 'synchronous' estuary (i.e., where the sea surface slope resulting from the axial gradient in phase of tidal elevation significantly exceeds the gradient in tidal amplitude ), an expression is derived for the slope of the sea bed. Thence, by integration we derive expressions for the axial depth profile and estuarine length, L, as a function of and D, the prescribed depth at the mouth. Calculated values of L are broadly consistent with observations. The synchronous estuary approach enables a number of dynamical parameters to be directly calculated and conveniently illustrated as functions of and D, namely: current amplitude Û, ratio of friction to inertia terms, estuarine length, stratification, saline intrusion length, flushing time, mean suspended sediment concentration and sediment in-fill times.Four separate derivations for the length of saline intrusion, L_I, all indicate a dependency on (U_o is the residual river flow velocity and f is the bed friction coefficient). Likely bathymetries for `mixed' estuaries can be delineated by mapping, against and D, the conditions L<sub>I</sub>/L<1,E<sub>X</sub>/L<1 (E<sub>X</sub> is the tidal excursion) alongside the Simpson-Hunter criteria D/U<sup>3</sup><50 m<sup>−2</sup> s<sup>3</sup>. This zone encompasses 24 out of 25 `randomly' selected UK estuaries.However, the length of saline intrusion in a funnel-shaped estuary is also sensitive to axial location. Observations suggest that this location corresponds to a minimum in landward intrusion of salt. By combining the derived expressions for L and L_I with this latter criterion, an expression is derived relating D_i, the depth at the centre of the intrusion, to the corresponding value of U_o. This expression indicates U_o is always close to 1 cm s⁻¹, as commonly observed. Converting from U_o to river flow, Q, provides a morphological expression linking estuarine depth to Q (with a small dependence on side slope gradients).These dynamical solutions are coupled with further generalised theory related to depth and time-mean, suspended sediment concentrations (as functions of and D). Then, by assuming the transport of fine marine sediments approximates that of a dissolved tracer, the rate of estuarine supply can be determined by combining these derived mean concentrations with estimates of flushing time, F_T, based on L_I. By further assuming that all such sediments are deposited, minimum times for these deposition rates to in-fill estuaries are determined. These times range from a decade for the shortest, shallowest estuaries to upwards of millennia in longer, deeper estuaries with smaller tidal ranges.     

Simulating a heavy metal spill under estuarine conditions: Effects on the clam Scrobicularia plana

We describe the effect of heavy metals Zn, Cd, Pb and Cu on the induction of methallothioneins on the clam Scrobicularia plana along a salinity gradient simulated under laboratory conditions. The clams were exposed to constant heavy metal concentrations in a dynamic estuary simulator during a 15-day assay to investigate possible induction of metal-binding proteins in them. The concentration of heavy metals in water was analysed. Clams were analysed for methallothionein concentrations. The speciation of Zn, Cd, Pb and Cu along the salinity gradient was modelled. Zn showed the highest concentrations and its prevalent species was the free ion. Intersite differences have been observed in methallothionein concentration and related to the salinity gradient. It seems that synthesis of methallothioneins is the result of physiological forces acting in concert with the changes in the chemical speciation of metals, owing to the trace metals uptake is controlled by means of an interaction of physiology and physicochemistry.     

Modelling the equilibrium speciation of nickel in the Tweed Estuary, UK: Voltammetric determinations and simulations using WHAM

The complexation of dissolved Ni has been evaluated in a rapidly-flushed, rural estuary (Tweed, UK) by ligand exchange-adsorptive cathodic stripping voltammetry. Results suggest the presence of strongly binding ligands, L, throughout, with average stability constants of about 10¹⁹ and which are saturated by ambient Ni concentrations. Equilibrium speciation calculations incorporating these constants in WHAM, version 6, predict an increase in Ni complexation (as NiL) from about 50% of total dissolved Ni in fresh water to over 90% in sea water. Equivalent calculations using the default-mode fulvic and humic substances (FS and HS, respectively) encoded in the WHAM database predict a reduction in complexation (as NiFS + NiHS) from about 20% in fresh water to less than 1% in sea water. Discrepancies arising from the two approaches are largely attributed to the different analytical detection windows employed. Thus, a better representation of Ni complexation is derived from including both types of complexant in the speciation calculations, resulting in estimates of net complexation in excess of 60% of total dissolved Ni throughout the estuary. The uncertainties and assumptions inherent in all computations illustrate the difficulty in measuring or predicting metal complexation in estuaries.     

Preservation conditions and the use of sediment pigments as a tool for recent ecological reconstruction in four Northern European estuaries

Down-core sediment pigment concentrations from four Northern European estuaries were measured using high-performance liquid chromatography (HPLC) to investigate phytoplankton community structure and preservation conditions over the last ca. 100 years where all sites have experienced different levels of eutrophication. Phytoplankton pigments have been shown to be useful biomarkers for phytoplankton community structure and abundance due to their taxonomic specificity. The pigment concentrations and sediment pigment inventory showed large variation between the four sites. Concentrations ranged from more than 6000 nmol/g OC to less than 100 nmol/g OC and the inventory integrated over the top 10 cm from more than 300 nmol/cm² to less than 30 nmol/cm² for total identified pigments. Good pigment preservation in Mariager Fjord (Denmark) reflected the almost permanently anoxic conditions. Pigments in Laajalahti (Finland) showed peak concentrations around the time of highest nitrogen loading events known from historical and modelled records over the past 100 years. In contrast, poor down-core preservation of pigments (especially carotenoids) was observed in the Ems-Dollard (The Netherlands) and Himmerfjärden (Sweden) estuaries. The Ems-Dollard site is an intertidal mudflat that experiences daily exposure to light and air, which enhances pigment degradation. In Himmerfjärden, resuspension is an important process affecting both the sedimentation rate and degradation properties. The different preservation conditions at the four sites were supported by the differences in two degradation indicators; the ratio of pheopigment-a to chlorophyll-a and total carotenoids to total pigments. Class-specific carotenoid pigments represented the dominant algal groups reported from each site, however, no distinct down-core changes in the pigment composition were observed at any of the four sites. This indicated that changes in plankton community structure on the group level have been limited over this time period or masked by low preservation of pigments.     

Spectrophotometric pH measurement in estuaries using thymol blue and m-cresol purple

The conditional acid dissociation constants (pK_a′) of two sulfonephthalein dyes, thymol blue (TB) and m-cresol purple (mCP), were assessed throughout the estuarine salinity range (0<S<40) using a tris/tris–HCl buffer and spectrophotometric measurement. The salinity dependence of the pK_a′ of both dyes was fitted to the equations (25 °C, total proton pH scale, mol kg soln⁻¹):

The estimated accuracy of pH measurements using these calculated pK_a′ values is considered to be comparable to that possible with careful use of a glass electrode (±0.01 pH unit) but spectrophotometric measurements in an estuary have the significant advantage that it is not necessary to calibrate an electrode at different salinities. pH was measured in an estuary over a tidal cycle with a precision of ±0.0005 pH unit at high (S>30) salinity, and ±0.002 pH unit at low (S<5) salinity. The pH increased rapidly in the lower salinity ranges (0<S<15) but less rapidly at higher salinities.     

New typologies for estuarine morphology

New theories to explain how tides and river flow determine estuarine bathymetries have been developed. The range of validity of these theories has been subsequently assessed against a recently available observational data set.Success in reproducing the evolution of estuarine bathymetries that has occurred over the last 10 000 years provides the basis for the translation, here, of these theories into new typologies. These then provide reference frameworks to enhance our perspective on existing morphologies and identify ‘anomalous’ estuaries. Further refinement of these new theories requires analyses of causal factors responsible for such anomalies drawing on classical morphological experience. These frameworks can also be used for calculating likely future bathymetric adjustments for prescribed changes in mean sea level, tidal amplitudes, river flow and sediment supply.     

Sediment origin and budget in Sepetiba Bay (Brazil) - an approach based on multielemental analysis

 Sediment origin and transport were determined in a shallow 447 km² coastal lagoon in Brazil, using the distribution of major elements in bottom and suspended sediments. Applying multivariate analysis, the sum of normalized concentrations of Ti, Fe, and Al in sediments was selected to trace the terrigenous influence in bottom sediments, whereas Ca, Si, and K were used as indicators of marine influence. Sepetiba Bay bottom sediments are dominated by up to 80% terrigenous inputs. The open bay sediments are enriched in P, Mg, K, and Ca because of the intensive water column primary production. The inner bay forms an independent circulation cell with a predominance of terrigenous sediments. The results suggest a long residence time for particles in the bay, demonstrated by the magnitude of resuspension flux compared to the total sediment input to the bay. Received: 29 May 1996 · Accepted: 17 December 1996     

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     

SELFE: A semi-implicit Eulerian–Lagrangian finite-element model for cross-scale ocean circulation

Unstructured-grid models grounded on semi-implicit, finite-volume, Eulerian–Lagrangian algorithms, such as UnTRIM and ELCIRC, have enjoyed considerable success recently in simulating 3D estuarine and coastal circulation. However, opportunities for improving the accuracy of this type of models were identified during extensive simulations of a tightly coupled estuary–plume–shelf system in the Columbia River system. Efforts to improve numerical accuracy resulted in SELFE, a new finite-element model for cross-scale ocean modeling. SELFE retains key benefits, including computational efficiency of existing semi-implicit Eulerian–Lagrangian finite-volume models, but relaxes restrictions on grids, uses higher-order shape functions for elevation, and enables superior flexibility in representing the bathymetry. Better representation of the bathymetry is enabled by a novel, “localized” vertical grid that resembles unstructured grids. At a particular horizontal location, SELFE uses either S coordinates or SZ coordinates, but the equations are consistently solved in Z space. SELFE also performs well relative to volume conservation and spurious oscillations, two problems that plague some finite-element models. This paper introduces SELFE as an open-source code available for community use and enhancement. The main focus here is on describing the formulation of the model and on showing results for a range of progressively demanding benchmark tests. While leaving details to separate publications, we also briefly illustrate the superior performance of SELFE over ELCIRC in a field application to the Columbia River estuary and plume.     

Kinetic and equilibrium studies of copper-dissolved organic matter complexation in water column of the stratified Krka River estuary (Croatia)

An interaction of dissolved natural organic matter (DNOM) with copper ions in the water column of the stratified Krka River estuary (Croatia) was studied. The experimental methodology was based on the differential pulse anodic stripping voltammetric (DPASV) determination of labile copper species by titrating the sample using increments of copper additions uniformly distributed on the logarithmic scale. A classical at-equilibrium approach (determination of copper complexing capacity, CuCC) and a kinetic approach (tracing of equilibrium reconstitution) of copper complexation were considered and compared. A model of discrete distribution of organic ligands forming inert copper complexes was applied. For both approaches, a home-written fitting program was used for the determination of apparent stability constants (K_i^equ), total ligands concentration (L_iT) and association/dissociation rate constants (k_i¹,k_i^- 1).A non-conservative behaviour of dissolved organic matter (DOC) and total copper concentration in a water column was registered. An enhanced biological activity at the freshwater–seawater interface (FSI) triggered an increase of total copper concentration and total ligand concentration in this water layer. The copper complexation in fresh water of Krka River was characterised by one type of binding ligands, while in most of the estuarine and marine samples two classes of ligands were identified. The distribution of apparent stability constants (log K₁^equ: 11.2–13.0, log K₂^equ:8.8–10.0) showed increasing trend towards higher salinities, indicating stronger copper complexation by autochthonous seawater organic matter.Copper complexation parameters (ligand concentrations and apparent stability constants) obtained by at-equilibrium model are in very good accordance with those of kinetic model. Calculated association rate constants (k₁¹:6.1–20 × 10³ (M s)^− 1, k₂¹: 1.3–6.3 × 10³ (M s)^− 1) indicate that copper complexation by DNOM takes place relatively slowly. The time needed to achieve a new pseudo-equilibrium induced by an increase of copper concentration (which is common for Krka River estuary during summer period due to the nautical traffic), is estimated to be from 2 to 4 h.It is found that in such oligotrophic environment (dissolved organic carbon content under 83 µM_C, i.e. 1 mg_CL^− 1) an increase of the total copper concentration above 12 nM could enhance a free copper concentration exceeding the level considered as potentially toxic for microorganisms (10 pM).