Internal tide modelling and the influence of wind effects

Initially the development of shallow sea three-dimensional barotropic tidal models is briefly reviewed with a view to determining what were the key measurements that allowed progress in this field and rigorous model validation. Subsequently this is extended to a brief review of baroclinic tidal models to try to determine a “way forward” for baroclinic model development. The difficulty of high spatial variability, and wind influence are identified as possibly important issues that must be considered in validating baroclinic tidal models. These are examined using a three-dimensional unstructured grid model of the M₂ internal tide on the shelf edge region off the west coast of Scotland. The model is used to investigate the spatial variability of the M₂ internal tide, and associated turbulence energy and mixing in the region. Initial calculations are performed with tidal forcing only, with subsequent calculations briefly examining how the tidal distribution is modified by down-welling and up-welling favourable winds. Calculations with tidal forcing only, show that there is significant spatial variability in the internal tide and associated mixing in the region. In addition, these are influenced by wind effects which may have to be taken into account in any model validation exercise. The paper ends by discussing the comprehensive nature of data sets that need to be collected to validate internal tidal models to the same level currently attained with three dimensional barotropic tidal models.     

Modeling internal tides over Fieberling Guyot: resolution, parameterization, performance

Terrain-following ocean models are being used to simulate baroclinic tides and provide estimates of the tidal fields for circulation and mixing studies. These models have successfully reproduced elevations with most of the remaining inaccuracies attributed to topographic errors; however, the replication of barotropic and baroclinic velocity fields has not been as robust. Part of the problem is the lack of an adequate observational dataset in the simulated regions to compare the models. This problem was addressed using a dataset collected during the Flow over Abrupt Topography initiative at Fieberling Guyot. To evaluate the capability of the Regional Ocean Model System (ROMS) to simulate baroclinic tidal velocities, the combined tides for four constituents, M₂, S₂, K₁, and O₁, were modeled over Fieberling Guyot. Model inputs, numerical schemes, and parameterizations were varied to improve agreement with observations. These included hydrography, horizontal resolution, and the vertical mixing parameterization. Other factors were evaluated but are not included in this paper. With the best case, semidiurnal baroclinic tides were well replicated with RMS differences between the model estimates and the observations of 1.85 and 0.60 cm s⁻¹ for the major axes of the tidal ellipses for M₂ and S₂, respectively. However, diurnal K₁ baroclinic tides were poorly simulated with RMS differences of 4.49 cm s⁻¹. In the simulations, the K₁ baroclinic tides remained bottom-trapped unlike the observed fields, which had free waves due to the contribution of the mean velocity to the potential vorticity. The model did not adequately simulate the mean velocity, and the K₁ tides remained trapped. A resolution of 1 km most accurately reproduced the major axes and mean velocities; however, a 4-km resolution was sufficient for a qualitative estimate of where baroclinic tidal generation occurred. Nine vertical mixing parameterizations were compared. The vertical mixing parameterization was found to have minor effects on the velocity fields, with most effects occurring over the crown of guyot and in the lower water column; however, it had dramatic effects on the estimation of vertical diffusivity of temperature. Although there was no definitive best performer for the vertical mixing parameterization, several parameterizations could be eliminated based on comparison of the vertical diffusivity estimates with observations. The best performers were Mellor–Yamada and three generic length scale schemes.     

Intense diurnal surface currents in the Bay of La Paz,Mexico

