Submesoscale Eddies Research Articles

Modulation of Phytoplankton Uptake by Mesoscale and Submesoscale Eddies in the California Current System

AbstractEddies play a crucial role in shaping ocean dynamics by affecting material transport, and generating spatio‐temporal heterogeneity. However, how eddies at different scales modulate biogeochemical transformation rates remains an open question. Applying a multi‐scale decomposition to a numerical simulation, we investigate the respective impact of mesoscale and submesoscale eddies on nutrient transport and biogeochemical cycling in the California Current System. First, the non‐linear nature of nutrient uptake by phytoplankton results in a 50% reduction in primary production in the presence of eddies. Second, eddies shape the vertical transport of nutrients with a strong compensation between mesoscale and submesoscale. Third, the eddy effect on nutrient uptake is controlled by the covariance of temperature, nutrient and phytoplankton fluctuations caused by eddies. Our results shed new light on the tight interaction between non‐linear fluid dynamics and ecosystem processes in realistic eddy regimes, which remain largely under‐resolved by global Earth system models.

Different vertical heat transport induced by submesoscale motions in the shelf and open sea of the northwestern South China Sea

Submesoscale currents are known to be associated with strong vertical velocities (O (10) m/day), regulating the redistributions of energy and matter balances. The northern South China Sea (SCS) is fulfilled with submesoscale motions, which might induce strong vertical heat transport (VHT). We set up a 1-km horizontal resolution Massachusetts Institute of Technology General Circulation Model (MITgcm) to study the seasonal variations in submesoscale vertical heat transport in shelf regions and open seas. Spectrum analysis shows that the spatial scale separating submesoscale and mesoscale motions are 14 and 30 km for the shelf and open regions, respectively. The submesoscale VHT in the shelf region is one order of magnitude larger than that in the open ocean. The former has the largest value in summer and winter, which might be induced by summer upwelling and winter downwelling, while the latter is strongest in winter and weakest in summer in open regions. The submesoscale VHT also appears to have intra-seasonal variations and might be attributed to the disturbances of tropical cyclones and life stages of submesoscale eddies. The submesoscale VHT is strongest in the pregeneration phase of the eddies, and the maximum VHT belt has an entrainment type at the developing and mature stages. The chlorophyll-a concentration also has the same temporal variation as the different life-stage of eddies. This study provides local VHT induced by submesoscale motions, which is expected to improve our understanding of submesoscale air–sea interactions and their biological effects.

Surface and subsurface oceanographic features drive forage fish distributions and aggregations: Implications for prey availability to top predators in the US Northeast Shelf ecosystem.

Forage fishes are a critical food web link in marine ecosystems, aggregating in a hierarchical patch structure over multiple spatial and temporal scales. Surface-level forage fish aggregations (FFAs) represent a concentrated source of prey available to surface- and shallow-foraging marine predators. Existing survey and analysis methods are often imperfect for studying forage fishes at scales appropriate to foraging predators, making it difficult to quantify predator-prey interactions. In many cases, general distributions of forage fish species are known; however, these may not represent surface-level prey availability to predators. Likewise, we lack an understanding of the oceanographic drivers of spatial patterns of prey aggregation and availability or forage fish community patterns. Specifically, we applied Bayesian joint species distribution models to bottom trawl survey data to assess species- and community-level forage fish distribution patterns across the US Northeast Continental Shelf (NES) ecosystem. Aerial digital surveys gathered data on surface FFAs at two project sites within the NES, which we used in a spatially explicit hierarchical Bayesian model to estimate the abundance and size of surface FFAs. We used these models to examine the oceanographic drivers of forage fish distributions and aggregations. Our results suggest that, in the NES, regions of high community species richness are spatially consistent with regions of high surface FFA abundance. Bathymetric depth drove both patterns, while subsurface features, such as mixed layer depth, primarily influenced aggregation behavior and surface features, such as sea surface temperature, sub-mesoscale eddies, and fronts influenced forage fish diversity. In combination, these models help quantify the availability of forage fishes to marine predators and represent a novel application of spatial models to aerial digital survey data.

