We introduce a simple variational principle for coherent material vortices in two-dimensional turbulence. Vortex boundaries are sought as closed stationary curves of the averaged Lagrangian strain. Solutions to this problem turn out to be mathematically equivalent to photon spheres around black holes in cosmology. The fluidic photon spheres satisfy explicit differential equations whose outermost limit cycles are optimal Lagrangian vortex boundaries. As an application, we uncover super-coherent material eddies in the South Atlantic, which yield specific Lagrangian transport estimates for Agulhas rings.

Reliable forecasts for the dispersion of oceanic contamination are important for coastal ecosystems, society, and the economy as evidenced by the Deepwater Horizon oil spill in the Gulf of Mexico in 2010 and the Fukushima nuclear plant incident in the Pacific Ocean in 2011. Accurate prediction of pollutant pathways and concentrations at the ocean surface requires understanding ocean dynamics over a broad range of spatial scales. Fundamental questions concerning the structure of the velocity field at the submesoscales (100 m to tens of kilometers, hours to days) remain unresolved due to a lack of synoptic measurements at these scales. Using high-frequency position data provided by the near-simultaneous release of hundreds of accurately tracked surface drifters, we study the structure of submesoscale surface velocity fluctuations in the Northern Gulf of Mexico. Observed two-point statistics confirm the accuracy of classic turbulence scaling laws at 200-m to 50-km scales and clearly indicate that dispersion at the submesoscales is local, driven predominantly by energetic submesoscale fluctuations. The results demonstrate the feasibility and utility of deploying large clusters of drifting instruments to provide synoptic observations of spatial variability of the ocean surface velocity field. Our findings allow quantification of the submesoscale-driven dispersion missing in current operational circulation models and satellite altimeter-derived velocity fields.T he Deepwater Horizon (DwH) incident was the largest accidental oil spill into marine waters in history with some 4.4 million barrels released into the DeSoto Canyon of the northern Gulf of Mexico (GoM) from a subsurface pipe over ∼84 d in the spring and summer of 2010 (1). Primary scientific questions, with immediate practical implications, arising from such catastrophic pollutant injection events are the path, speed, and spreading rate of the pollutant patch. Accurate prediction requires knowledge of the ocean flow field at all relevant temporal and spatial scales. Whereas ocean general circulation models were widely used during and after the DwH incident (2-6), such models only capture the main mesoscale processes (spatial scale larger than 10 km) in the GoM. The main factors controlling surface dispersion in the DeSoto Canyon region remain unclear. The region lies between the mesoscale eddy-driven deep water GoM (7) and the winddriven shelf (8) while also being subject to the buoyancy input of the Mississippi River plume during the spring and summer months (9). Images provided by the large amounts of surface oil produced in the DwH incident revealed a rich array of flow patterns (10) showing organization of surface oil not only by mesoscale straining into the loop current "Eddy Franklin," but also by submesoscale processes. Such processes operate at spatial scales and involve physics not currently captured in operational circulation models. Submesoscale motions, where they exist, can directly influence the local transport of biogeochemical tracers (11, 12) ...

No abstract

Mesoscale oceanic eddies are routinely detected from instantaneous velocities derived from satellite altimetry data. While simple to implement, this approach often gives spurious results and hides true material transport. Here it is shown how geodesic transport theory, a recently developed technique from nonlinear dynamical systems, uncovers eddies objectively. Applying this theory to altimetry-derived velocities in the South Atlantic reveals, for the first time, Agulhas rings that preserve their material coherence for several months, while ring candidates yielded by other approaches tend to disperse or leak within weeks. These findings suggest that available velocity-based estimates for the Agulhas leakage, as well as for its impact on ocean circulation and climate, need revision.

No abstract

[1] We demonstrate the feasibility of using dynamical systems tools to unambiguously identify mesoscale oceanic eddies from surface ocean currents derived using climatological hydrography and altimetry. Specifically, our analysis is based on extracting Lagrangian coherent structures (LCSs) from finite-time Lyapunov exponent (FTLE) fields. The FTLE fields reveal with unprecedented detail an intricate tangle of LCSs, which are hidden in ocean surface topography maps but sometimes are apparent in ocean color images. These LCSs delineate fluid domains with very different advective properties, and thus their detection provides an objective (i.e., frame-independent) means of identifying eddy boundaries. The importance of considering LCSs in quantifying transport by eddies is highlighted. Such a quantification does not rely on the common assumption-which is shown to be generally not valid-that transport is largely effected by the trapping and subsequent translation of water slugs inside eddies defined as the regions enclosed by sea height (streamfunction) contours within which rotation dominates over strain. LCSs are calculated for the whole globe and compared with satellite-tracked drogue drifter trajectories within a selected region of the South Atlantic. Citation: Beron-Vera, F. J., M. J.Olascoaga, and G. J. Goni (2008), Oceanic mesoscale eddies as revealed by Lagrangian coherent structures, Geophys. Res. Lett., 35, L12603,

[1] Analysis of drifter trajectories in the Gulf of Mexico has revealed the existence of a region on the southern portion of the West Florida Shelf (WFS) that is not visited by drifters that are released outside of the region. This socalled ''forbidden zone'' (FZ) suggests the existence of a persistent cross-shelf transport barrier on the southern portion of the WFS. In this letter a year-long record of surface currents produced by a Hybrid-Coordinate Ocean Model simulation of the WFS is used to compute Lagrangian coherent structures (LCSs), which reveal the presence of a persistent cross-shelf transport barrier in approximately the same location as the boundary of the FZ. The location of the cross-shelf transport barrier undergoes a seasonal oscillation, being closer to the coast in the summer than in the winter. A month-long record of surface currents inferred from high-frequency (HF) radar measurements in a roughly 60 km  80 km region on the WFS off Tampa Bay is also used to compute LCSs, and these also reveal the presence of transient transport barriers. While the HF-radarderived transport barriers cannot be unambiguously linked to the boundary of the FZ, this analysis does demonstrate the feasibility of monitoring transport barriers on the WFS using a HF-radar-based measurement system. The implications of a persistent cross-shelf transport barrier on the WFS for the development of harmful algal blooms on the shoreward side of the barrier are considered.

A framework for the study of surface ocean inertial particle motion is built from the Maxey-Riley set. A new set is obtained by vertically averaging each term of the original set, adapted to account for Earth's rotation effects, across the extent of a sufficiently small spherical particle that floats at an assumed unperturbed airsea interface with unsteady nonuniform winds and ocean currents above and below, respectively. The inertial particle velocity is shown to exponentially decay in time to a velocity that lies close to an average of seawater and air velocities, weighted by a function of the seawater-to-particle density ratio. Such a weighted average velocity turns out to fortuitously be of the type commonly discussed in the search-and-rescue literature, which alone cannot explain the observed role of anticyclonic mesoscale eddies as traps for marine debris or the formation of great garbage patches in the subtropical gyres, phenomena dominated by finite-size effects. A heuristic extension of the theory is proposed to describe the motion of nonspherical particles by means of a simple shape factor correction, and recommendations are made for incorporating wave-induced Stokes drift, and allowing for inhomogeneities of the carrying fluid density. The new Maxey-Riley set outperforms an ocean adaptation that ignored wind drag effects and the first reported adaption that attempted to incorporate them.