Currents in the northern Bay of La Paz were examined using an 8-month Acoustic Doppler Current Profiler (ADCP) record collected in the upper 185 m of the water column during 2007. Flow variability was dominated by tidal motions, which accounted for 43% (33% diurnal, 10% semidiurnal) of the total kinetic energy. The tidal motions had a pronounced vertical structure dominated within a shallow (∼30 m thick) surface layer by intense counterclockwise (CCW) rotary S₁ diurnal radiational currents that were highly coherent with the counterclockwise seabreeze. Motions within the semidiurnal frequency band were primarily associated with significant counterclockwise S₂ radiational tidal currents, which were also coherent with the seabreeze. Both S₁ and S₂ tidal ellipses in the upper layer were aligned perpendicular to the bay entrance with mean semi-major axes of 55 and 20 cm/s, respectively. Below the surface layer, tidal currents decayed rapidly to relatively weak, clockwise rotary barotropic motions. In contrast to those for radiational harmonics, tidal ellipses of the gravitational constituents (M₂, K₁ and O₁) were oriented cross-bay. Energy within the diurnal frequency band in the surface layer was dominated by a coherent component (barotropic, phase-locked baroclinic and radiational), which accounted for roughly 65% (59% from S₁ alone) of the total diurnal kinetic energy. Of the remaining diurnal band energy, 18% was associated with an incoherent baroclinic component and 17% with a background noise component. Below 30 m depth, the corresponding estimates are 40%, 32% and 28%, respectively. The persistent, surface-intensified CCW rotary currents observed at the mooring site are assumed to be forced by strong CCW seabreeze winds in the presence of a “slippery” low-density surface layer. This response may be further augmented by topographic narrowing at the bay entrance and by the close proximity of the diurnal and inertial frequency bands in the region.     

Effect of stratification on tidal circulation over the Scotian Shelf and Gulf of St. Lawrence: a numerical study using a three-dimensional shelf circulation model

A three-dimensional shelf circulation model is used to examine the effect of seasonal changes in water-column stratification on the tidal circulation over the Scotian Shelf and Gulf of St. Lawrence. The model is driven by tidal forcing specified at the model’s lateral open boundaries in terms of tidal sea surface elevations and depth-averaged currents for five major tidal constituents (M₂, N₂, S₂, K₁, and O₁). Three numerical experiments are conducted to determine the influence of baroclinic pressure gradients and changes in vertical mixing, both associated with stratification, on the seasonal variation of tidal circulation over the study region. The model is initialized with climatological hydrographic fields and integrated for 16 months in each experiment. Model results from the last 12 months are analyzed to determine the dominant semidiurnal and diurnal tidal components, M₂ and K₁. Model results suggest that the seasonal variation in the water-column stratification affects the M₂ tidal circulation most strongly over the shelf break and over the deep waters off the Scotian Shelf (through the development of baroclinic pressure gradients) and along Northumberland Strait in the Gulf of St. Lawrence (through changes in vertical mixing and bottom stress). For the K₁ constituent, the baroclinic pressure gradient and vertical mixing have opposing effects on the tidal circulation over several areas of the study region, while near the bottom, vertical mixing appears to play only a small role in the tidal circulation.     

Observations of semidiurnal internal tidal currents off central Chile (36.6°S)

Observations of semidiurnal internal tidal currents from three moorings deployed on the continental shelf off central Chile during summer and winter of 2005 are reported. The spectra of the baroclinic currents showed large peaks at the semidiurnal band with a dominant counterclockwise rotation, which was consistent with internal wave activity. The amplitude of the barotropic tidal currents varied according to the spring–neap cycle following the sea level fluctuations. In contrast, the amplitudes of the internal tide showed high spatial-temporal variability not directly related to the spring–neap modulation. Near the middle of the continental shelf and near the coast (San Vicente Bay) the variance of the semidiurnal baroclinic current is larger than the variance of its barotropic counterpart. The vertical structure of the baroclinic tidal current fluctuations was similar to the structure of the first baroclinic internal wave mode. In general, in the three study sites the variance of the baroclinic current was larger near the surface and bottom and tended to show a minimum value at mid depths. Kinetic energy related to semidiurnal internal waves was larger in winter when stratification of the water column was stronger. During summer, upwelling and the decrease of freshwater input from nearby rivers reduced the vertical density stratification. The amplitude of the semidiurnal internal tide showed a tendency to be enhanced with increasing stratification as observed in other upwelling areas. The continental shelf break and submarine canyons, which limit the continental shelf in the alongshore direction, represent near-critical slopes for the semidiurnal period and are suggested to be the main internal tide generation sites in the study region.     

Tidal currents in the central Great Barrier Reef

Observations at 8 sites in the outer central Great Barrier Reef show M₂, S₂, K₁, and O₁ tidal currents flow directly off-shelf (northeast), when the corresponding tide at Townsville is at zero height and falling, with typical amplitudes of 12, 6, 3, and 2 cm s^?1. On the slope (at 300 m depth), the vertically averaged long-shelf component was small. On the shelf, the eccentricity of the tidal ellipses decreases shoreward and the tidal ellipses rotate anticlockwise. The major axes of the tidal ellipses tilt left of cross-shelf, especially for the diurnal constituents. There is satisfactory agreement between the observed and modelled cross-shelf currents. The long-shelf velocity is sensitive to the long-shelf changes in amplitude and phase of the tide heights and high quality tidal data for open boundary conditions will be required if numerical models are to model these currents satisfactorily.     

