A new barotropic model of the wind-driven circulation

Rationalized by the observational circulation pattern in the upper ocean of the North Pacific, meridional friction term is first incorporated in a barotropic theoretical model of the wind-driven circulation. The governing potential vorticity equation thence has β term and wind stress curl term (the two of the Sverdrup balance), zonal friction term and meridional friction term. The analytical solution satisfactorily captures many important features of the wind-driven circulation in the North Pacific: Kuroshio, Oyashio, Kuroshio extension, North Equatorial Current, and especially the eastern boundary currents in the North Pacific, i.e. California current and Alaska current.     

Vertical structure of residual slope circulation driven by JEBAR and tides: an idealised model

The residual circulation over the continental slope, and in particular, its vertical structure, is analysed by means of an idealised hydrodynamic model. The model is based on the depth-dependent shallow-water equations under uniform along-isobath conditions and is forced by a prescribed meridional density gradient and tidal velocities. By means of expansion in the small Rossby number solutions are analysed for conditions representative for the continental slopes off the Hebrides and in the Bay of Biscay. The steady solution at zeroth order consists of a linear density-driven flow. At order a tidally rectified flow is found and a stationary flow due to self-interaction of the zeroth-order density-driven flow. At order ² the leading-order effect of the interaction between the zeroth-order density-driven flow and the tides is found: the ‘interaction current’. The solutions up to and including order ² constitute an along-isobath steady slope current which is comparable to field data. The slope current and the accompanying cross-shelf circulation depend strongly on the shelf and flow characteristics. For the Hebridean case the density forcing predominates, but for the Biscay case the tidal effects are of the same order of magnitude as the density effects. Under those conditions the interaction current is significant which implies that linear superposition of density and tidal effects differs from the non-linear combination of both. It is also shown that the depth-average of the interaction current differs essentially from the solution obtained from a depth-averaged model.     

Dynamics,sediment transport,and morphology in a tide-dominated embayment

