Coastal currents and eddies and their interaction with topography

The nonlinear interaction of vorticity driven coastal currents and eddies with topography is studied. The topography is either a semi-infinite escarpment perpendicular to the coast (such that topographic waves propagate toward the coast) or a semi-circular canyon or seamount attached to the coast. Assuming a piecewise constant potential vorticity distribution, the quasigeostrophic equations are solved using contour dynamics. Offshore propagating dipole eddies occur, whenever a coastal current or eddy interacts with escarpment and canyon topographies. The size and frequency at which dipoles form are found to depend on the vorticity of the current and amplitude of the topography. However, for a seamount, little eddy shedding is observed and the coastal current or eddy skirts around topography.     

A laboratory simulation of mesoscale flow interaction with the Alps

A series of laboratory experiments, aimed at the simulation of some aspects of Alpine lee cyclogenesis has been carried out in the rotating tank of the Coriolis Laboratory of LEGI-IMG in Grenoble. Dynamic and thermodynamic processes, typical of baroclinic development triggered by the orography, were simulated. The background flow simulating the basic state of the atmosphere consisted of a stream of intermediate density fluid introduced at the interface between two fluid layers. The structure of the intermediate current was established by mixing fluid obtained from the upper layer of fresh water with fluid removed from the heavier salty layer below.The dynamical similarity parameters are the Rossby (Ro), Burger (Bu) and Ekman (Ek) numbers, although this last, owing to its small values, need not be matched between model and prototype, since viscous effects are not important for small time scales. The flow in both the prototype and laboratory simulation is characterized by hydrostatics; this requires (Ro²δ²/Bu)1 (where δ=H/L is the aspect ratio of the obstacle) which is clearly satisfied, in the atmosphere and oceans, and for the laboratory experiment.A range of experiments for various Rossby and Burger numbers were conducted which delimited the region of parameter space for which background flows akin to that found to the northwest of the Alps prior to baroclinic cyclogenesis events, were observed.One such experiment was carried out by placing a model of the Alps at the appropriate place in the flow field. The subsequent motion in the laboratory was observed and dye tracer motions were used to obtain the approximate particle trajectories. The density field was also analyzed to provide the geopotential field of the simulated atmosphere. Using standard transformations from the similarity analysis, the laboratory observations were related to the prototype atmosphere. The flow and the geopotential fields gave results compatible with the particular atmospheric event presented.     

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.     

Instabilities of the tidally induced bottom boundary layer in the rotating frame and their mixing effect

To investigate the stability of the bottom boundary layer induced by tidal flow (oscillating flow) in a rotating frame, numerical experiments have been carried out with a two-dimensional non-hydrostatic model. Under homogeneous conditions three types of instability are found depending on the temporal Rossby number Ro_t, the ratio of the inertial and tidal periods. When Ro_t < 0.9 (subinertial range), the Ekman type I instability occurs because the effect of rotation is dominant though the flow becomes more stable than the steady Ekman flow with increasing Ro_t. When Ro_t > 1.1 (superinertial range), the Stokes layer instability is excited as in the absence of rotation. When 0.9 < Ro_t < 1.1 (near-inertial range), the Ekman type I or type II instability appears as in the steady Ekman layer. Being much thickened (100 m), the boundary layer becomes unstable even if tidal flow is weak (5 cm/s). The large vertical scale enhances the contribution of the Coriolis effect to destabilization, so that the type II instability tends to appear when Ro_t > 1.0. However, when Ro_t < 1.0, the type I instability rather than the type II instability appears because the downward phase change of tidal flow acts to suppress the latter. To evaluate the mixing effect of these instabilities, some experiments have been executed under a weak stratification peculiar to polar oceans (the buoyancy frequency N²  10⁻⁶ s⁻²). Strong mixing occurs in the subinertial and near-inertial ranges such that tracer is well mixed in the boundary layer and an apparent diffusivity there is evaluated at 150–300 cm²/s. This suggests that effective mixing due to these instabilities may play an important role in determining the properties of dense shelf water in the polar regions.     

Polytropic atmospheres and the scaling of potential vorticity

Summary Potential vorticity P is in widespread use as a diagnostic in atmospheric dynamics. Here, we employ the concepts of polytropic change of state and polytropic temperature for investigating the problem of scaling potential vorticity. For many applications, it has proven useful to scale P with a dimensionless scale factor to remove its exponential increase with altitude. This scale factor may be derived from the assumption of a polytropic reference atmosphere which includes the classical scaling using isothermal reference atmospheres as a special case.     

