风驱动下f-平面准地转三维海洋环流的形成及总质量守恒的应用

将理想化的南中国海海盆在垂直方向上划分为Ekman层、惯性层和摩擦层. Ekman层中的运动由大气风应力驱动,其底部的扰动压力将作为其下惯性层中运动的上边界条件. 惯性层中的运动是由f 平面三维非线性方程在准地转近似下位势涡度守恒控制,由此得到控制惯性层中运动关于扰动压力的三维椭圆型方程. 在惯性层以下考虑到深层的海盆水平尺度很小,由此引进带有底部摩擦的线性控制方程,方程的边界条件为惯性层和摩擦层交界面上的扰动压力连续,沿海盆边界假定海水与相邻的固壁间无热量交换,由此设在海盆边界上扰动温度为零. 在此基础上分别利用惯性层和摩擦层中的椭圆型控制方程计算了相应层次上冬、夏季的扰动压力和准地转流. 结果表明冬季各层上以气旋式环流为主,且随深度的增加流速减小;夏季各层上以反气旋式环流为主,流速也随深度增加而减小. 这在一定程度上和观测事实相符.     

Upwelling and surface cold patches in the Yellow Sea in summer: Effects of tidal mixing on the vertical circulation

A three-dimensional, prognostic, wave–tide–circulation coupled numerical model is developed to study the effects of tidal mixing on the summertime vertical circulation in the Yellow Sea (YS). The distribution and mechanisms of upwelling are investigated by numerical means. Validated by historical tide gauge data, satellite sea surface temperature (SST) data, and cruise observation data, the model shows satisfactory performances in reproducing the dominant tidal system and three-dimensional sea temperature structure. Model results suggest that strong tidal mixing plays an important role in the formation of the vertical circulation in the YS. The Yellow Sea Cold Water Mass (YSCWM) is fringed by typical tidal mixing fronts (TMFs), which separate the cold, stratified water at the offshore side from the warm, well-mixed, shallow water at the other side. Considerable baroclinic gradient across the TMF makes the frontal zone the spot where the most active vertical circulation occurs; a secondary circulation is triggered with a distinct upwelling branch occurring mainly on the mixed side of the front. The numerical model produces systematic upwelling belts surrounding the YSCWM, and the upwelling is essentially induced by the TMF over sloping topography. The relative importance of tidal mixing and wind forcing for upwelling is further examined in numerical experiments. The southerly wind enhances the upwelling off the western coasts, but its overall influences in the whole YS are less important than tidal mixing. As shown by both satellite data and numerical modeling, the summertime SST field in the YS is featured by the stable existence of several site-selective surface cold patches (SCPs), most of which scatter in the waters off convex coastlines. One of the SCPs is found off Subei Bank, and the others are located off the eastern tip of Shandong Peninsula and off the three tips of Korean Peninsula. Two processes give rise to the SCP: on the one hand, TMF-induced upwelling supplies cold water from the deep layer; on the other hand, tidal mixing itself can stir the bottom water upward and homogenize the water column vertically. In the waters around the tips of peninsula in the YS, the tidal currents are extraordinarily strong, which provides a possible explanation for the site-selectivity of the SCPs.     

Sediment dynamics in an offshore tidal channel in the southern Yellow Sea

The geomorphology of the southern Yellow Sea（SYS） is characterized by offshore radial sand ridges（RSR）.An offshore tidal channel（KSY Channel） is located perpendicular to the coast,comprised of a main and a tributary channel separated by a submarine sand ridge（KSY Sand Ridge） extending seaward.In order to investigate the interactions among water flow,sediment transport,and topography,current velocity and suspended sediment concentration（SSC） were observed at 11 anchor stations along KSY Channel in RSR during a spring tide cycle.High resolution bottom topography was also surveyed.Residual currents and tidally averaged suspended sediment fluxes were calculated and analyzed by using the decomposition method.Results suggested that the water currents became stronger landward but with asymmetrical current speed and temporal duration of flood and ebb tides.Residual currents showed landward water transport in the nearshore channel and a clockwise circulation around the KSY Sand Ridge.Tidally-averaged SSC also increased landward along the channel.The main mechanisms controlling SSC variations were resuspension and horizontal advection,with spatial and temporal variations in the channel,which also contributed to sediment redistribution between channels and sand ridges.Residual flow transport and the tidal pumping effect dominated the suspended sediment flux in the KSY Channel.The KSY Sand Ridge had a potential southward migration due to the interaction between water flow,sediment transport,and topography.     

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.     

