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.     

Deep circulation driven by strong vertical mixing in the Timor Basin

The importance of deep mixing in driving the deep part of the overturning circulation has been a long debated question at the global scale. Our observations provide an illustration of this process at the Timor Basin scale of ～1000 km. Long-term averaged moored velocity data at the Timor western sill suggest that a deep circulation is present in the Timor Basin. An inflow transport of ～0.15 Sv is observed between 1600 m and the bottom at 1890 m. Since the basin is closed on its eastern side below 1250 m depth, a return flow must be generated above 1600 m with a ～0.15 Sv outflow. The vertical turbulent diffusivity is inferred from a heat and transport balance at the basin scale and from Thorpe scale analysis. Basin averaged vertical diffusivity is as large as 1 × 10^?3 m² s^?1. Observations are compared with regional low-resolution numerical simulations, and the deep observed circulation is only recovered when a strong vertical diffusivity resulting from the parameterization of internal tidal mixing is considered. Furthermore, the deep vertical mixing appears to be strongly dependent on the choice of the internal tide mixing parameterization and also on the prescribed value of the mixing efficiency.     

The role of suspended sediments and phytoplankton in the partitioning and transport of silver in estuaries

Silver sorbs readily at low salinities to both phytoplankton and suspended sediments. As salinity increases, the degree of sorption decreases. Nearly 80% of silver sorbed to suspended sediments at low salinities will desorb at higher salinities, but desorption does not occur when silver is associated with phytoplankton. Thus, silver incorporation onto/into cellular material will increase the retention of silver within the estuary, reducing its rate of transport.     

Comparison of longitudinal equilibrium profiles of estuaries in idealized and process-based models

The process-based morphodynamic model Delft3D-MOR and the idealized model of schuttelaars and De Swart (2000) are compared with each other. The differences between the two models in their mathematical-physical formulation as well as the boundary conditions are identified. Their effect on producing cross-sectionally averaged morphological equilibria of tidal inlets with arbitrary length and forced at the seaward boundary by a prescribed M₂ and M₄ sea- surface elevation is studied and an inventory is made of all relevant differences. The physical formulations in the source code of Delft3D-MOR are modified in various steps to resemble the formulations in the idealized model. The effect of each of the differences between the idealized and process-based model are studied by comparing the results of the idealized model to those of the adapted process-based model. The results of the idealized model can be qualitatively reproduced by the process-based model as long as the same morphological boundary condition is applied at the open sea end. This means that the simplifications concerning the mathematical formulation of the physical processes in the idealized model can be justified. Furthermore, it can be inferred that the morphological boundary condition at the open sea end is an essential element in controlling the behaviour of morphodynamic models for tidal inlets and estuaries.Responsible Editor: Jens Kappenberg     

Observations of the spring–neap modulation of the gravitational circulation in a partially mixed estuary

The long-term variability of the non-tidal circulation in Southampton Water, a partially mixed estuary, was investigated using 71-day acoustic Doppler current profiler (ADCP) time series. The data show evidence that the spring–neap tidal variability of the turbulent mixing modulates the strength of the non-tidal residual circulation, with subtidal neap tide surface flows reaching 0.12 m s^–1 compared to <0.05 m s^–1 at spring tides. The amplitude of the neap-tide events in this non-tidal circulation is shown to be related to a critical value of the tidal currents, illustrating the strong dependence on tidal mixing. The results suggest that the dominant mechanism for generating these neap-tide circulation events is the baroclinic forcing of the horizontal density gradient, rather than barotropic forcing associated with ebb-induced periodic stratification. While tidal turbulence is thought to be the dominant control on this gravitational circulation, there is evidence of the additional effect of wind-driven mixing, including the effects of wind fetch and possibly wave development with along-estuary winds being more efficient at mixing the estuary than across-estuary winds. Rapid changes in atmospheric pressure also coincided with fluctuations in the gravitational circulation. The observed subtidal flows are shown to be capable of rapidly flushing buoyant material out of the estuary and into the coastal sea at neap tides.Responsible Editor: Iris Grabemann     

