We employ a recently developed single-trajectory Lagrangian diagnostic tool, the trajectory rotation average , to visualize oceanic vortices (or eddies) from sparse drifter data. We apply the to two drifter data sets that cover various oceanographic scales: the Grand Lagrangian Deployment and the Global Drifter Program. Based on the , we develop a general algorithm that extracts approximate eddy boundaries. We find that the outperforms other available single-trajectory-based eddy detection methodologies on sparse drifter data and identifies eddies on scales that are unresolved by satellite-altimetry.
Meso and submesoscale vortices (or eddies) can trap and transport material over large distances, thereby playing a crucial role in the dynamics of our ecosystem. In order to expand our understanding of the transport of marine tracers, we need to accurately and reliably track the evolution of vortical flow structures. Drifter trajectories represent a valuable but sparse source of information for this purpose. We utilize this information here to evaluate a single-trajectory Lagrangian diagnostic tool that approximates the local material rotation in the flow in a quasi-objective fashion. Our findings on two distinct data sets suggest that this diagnostic accurately highlights material vortices from sparse drifter data.
I. INTRODUCTION
Oceanic vortices (or eddies) are highly energetic coherent flow structures that can transport material over large distances. Studying the dynamics of eddies is key to understanding the dispersion of marine tracers, such as biological nutrients and pollutants.1–5 Lagrangian eddies, generally referred to as elliptic Lagrangian coherent structures (LCSs) in the dynamical systems literature,6 are material objects that trap and transport floating particles over large distances in the ocean in a coherent fashion. The size of such eddies ranges from a few kilometers (submesoscale) to hundreds of kilometers (mesoscale).
Mesoscale eddies have predominantly been inferred from the sea-surface-height (SSH) field derived from satellite-altimetry data.7–9 There is, however, increasing evidence that submesoscale currents on the order of a few kilometers influence the marine ecosystem at least as strongly as mesoscale eddies, especially along coastlines and oceanic fronts.10,11 Despite this observation, submesoscale eddies are rarely studied in detail because their footprint in the SSH field is weak or nonexistent. In contrast to the SSH field, however, drifters follow ocean currents closely and, hence, resolve small-scale features accurately.12 In addition, surface drifters provide information about the local ocean surface velocity field at a very high temporal resolution. Ocean drifters are, however, sparsely distributed in space, rendering most LCS diagnostics inapplicable to their trajectories. Indeed, most mathematically justifiable algorithms for the detection of elliptic LCSs require differentiation with respect to initial conditions,6 which is unfeasible for sparse drifter data.
Alternatively, elliptic LCSs can be viewed as coherently evolving sets of material trajectories in space. Based on this view, a wide range of clustering methods for extracting elliptic LCSs have been proposed (see Ref. 13 for a review). These objective methods differ depending on the specific clustering algorithm and the distance metric they employ to define clusters. Common algorithms include fuzzy clustering,14 spectral graph methods,15–17 and density-based clustering methods.18 Their results, however, rely on user defined input parameters such as the number of clusters, which a priori determines the possible elliptic LCS to be detected by the method in the domain. This limitation becomes even more pronounced for sparse drifters for which the number of eddies is a priori unknown.
To extract elliptic LCSs from drifter data, one is forced to rely on features of a single trajectory, such as trajectory length, velocity, acceleration, and curvature. These features, however, are all inherently non-objective (frame-dependent) quantities whereas LCSs, as material objects, are objective, i.e., indifferent to coordinate frame changes.6 Common single-trajectory methods, such as the absolute dispersion,19 trajectory length,20,21 Lagrangian spin,22 maximal trajectory extent,23 trajectory complexity,24 and wavelet ridge analysis,25,26 are limited in capturing elliptic LCSs as they are either non-objective or lack direct physical connection to material deformation in the fluid. Accordingly, while most of these methods were originally introduced to visualize LCSs from truly sparse trajectories, their application to drifter data has remained rare (see Ref. 27 for a review). Notable exceptions include the fully automated looper detection algorithms based on the Lagrangian spin12,28–30 and the wavelet ridge regression.31 Both methods seek to extract oscillatory motions from a time-series and lead to comparable results. The methodology employed by the authors in Refs. 12 and 28–30 is an aggregate measure of rotation within a time-series of the Lagrangian spin, whereas Lilly and Pérez-Brunius31 quantify the instantaneous oscillatory motion of the velocity using signal processing techniques. As a consequence, elliptic LCSs are visualized through looping (or oscillatory) trajectory segments.
In order to distinguish looping from non-looping trajectory segments, a user-defined threshold is inevitably required. In a practical setting, however, differentiating between looping and non-looping trajectories proves to be challenging and hence finding an appropriate parameter is far from trivial. A common approach is to manually tune a threshold value for looping based on a subset of trajectories. Although this provides a natural and intuitive way to characterize elliptic LCSs, a considerable amount of information is lost when one discards non-looping trajectory segments. It would be more desirable to retain all trajectory information and associate to each trajectory a Lagrangian diagnostic related to material rotation in the flow. This rotational diagnostic could then be plotted over evolving trajectory positions. This approach would provide a qualitative overview of individual vortical flow structures, without relying on any chosen threshold.
Here, we propose to identify eddies from sparse drifter data using the recently introduced adiabatically quasi-objective single-trajectory diagnostics in Ref. 32. These diagnostics closely approximate appropriate objective LCS diagnostics in frames satisfying conditions that typically hold in geophysical flows. We show that the adiabatically quasi-objective trajectory rotation average ( ) reveals elliptic LCSs (material eddies) at meso- and submesoscales from sparse drifter trajectory data. We additionally compare the vortical flow features extracted from the drifter-based computation with those obtained from other available Lagrangian single-trajectory diagnostics such as the trajectory length20,21 and the Lagrangian spin.12,22 We further validate the extracted features with respect to Lagrangian averaged vorticity deviation (LAVD)33 computed from geostrophic ocean velocity fields derived from satellite-altimetry data (AVISO).
II. DATA
A. Satellite-altimetry ocean-surface current product (AVISO)
B. Drifter data sets
While satellite-based altimetry yields mesoscale velocity fields, surface drifter observations provide direct estimates for the local surface velocity field. In order to illustrate the range of applications of the adiabatically quasi-objective single-trajectory rotation measure introduced in Ref. 32, we will focus on two drifter data sets: the Grand Lagrangian Deployment (GLAD) and the Global Drifter Program (GDP).
1. Grand Lagrangian Deployment (GLAD)
2. Global Drifter Program (GDP)
3. Drifter data preprocessing
At mid-latitudes, inertial oscillations have time-scales of 1–2 days, whereas the dominant period of submesoscale and mesoscale eddy motion is above 5 days.38 The characteristic time of submesoscale and mesoscale features, thus, greatly exceeds the inertial period. Hence, the anticyclonic looping arising from high-frequency diurnal inertial oscillations does not influence the relative dispersion of drifters in the submesoscale and mesoscale regimes.41 In this paper, we target coherent structures whose time-scales are far above the period of inertial oscillations.
The time interval of definition of a finite-time, temporally aperiodic dynamical system inherently determines the range of features one can identify in that system. Submesoscale eddies have characteristic time-scales of around one week, whereas mesoscale structures arise over months. In contrast, coherent structures tied to inertial oscillations have a temporal range of a few hours. Therefore, in the exploration of short term features, inertial oscillations should not be filtered out of the drifter trajectories.
III. METHODS
A. Trajectory rotation average
We associate to the trajectory over the time interval its corresponding value, where . Reconstructing the -field from sparse drifter trajectories via scattered interpolation allows visualization of vortical flow structures. The only parameter involved in the reconstruction of the -field is the choice of the scattered interpolation method. As vortices tend to be elliptic geometric objects, we use linear radial basis function interpolation that favors such structures (see Appendix B for further details). We then apply a spatial average filter of size in order to suppress noise but still retain small-scale features. The size of this spatial filter coincides with the spatial resolution of the AVISO data set. Local maxima in the -field mark centers of high material rotation and, hence, indicate underlying rotational flow features.
B. Trajectory length
C. Lagrangian spin
D. Lagrangian averaged vorticity deviation
E. Validation of eddy boundary extraction algorithm on numerical ocean model
and evaluated on the full and progressively subsampled AVISO data sets in the Agulhas region. (a) -field with full resolution ( ). The red curves denote the boundaries of the three -eddies. [(b)–(d)] Reconstructed -field from 2000, 200, and 20 randomly selected trajectories. The white curves indicate boundaries of the -eddies. The black triangles indicate the final trajectory positions.
and evaluated on the full and progressively subsampled AVISO data sets in the Agulhas region. (a) -field with full resolution ( ). The red curves denote the boundaries of the three -eddies. [(b)–(d)] Reconstructed -field from 2000, 200, and 20 randomly selected trajectories. The white curves indicate boundaries of the -eddies. The black triangles indicate the final trajectory positions.
We vary the number of trajectories from to and randomly subsample multiple times in order to obtain a statistical estimate of the DICE coefficient for each trajectory density. In Fig. 2, we plot the averaged DICE coefficient and its error bar as a function of two quantities: the averaged trajectory spacing normalized by the eddy diameter [Fig. 2(a)] and the trajectory density [Fig. 2(b)]. The eddy diameter of the three -eddies in the Agulhas region is around . The averaged trajectory spacing per eddy diameter relates the averaged drifter spacing to the eddy length scales we seek to resolve. Decreasing the drifter density (or, equivalently, increasing the averaged spacing) leads to a drop in the DICE coefficient. Even for very sparse data sets (e.g., ), it is possible to approximately identify coherent structures with maximum DICE coefficients of around . Note, however, that the results inevitably depend on the drifter distribution. If no drifters are inside an eddy, then we do not expect the eddy to be visible in the -field.
Averaged DICE coefficient plotted as a function of the averaged trajectory spacing per eddy diameter (a) and trajectory density (b). The error bars indicate the minimum and maximum DICE coefficient attained within the corresponding set of trajectories.
Averaged DICE coefficient plotted as a function of the averaged trajectory spacing per eddy diameter (a) and trajectory density (b). The error bars indicate the minimum and maximum DICE coefficient attained within the corresponding set of trajectories.
GLAD drifters (to be analyzed in Sec. IV A) have an average spacing of and a trajectory density of around . The diameter of the smallest submesoscale eddies observed during the GLAD experiment is around and the averaged drifter spacing per eddy diameter is . GDP drifters (to be analyzed in Sec. IV B) are on average separated by with an estimated trajectory density of . Mesoscale eddies have characteristic widths ranging from to 12 and, hence, the averaged spacing per eddy diameter ranges from to . Based on the statistical estimates of the DICE coefficient in Fig. 2, we expect to detect submesoscale and mesoscale eddies in the GLAD and GDP data sets with reasonable confidence.
IV. RESULTS
A. Grand Lagrangian deployment (GLAD)
Drifters released in the GLAD experiment have mostly been used to study dispersion in the ocean over a range of scales.38,39,41 Furthermore, the authors in Refs. 34 and 35 found a correlation between the evolution of the drifters and nearby attracting LCSs extracted from the geostrophic velocity field. Specifically, drifter behavior agreed with the tiger tail pattern inferred from the chlorophyll distribution shown in Fig. 3. A chlorophyll plume extended over more than 100 km from the outlet of the Mississippi River into the open sea and coincided with an attracting LCS.34 The attracting LCS (continuous black line in Fig. 3) is computed from the geostrophic velocity field using backward trajectories over the time interval according to Ref. 56. Drifters in the proximity of the tip of the tiger tail organized themselves into long filaments along the attracting LCS. Some of the drifters (in red) additionally exhibited some clustering, suggesting the presence of an elliptic LCS close to the chlorophyll front.
Evolution of GLAD-drifters overlayed on the 8-day composite chlorophyll concentration (12 August 2012). Most drifters (white triangles) formed a long filament along the chlorophyll plume (attracting LCS) extending deep into the Gulf of Mexico. Drifters additionally forming a visible tight cluster are colored red. The blue cross denotes the outlet of the Mississippi River. The black line corresponds to the attracting LCS. (a) . (b) . (c) . (d) .
Evolution of GLAD-drifters overlayed on the 8-day composite chlorophyll concentration (12 August 2012). Most drifters (white triangles) formed a long filament along the chlorophyll plume (attracting LCS) extending deep into the Gulf of Mexico. Drifters additionally forming a visible tight cluster are colored red. The blue cross denotes the outlet of the Mississippi River. The black line corresponds to the attracting LCS. (a) . (b) . (c) . (d) .
Here, we seek to visualize vortices by reconstructing the -field from drifter data over the interval using linear radial basis function interpolation. The , shown in Fig. 4, is plotted with respect to the final drifter positions at time . This plot reveals multiple rotational features marked by local maxima of the surrounded by a dense set of closed and almost convex contours. For comparison, we include three further Lagrangian eddy diagnostics: the LAVD computed from geostrophic velocity currents, the trajectory length (also called -function), and the Lagrangian spin computed from drifter data.
Different material eddy diagnostics computed from the GLAD data set or the AVISO data set. The dashed black lines correspond to the isobaths and highlight the coastal area. (a) Reconstructed -field plotted with respect to the position of the drifters at time . White crosses indicate local maxima. (b) -field computed from AVISO plotted at time . The closed red curve denotes the eddy boundary extracted from the -field. (c) Reconstructed -field plotted with respect to the position of the drifters at time . Red crosses indicate local minima. (d) Cyclonic (red) and anticyclonic (blue) loopers.
Different material eddy diagnostics computed from the GLAD data set or the AVISO data set. The dashed black lines correspond to the isobaths and highlight the coastal area. (a) Reconstructed -field plotted with respect to the position of the drifters at time . White crosses indicate local maxima. (b) -field computed from AVISO plotted at time . The closed red curve denotes the eddy boundary extracted from the -field. (c) Reconstructed -field plotted with respect to the position of the drifters at time . Red crosses indicate local minima. (d) Cyclonic (red) and anticyclonic (blue) loopers.
1. Coastal area
We start by focusing on the area at the outlet of the Mississippi River on the continental shelf. Coastal flows are regions of intense mixing,45 often characterized by high potential vorticity and strong horizontal shear, which are responsible for the formation of small-scale eddies.57,58 Anticyclonic looping segments extracted from the Lagrangian spin parameter confirm the presence of small-scale vortices at the outlet of the Mississippi River [Fig. 4(d)]. Evidence for the existence of submesoscale elliptic features and strong mixing areas on the continental shelf is also found in the -field [Fig. 4(a)]. Vortices are revealed in the as local maxima surrounded by a dense set of closed curves, indicating abrupt spatial variations. The local maximum of the at (88.6 W, 28.1 N) is surrounded by a dense set of closed contours, thereby indicating high spatial gradients in the -field. Close to the vortical flow features identified by the in the coastal area, the -function displays two local minima, respectively, at (88.6 W, 28.2 N) and (88.5 W,27.8 N) [Fig. 4(c)]. Similarly to the and the Lagrangian spin parameter, the -function also indicates the existence of coastal eddies. However, compared to the , which displays a unique local maximum surrounded by a dense set of closed and convex curves, the -function displays multiple local minima surrounded by sparse level sets, thereby indicating a region without sharp minima.
2. Open sea
The highlights another family of elliptic LCSs close to the chlorophyll plume extending from the outlet of the Mississippi River to the open ocean [Fig. 4(a)]. The two local maxima marking elliptic LCSs are located at (87.5 W, 26.8 N) and (87.3 W, 27.4 N). Due to the close proximity of these extrema, it is, however, a priori unclear whether they highlight two separate eddies or whether they are part of the same mesoscale eddy. The cyclonic looping segments confirm the presence of vortical flow features in this area [red trajectories in Fig. 4(d)], but do not specifically highlight eddy boundaries. Inspection of the cyclonic trajectories suggests that they are originally part of two distinct eddies as they come from two separate regions. Extracting looping trajectory segments from the time-series of the Lagrangian spin is a powerful methodology. However, discriminating between loopers and non-loopers requires an ad hoc choice of the minimum sustained looping period of a trajectory. In contrast, the retains all trajectory information and allows visualization of vortical flow structures from sparse drifter data using a scalar diagnostic.
The investigated elliptic LCSs highlighted by the and the Lagrangian spin are also visible in the , which displays a nearly convex vortical flow feature [closed red curve in Fig. 4(b)]. Compared to the , however, which indicates two separate albeit closely located local maxima, the computation clearly reveals a single mesoscale eddy.
In contrast to the aforementioned methods, the features resulting from the trajectory length diagnostic are far less pronounced [Fig. 4(c)]. Indeed, the -function displays multiple local minima, which are only partially correlated with the eddies suggested by the other methods. Furthermore, local minima in the -function are not sharp, as indicated by the sparsity of the surrounding level curves. Hence, extracting eddy-related features from the -function is challenging as there exist no distinguishable sharp and closed boundaries surrounding the local minima.
3. Eddy dynamics
In order to investigate the formation and evolution of the elliptic LCSs inferred from the drifter-based , we proceed by quasi-materially advecting the distribution over the time interval (Fig. 5). True material advection would require a spatiotemporally well-resolved velocity field. In our setting, however, the velocity field is only sparsely known and, hence, the advected structures are inherently non-material: At every time step, the -field must be approximated from the current drifter distribution. As a consequence, the evolution of the extracted eddy boundaries is not exactly Lagrangian. This implies that a set of drifters may not stay inside an eddy over the full time interval. For the same reason, the extracted eddies can potentially merge.
Quasi-materially advected for the GLAD data set over t [222 doy, 229 doy] [(a)–(h)]. The black line denotes the AVISO-based attracting LCS. Closed white curves mark eddy boundaries extracted from the drifter-based using the algorithm described in Appendix A. Closed red curves indicate eddy boundaries obtained from the AVISO-based computation. Triangles indicate the positions of the drifters.
Quasi-materially advected for the GLAD data set over t [222 doy, 229 doy] [(a)–(h)]. The black line denotes the AVISO-based attracting LCS. Closed white curves mark eddy boundaries extracted from the drifter-based using the algorithm described in Appendix A. Closed red curves indicate eddy boundaries obtained from the AVISO-based computation. Triangles indicate the positions of the drifters.
The closed white curves in Fig. 5 indicate eddy boundaries extracted from the using the algorithm described in Appendix A. The red closed curve denotes the truly materially advected vortex boundary extracted from the at time . The -based eddy at is materially advected using the geostrophic velocity , whereas the eddy boundaries inferred from the drifter-based are quasi-materially advected along drifter trajectories.
At , the suggests the presence of several small-scale vortices at the outlet of the Mississippi River and at open sea [Fig. 5(a)]. The submesoscale eddies close to the outlet of the Mississippi River remain trapped in coastal areas and eventually merge into a larger vortical flow feature. Over the time interval , the -based eddy does not coincide with any of the eddies inferred from the [Figs. 5(a)–5(d)]. The white eddy initially located at approximately (88.5 W, 27.5 N) is associated with the clustered red drifters identified in Fig. 3, thereby confirming the existence of the submesoscale eddy along the chlorophyll front, which agrees with the observation put forward in Refs. 35 and 34. The eddy develops along the attracting LCS and then slowly detaches away from the oceanic front [Figs. 5(a)–5(f)]. Drifters are thereby coherently transported from the coastline into the open sea. The drifter-based eddy evolving along the AVISO-based attracting LCS eventually merges with the submesoscale eddy originally located at (86.5 W, 27.25 N) to form a larger mesoscale eddy at [Fig. 5(g)]. Hence, elliptic LCSs generated along oceanic fronts offer a transport mechanism for particles, carrying material over long distances away from the coastline into the open sea. We point out that intersections between attracting and elliptic LCS (such as at ) are physically inconsistent as they imply contradicting material response. We attribute this inconsistency to the fact that the attracting LCS is computed from whereas the white elliptic LCSs are computed from drifter velocities.
Toward the end of the advection process, the -based elliptic LCSs [red curve in Fig. 5(h)] approximates the mesoscale eddy inferred from the [white curve centered at approximately (87.5 W, 27 N) in Fig. 5(h)]. The red eddy shows no degree of filamentation and remains coherent over the full time interval. The white eddy results from the vortex merger between two smaller eddies and is significantly larger than the red eddy. At , the circumference of the -based white eddy and the -based red eddy is 298.5 and . As the -based eddy is a materially advected closed curve, it can neither split nor merge with any other materially advected curve. Hence, by construction, the is unable to capture vortex mergers. In contrast, the quasi-materially advected, -based eddies can merge into larger eddies. This follows because the evolving eddies are not purely Lagrangian, given that the computation of the eddy boundary from the -field is independently carried out at each time step.
evaluated both on GLAD-drifters and on AVISO. The white and red eddies, respectively, indicate the eddies extracted from the drifter-based -field and the AVISO-based -field. The dashed black lines correspond to the isobaths . The black line denotes the AVISO-based attracting LCS. (a) Drifter-based -field at time . White, dashed curves denote a subset of GLAD drifter trajectories contained currently in the same -eddy. (b) AVISO-based -field at time . Red dashed curves denote a subset of AVISO trajectories contained currently in the same -eddy.
evaluated both on GLAD-drifters and on AVISO. The white and red eddies, respectively, indicate the eddies extracted from the drifter-based -field and the AVISO-based -field. The dashed black lines correspond to the isobaths . The black line denotes the AVISO-based attracting LCS. (a) Drifter-based -field at time . White, dashed curves denote a subset of GLAD drifter trajectories contained currently in the same -eddy. (b) AVISO-based -field at time . Red dashed curves denote a subset of AVISO trajectories contained currently in the same -eddy.
The quasi-material advection of the -based eddies sheds new light on the origin and the formation of the mesoscale eddy detected from the AVISO-based -field. Despite the minimal amount of data, the reveals eddy dynamics hidden in the and .
B. Global Drifter Program (GDP)
In our second example, we focus on a set of drifters in the western North-Atlantic. This oceanic region is characterized by a strong and persistent formation of eddies arising from the meanders of the Gulf Stream.62,63 On 4 October 2006, a floating sargassum patch was detected by the Medium Resolution Imaging Spectrometer (MERIS) on Envisat. This feature has a spiralling shape that is also visible from satellite-altimetry data.2 This floating sargassum patch is visualized in Fig. 7(b) using the Maximum Chlorophyll Index (MCI).64 Due to frequent cloud coverage, such clear snapshots of floating material in the ocean are very rare.
Lagrangian diagnostics evaluated on a subset of the GDP data set. (a) Reconstructed -field plotted with respect to the position of the drifters (triangles) at time . White cross indicates a unique local maximum in the -field. (b) Eddy boundary superimposed on the floating sargassum concentration inferred from the Maximum Chlorophyll Index (MCI) at . (c) Reconstructed -field plotted with respect to the position of the drifters (triangles) at time . White cross indicates a unique local maximum in the -field. (d) Cyclonic (red) loopers.
Lagrangian diagnostics evaluated on a subset of the GDP data set. (a) Reconstructed -field plotted with respect to the position of the drifters (triangles) at time . White cross indicates a unique local maximum in the -field. (b) Eddy boundary superimposed on the floating sargassum concentration inferred from the Maximum Chlorophyll Index (MCI) at . (c) Reconstructed -field plotted with respect to the position of the drifters (triangles) at time . White cross indicates a unique local maximum in the -field. (d) Cyclonic (red) loopers.
We use the to extract the eddy highlighted by the spiralling sargassum patch described in Ref. 2. We also compute two further single trajectory metrics: the trajectory length diagnostic [Fig. 7(c)] and the looping segments derived from the Lagrangian spin [Fig. 7(d)]. We include a snapshot of the floating sargassum patch that reveals the presence of an elliptic LCS [Fig. 7(b)]. Similarly to the MCI, the reveals a Lagrangian eddy centered at (67.6 W, 37.1 N). Indeed, this eddy is visible as a distinguished local maximum in the -field surrounded by a dense set of closed and convex contours [the zoomed inset of Fig. 7(a)]. The white eddy boundary is extracted from the using the algorithm proposed in Appendix A with the same parameters presented in Sec. IV A 3. The extracted eddy boundary underestimates the size of the sargassum patch, but correctly approximates the location [Fig. 7(b)]. The cyclonic looping exhibited by the two trajectories inside the eddy additionally confirms the presence of an elliptic LCS [Fig. 7(d)] but does not provide a specific estimate for the eddy boundary. The trajectory length diagnostic shows a nearby maximum [the zoomed inset of Fig. 7(c)], which is inconsistent with the generally suggested footprint (local minima) for a coherent eddy in the -field (see Sec. III B).
Figure 8 shows the -field with AVISO trajectories (red) launched backwards in time from true drifter positions at time . The mesoscale eddy inferred from the sparse drifter distribution (white curve in the zoomed inset of Fig. 8) is in close agreement with the -eddy (red curve in the zoomed inset of Fig. 8). As expected, drifter and AVISO trajectories show similar looping behavior inside the identified coherent structure, while they largely differ in other regions.
AVISO-based trajectories (red) and true drifter trajectories (white) superimposed on LAVD-field. The closed white and red curves in the zoomed inset, respectively, indicate the -and -eddy. Triangles denote drifter positions.
AVISO-based trajectories (red) and true drifter trajectories (white) superimposed on LAVD-field. The closed white and red curves in the zoomed inset, respectively, indicate the -and -eddy. Triangles denote drifter positions.
V. CONCLUSION
Lagrangian eddies (elliptic LCSs) are material objects responsible for the transport of floating particles over large distances in the ocean. They are, by their definition, frame-indifferent and, thus, can only be reliably deduced from objective feature extraction methods. The local ocean velocity is most accurately observed from float trajectory data, which, however, is inherently non-objective, representing an inconsistency that available eddy detection methods for sparse trajectories do not address. Those methods typically describe eddies by extracting the looping segments of a trajectory, but their definition of looping depends on the frame of the observer. Furthermore, looping segments of a trajectory are most commonly described by these methods in a statistical sense and, hence, are not geared toward highlighting individual Lagrangian eddy boundaries with high accuracy.
In this paper, we have proposed to tackle this inconsistency from a dynamical systems perspective by applying the adiabatically quasi-objective diagnostic32,46 to sparse drifter data sets. The approximates an objective measure of material rotation in frames satisfying specific conditions that generally hold in the ocean. We have found that vortical flow features are related to regions of high local material rotation identified as local maxima in the -field. The highlights both submesoscale and mesoscale vortices from sparse drifter data, as demonstrated in our two examples. Furthermore, it also succeeds in characterizing the mixing and stirring processes in coastal flow regions and captures the merger of two originally distinct eddies. Contrary to other single trajectory diagnostics, both the and the are physically related to the local material rotation in the flow. In contrast to the , which correctly highlights vortical flow features in sufficiently dense velocity data, the can be applied to arbitrarily sparse drifter data, given its lack of dependence on nearby drifter trajectories. Hence, it incorporates valuable drifter data into the analysis of oceanic coherent structures in a physically and mathematically justifiable way. This proves to be especially useful in ocean regions where satellite-altimetry data do not unravel the true ocean dynamics.
Importantly, the looper segments extracted from the Lagrangian spin coincide with the features identified in the . However, a spaghetti plot of looping trajectory segments does not immediately reveal transport barriers and eddy boundaries. Furthermore, potentially valuable information is lost when we discard non-looping trajectory segments based on a manually tuned threshold parameter. The trajectory length diagnostic20,21 also tended to show either minima or maxima near the material eddies highlighted by the local maxima of the -field. This variation in extremum types creates ambiguity in using the trajectory length as a stand-alone indicator for detecting elliptic LCSs from sparse data sets.
Apart from the visual inspection of the reconstructed -field, we have additionally presented an algorithm to extract approximate eddy boundaries from sparse drifter data. As vortical flow features are indicated by blobs close to local maxima in the -field, the proposed method resembles a blob detection algorithm. Passing from a continuous scalar diagnostic field to a set of closed curves inevitably requires introducing user-defined parameters. All in all, however, the number of free parameters is comparable to other multi-trajectory Lagrangian eddy detection methods.33,65 This is noteworthy as these algorithms were originally designed assuming knowledge of the underlying velocity field.
ACKNOWLEDGMENTS
The authors acknowledge financial support from Priority Program SPP 1881 (Turbulent Superstructures) of the German National Science Foundation (DFG).
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Alex P. Encinas Bartos: Data curation (equal); Investigation (equal); Writing – original draft (equal). Nikolas O. Aksamit: Supervision (equal); Writing – review & editing (equal). George Haller: Conceptualization (equal); Funding acquisition (equal); Supervision (equal); Writing – review & editing (equal).
DATA AVAILABILITY
The AVISO geostrophic current velocity product used in this study “Global Ocean Gridded L4 Sea Surface Heights and Derived Variables Reprocessed” is freely available and is hosted by the Copernicus Marine Environment Monitoring Service.36 The GDP drifter data is openly available in the NOAA Global Drifter Program data set.66 The GLAD drifter data is distributed by the Gulf of Mexico Research Initiative Information and Data Cooperative (GRIIDC).67 The attracting LCS and LAVD computations have been performed with the software package TBarrier. Jupyter notebooks (together with drifter data) implementing the methods used here are available under Github/EncinasBartos. These notebooks can readily be applied to any sparse trajectory data set.
APPENDIX A: EDDY BOUNDARY EXTRACTION ALGORITHM
While various eddy extraction algorithms have been presented in Refs. 16, 17, 50, and 68, these methods assume a trajectory density that is generally unavailable for drifter trajectories in the ocean. Here, we propose an algorithm that extracts approximate eddy boundaries from the topology of the reconstructed -field. Vortices are then identified by this algorithm as local maxima of the -field surrounded by a dense set of closed and convex curves characterized by high spatial gradients.
Extraction of approximate eddy boundaries from -field
Input: Trajectories over the time interval . |
1. Reconstruct -field at time using linear radial basis interpolation. It is recommended to additionally filter the resulting using a spatial average filter to reduce noise. |
2. Find local maxima of which are above a threshold . |
3. Compute for each closed level set surrounding a local maximum, the averaged along the level set. |
4. Find closed level set with the highest averaged , which additionally |
(a) has at least one local maximum of in its interior. |
(b) contains at least -trajectories. |
5. Take the convex hull of all the selected closed curves. |
6. If two or more convex closed curves intersect, then take the union of these curves. |
7. Take the convex hull of the resulting closed curves. |
Output: Approximate eddy boundaries at time . |
Input: Trajectories over the time interval . |
1. Reconstruct -field at time using linear radial basis interpolation. It is recommended to additionally filter the resulting using a spatial average filter to reduce noise. |
2. Find local maxima of which are above a threshold . |
3. Compute for each closed level set surrounding a local maximum, the averaged along the level set. |
4. Find closed level set with the highest averaged , which additionally |
(a) has at least one local maximum of in its interior. |
(b) contains at least -trajectories. |
5. Take the convex hull of all the selected closed curves. |
6. If two or more convex closed curves intersect, then take the union of these curves. |
7. Take the convex hull of the resulting closed curves. |
Output: Approximate eddy boundaries at time . |
Passing from a continuous scalar field to a set of discrete closed curves representing eddy boundaries inevitably requires introducing threshold parameters. There are two main parameters involved in Algorithm 1. The first user-defined quantity aids the identification of the local maxima in the -field. As we identify vortical flow features with regions of high , local maxima below a predefined threshold are neglected. Additionally, we also need to specify the minimum number of drifters inside an eddy. As elliptic LCSs are often observed via a dense clustering of multiple drifters, this parameter is generally set to be greater than 1. In this work, the parameters are consistently set to , where is the global maximum in the -field and , as we expect multiple drifters to be part of the same coherent structure. Hence, our eddy extraction algorithm detects elliptic LCSs if at least two trajectories are contained inside an eddy. By construction, eddy boundaries are closed convex curves characterized by sharp gradients. These curves do not necessarily coincide with level sets of the .
APPENDIX B: SENSITIVITY ANALYSIS WITH RESPECT TO THE INTERPOLATION METHOD
The features of the reconstructed -field clearly depend on the employed interpolation scheme, but we expect the topology of the -field to be robust with respect to the interpolation method in the vicinity of sharp maxima. Here, we verify the persistence of such local maxima and the extracted -eddies with respect to common interpolants: linear radial basis function (rbf), linear ( -interpolant), and natural neighbor ( -interpolant). The reconstructed scalar fields will all be equally pre-and post-processed: inertial oscillations are removed from the drifter trajectories and a spatial averaging filter of size is applied to the reconstructed -field.
Figure 9 shows the reconstructed -field using linear rbf, linear, and natural neighbor interpolation for the GLAD drifters. For comparison, we also included the raw plot. Two major eddies are detected independent of the interpolation scheme [see Figs. 9(a)–9(c)]. Both the linear rbf and natural neighbor interpolation reveal three similar local maxima [white crosses/squares in Figs. 9(a) and 9(c)]. On average, the white crosses and squares are only separated by ( ) and the average distance of the local maxima to the closest drifter is around ( ) in both cases.
Comparison of -fields obtained from different interpolation functions with the raw plot for the GLAD data set at time . Closed white curves indicate eddies extracted from the reconstructed -field. White crosses, circles, and squares denote local maxima. Triangles indicate the position of the drifters. (a) Linear rbf, (b) linear ( -interpolation), (c) natural neighbor ( -interpolation), and (d) scattered, raw . Color and size of the triangles are related to the -value associated with each drifter.
Comparison of -fields obtained from different interpolation functions with the raw plot for the GLAD data set at time . Closed white curves indicate eddies extracted from the reconstructed -field. White crosses, circles, and squares denote local maxima. Triangles indicate the position of the drifters. (a) Linear rbf, (b) linear ( -interpolation), (c) natural neighbor ( -interpolation), and (d) scattered, raw . Color and size of the triangles are related to the -value associated with each drifter.
The linear interpolant reveals three further local maxima closely located to those suggested by the linear rbf and natural neighbor interpolation scheme [white circles in Fig. 9(b)]. These additional local maxima, however, are not strong enough to lead to the emergence of new eddies. Furthermore, local maxima are located in regions with dense -values [see Fig. 9(d)]. All three interpolation schemes suggest the emergence of two robustly persisting eddies.
We perform an analogous sensitivity analysis on the reconstructed -field for the GDP drifters at time (see Fig. 10). Local maxima of the linear rbf and natural neighbor interpolation are in close agreement with each other, whereas the linear interpolant introduces two further local maxima. Again, these additional local maxima are weak and do not lead to the detection of additional eddies. In summary, the eddy boundary extraction algorithm consistently reveals a mesoscale eddy independent of the interpolation scheme.
Overall, we find that the eddy location does not strongly depend on the interpolation scheme at least in the two data sets considered in this paper. However, the exact shape and size of the eddies vary as a function of the interpolation method. In very sparse data sets, such as the GDP-data set, the identified eddy area may vary substantially under changes in the interpolation method. The location of the largest local maximum signaling the presence of an eddy, however, is fairly robust with respect to the interpolation method. Due to its inherent radial symmetry, linear rbf interpolation tends to favor elliptic structures. In contrast, linear and natural neighbor interpolation leads to sharper and non-smooth eddy boundaries.