Influence of bottom frictional effects in sill regions upon lee wave generation and implications for internal mixing

A cross-sectional nonhydrostatic model using idealized sill topography is used to examine the influence of bottom friction upon unsteady lee wave generation and flow in the region of a sill. The implications of changes in shear and lee wave intensity in terms of local mixing are also considered. Motion is induced by a barotropic tidal flow which produces a hydraulic transition, associated with which are convective overturning cells, wave breaking, and unsteady lee waves that give rise to mixing on the lee side of the sill. Calculations show that, as bottom friction is increased, current profiles on the shallow sill crest develop a highly sheared bottom boundary layer. This enhanced current shear changes the downwelling of isotherms downstream of the sill with an associated increase in the hydraulic transition, wave breaking, and convective mixing in the upper part of the water column. Both short and longer time calculations with wide and narrow sills for a number of sill depths and buoyancy frequencies confirm that increasing bottom friction modifies the flow and unsteady lee wave distribution on the downstream side of a sill. Associated with this increase in bottom friction coefficient, there is increased mixing in the upper part of the water column with an associated decrease in the vertical temperature gradient. However, this increase in mixing and decrease in temperature gradient in the upper part of the water column is very different from the conventional change in near-bed temperature gradient produced by increased bottom mixing that occurs in shallow sea regions as the bottom drag coefficient is increased.     

The baroclinic circulation structure of Yellow Sea Cold Water Mass

The Yellow Sea is a semi-enclosed shallow sea with a deep trough of about 80 m. On the hy-drographic condition in the Yellow Sea, Lie[1] pointed out that it is strongly associated with winter cooling and summer heating, fresh input from rivers into the co…     

Circulation induced by subglacial discharge in glacial fjords: Results from idealized numerical simulations

The flow caused by the discharge of freshwater underneath a glacier into an idealized fjord is simulated with a 2D non-hydrostatic model. As the freshwater leaves horizontally the subglacial opening into a fjord of uniformly denser water it spreads along the bottom as a jet, until buoyancy forces it to rise. During the initial rising phase, the plume meanders into complex flow patterns while mixing with the surrounding fluid until it reaches the surface and then spreads horizontally as a surface seaward flowing plume of brackish water. The process induces an estuarine-like circulation. Once steady-state is reached, the flow consists of an almost undiluted buoyant plume rising straight along the face of the glacier that turns into a horizontal surface layer thickening as it flows seaward. Over the range of parameters examined, the estuarine circulation is dynamically unstable with gradient Richardson number at the sheared interface having values of <1/4. The surface velocity and dilution factors are strongly and non-linearly related to the Froude number. It is the buoyancy flux that primarily controls the resulting circulation with the momentum flux playing a secondary role.     

A robust well-balanced finite volume model for shallow water flows with wetting and drying over irregular terrain

An unstructured Godunov-type finite volume model is developed for the numerical simulation of geometrically challenging two-dimensional shallow water flows with wetting and drying over convoluted topography. In the framework of sloping bottom model, a modified formulation of shallow water equations is used to preserve mass conservation during flooding and recession. The key ingredient of the model is the use of this combination of the sloping bottom model and the modified shallow water equations to provide a robust technique for wet/dry fronts tracking and, together with centered discretization of the bed slope source term, to exactly preserve the static flow on irregular topographies. The variable reconstruction technique ensures nonnegative reconstructed water depth and reasonable reconstructed velocity, and the friction terms are solved by semi-implicit scheme that does not invert the direction of velocity components. The robustness and accuracy of the proposed model are assessed by comparing numerical and reference results of extensive test cases. Moreover, the results of a dam-break flooding over real topography are presented to show the capability of the model on field-scale application.     

Observation of the seasonal evolution of the Yellow Sea Cold Water Mass in 1996–1998