Selenium in sediments and biota from estuaries of southwest England

Selenium concentrations have been measured in sediment, fucoid macroalgae and macroinvertebrates from four estuaries of SW England (Yealm, Plym, Looe, Fal). Sediment concentrations ranged from about 0.4 μg g⁻¹ in the Yealm to 1.49 μg g⁻¹ at one site in the Plym. Concentrations in Fucus vesiculosus (0.05–0.31 μg g⁻¹) and F. ceranoides (0.05–0.51 μg g⁻¹) were significantly lower than corresponding concentrations in sediment but there was no correlation between algal and sediment concentrations. Selenium concentrations in Littorina littorea (∼4 μg g⁻¹), Hediste diversicolor (2.82–12.68 μg g⁻¹), Arenicola marina (∼17 μg g⁻¹) and Scrobicularia plana (1.18–6.85 μg g⁻¹) were considerably higher than concentrations in macroalga or sediment, suggesting that Se is effectively accumulated from the diet. Although Se concentrations in some invertebrates exceed toxicity thresholds for the diet of predacious birds and fish, no specific evidence for Se toxicity exists in these estuaries.     

Lateral friction in the oceans as an effect of potential vorticity mixing

Abstract

Applying a mixing-length calculation to potential vorticity rather than to momentum a new type of lateral friction appears in the oceanic mass transport equations. This friction is evaluated for the special case of horizontally homogeneous, quasi-geostrophic turbulence. The main effect is a westward force arising from the so-called β-term. This produces an additional southward interior transport and a strengthening of the western boundary current. A turbulent exchange coefficient K^H = 10⁸ em²s^?1 is sufficient to give a Gulf Stream transport twice that obtained by the classical Sverdrup model.     

Anthropogenic influence on the organic fraction of sediments in two contrasting estuaries: a biochemical approach

Abundance and biochemical composition of organic materials in sediments of the estuaries of Mundaka and Bilbao (SE Bay of Biscay) were analysed to assess the effect of organic wastes released to these systems. Organic and labile contents were higher in the sediments of Bilbao, denoting organic enrichment with poorly decomposed materials by human dumping. Spatially, organic matter distribution skewed seaward in Bilbao and upward in Mundaka, in agreement with the location of major point sources of sewage. Labile material, proteins and protein:carbohydrate ratio showed upward increases in both estuaries, attributed to the oxygen restrictive conditions to decompose organic materials. Vertical distribution of organic components into sediments reflected the history of human activities. In Bilbao, the significant reduction of organic and labile compounds, and the decrease of protein:carbohydrate values in top layers seem a result of recent remedial actions to reduce urban sewage inputs and improve environmental conditions. Higher contributions of proteins and lipids were associated with anthropogenic materials, and higher contributions of carbohydrates with autochthonous materials. High protein:carbohydrate values indicate poorly decomposed materials of human origin. The close relationship of carbohydrates with chloropigments in Mundaka suggests that phytobenthic populations and derived detritus contributed to a greater extent to the organic fraction in this estuary.     

On the tidally driven circulation in the South China Sea: modeling and analysis

The South China Sea is a large marginal sea surrounded by land masses and island chains, and characterized by complex bathymetry and irregular coastlines. An unstructured-grid SUNTANS model is employed to perform depth-averaged simulations of the circulation in the South China Sea. The model is tidally forced at the open ocean boundaries using the eight main tidal constituents as derived from the OSU Tidal Prediction Software. The model simulations are performed for the year 2005 using a time step of 60 s. The model reproduces the spring-neap and diurnal and semidiurnal variability in the observed data. Skill assessment of the model is performed by comparing model-predicted surface elevations with observations. For stations located in the central region of the South China Sea, the root mean squared errors (RMSE) are less than 10 % and the Pearson’s correlation coefficient (r) is as high as 0.9. The simulations show that the quality of the model prediction is dependent on the horizontal grid resolution, coastline accuracy, and boundary locations. The maximum RMSE errors and minimum correlation coefficients occur at Kaohsiung (located in northern South China Sea off Taiwan coast) and Tioman (located in southern South China Sea off Malaysia coast). This may be explained with spectral analysis of sea level residuals and winds, which reveal dynamics at Kaohsiung and Tioman are strongly influenced by the seasonal monsoon winds. Our model demonstrates the importance of tidally driven circulation in the central region of the South China Sea.     

