Near-surface circulation of the northeast Pacific Ocean derived from WOCE-SVP satellite-tracked drifters

Statistics of the near-surface circulation in the northeast Pacific Ocean were derived from the trajectories of nearly 100 surface drifters tracked between August 1990 and December 1995 as part of the World Ocean Circulation Experiment's (WOCE) Surface Velocity Program (SVP). Drifters were drogued within the mixed layer (15 m drogue depth) or near the top of the permanent halocline (120 m). All branches of the Alaskan Gyre were well-sampled at both depths, revealing a weak Subarctic Current, a bifurcation of the Subarctic Current near 48°N, 130°W at 15 m depth, and strong, variable flow in the Alaska Current and Alaskan Stream. At 120 m depth, northward flow in the Alaska Current occurred much farther offshore than within the mixed layer. The drifter trajectories revealed interannual variability, with evidence of an intensified Alaskan Gyre during the winters of 1991–92 and 1992–93 and more southerly transport during winter 1994–95. A minimum in eddy kinetic energy was found at both depths within the northern branch of the Subtropical Gyre. Eddy kinetic energies were nearly twice as high in the mixed layer compared to below, and were 2–3 times larger in winter than in summer throughout most of the near-surface Alaskan Gyre. High eddy energies observed near the eastern perimeter of the Alaskan Gyre may be due to the offshore intrusion of eddies formed by coastal current instabilities.Taylor's theory of single-particle dispersion was applied to the drifter ensembles to estimate Lagrangian decorrelation scales and eddy diffusivities. Both the initial dispersion and random walk regimes were identified in the dispersion time series computed for several regions of both ensembles. The integral time scales and eddy diffusivities computed from the dispersion scale linearly with r.m.s. velocity, which is consistent with drifter studies from the Atlantic. An exception is the meridional integral time scales, which were nearly constant throughout the study area and at both drogue depths. The magnitudes of the derived eddy statistics are comparable to those derived from surface drifters in other parts of the world ocean. These are the first Lagrangian estimates of particle dispersion over a broad region of the near-surface North Pacific, and the consistency of the results with previous studies from the Atlantic lends credence to the idea that the simplifying assumptions of Taylor (1921) (Proceedings of the London Mathematical Society Series A 20, 196–221) are reasonably valid throughout the upper ocean. This bodes well for the effective parameterization of near-surface diffusivities in general circulation models. Finally, the drifter-derived velocity statistics were used to speculate on the source regions of waters of possible coastal origin observed at offshore stations during the field studies of the Canadian Joint Global Ocean Flux Study.     

Lagrangian predictability assessed in the East China Sea

A group of 30 surface drifters, launched over a 4 day period as part of a US Navy exercise in early October 2007, are used to assess the predictability of trajectories in a confined geographic region at the northwestern edge of the Kuroshio north of Taiwan. Model trajectories were computed from archives of hourly hindcast velocities from the US Navy East Asian Seas (EAS16) model with 1/16° horizontal resolution. Three metrics are defined for comparing observed and modeled trajectories. All three metrics indicated that model trajectories separated from observations by roughly 15 km after the first 24 h on average. Because of the unique launch strategy for these drifters, with six repetitions of launches from four locations, the dependence of predictability on both launch time and launch location could be assessed separately. Predictive skill displayed only modest dependence on launch time, likely influenced by the passage of a typhoon near the experiment area a few days prior to the first drifter launch. Launch location was a much more reliable indicator of predictive skill, with trajectories for launches closest to the edge of the Kuroshio typically hardest to predict, and those for launches on the shelf, where currents tended to be weaker, predicted more accurately. Comparisons of skill metric statistics for modeled trajectories from hindcasts with and without tides suggested that tidal currents have only a small impact on predictive skill. The influence of archive time and space resolution was also studied using sets of model trajectories computed from hindcast archives that were systematically subsampled separately in space and time. Coarsening by up to a factor of eight in either space or time had little impact on predictive skill. Further coarsening degraded trajectory predictions, particularly when coarsening in time leads to an archive time step too large to adequately resolve the tides. While accurate trajectory predictions remain challenging for ocean models, skill assessments like the one presented here are important for developing error estimates for users of trajectory forecasts and for gaining new insight into potential sources of model errors.     

Agulhas leakage into the Atlantic estimated with subsurface floats and surface drifters

