Showing 1–45 of 45 results for author: Beron-Vera, F J

Dynamics of inertial particles on the ocean surface with unrestricted reserve buoyancy

Authors: F. J. Beron-Vera

Abstract: The purpose of this note is to present an enhancement to a Maxey-Riley theory proposed in recent years for the dynamics of inertial particles on the ocean surface. This model upgrade removes constraints on the reserve buoyancy, defined as the fraction of the particle volume above the ocean surface. The refinement results in an equation that correctly describes both the neutrally buoyant and fully… ▽ More The purpose of this note is to present an enhancement to a Maxey-Riley theory proposed in recent years for the dynamics of inertial particles on the ocean surface. This model upgrade removes constraints on the reserve buoyancy, defined as the fraction of the particle volume above the ocean surface. The refinement results in an equation that correctly describes both the neutrally buoyant and fully buoyant particle scenarios. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Submitted to Physics of Fluids

Properties and baroclinic instability of stratified thermal upper-ocean flow

Authors: F. J. Beron-Vera, M. J. Olascoaga

Abstract: We study the properties of, and perform a stability analysis of a baroclinic zonal current in, a thermal rotating shallow-water model, sometimes called \emph{Ripa's model}, featuring stratification for quasigeostrophic upper-ocean dynamics. The model has Lie--Poisson Hamiltonian structure. In addition to Casimirs, the model supports integrals of motion that neither form the kernel of the bracket n… ▽ More We study the properties of, and perform a stability analysis of a baroclinic zonal current in, a thermal rotating shallow-water model, sometimes called \emph{Ripa's model}, featuring stratification for quasigeostrophic upper-ocean dynamics. The model has Lie--Poisson Hamiltonian structure. In addition to Casimirs, the model supports integrals of motion that neither form the kernel of the bracket nor are related to explicit symmetries. The model sustains Rossby waves and a neutral model, whose spurious growth is prevented by a positive-definite integral, quadratic on the deviation from the motionless state. A baroclinic zonal jet with vertical curvature is found to be spectrally stable for specific configurations of the gradients of layer thickness, vertically averaged buoyancy, and buoyancy frequency. Only a subset of such states was found Lyapunov stable using the available integrals, except the newly reported ones, whose role in constraining stratified thermal flow remains elusive. The existence of Lyapunov-stable states enabled bounding the nonlinear growth of perturbations to spectrally unstable states. Our results do not support the generality of earlier numerical evidence on the suppression of submesoscale wave activity as a result of the inclusion of stratification in thermal shallow-water theory, which we supported with direct numerical simulations. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: Submitted to Proceeding of the Royal Society A

On a priori bounding the growth of thermal instability waves

Authors: F. J. Beron-Vera

Abstract: We have previously shown that the nonlinear growth of a finite-amplitude perturbation to a basic state given by a baroclinic zonal flow on the $β$-plane in a thermal quasigeostrophic reduced-gravity model can be a priori bounded. In this note we show that, unlike we stated earlier, Lyapunov stability can be proved even when buoyancy varies linearly with the meridional coordinate. In addition to re… ▽ More We have previously shown that the nonlinear growth of a finite-amplitude perturbation to a basic state given by a baroclinic zonal flow on the $β$-plane in a thermal quasigeostrophic reduced-gravity model can be a priori bounded. In this note we show that, unlike we stated earlier, Lyapunov stability can be proved even when buoyancy varies linearly with the meridional coordinate. In addition to rectifying our earlier results we expand them by deriving an instability saturation bound by making use of the existence of such a class of Lyapunov-stable basic states. This bound can be smaller than that one we estimated before, reinforcing our previous conclusions. We also present a numerical test of the accuracy of the derived bound. △ Less

Submitted 19 March, 2024; v1 submitted 20 February, 2024; originally announced February 2024.

Comments: To appear in Physics of Fluids

Physics-informed laboratory estimation of Sargassum windage

Authors: M. J. Olascoaga, F. J. Beron-Vera, R. T. Beyea, G. Bonner, M. Castellucci, G. J. Goni, C. Guigand, N. F. Putman

Abstract: A recent Maxey--Riley theory for \sarg raft motion, which models a raft as a network of elastically interacting finite-size, buoyant particles, predicts the carrying flow velocity to be given by the weighted sum of the water and air velocities $(1-α)\mathbf{v} + α\mathbf w$. The theory provides a closed formula for parameter $α$, referred to as \emph{windage}, depending on water-to-particle-densit… ▽ More A recent Maxey--Riley theory for \sarg raft motion, which models a raft as a network of elastically interacting finite-size, buoyant particles, predicts the carrying flow velocity to be given by the weighted sum of the water and air velocities $(1-α)\mathbf{v} + α\mathbf w$. The theory provides a closed formula for parameter $α$, referred to as \emph{windage}, depending on water-to-particle-density ratio or buoyancy ($δ$). From a series of laboratory experiments in an air--water stream flume facility under controlled conditions, we estimate $α$ ranging from 0.02 to 0.96\pct. On average, our windage estimates can be up to 9 times smaller than considered in conventional \emph{Sargassum} raft transport modeling, wherein it is customary to add a fraction of $\mathbf w$ to $\mathbf{v}$ chosen in an ad-hoc piecemeal manner. Using the formula provided by the Maxey--Riley theory, we estimate $δ$ ranging from 1.00 to 1.49. This is consistent with direct $δ$ measurements, ranging from 0.9 to 1.25, which provide support for our $α$ estimation. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Journal ref: Physics of Fluids 35, 111702 (2023)

Coherent spore dispersion via drop-leaf interactions

Authors: Zixuan Wu, Saikat Basu, Seungho Kim, Mark Sorrells, Francisco J. Beron-Vera, Sunghwan Jung

Abstract: The dispersion of plant pathogens, such as rust spores, is responsible for more than 20% of global yield loss annually, and poses a significant threat to human health. However, the release mechanics of pathogens from flexible plant surfaces into the canopy is not well understood. In this study, we investigated the interplay between leaf elasticity and raindrop momentum, revealing how it induces fl… ▽ More The dispersion of plant pathogens, such as rust spores, is responsible for more than 20% of global yield loss annually, and poses a significant threat to human health. However, the release mechanics of pathogens from flexible plant surfaces into the canopy is not well understood. In this study, we investigated the interplay between leaf elasticity and raindrop momentum, revealing how it induces flow coherence and enhances spore transport with 2-10 times greater energy compared to impacts on stationary surfaces. We observed that a flexible leaf generates vortex dipoles, leading to a super-diffusive stream flow. We then developed a theoretical model that accurately predicted the average air flux from leaf edges and the vortex strength to be proportional the vibration speed of the leaves. With Lagrangian diagnostics, we further revealed the presence of hyperbolic and elliptical coherent structures around fluttering leaves, providing the dynamical description of spore transport. Our model demonstrated that a leaf aspect ratio (length/width) negatively correlates with dispersion, indicating that shorter and wider leaves promote greater pathogen spread. Additionally, we found that leaf rigidity positively correlates with dispersion due to damping effects. These mechanistic insights would help the construction of physically informed analytical models for improve local crop disease management. △ Less

Submitted 24 July, 2023; v1 submitted 24 June, 2023; originally announced June 2023.

Comments: 12 pages, 4 figures

Stability of temporal statistics in Transition Path Theory with sparse data

Authors: G. Bonner, F. J. Beron-Vera, M. J. Olascoaga

Abstract: Transition Path Theory (TPT) provides a rigorous statistical characterization of the ensemble of trajectories connecting directly, i.e., without detours, two disconnected (sets of) states in a Markov chain, a stochastic process that undergoes transitions from one state to another with probability depending on the state attained in the previous step. Markov chains can be constructed using trajector… ▽ More Transition Path Theory (TPT) provides a rigorous statistical characterization of the ensemble of trajectories connecting directly, i.e., without detours, two disconnected (sets of) states in a Markov chain, a stochastic process that undergoes transitions from one state to another with probability depending on the state attained in the previous step. Markov chains can be constructed using trajectory data via counting of transitions between cells covering the domain spanned by trajectories. With sparse trajectory data, the use of regular cells is observed to result in unstable estimates of the total duration of transition paths. Using Voronoi cells resulting from k-means clustering of the trajectory data, we obtain stable estimates of this TPT statistic, which is generalized to frame the remaining duration of transition paths, a new TPT statistic suitable for investigating connectivity. △ Less