Morphological controls in sandy estuaries: the influence of tidal flats and bathymetry on sediment transport

The morphodynamics of shallow, vertically well-mixed estuaries, characterised by tidal flats and deeper channels, have been investigated. This paper examines what contributes to flood/ebb-dominant sediment transport in localised regions through a 2D model study (using the TELEMAC modelling system). The Dyfi Estuary in Wales, UK has been used as a case study and, together with idealised estuary shapes, shows that shallow water depths lead to flood dominance in the inner estuary whilst tidal flats and deep channels cause ebb dominance in the outer estuary. For medium sands and with an artificially ‘flattened’ bathymetry (i.e. no tidal flats), the net sediment transport switches from ebb-dominant to flood-dominant where the parameter a/h (local tidal amplitude ÷ local tidally averaged water depth) exceeds 1.2. Sea level rise will reduce this critical value of a/h and also reduce the ebb-directed sediment transport significantly, leading to a flood-dominated estuarine system. A similar pattern, albeit with greater transport, was simulated with tidal flats included and also with a reduced grain size. This suggests that analogous classifications for flood/ebb asymmetry of the tide in estuaries as a whole may not represent the local sediment transport in sufficient detail. Through the Dyfi simulations, the above criterion involving a/h is shown to be complicated further by augmented flow past a spit at the estuary mouth which gives rise to a self-maintaining scour hole. Simulations of one year of bed evolution in an idealised flat-bottomed estuary, including tidal flow past a spit, recreate the flood/ebb dominance on either side of the spit and the formation of a scour hole in between. The erosion rate at the centre of the hole is reduced as the hole deepens, suggesting the establishment of a self-maintaining equilibrium state.     

Plutonium and210Pb distributions in northeast Atlantic sediments: subsurface anomalies caused by non-local mixing

The depth distributions of²¹⁰Pb and^239,240Pu measured in a suite of box cores collected from water depths of 4000–5000 m in the northeast Atlantic Ocean exhibit pronounced subsurface maxima caused by sediment reworking by benthic infauna. Small-scale spatial heterogeneity in bioturbation rates is indicated by large differences in tracer profiles from duplicate cores separated only by a few centimeters.²¹⁰Pb and^239,240Pu activity distributions from each subcore exhibit a high degree of correlation, and most tracer profiles exhibit one or more subsurface maxima.One-dimensional, “biodiffusion” analogue models do not adequately simulate the principal features of this data set. However, an inverse “conveyer belt” mixing model which simulates subsurface egestion (or a functionally equivalent process) of surficial material which is enriched both in organic debris and radioactive tracers can reproduce the subsurface tracer maxima. Single-event and continuous subsurface egestion models have been formulated and solved for different “feeding rates” and background biodiffusive fields. The single-event model provides a better fit to the data and, in particular, ensures the observed, high degree of correlation between the²¹⁰Pb and^239,240Pu activity profiles, regardless of the different tracer input functions. The most likely candidate responsible for subsurface tracer egestion is a large infaunal worm of the phylum Sipunculida which dominates the biomass below a depth of 3 cm.     

Modeling seasonal circulation, upwelling and tidal mixing in the Arafura and Timor Seas

