The emergence of coherent Lagrangian swirls (CLSs) among submesoscale motions in the ocean is illustrated. This is done by applying recent nonlinear dynamics tools for Lagrangian coherence detection on a surface flow realization produced by a dataassimilative submesoscale-permitting ocean general circulation model simulation of the Gulf of Mexico. Both mesoscale and submesoscale CLSs are extracted. These extractions prove the relevance of coherent Lagrangian eddies detected in satellitealtimetry-based geostrophic flow data for the arguably more realistic ageostrophic multiscale flow.Lagrangian coherence | mesoscale eddies | submesoscale eddies | satellite altimetry | high-resolution ocean models A dvances in nonlinear dynamical systems theory over the past few years (1-5) have enabled the discovery of unexpectedly resilient mesoscale material eddies in ocean velocity fields derived from satellite-sensed sea-surface height (SSH) fields (6). Ranging between 50 km and 250 km in diameter, these supercoherent Lagrangian eddies show no breakup or disintegration over a period of several months, in some cases up to 2 y (3, 7-10). Detected in a fashion fully independent from the observer, such Lagrangian eddies are bounded by uniformly stretching fluid (i.e., material) loops that stationarize the material-line-averaged tangential strain functional (3, 4). As a result, these material boundaries or elliptic Lagrangian coherent structures (11) will nearly exactly reassume their initial arc length at the end of the coherence assessment interval. Coupled with area preservation in the incompressible case, the preservation of boundary length renders these special Lagrangian eddies coherent to a previously undocumented and unexpected degree.The presence of one such supercoherent Lagrangian eddy in the Gulf of Mexico has recently been confirmed by independent measurements of satellite-derived chlorophyll and satellitetracked drifter paths (12). On the other hand, while satellitebased measurements are widely used for monitoring mesoscale variability in the ocean (13), they are incapable of resolving submesoscales, ranging from 100 m to a few tens of kilometers (14). There is, therefore, reason to believe (15, 16) that unresolved submesoscale motions, if sufficiently energetic, will erode the boundaries of supercoherent Lagrangian eddies inferred from satellite data. For this reason, the ubiquitous presence of submesoscale turbulence in the ocean (17, 18) has casted some doubt on the practical relevance of supercoherent mesoscale eddies inferred from satellite altimetry.The goal of this paper is to demonstrate that, even in the presence of active submesoscale motions, mesoscale swirling structures will exist, even though their boundaries will show increased filamentation relative to the supercoherent boundaries inferred from satellite altimetry. Such material interfaces will constrain mixing in the presence of multiscale ocean motions, ranging from largely geostrophic (balanced, equilibrated) mesoscale motions to likely ageostrophi...