Submitted 30 January, 2023; originally announced February 2023.

Comments: Submitted to Chaos. Comments welcomed!

Dynamical geography and transition paths of Sargassum in the tropical Atlantic

Authors: F. J. Beron-Vera, M. J. Olascoaga, N. Putman, J. Trinanes, G. J. Goni, R. Lumpkin

Abstract: By analyzing a time-homogeneous Markov chain constructed using trajectories of undrogued drifting buoys from the NOAA's Global Drifter Program, we find that probability density can distribute in a manner that resembles very closely the recently observed recurrent belt of high Sargassum density in the tropical Atlantic between 5--10$^{\circ}$N, coined the \emph{Great Atlantic Sargassum Belt} (\emph… ▽ More By analyzing a time-homogeneous Markov chain constructed using trajectories of undrogued drifting buoys from the NOAA's Global Drifter Program, we find that probability density can distribute in a manner that resembles very closely the recently observed recurrent belt of high Sargassum density in the tropical Atlantic between 5--10$^{\circ}$N, coined the \emph{Great Atlantic Sargassum Belt} (\emph{GASB}). A spectral analysis of the associated transition matrix further unveils a forward attracting almost-invariant set in the northwestern Gulf of Mexico with a corresponding basin of attraction disconnected from the Sargasso Sea, but including the nutrient-rich regions around the Amazon and Orinoco rivers mouths and also the upwelling system off the northern coast of west Africa. This represents a data-based inference of potential remote sources of Sargassum recurrently invading the Intra-Americas Seas (IAS). By further applying Transition Path Theory (TPT) on the data-derived Markov chain model, two potential pathways for Sargassum into the IAS from the upwelling system off the coast of Africa are revealed. One TPT-inferred pathway takes place along the GASB. The second pathway is more southern and slower, first going through the Gulf of Guinea, then across the tropical Atlantic toward the mouth of the Amazon River, and finally along the northeastern South American margin. The existence of such a southern TPT-inferred pathway may have consequences for bloom stimulation by nutrients from river runoff. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: To appear in AIP Advances

Exploring the use of Transition Path Theory in building an oil spill prediction scheme

Authors: M. J. Olascoaga, F. J. Beron-Vera

Abstract: The Transition Path Theory (TPT) of complex systems has proven a robust means for statistically characterizing the ensemble of trajectories that connect any two preset flow regions, say $\mathcal A$ and $\mathcal B$, directly. More specifically, transition paths are such that they start in $\mathcal A$ and then go to $\mathcal B$ without detouring back to $\mathcal A$ or $\mathcal B$. This way, th… ▽ More The Transition Path Theory (TPT) of complex systems has proven a robust means for statistically characterizing the ensemble of trajectories that connect any two preset flow regions, say $\mathcal A$ and $\mathcal B$, directly. More specifically, transition paths are such that they start in $\mathcal A$ and then go to $\mathcal B$ without detouring back to $\mathcal A$ or $\mathcal B$. This way, they make an effective contribution to the transport from $\mathcal A$ to $\mathcal B$. Here, we explore its use for building a scheme that enables predicting the evolution of an oil spill in the ocean. This involves appropriately adapting TPT such that it includes a reservoir that pumps oil into a typically open domain. Additionally, we lift up the restriction of the oil not to return to the spill site en route to a region that there is interest to be protected. TPT is applied on oil trajectories available up to the present, e.g., as integrated using velocities produced by a data assimilative system or as inferred from high-frequency radars, to make a prediction of transition oil paths beyond, without relying on forecasted oil trajectories. As a proof of concept we consider a hypothetical oil spill in the Trion oil field, under development within the Perdido Foldbelt in the northwestern Gulf of Mexico, and the \emph{Deepwater Horizon} oil spill. This is done using trajectories integrated from climatological and hindcast surface velocity and winds as well as produced by satellite-tracked surface drifting buoys, in each case discretized into a Markov chain that provides a framework for the TPT-based prediction. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: Submitted to Frontiers

Sampling--Dependent Transition Paths of Iceland-Scotland Overflow Water

Authors: F. J. Beron-Vera, M. J. Olascoaga, L. Helfmann, P. Miron

Abstract: In this note, we apply Transition Path Theory (TPT) from Markov chains to shed light on the problem of Iceland-Scotland Overflow Water (ISOW) equatorward export. A recent analysis of observed trajectories of submerged floats demanded revision of the traditional abyssal circulation theory, which postulates that ISOW should steadily flow along a deep boundary current (DBC) around the subpolar North… ▽ More In this note, we apply Transition Path Theory (TPT) from Markov chains to shed light on the problem of Iceland-Scotland Overflow Water (ISOW) equatorward export. A recent analysis of observed trajectories of submerged floats demanded revision of the traditional abyssal circulation theory, which postulates that ISOW should steadily flow along a deep boundary current (DBC) around the subpolar North Atlantic prior to exiting it. The TPT analyses carried out here allow to focus the attention on the portions of flow from the origin of ISOW to the region where ISOW exits the subpolar North Atlantic and suggest that insufficient sampling may be biasing the aforementioned demand. The analyses, appropriately adapted to represent a continuous input of ISOW, are carried out on three time-homogeneous Markov chains modeling the ISOW flow. One is constructed using a high number of simulated trajectories homogeneously covering the flow domain. The other two use much fewer trajectories which heterogeneously cover the domain. The trajectories in the latter two chains are observed trajectories or simulated trajectories subsampled at the observed frequency. While the densely sampled chain supports a well-defined DBC, the more heterogeneously sampled chains do not, irrespective of whether observed or simulated trajectories are used. Studying the sampling sensitivity of the Markov chains, we can give recommendations for enlarging the existing float dataset to improve the significance of conclusions about time-asymptotic aspects of the ISOW circulation. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: Submitted to JPO. Comments welcomed!

Do eddies connect the tropical Atlantic Ocean and the Gulf of Mexico?

Authors: F. Andrade-Canto, F. J. Beron-Vera

Abstract: Consistent with satellite-tracked trajectories of drogued drifters, but at odds with Eulerian assessment of satellite-altimetry measurements of sea-surface height, we show that North Brazil Currents Rings (NBCRs) are incapable of bypassing the Lesser Antilles as structures that coherently transport material. The nature of the inability of the de-facto oceanographic Eulerian, streamline-based eddy… ▽ More Consistent with satellite-tracked trajectories of drogued drifters, but at odds with Eulerian assessment of satellite-altimetry measurements of sea-surface height, we show that North Brazil Currents Rings (NBCRs) are incapable of bypassing the Lesser Antilles as structures that coherently transport material. The nature of the inability of the de-facto oceanographic Eulerian, streamline-based eddy detection technique to produce a correct assessment is rooted in its lack of objectivity. We arrive at this conclusion by applying \emph{geodesic eddy detection} on the altimetric dataset over nearly its entire extent. While we detect northwestward translating NBCRs that can be classified as coherent Lagrangian eddies, they typically experience strong filamentation and complete loss of coherence prior to reaching the Lesser Antilles. Moreover, the filamented material hardly penetrates into the Caribbean Sea, let alone the Gulf of Mexico, and not without substantively mixing with the ambient fluid east of the archipelago. △ Less

Submitted 22 September, 2022; v1 submitted 10 May, 2022; originally announced May 2022.

Comments: Submitted to Geophysical Research Letters. Typos and references fixed. Appendices added showing robustness of geodesic boundaries under flow perturbations; support for the x-to-o behavior in the bottom panels of Fig.~1; and proof of the lack of objectivity of PV

Carriers of \emph{Sargassum} and mechanism for coastal inundation in the Caribbean Sea

Authors: F. Andrade-Canto, F. J. Beron-Vera, G. J. Goni, D. Karrasch, M. J. Olascoaga, J. Trinanes