Western Port, Victoria, Australia is a tide-dominated embayment with an unusual and complex shape. Bottom currents and circulation and their effects on sediment processes were examined using instrumented tripods to measure currents, tides and wave activity, and to estimate sediment transport at nine locations in the bay. Overall bottom water movement patterns were determined by use of seabed drifters. The characteristics of the bay reflect a small catchment and low freshwater and sediment input. A complex system of channels is flanked by extensive intertidal areas. Tidal range varies up to more than 3 m, generating bottom currents up to 70cm s^?1. Flow directions generally conform to channel alignments but major deviations are important. Net circulation in the bay is clockwise around the large central island. The ratio of tidal range to half-tidal period (ΔH/ΔT) when compared with measured currents gave a method of prediction of the annual frequency distribution of maximum bottom current velocities. Determination of threshold current velocities enabled prediction of annual frequency of bedload movement (generally 50–100 per cent of tide cycles). Bedload mass transport for all observed tide cycles was calculated, and estimates of annual mass transport capacity (between 10² and 10⁶ g cm^?1 a^?1( were obtained from a relationship between predicted tidal conditions and mass transport. Seabed drifters delineated the major bottom water movements in the bay and adjacent Bass Strait, and also detailed circulation patterns. This linked the data from the tripod stations, in particular patterns of ebb- and flood-dominance. A coherent picture of the processes operating in Western Port is presented by integrating these studies with corroborative studies of sediment distribution and morphology, hydrochemistry, and mathematical modelling. Some morphological characteristics related to tide-dominance are discussed.     

Global surface currents: a high-resolution product for investigating ocean dynamics

A global 1/4° resolution product of surface currents has been developed by the Centre de Topographie des Océans et de l’Hydrosphère. The surface current is calculated from a combination of Ekman currents derived from wind estimates from QuikSCAT satellite, geostrophic current anomalies derived from altimetry, and a mean geostrophic current derived from climatology. In the equatorial band, the currents are adjusted following the methodology proposed by Lagerloef et al. (J Geophys Res, 104(C10):22313–22326, 1999). These satellite-derived currents have been compared to different types of in situ current observations. A global validation is performed using Lagrangian surface drifting buoys and acoustic Doppler current profiler current observations along ship tracks. The comparison shows a very good agreement in the subtropical and mid-latitude bands. The correlation between the satellite-derived currents and the drifter currents in zonal mean bands is around 0.7 for most of the world oceans, both for the zonal and the meridional components. This correlation rises up to 0.8 in the regions of strong boundary currents. In the equatorial band, the correlation with the surface drifting buoys is reduced. A direct comparison with the TOGA/TAO moored current meter data at the equator shows that the low frequency currents are captured by the satellite current product, but there is a substantial high-frequency signal (<20 days), which is not reproduced. This is especially the case for the meridional component and is mainly related to the tropical instability waves. We also show that using daily QuikSCAT wind forcing improves the satellite current product, particularly in the high-latitude westerly wind belt and in the tropical Indian Ocean.     

An externally forced barotropic circulation model for the New York Bight

The sea level and the barotropic, frictional circulation response for the New York Bight are used to demonstrate the effects of external sea-level forcing, bathymetry, and variable friction. The governing equation is the steady, integrated vorticity equation and is computed by finite differencing over a curvilinear grid conforming to the 10- and 100-m isobaths and extending for 250 km alongshore. The boundary conditions are based on the hypothesis that the dynamics of the shelf are driven by the external sea-level gradient and the coastal no-flux condition; and consequently the conditions at the lateral boundaries are dependent thereon. Therefore, the external sea-level slope must be independently specified, and the lateral boundary conditions must be dependently generated. The diabathic component of the external sea slope forces the calm wind circulation by its effect on the transport through the upstream boundary; and the parabathic component has also an important modifying effect by forcing a shelf convergent transport. The parabathic sea slope at the coast is independent of its offshore value, being instead a direct product of the coastal boundary condition.The bottom friction is expressed as related to the sea level through a bottom length parameter and a veer angle, both of which are taken to increase shoreward. An additional bottom stress component, related to the surface stress, is determined for bottom depths less than the Ekman depth. Such bottom stress variability produces significant alterations in the nearshore flow field, over the constant bottom stress formulation, by reducing it and causing it to veer downgradient and downwind in the nearshore.The model is forced by different wind directions and the results are discussed. The circulations generally conform to the observed mean flow patterns, but with several smaller-scale features. The strong bathymetric feature of the Hudson Shelf Valley causes a polarized up- and downvalley flow for winds with an eastward or westward component, respectively. Under mean westerly winds, there is a divergence in the shelf valley flow at about the 60-m isobath. The Apex gyre existing off the western tip of Long Island becomes more extensive for winds changing from northeast to southwest. Mean flow reversals (to the northeast) occur off both Long Island and New Jersey for wind directions changing counterclockwise from northwest to southeast and from west to east, respectively. Southeastward transport over the outer New Jersey shelf tends to be enhanced by wind and external sea-level conditions; and the transport over the New Jersey midshelf, particularly in the lee of the shelf valley, tends to be weak and variable also under these mean conditions.     

Numerical models of continental and submarine erosion: application to the northern Ligurian Margin (Southern Alps,France/Italy)

We present a new numerical surface process model allowing us to take into account submarine erosion processes due either to submarine landslides or to hyperpycnal currents. A first set of models show that the frequency of hyperpycnal flows influences the development of submarine canyons at the mouth of continental rivers. Further experiments show that an increase in submarine slope leads to faster regressive canyon erosion and a more dentritic canyon network, whereas increasing the height of the unstable sediment pile located on the shelf break leads to wider and less dendritic canyons. The models are then applied to the western segment of the north Ligurian margin (northwestern Mediterranean), which displays numerous submarine canyons with various sizes and morphologies. From west to east, canyon longitudinal profiles as well as margin‐perpendicular profiles progressively change from moderately steep, concave‐up shapes to steeper linear to convex‐up shapes suggesting increasing eastward margin uplift. Moreover, the foot of the margin is affected by a marked slope increase with evidences of mass transport due to landslides. Numerical models which reproduce well the North Ligurian margin morphologic features indicate that the western part of the margin is submitted to rather low (i.e. 0.4 mm yr^?1) uplift and intense submarine erosion due to frequent hyperpycnal currents, whereas the eastern part bears more rapid (i.e. 0.7 mm yr^?1) uplift and has little or no hyperpycnal currents. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical simulation of tidal current and erosion and sedimentation in the Yangshan deep-water harbor of Shanghai

Yangshan near-shore sea area is the multi-island and multi-channel area with strong flow velocity and high suspended sediment concentration. Based on the characteristics of tidal currents, waves, and sediment in the Yangshan area, a two-dimensional numerical model of tidal currents, sediment transport, and sea bed deformation is developed. In the model, the effects of tidal currents and wind waves on sediment transport are considered. According to characteristics of the study area, unstructured grids are applied to fit the boundaries of the near-shore sea area. The results show that the calculated values are in good agreement with the measured data. The field of tidal currents, suspended sediment concentrations, and the deformation of the seabed can be successfully simulated.     

Interactions among waves,current, and mud: Numerical and laboratory studies

Interactions between waves, current, mud and turbulence are very complicated in the coastal and estuarine turbid waters. It is still necessary to improve our understanding of the fundamental physical processes governing the cohesive sediment transport in the coastal and estuarine waters. A numerical model is developed to study the interactions among waves, current, and mud. An eddy viscosity model for wave and current is proposed in order to close the equations of wave motion or of current motion in a combined flow, respectively. The equations of mud transport are derived based on the visco-elastic properties of mud. Coupling the equations of wave motion or of current motion for water layer with those of mud layer can give (1) wave height; (2) distributions of current velocities in the water layer; (3) distributions of transport velocities at the water–mud interface; and (4) distributions of mass transport velocities within the mud layer. These modeled results are in a reasonable agreement with experimental results. Results suggest that (1) the rate of wave attenuation increases in the opposing currents (currents against in the direction in which the waves propagate) and decreases in the following currents (currents in the same direction as that in which the waves propagate); (2) the opposing currents would have more significant effects on the rate of wave height attenuation than the following currents; (3) the effect of current on the rate of wave attenuation on the muddy bottom is larger than that on the rigid bottom; (4) mud transport rate increased in the following currents but decreased in the opposing currents; and (5) the rate of wave height attenuation on the mud bottom is one order of magnitude larger than that on the rigid bottom.     

First record of contourites from Lower Devonian Liptrap Formation in southeast Australia

Typical contourite deposits associated with submarine turbidite fan deposits are recognized for the first time from the Lower Devonian Liptrap Formation at Cape Liatrap, Victoria in southeast Australia. The contourites are well integrated within the turbidite fan deposits and are characterized by thin (5-8 cm), lenticular, well-sorted coarse-grained siltstones to fine-grained sandstones with current-ripples and cross beddings. The palaeocurrent directions of the turbidite fan and contourites are perpendicular to each other, with the former directed generally westward while the latter varying from 165° to 190° southward. In view of the facies types and architecture, we suggest that the turbidite fan was developed at the base of a westward inclined palaeo-slope, at the front of which the contourites were deposited as a result of southward flowing deep-sea contour (geostrophic) currents. The depositional setting inter- preted for the Liptrap Formation thus may provide a provisional model for the Lower Devonian conti- nental slope and abyssal basin environment in the southeastern part of the Melbourne Trough.     

Mechanisms of bottom boundary fluxes in a numerical model of the Shetland shelf

Across-slope bottom boundary layer (BBL) fluxes on the shelf-edge connect this region to deeper waters. Two proposed ways in which across-slope BBL fluxes can occur, in regions that have a slope current aligned to the bathymetry, are the frictional veering of bottom currents termed the ‘Ekman drain’ and through local wind-forced downwelling (wind-driven surface Ekman flow with an associated bottom flow). We investigate the variability, magnitude and spatial scale of BBL fluxes on the Shetland shelf, which has a prominent slope current, using a high-resolution (～2 km) configuration of the MITgcm model. Fluxes are analysed in the BBL at the shelf break near the 200 m isobath and are found to have a seasonal variability with high/low volume transport in winter/summer respectively. By using a multivariate regression approach, we find that the locally wind-driven Ekman transport plays no explicit role in explaining daily bottom fluxes. We can better explain the variability of the across-slope BBL flux as a linear function of the speed and across-slope component of the interior flow, corresponding to an Ekman plus mean-flow flux. We estimate that the mean-flow is a greater contributor than the Ekman flux to the BBL flux. The spatial heterogeneity of the BBL fluxes can be attributed to the mean-flow, which has a much shorter decorrelation length compared to the Ekman flux. We conclude that both the speed and direction of the interior current determines the daily BBL flux. The wind does not explicitly contribute through local downwelling, but may influence the interior current and therefore implicitly the BBL fluxes on longer timescales.     

Seasonal Variations of Currents over the Western Slope of the Middle Caspian Bed

Results of field observations of current dynamics in the frontal zone of the western Middle Caspian are given. The cyclonic circulation over the western slope in winter is shown to be a unidirectional intense current with velocities up to 100 cm/s. In summer, the current slows down and separates into branches—it turns southwestward and westward at the slope depth down to 150 m, southward and southeastward at the depth of ~100–350 m, and eastward at larger depths. In summer, shelf currents interact with the flow of Middle Caspian cyclonic circulation, resulting in that anticyclonic vortices reach the shelf.     

Coastal and shelf circulation in the vicinity of Camamu Bay (14°S), Eastern Brazilian Shelf

The Camamu Bay (CMB) is located on the narrowest shelf along the South American coastline and close to the formation of two major Western Boundary Currents (WBC), the Brazil/North Brazil Current (BC/NBC). These WBC flow close to the shelf break/slope region and are expected to interact with the shelf currents due to the narrowness of the shelf. The shelf circulation is investigated in terms of current variability based on an original data set covering the 2002-2003 austral summer and the 2003 austral autumn. The Results show that the currents at the shelf are mainly wind driven, experiencing a complete reversal between seasons due to a similar change in the wind field. Currents at the inner-shelf have a polarized nature, with the alongshore velocity mostly driven by forcings at the sub-inertial frequency band and the cross-shore velocity mainly supra-inertially forced, with the tidal currents playing an important role at this direction. The contribution of the forcing mechanisms at the mid-shelf changes between seasons. During the summer, forcings in the two frequency bands are important to drive the currents with a similar contribution of the tidal currents. On the other hand, during the autumn season, the alongshore velocity is mostly driven by sub-inertial forcings and tidally driven currents still remain important in both directions. Moreover, during the autumn when the stratification is weaker, the response of the shelf currents to the wind forcing presents a barotropic signature. The meso-scale processes related to the WBC flowing at the shelf/slope region also affect the circulation within the shelf, which contribute to cause significant current reversals during the autumn season. Currents at the shelf-estuary connection are clearly supra-inertially forced with the tidal currents playing a key role in the generation of the along-channel velocities. The sub-inertial forcings at this location act mainly to drive the weak ebb currents which were highly correlated with both local and remote wind forcing during the summer season.     

Reflection of equatorial Kelvin waves at eastern ocean boundaries Part II: Pacific and Atlantic Oceans

The effect of viscosity, non linearities, incident wave period and realistic eastern coastline geometry on energy fluxes are investigated using a shallow water model with a spatial resolution of 1/4 degree in both meridional and zonal directions. Equatorial and mid-latitude responses are considered. It is found that (1) the influence of the coastline geometry and the incident wave period is more important for the westward energy flux than for the poleward flux, and (2) the effect of the inclination of the eastern ocean boundary on the poleward energy flux, for the Pacific and Atlantic Oceans, decline as the period of the incident wave increases. Furthermore, the model simulations suggest that the poleward energy fluxes from meridional boundaries give plausible results for motions of seasonal and annual periods. For comparatively shorter periods, a realistic coastline geometry has to be included for more accurate results. It is recommended that any numerical model involving the reflection of baroclinic Rossby waves (of intraseasonal, seasonal or annual periods) on the eastern Pacific or Atlantic Oceans, should consider the effect of the coastline geometry in order to improve the accuracy of the results.     

Topographically generated steady currents in barotropic turbulence

Abstract

Steady currents develop in oceanic turbulence above topography even in the absence of steady forcing. Mesoscale steady currents are correlated with mesoscale topography with anticyclonic eddies above topographic bumps, and large scale westward flows develop when β is non-zero. The relationship between those two kinds of steady currents, as well as their dependence on various parameters, is studied using a barotropic quasi-geostrophic channel model. The percentage of steady energy is found to depend on the forcing, friction and topography in a non-monotonic fashion. For example, the percentage of steady currents grows with the energy level in the linear regime (low energies) and decreases when the energy level increases in the nonlinear regime (high energies). Mesoscale steady currents are the energy source for the steady westward flow U, and therefore U is the maximum when large scale and mesoscale currents are of the same order of magnitude. This happens when the ratio S of the large scale slope βH/f ₀ and the mesoscale rms topographic slope α is of order one. U decreases for both small and large values of S.     

Baroclinic dynamics of wind-driven circulation in a stratified bay: A numerical study using models of varying complexity

The baroclinic response of a stratified coastal embayment (Lunenburg Bay of Nova Scotia) to the observed wind forcing is examined using two numerical models. A linear baroclinic model based on the normal mode approach shows skill at reproducing the observed isotherm movements and sub-surface currents during a time of strong stratification in the bay. The linear model also shows that the isotherm movement in Lunenburg Bay is influenced by the wind forcing and propagation of baroclinic Kelvin waves from neighbouring Mahone Bay. The effects of nonlinearity and topography are investigated using a three-dimensional nonlinear coastal circulation model. The nonlinear model results demonstrate that the nonlinear advection terms generate a gyre circulation at the entrance of Lunenburg Bay, and the slope bottom topography at the mouth of the bay strengthens the sub-surface time-mean inflow on the southern side of the bay. A comparison of model-calculated currents in different numerical experiments clearly shows that baroclinicity plays a dominant role in the dynamics of wind-driven circulation in Lunenburg Bay.     

Effect of slope angle of an artificial pool on distributions of turbulence

Experiments were carried out over a 2-dimentional pool with a constant length of 1.5 m and four different slopes.The distributions of velocity,Reynolds stress and turbulence intensities have been studied in this paper.Results show that as flow continues up the exit slope,the flow velocity increases near the channel bed and decreases near the water surface.The flow separation was not observed by ADV at the crest of the bed-form.In addition,the length of the separation zone increases with the increasing of entrance and exit slopes.The largest slope angle causes the maximum normalized shear stress.Based on the experiments,it is concluded that the shape of Reynolds stress distribution is generally dependent on the entrance and exit slopes of the pool.Also,the shape of Reynolds stress distribution is affected by both decelerating and accelerating flows.Additionally,with the increase in the slope angle,secondary currents are developed and become more stable.Results of the quadrant analysis show that the momentum between flow and bed-form is mostly transferred by sweep and ejection events.     

Baroclinic and topographic influences on the transport in western boundary currents

Abstract

One of the central unsolved theoretical problems of the large scale ocean circulation is concerned with explaining the very large transports measured in western boundary currents such as the Gulf Stream and the Kuroshio. The only theory up to now that can explain the size of these transports is that of non-linear recirculation in which the advective terms in the momentum equations became important near the western boundary. In this paper an alternative explanation is suggested. When bottom topography and baroclinic effects are included in a wind-driven ocean model it is shown that the western boundary current can have a transport larger than that predicted from the wind stress distribution even when the nonlinear advective terms are ignored. The explanation lies in the presence of pressure torques associated with bottom topography which can contribute to the vorticity balance in the same sense as the wind stress curl.

Three numerical experiments have been carried out to explore the nature of this process using a three dimensional numerical model. The first calculation is done for a baroclinic ocean of constant depth, the second for a homogeneous ocean with an idealized continental slope topography, and the third for a baroclinic ocean with the same continental slope topography. The nature of the vorticity balance and of the circulation around closed paths is examined in each case, and it is shown that bottom pressure torques lead to enhanced transport in the western boundary current only for the baroclinic case with variable depth.     

On the meridional circulation and balance of momentum in the Southern Ocean of POP

 The circulation of the Southern Ocean is studied in the eddy-resolving model POP (Parallel Ocean Program) by an analysis of zonally integrated balances. The TEM formalism (Transformed Eulerian Mean) is extended to include topography and continental boundaries, thus deviations from a zonally integrated state involve transient and standing eddies. The meridional circulation is presented in terms of the Eulerian, eddy-induced, and residual streamfunctions. It is shown that the splitting of the meridional circulation into Ekman and geostrophic transports and the component induced by subgrid and Reynolds stresses is identical to a particular form of the zonally integrated balance of zonal momentum. In this balance, the eddy-induced streamfunctions represent the interfacial form stresses by transient and standing eddies and the residual streamfunction represents the acceleration of the zonal current by density fluxes in a zonally integrated frame. The latter acceleration term is directly related to the surface flux of density and interior fluxes due to the resolved and unresolved eddies. The eddy-induced circulation is extremely vigorous in POP. In the upper ocean a shallow circulation, reversed in comparison to the Deacon cell and mainly due to standing eddies, appears to the north of Drake Passage latitudes, and in the Drake Passage belt of latitudes a deep-reaching cell is induced by transient eddies. In the resulting residual circulation the Deacon cell is largely cancelled and the residual advection of the zonal mean potential density is balanced by diapycnal eddy and subgrid fluxes which are strong in the upper few hundred meters but small in the ocean interior. The balance of zonal momentum is consistent with other eddy-resolving models; a new aspect is the clear identification of density effects in the zonally integrated balance. We show that the wind stress and the stress induced by the residual circulation drive the eastward current, whereas both eddy species result in a braking. Finally, we extend the Johnson–Bryden model of zonal transport to incorporate all relevant terms from the zonal momentum balance. It is shown that wind stress and induction by the residual circulation carry an eastward transport while bottom form stress and the stress induced by standing eddies yield westward components of transport. Received: 26 June 2001 / Accepted: 2 November 2001     

Numerical prediction of turbidity currents on unsteady sediment-supply

A three-dimensional(3D) numerical model of unstable turbidity currents is developed based on the mechanism of sediment transport and turbulence theory.In this model,numerical simulation of turbidity currents without subsequent supply of muddy water was conducted using the same parameters as were used in the flume experiments.The evolution process of turbidity currents of completely losing supplies observed in the experiment was simulated by the model;validation of the numerical model and the algorithm was conducted.If momentarily interrupted process is regarded as a special case of the gradually interrupted,based on the preceding numerical simulation validity,it is feasible to simulate the motion law of turbidity currents under losing gradually supplies.By this method in this article,the characteristic of sediment-laden flow of losing gradually supplies was obtained,as well as its relationship between front velocity and sediment concentration.     

Understanding flow dynamics and density currents in a river-reservoir system under upstream reservoir releases

ABSTRACT

A three-dimensional flow and temperature model was applied for a 124 km river-reservoir system from Lewis Smith Dam tailrace to Bankhead Lock & Dam, Alabama. The model was calibrated against measured water levels, temperatures, velocities and flow rates from 4 May to 3 September 2011 under small constant release (2.83 m³/s) and large intermittent releases (~140 m³/s) from an upstream reservoir. Distributions of simulated flow and temperatures and particle tracking at various locations were analyzed which revealed the complex interactions of density currents, dynamic surface waves and solar heating. Flows in the surface and bottom layers moved in both upstream and downstream directions. If there was small constant release only from Smith Dam, simulated bottom temperatures at Cordova were on average 4.8°C higher than temperatures under actual releases. The momentum generated from large releases pushed bottom density currents downstream, but the released water took several days to reach Cordova.

Editor D. Koutsoyiannis; Associate editor B. Dewals