Surface drifters and subsurface floats drifting at depths near 800 m were used to study the pathways of warm, salty Indian Ocean water leaking into the South Atlantic that is a component of the upper limb of the Atlantic meridional overturning circulation (MOC). Four drifters and 5 floats drifted from the Agulhas Current directly into the Benguela Current. Others looped for various amounts of time in Agulhas rings and cyclones, which translated westward into the Atlantic, contributing a large part of Indian Ocean leakage. Agulhas rings translated into the Benguela Current, where they slowly decayed. Some large, blob-like Agulhas rings with irregular shapes were found in the southeastern Cape Basin. Drifter trajectories suggest these rings become more circular with time, eventually evolving into the circular rings observed west of the Walvis Ridge. Agulhas cyclones, which form on the north side of the Agulhas Current south of Africa, translated southwestward (to 6°E) and contributed water to the southern Cape Basin. A new discovery is a westward extension from the mean Agulhas retroflection measured by westward drifting floats near 41°S out to at least 5°W, with some floats as far west as 25°W. The Agulhas extension appears to split the South Atlantic Current (SAC) into two branches and to transport Agulhas water westward, where it is mixed and blended with eastward-flowing water from the western Atlantic. The blended mixture flows northeastward in the northern branch of the SAC and into the Benguela Current. Agulhas leakage transport was estimated from drifters and floats to be at least 15 Sv in the upper 1000 m, which is equivalent to the transport of the upper layer MOC. It is suggested that the major component of the upper layer overturning circulation in the Atlantic is Agulhas leakage in the form of Agulhas rings.     

Ocean current estimation using a Multi-Model Ensemble Kalman Filter during the Grand Lagrangian Deployment experiment (GLAD)

In the summer and fall of 2012, during the GLAD experiment in the Gulf of Mexico, the Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE) used several ocean models to assist the deployment of more than 300 surface drifters. The Navy Coastal Ocean Model (NCOM) at 1 km and 3 km resolutions, the US Navy operational NCOM at 3 km resolution (AMSEAS), and two versions of the Hybrid Coordinates Ocean Model (HYCOM) set at 4 km were running daily and delivering 72-h range forecasts. They all assimilated remote sensing and local profile data but they were not assimilating the drifter’s observations. This work presents a non-intrusive methodology named Multi-Model Ensemble Kalman Filter that allows assimilating the local drifter data into such a set of models, to produce improved ocean currents forecasts. The filter is to be used when several modeling systems or ensembles are available and/or observations are not entirely handled by the operational data assimilation process. It allows using generic in situ measurements over short time windows to improve the predictability of local ocean dynamics and associated high-resolution parameters of interest for which a forward model exists (e.g. oil spill plumes). Results can be used for operational applications or to derive enhanced background fields for other data assimilation systems, thus providing an expedite method to non-intrusively assimilate local observations of variables with complex operators. Results for the GLAD experiment show the method can improve water velocity predictions along the observed drifter trajectories, hence enhancing the skills of the models to predict individual trajectories.     

Year-long float trajectories in the Labrador Sea Water of the eastern North Atlantic Ocean

Year-long Lagrangian trajectories within the Labrador Sea Water of the eastern North Atlantic Ocean are analysed for basic flow statistics. Root-mean-square velocities at 1750 m depth are about 2 cm/s, except within the North Atlantic Current, where they are twice as large. These values are consistent with previous Eulerian measurements and extend those results to a much larger domain of the eastern basin. Mean flow estimates in boxes large enough to contain about 1 float-year of data indicate that Labrador Sea Water, having crossed the Mid- Atlantic Ridge (not resolved) near 50–55^°N, presumably with the North Atlantic Current, partially recirculates to the north in the subpolar gyre, as well as entering the subtropical gyre and continuing south and west. The circulation of this water mass, as defined by the 1 yr average velocities, is stronger than traditional models of deep circulation would suggest, with an interior flow of roughly 1 cm/s. Mean speeds up to 3 cm/s were observed, with the highest values near the Azores Plateau. North of 45°N–55°N, mean eastward speeds closer to 0.2 cm/s were observed. Wind-generated barotropic fluctuations may be responsible for some part of the transport at this depth.     

Surface circulation in the southwestern Japan/East Sea as observed from drifters and sea surface height