Abstract: We identify effective carriers of \emph{Sargassum} in the Caribbean Sea and describe a mechanism for coastal choking. Revealed from satellite altimetry, the carriers of \emph{Sargassum} are mesoscale eddies (vortices of 50-km radius or larger) with coherent material (i.e., fluid) boundaries. These are observer-independent -- \emph{unlike} eddy boundaries identified with instantaneously closed stre… ▽ More We identify effective carriers of \emph{Sargassum} in the Caribbean Sea and describe a mechanism for coastal choking. Revealed from satellite altimetry, the carriers of \emph{Sargassum} are mesoscale eddies (vortices of 50-km radius or larger) with coherent material (i.e., fluid) boundaries. These are observer-independent -- \emph{unlike} eddy boundaries identified with instantaneously closed streamlines of the altimetric sea-surface height field -- and furthermore harbor finite-time attractors for networks of elastically connected finite-size buoyant or "inertial" particles dragged by ocean currents and winds, a mathematical abstraction of \emph{Sargassum} rafts. The mechanism of coastal inundation, identified using a minimal model of surface-intensified Caribbean Sea eddies, is thermal instability in the presence of bottom topography. △ Less

Submitted 31 December, 2021; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: To appear in Physics Fluids

Transition paths of North Atlantic Deep Water

Authors: P. Miron, F. J. Beron-Vera, M. J. Olascoaga

Abstract: We use Transition Path Theory (TPT) to frame, in a statistically more robust fashion than earlier analyses, equatorward routes of North Atlantic Deep Water (NADW) in the subpolar North Atlantic. TPT is applied on all available RAFOS and Argo floats in the area by means of a discretization of the Lagrangian dynamics described by their trajectories. By considering floats at different depths, we inve… ▽ More We use Transition Path Theory (TPT) to frame, in a statistically more robust fashion than earlier analyses, equatorward routes of North Atlantic Deep Water (NADW) in the subpolar North Atlantic. TPT is applied on all available RAFOS and Argo floats in the area by means of a discretization of the Lagrangian dynamics described by their trajectories. By considering floats at different depths, we investigate transition paths of NADW in its upper (UNADW) and lower (LNADW) layers. We find that the majority of UNADW transition paths sourced in the Labrador and southwestern Irminger Seas reach the western side of a target arranged zonally along the southern edge of the subpolar North Atlantic domain visited by the floats. This is accomplished in the form of a well-organized deep boundary current (DBC). LNADW transition paths sourced west of the Reykjanes Ridge reveal a similar pattern, while those sourced east of the ridge are found to hit the western side of the target via a DBC and also several other places along it in a less organized fashion, indicating southward flow along the eastern and western flanks of the Mid-Atlantic Ridge. Naked-eye inspection of trajectories suggest generally much more diffusive equatorward NADW routes. A source-independent dynamical decomposition of the flow domain into analogous backward-time basins of attraction, beyond the reach of direct inspection of trajectories, reveals a much wider influence of the western side of the target for UNADW than for LNADW. For UNADW, the average expected duration of the pathways from the Labrador and Irminger Seas was found to be of 2 to 3 years. For LNADW, the duration was found to be influenced by the Reykjanes Ridge, being as long as 8 years from the western side of the ridge and of about 3 years on average from its eastern side. △ Less

Submitted 24 May, 2022; v1 submitted 31 August, 2021; originally announced August 2021.

Journal ref: Miron, P., F.J. Beron and M.J. Olascoaga (2022). Transition paths of North Atlantic Deep Water. The Journal of Atmospheric and Oceanic Technology

Extended shallow-water theories with thermodynamics and geometry

Authors: F. J. Beron-Vera

Abstract: Driven by growing momentum in two-dimensional geophysical flow modeling, this paper introduces a general family of "thermal" rotating shallow-water models. The models are capable of accommodating thermodynamic processes, such as those acting in the ocean mixed layer, by allowing buoyancy to vary in horizontal position and time as well as with depth, in a polynomial fashion up to an arbitrary degre… ▽ More Driven by growing momentum in two-dimensional geophysical flow modeling, this paper introduces a general family of "thermal" rotating shallow-water models. The models are capable of accommodating thermodynamic processes, such as those acting in the ocean mixed layer, by allowing buoyancy to vary in horizontal position and time as well as with depth, in a polynomial fashion up to an arbitrary degree. Moreover, the models admit Euler--Poincare variational formulations and possess Lie--Poisson Hamiltonian structure. Such a geometric property provides solid fundamental support to the theories described with consequences for numerical implementation and the construction of unresolved motion parametrizations. In particular, it is found that stratification halts the development of small-scale filament rollups recently observed in a popular model, which, having vertically homogeneous density, represents a special case of the models presented here. △ Less

Submitted 22 September, 2021; v1 submitted 20 August, 2021; originally announced August 2021.

Comments: Revised; to appear in Physics of Fluids

Geometry of shallow-water dynamics with thermodynamics

Authors: F. J. Beron-Vera

Abstract: We review the geometric structure of the IL$^0$PE model, a rotating shallow-water model with variable buoyancy, thus sometimes called ``thermal'' shallow-water model. We start by discussing the Euler--Poincaré equations for rigid body dynamics and the generalized Hamiltonian structure of the system. We then reveal similar geometric structure for the IL$^0$PE. We show, in particular, that the model… ▽ More We review the geometric structure of the IL$^0$PE model, a rotating shallow-water model with variable buoyancy, thus sometimes called ``thermal'' shallow-water model. We start by discussing the Euler--Poincaré equations for rigid body dynamics and the generalized Hamiltonian structure of the system. We then reveal similar geometric structure for the IL$^0$PE. We show, in particular, that the model equations and its (Lie--Poisson) Hamiltonian structure can be deduced from Morrison and Greene's (1980) system upon ignoring the magnetic field ($\vec{\mathrm B} = 0$) and setting $U(ρ,s) = \frac{1}{2}ρs$, where $ρ$ is mass density and $s$ is entropy per unit mass. These variables play the role of layer thickness ($h$) and buoyancy ($\t$) in the IL$^0$PE, respectively. Included in an appendix is an explicit proof of the Jacobi identity satisfied by the Poisson bracket of the system. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: This is review work (with a novel result on an overlooked connection with earlier work) carried out during the Aspen Center for Physics workshop Transport and Mixing of Tracers in Geophysics and Astrophysics, which has motivated original ongoing work

Nonlinear saturation of thermal instabilities

Authors: F. J. Beron-Vera

Abstract: Low-frequency simulations of a one-layer model with lateral buoyancy variations (i.e., thermodynamically active) have revealed circulatory motions resembling quite closely submesoscale observations in the surface ocean rather than indefinitely growing in the absence of a high-wavenumber instability cutoff. In this note it is shown that the existence of a convex pseudoenergy--momentum integral of m… ▽ More Low-frequency simulations of a one-layer model with lateral buoyancy variations (i.e., thermodynamically active) have revealed circulatory motions resembling quite closely submesoscale observations in the surface ocean rather than indefinitely growing in the absence of a high-wavenumber instability cutoff. In this note it is shown that the existence of a convex pseudoenergy--momentum integral of motion for the inviscid, unforced dynamics provides a mechanism for the nonlinear saturation of such thermal instabilities in the zonally symmetric case. The result is an application of \citet{Arnold-66} and \citet{Shepherd-88a} methods. △ Less

Submitted 22 January, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: Comments welcomed!

Transition paths of marine debris and the stability of the garbage patches

Authors: P. Miron, F. J. Beron-Vera, L. Helfmann, P. Koltai

Abstract: We used transition path theory (TPT) to infer "reactive" pathways of floating marine debris trajectories. The TPT analysis was applied on a pollution-aware time-homogeneous Markov chain model constructed from trajectories produced by satellite-tracked undrogued buoys from the NOAA Global Drifter Program. The latter involved coping with the openness of the system in physical space, which further re… ▽ More We used transition path theory (TPT) to infer "reactive" pathways of floating marine debris trajectories. The TPT analysis was applied on a pollution-aware time-homogeneous Markov chain model constructed from trajectories produced by satellite-tracked undrogued buoys from the NOAA Global Drifter Program. The latter involved coping with the openness of the system in physical space, which further required an adaptation of the standard TPT setting. Directly connecting pollution sources along coastlines with garbage patches of varied strengths, the unveiled reactive pollution routes represent alternative targets for ocean cleanup efforts. Among our specific findings we highlight: constraining a highly probable pollution source for the Great Pacific Garbage Patch; characterizing the weakness of the Indian Ocean gyre as a trap for plastic waste; and unveiling a tendency of the subtropical gyres to export garbage toward the coastlines rather than to other gyres in the event of anomalously intense winds. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: Submitted to Chaos