A numerical study of island wake generated by an elliptical tidal flow

An idealized numerical study of the influence of a tidal flow around an island has been undertaken with ROMS. The study focusses on coastal island wakes which are mainly controlled by elliptical tidal current flows on shallow shelves. This model is typical of some isolated continental shelf islands. The model is forced by a semi-diurnal barotropic inertia gravity wave imposed on the four open boundaries of a rectangular domain and its propagation results in an elliptical tidal flow within the domain in which the circular island lies. The influence of the surrounding island bathymetry and of the ellipse shape has been studied both in two and three dimensions. In the island vicinity, the residual circulation patterns over a tidal period show alongshore flow divergence along the major axis and convergence along the minor axis. A thin tidal ellipse (i.e. with a large ratio between major and minor axes) leads to strong eddy activity periods in the lee of the island during the flood and ebb phases, with eddy dissipation phases in between. By contrast, an almost round ellipse (axis ratio nearly 1) leads to vorticity filaments which continuously progress around the island without eddy shedding. The presence of a topographic slope in the vicinity of the island strengthens the eddy activity. This study suggests that the tidal current rotation favors the development of the eddy rotating in the same direction and weakens the development of the second eddy. In three dimensions with a surrounding bathymetry, an intense upwelling occurs in a large area in the lee of the island and the vertical velocities are stronger with thinner ellipses. With a flat bottom the vertical motions are almost fully generated by convergence and divergence of the secondary flow. With a varying bottom topography, the vertical motions come from a combination of this mechanism with convergence and divergence of the depth averaged flow.     

Variations in the schumann resonance polarization ellipse in the horizontal and vertical planes according to observations at the Barentsburg and Lovozero observatories

The diurnal variations in the amplitude and parameters of polarization ellipses of the first Schumann resonance according to three magnetic field components, observed on December 1–10, 2007, at the Barentsburg and Lovozero observatories, have been studied. Ellipses have been constructed in the (H, D) and (H, Z) planes. The value of the minor axis, inclination of the major axis, ellipticity, and the rotation direction have been estimated. The vertical magnetic component of the Schumann resonance is three to four times as small as the horizontal component. The difference in the diurnal variations in the ellipse parameters between both observatories has been found. The effect can be caused by a difference in the Earth’s conductivity in the vicinity of the observatories and by conductivity anisotropy. The major axis inclination and ellipticity have semidiurnal components. The polarization vector mostly rotates from D to H in the horizontal plane and from Z to H in the vertical plane at both stations.     

Changes in the dynamics/energetics of surface and internal tides in the White Sea on the coexistence of shore-fast and drifting ice covers

The results of modeling for _M2${\mathrm{M}}_2$ surface and internal tides in the White Sea are discussed. These results are obtained for the case when shore-fast and drifting ice covers are present concurrently. It is assumed that the interface between ice covers is of non-tidal origin (i.e., it is pre-assigned) and that ice rheology is viscous-elastic, representative of the low temperatures typical of winter conditions. Emphasis is placed on tidal energetics and, in particular, on the averaged (over a tidal cycle) values of the density and the dissipation rate of barotropic/baroclinic tidal energy. It is shown that in the White Sea, unlike in other marginal seas, the averaged (over a tidal cycle) and depth-integrated density of baroclinic tidal energy for the combined ice cover is much less than the same defined density of barotropic tidal energy. Similarly, the averaged and integrated (over the volume of the White Sea) rate of baroclinic tidal energy dissipation is much less than the same defined rate of barotropic tidal energy dissipation. The latter, in turn, is greater than for the shore-fast ice cover, but is smaller than for the drifting ice cover.     

Variability of baroclinic tidal currents on the Mackenzie Shelf,the Southeastern Beaufort Sea