The mean circulation of the surface layer of the southwestern Japan/East Sea (JES) was examined using current measurements collected at 15 m by satellite-tracked drifters and merged sea level anomalies from satellite altimeters. The study of circulation patterns in this paper focused on the inflow passing through the western channel of the Korea Strait from the East China Sea. Empirical Orthogonal Function (EOF) analysis of non-seasonal sea level anomalies revealed that significant energy in the circulation pattern of Ulleung Basin was controlled by the inflow conditions through the Korea Strait. Three circulation patterns were identified that depended on the initial relative vorticity of the inflow. When inflow had initially large negative vorticity, the flow gained more negative vorticity due to deepening of the bottom (stretching) and then turned right after entering the JES. The inflow then followed the path of the Tsushima Warm Current along the coast of Japan. When the inflow was strong, with a speed in excess of 55 cm/s and with a large positive vorticity, potential vorticity appeared to be conserved. In this case, the EKWC followed isobaths along the coast and then left the coast, following topographic features north of Ulleung-Do. The northward flowing jet developed inertial meandering after leaving the coast, which is a characteristic of many western boundary currents. The regular, bimonthly deployments of drifters in the western portion of the Korea Strait revealed that splitting or branching of the flow through the western channel of the Korea Strait occurred only 15% of the time. And splitting or branching rarely occurred during the fall and winter seasons, when the inflow splitting was previously reported in hydrographic surveys. The time-averaged circulation map of the EKWC and its seaward extension were considerably enhanced by using regularly sampled geostrophic velocities calculated from sea level anomalies to remove biases in the mean velocity that were caused by irregular spatial and temporal drifter observations. The East Korean Warm Current, a mean coastal current along the Korean coast, behaved like the simple model by Arruda et al. (2004) in which the generation of the Ulleung Warm Eddy and the meandering circulation pattern were well reproduced.     

A numerical study of nonlinear wave run-up on a vertical plate

A finite difference model based on a recently derived highly-accurate Boussinesq-type formulation is presented. Up to the third-order space derivatives in terms of the velocity variables are retained, and the horizontal velocity variables are re-formulated in terms of a velocity potential. This decreases the total number of unknowns in two horizontal dimensions from seven to five, simplifying the implementation, and leading to increased computational efficiency. Analysis of the embedded properties demonstrates that the resulting model has applications with errors of 2 to 3% for (wavenumber times depth) kh ≤ 10 in terms of dispersion and kh ≤ 4 in terms of internal kinematics. The stability and accuracy of the discrete linearised systems are also analysed for both potential and velocity formulations and the advantages and disadvantages of each are discussed. The velocity potential model is then used to study physically demanding problems involving highly nonlinear wave run-up on a bottom-mounted (surface-piercing) plate. New cases involving oblique incidence are considered. In all cases, comparisons with recent physical experiments demonstrate good quantitative accuracy, even in the most demanding cases, where the local wave steepness can exceed (waveheight divided by wavelength) H / L = 0.20. The velocity potential model is additionally shown to have numerical advantages when dealing with wave–structure interactions, requiring less smoothing around exterior structural corners.     

The swarm dynamics of northern krill (Meganyctiphanes norvegica) and pteropods (Cavolinia inflexa) during vertical migration in the Ligurian Sea observed by an acoustic Doppler current profiler

A ship-mounted 153 kHz narrow-band ADCP and 1 m² MOCNESS were deployed between 16 and 24 Sept. 1997 in the Ligurian central zone (∼43°20′N 7°48′E). Results from both instruments showed that the zooplankton community performed vertical migrations that conformed to the classical pattern of ascent at dusk (∼18:30 h) and descent at dawn (∼06:30 h). Depth-discrete net samples between 0 and 500 m showed that the community was dominated by two species, the euphausiid Meganyctiphanes norvegica (Northern krill) and the pteropod Cavolinia inflexa, which migrated in separate discrete bands that were detectable by the ADCP. Information from the ADCP was used to estimate vertical migration speed in two ways: (i) from the trajectory of the back-scattering bands over time and (ii) from the Doppler-shift vertical velocity measured within depth zones at the corresponding time and depth of these bands. Estimates of the migration speed of C. inflexa were between 2 and 7 cm s⁻¹ upwards and between 4 and 7 cm s⁻¹ downwards. M. norvegica was estimated to migrate at speeds between 7 and 8 cm s⁻¹ upwards and over 11 cm s⁻¹ downwards. The consistently lower migration speeds estimated from Doppler measurements as compared with estimates obtained from measuring trajectories of back-scattering bands over time was believed to result from a methodological artefact. The Doppler measurements were nevertheless useful in a relative sense in revealing the relative speed of individuals within swarms. It was shown that individuals at the front of the upwardly migrating band of M. norvegica moved more slowly than those at the rear. These results illustrate the extra biological information that can be obtained by ADCPs compared with conventional echo-sounders.     