The extensive shallow tropical seas off northern Australia, encompassing the Arafura and Timor Seas, have been identified as one of the most pristine marine environments on the planet. However, the remoteness and the absence of major industrial development that has contributed to this status have the additional consequence that relatively little is known about these systems. This study is the first to model oceanographic conditions across the tidally dominated Arafura and Timor Seas, and their seasonal variability. The results are based on a high-resolution (0.05°) ocean circulation model forced by realistic winds, waves and tides. The main focus of the study is on physical processes that influence the distributions of sediments and primary productivity across the system. Regions of high bottom stress and tidal mixing have been identified, including a large offshore area around Van Diemen Rise (Timor Sea). Lagrangian particle tracks have revealed a seasonal overturning cell that stretches across the Gulf of Carpentaria (Arafura Sea) with upwelling and downwelling on either side of the Gulf. The presence of coastal upwelling and downwelling is shown to provide a dynamically consistent explanation for the persistent turbid boundary layer observed around the shallow coastal waters of the Gulf.     

Vertical structure of residual slope circulation driven by JEBAR and tides: an idealised model

The residual circulation over the continental slope, and in particular, its vertical structure, is analysed by means of an idealised hydrodynamic model. The model is based on the depth-dependent shallow-water equations under uniform along-isobath conditions and is forced by a prescribed meridional density gradient and tidal velocities. By means of expansion in the small Rossby number solutions are analysed for conditions representative for the continental slopes off the Hebrides and in the Bay of Biscay. The steady solution at zeroth order consists of a linear density-driven flow. At order a tidally rectified flow is found and a stationary flow due to self-interaction of the zeroth-order density-driven flow. At order ² the leading-order effect of the interaction between the zeroth-order density-driven flow and the tides is found: the ‘interaction current’. The solutions up to and including order ² constitute an along-isobath steady slope current which is comparable to field data. The slope current and the accompanying cross-shelf circulation depend strongly on the shelf and flow characteristics. For the Hebridean case the density forcing predominates, but for the Biscay case the tidal effects are of the same order of magnitude as the density effects. Under those conditions the interaction current is significant which implies that linear superposition of density and tidal effects differs from the non-linear combination of both. It is also shown that the depth-average of the interaction current differs essentially from the solution obtained from a depth-averaged model.     

Heavy metals in surface sediments from nine estuaries along the coast of Bohai Bay,Northern China

Concentrations of heavy metals in river water and sediment were investigated in nine estuaries along the coast of Bohai Bay, Northern China. Multivariate statistical techniques such as principal component analysis and cluster analysis, in combination with metal concentration analysis and correlation analysis, were used to identify the possible sources of the metals and the pollution pattern in nine estuaries along the coast of Bohai Bay. The environmental risks of metals, evaluated by sediment quality guidelines and background values, revealed Hg contamination in the estuaries. However, levels of Cd in estuarine sediments were low, and they were less than those levels in river sediments, partly due to the high mobility and dilution of river or seawater. Cd did not contribute to sediment deposits in estuaries. High organic matter from effluents from large municipal sewage treatment plants was predominantly responsible for restricting Hg mobility from the river to Bohai Bay.     

Steady flows driven by sources of random strength in heterogeneous aquifers with application to partially penetrating wells

Average steady source flow in heterogeneous porous formations is modelled by regarding the hydraulic conductivity K(x) as a stationary random space function (RSF). As a consequence, the flow variables become RSFs as well, and we are interested into calculating their moments. This problem has been intensively studied in the case of a Neumann type boundary condition at the source. However, there are many applications (such as well-type flows) for which the required boundary condition is that of Dirichlet. In order to fulfill such a requirement the strength of the source must be proportional to K(x), and therefore the source itself results a RSF. To solve flows driven by sources whose strength is spatially variable, we have used a perturbation procedure similar to that developed by Indelman and Abramovich (Water Resour Res 30:3385–3393, 1994) to analyze flows generated by sources of deterministic strength. Due to the linearity of the mathematical problem, we have focused on the explicit derivation of the mean head distribution G _d (x) generated by a unit pulse. Such a distribution represents the fundamental solution to the average flow equations, and it is termed as mean Green function. The function G _d (x) is derived here at the second order of approximation in the variance σ² of the fluctuation (where K _A is the mean value of K(x)), for arbitrary correlation function ρ(x), and any dimensionality d of the flow domain. We represent G _d (x) as product between the homogeneous Green function G _d ⁽⁰⁾(x) valid in a domain with constant K _A , and a distortion term Ψ _d (x) = 1 + σ²ψ _d (x) which modifies G _d ⁽⁰⁾(x) to account for the medium heterogeneity. In the case of isotropic formations ψ _d (x) is expressed via one quadrature. This quadrature can be analytically calculated after adopting specific (e.g.. exponential and Gaussian) shape for ρ(x). These general results are subsequently used to investigate flow toward a partially-penetrating well in a semi-infinite domain. Indeed, we construct a σ²-order approximation to the mean as well as variance of the head by replacing the well with a singular segment. It is shown how the well-length combined with the medium heterogeneity affects the head distribution. We have introduced the concept of equivalent conductivity K ^eq(r,z). The main result is the relationship where the characteristic function ψ^(w)(r,z) adjusts the homogeneous conductivity K _A to account for the impact of the heterogeneity. In this way, a procedure can be developed to identify the aquifer hydraulic properties by means of field-scale head measurements. Finally, in the case of a fully penetrating well we have expressed the equivalent conductivity in analytical form, and we have shown that (being the effective conductivity for mean uniform flow), in agreement with the numerical simulations of Firmani et al. (Water Resour Res 42:W03422, 2006).     