The vorticity balance on the Southern British Columbia continental shelf

Abstract

An array of moored current meters was deployed over a canyon system on the British Columbia continental shelf in the spring of 1985 and maintained for 110 days. The array was specifically designed to monitor the various terms in the vorticity budget.

The observations indicate that the substantial variations in vorticity with time are being driven by vertical velocities originating from the head of a small canyon. These vertical velocities are in turn being directly driven by deep northward flows along the axis of the canyon and these latter flows are closely related to the southward shelf‐edge‐flow.     

Steering of upper ocean currents and fronts by the topographically constrained abyssal circulation

A two-layer theory is used to investigate (1) the steering of upper ocean current pathways by topographically constrained abyssal currents that do not impinge on the bottom topography and (2) its application to upper ocean – topographic coupling via flow instabilities where topographically constrained eddy-driven deep mean flows in turn steer the mean pathways of upper ocean currents and associated fronts. In earlier studies the two-layer theory was applied to ocean models with low vertical resolution (2–6 layers). Here we investigate its relevance to complex ocean general circulation models (OGCMs) with high vertical resolution that are designed to simulate a wide range of ocean processes. The theory can be easily applied to models ranging from idealized to complex OGCMs, provided it is valid for the application. It can also be used in understanding some persistent features seen in observed ocean frontal pathways (over deep water) derived from satellite imagery and other data. To facilitate its application, a more thorough explanation of the theory is presented that emphasizes its range of validity. Three regions of the world ocean are used to investigate its application to eddy-resolving ocean models with high vertical resolution, including one where an assumption of the two-layer theory is violated. Results from the OGCMs with high vertical resolution are compared to those from models with low vertical resolution and to observations. In the Kuroshio region upper ocean – topographic coupling via flow instabilities and a modest seamount complex are used to explain the observed northward mean meander east of Japan where the Kuroshio separates from the coast. The Japan/East Sea (JES) is used to demonstrate the impact of upper ocean – topographic coupling in a relatively weak flow regime. East of South Island, New Zealand, the Southland Current is an observed western boundary current that flows in a direction counter to the demands of Sverdrup flow and counter to the direction simulated in nonlinear global flat bottom and reduced gravity models. A model with high vertical resolution (and topography extending through any number of layers) and a model with low vertical resolution (and vertically compressed but otherwise realistic topography confined to the lowest layer) both simulate a Southland Current in the observed direction with dynamics depending on the configuration of the regional seafloor. However, the dynamics of these simulations are very different because the Campbell Plateau and Chatham Rise east and southeast of New Zealand are rare features of the world ocean where the topography intrudes into the stratified water column over a relatively broad area but lies deeper than the nominal 200 m depth of the continental shelf break, violating a limitation of the two-layer theory. Observations confirm the results from the high vertical resolution model. Overall, the model simulations show increasingly widespread upper ocean – topographic coupling via flow instabilities as the horizontal resolution of the ocean models is increased, but fine resolution of mesoscale variability and the associated flow instabilities are required to obtain sufficient coupling. As a result, this type of coupling is critical in distinguishing between eddy-resolving and eddy-permitting ocean models in regions where it occurs.     

Stratified flow over three-dimensional topography

In order to investigate flows over topography in an atmospheric context, we have studied experimentally the wake structure of axi-symmetric Gaussian obstacles towed through a linearly stratified fluid. Three dimensionless parameters govern the flow dynamics: F, the Froude number based on the topography height h; Re, the Reynolds number and the aspect ratio r = h/L, where L is the topography horizontal scale. Two-dimensional (2-D), saturated lee wave (SLW) and three-dimensional (3-D) regimes, as defined in Chomaz et al. (1993), are found to be functions of F and r only (Fig. 1) as soon as Re is larger than Re_c ≈ 2000. For F < 0.7 the flow goes around the obstacle and the motion in the wake is quasi-two-dimensional. This 2-D layer is topped by a region affected by lee wave motions with amplitude increasing with r and F. For 0.7 < F < 1/r, the flow is entirely dominated by a lee wave of saturated amplitude which suppresses the separation of the boundary layer from the obstacle. Above the critical value 1/r, the lee wave amplitude decreases with F and a recirculating zone appears behind the obstacle. Simultaneously, coherent large-scale vortices start to be shed periodically from the wake at a Strouhal number which decreases as 1/F until it reaches its neutral asymptotic value.     