Flux comparison of Eulerian and Lagrangian estimates of Agulhas leakage: A case study using a numerical model

Estimating the magnitude of Agulhas leakage, the volume flux of water from the Indian to the Atlantic Ocean, is difficult because of the presence of other circulation systems in the Agulhas region. Indian Ocean water in the Atlantic Ocean is vigorously mixed and diluted in the Cape Basin. Eulerian integration methods, where the velocity field perpendicular to a section is integrated to yield a flux, have to be calibrated so that only the flux by Agulhas leakage is sampled. Two Eulerian methods for estimating the magnitude of Agulhas leakage are tested within a high-resolution two-way nested model with the goal to devise a mooring-based measurement strategy. At the GoodHope line, a section halfway through the Cape Basin, the integrated velocity perpendicular to that line is compared to the magnitude of Agulhas leakage as determined from the transport carried by numerical Lagrangian floats. In the first method, integration is limited to the flux of water warmer and more saline than specific threshold values. These threshold values are determined by maximizing the correlation with the float-determined time series. By using the threshold values, approximately half of the leakage can directly be measured. The total amount of Agulhas leakage can be estimated using a linear regression, within a 90% confidence band of 12 Sv. In the second method, a subregion of the GoodHope line is sought so that integration over that subregion yields an Eulerian flux as close to the float-determined leakage as possible. It appears that when integration is limited within the model to the upper 300 m of the water column within 900 km of the African coast the time series have the smallest root-mean-square difference. This method yields a root-mean-square error of only 5.2 Sv but the 90% confidence band of the estimate is 20 Sv. It is concluded that the optimum thermohaline threshold method leads to more accurate estimates even though the directly measured transport is a factor of two lower than the actual magnitude of Agulhas leakage in this model.     

Optimisation of an idealised ocean model,stochastic parameterisation of sub-grid eddies

An optimisation scheme is developed to accurately represent the sub-grid scale forcing of a high dimensional chaotic ocean system. Using a simple parameterisation scheme, the velocity components of a 30 km resolution shallow water ocean model are optimised to have the same climatological mean and variance as that of a less viscous 7.5 km resolution model. The 5 day lag-covariance is also optimised, leading to a more accurate estimate of the high resolution response to forcing using the low resolution model.The system considered is an idealised barotropic double gyre that is chaotic at both resolutions. Using the optimisation scheme, we find and apply the constant in time, but spatially varying, forcing term that is equal to the time integrated forcing of the sub-grid scale eddies. A linear stochastic term, independent of the large-scale flow, with no spatial correlation but a spatially varying amplitude and time scale is used to represent the transient eddies. The climatological mean, variance and 5 day lag-covariance of the velocity from a single high resolution integration is used to provide an optimisation target. No other high resolution statistics are required. Additional programming effort, for example to build a tangent linear or adjoint model, is not required either.The focus of this paper is on the optimisation scheme and the accuracy of the optimised flow. However the forcing can provide insights in the design of deterministic and stochastic parameterisations. In the present study, we found that the stochastic parameterisation correcting the model variance is associated with the spatial pattern of eddy-decorrelation timescales rather than the spatial pattern of the amplitude of the variance. The method can be applied in future investigations into the physical processes that govern barotropic turbulence and it can perhaps be applied to help understand and correct biases in the mean and variance of a more realistic coarse or eddy-permitting ocean model. The method is complementary to current parameterisations and can be applied at the same time without modification.     

Hydrographic and hydrodynamic characteristics of the eastern basin of the Bransfield Strait (Antarctica)