We use the hydrographic data obtained during the joint survey of the Yellow Sea by the First Institute of Oceanography, China and the Korea Ocean Research and Development Institute, Korea, to quantify the spatial structures and temporal evolution of the southern Yellow Sea Cold Water Mass (YSCWM). It is indicated that the southern YSCWM is a water mass that develops in summer and decays in fall. In winter, due to the intrusion of the Yellow Sea Warm Current (YSWC), the central area (approximately between 34°N and 35°N, 122°E and 124°E) of the Yellow Sea is mainly occupied by relatively high temperature water (T>10 °C). By contrast, from early summer to fall, under the seasonal thermocline, the central area of Yellow Sea is occupied by cold water (T<10 °C). In summer, the southern YSCWM has two cold cores. One is formed locally southeast of Shandong Peninsula, and the other one has a tongue-like feature occupying the area approximately between 34°N and 37°N, 123°E and 126°E. The bottom layer temperature anomalies from February to July in the cold tongue region, along with the trajectories of the bottom floaters, suggest that the cold water mass in the northeast region has a displacement from the north to the central area of the Yellow Sea during the summer.     

A three-dimensional numerical shelf-sea front model with variable eddy viscosity and diffusivity

A three-dimensional numerical model of circulation and eddy development in shelf-sea fronts is applied to three frontal structures, with two parameterization schemes for vertical eddy viscosity and diffusivity. The three fronts resemble those in the German Bight (a front between relatively fresh coastal water and saltier water offshore, with an interface extending from surface to bottom), the Norwegian Coastal Current (also formed by fresh coastal water but with a thermocline on one side), and the Celtic Sea (a front between water which is stratified in summer and water which is well mixed throughout the year). The two mixing assumptions, modelling the reduction of turbulence in stratified zones, are based on the Munk-Anderson scheme and the turbulent energy equation. Many features of frontal dynamics are common to all the results: strong surface currents along the front, cross-frontal circulation cells, a considerable enhancement of vertical velocities when eddies are formed, and development of eddies into cyclonic-anticyclonic vortex pairs. Cross-frontal circulation and frontal sharpening are the variables most sensitive to the different mixing assumptions. The German Bight front is the one most affected by changing these assumptions. The comparisons suggest that realistic results may be obtained from models despite the present uncertainty about vertical mixing in stratified shelf seas.     

Analytical modelling of wintertime coastal jets in the Adriatic Sea

Two diagnostic models, reproducing circulation generated in a marginal sea by variable density, have been developed. The models’ domain is a 2D transverse section for which analytical solutions have been obtained. They describe the winter situation in the northern Adriatic, with a strong vertical mixing present and the density maximum dominating the centre of the basin. Both models employ Boussinesq-type parametrisation of friction and linear slip at the bottom. The first model allows for frictional departure from hydrostatic equilibrium and includes vertical friction only. The second one is hydrostatic but allows for lateral friction as well. The results obtained by the two models are similar and to some extent dependent on the vertical and bottom friction. They reproduce several well known characteristics of the Adriatic circulation (cyclonic surface flow, downwelling in the central and larger part of the basin compensated by upwelling in the coastal zone) but also predict some phenomena that are still not well understood. A conspicuous feature of the model results are coastal jets, which were observed in the Adriatic on several occasions. The present models show that the distance of jets from the coasts depends on lateral friction: it is found to vary from 1 up to 10 km on the Italian side and between 2 and 15 km on the Croatian side. Both models reproduce the west–east asymmetry, with the wider current on the east side of the basin. The asymmetry is a subject on which conflicting empirical results exist in the Adriatic. In the two models cyclonic flow occupies the whole water column, which disagrees with some recent theoretical findings of the near-bottom anticyclonic flow and thus leaves the issue open.     

Modelling the Gibraltar Strait/Western Alboran Sea ecohydrodynamics

