Geostrophic Spirals Generated by the Horizontal Diffusion of Vortex Stretching in the Yellow Sea

Horizontal velocity spirals with a clockwise rotation(downward looking) rate of 1.7?m^-1, on average, were observed in the western and northern Yellow Sea from December 2006 to February 2007. With the observed thermal wind relation,the beta-spiral theory was used to explain the dynamics of spirals. It was found that the horizontal diffusion of geostrophic vortex stretching is likely to be a major mechanism for generating geostrophic spirals. Vertical advection associated with surface/bottom Ekman pumping and topography-induced upwelling is too weak to support these spirals. Strong wind stirring and large heat loss in wintertime lead to weak stratification and diminish the effects of vertical advection. The cooling effect and vertical diffusion are offset by an overwhelming contribution of horizontal diffusion in connection with vortex stretching. The Richardson number-dependent vertical eddy diffusivity reaches a magnitude of 10^-4 m^2 s^-1 on average. An eddy diffusivity of 2870 m^2 s^-1 is required for dynamic balance by estimating the residual term. This obtained value of 10-4 m^2 s^-1 is in good agreement with the estimation in terms of observed eddy activities. The suppressed and unsuppressed diffusivities in the observation region are 2752 and 2881 m^2 s^-1, respectively, which supports a closed budget for velocity rotation.     

The turbulent bottom boundary layer and its influence on local dynamics over the continental shelf

We report on measurements of the structure of the bottom boundary layer on the continental shelf off Vancouver Island. A time series of vertical profiles obtained with the microstructure profiler FLY II revealed large temporal variations in the dissipation rate and in the density structure. The near-bottom current structure was simultaneously measured at fixed heights with conventional current meters.The data reveal the association between the predominantly diurnal tide and the variations in the structure of the bottom boundary layer. A clear distinction appears between the turbulent bottom boundary layer (8–40 m deep in a total water depth of 138 m) and the well-mixed layer (20–40 m deep). The two layers vary independently, with horizontal advection dominating the fluctuations in the thickness of the well-mixed layer while local dissipation is more closely related to the thickness of the turbulent layer. Variations in the density structure of the bottom layer are related to the strength and direction of the vertical shear and to the regional distribution of isopycnals. Current veering is commonly concentrated above the well-mixed layer. The evolution of the characteristics of the bottom layers is followed through a tidal cycle and related to local dissipation and other variables.Microstructure measurements from six locations over the southern portion of the Vancouver Island continental shelf are used to estimate the influence of turbulent energy dissipation on regional-scale flows. That fraction of the dissipation taking place in the bottom boundary layer is attributed to barotropic tidal flows, while that occurring above it is associated with nearly geostrophic baroclinic flows. Results give a lower limit of ∼ 1070 km for the length scale of shelf wave decay, in good agreement with current models; the estimates of tidal friction based on our dissipation measurements are however much lower than required by astronomical observations. An estimate of 230 h is obtained for the spin-down time of the local Tully eddy, rather longer than the decay time of 68 h obtained from observations.     

The effect of diapycnal mixing on the ventilation and CFC-11 uptake in the Southern Ocean

1. Introduction Ocean General Circulation Models (OGCMs) arekey tools in the assessment of the future ocean up-take of atmospheric greenhouse gases and heat. Fur-thermore, whereas nature experiences one realisationof the climate state, climate models can be used as alaboratory to produce a multitude of climate realisa-tions, and by that contribute to the understanding ofthe variability and stability properties of the system.It is, in this respect, crucial to evaluate the climatemodels ag…     

Winter‐time near‐surface currents in the strait of Juan de Fuca

Abstract