Seasonality of Upper-Ocean Vertical Eddy Buoyancy Flux in the Kuroshio Extension

Abstract Oceanic mesoscale and submesoscale eddies produce a pronounced vertical buoyancy flux, playing an important role in ocean restratification. This study used a 1-km ocean simulation to investigate the seasonality of the vertical eddy buoyancy flux (VEBF) in the Kuroshio Extension as well as its underlying dynamics. The simulated VEBF in the upper 200 m over the Kuroshio Extension has a pronounced seasonal cycle. The winter VEBF peaks in the mixed layer, whereas the summer VEBF has a much smaller magnitude but a more complicated vertical structure with a narrow peak in the shallow mixed layer and a broader and stronger peak in the seasonal thermocline. The baroclinic instability (BCI), frontogenesis, and turbulent thermal wind (TTW) balance all contribute to the VEBF seasonal cycle. In winter, large surface heat loss and intense winds destroy stratification and enhance turbulent vertical mixing in the upper ocean. These phenomena intensify VEBF by promoting its components induced by the frontogenesis and TTW balance and by triggering mixed layer instability (MLI). In summer, strong stratification associated with suppressed turbulent vertical mixing weakens the contributions of the frontogenesis and TTW balance to VEBF and shifts the dominant BCI type from the MLI to the surface Charney- and Philips-like types with greatly reduced growth rate compared with that of MLI in winter. The shallow peak of the VEBF in summer is mainly attributed to the TTW balance, whereas the BCI and frontogenesis account primarily for its deep peak.

Characteristics of submesoscale eddy structures within mesoscale eddies in the Gulf of Mexico from 1/48° ECCO estimates

Submesoscale oceanic structures (<10-20 km) such as eddies and fronts are often difficult to describe given the influence of the mesoscale. In order to characterize the surface signatures of submesoscale structures, we utilize a custom spatial filtering function to separate the meso- and large-scale sea surface height (SSH) signal from the small scale SSH signal of 1/48° high resolution estimates provided by NASA’s Estimating the Circulation and Climate of the Oceans (ECCO) project. In this study, we use ECCO estimates from a 14-month global simulation between September 2011 and November 2012 with a 2 km horizontal grid spacing in the Gulf of Mexico. We then use an eddy detection and tracking algorithm to identify persistent circular features on both scales, giving rise to an atlas of submesoscale eddy-like variabilities (SEVs). We briefly investigate the geographic and temporal variability of SEVs as a whole before collocating SEVs inside mesoscale eddies, allowing us to evaluate the characteristics of internal SEVs and the impact of SEVs on mesoscale eddies. We find that SEVs, both anticyclonic and cyclonic, are ubiquitous inside mesoscale eddies with lifetimes longer than a week, accounting for an average of 10-20% of the spatial area and eddy kinetic energy of mesoscale eddies. We also show that internal SEVs are persistently associated with temperature and salinity anomalies in both eddy centers and edges of up to 0.1 °C and 0.05 psu, with anticyclonic internal SEVs being warmer and fresher while cyclonic internal SEVs are colder and saltier. Finally, we examine the life cycle of an anticyclonic Loop Current eddy, demonstrating that the number and intensity of internal SEVs within increases as the eddy approaches separation from the Loop Current until a maximum is obtained just after separation. In light of forthcoming submesoscale SSH observations from NASA’s Surface Water and Ocean Topography (SWOT) mission, our results showcase the variability of submesoscale eddy structures and their possible implications for biogeochemical cycling, the inverse energy cascade, and Loop Current prediction techniques.

An advanced towed CTD chain for physical-biological high resolution in situ upper ocean measurements