Genesis, evolution, and apocalypse of Loop Current rings

Authors: F Andrade-Cano, D. Karrasch, F. J. Beron-Vera

Abstract: We carry out assessments of the life cycle of Loop Current vortices, so-called rings, in the Gulf of Mexico by applying three objective (i.e., observer-independent) coherent Lagrangian vortex detection methods on velocities derived from satellite altimetry measurements of sea-surface height (SSH). The methods reveal material vortices with boundaries that withstand stretching or diffusion, or whose… ▽ More We carry out assessments of the life cycle of Loop Current vortices, so-called rings, in the Gulf of Mexico by applying three objective (i.e., observer-independent) coherent Lagrangian vortex detection methods on velocities derived from satellite altimetry measurements of sea-surface height (SSH). The methods reveal material vortices with boundaries that withstand stretching or diffusion, or whose fluid elements rotate evenly. This involved a technology advance that enables framing vortex genesis and apocalypse robustly and with precision. We find that the stretching- and diffusion-withstanding assessments produce consistent results, which show large discrepancies with Eulerian assessments that identify vortices with regions instantaneously filled with streamlines of the SSH field. The even-rotation assessment, which is vorticity-based, is found to be quite unstable, suggesting life expectancies much shorter than those produced by all other assessments. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: Submitted to Physics of Fluids

Nonlinear dynamics of inertial particles in the ocean: From drifters and floats to marine debris and Sargassum

Authors: F. J. Beron-Vera

Abstract: Buoyant, finite-size or inertial particle motion is fundamentally unlike neutrally buoyant, infinitesimally small or Lagrangian particle motion. The de-jure fluid mechanics framework for the description of inertial particle dynamics is provided by the Maxey-Riley equation. Derived from first principles - a result of over a century of research since the pioneering work by Sir George Stokes - the Ma… ▽ More Buoyant, finite-size or inertial particle motion is fundamentally unlike neutrally buoyant, infinitesimally small or Lagrangian particle motion. The de-jure fluid mechanics framework for the description of inertial particle dynamics is provided by the Maxey-Riley equation. Derived from first principles - a result of over a century of research since the pioneering work by Sir George Stokes - the Maxey-Riley equation is a Newton-type-law with several forces including (mainly) flow, added mass, shear-induced lift, and drag forces. In this paper we present an overview of recent efforts to port the Maxey-Riley framework to oceanography. These involved: 1) including the Coriolis force, which was found to explain behavior of submerged floats near mesoscale eddies; 2) accounting for the combined effects of ocean current and wind drag on inertial particles floating at the air-sea interface, which helped understand the formation of great garbage patches and the role of anticyclonic eddies as plastic debris traps; and 3) incorporating elastic forces, which are needed to simulate the drift of pelagic Sargassum. Insight on the nonlinear dynamics of inertial particles in every case was possible to be achieved by investigating long-time asymptotic behavior in the various Maxey-Riley equation forms, which represent singular perturbation problems involving slow and fast variables. △ Less

Submitted 8 November, 2020; v1 submitted 21 July, 2020; originally announced July 2020.

Journal ref: Nonlinear Dynamics (2020)

Clustering of marine-debris- and \emph{Sargassum}-like drifters explained by inertial particle dynamics}

Authors: P. Miron, M. J. Olascoaga, F. J. Beron-Vera, N. F. Putman, J. Trinanes, R. Lumpkin, G. J. Goni

Abstract: Drifters designed to mimic floating marine debris and small patches of pelagic \emph{Sargassum} were satellite tracked in four regions across the North Atlantic. Though subjected to the same initial conditions at each site, the tracks of different drifters quickly diverged after deployment. We explain the clustering of drifter types using a recent Maxey-Riley theory for surface ocean inertial part… ▽ More Drifters designed to mimic floating marine debris and small patches of pelagic \emph{Sargassum} were satellite tracked in four regions across the North Atlantic. Though subjected to the same initial conditions at each site, the tracks of different drifters quickly diverged after deployment. We explain the clustering of drifter types using a recent Maxey-Riley theory for surface ocean inertial particle dynamics applied on multidata-based mesoscale ocean currents and winds from reanalysis. Simulated trajectories of objects at the air-sea interface are significantly improved when represented as inertial (accounting for buoyancy and size), rather than as perfectly Lagrangian (fluid following) particles. Separation distances between simulated and observed trajectories were substantially smaller for debris-like drifters than for \emph{Sargassum}-like drifters, suggesting that additional consideration of its physical properties relative to fluid velocities may be useful. Our findings can be applied to model variability in movements and distribution of diverse objects floating at the ocean surface. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Comments: Submitted to Geophy. Res. Lett

A minimal Maxey--Riley model for the drift of \emph{Sargassum} rafts

Authors: Francisco J Beron-Vera, Philippe Miron

Abstract: Inertial particles (i.e. with mass and of finite size) immersed in a fluid in motion are unable to adapt their velocities to the carrying flow and thus they have been the subject of much interest in fluid mechanics. In this paper we consider an ocean setting with inertial particles elastically connected forming a network that floats at the interface with the atmosphere. The network evolves accordi… ▽ More Inertial particles (i.e. with mass and of finite size) immersed in a fluid in motion are unable to adapt their velocities to the carrying flow and thus they have been the subject of much interest in fluid mechanics. In this paper we consider an ocean setting with inertial particles elastically connected forming a network that floats at the interface with the atmosphere. The network evolves according to a recently derived and validated Maxey--Riley equation for inertial particle motion in the ocean. We rigorously show that, under sufficiently calm wind conditions, rotationally coherent quasigeostrophic vortices (which have material boundaries that resist outward filamentation) always possess finite-time attractors for elastic networks if they are anticyclonic, while if they are cyclonic provided that the networks are sufficiently stiff. This result is supported numerically under more general wind conditions and, most importantly, is consistent with observations of rafts of pelagic \emph{Sargassum}, for which the elastic inertial networks represent a minimal model. Furthermore, our finding provides an effective mechanism for the long range transport of \emph{Sargassum}, and thus for its connectivity between accumulation regions and remote sources. △ Less

Submitted 31 July, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: To appear in J. Fluid Mechanics

Observation and quantification of inertial effects on the drift of floating objects at the ocean surface

Authors: M. J. Olascoaga, F. J. Beron-Vera, P. Miron, J. Triñanes, N. F. Putman, R. Lumpkin, G. J. Goni

Abstract: We present results from an experiment designed to better understand the mechanism by which ocean currents and winds control flotsam drift. The experiment consisted in deploying in the Florida Current and subsequently satellite tracking specially designed drifting buoys of varied sizes, buoyancies, and shapes. We explain the differences in the trajectories described by the special drifters as a res… ▽ More We present results from an experiment designed to better understand the mechanism by which ocean currents and winds control flotsam drift. The experiment consisted in deploying in the Florida Current and subsequently satellite tracking specially designed drifting buoys of varied sizes, buoyancies, and shapes. We explain the differences in the trajectories described by the special drifters as a result of their inertia, primarily buoyancy, which constrains the ability of the drifters to adapt their velocities to instantaneous changes in the ocean current and wind that define the carrying flow field. Our explanation of the observed behavior follows from the application of a recently proposed Maxey--Riley theory for the motion of finite-size particles floating at the surface ocean. The nature of the carrying flow and the domain of validity of the theory are clarified, and a closure proposal is made to fully determine its parameters in terms of the carrying fluid system properties and inertial particle characteristics. △ Less

Submitted 17 November, 2019; v1 submitted 10 November, 2019; originally announced November 2019.

Comments: Physics of Fluids, submitted 17-Nov-2019

Building a Maxey--Riley framework for surface ocean inertial particle dynamics

Authors: F. J. Beron-Vera, M. J. Olascoaga, P. Miron

Abstract: 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 air--sea interface with unsteady nonuniform winds and ocean currents ab… ▽ More 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 air--sea 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. △ Less

Submitted 10 August, 2019; v1 submitted 28 May, 2019; originally announced June 2019.

Comments: To appear in Phys. Fluids

Markov-chain-inspired search for MH370

Authors: P. Miron, F. J. Beron-Vera, M. J. Olascoaga, P. Koltai