During November 1976 to February 1977 near‐surface wind, current and temperature measurements were made at three sites along the Strait of Juan de Fuca. Strong tidal currents and major intrusions of warmer, fresher offshore coastal water were superimposed upon the estuarine circulation of near‐surface seaward flow. The r.m.s. amplitudes of the diurnal and semidiurnal tidal currents were ～30 cms^‐1 and 30–47 cm s^‐1, respectively. The vector‐mean flow at 4 m‐depth was seaward and decreased in speed from 28 cm s^‐1 at 74 km from the entrance to 9 cm s^‐1 at 11 km from the entrance. On five occasions intrusions of 1–3 C warmer northeast Pacific coastal water occurred for durations of 1–10 days. The 25 cm s^‐1 up‐strait speed of the intrusive lens agreed to within 20% of the gravity current speed computed from Benjamin's (1968) hydraulic model. The near‐surface currents associated with the intrusions and the southerly coastal winds were significantly correlated, indicating that the intrusions were initiated when shoreward Ekman currents advected Pacific coastal water into the Strait. The reversals were not significantly coherent with the along‐strait sea surface slope measured along the north side of the Strait nor were they strongly related to local wind forcing.     

Turbulence and mixing: Sources of nutrients on the Vancouver island continental shelf

Abstract

Estimates of the rate of dissipation of turbulent energy,?_v, were made on the continental shelf off the southwest coast of Vancouver Island throughout the diurnal tidal cycle at six stations. If turbulent mixing takes place away from boundaries, and a constant fraction of the energy supplied to turbulence is converted into potential energy, this dissipation rate per unit volume, ?_v, is shown to be proportional to the rate of turbulent mass and nutrient flux across the isopycnals.

We compute these fluxes to determine the source of nutrients in the upper mixed layer on the shelf. It is found that tidal mixing on the shelf throughout the water column contributes less than 10% of the flux of nutrients supplied by the estuarine outflow out of Juan de Fuca Strait. Strong winds during upwelling events supply nutrients at rates greater than those due to tidal mixing, but at rates that are likely smaller than the Juan de Fuca source. Therefore, the nutrient‐rich waters observed in the euphotic zone in spring and summer on this shelf derive mainly from Juan de Fuca Strait. Nutrients in bottom water derive from upwelling along shelf‐break canyons.     

Dynamics of the Leeuwin Current: Part 2. Impacts of mixing,friction, and advection on a buoyancy-driven eastern boundary current over a shelf

The boundary currents over the Western Australian continental shelf and slope consist of the poleward flowing Leeuwin Current (LC) and the equatorward flowing Leeuwin Undercurrent (LUC). Key properties of the LC are its poleward strengthening, deepening to the south, and shelfbreak intensification. The alongshore flow reverses direction below about 300 m, forming the LUC at greater depths. To investigate the processes that cause these features, we obtain solutions to an idealized, regional ocean model of the South Indian Ocean. Solutions are forced by relaxing surface density to a prescribed, meridionally varying density profile ρ^*(y) with a timescale of δt. In addition, vertical diffusion is intensified near the ocean surface. This diffusion establishes the minimum thickness over which density is well-mixed. We define this thickness as the “upper layer”. Solutions are obtained with and without a continental shelf and slope off Western Australia and for a range of values of δt and mixing parameters. Within this upper layer, there is a meridional density gradient that balances a near-surface, eastward geostrophic flow. The eastward current downwells near the eastern boundary, leading to westward flow at depth. The upper layer's meridional structure and zonal currents crucially depend on coastal processes, including the presence of topography near the eastern boundary. Kelvin waves inhibit the upper layer from deepening at the coast. Rossby waves propagate the coastal density structure offshore, hence modifying the interior currents. A comparison of the solutions with or without a continental shelf and slope demonstrate that topographic trapping of Rossby waves is a necessary process for maintaining realistic eastern boundary current speeds. Significant poleward speeds occur only onshore of where the upper layer intersects the slope, that is, at a grounding line. Its poleward transport increases when surface-enhanced vertical mixing is applied over a greater depth. When the timescale δt is sufficiently short, the poleward current is nearly barotropic. The current's spatial structure over the shelf is controlled by horizontal mixing, having the structure of a Munk layer. Increasing vertical diffusion deepens the upper layer thickness and strengthens the alongshore current speed. Bottom drag leads to an offshore flow along the bottom, reducing the net onshore transport and weakening the current's poleward acceleration. When δt is long, poleward advection of buoyancy forms a density front near the shelf break, intensifying poleward speeds near the surface. With bottom drag, a bottom Ekman flow advects density offshore, shifting the jet core offshore of the shelf break. The resulting cross-shelf density gradient reverses the meridional current's direction at depth, leading to an equatorward undercurrent.     