Semidiurnal tidal currents on the outer shelf of the Mackenzie Shelf in the Beaufort Sea were found to be strongly influenced by the locally generated baroclinic tide. Two primary factors are involved in this process: (1) the sharp shelf break along the northeastern Mackenzie Shelf, promoting the generation of vigorous internal tidal waves; and (2) the proximity to critical latitudes for M₂ and N₂ motions locking these waves and preventing them from leaving the source region. As a result, internal tides are resonantly trapped between the shelf and critical latitudes. The physical properties and temporal variations of tidal motions were examined using current meter measurements obtained from 1987–1988 at four sites (SS1, SS2, SS3, and SS4) offshore of the shelf break at depths of ∼200 m. Each mooring had Aanderaa RCM4s positioned at ∼35 m below the surface and ∼50 m above the bottom. Complex demodulation was used to compute the envelopes (amplitude modulation) of these components. A striking difference in the variability of clockwise (CW) and counterclockwise (CCW) tidal currents was found. The CW tides are highly variable, have greater amplitude, exhibit a burst-like character associated with wind events and contain about 80% of the total energy of the semidiurnal tidal currents. In contrast, the CCW components have a more regular temporal regime with distinct monthly, fortnightly and 10-day modulation at astronomical periodicities associated with frequency differences M₂–N₂ (0.03629 cpd), S₂–M₂ (0.06773 cpd), and S₂–N₂ (0.10402 cpd). Significant horizontal correlation of the CW current envelopes was found only between stations near the northeast Mackenzie Shelf, indicating this to be the main area of baroclinic internal wave generation.     

Tidal and wind-driven surface currents in the German Bight: HFR observations versus model simulations

Tidal and wind-driven surface currents in the German Bight between shallow mudflats of the North Frisian islands and the island of Helgoland are studied using coastal high-frequency radar (HFR) observations and hindcasts from a primitive equation numerical model. The setup of the observational system is described, and estimates of expected measurement errors are given. A quantitative comparison of numerical model results and observations is performed. The dominant tidal components are extracted from the two data sources using tidal harmonic analysis and the corresponding tidal ellipses are defined. Results show that the spatial patterns of different tidal ellipse parameters are consistent in the two data sets. Model sensitivity studies with constant and variable salinity and temperature distributions are used to study density-related mechanisms of circulation. Furthermore, the role of the surface wind field in driving the German Bight circulation is investigated using the complex correlation between wind and surface current vectors. The observed change of the respective correlation patterns from the coastal to open ocean is shown to be due to a combination of density effects, the coastline and topography. The overall conclusion is that HFR observations resolve the small-scale and rapidly evolving characteristics of coastal currents well in the studied area and could present an important component for regional operational oceanography when combined with numerical modelling. Some unresolved issues associated with the complex circulation and large instability of circulation in front of the Elbe River Estuary justify further considerations of this area using dedicated surveys and modelling efforts.     

Modeling the influence of river discharge on salt intrusion and residual circulation in Danshuei River estuary,Taiwan

A 3-D, time-dependent, baroclinic, hydrodynamic and salinity model was implemented and applied to the Danshuei River estuarine system and the adjacent coastal sea in Taiwan. The model forcing functions consist of tidal elevations along the open boundaries and freshwater inflows from the main stream and major tributaries in the Danshuei River estuarine system. The bottom friction coefficient was adjusted to achieve model calibration and verification in model simulations of barotropic and baroclinic flows. The turbulent diffusivities were ascertained through comparison of simulated salinity time series with observations. The model simulation results are in qualitative agreement with the available field data.     

Numerical investigations of the turbulent kinetic energy dissipation rate in the Rhine region of freshwater influence

The turbulent kinetic energy dissipation rate, ε, in tidal seas is maximum at the bottom during full flood and during full ebb, i.e. when tidal currents are strongest. In coastal regions with tides similar to a Kelvin wave, this coincides with high water and low water. If there is a freshwater source at the coast, stratification in such a region will be most stable at high water and least at low water. Measurements of ε in the Rhine region of freshwater influence performed by previous studies have revealed bottom maxima at both high and low water. In addition, a maximum in the upper half of the water column was found around high water, which cannot be explained by tidal shear at the bottom, convective instabilities or wind mixing. This study investigates the dissipation rate and relevant physical properties in the Rhine region of freshwater influence by means of three-dimensional numerical simulations using the General Estuarine Transport Model and idealised conditions. The measurements are well reproduced; two distinct peaks of ε are evident in the upper layer shortly before and after high water. These maxima turn out to be due to strong peaks in the alongshore shear occurring when the fore- and the back-front of the plume transit the water column.     

