Modeling the Effects of Tidal Energy Extraction on Estuarine Hydrodynamics in a Stratified Estuary

A 3-D coastal ocean model with a tidal turbine module was used in this paper to study the effects of tidal energy extraction on temperature and salinity stratification and density-driven two-layer estuarine circulation. Numerical experiments with various turbine array configurations were carried out to investigate the changes in tidally averaged temperature, salinity, and velocity profiles in an idealized stratified estuary that connects to coastal water through a narrow tidal channel. The model was driven by tides, river inflow, and sea surface heat flux. To represent the realistic size of commercial tidal farms, model simulations were conducted based on a small percentage (less than 10 %) of the total number of turbines that would generate the maximum extractable energy in the system. Model results show that extraction of tidal in-stream energy will increase the vertical mixing and decrease the stratification in the estuary. Installation of in-stream tidal farm will cause a phase lag in tidal wave, which leads to large differences in tidal currents between baseline and tidal farm conditions. Extraction of tidal energy in an estuarine system has stronger impact on the tidally averaged salinity, temperature, and velocity in the surface layer than the bottom layer even though the turbine hub height is close to the bottom. Finally, model results also indicate that extraction of tidal energy weakens the two-layer estuarine circulation, especially during neap tides when tidal mixing is weakest and energy extraction is smallest.     

Hydrography and circulation of the Chubut River estuary (Argentina)

The hydrography and circulation of the Chubut River were investigated under exceptionally low river discharge. The frontal zone formed by the entrance of the tide in the estuary may be observed as far as 4.5 km from the mouth, showing that the salt intrusion due to tidal effects reaches further inland than during normal river discharge. Based on the classification of Hansen and Rattray (1966), the estuary corresponds to Type 1 with some vertical stratification observed on the seaward side of the frontal zone. A lateral salinity gradient was found, which was not the result of Coriolis force. The general morphology of the estuary and the consequent secondary circulation due to meanders and interchannel bars may explain the lateral variation. Wind effect is a major component of the circulation and mixing of this shallow estuary.     

The calculation of estuarine turnover times using freshwater fraction and tidal prism models: A critical evaluation

Freshwater fraction and tidal prism models are simple methods for estimating the turnover time of estuarine water. The freshwater fraction method prominently features flushing by freshwater inflow and has sometimes been criticized because it appears not to include flushing by seawater, but this is accounted for implicitly because the average estuary salinity used in the calculation reflects all the processes that bring seawater into the estuary, including gravitational circulation and tidal processes. The model relies on measurable salinity differences among water masses and so must be used for estuaries with substantial freshwater inflow. Tidal prism models are based on flushing by flood tide inflow and ignore seawater inflow due to gravitational circulation. These models should only be applied to estuaries with weak or nonexistent gravitational circulation, which are generally those with little freshwater inflow. Using a framework that is less ambioguous and more directly applicable to the estimation of turnover times than those used previously, this paper critically examines the application of tidal prism models in well-mixed estuaries with complete tidal exchange, partial ebb return, or incomplete flood mixing and in partially mixed estuaries. Problems with self-consistency in earlier versions of these models also apply to the budgeting procedure used by the LOICZ (Land-Ocean Interactions in the Coastal Zone) program. Although freshwater fraction and tidal prism models are different approaches to estimating turnover times in systems with very different characteristics, consistent derivation shows that these models have much in common with each other and that they yield equivalent values that can be used to make comparisons across systems.     

Computed and observed currents,elevations, and salinity in a branching estuary

A one-dimensional, hydrodynamical model of the Tamar Estuary shows good agreement with measured tidal elevations and currents. Computed currents are used to drive a one-dimensional moving-element model of the salt balance. The moving-element model overcomes the numerical difficulties associated with strong tidal advection. Axial distributions of salinity at high water, computed using the moving-element model, compare well with measurements. The modelled and observed high water salinity distributions in this macrotidal estuary show little dependence on tidal range. The major variability in salinity is due to runoff. This strong and rapid dependence on runoff is a consequence of short residence (or flushing) times. Typically, residence times are less than one day throughout the year in the upper 10 km of estuary. The residence times maximize in summer, reaching 14 d for the whole estuary. During high runoff winter periods residence times are less than 5 d. Mixing coefficients for the moving-element salinity model are deduced from salinity measurements. Dispersion coefficients at fixed locations along the estuary are deduced from solutions of the salinity model. The spatially-averaged coefficients at mean spring and neap tides are 180 and 240 m² s^?1, respectively, for average runoff. Therefore, spring-neap variations in dispersion are fairly small and show a negative correlation with tidal range. The spatially-averaged dispersion coefficients at mean tides vary from 150 to 300 m² s^?1 for typical summer and winter runoff, respectively. The increase in dispersion with runoff and the decrease with tidal range implies that buoyancy-driven currents generate an important component of the shear dispersion in this estuary.     