Steering of upper ocean currents and fronts by the topographically constrained abyssal circulation

A two-layer theory is used to investigate (1) the steering of upper ocean current pathways by topographically constrained abyssal currents that do not impinge on the bottom topography and (2) its application to upper ocean – topographic coupling via flow instabilities where topographically constrained eddy-driven deep mean flows in turn steer the mean pathways of upper ocean currents and associated fronts. In earlier studies the two-layer theory was applied to ocean models with low vertical resolution (2–6 layers). Here we investigate its relevance to complex ocean general circulation models (OGCMs) with high vertical resolution that are designed to simulate a wide range of ocean processes. The theory can be easily applied to models ranging from idealized to complex OGCMs, provided it is valid for the application. It can also be used in understanding some persistent features seen in observed ocean frontal pathways (over deep water) derived from satellite imagery and other data. To facilitate its application, a more thorough explanation of the theory is presented that emphasizes its range of validity. Three regions of the world ocean are used to investigate its application to eddy-resolving ocean models with high vertical resolution, including one where an assumption of the two-layer theory is violated. Results from the OGCMs with high vertical resolution are compared to those from models with low vertical resolution and to observations. In the Kuroshio region upper ocean – topographic coupling via flow instabilities and a modest seamount complex are used to explain the observed northward mean meander east of Japan where the Kuroshio separates from the coast. The Japan/East Sea (JES) is used to demonstrate the impact of upper ocean – topographic coupling in a relatively weak flow regime. East of South Island, New Zealand, the Southland Current is an observed western boundary current that flows in a direction counter to the demands of Sverdrup flow and counter to the direction simulated in nonlinear global flat bottom and reduced gravity models. A model with high vertical resolution (and topography extending through any number of layers) and a model with low vertical resolution (and vertically compressed but otherwise realistic topography confined to the lowest layer) both simulate a Southland Current in the observed direction with dynamics depending on the configuration of the regional seafloor. However, the dynamics of these simulations are very different because the Campbell Plateau and Chatham Rise east and southeast of New Zealand are rare features of the world ocean where the topography intrudes into the stratified water column over a relatively broad area but lies deeper than the nominal 200 m depth of the continental shelf break, violating a limitation of the two-layer theory. Observations confirm the results from the high vertical resolution model. Overall, the model simulations show increasingly widespread upper ocean – topographic coupling via flow instabilities as the horizontal resolution of the ocean models is increased, but fine resolution of mesoscale variability and the associated flow instabilities are required to obtain sufficient coupling. As a result, this type of coupling is critical in distinguishing between eddy-resolving and eddy-permitting ocean models in regions where it occurs.     

Island Wakes in Shallow Water

Island wakes in shallow water are investigated using the Regional Ocean Modeling System (ROMS). In contrast to deep water where bottom stress can be neglected in island wakes, shallow water implies that inhomogeneity in the bottom stress plays an important role in the wake vorticity generation. A series of numerical experiments are conducted to investigate wake formation and evolution in shallow water. It is found that the vertical structures of shallow-water and deep-water wakes are significantly different because of the presence of a density frontal jet, which results from the interaction between stratification and bottom topography. The frontal jet reaches its maximum within the bottom boundary layer over the shelf, giving rise to vorticity. The potential vorticity (PV) balance analysis reveals that frictional and diapycnal processes play different roles in the PV anomalies. With the absence of lateral stress (i.e., a sea mountain case), the surface vorticity becomes much weaker than that in the presence of an island.     

Controlling the open ocean forcing of baroclinic coastal jets with a 4D-var method

An optimal control approach was carried out to investigate the influence of perturbations of baroclinic coastal jets which are located upstream of the observation site. Two academic numerical experiments were considered: the intrusion of the buoyancy anomaly within an established density current along a straight coast (experiments A), and the interaction of a coastal current with a shelf topography (experiments B). We explored the sensitivity of the forcing retrieval to the spatial distribution of observations (experiments A) and to the temporal sampling of the forcing (experiments B), using a 4D-var formalism. The sensitivity of numerical simulations is discussed in the framework of analytical solutions. The results provide insight into how observational network and data assimilation method would impact the realistic modelling of coastal jets and the shelf slope exchanges.     