An experimental study of global instabilities due to the tidal (elliptical) distortion of a rotating elastic cylinder

Abstract

Broad band secondary instability of elliptical vortex motion has been proposed as a principal source of shear-flow turbulence. Here experiments on such instability in an elliptical flow with no shear boundary layer are described. This is made possible by the mechanical distortion in the laboratory frame of a rotating fluid-filled elastic cylinder. One percent ellipticity of a 10 cm diameter cylinder rotating once each second can give rise to an exponentially-growing mode stationary in the laboratory frame. In first order this mode is a sub-harmonic parametric Faraday instability. The finite-amplitude equations represent angular momentum transfer on an inertial time scale due to Reynolds stresses. The growth of this mode is not limited by boundary friction but by detuning and centrifugal stabilization. On average, a generalized Richardson number achieves a marginal value through much of the evolved flow. However, the characteristic flow is intermittent with the cycle: rapid growth, stabilizing momentum transfer from the mean flow, interior re-spin up, and then again. Data is presented in which, at large Reynolds numbers, seven percent ellipticity causes a fifty percent reduction in the kinetic energy of the rotating fluid. In the geophysical setting, this tidal instability in the earth's interior could be inhibited by sub-adiabatic temperature gradients. A near adiabatic region greater than 10 km in height would permit the growth of tidally destabilized modes and the release of energy to three-dimensional disturbances. Such disturbances might play a central role in the geodynamo and add significantly to overall tidal dissipation.     

Strong tidal currents over two shallow banks in the western subarctic Pacific

The south-flowing waters of the Kamchatka and Oyashio currents are key components of the western subarctic Pacific gyre. The dissipation of tidal energy in shallow and coastal regions of these currents and the attendant mixing are the important processes that affect the upper layer temperature and salinity. Examples of the impact of tidal currents on water temperature and salinity are the persistent tide-driven mixing around the Kashevarov and Kruzenshtern banks. The Kruzenshtern Bank is a shallow submarine bank stretching along the eastern continental slope of the Kuril Islands with the minimum depth of 86 m. Surface drifters observations are used to determine the characteristics of tidal currents and the circulation over these banks. New software that allows more versatility in the harmonic analysis is used for drifter’s data. The two banks have similar features. The variations in current velocities are dominated by the diurnal signals. The K₁ and O₁ tidal ellipses over the banks are the largest and clockwise. The enhanced tidal currents suggest that the formation of cold and saline water in summer is due to mixing of water column over the banks with intermediate waters. Variations of tidal ellipses over the bank may explain the formation of polynya at the western end of the Kashevarov Bank. We found that the 18.6-year lunar nodal cycle is a significant characteristic of salinity variation at the temperature minimum in the eddies eastward of the Boussole Strait over the period 1990–2015.     

ANISOTROPIC TRAVELTIME TOMOGRAPHY

Velocity estimation technique using seismic data is often based on time/distance equations which are independent of direction, and even though we now know that many rocks are quite anisotropic, useful results have been obtained over the years from these isotropic estimates. Nevertheless, if velocities are significantly direction-dependent, then the isotropic assumption may lead to serious structural interpretation errors. Additionally, information on angle-dependence may lead to a better understanding of the lithology of the rocks under measurement. VSP and cross-well data may each lack the necessary aperture to estimate more than two velocity parameters for each wave type, and if the data straddle a symmetry axis, then these may be usefully chosen to be the direct velocities (from time-and-distance measurements along the axis) and NMO velocities (from differential time-offset measurements). These sets of two parameters define ellipses, and provide intermediate models for the variation of velocity with angle which can later be assembled and translated into estimates of the elastic moduli of the rocks under scrutiny. If the aperture of the measurements is large enough e.g. we have access to both VSP and cross-well data, we divide the procedure into two independent steps, first fitting best ellipses around one symmetry axis and then fitting another set around the orthogonal axis. These sets of four elliptical parameters are then combined into a new, double elliptical approximation. This approximation keeps the useful properties of elliptical anisotropy, in particular the simple relation between group and phase velocities which simplifies the route from the traveltimes measurements to the elastic constants of the medium. The inversion proposed in this paper is a simple extension of well-known isotropic schemes and it is conceptually identical for all wave types. Examples are shown to illustrate the application of the technique to cross-well synthetic and field P-wave data. The examples demonstrate three important points: (a) When velocity anisotropy is estimated by iterative techniques such as conjugate gradients, early termination of the iterations may produce artificial anisotropy. (b) Different components of the velocity are subject to different type of artifacts because of differences in ray coverage, (c) Even though most rocks do not have elliptical dispersion relations, our elliptical schemes represent a useful intermediate step in the full characterization of the elastic properties.     

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.     