The ecohydrodynamics of the Gibraltar Strait and the Western Alboran Sea is investigated using a 3-D, two-way nested, coupled hydrodynamic/plankton ecosystem model, exploiting the MEDATLAS climatological database. A high-resolution model (~1 km) of the Gibraltar/Western Alboran region embedded within a coarse-resolution model of the West Mediterranean (~5 km) is implemented. The model seasonal climatology of the 3-D circulation and the flow characteristics at the Gibraltar Strait and the Alboran Sea are discussed, and their impact on the plankton ecosystem evolution is explored. An important ecohydrodynamic feature produced by the model is a permanent upwelling zone in the northwestern part of the Alboran Sea in agreement with observations. Model results show that both horizontal and vertical current intensity of the Atlantic Jet increases progressively at the strait to obtain maximum values in the northeastern Mediterranean entrance, inducing an upward displacement of the nitracline. The nutrient-rich water transport through the strait along with the generation of cyclonic vorticity in the northwestern Alboran Sea result in the accumulation of nutrients there and thus induce a permanent fertilisation of this area.     

Dense water generation on a shelf: the case of the Adriatic Sea

The paper overviews recent and past studies of preconditioning, generation and spreading of North Adriatic Dense Water (NAdDW), by analysing both transient episodes and climatological data. The importance of wind stress, heat and water fluxes, and particularly river discharges during the preconditioning and generative period is emphasized, as well as the advection of saline levantine intermediate water from the southeast. After the generation, NAdDW affects deep and bottom layers of the middle and south Adriatic Sea, flowing as a dense current and mixing with the adjacent waters; it can be traced even in the Otranto Strait, contributing to the formation of deep water in the Eastern Mediterranean. Objective shortcomings and the projections in NAdDW investigations are also discussed, due to their high importance in the circulation of the Adriatic Sea, which may relate other similar basins around the world.     

Recent changes of the thermohaline circulation in the subpolar North Atlantic

Time series of hydrographic sections in the northern North Atlantic from the period 1990 to 2004 are analyzed for changes in the characteristics and distribution of water masses that are involved in the thermohaline circulation (THC). During the 1990s, the North Atlantic Oscillation (NAO) alternates from a positive phase (strong westerlies) to a negative phase (weak westerlies). The reduced ocean heat loss confined the convection in the Labrador Sea to the upper 1,200 m, generating a new salinity minimum layer characterizing the Upper Labrador Sea Water (ULSW), and led to a warming and salinization of the older LSW below due to lateral mixing. The Lower LSW, formed in the first half of the 1990s, spread in the subpolar gyre and reached the Newfoundland and Irminger basins after about 1 to 2 years, where the associated isopycnal doming contributed to eastward frontal shifts in the upper layer. After 5 and 6 years, it arrived in the Iceland and West European basins, respectively. The collapse of the isopycnal dome in the Labrador Sea, associated with the drainage of the Lower LSW, resulted in a slowing of the cyclonic circulation of the subpolar gyre. This was accompanied in the upper layer by a westward shift of the southeastern extension of the gyre and a northward advection of warm and saline subtropical water in its eastern part, which finally reached the Labrador Sea after about 7 years. In the upper layer of the Labrador Sea, the advection of warm and saline water dominated over the heat loss to the atmosphere and the freshwater gain from melting ice and precipitation in the NAO-low period, so that no accumulation of freshwater but an increase of the heat and salt contents were observed, as in the whole eastern part of the subpolar gyre. Within 1 to 2 years after the drop of the NAO in the winter of 1995/1996, the Subarctic (Subpolar) Front shifted northward and westward north of about 50°N, favored by the retreat of the low-salinity tongue extending eastward from the southern Labrador Sea, and it shifted southward and eastward in the Newfoundland Basin. Therefore, the enhanced northward advection of subtropical waters in the northeastern North Atlantic is balanced by the enhanced southward advection of subarctic waters, including Lower LSW in the Newfoundland Basin, indicating a strong response of the gyre component of the THC.     

Pb in the western Indian Ocean: distribution, disequilibrium, and partitioning between dissolved and particulate phases