Submesoscale eddies, fronts, and filaments are ubiquitous in the upper ocean and play an important role in biogeochemical and mixing processes as well as in the energy budget. To capture the high spatial variability of submesoscale processes, it is desirable to simultaneously resolve the vertical and horizontal gradients of hydrographic properties on scales of 10 m to 10 km. We present a revised towed CTD chain, for rapid quasi-synoptic in situ measurements of submesoscale oceanographic features, that is lighter, more robust and scientifically more useful than previous towed CTD chains. This new instrument provides a horizontal resolution of O(1 m) and can be towed at speeds of up to 5 ms-1 for measurements of the upper 100 m of the water column while providing a reasonable vertical resolution of O(1 m – 10 m). Individual CTD probes are equipped with temperature, conductivity, pressure and either rapid response dissolved oxygen or fluorescence sensors at multiple depths, enabling both hydrographic and biogeochemical studies at high resolution. A flexible probe hardware allows either real-time data collection or internal data logging for offline post-processing. Finally, we outline the necessary post-processing steps and provide data examples. With the presented data examples we show and conclude that the advanced towed CTD chain is a flexible and lightweight take on the towed CTD chain concept. It can easily be adapted to scientific needs and provides high quality very high resolution oceanographic data.

Detection and characterization of submesoscale eddies from optical images: a case study in the Argentine continental shelf

ABSTRACT Eddies are structures of recognized physical and biogeochemical relevance. Because of their small space and time scales, submesoscale eddy detection is a challenge, and detection methods for these structures are scarce. In this study, we detect submesoscale eddies (SME) by means of a method implemented on satellite chlorophyll-a (CSAT) and sea surface temperature (SST) data from Visible Infrared Imaging Radiometer Suite (VIIRS) and Moderate Resolution Imaging Spectroradiometer (MODIS) images. The performance of this method is shown through the characterization of the shape and size of the detected eddies. The trajectory and characteristics of one particular eddy could be observed for 6 days over the Argentine continental shelf. The eddy reached the highest CSAT concentration (1.52 mg m−3, decreasing to less than 0.8 mg m−3 in about 3.5 km), and the lowest SST on the same day it developed its largest size. All SME detected within the study period are in agreement with cyclonic eddies and have a radius smaller than the baroclinic Rossby radius of deformation. This paper presents a workflow for the automated recognition of SME that could be applied to other regions with similar characteristics, or different ones with previous changes in the parameters.

Numerical simulation of upper ocean responses to the passage of a submesoscale eddy

In this study, a large eddy simulation model (LES) is used to investigate the effect of an idealized westward propagating submesoscale eddy on ocean surface boundary layer turbulence. Large scale forcing (LSF) terms are introduced to the model to represent the effects of submesoscale eddies using the scale separation approach.Although re-stratification is observed at the arrival of the eddy center, consistent with previous studies, strong mixing and deepening is detected before/after the arrival of the eddy center. LES experiments are conducted to investigate the mechanism behind these dynamical responses, and explore the relative importance of the LSF terms. The interaction among these forcing terms is shown to be highly nonlinear, and the combined effects of buoyancy and momentum eddy forcing dominate the boundary layer response to the submesoscale eddy. While buoyancy flux is responsible for the enhanced mixing in the boundary layer, our analysis suggests that the momentum eddy forcing can significantly alter the timing of the enhanced mixing and boost its strength through generating strong mean current in the water column that accelerates westward density advection generated by the buoyancy eddy forcing. The lighter/heavier fluid intrusion brought by the mean flow leads to strong upwelling/downwelling, enhanced mixing and deepening of the mixed layer before the arrival of the eddy. Simulations using a general circulation model do not show similar response in the upper ocean because the K profile parameterization used in the model is a function of vertical shear only, and cannot represent the enhanced turbulence due to lateral mixing brought by the submesoscale eddy.Eddy distortion due to its asymmetric decay with time and northward Ekman transport generated by the wind stress also leads to substantial differences in boundary layer responses at different locations relative to the eddy center.

Towards the Upper-Ocean Unbalanced Submesoscale Motions in the Oleander Observations