Tidal wave transformations in the German Bight

Mesoscale and submesoscale dynamics associated with tidal wave transformations were addressed in the German Bight using numerical simulations. Tidal gauge and velocity observations in several locations were used to validate the numerical model. A downscaling approach included analysis of simulations with horizontal resolutions of 1, 0.4, and 0.2 km. It was shown that the modified tidal wave lost most of its energy after reflection or refraction over the eastern part of the German Bight. Energy loss resulted in a pronounced change of the wave’s spectral composition and generation of overtides. Tidal oscillations were modified by mesoscale processes associated with bathymetric channels. Semidiurnal and quarterdiurnal tides revealed very different spatial patterns. The former were aligned with the bathymetric channels, while the latter were rather “patchy” and had about half the spatial scales. In numerous areas around the bathymetric channels, the major axis of the M₄ ellipses was normal or at some angle with the major axis of the M₂ ellipses. Thus, higher harmonics developed “orthogonal” patterns that drove secondary circulations. Moreover, the ratio between spring and neap tidal amplitudes was relatively low in the Wadden Sea, showing reduced sensitivity of this very shallow area to fortnightly tidal variations. It was demonstrated that simulated hydrodynamics patterns help explain the physical mechanism shaping the median grain size distribution in the German Bight.     

Three-dimensional modelling of estuarine turbidity maxima in a tidal estuary

A new numerical model for simulating estuarine dynamics is introduced here. This model, called General Estuarine Transport Model (GETM), has been specifically designed for reproducing baroclinic, bathymetry-guided flows where the tidal range may exceed the mean water depth in large parts of the domain such that drying and flooding processes are relevant. Several physical and numerical features of the model support exact and stable results for such domains. For the physics, high-order turbulence closure schemes guarantee proper reproduction of vertical exchange processes. Among the specific numerical features, generalised vertical coordinates, orthogonal curvilinear horizontal coordinates, high-order TVD advection schemes and stable drying and flooding algorithms have been implemented into GETM. The model is applied here to simulate the dynamics of estuarine turbidity maxima (ETMs), a complex feature present in most tidal estuaries. First, idealised simulations for a two-dimensional domain in the xz space will be shown to reproduce the basic generation mechanisms for ETMs. Then, a realistic three-dimensional simulation of the Elbe estuary in Northern Germany will be carried out. It is demonstrated that for a given forcing situation the model reproduces a stable ETM at the correct location.Responsible Editor: Phil Dyke     

Tidal-induced buoyancy flux and mean transverse circulation

A weakly nonlinear model is used to examine the mean transverse circulation (cross-isobath) driven by tidal-induced buoyancy flux. The mean Eulerian flows driven by both the barotropic and baroclinic tide are presented for a semi-infinite wedge. The mean flow driven by the barotropic tide is significant only near the apex where the thickness of the frictional boundary layer is comparable to the water depth. The mean flow there is characterized by a single-cell circulation with offshore flow near the bottom, and its magnitude can reach a few percentage or a significant fraction of the tidal velocity in oceanic applications. The mean flow driven by the baroclinic tide, on the other hand, is characterized by pairs of half-open (on the seaward side) counter-rotating cells, the number of which equals the vertical mode number. For a baroclinic tide propagating onshore, the mean flow near the top and bottom surfaces is always directed offshore and its magnitude can reach a large fraction of the tidal velocity. Taken together, the model thus predicts a mean offshore flow near the bottom while higher up in the water column the mean flow direction is less definite due to the contribution from different tidal components. The model results are consistent with some current measurements over the Georges Bank.