Flood–ebb and spring–neap variations of mixing,stratification and circulation in Chesapeake Bay

A three-dimensional hydrodynamic model is used to investigate intra-tidal and spring–neap variations of turbulent mixing, stratification and residual circulation in the Chesapeake Bay estuary. Vertical profiles of salinity, velocity and eddy diffusivity show a marked asymmetry between the flood and ebb tides. Tidal mixing in the bottom boundary layer is stronger and penetrates higher on flood than on ebb. This flood–ebb asymmetry results in a north–south asymmetry in turbulent mixing because tidal currents vary out of phase between the lower and upper regions of Chesapeake Bay. The asymmetric tidal mixing causes significant variation of salinity distribution over the flood–ebb tidal cycle but insignificant changes in the residual circulation. Due to the modulation of tidal currents over the spring–neap cycle, turbulent mixing and vertical stratification show large fortnightly and monthly fluctuations. The stratification is not a linear function of the tidal-current amplitude. Strong stratification is only established during those neap tides when low turbulence intensity persists for several days. Residual circulation also shows large variations over the spring–neap cycle. The tidally averaged residual currents are about 50% stronger during the neap tides than during the spring tides.     

Regularities of hydrological processes in the Chesapeake Bay (USA): Case study of a classical estuary

The main regularities of hydrological and hydrological-environmental processes occurring within the complex estuary, the Chesapeake Bay and the mouths of its tributaries, are discussed. The peculiarities of the estuary morphological structure, including the structures of tidal and net currents, salinity and water turbidity fields and their variability, the environmental conditions, and their human-induced changes. Using the Chesapeake Bay as an example, it became possible to reveal the basic features of classical estuaries subject to a considerable impact of river runoff and featuring mixing of river and sea water and moderate stratification of the water mass. It is shown that the regularities of hydrological processes in the Chesapeake Bay are typical of many mouth water bodies of estuarine type (inlets, drowned river valleys, lagoons, and tidal estuaries proper).     

Resonance and sea level variability in Chesapeake Bay

A numerical model is used to determine the resonant period and quality factor Q of Chesapeake Bay and explore physical mechanisms controlling the resonance response in semi-enclosed seas. At the resonant period of 2 days, the mouth-to-head amplitude gain is 1.42 and Q is 0.9, indicating that Chesapeake Bay is a highly dissipative system. The modest amplitude gain results from strong frictional dissipation in shallow water. It is found that the spatial distribution of energy dissipation varies with forcing frequency. While energy at tidal frequencies is dissipated around topographic hotspots distributed throughout the Bay, energy dissipation at subtidal frequencies is mainly concentrated in the shallow-water lower Bay. An analytic calculation shows that the bottom friction parameter is much larger in Chesapeake Bay than in other coastal systems with strong resonance response. The model-predicted amplitude gains and phase changes agree well with the observations at semidiurnal and diurnal tidal frequencies. However, the predicted amplitude gain in the resonant frequency band (34–54 h period) falls below that inferred from band-passed sea level observations. This discrepancy can be attributed to the local wind forcing which amplifies the sea level response in the upper Bay. The model is also used to show that rising sea levels associated with global warming will shift the resonance period of Chesapeake Bay closer to the diurnal tides and thus exacerbate flooding problems by causing an increase in tidal ranges.     

Horizontal patterns of water temperature and salinity in an estuarine tidal channel: Ria de Aveiro

This work presents results from two complementary and interconnected approaches to study water temperature and salinity patterns in an estuarine tidal channel. This channel is one of the four main branches of the Ria de Aveiro, a shallow lagoon located in the Northwest coast of the Iberian Peninsula. Longitudinal and cross-sectional fields of water temperature and salinity were determined by spatial interpolation of field measurements. A numerical model (Mohid) was used in a 2D depth-integrated mode in order to compute water temperature and salinity patterns. The main purpose of this work was to determine the horizontal patterns of water temperature and salinity in the study area, evaluating the effects of the main forcing factors. The field results were depth-integrated and compared to numerical model results. These results obtained using extreme tidal and river runoff forcing, are also presented. The field results reveal that, when the river flow is weak, the tidal intrusion is the main forcing mechanism, generating saline and thermal fronts which migrate with the neap/spring tidal cycle. When the river flow increases, the influence of the freshwater extends almost as far as the mouth of the lagoon and vertical stratification is established. Results of numerical modelling reveal that the implemented model reproduces quite well the observed horizontal patterns. The model was also used to study the hydrology of the study area under extreme forcing conditions. When the model is forced with a low river flow (1 m³ s⁻¹) the results confirm that the hydrology is tidally dominated. When the model is forced with a high river flow (1,000 m³ s⁻¹) the hydrology is dominated by freshwater, as would be expected in such an area.     