Baroclinic and topographic influences on the transport in western boundary currents

Abstract

One of the central unsolved theoretical problems of the large scale ocean circulation is concerned with explaining the very large transports measured in western boundary currents such as the Gulf Stream and the Kuroshio. The only theory up to now that can explain the size of these transports is that of non-linear recirculation in which the advective terms in the momentum equations became important near the western boundary. In this paper an alternative explanation is suggested. When bottom topography and baroclinic effects are included in a wind-driven ocean model it is shown that the western boundary current can have a transport larger than that predicted from the wind stress distribution even when the nonlinear advective terms are ignored. The explanation lies in the presence of pressure torques associated with bottom topography which can contribute to the vorticity balance in the same sense as the wind stress curl.

Three numerical experiments have been carried out to explore the nature of this process using a three dimensional numerical model. The first calculation is done for a baroclinic ocean of constant depth, the second for a homogeneous ocean with an idealized continental slope topography, and the third for a baroclinic ocean with the same continental slope topography. The nature of the vorticity balance and of the circulation around closed paths is examined in each case, and it is shown that bottom pressure torques lead to enhanced transport in the western boundary current only for the baroclinic case with variable depth.     

Asymptotic solutions of differential rotation driven by convection in rapidly rotating fluid spheres with the non-slip boundary condition

A strong Coriolis force in a rapidly rotating planet and star not only enforces two-dimensionality of fluid motion driven by thermal instabilities but also generates strong differential rotation even in the vicinity of the onset of instabilities. We derive an asymptotic solution describing convection-driven differential rotation in rotating, self-gravitating Boussinesq fluid spheres with the no-slip boundary condition, taking into account full spherical curvature and being valid for asymptotically small Ekman numbers. For the purpose of validating the asymptotic solution, the corresponding numerical analysis valid for large or small Ekman numbers is also carried out, showing a satisfactory agreement between the asymptotic and numerical solutions.     

Interstitial transport of solutes in non-steady state accumulating and compacting sediments

A general one-dimensional equation for interstitial transport in accumulating and compacting sediments under non-steady state conditions is derived. As a consequence of compaction the metric along the path of a given horizon, i.e. the spatial distance between neighbouring particles, changes continually. Special emphasis is put on the treatment of advection caused by compaction. The resulting partial differential equation for the interstitial concentration of a given solute contains terms which can be evaluated based on data from a single sediment core. In addition, an integral over the time-derivative of porosity appears which would make it necessary to compare cores from the same site but at different times. Under quite general assumptions this last term may, however, be transformed into a form for which evaluation from a single core becomes possible. Several special solutions are treated such as total steady state, steady state at the surface, non-constant sedimentation rate with steady state compaction, and non-steady state with steady state compaction. The last case applies, e.g., to diffusion under the influence of changing boundary conditions at the water/sediment interface while the accumulation process remains in steady state.