Water exchange in fjords induced by tidally generated internal lee waves

Tidal processes are examined that control the water exchange between two basins of the Trondheimsfjord through a narrow channel with sills. For this purpose, a non-hydrostatic numerical model based on the laterally averaged Reynolds equations in the Boussinesq approximation was developed. The model takes into account the real vertical fluid stratification, variable bottom topography and variable cross-section of the fjord. Numerical experiments were performed to investigate tidally generated internal waves and their influence on the water exchange.The model produces both baroclinic tides and tidally generated lee waves. It was found that, for the Skarnsund strait which connects the Middle Fjord and the Beitstadfjord, the internal tides generated over the Skarnsund sills are very weak. Their amplitudes do not exceed 1 m.The intense short internal waves, which are identified as unsteady lee waves, comprise the basic input of the total internal wave field. These waves are generated by tidal currents at sill breaks, are trapped by topography in the generation area and grow by continuing feedback into large-amplitude waves. As the tidal flow slackens, they move upstream as freely propagating waves.As essentially nonlinear responses, the lee waves cause a nonlinear water transport. The detailed analysis of the residual currents produced by unsteady lee waves (which are propagating in both directions from the Scarnsund sills) has shown, in particular, that the residual currents can reach values as high as 0.27 m s⁻¹.It was also found that such currents exert a considerable effect on the water exchange through the Skarnsund strait between the adjacent basins. This mechanism can play an important role in water renewal and formation of the Beitasdfjord waters.     

A tidal model of the northeast pacific

Abstract

The development of a tidal model for the west coast of Canada is described. The model is intermediate in resolution between coarse‐gridded global models and fine‐gridded local models; it provides a good representation of the main shelf regions and also includes a substantial area of the neighbouring ocean. The physical processes relevant to tides in both deep and shallow water are included. Calculations have been carried out for the M₂ and K₁ constituents and the model results were compared with extensive tide‐gauge observations and empirically based charts. For M₂, the agreement between model results and observations is generally excellent, but for K₁, which contains more small‐scale variability, the model results are not quite so good. The variability in K₁ is associated with tidally generated continental shelf waves. Examination of the computed currents and energy fluxes suggests that shelf‐wave components are present in the model solution but, for the Vancouver Island shelf, their propagation is not reproduced accurately. This may be due to deficiencies in the model and/or to the influences of stratification and mean currents, which are neglected here. The model predicts that shelf‐wave components should also occur in diurnal tides on the Alaskan shelf.

The significance of the tide‐generating potential and advection are also examined and further work proposed.     

Aircraft-Based Measurements Of Turbulence Structures In The Katabatic Flow Over Greenland

Turbulence structures in the katabatic flow in the stable boundary layer (SBL) over the ice sheet are studied for two case studies with high wind speeds during the aircraft-based experiment KABEG (Katabatic wind and boundary layer front experiment around Greenland) in the area of southern Greenland. The aircraft data allow the direct determination of turbulence structures in the katabatic flow. For the first time, this allows the study of the turbulence structure in the katabatic wind system over the whole boundary layer and over a horizontal scale of 80 km.The katabatic flow is associated with a low-level jet (LLJ), with maximum wind speeds up to 25 m s^-1. Turbulent kinetic energy (TKE) and the magnitude of the turbulent fluxes show a strong decrease below the LLJ. Sensible heat fluxes at the lowest level have values down to -25 W m^-2. Latent heat fluxes are small in general, but evaporation values of up to +13 W m^-2 are also measured. Turbulence spectra show a well-defined inertial subrange and a clear spectral gap around 250-m wavelength. While turbulence intensity decreases monotonously with height above the LLJ for the upper part of the slope, high spectral intensities are also present at upper levels close to the ice edge. Normalized fluxes and variances generally follow power-law profiles in the SBL.Terms of the TKE budget are computed from the aircraft data. The TKE destruction by the negative buoyancy is found to be very small, and the dissipation rate exceeds the dynamical production.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅰ: Retrieval of Microphysical Elements and Vertical Velocity