Long-term trends in the macrobenthos and water quality of the lower Chesapeake Bay (1985–1991)

Long-term trends in macrobenthic communities of the lower Chesapeake Bay, USA, were examined using data collected quarterly (March, June, September and December) from 1985 to 1991 at 16 stations along a salinity gradient from tidal freshwater regions of the major tributaries (James, York and Rappahannock rivers) to the polyhaline region of the main-stem of Chesapeake Bay. A non-parametric trend analysis procedure was applied to five parameters characterizing macrobenthic community structure: community biomass, species richness, abundance of individuals, proportion of biomass composed of opportunistic species (opportunistic biomass composition) and proportion of biomass composed of equilibrium species (equilibrium biomass composition). For the parameters tested 36 trends were detected. For community biomass, five trends were significant; all had positive slopes and occurred in the James and York rivers. For species richness, six trends were significant; all had positive slopes with three trends in the James River, two trends in the York River and one trend in the main-stem of Chesapeake Bay. For abundance of individuals, 17 trends were detected; all abundance trends were seasonally dependent, had positive slopes and occurred at 12 of the 16 stations. For opportunistic biomass composition, four trends were significant; all had positive slopes with one trend in the lower Rappahannock River and three trends in the main-stem of Chesapeake Bay. For equilibrium biomass composition four trends were significant; two trends had positive slopes (one in the James River and one in the York River) and two trends had negative slopes (one in the Rappahannock River and one in the main-stem of Chesapeake Bay). Trends in the James and York rivers were considered to indicate improving conditions for the benthos, while trends in the lower Rappahannock River and the main-stem of the Chesapeake Bay were considered to indicate deteriorating conditions. Deteriorating conditions for the benthos were associated with regions exposed to summer, low dissolved oxygen events. The trends in the indicators of benthic biological community health were inferentially related to trends observed in water quality conditions in the tributaries and main-stem of Chesapeake Bay. All major water quality and biotic trends appeared to correspond in an ecologically meaningful manner.     

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.     

Analysis and modeling of the seasonal South China Sea temperature cycle using remote sensing

The present paper describes the analysis and modeling of the South China Sea (SCS) temperature cycle on a seasonal scale. It investigates the possibility to model this cycle in a consistent way while not taking into account tidal forcing and associated tidal mixing and exchange. This is motivated by the possibility to significantly increase the model’s computational efficiency when neglecting tides. The goal is to develop a flexible and efficient tool for seasonal scenario analysis and to generate transport boundary forcing for local models. Given the significant spatial extent of the SCS basin and the focus on seasonal time scales, synoptic remote sensing is an ideal tool in this analysis. Remote sensing is used to assess the seasonal temperature cycle to identify the relevant driving forces and is a valuable source of input data for modeling. Model simulations are performed using a three-dimensional baroclinic-reduced depth model, driven by monthly mean sea surface anomaly boundary forcing, monthly mean lateral temperature, and salinity forcing obtained from the World Ocean Atlas 2001 climatology, six hourly meteorological forcing from the European Center for Medium range Weather Forecasting ERA-40 dataset, and remotely sensed sea surface temperature (SST) data. A sensitivity analysis of model forcing and coefficients is performed. The model results are quantitatively assessed against climatological temperature profiles using a goodness-of-fit norm. In the deep regions, the model results are in good agreement with this validation data. In the shallow regions, discrepancies are found. To improve the agreement there, we apply a SST nudging method at the free water surface. This considerably improves the model’s vertical temperature representation in the shallow regions. Based on the model validation against climatological in situ and SST data, we conclude that the seasonal temperature cycle for the deep SCS basin can be represented to a good degree. For shallow regions, the absence of tidal mixing and exchange has a clear impact on the model’s temperature representation. This effect on the large-scale temperature cycle can be compensated to a good degree by SST nudging for diagnostic applications.     