The application of a segmented tidal mixing model to the Great Bay Estuary,N.H.

Measurements show that in general salt is vertically well-mixed everywhere in the Great Bay Estuary, New Hampshire except near the river entrances at the head of the estuary. Dyer and Taylor’s (1973) modified version of Ketchum’s segmented tidal prism model has been applied to the Great Bay Estuarine System in order to predict high and low water salinity distribution for a specified river flow. The theory has been modified here to account for the mixing which occurs at the junction of two branches of an estuary. The mixing parameter, which in this model is related to the tidal excursion of water in the estuary, has been determined for different segments in the estuary on the basis of a comparison between predictions and a comprehensive data set obtained for a low river flow period. Using a mixing parameter distribution based on the low river flow calibration procedure the salinity distribution has been predicted for high river flow. The resulting salinity distribution compares favorably with observations for most of the estuary. The corresponding flushing times for water parcels entering at the head of the estuary during periods of low and high river flow is 54.5 and 45.9 tidal cycles respectively.     

Transport of Crustacean Larvae Between a Low-Inflow Estuary and Coastal Waters

The effectiveness of larval behavior in regulating transport between well-mixed, low-inflow estuaries and coastal waters in seasonally arid climates is poorly known. We determined the flux of an assemblage of benthic crustacean larvae relative to physical conditions between a shallow estuary and coastal waters on the upwelling coast of northern California (38°18′N, 123°03′W) from 29 to 31 March 2006. We detected larval behaviors that regulate transport in adjacent coastal waters and other estuaries for only two taxa in the low-inflow estuary, but they were apparent for taxa outside the estuary. Vertical mixing in the shallow estuary may have overwhelmed larvae of some species, or salinity fluctuations may have been too slight to cue tidal vertical migrations. Nevertheless, all larval stages of species that complete development in nearshore coastal waters were present in the estuary, because they remained low in the water column reducing seaward advection or they were readily exchanged between the estuary and open coast by tidal flows. Weak tidal flows and gravitational circulation at the head of the estuary reduced seaward transport during development for species that completed development nearshore, whereas larval release during nocturnal ebb tides enhanced seaward transport for species that develop offshore. Thus, nonselective tidal processes dominated larval transport for most species back and forth between the low-inflow estuary and open coastal waters, whereas in adjacent open coastal waters, larval behavior in the presence of wind-induced shear was more important in regulating migrations between adult and larval habitats along this upwelling coast.     

Impact of Channel Deepening on Tidal and Gravitational Circulation in a Highly Engineered Estuarine Basin

Deepening of estuarine channels is a common practice to ensure navigation. Here, we investigate whether such deepening impacts physical processes such as the strength of the estuarine exchange flow, the horizontal salinity gradient, and tidal dynamics. We analyze recent and historical hydrodynamic observations in Newark Bay, New Jersey, to assess the effect of channel deepening on tides, circulation, and salinity. The Bay’s navigational channel has undergone significant deepening, from 3 to 10 m in the nineteenth century to ~16 m today. Observations presented here include sea-level data from the nineteenth, twentieth, and twenty-first century, and moored Doppler current data and bottom salinity measurements made over the past 20 years. Results show a doubling of the estuarine exchange flow, a slight increase in salinity and in the horizontal salinity gradient, a decrease in tidal current amplitude, and a spatially variable change in the tidal range. The doubling of the exchange flow is consistent with the Hansen and Rattray scaling provided that the horizontal salinity gradient is unable to fully adjust landward because the dredging is limited to a short reach of the estuary. However, uncertainty in channel depth leaves open the possibility that the exchange flow is also augmented by an increase in the horizontal salinity gradient and/or a reduction in vertical mixing. Nevertheless, results demonstrate that a relatively small (15%) increase in depth appears to have doubled the exchange flow. We believe that this result is relevant to other systems where dredging is limited to a short reach of an estuary.     