Hydrographic station and current meter data are used to estimate circulation and transport in the eastern basin of the Bransfield Strait. The short distance between adjacent hydrographic stations (20 km) allows evaluation of structures at scales seldom addressed in previous studies. The main feature of the derived circulation is the Bransfield Front and its associated baroclinic jet (the Bransfield Current). This frontal current crosses the northern half of the basin in a generally SW–NE direction, has maximum geostrophic speeds of 22 cm s^−l (at the jet entrance), and has geostrophic transport relative to 500 dbar estimated to be 1 Sv. Dynamically significant mesoscale features associated with the Bransfield Current are seen to be relevant down to 500 dbar. Specific aspects inferred from our analysis are the apparent high degree of stationarity of the described circulation, the shallow intrusions of Circumpolar Deep Water through the northern boundary of the domain (from the Drake Passage), and the northward sinking of Weddell Sea water over most of the domain.     

On the nature of oceanic eddies shed by the Island of Gran Canaria

We use hydrographic and buoy data to compare the initial temperature fields and Lagrangian evolution of water parcels in two vortices generated by the southward flowing Canary Current passing around the island of Gran Canaria Island. One vortex is anticyclonic, shed in June 1998 as the result of an incident current of about 0.05 m s⁻¹, and the second one is cyclonic, shed in June 2005 with the impinging current estimated as 0.03 m s⁻¹. The two vortices exhibit contrasting characteristics yet display some important similarities. The isopycnals are depressed in the core of the anticyclonic vortex, at least down to a depth of 700 m, whilst they dome up in the core of the cyclonic vortex but only down to 450 m. In the top 300 m the depression/doming of the isotherms is similar for both vortices, with a maximum vertical displacement of the isotherm of about 80 m, which correspond to temperature anomalies of some 2.5 °C at a given depth. A simple method is developed to obtain the initial orbital velocity field from the temperature data, from which we estimate peak values of 0.7 and 0.5 m s⁻¹ for the anticyclonic and cyclonic vortices, respectively. The buoys, three for the anticyclonic vortex and two for the cyclonic one, were drougued at 100 m depth, below the surface mixed layer, and their initial velocities are consistent with the above values. In both vortices, the buoys revolve either within a central core, where the rotation rate remains stable and large for several weeks, or in an outer ring, where the rotation rate is significantly smaller and displays large radial fluctuations. Within the inner core the anticyclonic vortex has significant inward radial velocity, while the cyclonic vortex has near-zero radial mean motions. The cyclonic vortex rotates more slowly than the anticyclonic, their initial periods being 4.5 and 2.5 days, respectively. A simple axisymmetric model with radial diffusion (coefficient K_h≅25 m² s⁻¹) and advection reproduces the observations reasonably well, the diffusive effect being more important than that resulting from the observed radial advection. The model also supports the hypothesis that the rotation rate of cyclonic vortices is less than that of anticyclonic vortices, as otherwise they would become inertially unstable. Both the buoys data and sea surface temperature images confirm that the vortices evolve from youth to maturity, as the cores shrink and the outer rings expands, and then to a decay stage, as the core rotation rates decrease, though frequent interactions with other mesoscale structures result in more accelerated aging. Despite these interaction they last many months as coherent structures south of the Canary Islands.     

Hydrography and circulation of the West Antarctic Peninsula Continental Shelf

The water mass structure and circulation of the continental shelf waters west of the Antarctic Peninsula are described from hydrographic observations made in March–May 1993. The observations cover an area that extends 900 km alongshore and 200 km offshore and represent the most extensive hydrographic data set currently available for this region. Waters above 100–150 m are composed of Antarctic Surface Water and its end member Winter Water. Below the permanent pycnocline is a modified version of Circumpolar Deep Water, which is a cooled and freshened version of Upper Circumpolar Deep Water. The distinctive signature of cold and salty water from the Bransfield Strait is found at some inshore locations, but there is little indication of significant exchange between Bransfield Strait and the west Antarctic Peninsula shelf. Dynamic topography at 200 m relative to 400 m indicates that the baroclinic circulation on the shelf is composed of a large, weak, cyclonic gyre, with sub-gyres at the northeastern and southwestern ends of the shelf. The total transport of the shelf gyre is 0.15 Sv, with geostrophic currents of order 0.01 m s^-1. A simple model that balances across-shelf diffusion of heat and salt from offshore Upper Circumpolar Deep Water with vertical diffusion of heat and salt across the permanent pycnocline into Winter Water is used to explain the formation of the modified Circumpolar Deep Water that is found on the shelf. Model results show that the observed thermohaline distributions across the shelf can be maintained with a coefficient of vertical diffusion of 10^-4 m² s^-1 and horizontal diffusion coefficients for heat and salt of 200 and 1200 m² s^-1, respectively. When the effects of double diffusion are included in the model, the required horizontal diffusion coefficients for heat and salt are 200 and 400 m² s^-1, respectively.     