Assuming that cloud reaches static state in the warm microphysical processes, water vapor mixing ratio(qv), cloud water mixing ratio (qc), and vertical velocity (w) can be calculated from rain water mixing ratio (qr)- Through relation of Z-qr, qr can be retrieved by radar reflectivity factor (Z). Retrieval results indicate that the distributions of mixing ratios of vapor, cloud, rain, and vertical velocity are consistent with radar images, and the three-dimensional spatial structure of the convective cloud is presented. Treating q,v saturated at the echo area, the retrieved qr is about 0.1 g kg-1, qc is always less than 0.3 g kg-1, w is usually below 0.5 m s-1, and rain droplet terminal velocity (vr) is around 5.0 m s-1 in the place where radar reflectivity factor is about 25 dBz; in the place where echo is 45 dBz, the retrieved qr and qc are always about 3.0 g kg-1, w is greater than 5.0 m s-1, and vr is around 7.0 m s-1. In the vertical, the maximum updraft velocity is greater than 3.0 m s-1 at the height of around 5.0 kin, the maximum cloud water content is about 3.0 g kg-1 above 5 km and the maximum rain water content is about 3.0 g kg-1 below 6 kin. Due to the assumption that the cloud is in static state, there will be some errors in the retrieved variables within the clouds which axe rapidly growing or dying-out, and in such cases, more sophisticated radar data control technique will help to improve the retrieval results.     

On the mean and tidal currents in Hudson strait

Abstract

The vertical structures of the mean and tidal flows in Hudson Strait are described from moored current‐meter data collected during an 8‐week period in August to October of 1982. The residual flow in the strongly stratified waters off Quebec is directed along the Strait to the southeast, is highly baroclinic and is concentrated near shore (within an offshore length scale of approximately an internal Rossby radius). Maximum mean speeds of 0.3 m s^?1 were observed near‐surface (30 m). In the weakly stratified waters on the northern side of the Strait along Baffin Island the mean flow is northwestward. The maximum speeds are 0.1 m s^?1 near‐surface (30 m) and the current amplitudes decrease to 0.05 m s^?1 at 100 m. The mean southeastward transport is estimated to be 0.93 ±0.23 × 10⁶ m³ s^?1 with a northwestward transport of 0.82 ± 0.24 × 10⁶ m³ s^?1. Over most of the Strait the across‐channel residual currents are directed towards the Quebec shore with velocities ranging from 0.02 to 0.1 ms^?1. Current variability is dominated by the tides, the M₂ being the major tidal constituent. In the vicinity of the mooring the M₂ tide is primarily barotropic, progressive in nature, and has along‐channel current amplitudes varying across the Strait from 0.20 to 0.45 m s^?1. Observed differences in tidal sea‐level elevations across the Strait can be accounted for by the cross‐channel variations characteristic of Kelvin waves.     

Sensitivity study of the seasonal mean circulation in the northern South China Sea