Seasonal Variability and Estuary–Shelf Interactions in Circulation Dynamics of a River-dominated Estuary

The long-term response of circulation processes to external forcing has been quantified for the Columbia River estuary using in situ data from an existing coastal observatory. Circulation patterns were determined from four Acoustic Doppler Profilers (ADP) and several conductivity–temperature sensors placed in the two main channels. Because of the very strong river discharge, baroclinic processes play a crucial role in the circulation dynamics, and the interaction of the tidal and subtidal baroclinic pressure gradients plays a major role in structuring the velocity field. The input of river flow and the resulting low-frequency flow dynamics in the two channels are quite distinct. Current and salinity data were analyzed on two time scales—subtidal (or residual) and tidal (both diurnal and semidiurnal components). The residual currents in both channels usually showed a classical two-layer baroclinic circulation system with inflow at the bottom and outflow near the surface. However, this two-layer system is transient and breaks down under strong discharge and tidal conditions because of enhanced vertical mixing. Influence of shelf winds on estuarine processes was also observed via the interactions with upwelling and downwelling processes and coastal plume transport. The transient nature of residual inflow affects the long-term transport characteristics of the estuary. Effects of vertical mixing could also be seen at the tidal time scale. Tidal velocities were separated into their diurnal and semidiurnal components using continuous wavelet transforms to account for the nonstationary nature of velocity amplitudes. The vertical structure of velocity amplitudes were considerably altered by baroclinic gradients. This was particularly true for the diurnal components, where tidal asymmetry led to stronger tidal velocities near the bottom.     

Seasonal variation of the salinity in the Zuari estuary,Goa, India

The annual salt budget of the Zuari is examined. The characteristics of the estuary differ markedly from the low run off season during November–May to the heavy run off period of the southwest monsoon from June to October. During November–May the estuary is vertically mixed and the two processes controlling the transport of salt are run off induced advective transport out of the estuary, and tidally induced diffusive transport into the estuary. The magnitude of the latter is about 20% larger, leading to a salinity rise in the estuary. The diffusion coefficient has been estimated to be 233 ± 101 m²/sec. With the onset of the southwest monsoon, the run off increases dramatically, and the estuary loses about 75% of its salt during the first two months of the season. About 2/3 of this loss is recovered in the next two months when the run off decreases. Because the estuary is partially stratified during June–October, gravitational circulation is expected to play a role in addition to tidal diffusion and run off. The magnitude of its contribution has, however, not yet been determined.     

Evidence of dynamics reversal in tropical estuaries, geomorphological and sedimentological consequences (Salum and Casamance Rivers, Senegal)

South Dakar Senegambian estuaries are subject to an unusual hydrodynamical regime caused by weak or absent run-off. In the Salum delta, each distributary lacks fresh water during most of the year. Only the tidal flows are responsible for geomorphological and sedimentological effects. The current distribution shows a net discharge upstream due to the extensive evaporation and evapotranspiration in mangrove swamps and tidal flats. Consequently the salinity is always higher towards the river than near the sea. A high salinity bottom layer suggests the occurrence of a supersaline wedge of reverse sense to the salt wedge of a normal estuary. Such an inverse pattern is similarly displayed by sedimentological features (double upstream turned spits) and by the external location of the turbidity maximum. A coherent reverse estuary model is suggested from our field observations.     

Time-dependent stratification in the Gauthami-Godavari estuary

The time-dependent salinity stratification in an environment of highly variable turbulence of Gauthami-Godavari estuary is described. From time series observations at three stations over two tidal cycles each in the winter (December 1989) and the dry (April 1990) seasons in the Gauthami-Godavari estuary, the development of stratification is correlated with periods of substantially reduced velocity shear. Stratification is observed to be greatest around the low water slacks and least around the high water slacks. The formation of stratification relaxes viscous constraints and a buoyancy circulation rapidly develops. The breakdown of stratification drastically modifies the circulation and largely removes the vertical shear associated with the density driven flow. Destratification occurs around the high water slacks in the lower reaches fairly close to the mouth of the estuary. The variations in the fields of mass will strongly affect the response of the velocity field.     

