Oceanic mesoscale eddies with typical sizes of 30–200 km contain more than half of the kinetic energy of the ocean. With an average lifespan of several months, they are major contributors to the transport of heat, nutrients, plankton, dissolved oxygen and carbon in the ocean. Mesoscale eddies have been observed and studied over the past 50 years, nonetheless our understanding of the details of their structure remains incomplete due to lack of systematic high-resolution measurements. To bridge this gap, a survey of a mesoscale anticyclone was conducted in early 2014 in the South China Sea capturing its structure at submesoscale resolution. By modeling an anticyclone of comparable size and position at three horizontal resolutions the authors verify the resolution requirements for capturing the observed variability in dynamical quantities, and quantify the role of ageostrophic motions on the vertical transport associated with the anticyclone. Results indicate that different submesoscale processes contribute to the vertical transport depending on depth and distance from the eddy center, with frontogenesis playing a key role. Vertical transport by anticyclones cannot be reliably estimated by coarse-resolution or even mesoscale-resolving models, with important implications for global estimates of the eddy-driven vertical pumping of biophysical and chemical tracers.

Direct microstructure observations across three warm mesoscale eddies were conducted in the northern South China Sea during the field experiments in July 2007, December 2013, and January 2014, respectively, along with finestructure measurements. An important finding was that turbulent mixing in the mixed layer was considerably elevated in the periphery of each of these eddies, with a mixing level 5–7 times higher than that in the eddy center. To explore the mechanism behind the high mixing level, this study carried out analyses of the horizontal wavenumber spectrum of velocities and spectral fluxes of kinetic energy. Spectral slopes showed a power law of k−2 in the eddy periphery and of k−3 in the eddy center, consistent with the result that the kinetic energy of submesoscale motion in the eddy periphery was more greatly energized than that in the center. Spectral fluxes of kinetic energy also revealed a forward energy cascade toward smaller scales at the wavelength of kilometers in the eddy periphery. This study illustrated a possible route for energy cascading from balanced mesoscale dynamics to unbalanced submesoscale behavior, which eventually furnished turbulent mixing in the upper ocean.

The spatial scale of submesoscales is an important parameter for studies of submesoscale dynamics and multi-scale interactions. The horizontal spatial scales of baroclinic, geostrophic-branch mixed layer instabilities (MLI) are investigated globally (without the equatorial or Arctic oceans) based on observations and simulations in the surface and bottom mixed layers away from significant topography. Three high vertical-resolution boundary layer schemes driven with profiles from a MITgcm global submesoscale-permitting model improve robustness. The fastest-growing MLI wavelength decreases toward the poles. The zonal median surface MLI wavelength is 51 km to 2.9 km when estimated from the observations and 32, 25, and 27 km to 2.5, 1.2, and 1.1 km under the K-profile Parameterization (KPP), Mellor-Yamada (MY), and k-ε schemes, respectively. The surface MLI wavelength has a strong seasonality with a median value 1.6 times smaller in summer (10 km) than winter (16 km) globally from the observations. The median bottom MLI wavelengths estimated from simulations are 2.1, 1.4, and 0.41 km globally under the KPP, MY and k-ε schemes with little seasonality. The estimated required ocean model grid spacings to resolve wintertime surface mixed layer eddies are 1.9 km (50% of regions resolved) and 0.92 km (90%) globally. To resolve summertime eddies or MLI seasonality requires grids finer than 1.3 km (50%) and 0.55 km (90%). To resolve bottom mixed layer eddies grids finer than 257, 178, and 51 m (50%) and 107, 87, and 17 m (90%) are estimated under the KPP, MY, and k-ε schemes.

Submesoscale processes in the upper ocean vary seasonally, in tight correspondence with mixed layer thickness variability. Based on a global high‐resolution MITgcm simulation, seasonal evaluation of strong vorticity and spectral analysis of the kinetic energy in the Kuroshio Extension System show the strongest submesoscales occur in March, implying a lag of about a month behind mixed layer thickness maximum in February. An analysis of spectral energy sources and transfers indicates that the seasonality of the submesoscale energy content is a result of the competition between the conversion of available potential energy into submesoscale kinetic energy via a buoyancy production/vertical buoyancy flux associated with mixed layer instability and nonlinear energy transfers to other scales associated with an energy cascade. The buoyancy production is seasonally in phase with the mixed layer depth, but the transfers of energy across scales makes energizing the reservoir of submesoscale kinetic energy lag behind by a month.