Importance of cross-channel bathymetry and eddy viscosity parameterisation in modelling estuarine flow

For a proper understanding of flow patterns in curved tidal channels, quantification of contributions from individual physical mechanisms is essential. We study quantitatively how such contributions are affected by cross-channel bathymetry and three alternative eddy viscosity parameterisations. Two models are presented for this purpose, both describing flow in curved but otherwise prismatic channels with an (almost) arbitrary transverse bathymetry. One is a numerical model based on the full three-dimensional shallow water equations. Special feature of this diagnostic model is that assumptions regarding the relative importance of particular physical mechanisms can be incorporated in the computations by switching corresponding terms in the model equations on or off. We also present an idealized model that provides semi-analytical approximate solutions of the shallow water equations for all three considered alternative eddy viscosity parameterisations. It forms an aid in explaining and theorising about results obtained with the numerical model. Observations regarding Chesapeake Bay serve as a reference case for the present study. We find that the relative importance of both along-channel advective forcing and transverse diffusive forcing depends on local characteristics of the cross-sectional bottom profile rather than global ones. In our reference case, tide-residual along-channel flow induced by these forcings is not small compared to the total tidal residual. Building on this observation, we present an indicative test to judge whether advective processes should be included in leading order in modelling tide-dominated estuarine flow. Furthermore, depending on the applied eddy viscosity parameterisation (uniformly or parabolically distributed over the vertical), we find qualitatively different spatial patterns for the along-channel advective forcing.     

Far-field effects of tidal energy extraction in the Minas Passage on tidal circulation in the Bay of Fundy and Gulf of Maine using a nested-grid coastal circulation model

The Bay of Fundy in eastern Canada has the highest tides in the world. Harnessing the tidal energy in the region has long been considered. In this study, the effects of tidal in-stream energy extraction in the Minas Passage on the three-dimensional (3D) tidal circulation in the Bay of Fundy (BoF) and the Gulf of Maine (GoM) are examined using a nested-grid coastal ocean circulation model based on the Princeton Ocean Model (POM). The nested-grid model consists of a coarse-resolution (~4.5 km) parent sub-model for the GoM and a high-resolution (~1.5 km) child sub-model for the BoF. The tidal in-stream energy extraction in the model is parameterized in terms of nonlinear Rayleigh friction in the momentum equation. A suite of numerical experiments are conducted to determine the ranges of extractable tidal in-stream energy and resulting effects on the 3D tidal circulation over the Bay of Fundy and the Gulf of Maine (BoF-GoM) in terms of the Rayleigh friction coefficients. The 3D model results suggest that the maximum energy extraction in the Minas Passage increases tidal elevations and tidal currents throughout the GoM and reduces tidal elevations and circulation in the upper BoF, especially in the Minas Basin. The far-field effect of tidal energy extraction in the Passage on the 3D tidal circulation in the BoF-GoM is examined in two cases of harnessing tidal in-stream energy from (a) the entire water column and (b) the lower water column within 20 m above the bottom in the Passage. The 3D model results demonstrate that tidal in-stream energy extraction from the lower water column has less impact on the tidal elevations and circulation in the BoF-GoM than the energy extraction from the whole water column in the Minas Passage.     

Tidal effects on estuarine circulation and outflow plume in the Chesapeake Bay