Circulation and mixing due to tidal forcing in the Bertioga Channel, São Paulo, Brazil

Bertioga Channel is a partially mixed (type 2) tidal estuary on the coastal plain of São Paulo, Brazil. Hourly current and salinity measurements during neap and spring tides in July 1991 yielded information about the physical structure of the system. Peak along-channel velocities varied from 40 cm s^?1 to 60 cm s^?1 during flood tides and from 70 cm s^?1 to 100 cm s^?1 during ebb tides. Net vertical velocity profiles indicate that the net current reverses directions at a depth of 2.5–3.0 m in the halocline. Due to appreciable fortnightly tidal modulation, the estuary alternates from being highly stratified (type 2b) during neap tides, with advection and diffusion contributing equally to the net upstream salt flux, to being moderately stratified (type 2a) during spring tides, when 90% of the net upstream salt transport is the result of effective tidal diffusion. Decomposition of the salt flux indicates that the relative contribution to the upstream salt transport by gravitational circulation shear is greater than the oscillatory tidal flux by a factor of 2.6 during neap tides. The oscillatory tidal flux is generated by the correlation of the tidal components of the u-velocity and salinity and is responsible for approximately the same amount of upstream salt transport, during neap and spring tides. However, during spring tides, this oscillatory term is greater than the other salt flux terms by a factor of 1.4. The total salt transport, through a unit width of the section perpendicular to the flow, was within 2% of the sum of the seven major decomposed, advective and dispersive terms. On the assumption that the Bertioga Channel is laterally homogeneous, the results also indicate that the estuary is not in steady state with respect to salt flux.     

Transport in the Hudson estuary: A modeling study of estuarine circulation and tidal trapping

The effects of estuarine circulation and tidal trapping on transport in the Hudson estuary were investigated by a large-scale, high-resolution numerical model simulation of a tracer release. The modeled and measured longitudinal profiles of surface tracer concentrations (plumes) differ from the ideal Gaussian shape in two ways: on a large scale the plume is asymmetric with the downstream end stretching out farther, and small-scale (1–2 km) peaks are present at the upstream and downstream ends of the plume. A number of diagnostic model simulations (e.g., remove freshwater flow) were performed to understand the processes responsible for these features. These simulations show that the large-scale asymmetry is related to salinity. The salt causes an estuarine circulation that decreases vertical mixing (vertical density gradient), increases longitudinal dispersion (increased vertical and lateral gradients in longitudinal velocities), and increases net downstream velocities in the surface layer. Since salinity intrusion is confined to the downstream end of the tracer plume, only that part of the plume is effected by those processes, which leads to the largescale asymmetry. The small-scale peaks are due to tidal trapping. Small embayments along the estuary trap water and tracer as the plume passes by in the main channel. When the plume in the main channel has passed, the tracer is released back to the main channel, causing a secondary peak in the longitudinal profile.     

The behavior of trace metals in the Geum Estuary,Korea

The distributions of trace metals in the Geum Estuary of western Korea were studied with regard to changes in other estuarine chemical parameters. Dissolved oxygen, pH, and alkalinity increased with increasing salinity. Dissolved aluminum concentrations increased at low salinities and were perhaps influenced by the solubility of particulate aluminosilicate phases. Iron, managanese, cobalt, and zinc are removed from solution in the low salinity end of the estuary. Cobalt and nickel have mid-estuary concentration maxima that may be due to an anthropogenic source. Cadmium, copper, lead, and zinc concentrations also increased in the estuary, possibly as the result of remobilization in the sediments. Cadmium increases are also linked to remineralization from tidal flat sediments in the outer estuary. The source of an increase in dissolved lead at low salinity is unclear, but may be due to release from particles.     

Axial zonation patterns of subtidal macrozoobenthos in the Gamtoos estuary,South Africa