The cross‐channel flow at the entrance of Lancaster Sound

Abstract

A westward current flows along the northern side of Lancaster Sound and an eastward current flows along the southern side. A cross‐channel flow is commonly observed to link them near the eastern entrance of Lancaster Sound; this flow is modelled assuming inviscid flow and conservation of potential vorticity. It is shown that the westward decrease in depth is sufficient to cause a cross‐channel flow that couples the inflow to the outflow. The modelled cross‐channel flow takes place at a distance inside the entrance that is less than that observed for the surface current. Obviously stratification reduces the coupling of the surface current to the bathymetry. A more realistic result is obtained with the barotropic model if the bottom slope is halved.

An inviscid mean barotropic flow out of the channel is also modelled and found to be concentrated on the southern bank in order to conserve potential vorticity. It seems that barotropic instability and friction would limit the narrowing of the flow.     

Weakly nonlinear oscillatory flow over an isolated seamount

Abstract

Numerical models can show large errors in the modelling of flow around tall topography. We present a weakly nonlinear analytical solution of barotropic oscillatory flow over a tall axisymmetric seamount and compare the results to a numerical tidal model. The comparison is formulated for a flat topped seamount with parabolic sides and for a standing wave boundary condition in order to make the problem analytically tractable and to reduce it to a series of ordinary differential equations. The comparison shows reasonable agreement within the confines of the analytical formulation confirming that the numerical model can exploit tall topography. Comparisons to a Gaussian seamount and Kelvin wave forcing illustrate the generality of the results.

The analytical formulation produces an equation relating relative errors in the comparison to discontinuities in topographical slope; the numerical model should have smaller errors over smoother topography. Residual velocities are found to be proportional to the aspect ratio of the topography.     

A STUDY ON MESOSCALE VORTICITY PROPAGATION IN A TYPHOON-LIKE VORTEX*

By using an f-plane barotropic quasi geostrophic model in the rectangular coordinates with a grid spacing of 5 km,ten experiments whose integration time is 36 hours are performed in order to study the interaction between a typhoon vortex and a mesoscale vortex whose initial center position is located at 2 r_m northwest to the typhoon center,where r_m is the radius of maximum wind of the typhoon vortex. Results show that the interaction can create a pair of smaller scale vortices or lumps,which extend from the outside region of the typhoon to near its center,resulting in the inward propagation of mesoscale vorticity.In this process,the vorticity concentration of the mesoscale vortex may appear.The coexistence of the propagation and the concentration makes the increase of vorticity in the inside region i.e.a more intensive typhoon.Meanwhile,the intensity of the lump with positive vorticity oscillates with time,with the oscillation period being several hours,the distance from the typhoon center to the lump center also has a similar oscillation period,which reduces the oscillation of typhoon intensity.In the case of stronger circular basic flow,the interaction can make the intensification of typhoon more obviously. In addition,in some parametric conditions,the interaction may break down the continuous vorticity zone,exhibiting a cluster of smaller vorticity lumps.     

Internal tide generation at a continental shelf/slope junction: A comparison between theory and a laboratory experiment

Observations of internal tide generation over continental slopes in a laboratory experiment have been carried out, with the objectives of making comparisons with linear generation theory and investigating its limitations. Both continuous and layered stratification have been considered. A measure of the amplitude of the barotropic tidal forcing (and hence of non-linearity) is given by the Froude Number F = u_sb/c_w, where u_sb is the maximum barotropic tidal velocity at the shelf break, and c_w is the long-wave speed of the lowest internal mode.For continuous stratification, good agreement was obtained for “steep” slopes (α/c > 1), where α is the slope at the continental slope and c is the slope of the internal wave rays of tidal frequency), even for quite large amplitude motions (F < 1.6), and the upper limit of its quantitative usefulness was not reached. For “flat” slopes (α/c < 1) reasonable agreement was also obtained, even up to quite amplitudes (F < 3.1), although some departure from linear theory was apparent.For two-layer flows the applicability of linear theory was much more restricted. For F 0.5 there was only qualitative agreement and for larger F (>1) significant differences were observed. The latter were due to the substantial advection and associated hydraulic jumps which occured seaward of the shelf-break during the ebb-phase of the barotropic tide. Shelf-break values of F > 1 are common in the ocean.     

Rectified flow along a vertical coastline