Float trajectories in the deep western boundary current and deep equatorial jets of the tropical Atlantic

Fourteen neutrally buoyant SOFAR floats at a nominal depth of 1800 m were tracked acoustically for 3.7 yr in the vicinity of the western boundary and the equator of the Atlantic Ocean. The trajectories revealed a swift, narrow, southward-flowing deep western boundary current (DWBC) extending from 7N across the equator. Two floats crossed the equator in the DWBC and went to 10S. Two other floats left the DWBC and drifted eastward in the equatorial band (3S–3N). Three floats entered the DWBC from the equatorial current system and drifted southward. These results suggest that at times the DWBC flows directly southward across the equator with a mean velocity of 8–9 cm/s averaged over long distances (∼2800 km). At other times DWBC water is diverted eastward near the equator for long periods (2–3 yr), which can reduce the mean along-boundary velocity to 1–2 cm/s. This is much less than the instantaneous along-boundary velocities in the DWBC, which are often above 25 cm/s and occasionally exceed 50 cm/s. Mean eastward-flowing jets were observed near 2N and 2S bounding a mean westward jet centered on the equator (1S–1N). The southern jet at 2S coincides with a CFC-rich plume centered south of the equator. The CFC plume is inferred to have been advected by the southern jet across the Atlantic and into the Gulf of Guinea.     

Tuning a physically-based model of the air–sea gas transfer velocity

Air–sea gas transfer velocities are estimated for one year using a 1-D upper-ocean model (GOTM) and a modified version of the NOAA–COARE transfer velocity parameterization. Tuning parameters are evaluated with the aim of bringing the physically based NOAA–COARE parameterization in line with current estimates, based on simple wind-speed dependent models derived from bomb-radiocarbon inventories and deliberate tracer release experiments. We suggest that A = 1.3 and B = 1.0, for the sub-layer scaling parameter and the bubble mediated exchange, respectively, are consistent with the global average CO₂ transfer velocity k. Using these parameters and a simple 2nd order polynomial approximation, with respect to wind speed, we estimate a global annual average k for CO₂ of 16.4 ± 5.6 cm h^?1 when using global mean winds of 6.89 m s^?1 from the NCEP/NCAR Reanalysis 1 1954–2000. The tuned model can be used to predict the transfer velocity of any gas, with appropriate treatment of the dependence on molecular properties including the strong solubility dependence of bubble-mediated transfer. For example, an initial estimate of the global average transfer velocity of DMS (a relatively soluble gas) is only 11.9 cm h^?1 whilst for less soluble methane the estimate is 18.0 cm h^?1.     

Observations on the Rim Current structure,CIW formation and transport in the western Black Sea

CTD and ADCP measurements together with a sequence of satellite images indicate pronounced current meandering and eddy activity in the western Black Sea during April 1993. The Rim Current is identified as a well-defined meandering jet stream confined over the steepest topographic slope and associated cyclonic–anticyclonic eddy pairs located on both its sides. It has a form of highly energetic and unstable flow system, which, as it propagates cyclonically along the periphery of the basin, is modified in character. It possesses a two-layer vertical structure with uniform upper layer speed in excess of 50 cm/s (maximum value ∼100 cm/s), followed by a relatively sharp change across the pycnocline (between 100 and 200 m) and the uniform sub-pycnocline currents of 20 cm/s (maximum value ∼40 cm/s) observed up to the depth of ∼350 dbar, being the approximate limit of ADCP measurements. The cross-stream velocity structure exhibits a narrow core region (∼30 km), flanked by a narrow zone of anticyclonic shear on its coastal side and a broader region of cyclonic shear on its offshore side. The northwestern shelf circulation is generally decoupled from the influence of the basinwide circulation and is characterized by much weaker currents, less than 10 cm/s. The southward coastal flow associated with the Danube and Dinepr Rivers is weak during the measurement period and is restricted to a very narrow coastal zone.The data suggest the presence of temperature-induced overturning prior to the measurements, and subsequent formation of the Cold Intermediate Water mass (CIW) within the Northwestern Shelf (NWS) and interior of the western basin. The newly formed shelf CIW is transported in part along the shelf by the coastal current system, and in part it flows downslope across the shelf and intrudes into the Rim Current convergence zone. A major part of the cold water mass, however, seems to be trapped within the northwestern shelf. The CIW mass, injected into the Rim Current zone from the shelf and the interior region, is then circulated around the basin.     

Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.     