Abstract: Markov-chain models are constructed for the probabilistic description of the drift of marine debris from Malaysian Airlines flight MH370. En route from Kuala Lumpur to Beijing, the MH370 mysteriously disappeared in the southeastern Indian Ocean on 8 March 2014, somewhere along the arc of the 7th ping ring around the Inmarsat-3F1 satellite position when the airplane lost contact. The models are obt… ▽ More Markov-chain models are constructed for the probabilistic description of the drift of marine debris from Malaysian Airlines flight MH370. En route from Kuala Lumpur to Beijing, the MH370 mysteriously disappeared in the southeastern Indian Ocean on 8 March 2014, somewhere along the arc of the 7th ping ring around the Inmarsat-3F1 satellite position when the airplane lost contact. The models are obtained by discretizing the motion of undrogued satellite-tracked surface drifting buoys from the global historical data bank. A spectral analysis, Bayesian estimation, and the computation of most probable paths between the Inmarsat arc and confirmed airplane debris beaching sites are shown to constrain the crash site, near 25$^{\circ}$S on the Inmarsat arc. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: Submitted to Chaos

Journal ref: Chaos 29, 041105 (2019)

Persistent Lagrangian transport patterns in the northwestern Gulf of Mexico

Authors: M. K. Gough, F. J. Beron-Vera, M. J. Olascoaga, J. Sheinbaum, J. Juoanno, R. Duran

Abstract: Persistent Lagrangian transport patterns at the ocean surface are revealed from Lagrangian Coherent Structures (LCSs) computed from daily climatological surface current velocities in the northwestern Gulf of Mexico (NWGoM). The velocities are produced by a submesoscale permitting regional ocean model of the Gulf of Mexico. The significance of the climatological LCSs (cLCSs) is supported with ensem… ▽ More Persistent Lagrangian transport patterns at the ocean surface are revealed from Lagrangian Coherent Structures (LCSs) computed from daily climatological surface current velocities in the northwestern Gulf of Mexico (NWGoM). The velocities are produced by a submesoscale permitting regional ocean model of the Gulf of Mexico. The significance of the climatological LCSs (cLCSs) is supported with ensemble-mean drifter density evolutions from simulated and historical satellite-tracked drifter trajectories. A persistent attracting barrier between the NWGoM shelf and the deep ocean is effectively identified by a hook-like pattern associated with groups of overall strongly attracting cLCSs that extend along the shelf break. Localized reductions in the attraction rate along these overall strongly attracting cLCSs proximal to cross-shore oriented cLCSs identify a pathway for potential transport across the shelf break. Groups of overall weakly-attracting cLCSs are not seen to strongly constrain material transport. Tracers originating over the shelf tend to be trapped there by the hook-like pattern as they spread cyclonically. Tracers originating beyond the shelf tend to be initially attracted to the hook-like pattern as they spread anti-cyclonically and eventually over the deep ocean. The findings have important implications for the mitigation of contaminant accidents such as oil spills. △ Less

Submitted 11 October, 2017; originally announced October 2017.

Comments: Submitted to J. Phys. Oceanography

Enduring Lagrangian coherence of a Loop Current ring assessed using independent observations

Authors: M. J. Olascoaga, F. J. Beron-Vera, Y. Wang, J. Triñanes, P. Pérez-Brunius

Abstract: Ocean flows are routinely inferred from low-resolution satellite altimetry measurements of sea surface height assuming a geostrophic balance. Recent nonlinear dynamical systems techniques have revealed that surface currents derived from altimetry can support mesoscale eddies with material boundaries that do not filament for many months, thereby representing effective transport mechanisms. However,… ▽ More Ocean flows are routinely inferred from low-resolution satellite altimetry measurements of sea surface height assuming a geostrophic balance. Recent nonlinear dynamical systems techniques have revealed that surface currents derived from altimetry can support mesoscale eddies with material boundaries that do not filament for many months, thereby representing effective transport mechanisms. However, the long-range Lagrangian coherence assessed for mesoscale eddy boundaries detected from altimetry is constrained by the impossibility of current altimeters to resolve ageostrophic submesoscale motions. These may act to prevent Lagrangian coherence from manifesting in the rigorous form described by the nonlinear dynamical systems theories. Here we use a combination of satellite ocean color and surface drifter trajectory data, rarely available simultaneously over an extended period of time, to provide observational evidence for the enduring Lagrangian coherence of a Loop Current ring detected from altimetry. We also seek indications of this behavior in the flow produced by a data-assimilative system which demonstrated ability to reproduce observed relative dispersion statistics down into the marginally submesoscale range. However, the simulated flow, total surface and subsurface or subsampled emulating altimetry, is not found to support the long-lasting Lagrangian coherence that characterizes the observed ring. This highlights the importance of the Lagrangian metrics produced by the nonlinear dynamical systems tools employed here in assessing model performance. △ Less

Submitted 13 July, 2018; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: In press in nature.com/Scientific Reports

Extracting quasi-steady Lagrangian transport patterns from the ocean circulation: An application to the Gulf of Mexico

Authors: Rodrigo Duran, Francisco J. Beron-Vera, Maria Josefina Olascoaga

Abstract: We construct a climatology of Lagrangian coherent structures (LCSs), the concealed skeleton that shapes transport, with a twelve-year-long data-assimilative simulation of the sea-surface circulation in the Gulf of Mexico (GoM). Computed as time-mean Cauchy-Green strain tensorlines of the climatological velocity, the climatological LCSs (cLCSs) unveil recurrent Lagrangian circulation patterns. cLCS… ▽ More We construct a climatology of Lagrangian coherent structures (LCSs), the concealed skeleton that shapes transport, with a twelve-year-long data-assimilative simulation of the sea-surface circulation in the Gulf of Mexico (GoM). Computed as time-mean Cauchy-Green strain tensorlines of the climatological velocity, the climatological LCSs (cLCSs) unveil recurrent Lagrangian circulation patterns. cLCSs strongly constrain the ensemble-mean Lagrangian circulation of the instantaneous model velocity, thus we show that a climatological velocity may preserve meaningful transport information. Also, the climatological transport patterns we report agree well with GoM kinematics and dynamics, as described in several previous observational and numerical studies. For example, cLCSs identify regions of persistent isolation, and suggest that coastal regions previously identified as high-risk for pollution impact, are regions of maximal attraction. Also, we show examples where cLCSs are remarkably similar to transport patterns observed during the Deepwater Horizon and Ixtoc oil spills, and during the Grand LAgrangian Deployment (GLAD) experiment. Thus, it is shown that cLCSs are an efficient way of synthesizing vast amounts of Lagrangian information. The cLCS method confirms previous GoM studies, and contributes to our understanding by revealing the persistent nature of the dynamics and kinematics treated therein. △ Less

Submitted 21 October, 2017; v1 submitted 7 April, 2017; originally announced April 2017.

Comments: To be submitted

Lagrangian dynamical geography of the Gulf of Mexico

Authors: P. Miron, F. J. Beron-Vera, M. J. Olascoaga, J. Sheinbaum, P. Perez-Brunius, G. Froyland

Abstract: We construct a Markov-chain representation of the surface-ocean Lagrangian dynamics in a region occupied by the Gulf of Mexico (GoM) and adjacent portions of the Caribbean Sea and North Atlantic using satellite-tracked drifter trajectory data, the largest collection so far considered. From the analysis of the eigenvectors of the transition matrix associated with the chain, we identify almost-invar… ▽ More We construct a Markov-chain representation of the surface-ocean Lagrangian dynamics in a region occupied by the Gulf of Mexico (GoM) and adjacent portions of the Caribbean Sea and North Atlantic using satellite-tracked drifter trajectory data, the largest collection so far considered. From the analysis of the eigenvectors of the transition matrix associated with the chain, we identify almost-invariant attracting sets and their basins of attraction. With this information we decompose the GoM's geography into weakly dynamically interacting provinces, which constrain the connectivity between distant locations within the GoM. Offshore oil exploration, oil spill contingency planning, and fish larval connectivity assessment are among the many activities that can benefit from the dynamical information carried in the geography constructed here. △ Less

Submitted 6 June, 2017; v1 submitted 30 March, 2017; originally announced March 2017.

Comments: Submitted to Scientific Reports

Inertia-induced accumulation of flotsam in the subtropical gyres

Authors: Francisco J. Beron-Vera, Maria J. Olascoaga, Rick Lumpkin