Abstract Oceanic submesoscale dynamics with horizontal scales < 20 km have similar temporal and spatial scales as internal gravity waves (IGWs), but they differ dynamically and have distinct impacts on the ocean. Separating unbalanced submesoscale motions (USMs), quasi-balanced submesoscale motions (QBMs), and IGWs in observations remains a great challenge. Based on the wave–vortex decomposition and the vertical scale separation approach for distinguishing IGWs and USMs, the long-term repeat Oleander observations in the Gulf Stream region provide an opportunity to quantify these processes separately. Here in this study, the role of USMs in the divergence is emphasized, which has confounded the wave–vortex decomposition of wintertime data in previous analyses. We also adopt the vertical filtering approach to identify the USMs by applying a high-pass filter to the vertical scales, as USMs are characterized by smaller vertical scales. This approach is tested with submesoscale-permitting model data to confirm its effectiveness in filtering the submesoscale velocity perturbations, before being applied to the compiled velocity data of the Oleander dataset (years 2005–18). The results show that the averaged submesoscale eddy kinetic energy by USMs can reach ∼1 × 10−3 m2 s−2 at z = −30 m in winter, much stronger than found in other seasons. Importantly, this study exemplifies the possibility of obtaining USMs from existing ADCP observations and reveals the seasonal dynamical regimes for the submesoscales.

Toward Earth system modeling with resolved clouds and ocean submesoscales on heterogeneous many-core HPCs.

With the aid of the newly developed 'Sunway' heterogeneous-architecture supercomputer, which has world-leading HPC (high-performance computer) capability, a series of high-resolution coupled Earth system models (SW-HRESMs) with up to 5km of atmosphere and 3km of ocean have been developed. These models can meet the needs of multiscale interaction studies with different computational costs. Here we describe the progress of SW-HRESMs development, with an overview of the major advancements made by the international Earth science community in HR-ESMs. We also show the preliminary results of SW-HRESMs with regard to capturing major weather-climate extremes in the atmosphere and ocean, stressing the importance of permitted clouds and ocean submesoscale eddies in modeling tropical cyclones and eddy-mean flow interactions, and paving the way for further model development to resolve finer scales with even higher resolution and more realistic physics. Finally, in addition to increasing model resolution, the development procedure for a non-hydrostatic cloud and ocean submesoscale resolved ESM is discussed, laying out the major scientific directions of such a huge modeling advancement.

Submesoscale inverse energy cascade enhances Southern Ocean eddy heat transport

Oceanic eddy-induced meridional heat transport (EHT) is an important process in the Southern Ocean heat budget, the variability of which significantly modulates global meridional overturning circulation (MOC) and Antarctic sea-ice extent. Although it is recognized that mesoscale eddies with scales of ~40–300 km greatly contribute to the EHT, the role of submesoscale eddies with scales of ~1–40 km remains unclear. Here, using two state-of-the-art high-resolution simulations (resolutions of 1/48° and 1/24°), we find that submesoscale eddies significantly enhance the total poleward EHT in the Southern Ocean with an enhancement percentage reaching 19–48% in the Antarctic Circumpolar Current band. By comparing the eddy energy budgets between the two simulations, we detect that the primary role of submesoscale eddies is to strengthen mesoscale eddies (and thus their heat transport capability) through inverse energy cascade rather than directly through submesoscale heat fluxes. Due to the submesoscale-mediated enhancement of mesoscale eddies in the 1/48° simulation, the clockwise upper cell and anti-clockwise lower cell of the residual-mean MOC in the Southern Ocean are weakened and strengthened, respectively. This finding identifies a potential route to improve the mesoscale parameterization in climate models for more accurate simulations of the MOC and sea ice variability in the Southern Ocean.

Equatorial Submesoscale Eddies Contribute to the Asymmetry in ENSO Amplitude

AbstractEl Niño/Southern Oscillation (ENSO) exhibits evident amplitude asymmetry with stronger El Niño than La Niña events. Equatorial submesoscale eddies act as an important damper of ENSO via their induced upward heat flux () from the subsurface to the surface ocean. Yet their effect on ENSO amplitude asymmetry remains unexplored. Using a high‐resolution global climate simulation, we found that strengthening of in the Niño3.4 region during La Niña is more evident than weakening of during El Niño, resulting in stronger damping of La Niña than El Niño events. This asymmetric response of to ENSO is primarily ascribed to that of the temperature fronts associated with larger‐scale processes. Using a recharge oscillator model modified to include the effects of , we show that the asymmetric damping of ENSO by contributes importantly to the ENSO amplitude asymmetry.

ICON-Sapphire: simulating the components of the Earth system and their interactions at kilometer and subkilometer scales