A study of the circulation in the northern South China Sea （SCS） is carried out with the aid of a three-dimensional, high-resolution regional ocean model. One control and two sensitivity experiments are performed to qualitatively investigate the effects of surface wind forcing, Kuroshio intrusion, and bottom topographic influence on the circulation in the northern SCS. The model results show that a branch of the Kuroshio in the upper layer can intrude into the SCS and have direct influence on the circulation over the continental shelf break in the northern SCS. There are strong southward pressure gradients along a zonal belt largely seaward of the continental slope. The pressure gradients are opposite in the southern and northern parts of the Luzon Strait, indicating inflow and outflow through the strait, respectively. The sensitivity experiments suggest that the Kuroshio intrusion is responsible for generating the imposed pressure head along the shelf break and has no obvious seasonal variations. The lateral forcing through the Luzon Strait and Taiwan Strait can induce the southwestward slope current and the northeastward SCS Warm Current in the northern SCS. Without the lateral forcing, there is the continental slope. The wind forcing mainly causes the The wind-induced water pile-up results in the southward no high-pressure-gradient zonal belt seaward of seasonal variation of the circulation in the SCS. high pressure gradient along the northwestern boundary of the basin. Without the blocking of the plateau around Dongsha Islands, the intruded Kuroshio tends to extend northwest and the SCS branch of the Kuroshio becomes wider and stronger. The analyses presented here are qualitative in nature but should lead to a better understanding of the oceanic responses in the northern SCS to these external influence factors.     

Exploring the vertical extent of breaking internal wave turbulence above deep-sea topography

A mooring equipped with 200 high-resolution temperature sensors between 6 and 404 m above the bottom was moored in 1890 m water depth above a steep, about 10° slope of Mount Josephine, NE-Atlantic. The sensors have a precision of less than 0.5 mK. They are synchronized via induction every 4 h so that the 400 m range is measured to within 0.02 s, every 1 s. Thin cables and elliptical buoyancy assured vertical mooring motions to be smaller than 0.1 m under maximum 0.2 m s⁻¹ current speeds. The local bottom slope is supercritical for semidiurnal internal tides by a factor of two. Exploring a one-month record in detail, the observations show: 1/semidiurnal tidal dominance in variations of dissipation rate ε, eddy diffusivity K_z and temperature, but no significant correlation between the records of ε and total kinetic energy, 2/a variation with time over four orders of magnitude of 100-m vertically averaged ε, 3/a local minimum in density stratification between 50 and 100 m above the bottom, 4/a gradual decrease in daily or longer averaged ε and K_z by one order of magnitude over a vertical distance of 250 m, upwards from 150 m above the bottom, 5/monthly mean values of <[ε]> = 2 ± 0.5 × 10⁻⁷ m² s⁻³, <[K_z]> = 8 ± 3 × 10⁻³ m² s⁻¹ averaged over the lower 150 m above the bottom.     

Tidal currents and mixing at the INSTANT mooring locations

Tides affect transport and mixing in the Indonesian Seas, impacting the throughflow and the return flow of the global thermohaline circulation. In a previous study, barotropic and baroclinic tides were simulated for the Indonesian Seas at 5 km resolution in order to characterize the tides of the region and to identify and quantify locations of tidal mixing. Baroclinic tidal velocities exceeded barotropic velocities except in shallow regions and their variability was on smaller scales. Model results agreed reasonably with observations and are consistent with the resolution. However, only four mooring locations were available for comparison. The new International Nusantara Stratification (INSTANT) data set enables a more comprehensive comparison. With the exception of Lombok Strait, the model replicated the observed INSTANT velocity spectra, falling within the 90% confidence limits of the observed spectra, both in regions of high and low baroclinic tidal activity for the band of frequencies from 0.02 cph to 0.33 cph (periods of 50–3 h, respectively), which includes the major semidiurnal and diurnal tides and several of their first harmonics. The model overestimated the semidiurnal baroclinic tides in the narrow Lombok Strait, which is not well resolved in the model. Comparisons of vertical profiles of the major axes of the tidal ellipses at the mooring sites generally reproduced the vertical pattern, although there were exceptions, such as Lombok and Ombai Straits. Rms differences between the model estimates and hourly observations for the major axes of the tidal ellipses were typically 1–8 cm s⁻¹ in regions of high tidal activity, 1–5 cm s⁻¹ in regions of low tidal activity, and 1–20 cm s⁻¹ for the semidiurnal tides in Lombok and Ombai Straits. Rms errors of 1–6 cm s⁻¹ are typical in regions of moderate baroclinic tidal activity at this model resolution (5 km). Many of the larger rms differences result from vertical discrepancies in the depths of the internal tidal beams. The local nature of the internal tides generation and beam propagation results in large differences from small vertical shifts in the beams or generation due to topographic differences between the model topography and the actual topography. In addition, the moorings experienced severe blowdown. The blowdown adds uncertainty to the depths of the instruments and introduces errors in the observational tidal analysis in magnitude of the tidal constituents, both of which contribute to rms differences. Tidal mixing was found to occur in intense local regions with strong internal tidal shear. The local regions of mixing were typically along the bottom in steep slopes and over sills. In conclusion, the tidal model was found to reproduce the kinetic energy distribution and transfer of energy from tides to other frequencies in the Indonesian Seas and to roughly replicate the observed structure and magnitude of the tidal currents. Improvements in the tidal simulations in reproducing observations are expected with increased resolution.     