The response of the Chesapeake Bay to river discharge under the influence and absence of tide is simulated with a numerical model. Four numerical experiments are examined: (1) response to river discharge only; (2) response to river discharge plus an ambient coastal current along the shelf outside the bay; (3) response to river discharge and tidal forcing; and (4) response to river discharge, tidal forcing, and ambient coastal current. The general salinity distribution in the four cases is similar to observations inside the bay. Observed features, such as low salinity in the western side of the bay, are consistent in model results. Also, a typical estuarine circulation with seaward current in the upper layer and landward current in the lower layer is obtained in the four cases. The two cases without tide produce stronger subtidal currents than the cases with tide owing to greater frictional effects in the cases with tide. Differences in salinity distributions among the four cases appear mostly outside the bay in terms of the outflow plume structure. The two cases without tide produce an upstream (as in a Kelvin wave sense) or northward branch of the outflow plume, while the cases with tide produce an expected downstream or southward plume. Increased friction in the cases with tide changes the vertical structure of outflow at the entrance to the bay and induces large horizontal variations in the exchange flow. Consequently, the outflow from the bay is more influenced by the bottom than in the cases without tide. Therefore, a tendency for a bottom-advected plume appears in the cases with tide, rather than a surface-advected plume, which develops in the cases without tide. Further analysis shows that the tidal current favors a salt balance between the horizontal and vertical advection of salinity around the plume and hinders the upstream expansion of the plume outside the bay.     

Predicting dissolved oxygen in the Chesapeake Bay: applications and implications

Eutrophic depletion of dissolved oxygen (DO) and its consequences for ecosystem dynamics have been a central theme of research, assessment and management policies for several decades in the Chesapeake Bay. Ongoing forecast efforts predict the extent of the summer hypoxic/anoxic area due to nutrient loads from the watershed. However, these models neither predict DO levels nor address the intricate interactions among various ecological processes. The prediction of spatially explicit DO levels in the Chesapeake Bay can eventually lead to a reliable depiction of the comprehensive ecological structure and functioning, and can also allow the quantification of the role of nutrient reduction strategies in water quality management. In this paper, we describe a three dimensional empirical model to predict DO levels in the Chesapeake Bay as a function of water temperature, salinity and dissolved nutrient concentrations (TDN and TDP). The residual analysis shows that predicted DO values compare well with observations. Nash–Sutcliffe efficiency (NSE) and root mean square error-observations standard deviation ratio (RSR) are used to evaluate the performance of the empirical model; the scores demonstrate the usability of model predictions (NSE, surface layer = 0.82–0.86; middle layer = 0.65–0.82; bottom layer = 0.70–0.82; RSR surface layer = 0.37–0.44; middle layer = 0.43–0.58 and bottom layer = 0.43–0.54). The predicted DO values and other physical outputs from downscaling of regional weather and climate predictions, or forecasts from hydrodynamic models, can be used to forecast various ecological components. Such forecasts would be useful for both recreational and commercial users of the Chesapeake Bay.     

Mixing,hypersalinity and gradients in Hervey Bay,Australia

Hervey Bay, a large coastal embayment situated off the central eastern coast of Australia, is a shallow tidal area (average depth = 15 m), close to the continental shelf. It shows features of an inverse estuary, due to the high evaporation rate (approx. 2 m/year), low precipitation (less than 1 m/year) and on average almost no freshwater input from rivers that drain into the bay. The hydro- and thermodynamical structures of Hervey Bay and their variability are presented here for the first time, using a combination of four-dimensional modelling and observations from field studies. The numerical studies are performed with the Coupled Hydrodynamical Ecological Model for Regional Shelf Seas (COHERENS). Due to the high tidal range (>3.5 m), the bay is considered as a vertically well-mixed system, and therefore, only horizontal fronts are likely. Recent field measurements, but also the numerical simulations, indicate characteristic features of an inverse/hypersaline estuary with low salinity (35.5 psu) in the open ocean and peak values (>39.0 psu) in the head water of the bay. The model further predicts a nearly persistent mean salinity gradient of 0.5 psu across the bay (with higher salinities close to the shore). The investigation further shows that air temperature, wind direction and tidal regime are mainly responsible for the stability of the inverse circulation and the strength of the salinity gradient across the bay. Due to an ongoing drying trend, the occurrence of severe droughts at the central east coast of Australia and, therefore, a reduction in freshwater supply, the salinity flux out of the bay has increased, and the inverse circulation has also strengthened.     

Analytical and numerical analysis of tides and salinities in estuaries; part II: salinity distributions in prismatic and convergent tidal channels