Abstract. State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

GENERATION OF INTERNAL WAVES BY A SUB-MESOSCALE EDDY AT CAPE SVYATOY NOS IN THE BARENTS SEA

Data on generation of soliton-like internal waves by a mesoscale eddy in the region of the Barents Sea near Cape Svyatoy Nos, widely known for its dynamic activity, are presented here. The results are based on the analysis of satellite optical (Sentinel-2 MSI) and radar (Sentinel-1 SAR) images obtained sequentially on August 8 and 9, 2021. A detailed analysis of the dynamics of the eddy and internal waves was made possible due, among other things, to the phytoplankton bloom, which created the necessary optical contrasts that make it possible to monitor the dynamics of currents in the near-surface layer of the sea. Significant parameters of recorded internal waves were measured using satellite methods and analyzed using the TPXO9.0 model and numerical solution of the internal wave equation. It has been shown, that internal waves were caused by a submesoscale eddy formed near the cape. Waves had quite low velocity (0.10–0.16 m/s) what was due to weak stratification of the water environment.

Possible Mechanism for the Generation of Internal Waves by Submesoscale Eddies in the Mozambique Channel

Possible Mechanism for the Generation of Internal Waves by Submesoscale Eddies in the Mozambique Channel

Enhancement of Oceanic Eddy Activity by Fine‐Scale Orographic Winds Drives High Productivity, Low Oxygen, and Low pH Conditions in the Santa Barbara Channel

AbstractThe Santa Barbara Channel is one of the most productive regions of the California Current System. Yet, the physical processes that sustain this high productivity remain unclear. We use a high‐resolution physical‐biogeochemical model to show that submesoscale eddies generated by islands are energized by orographic effects on the wind, with significant impacts on nutrient, carbon, and oxygen cycles. These eddies are modulated by two co‐occurring air‐sea‐land interactions: transfer of wind energy to ocean currents that intensifies ocean eddies, and a wind‐current feedback that tends to dampen them. Here we show that the dampening is overwhelmed by fine scale wind patterns induced by the presence of surrounding capes and islands. The fine‐scale winds cause an additional transfer of momentum from the atmosphere to the ocean that energizes submesoscale eddies. This drives upward doming of isopycnals in the center of the channel, allowing a more efficient injection of nutrients to the surface, and triggering intense phytoplankton blooms that nearly double productivity relative to the case without fine‐scale winds. The intensification of the doming effect by the wind‐curl and submesoscale eddies pumps deep low oxygen, acidic waters to the center of the cyclonic eddies. These eddies are then transported away from the Channel into the California Current, where they impact a wider area along the central coast, with potential ecological consequences. Our study highlights the important role of air‐sea‐land interactions in modulating coastal processes, and suggests that submesoscale resolving models are required to correctly represent coastal processes and their ecological impacts.

Submesoscale eddies in eastern Guangdong identified using high-frequency radar observations

The statistical characteristics and mechanism underlying the generation of submesoscale eddies in eastern Guangdong were examined using high-frequency radar (HFR) observations of surface currents during 2019_2021. By applying the vector geometry-based (VG) eddy detection algorithm to HFR datasets, 313 submesoscale eddies with 229 cyclonic and 84 anticyclonic eddies were identified. Cyclonic eddies were significantly more common than the anticyclonic eddies, but they had similar distributions in terms of eddy lifespan and radius. The trajectories showed that the submesoscale eddies within the study area tended to move northward and westward. Submesoscale eddies revealed annual variations with decreasing cyclonic eddies and increasing anticyclonic eddies during the observation period. Submesoscale eddies displayed a remarkable spatial distribution. Most of cyclonic (anticyclonic) eddies occurred at the eastern and southeastern regions of the Huilai Headland (south of the Nanao Island). Eddy evolution was linked to the production of vorticity in a narrow layer along the headland, the injection of vorticity into the interior at the point of flow separation, and vorticity damping by bottom friction.

Internal tides off the Amazon shelf during two contrasted seasons: interactions with background circulation and SSH imprints