The distribution of macroinfauna was quantified in subtidal, soft-bottom habitats, extending from the estuarine mouth to the tidal head of the Gamtoos—a small, shallow, temperate estuary situated on the south coast of South Africa. Sampling covered the full salinity gradient from fresh to marine waters, and all sediment types from marine sands to fluvial silts. A total of 35 taxa was recorded, of which 22 occurred throughout the year. Species richness and diversity declined from the seawater-dominated mouth region toward the fresh water section at the tidal head of the estuary. Sediment type generally bore no clear relation to biotic diversity. A marked drop in salinity between winter and summer sample series (Δ 0.2‰ to 24‰) coincided with a reduction of mean macrofaunal density by 70%, a more seaward relocation, and a compression of axial ranges of most taxa. Numerical classification and ordination of faunistically similar regions and of co-occurring species delineated four habitat zones along the longitudinal axis of the estuary which harbour four distinct macrofaunal assemblages: 1) A tidal inlet area with salinities close to seawater; clean, coarse, marine sands, rich in CaCO₃ harbour a stenohaline fauna normally found on adjacent, marine sandy beaches. 2) In the lower reaches, where fine, fluvial silts of high organic content prevail, euryhaline polychaetes dominate the macrozoobenthic community; bottom salinities in this zone seldom dropped below 25‰ 3) The middle reaches, characterized by oligohaline- to polyhaline waters, stretch over sandy sediments of intermediate carbonate, silt, and organic fractions; the fauna comprises typical estuarine forms, which occurred throughout most of the estuary except at its seaward and landward limits. 4) The upper reaches encompass the limnetic waters near the tidal head of the estuary with sediments in this zone being composed mostly of coarse, clean sands, low in CaCO₃; the macrobenthos in this region is dominated by taxa of freshwater origin, which generally do not penetrate seaward beyond the oligohaline waters, and by exceptionally euryhaline estuarine species. Salinity appears as the main factor in controlling faunal assemblages at both extremes of the estuarine gradient (i.e., tidal inlet and head), whereas sediment type delineates between communities in the mesohaline to polyhaline reaches. Axial (i.e., from tidal inlet to tidal head of the estuary) zonation patterns of macroinfauna broadly matched those of mesozooplankton and fishes, supporting the notion of a general structure underlying species distribution patterns in the Gamtoos estuary.     

Effect of a broad shallow sill on tidal circulation and salt transport in the entrance to a coastal plain estuary (Mira—Vila Nova de Milfontes, Portugal)

We describe the tidal circulation and salinity regime of a coastal plain estuary that connects to the ocean through a flood tide delta. The delta acts as a sill, and we examine the mechanisms through which the sill affects exchange of estuarine water with the ocean. Given enough buoyancy, the dynamics of tidal intrusion fronts across the sill and selective withdrawal (aspiration) in the deeper channel landward appear to control the exchange of seawater with estuarine water. Comparison of currents on the sill and stratification in the channel reveals aspiration depths smaller than channel depth during neap tide. During neap tide and strong vertical stratification, seawater plunges beneath the less dense estuarine water somewhere on the sill. Turbulence in the intruding bottom layer on the sill promotes entrainment of fluid from the surface layer, and the seawater along the sill bottom is diluted with estuarine water. During ebb flow, salt is effectively trapped landward of the sill in a stagnant zone between the aspiration depth and the bottom where it can be advected farther upstream by flood currents. During spring tide, the plunge point moves landward and off the sill, stratification is weakened in the deep channel, and aspiration during ebb extends to the bottom. This prevents the formation of stagnant water near the bottom, and the estuary is flooded with high salinity water far inland. The neapspring cycle of tidal intrusion fronts on flood coupled with aspiration during ebb interacts with the sill to play an important role in the transport and retention of salt within the estuary.     

Modeling the transport and distribution of lead in tidal Keelung River estuary

A vertical, two-dimensional heavy metal (lead) transport model incorporated into the hydrodynamic, salt, and sediment transport modules was developed to simulate the lead concentration in the tidal Keelung River estuary of northern Taiwan. We validated the developed model with measured data, including longitudinal velocity, salinity, suspended sediment, and heavy metal (lead) concentration, obtained in 1998. An exponential relationship relating the salinity and suspended-sediment concentrations was established to calculate the partition coefficient of lead in the estuary. The simulated results of dissolved, particulate, and total lead concentrations agreed well with the measured data. A model sensitivity analysis indicated that the partition coefficient plays an important role in the distribution of dissolved and particulate lead concentrations along the tidal Keelung River estuary.     