Estuaries, commonly, are densely populated areas serving the needs of the inhabitants in multiple ways. Often the interests are conflicting and decisions need to be made by the local managers. Intake of fresh water for consumption, agricultural purposes or use by industries may take place within a region not far landward of the limit of salt intrusion. Human interventions (e.g. deepening of the navigation channels) or climate changes (sea level rise, reduction of the river discharge) can bring these intake locations within the reach of saline or brackish water and consequently endanger their function. To support policy and managerial decisions, a profound knowledge of processes associated with the salinity structure in estuaries is required. Although nowadays advanced numerical three-dimensional models are available that are able to cope with the complexity of the physics there is still a need for relatively simple tools for quick-scan actions in a pre-phase of a project or for instructive purposes. The analytical model described in this paper may serve these needs. It computes the maximum salinity distribution using the dispersion coefficient in the mouth as the only model parameter. The model has been calibrated using observational data in a large number of estuaries and experimental data in a tidal flume. The dispersion coefficient was successfully related to geometric and hydrodynamic parameters resulting in an expression that can be used for convergent estuaries as well as prismatic channels, see Eqs. 25a and 25b. Application of the model in a predictive mode showed its promising capabilities. Comparison with three-dimensional numerical models indicates that the channel geometry in the estuary mouth largely influences dispersive processes. The analytical model for salt intrusion may be used in combination with the analytical model for tidal propagation in convergent estuaries and tidal channels by Van Rijn (part I). In this way, input is obtained on the tidal velocity amplitude and the Chézy roughness following calibration of this model on tidal amplitudes along the estuary.     

滨海地区基于避咸蓄淡模式的供水安全临界流量研究——以粤港澳大湾区珠澳供水为例

枯水期咸潮入侵已经严重威胁到了感潮河流区域供水安全.本文通过构建避咸蓄淡供水模型,耦合了咸度预测、河库联合供水调度和供水安全分析模块,为依赖感潮河流为水源地的区域供水安全管理提供了一种整体思路和决策方法.以面向粤港澳大湾区珠海东部及澳门的珠江三角洲磨刀门水道取供水为例,基于潮汐、径流和风等因子及咸度实测数据,较好地拟合了基于BP神经网络的咸度预测模型,及含氯度与超标时间的曲线函数,建立了上游来水和咸度超标时间之间的联系,得到了避咸蓄淡取水时机.咸度预测与当地河库联合供水调度相结合,得到了上游枯水期来水过程的当地区域供需平衡状况.枯水期不考虑水库调蓄的资源性缺水临界需水量为3.22亿m³,水库参与调蓄的工程性缺水临界需水量为3.75亿m³.通过供水安全分析模块,基于设定的风险阈值和临界阈值识别出了不同需水规模的上游来水临界流量特征.对于当地规划的需水规模4.23亿m³,期望上游枯水期临界流量均值约为3372 m³/s.整体上来说,需水规模越大,对上游来水期望的临界流量越大,但同时还与枯水期流量分布有关.     

Impact of historic land‐use change on sediment delivery to a Chesapeake Bay subestuarine delta

Historic land use in the Chesapeake Bay drainage basin induced large fluxes of fluvial sediment to subestuarine tributaries. Stratigraphic and palaeoecologic analyses of deltaic deposits may be used to infer changes on the landscape, but are not sufficient to quantify past sediment supply. When viewed as an inverse boundary‐value problem, reconstruction of the sediment supply function may be achieved by combining deltaic sedimentation chronologies with an equation governing delta progradation. We propose that the diffusion equation is appropriate for simulating delta progradation and obtaining the sediment supply function provided a suitable diffusion constant (D) can be determined. Three new methods for estimating D are presented for the case of estuarine deltas. When the inverse boundary‐value technique was applied to Otter Point Creek, a tidal freshwater delta at the head of Bush River in upper Chesapeake Bay, D values ranged from 3763 to 6199 m² a^?1. Delta growth simulations showed a 1740–1760 initial pulse, a 1760–1780 erosive/redistributive interval, a 1780–1920 growth period, and a 1920‐present erosive/redistributive era. Coupling of simulated delta elevations with an empirical plant habitat predictive equation allowed for comparison of predicted versus actual relative habitat areas. Also, the model yielded reconstructed watershed erosion rates and stream suspended sediment concentrations that could be useful for development of water quality regulations. Copyright © 2001 John Wiley & Sons, Ltd.     