Abstract. The Amazon shelf break is a key region for internal tide (IT) generation. It also shows a large seasonal variation in circulation and associated stratification. This study, based on a high-resolution model (1/36∘) explicitly forced by tide, aims to better characterize how the ITs vary between two contrasted seasons. During the season from March to July (MAMJJ) the currents and mesoscale eddies are weak while the pycnocline is shallower and stronger. From August to December (ASOND) mean currents and mesoscale eddies are strong, and the pycnocline is deeper and weaker than in MAMJJ. For both seasons, semi-diurnal M2 ITs are generated on the shelf break mainly between the 100 and 1000 m isobath in the model. South of 2∘ N, the conversion from barotropic to baroclinic tide is more efficient in MAMJJ than in ASOND. Local dissipation of the coherent M2 at the generation sites is higher in MAMJJ (30 %) than in ASOND (22 %), because higher modes are favourably generated (modes 2 and 3), making the internal wave packet more dispersive. The remaining fraction (70 %–80 %) propagates away from the generation sites and mainly dissipates locally every ∼ 100 km, which corresponds to the mode 1 reflection beams. About 13 %, 30 %, and 40 % of the M2 coherent IT dissipates at the first, second, and third beams. M2 coherent baroclinic flux propagates more northward during MAMJJ while it seems to be blocked at 6∘ N during ASOND. There is no intensified dissipation of the coherent M2 that could explain the disappearance of the coherent flux. In fact, the flux at this location becomes more incoherent because of strong interaction with the currents. This has been shown in the paper using 25 h mean snapshots of the baroclinic flux that shows branching and stronger eastward deviation of the IT when interacting with mesoscale eddies and stratification during ASOND. Finally, we evaluated sea surface height (SSH) frequency and wavenumber spectra for subtidal (f<1/28 h−1), tidal (1/28 h−1 < f<1/11 h−1), and supertidal (f>1/11 h−1) frequencies. Tidal frequencies explain most of the SSH variability for wavelengths between 250 and 70 km. Below 70 km, the SSH is mainly incoherent and supertidal. The length scale at which the SSH becomes dominated by unbalanced (non-geostrophic) IT was estimated to be around 250 km. Our results highlight the complexity of correctly predicting IT SSH in order to better observe mesoscale and sub-mesoscale eddies from existing and upcoming altimetric missions, notably the Surface Water Ocean Topography (SWOT) mission.

Weakened Submesoscale Eddies in the Equatorial Pacific Under Greenhouse Warming

AbstractSubmesoscale eddies in the equatorial Pacific induce a pronounced upward heat flux from the subsurface to the surface ocean. Yet their response to greenhouse warming and feedback on the mean‐state change remain unknown. Using a long‐term high‐resolution climate simulation under a high carbon emission scenario, we show that submesoscale eddies in the central and eastern equatorial Pacific weaken, leading to an ∼50% reduction in their upward heat flux at the end of the 21st century, compared to its mid‐20th century level. This reduction is caused by an intensified stratification and a weakened cold tongue thermal front associated with a projected El Niño‐like warming pattern. The reduced vertical submesoscale eddy heat flux plays an important role in the anomalous heat budget of the surface layer over the central and eastern equatorial Pacific under greenhouse warming, acting as negative feedback on the El Niño‐like warming pattern.

Topological Feature Tracking for Submesoscale Eddies

AbstractCurrent state‐of‐the art procedures for studying modeled submesoscale oceanographic features have made a strong assumption of independence between features identified at different times. Therefore, all submesoscale eddies identified in a time series were studied in aggregate. Statistics from these methods are illuminating but oversample identified features and cannot determine the lifetime evolution of the transient submesoscale processes. To this end, the authors apply the Topological Feature Tracking (TFT) algorithm to the problem of identifying and tracking submesoscale eddies over time. TFT identifies critical points on a set of time‐ordered scalar fields and associates those points between consecutive timesteps. The procedure yields tracklets which represent spatio‐temporal displacement of eddies. In this way we study the time‐dependent behavior of submesoscale eddies, which are generated by a 1‐km resolution submesoscale‐permitting model. We summarize the submesoscale eddy data set produced by TFT, which yields unique, time‐varying statistics.