Symmetric instability (SI) extracts kinetic energy from fronts in the surface mixed layer (SML), potentially affecting the SML structure and dynamics. Here, a global submesoscale-permitting ocean model named MITgcm LLC4320 simulation is used to examine the Stone (1966) linear prediction of the maximum SI scale to estimate grid spacings needed to begin resolving SI. Furthermore, potential effects of SI on the usable wind-work are estimated roughly: this estimate of SI “activity” is useful for assessing if these modes should be resolved or parameterized. The maximum SI scale varies by latitude with median values of 568 m to 23 m. Strong seasonality is observed in the SI scale and activity. The median scale in winter is 188 m globally, 2.5 times of that of summer (75 m). SI is more active in winter: 15% of the time compared with 6% in summer. The strongest SI activity is found in the western Pacific, western Atlantic, and Southern Oceans. The required grid spacings for a global model to begin resolving SI eddies in the SML are 24 m (50% of regions resolved) and 7.9 m (90%) in winter, decreasing to 9.4 m (50%) and 3.6 m (90%) in summer. It is also estimated that SI may reduce usable wind-work by an upper bound of 0.83 mW m−2 globally, or 5% of the global magnitude. The sensitivity of these estimates to empirical thresholds is provided in the text.

From the analysis of oceanic eddies detected in the drifter trajectories of the Global Drifter Program (GDP) data set, it was found that oceanic eddies are asymmetrically distributed across the Kuroshio in the East China Sea: predominant cyclonic (anticyclonic) eddies are on the western (eastern) sides of Kuroshio. This distribution is confirmed by high-resolution numerical modeling output as well. Most of these eddies are 5~20 km in radius, less than the local first baroclinic deformation radius, thus categorized as submesoscale. The generation mechanism of these submesoscale eddies is speculated to be related to the horizontal velocity shear of the Kuroshio when it flows northeastward along the shelf break in the East China Sea. The budget analysis of eddy kinetic energy shows that both the horizontal shear and vertical buoyancy flux are important energy sources for eddy generation on the two sides of Kuroshio axis. The finding highlights the unique feature of oceanic eddies along the western boundary currents.

Physical processes in the ocean surface mixed layer (SML) modulate the SML, impacting the energy budget of global ocean and thereby climate change. As SML processes including surface waves, Langmuir circulation and submesoscale currents have spatial scales far from being resolved by climate ocean models (from tens of meters to kilometers), parameterizations are a key avenue to improve climate and weather models (e.g., Fox

Critical latitudes are a significant area of tidal dissipation. Generally, critical latitudes are taken to be the exact latitude where the tidal frequency equals the inertial frequency. However, the key is really where the tidal frequency equals the combination of planetary vorticity and relative vorticity from background currents. Although the influence of background currents on critical latitude effects and nonlinear interactions have been noted for many years, their exact impacts are not well known. The latitude dependence of critical latitude impacts on the tides, internal tides, and internal waves in the presence of background currents was investigated using the Regional Ocean Modeling System by shifting a small domain including a seamount from 20.6° to 38.6°S and comparing simulations with and without background currents. The diurnal kinetic energy with mesoscale currents was relatively unchanged for most latitudes, except for a slight decrease 1–4° poleward of the critical latitude. However, the semidiurnal and high‐frequency (≥3 cycles per day) kinetic energy increased with the presence of mesoscale currents, especially within the diurnal critical latitude range. Spectral and nonlinear analyses indicated mesoscale currents broadened the range of critical latitude effects and enhanced energy transferring from diurnal frequencies to semidiurnal and high frequencies and from low to high mode waves. Local diffusivities increased, roughly an order of magnitude, when mesoscale currents were present. The impacts of mesoscale currents on the broadening of the critical latitude range and enhancement of nonlinear interactions were attributed to the additional relative vorticity and near‐inertial internal waves generated by mesoscale currents.