Physical and dynamical oceanography of Liverpool Bay

The UK National Oceanography Centre has maintained an observatory in Liverpool Bay since August 2002. Over 8 years of observational measurements are used in conjunction with regional ocean modelling data to describe the physical and dynamical oceanography of Liverpool Bay and to validate the regional model, POLCOMS. Tidal dynamics and plume buoyancy govern the fate of the fresh water as it enters the sea, as well as the fate of its sediment, contaminants and nutrient loads. In this context, an overview and summary of Liverpool Bay tidal dynamics are presented. Freshwater forcing statistics are presented showing that on average the bay receives 233 m³ s^− 1. Though the region is salinity controlled, river input temperature is shown to significantly modulate the plume buoyancy with a seasonal cycle. Stratification strongly influences the region’s dynamics. Data from long-term moored instrumentation are used to analyse the stratification statistics that are representative of the region. It is shown that for 65% of tidal cycles, the region alternates between being vertically mixed and stratified. Plume dynamics are diagnosed from the model and are presented for the region. The spring–neap modulation of the plume’s westward extent, between 3.5 °W and 4°W, is highlighted. The rapid eastward erosion of the plume during spring tides is identified as a potentially important freshwater mixing mechanism. Novel climatological maps of temperature, salinity and density from the CTD surveys are presented and used to validate numerical simulations. The model is found to be sensitive to the freshwater forcing rates, temperature and salinities. The existing CTD survey grid is shown to not extend sufficiently near the coast to capture the near coastal and vertically mixed component the plume. Instead the survey grid captures the westward spreading, shallow and transient, portion of the plume. This transient plume feature is shown in both the long-term averaged model and observational data as a band of stratified fluid stretching between the mouth of the Mersey towards the Isle of Man. Finally the residual circulation is discussed. Long-term moored ADCP data are favourably compared with model data, showing the general northward flow of surface water and southward trajectory of bottom water.     

Wind effects on the lateral structure of density-driven circulation in Chesapeake Bay

The response of the density-driven circulation in the Chesapeake Bay to wind forcing was studied with numerical experiments. A model of the bay with realistic bathymetry was first applied to produce the density-driven flow under average river discharge and tidal forcing. Subsequently, four spatially uniform wind fields (northeasterly, northwesterly, southwesterly, and southeasterly) were imposed to examine the resulting cross-estuary structure of salinity and flow fields. In general, northeasterly and northwesterly winds intensified the density-driven circulation in the upper and middle reaches of the bay, whereas southeasterly and southwesterly winds weakened it. The response was different in the lower bay, where downwind flow from the upper and middle reaches of the bay competed with onshore/offshore coastal flows. Wind remote effects were dominant, over local effects, on volume transports through the bay entrance. However, local effects were more influential in establishing the sea-level slopes that drove subtidal flows and salinity fields in most of the bay. The effect of vertical stratification on wind-induced flows was also investigated by switching it off. The absence of stratification allowed development of Ekman layers that reached depths of the same order as the water depth. Consequently, bathymetric effects became influential on the homogeneous flow structure causing the wind-induced flow inside the bay to show a marked transverse structure: downwind over the shallow areas and upwind in the channels. In the presence of stratification, Ekman layers became shallower and the wind-induced currents showed weaker transverse structure than those that developed in the absence of stratification. In essence, the wind-driven flows were horizontally sheared under weak stratification and vertically sheared under stratified conditions.     

On the vertical structure of the Rhine region of freshwater influence