Estimation of net physical transport and hydraulic residence times for a coastal plain estuary using box models

A box model based on salinity distributions and freshwater inflow measurements was developed and used to estimate net non-tidal physical circulation and hydraulic residence times for Patuxent River estuary, Maryland, a tributary estuary of Chesapeake Bay. The box model relaxes the usual assumption that salinity is at steady-state, an important improvement over previous box model studies, yet it remains simple enough to have broad appeal. Average monthly 2-dimensional net non-tidal circulation and residence times for 1986–1995 are estimated and related to river flow and salt water inflow as estimated by the box model. An important result is that advective exchange at the estuary mouth was not correlated with Patuxent River flow, most likely due to effects of offshore salinity changes in Chesapeake Bay. The median residence time for freshwater entering at the head of the estuary was 68 d and decreased hyperbolically with increasing river flow to 30 d during high flow. Estimates of residence times for down-estuary points of origin showed that, from the head of the estuary to its mouth, control of flushing changed from primarily river flow to other factors regulating the intensity of gravitational circulation.     

Alteration of Residual Circulation Due to Large-Scale Infrastructure in a Coastal Plain Estuary

Large-scale human-built infrastructure is shown to alter the salinity and subtidal residual flow in a realistic numerical simulation of hydrodynamic circulation in a coastal plain estuary (Tampa Bay). Two model scenarios are considered. The first uses a modern bathymetry and boundary conditions from the years 2001–2003. The second is identical to the first except that the bathymetry is based on depth soundings from the pre-construction year 1879. Differences between the models' output can only result from changes in bay morphology, in particular built infrastructure such as bridges, causeways, and dredging of the shipping channel. Thirty-day means of model output are calculated to remove the dominant tidal signals and allow examination of the subtidal dynamics. Infrastructure is found to steepen the mean axial salinity gradient $ \partial \overline{s}/ dx $ by ~40% when there is low freshwater input but flatten $ \partial \overline{s}/ dx $ by ~25% under more typical conditions during moderate freshwater inflow to the estuary. Deepening of the shipping channel also increases the magnitude of the residual Eulerian circulation, allowing for larger up-estuary salt transport. Local bathymetry and morphology are important. Some regions within the estuary show little change in residual circulation due to infrastructure. In others, the residual circulation can vary by a factor of 4 or more. Major features of the circulation and changes due to infrastructure can be partially accounted for with linear theory.     

A modeling study of the Satilla River estuary,Georgia. I: Flooding-drying process and water exchange over the salt marsh-estuary-shelf complex

The flooding-drying process over the intertidal zone of the Satilla River estuary of Georgia was examined using a three-dimensional (3-D) primitive equations numerical model with Mellor and Yamada's (1982) level 2.5 turbulent closure scheme. The model was forced by the semi-diurnal M₂, S₂, and N₂ tides and freshwater discharge at the upstream end of the estuary. The intertidal salt marsh was treated using a 3-D wet-dry point treatment technique that was developed for the σ-coordinate transformation estuary model. Good agreement was found between model-data comparison at anchor monitoring sites and also along the estuary that suggested that the model provided a reasonable simulation of the temporal and spatial distribution of the 3-D tidal current and salinity in the Satilla River estuary. Numerical experiments have shown that the flooding-drying process plays a key role in the simulation of tidal currents in the main river channel and in water transport over the estuarine-salt marsh complex. Ignoring this process could lead to a 50% under-estimation of the amplitude of tidal currents. The model results also revealed a complex spatial structure of the residual flow in the main channel of the river, with characteristics of multiple eddy-like cell circulations. These complicated residual currents are formed due to tidal rectification over variable topography with superimposition of inertial effects, asymmetry of tidal currents, and baroclinic pressure gradients. Water exchanges over the estuary-intertidal salt marsh complex are asymmetric across the estuary, and tend to vary periodically on the northern side while quickly washing out of the marsh zone on the southern side. Strong Stokes’ drifting velocity was predicted in the estuary, so that the Lagrangian trajectories of particles were characterized by strong nonlinear processes that differ significantly from those estimated by the Eulerian residual currents.