Direct observation of subtropical mode water circulation in the western North Atlantic Ocean

Eighteen Degree Water (EDW) is the dominant subtropical mode water of the North Atlantic subtropical gyre and is hypothesized as an interannual reservoir of anomalous heat, nutrients and CO₂. Although isolated beneath the stratified upper-ocean at the end of each winter, EDW may re-emerge in subsequent years to influence mixed layer properties and consequently air–sea interaction and primary productivity. Here we report on recent quasi-Lagrangian measurements of EDW circulation and stratification in the western subtropical gyre using an array of acoustically-tracked, isotherm-following, bobbing profiling floats programmed to track and intensively sample the vertically homogenized EDW layer and directly measure velocity on the 18.5 °C isothermal surface.The majority of the CLIVAR Mode Water Dynamics Experiment (CLIMODE) bobbers drifted within the subtropical gyre for 2.5–3.5 years, many exhibiting complex looping patterns indicative of an energetic eddy field. Bobber-derived Lagrangian integral time and length scales (3 days, 68 km) associated with motion on 18.5 °C were consistent with previous measurements in the Gulf Stream extension region and fall between previous estimates at the ocean surface and thermocline depth. Several bobbers provided evidence of long-lived submesoscale coherent vortices associated with substantial EDW thickness. While the relative importance of such vortices remains to be determined, our observations indicate that these features can have a profound effect on EDW distribution. EDW thickness (defined using a vertical temperature gradient criterion) exhibits seasonal changes in opposition to a layer bounded by the 17 °C and 19 °C isotherms. In particular, EDW thickness is generally greatest in winter (as a result of buoyancy-forced convection), while the 17°–19 °C layer is thickest in summer consistent with seasonal Ekman pumping. Contrary to previous hypotheses, the bobber data suggest that a substantial fraction of subducted EDW is isolated from the atmosphere for periods of less than 24 months. Seasonal-to-biennial re-emergence (principally within the recirculation region south of the Gulf Stream) appears to be a common scenario which should be considered when assessing the climatic and biogeochemical consequences of EDW.     

Diagnosis of vertical velocities with the QG omega equation: an examination of the errors due to sampling strategy

Vertical motion at the mesoscale plays a key role in ocean circulation, ocean-atmosphere interaction, and hence climate. It is not yet possible to make direct Eulerian measurements of vertical velocities less than 1000 m day⁻¹. However, by assuming quasi-geostrophic (QG) balance, vertical velocities O (10 m day⁻¹) can be diagnosed from the geostrophic velocity field and suitable boundary conditions. Significant errors in the accuracy of this diagnosis arise from the necessary compromise between spatial resolution and synopticity of a hydrographic survey. This problem has been addressed by sampling the output of a numerical ocean model to simulate typical oceanographic surveys of mesoscale fronts. The balance between the number of observations and the synopticity of observations affects the apparent flow and in particular the diagnosed vertical motion. A combination of effects can typically lead to errors of 85% in the estimation of net vertical heat flux. An analytical two-layer model is used to understand components of this error and indicate the key parameters for the design of mesoscale sampling.     

The pattern of cross-slope depositional fluxes

A simple model with horizontal and vertical diffusivities and settling velocity is used to calculate expected distributions of suspended particulate matter in a section across the continental shelf and slope. Dependencies on the shelf and slope profile, diffusivities, settling velocity, cross-slope advection and boundary sources/sinks are explored. It is found that the strongest factors are relative values of diffusivities and settling velocity, and the distribution of sources and sinks – including bottom deposition or resuspension. The latter is the principal means whereby an increased concentration near the bottom is likely, and is suggested as the usual reason for increased deposition recorded by sediment traps nearer the bottom. Observed thin, near-horizontal intermediate nepheloid layers put bounds on the vertical diffusivity and settling velocity, e.g. O(10^-4 m² s^-1, 10^-5 m s^-1) over Goban Spur in OMEX.