Abstract: Recent surveys of marine plastic debris density have revealed high levels in the center of the subtropical gyres. Earlier studies have argued that the formation of great garbage patches is due to Ekman convergence in such regions. In this work we report a tendency so far overlooked of drogued and undrogued drifters to accumulate distinctly over the subtropical gyres, with undrogued drifters accumu… ▽ More Recent surveys of marine plastic debris density have revealed high levels in the center of the subtropical gyres. Earlier studies have argued that the formation of great garbage patches is due to Ekman convergence in such regions. In this work we report a tendency so far overlooked of drogued and undrogued drifters to accumulate distinctly over the subtropical gyres, with undrogued drifters accumulating in the same areas where plastic debris accumulate. We show that the observed accumulation is too fast for Ekman convergence to explain it. We demonstrate that the accumulation is controlled by finite-size and buoyancy (i.e., inertial) effects on undrogued drifter motion subjected to ocean current and wind drags. We infer that the motion of flotsam in general is constrained by similar effects. This is done by using a newly proposed Maxey--Riley equation which models the submerged (surfaced) drifter portion as a sphere of the fractional volume that is submerged (surfaced). △ Less

Submitted 4 October, 2016; v1 submitted 8 June, 2016; originally announced June 2016.

Comments: Submitted to Geophys. Res. Letters

The life cycle of a coherent Lagrangian Agulhas ring

Authors: Y. Wang, F. J. Beron-Vera, M. J. Olascoaga

Abstract: We document the long-term evolution of an Agulhas ring detected from satellite altimetry using a technique from nonlinear dynamical systems that enables objective (i.e., observer-independent) eddy framing. Such objectively detected eddies have Lagrangian (material) boundaries that remain coherent (unfilamented) over the detection period. The ring preserves a quite compact material entity for a per… ▽ More We document the long-term evolution of an Agulhas ring detected from satellite altimetry using a technique from nonlinear dynamical systems that enables objective (i.e., observer-independent) eddy framing. Such objectively detected eddies have Lagrangian (material) boundaries that remain coherent (unfilamented) over the detection period. The ring preserves a quite compact material entity for a period of about 2 years even after most initial coherence is lost within 5 months after detection. We attribute this to the successive development of short-term coherent material boundaries around the ring. These boundaries provide effective short-term shielding for the ring, which prevents a large fraction of the ring's interior from being mixed with the ambient turbulent flow. We show that such coherence regain events cannot be inferred from Eulerian analysis. This process is terminated by a ring-splitting event which marks the ring demise, near the South American coast. The genesis of the ring is characterized by a ring-merging event away from the Agulhas retroflection, followed by a 4-month-long partial coherence stage, scenario that is quite different than a simple current occlusion and subsequent eddy pinch off. △ Less

Submitted 7 January, 2016; originally announced January 2016.

Comments: Submitted to J. Geophys. Res. Oceans

Statistics of simulated and observed pair separations in the Gulf of Mexico

Authors: F. J. Beron-Vera, J. H. LaCasce

Abstract: Pair-separation statistics of in-situ and synthetic surface drifters deployed near the \emph{Deepwater Horizon} site in the Gulf of Mexico are investigated. The synthetic trajectories derive from a 1-km-resolution data-assimilative Navy Coastal Ocean Model (NCOM) simulation. The in-situ drifters were launched in the Grand LAgrangian Deployment (GLAD). Diverse measures of the dispersion are calcula… ▽ More Pair-separation statistics of in-situ and synthetic surface drifters deployed near the \emph{Deepwater Horizon} site in the Gulf of Mexico are investigated. The synthetic trajectories derive from a 1-km-resolution data-assimilative Navy Coastal Ocean Model (NCOM) simulation. The in-situ drifters were launched in the Grand LAgrangian Deployment (GLAD). Diverse measures of the dispersion are calculated and compared to theoretical predictions. For the NCOM pairs, the measures indicate nonlocal pair dispersion at the smallest sampled scales. At separations exceeding 100 km, pair motion is uncorrelated, indicating absolute rather than relative dispersion. With the GLAD drifters however the statistics suggest local dispersion (in which pair separations exhibit power law growth), in line with previous findings. The disagreement stems in part from inertial oscillations, which affect the energy levels at small scales without greatly altering the net particle displacements. They were significant in GLAD but much weaker in the NCOM simulations. In addition the GLAD drifters were launched close together, producing few independent realizations and hence weaker statistical significance. Restricting the NCOM set to those launched at the same locations yields very similar statistics. △ Less

Submitted 19 April, 2016; v1 submitted 13 May, 2015; originally announced May 2015.

Comments: In press, Journal of Physical Oceanography, 2016

Coherent water transport across the South Atlantic

Authors: Y. Wang, M. J. Olascoaga, F. J. Beron-Vera

Abstract: The role of mesoscale eddies in transporting Agulhas leakage is investigated using a recent technique from nonlinear dynamical systems theory applied on geostrophic currents inferred from the over two-decade-long satellite altimetry record. Eddies are found to acquire material coherence away from the Agulhas retroflection, near the Walvis Ridge in the South Atlantic. Yearly, 1 to 4 coherent materi… ▽ More The role of mesoscale eddies in transporting Agulhas leakage is investigated using a recent technique from nonlinear dynamical systems theory applied on geostrophic currents inferred from the over two-decade-long satellite altimetry record. Eddies are found to acquire material coherence away from the Agulhas retroflection, near the Walvis Ridge in the South Atlantic. Yearly, 1 to 4 coherent material eddies are detected with diameters ranging from 40 to 280 km. A total of 23 eddy cores of about 50 km in diameter and with at least 30\pct of their contents traceable into the Indian Ocean were found to travel across the subtropical gyre with minor filamentation. Only 1 eddy core was found to pour its contents on the North Brazil Current. While ability of eddies to carry Agulhas leakage northwestward across the South Atlantic is supported by our analysis, this is more restricted than suggested by earlier ring transport assessments. △ Less

Submitted 1 May, 2015; v1 submitted 18 April, 2015; originally announced April 2015.

Comments: In press

Journal ref: Geophysical Research Letters, 2015

Dissipative inertial transport patterns near coherent Lagrangian eddies in the ocean

Authors: F. J. Beron-Vera, M. J. Olascoaga, G. Haller, M. Farazmand, J. Trinanes, Y. Wang

Abstract: Recent developments in dynamical systems theory have revealed long-lived and coherent Lagrangian (i.e., material) eddies in incompressible, satellite-derived surface ocean velocity fields. Paradoxically, observed drifting buoys and floating matter tend to create dissipative-looking patterns near oceanic eddies, which appear to be inconsistent with the conservative fluid particle patterns created b… ▽ More Recent developments in dynamical systems theory have revealed long-lived and coherent Lagrangian (i.e., material) eddies in incompressible, satellite-derived surface ocean velocity fields. Paradoxically, observed drifting buoys and floating matter tend to create dissipative-looking patterns near oceanic eddies, which appear to be inconsistent with the conservative fluid particle patterns created by coherent Lagrangian eddies. Here we show that inclusion of inertial effects (i.e., those produced by the buoyancy and size finiteness of an object) in a rotating two-dimensional incompressible flow context resolves this paradox. Specifically, we obtain that anticyclonic coherent Lagrangian eddies attract (repel) negatively (positively) buoyant finite-size particles, while cyclonic coherent Lagrangian eddies attract (repel) positively (negatively) buoyant finite-size particles. We show how these results explain dissipative-looking satellite-tracked surface drifter and subsurface float trajectories, as well as satellite-derived \emph{Sargassum} distributions. △ Less

Submitted 23 February, 2015; v1 submitted 27 August, 2014; originally announced August 2014.

Comments: Submitted to \emph{Chaos} Focus Issue on Objective detection of Lagrangian Coherent Structures. Revised 23-Feb-15

Submesoscale dispersion in the vicinity of the Deepwater Horizon spill

Authors: Andrew C. Poje, Tamay M. Özgökmen, Bruce Lipphardt, Jr., Brian K. Haus, Edward H. Ryan, Angelique C. Haza, A. J. H. M. Reniers, Josefina Olascoaga, Guillaume Novelli, Francisco J. Beron-Vera, Shuyi Chen, Arthur J. Mariano, Gregg Jacobs, Pat Hogan, Emanuel Coelho, A. D. Kirwan, Jr., Helga Huntley, Annalisa Griffa