²²⁶Ra profiles have been measured in the western Indian Ocean as part of the 1977–1978 Indian Ocean GEOSECS program. These profiles show a general increase in deep and bottom water Ra concentration from the Circumpolar region to the Arabian Sea. A deep Ra maximum which originates in the Arabian Sea and in the Somali basin at about 3000 m depth spreads southward into the Mascarene basin and remains discernible in the Madagascar and Crozet basins. In the western Indian Ocean, the cold Antarctic Bottom Water spreads northward under the possibly southward-flowing deep water, forming a clear benthic front along the Crozet basin across the Southwest Indian Ridge into the Madagascar and Mascarene basins. The Antarctic Bottom Water continues to spread farther north to the Somali basin through the Amirante Passage at 10°S as a western boundary current. The benthic front and other characteristic features in the western Indian Ocean are quite similar to those observed in the western Pacific where the benthic front as a distinctive feature was first described by Craig et al. [15]. Across the Mid-Indian Ridge toward the Ceylon abyssal plain near the triple junction, Ra profiles display a layered structure, reflecting the topographic effect of the mid-ocean ridge system on the mixing and circulation of the deep and bottom waters. Both Ra and Si show a deep maximum north of the Madagascar basin. Linear relationships between these two elements are observed in the deep and bottom water with slopes increasing northward. This suggests a preferential input of Ra over Si from the bottom sediments of the Arabian Sea and also from the flank sediments of the Somali basin.     

On the influence of adequate Weddell Sea characteristics in a large-scale global ocean circulation model

Global ocean circulation models usually lack an adequate consideration of high-latitude processes due to a limited model domain or insufficient resolution. Without the processes in key areas of the global thermohaline circulation, the characteristics and flow of deep and bottom waters cannot be modeled realistically. In this study, a high-resolution (~20 km) ocean model focused on the Weddell Sea sector of the Southern Ocean is combined with a low-resolution (2^° × 2^°) global ocean model applying the state estimation technique. Temperature, salinity, and velocity data on two Weddell Sea sections from the regional model are used as constraints for the large-scale model in addition to satellite altimetry and sea-surface temperatures. The differences between the model with additional constraints and without document that the Weddell Sea circulation exerts significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. Furthermore, a warming trend in the period 1993–2001 was found in the Weddell Sea and adjacent basins in agreement with float measurements in the upper Southern Ocean. Teleconnections to the North Atlantic are suggested but need further studies to demonstrate their statistical significance.     

Oceanography of coastal and riverine waters around Timika, West Central Irian Jaya, Arafura Sea

Oceanographic studies have been carried out in coastal and riverine waters of the area around Timika, West Papua in November 1999, March–April, July and November 2000. The temperature of the seawater along the coast is around 28 °C in winter (November 99), rising to 30.0 °C (November 00). In the open sea, 30 miles off the coast at 40 m water depth, the temperature is >30 °C with no stratification. Water temperature near the coast is consistently lower than in the open sea. This is thought to be due the cooling effect of the land, being densely covered by mangrove forest. In the upper parts of the Kamora, West Tipuka, East Tipuka, Ajkwa, Minajerwi, Mawati and Otakwa Rivers, at salinity zero psu, water temperature varies between 24.6 and 26.2 °C, which is as cold as the temperature in the upwelling Banda Sea to the NW. Some of these rivers are fed by glacial melt water from the high mountains to the east. At mid estuary, warm seawater is found under the cooler river water.Salinity near this coast varied between 24 and 30, and offshore salinity was 31–33 with no stratification. Inshore surface waters were turbid (11–14 ntu), and near bottom waters were generally much more turbid from river sediment supply and tidal resuspension. The Ajkwa River estuary has the highest turbidity (750 ntu) at zero salinity. Offshore waters were very clear (5.0–6.0 ntu), and there was no increase in turbidity near the bottom.     

On the variability of the flow along the Meso-American Barrier Reef system: a numerical model study of the influence of the Caribbean current and eddies