Laboratory experiments are conducted on a physical system in which an oscillatory, along-shore, free stream flow of a homogeneous fluid occurs in the vicinity of a long coastline with vertical slope; the model sea-floor is horizontal. Particular attention is given to the resulting rectified (mean) current which is along the coastline with the shore on the right, facing downstream. In the lateral far field region defined by (1), where y is the offshore coordinate and H is the depth of the fluid, the motion field is approximately independent of the lateral distance from the coast. The vertical structure of the cross-stream motion in this region consists of Ekman layers near the sea-floor and interior adjustment flows, both periodic in time. In the near field, defined by (1), the motion is strongly dependent on the cross-stream coordinate as well as time, and rectified currents are observed. The mechanism responsible for the rectification is a complex nonlinear coupling between laterally directed adjustment flows driven by the transport in the bottom Ekman layers, and the free stream motion field. The rectified current is found to be substantially wider than the Stewartson layer thickness but much narrower than the Rossby deformation radius. The characteristic width, δ_y, of the rectified current is shown to scale as , where Ro is the Rossby number Ro_t is the temporal Rossby number and E is the Ekman number. Experiments are presented which support this scaling.     

Surface-layer flow in complex terrain: Comparison of models and full-scale observations

Four models of surface boundary-layer flow in complex terrain are compared with observations made at Blashaval Hill, North Uist, Scotland. The field experiment is described by Mason and King (1985). Three of the models are derived from the two-dimensional theory of Jackson and Hunt (1975) and are described in Mason and King (1985), Walmsley et al. (1986) and Troen and Petersen (1989). The fourth is a mass-consistent code based on Traci et al. (1979). The model results are in good agreement with each other and are generally within the observed range of variation ( ~ ± 16%) in normalized wind speed. For most wind direcions (7 of 11), model results of normalized wind speed at the summit were within 7% of the observed mean values. For some wind directions, calculations using the Guidelines of Walmsley et al. (1989) suggested that variations in surface roughness were important. This led us to apply one of our models incorporating nonuniform surface roughness. The lack of significant improvement for cases when water lay upstream of Blashaval Hill is attributed to compensating changes at summit and reference sites and to very local effects on the wind data. Sensitivity to topography lying to the west and northwest of Blashaval was also investigated. Results suggested an influence from those distant topographic features for some wind directions. When those features were incorporated, maximum errors in normalized wind speed at the summit were reduced from 18 to 13%.     

A note on the overall vorticity balance in wind-driven ocean circulation

Some constraints resulting from a required overall vorticity balance in an ocean forced by a net wind-stress curl are considered. It is pointed out that the Sverdrup balance holds for Ekman-geostrophic flow in an area-averaged sense even in the presence of topography. It is suggested that nonlinear stretching and twisting of vorticity in an inertial boundary current may provide a source of vorticity which may balance a net wind-stress curl, without help of frictional effects at the bottom or sides. It is also suggested that potential vorticity generated by the wind-stress curl can be balanced by other local potential vorticity sources in the top layer. At this stage no proof for the need of bottom or side friction in a generally stratified ocean exists, and a counterexample is also lacking, leaving the question open for further studies.     

"7.21"北京大暴雨过程的地形作用分析和数值试验研究

利用NCEP 1°×1°的6 h再分析资料和常规气象观测资料,对2012年7月21-22日北京地区的大暴雨过程中的地形作用进行诊断分析,研究其对暴雨过程前后的热力、动力及水汽条件的影响。并利用WRFV3.2模式尝试进行一些地形敏感性数值试验,探讨地形改变后对此次暴雨过程的影响,初步得出了一些结论：高空槽配合中层切变线,在低层形成低涡,低涡的东移并加强引发了此次暴雨。暴雨发生前,中低层强的正涡度中心和负散度中心配合很好,它们在东移过程中,组织起了一个很好的垂直环流圈。此环流圈加强了迎风坡气流的上升运动;对流层中低层东南和西南两条水汽通道提供了强大的偏南气流,不断往北京地区上空输送暖湿空气,从而在低层积累了大量的不稳定能量;对流层中低层的位涡从高处往下传,且正位涡中心在东移过程中加强,对应山前迎风坡出现了局地大暴雨,由此说明位涡中心与暴雨落区具有良好的对应关系。地形敏感性数值试验验证了地形虽然对于整体降雨带的影响不大,但对于局地暴雨落区和中心强度的影响却很大。