2021年黑龙江两次超级单体龙卷过程多尺度特征

2021年6月1日和6月9日黑龙江省哈尔滨尚志市及阿城区和齐齐哈尔梅里斯区分别发生双龙卷事件。利用常规气象观测、多普勒天气雷达等资料对比分析二者的多尺度特征。结果表明:两次龙卷均发生在东北冷涡的东南象限，高空急流出口区左侧，中低层偏南气流有利于暖湿气流输送和垂直运动发展。6月1日和6月9日分别以短时强降水和雷暴大风天气为主，6月1日水汽条件、垂直运动、0~1 km高度垂直风切变和抬升凝结高度更有利于产生强龙卷，且中尺度气旋维持时间更长。干线与地面辐合线为中尺度触发机制。雷暴冷池出流与中尺度暖锋形成的伪冷锋有利于龙卷的发展和维持。龙卷出现在地面伪冷锋与干线交界处的湿区一侧，冷池前沿，龙卷母云为超级单体。暖湿气流产生的入流缺口是钩状回波发展的前兆，中等到高强度的中尺度气旋在3 km高度产生并发展，5~10 min后触地，当钩状回波与中尺度气旋同时出现时龙卷产生。     

A numerical world ocean general circulation model: Part II. A baroclinic experiment

A new six-layer world ocean general circulation model based on the primitive system of equations is described in detail and its performance in the case of a homogeneous ocean is described. These test integrations show that the model is capable of reproducing the observed mean barotropic or vertically-integrated transport, as well as the seasonal variability of the major ocean gyres. The surface currents, however, are dominated by the Ekman transport, and such non-linear features as the western boundary currents and the equatorial countercurrents are poorly represented. The abyssal boundary countercurrents are also absent due to the lack of thermohaline forcing. The most conspicuous effect of the bottom topography on a homogeneous ocean is seen in the Southern ocean where the calculated Antarctic circumpolar transport through the Drake passage ( ≈ 10 Sv, with bathymetry included) greatly underestimates the observed transport (≈ 100 Sv).     

A numerical world ocean general circulation model Part I. Basic design and barotropic experiment

A new six-layer world ocean general circulation model based on the primitive system of equations is described in detail and its performance in the case of a homogeneous ocean is described. These test integrations show that the model is capable of reproducing the observed mean barotropic or vertically-integrated transport, as well as the seasonal variability of the major ocean gyres. The surface currents, however, are dominated by the Ekman transport, and such non-linear features as the western boundary currents and the equatorial countercurrents are poorly represented. The abyssal boundary countercurrents are also absent due to the lack of thermohaline forcing. The most conspicuous effect of the bottom topography on a homogeneous ocean is seen in the Southern ocean where the calculated Antarctic circumpolar transport through the Drake passage ( ≈ 10 Sv, with bathymetry included) greatly underestimates the observed transport (≈ 100 Sv).     