An idealised three-dimensional numerical model of the Rhine region of fresh water influence (ROFI) was set up to explore the effect of stratification on the vertical structure of the tidal currents. Prandle’s dynamic Ekman layer model, in the case of zero-depth-averaged, cross-shore velocities, was first used to validate the response of the numerical model in the case of barotropic tidal flow. Prandle’s model predicted rectilinear tidal currents with an ellipse veering of up to 2%. The behaviour of the Rhine ROFI in response to both a neap and a spring tide was then investigated. For the given numerical specifications, the Rhine plume region was well mixed over the vertical on spring tide and stratified on neap tide. During spring conditions, rectilinear tidal surface currents were found along the Dutch coast. In contrast, during neap conditions, significant cross-shore currents and tidal straining were observed. Prandle’s model predicted ellipse veering of 50%, and was found to be a good indicator of ellipticity magnitude as a function of bulk vertical eddy viscosity. The modelled tidal ellipses showed that surface currents rotated anti-cyclonically whereas bottom currents rotated cyclonically. This caused a semi-diurnal cross-shore velocity shearing which was 90° out of phase with the alongshore currents. This cross-shore shear subsequently acted on the horizontal density gradient in the plume, thereby causing a semi-diurnal stratification pattern, with maximum stratification around high water. The same behaviour was exhibited in simulations of a complete spring–neap tidal cycle. This showed a pattern of recurring stratification on neaps and de-stratification on springs, in accordance with observations collected from field campaigns in the 1990’s. To understand the increase in ellipticities to 30% during neaps and the precise shape of the vertical ellipse structure, stratification has to be taken into account. Here, a full three-dimensional numerical model was employed, and was found to represent the effect of de-coupling of the upper and lower layers due to a reduction of mixing at the pycnocline.     

Numerical model studies on the generation and dissipation of the diurnal coastal-trapped waves over the Sakhalin shelf in the Sea of Okhotsk

To clarify the generation and dissipation mechanisms of diurnal coastal-trapped waves (CTWs) over the Sakhalin shelf, a series of numerical experiments were conducted using a three-dimensional tidal model of the Okhotsk Sea with density stratification. The tidal model used has good reproduction owing to the careful fitting to the recent observations. The numerical experiments suggested that diurnal CTWs are primarily (~60%) generated by the conversion of tidal energy at the northern corner of the Sakhalin shelf, and further amplified by vorticity generation due to the water column oscillation from Sakhalin Bay and the influence of Kashevarov Bank. From the observations, it was found that diurnal CTWs are effectively dissipated by the strong spin-down due to bottom friction. The conventional turbulent closure model cannot reproduce the observed damping of diurnal CTWs, which raises a caution in modeling the tidal fields in high-latitude regions where diurnal CTWs exist. To resolve this underestimation of the damping, the vertical eddy viscosity was parameterized using its dependence on the observed major axis length of the diurnal tidal current ellipses, which improves the model reproduction on the damping of diurnal CTWs. The model also suggests that the spin-down effects due to friction associated with the sea-ice cover play an important role in the tidal current reduction in the region where diurnal CTWs exist, as the observations suggested.     

Simulation of submarine ground water discharge to a marine estuary: Biscayne Bay, Florida

Variable density ground water flow models are rarely used to estimate submarine ground water discharge because of limitations in computer speed, data availability, and availability of a simulation tool that can minimize numerical dispersion. This paper presents an application of the SEAWAT code, which is a combined version of MODFLOW and MT3D, to estimate rates of submarine ground water discharge to a coastal marine estuary. Discharge rates were estimated for Biscayne Bay, Florida, for the period from January 1989 to September 1998 using a three-dimensional, variable density ground water flow and transport model. Hydrologic stresses in the 10-layer model include recharge, evapotranspiration, ground water withdrawals from municipal wellfields, interactions with surface water (canals in urban areas and wetlands in the Everglades), boundary fluxes, and submarine ground water discharge to Biscayne Bay. The model was calibrated by matching ground water levels in monitoring wells, baseflow to canals, and the position of the 1995 salt water intrusion line. Results suggest that fresh submarine ground water discharge to Biscayne Bay may have exceeded surface water discharge during the 1989, 1990, and 1991 dry seasons, but the average discharge for the entire simulation period was only approximately 10% of the surface water discharge to the bay. Results from the model also suggest that tidal canals intercept fresh ground water that might otherwise have discharged directly to Biscayne Bay. This application demonstrates that regional scale variable density models are potentially useful tools for estimating rates of submarine ground water discharge.