A high resolution (3–8 km grid), 3D numerical ocean model of the West Caribbean Sea (WCS) is used to investigate the variability and the forcing of flows near the Meso-American Barrier Reef System (MBRS) which runs along the coasts of Mexico, Belize, Guatemala and Honduras. Mesoscale variations in velocity and temperature along the reef were found in seasonal model simulations and in observations; these variations are associated with meandering of the Caribbean current (CC) and the propagation of Caribbean eddies. Diagnostic calculations and a simple assimilation technique are combined to infer the dynamically adjusted flow associated with particular eddies. The results demonstrate that when a cyclonic eddy (negative sea surface height anomaly (SSHA)) is found near the MBRS the CC shifts offshore, the cyclonic circulation in the Gulf of Honduras (GOH) intensifies, and a strong southward flow results along the reef. However, when an anticyclonic eddy (positive SSHA) is found near the reef, the CC moves onshore and the flow is predominantly westward across the reef. The model results help to explain how drifters are able to propagate in a direction opposite to the mean circulation when eddies cause a reversal of the coastal circulation. The effect of including the Meso-American Lagoon west of the Belize Reef in the model topography was also investigated, to show the importance of having accurate coastal topography in determining the variations of transports across the MBRS. The variations found in transports across the MBRS (on seasonal and mesoscale time scales) may have important consequences for biological activities along the reef such as spawning aggregations; better understanding the nature of these variations will help ongoing efforts in coral reef conservation and maintaining the health of the ecosystem in the region.     

Modelling tide–surge interaction effects using finite volume and finite element models of the Irish Sea

An unstructured mesh model of the west coast of Britain, covering the same domain and using topography and open boundary forcing that are identical to a previous validated uniform grid finite difference model of the region, is used to compare the performance of a finite volume (FV) and a finite element (FE) model of the area in determining tide–surge interaction in the region. Initial calculations show that although qualitatively both models give comparable tidal solutions in the region, comparison with observations shows that the FV model tends to under-estimate tidal amplitudes and hence background tidal friction in the eastern Irish Sea. Storm surge elevations in the eastern Irish Sea due to westerly, northerly and southerly uniform wind stresses computed with the FV model tend to be slightly higher than those computed with the FE model, due to differences in background tidal friction. However, both models showed comparable non-linear tide–surge interaction effects for all wind directions, suggesting that they can reproduce the extensive tide–surge interaction processes that occur in the eastern Irish Sea. Following on from this model comparison study, the physical processes contributing to surge generation and tide–surge interaction in the region are examined. Calculations are performed with uniform wind stresses from a range of directions, and the balance of various terms in the hydrodynamic equations is examined. A detailed comparison of the spatial variability of time series of non-linear bottom friction and non-linear momentum advection terms at six adjacent nodes at two locations in water depths of 20 and 6 m showed some spatial variability from one node to another. This suggests that even in the near coastal region, where water depths are of the order of 6 m and the mesh is fine (of order 0.5 km), there is significant spatial variability in the non-linear terms. In addition, distributions of maximum bed stress due to tides and wind forcing in nearshore regions show appreciable spatial variability. This suggests that intensive measurement campaigns and very high-resolution mesh models are required to validate and reproduce the non-linear processes that occur in these regions and to predict extreme bed stresses that can give rise to sediment movement. High-resolution meshes will also be required in pollution transport problems.     

Effect of water depth and the bottom boundary layer upon internal wave generation over abrupt topography

The role of water depth and bottom boundary layer turbulence upon lee-wave generation in sill regions is examined. Their effect upon vertical mixing is also considered. Calculations are performed using a non-hydrostatic model in cross-section form with a specified tidal forcing. Initial calculations in deeper water and a sill height such that the sill top is well removed from the surrounding bed region showed that downstream lee-wave generation and associated mixing increased as bottom friction coefficient k increased. This was associated with an increase in current shear across the sill. However, for a given k, increasing vertical eddy viscosity A _v reduced vertical shear in the across sill velocity, leading to a reduction in lee-wave amplitude and associated mixing. Subsequent calculations using shallower water showed that for a given k and A _v, lee-wave generation was reduced due to the shallower water depth and changes in the bottom boundary layer. However, in this case (unlike in the deepwater case), there is an appreciable bottom current. This gives rise to bottom mixing which in shallow water extends to mid-depth and enhances the mid-water mixing that is found on the lee side of the sill. Final calculations with deeper water but small sill height showed that lee waves could propagate over the sill, thereby reducing their contribution to mixing. In this case, bottom mixing was the major source of mixing which was mainly confined to the near bed region, with little mid-water mixing.