Abstract: 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… ▽ More 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 meters to tens of kilometers, hours to days) remain unresolved due to a lack of synoptic measurements at these scales. \textcolor{black} {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 Mexico. Observed two-point statistics confirm the accuracy of classic turbulence scaling laws at 200m$-$50km scales and clearly indicate that dispersion at the submesoscales is \textit{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. △ Less

Submitted 18 July, 2014; v1 submitted 11 July, 2014; originally announced July 2014.

Comments: 9 pages, 6 figures

MSC Class: 76-XX

Journal ref: Proceedings of the National Academy of Sciences, vol. 111, no. 35, pp. 12693-12698, 2014

Addendum to "Coherent Lagrangian vortices: The black holes of turbulence"

Authors: G. Haller, F. J. Beron-Vera

Abstract: In Haller and Beron-Vera (2013) we developed a variational principle for the detection of coherent Lagrangian vortex boundaries. The solutions of this variational principle turn out to be closed null-geodesics of the Lorentzian metric associated with a generalized Green-Lagrange strain tensor family. This metric interpretation implies a mathematical analogy between coherent Lagrangian vortex bound… ▽ More In Haller and Beron-Vera (2013) we developed a variational principle for the detection of coherent Lagrangian vortex boundaries. The solutions of this variational principle turn out to be closed null-geodesics of the Lorentzian metric associated with a generalized Green-Lagrange strain tensor family. This metric interpretation implies a mathematical analogy between coherent Lagrangian vortex boundaries and photon spheres in general relativity. Here we give an improved discussion on this analogy. △ Less

Submitted 8 July, 2014; v1 submitted 4 June, 2014; originally announced June 2014.

Comments: Revised 27 June 2014

Coherent Lagrangian vortices: The black holes of turbulence

Authors: G. Haller, F. J. Beron-Vera

Abstract: 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 l… ▽ More 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. △ Less

Submitted 10 August, 2013; originally announced August 2013.

Comments: To appear in JFM Rapids

Zonal Jets as Transport Barriers in Planetary Atmospheres

Authors: F. J. Beron-Vera, M. G. Brown, M. J. Olascoaga, I. I. Rypina, H. Kocak, I. A. Udovydchenkov

Abstract: The connection between transport barriers and potential vorticity (PV) barriers in PV-conserving flows is investigated with a focus on zonal jets in planetary atmospheres. A perturbed PV-staircase model is used to illustrate important concepts. This flow consists of a sequence of narrow eastward and broad westward zonal jets with a staircase PV structure; the PV-steps are at the latitudes of the… ▽ More The connection between transport barriers and potential vorticity (PV) barriers in PV-conserving flows is investigated with a focus on zonal jets in planetary atmospheres. A perturbed PV-staircase model is used to illustrate important concepts. This flow consists of a sequence of narrow eastward and broad westward zonal jets with a staircase PV structure; the PV-steps are at the latitudes of the cores of the eastward jets. Numerically simulated solutions to the quasigeostrophic PV conservation equation in a perturbed PV-staircase flow are presented. These simulations reveal that both eastward and westward zonal jets serve as robust meridional transport barriers. The surprise is that westward jets, across which the background PV gradient vanishes, serve as robust transport barriers. A theoretical explanation of the underlying barrier mechanism is provided. It is argued that transport barriers near the cores of westward zonal jets, across which the background PV gradient is small, are found in Jupiter's midlatitude weather layer and in the Earth's summer hemisphere subtropical stratosphere. △ Less

Submitted 19 March, 2008; originally announced March 2008.

Comments: Accepted for publication in JAS

Tracing the early development of harmful algal blooms with the aid of Lagrangian coherent structures

Authors: M. J. Olascoaga, F. J. Beron-Vera, L. E. Brand, H. Koçak

Abstract: Several theories have been proposed to explain the development of harmful algal blooms (HABs) produced by the toxic dinoflagellate \emph{Karenia brevis} on the West Florida Shelf. However, because the early stages of HAB development are usually not detected, these theories have been so far very difficult to verify. In this paper we employ simulated \emph{Lagrangian coherent structures} (LCSs) to… ▽ More Several theories have been proposed to explain the development of harmful algal blooms (HABs) produced by the toxic dinoflagellate \emph{Karenia brevis} on the West Florida Shelf. However, because the early stages of HAB development are usually not detected, these theories have been so far very difficult to verify. In this paper we employ simulated \emph{Lagrangian coherent structures} (LCSs) to trace the early location of a HAB in late 2004 before it was transported to an area where it could be detected by satellite imagery, and then we make use of a population dynamics model to infer the factors that may have led to its development. The LCSs, which are computed based on a surface flow description provided by an ocean circulation model, delineate past and future histories of boundaries of passively advected fluid domains. The population dynamics model determines nitrogen in two components, nutrients and phytoplankton, which are assumed to be passively advected by the simulated surface currents. Two nearshore nutrient sources are identified for the HAB whose evolution is found to be strongly tied to the simulated LCSs. While one nutrient source can be associated with a coastal upwelling event, the other is seen to be produced by river runoff, which provides support to a theory of HAB development that considers nutrient loading into coastal waters produced by human activities as a critical element. Our results show that the use of simulated LCSs and a population dynamics model can greatly enhance our understanding of the early stages of the development of HABs. △ Less

Submitted 30 August, 2007; originally announced August 2007.

Comments: Submitted to JGR-Oceans

Persistent Transport Barrier on the West Florida Shelf

Authors: M. J. Olascoaga, I. I. Rypina, M. G. Brown, F. J. Beron-Vera, H. Koçak, L. E. Brand, G. R. Halliwell, L. K. Shay

Abstract: 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 so-called ``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 reco… ▽ More 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 so-called ``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 identify Lagrangian coherent structures (LCSs), which reveal the presence of a robust and 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 $\times$ 80 km region on the WFS off Tampa Bay is also used to identify LCSs, which reveal the presence of robust transient transport barriers. While the HF-radar-derived 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. △ Less

Submitted 7 August, 2006; originally announced August 2006.

Comments: Submitted to Geophysical Research Letters

On the Lagrangian Dynamics of Atmospheric Zonal Jets and the Permeability of the Stratospheric Polar Vortex

Authors: I. I. Rypina, F. J. Beron-Vera, M. G. Brown, H. Kocak, M. J. Olascoaga, I. A. Udovydchenkov

Abstract: The Lagrangian dynamics of zonal jets in the atmosphere are considered, with particular attention paid to explaining why, under commonly encountered conditions, zonal jets serve as barriers to meridional transport. The velocity field is assumed to be two-dimensional and incompressible, and composed of a steady zonal flow with an isolated maximum (a zonal jet) on which two or more travelling Ross… ▽ More The Lagrangian dynamics of zonal jets in the atmosphere are considered, with particular attention paid to explaining why, under commonly encountered conditions, zonal jets serve as barriers to meridional transport. The velocity field is assumed to be two-dimensional and incompressible, and composed of a steady zonal flow with an isolated maximum (a zonal jet) on which two or more travelling Rossby waves are superimposed. The associated Lagrangian motion is studied with the aid of KAM (Kolmogorov--Arnold--Moser) theory, including nontrivial extensions of well-known results. These extensions include applicability of the theory when the usual statements of nondegeneracy are violated, and applicability of the theory to multiply periodic systems, including the absence of Arnold diffusion in such systems. These results, together with numerical simulations based on a model system, provide an explanation of the mechanism by which zonal jets serve as barriers to meridional transport of passive tracers under commonly encountered conditions. Causes for the breakdown of such a barrier are discussed. It is argued that a barrier of this type accounts for the sharp boundary of the Antarctic ozone hole at the perimeter of the stratospheric polar vortex in the austral spring. △ Less

Submitted 17 May, 2006; originally announced May 2006.

Comments: Submitted to Journal of the Atmospheric Sciences

Linear waves and baroclinic instability in an inhomogeneous-density layered primitive-equation ocean model

Authors: F. J. Beron-Vera, M. J. Olascoaga, J. Zavala-Garay

Abstract: We consider a multilayer generalization of Ripa's inhomogeneous-density single-layer primitive-equation model. In addition to vary arbitrarily in horizontal position and time, the horizontal velocity and buoyancy fields are allowed to vary linearly with depth within each layer of the model. Preliminary results on linear waves and baroclinic instability suggest that a configuration involving a fe… ▽ More We consider a multilayer generalization of Ripa's inhomogeneous-density single-layer primitive-equation model. In addition to vary arbitrarily in horizontal position and time, the horizontal velocity and buoyancy fields are allowed to vary linearly with depth within each layer of the model. Preliminary results on linear waves and baroclinic instability suggest that a configuration involving a few layers may set the basis for a quite accurate and numerically efficient ocean model. △ Less

Submitted 11 June, 2004; v1 submitted 16 January, 2004; originally announced January 2004.

Comments: To appear in ICTAM04 Abstract Book and CD-ROM Proceedings, ISBN 83-89687-01-1, IPPT PAN, Warsaw

A Didactic Approach to Linear Waves in the Ocean

Authors: F. J. Beron-Vera

Abstract: The general equations of motion for ocean dynamics are presented and the waves supported by the (inviscid, unforced) linearized system with respect to a state of rest are derived. The linearized dynamics sustains one zero frequency mode (called buoyancy mode) in which salinity and temperature rearrange in such a way that seawater density does not change. Five nonzero frequency modes (two acousti… ▽ More The general equations of motion for ocean dynamics are presented and the waves supported by the (inviscid, unforced) linearized system with respect to a state of rest are derived. The linearized dynamics sustains one zero frequency mode (called buoyancy mode) in which salinity and temperature rearrange in such a way that seawater density does not change. Five nonzero frequency modes (two acoustic modes, two inertia--gravity or Poincaré modes, and one planetary or Rossby mode) are also sustained by the linearized dynamics, which satisfy an asymptotic general dispersion relation. The most usual approximations made in physical oceanography (namely incompressibility, Boussinesq, hydrostatic, and quasigeostrophic) are also consider, and their implications in the reduction of degrees of freedom (number of independent dynamical fields or prognostic equations) of, and compatible waves with, the linearized governing equations are particularly discussed and emphasized. △ Less

Submitted 8 January, 2004; originally announced January 2004.

Comments: Submitted to the American Journal of Physics

Multilayer shallow-water model with stratification and shear

Authors: F. J. Beron-Vera

Abstract: The purpose of this paper is to present a shallow-water-type model with multiple inhomogeneous layers featuring variable linear velocity vertical shear and startificaion in horizontal space and time. This is achieved by writing the layer velocity and buoyancy fields as linear functions of depth, with coefficients that depend arbitrarily on horizontal position and time. The model is a generalizatio… ▽ More The purpose of this paper is to present a shallow-water-type model with multiple inhomogeneous layers featuring variable linear velocity vertical shear and startificaion in horizontal space and time. This is achieved by writing the layer velocity and buoyancy fields as linear functions of depth, with coefficients that depend arbitrarily on horizontal position and time. The model is a generalization of Ripa's (1995) single-layer model to an arbitrary number of layers. Unlike models with homogeneous layers the present model is able to represent thermodynamics processes driven by heat and freshwater fluxes through the surface or mixing processes resulting from fluid exchanges across contiguous layers. A model configuration with only one layer has been previously shown to provide: a very good representation of the exact vertical normal modes up to the first internal mode; an exact representation of long-perturbation (free boundary) baroclinic instability; and a very reasonable representation of short-perturbation (classical Eady) baroclinic instability. Here it is shown that substantially more accurate overall results with respect to single-layer calculations can be achieved by considering a stack of only a few layers. A similar behavior is found in ageostrophic (classical Stone) baroclinic instability by describing accurately the dependence of the solutions on the Richardson number with only two layers. △ Less

Submitted 3 December, 2020; v1 submitted 11 December, 2003; originally announced December 2003.

Comments: To appear in Revista Mexicana de Fisica

Deep ocean influence on upper ocean baroclinic instability saturation

Authors: M. J. Olascoaga, F. J. Beron-Vera, J. Sheinbaum

Abstract: In this paper we extend earlier results regarding the effects of the lower layer of the ocean (below the thermocline) on the baroclinic instability within the upper layer (above the thermocline). We confront quasigeostrophic baroclinic instability properties of a 2.5-layer model with those of a 3-layer model with a very thick deep layer, which has been shown to predict spectral instability for b… ▽ More In this paper we extend earlier results regarding the effects of the lower layer of the ocean (below the thermocline) on the baroclinic instability within the upper layer (above the thermocline). We confront quasigeostrophic baroclinic instability properties of a 2.5-layer model with those of a 3-layer model with a very thick deep layer, which has been shown to predict spectral instability for basic state parameters for which the 2.5-layer model predicts nonlinear stability. We compute and compare maximum normal-mode perturbation growth rates, as well as rigorous upper bounds on the nonlinear growth of perturbations to unstable basic states, paying particular attention to the region of basic state parameters where the stability properties of the 2.5- and 3-layer model differ substantially. We found that normal-mode perturbation growth rates in the 3-layer model tend to maximize in this region. We also found that the size of state space available for eddy-amplitude growth tends to minimize in this same region. Moreover, we found that for a large spread of parameter values in this region the latter size reduces to only a small fraction of the total enstrophy of the system, thereby allowing us to make assessments of the significance of the instabilities. △ Less

Submitted 29 July, 2003; originally announced July 2003.

Comments: To appear \emph{in} O. U. Velasco-Fuentes et al. (eds.), \textit{Nonlinear Processes in Geophysical Fluid Dynamics}, Kluwer Academic

Constrained-Hamiltonian Shallow-Water Dynamics on the Sphere

Authors: F. J. Beron-Vera

Abstract: Salmon's nearly geostrophic model for rotating shallow-water flow is derived in full spherical geometry. The model, which results upon constraining the velocity field to the height field in Hamilton's principle for rotating shallow-water dynamics, constitutes an important prototype of Hamiltonian balanced models. Instead of Salmon's original approach, which consists in taking variations of parti… ▽ More Salmon's nearly geostrophic model for rotating shallow-water flow is derived in full spherical geometry. The model, which results upon constraining the velocity field to the height field in Hamilton's principle for rotating shallow-water dynamics, constitutes an important prototype of Hamiltonian balanced models. Instead of Salmon's original approach, which consists in taking variations of particle paths at fixed Lagrangian labels and time, Holm's approach is considered here, namely variations are taken on Lagrangian particle labels at fixed Eulerian positions and time. Unlike the classical quasigeostrophic model, Salmon's is found to be sensitive to the differences between geographic and geodesic coordinates. One consequence of this result is that the $β$ plane approximation, which is included in Salmon's original derivation, is not consistent for this class of model. △ Less

Submitted 29 July, 2003; originally announced July 2003.

Comments: To appear \textit{in} O. U. Velasco-Fuentes et al. (eds.), \textit{Nonlinear Processes in Geophysical Fluid Dynamics}, Kluwer Academic.}

Passive tracer patchiness and particle trajectory stability in incompressible two-dimensional flows

Authors: Francisco J. Beron-Vera, María J. Olascoaga, Michael G. Brown

Abstract: Particle motion is considered in incompressible two-dimensional flows consisting of a steady background gyre on which an unsteady wave-like perturbation is superimposed. A dynamical systems point of view that exploits the action--angle formalism is adopted. It is argued and demonstrated numerically that for a large class of problems one expects to observe a mixed phase space, i.e., the occurrenc… ▽ More Particle motion is considered in incompressible two-dimensional flows consisting of a steady background gyre on which an unsteady wave-like perturbation is superimposed. A dynamical systems point of view that exploits the action--angle formalism is adopted. It is argued and demonstrated numerically that for a large class of problems one expects to observe a mixed phase space, i.e., the occurrence of ``regular islands'' in an otherwise ``chaotic sea.'' This leads to patchiness in the evolution of passive tracer distributions. Also, it is argued and demonstrated numerically that particle trajectory stability is largely controlled by the background flow: trajectory instability, quantified by various measures of the ``degree of chaos,'' increases on average with increasing $|\mathrm{d}ω/\mathrm{d}I|$, where $ω(I)$ is the angular frequency of the trajectory in the background flow and $I$ is the action. △ Less

Submitted 5 July, 2003; originally announced July 2003.

Comments: Submitted to NPG