极端雷暴大风的环境参量特征

为了研究极端雷暴大风天气环境要素特点,选取2002—2017年中国各地区极端雷暴大风个例95个和不伴随强对流的普通雷暴个例95个,通过两者间关键环境参数的对比,揭示极端雷暴大风事件的关键环境参数特征。结果表明:极端雷暴大风天气发生在对流层中层相对干的环境下,表现为400~700 hPa极端雷暴大风对应的单层最大温度露点差和平均温度露点差平均值分别为25.7℃和13.6℃,而普通雷暴的相应值分别为16.2℃和6.5℃。统计结果表明:尽管产生极端雷暴大风的对流风暴和普通雷达对应的地面露点差异并不大,但前者相应的大气可降水量（平均值为37 mm）明显低于后者（平均值为51 mm）,差异突出表现在两者湿层厚度的不同上;相对于普通雷暴事件,极端雷暴大风事件对应的对流有效位能值（平均值为1820 J·kg-1）明显高于普通雷暴事件的对应值（平均值为470 J·kg-1）;此外,极端雷暴大风事件对应的对流层中下层垂直温度递减率、下沉有效位能、夹卷层平均风速和0~6 km,0~3 km垂直风切变均明显大于普通雷暴事件对应的相应值。     

Experiments with density currents on a sloping bottom in a rotating fluid

The behaviour of a density current on a sloping bottom in a rotating system is investigated by laboratory experiments. The main result is that the dense bottom outflow induces cyclonic vortices in the upper fluid layer, which are formed periodically and move to the west parallel to coast. Two regimes of vortex formation have been identified. For strong density currents and weak rotation, vortices are formed by stretching of the upper layer near the source as found also in the experiments by Lane-Serff and Baines (1998) [Lane-Serff, G.F., Baines, P.G., 1998. Eddy formation by dense flows on slopes in a rotating fluid. J. Fluid Mech. 363, 229–253]. For weak density currents and strong rotation vortices are due to instability of the bottom plume itself as found in the numerical simulations of Jiang and Garwood (1996) [Jiang, L., Garwood, W. Jr., 1996. Three-dimensional simulations of overflows on continental slopes. J. Phys. Oceanogr. 26, 1224–1233].     

Tidal circulation and buoyancy effects in the St. Lawrence Estuary

Abstract

The action of tides on density‐driven circulation, internal gravity waves, and mixing was investigated in the St. Lawrence Estuary between Rimouski and Québec City. Time‐varying fields of water level, currents and density were computed under typical summer conditions using a three‐dimensional hydrostatic coastal ocean model that incorporates a second order turbulence closure submodel. These results are compared with current meter records and other observations. The model and the observations reveal buoyancy effects produced by tidal forcing. The semi‐diurnal tide raises the isopycnals over the sills at the head of the Laurentian Trough and English Bank, producing internal tides radiating seaward. Relatively dense intermediate waters rise from below 75‐m depth to the near surface over the sills, setting up gravity currents on the inner slopes. Internal hydraulic controls develop over the outer sills; during flood, surface flow separation occurs at the entrances of the Saguenay Fjord and the upper estuary west of Ilet Rouge Bank. Early during ebb flow (restratification), the surface layer deepens to encompass the tops of the sills. As the ebb current intensifies, the model predicts the formation of seaward internal jumps over the outer sills, which were confirmed from acoustic reflection observations. As the internal Froude number increases further, flow separation migrates up to sill height. As a result of these transitions, internal bores emanate from the head region one to two hours before low water. We find that the mixing of oceanic and surface waters near the sills is driven by the vertical shear produced during ebb in the channel south of Ilet Rouge, the shear produced in the bottom gravity flood currents, and, to a lesser extent, the processes over the sills.