Spectra and fluxes of enstrophy and kinetic energy (KE) are estimated in different parts of the mid‐latitudinal oceans using geostrophic currents derived from altimetry data. The presence of a strong inverse flux of surface KE is confirmed at scales larger than approximately 200 km, whereas a robust enstrophy cascading regime, accompanied by an approximate k^‐3 KE spectrum, is observed from about 200 to 100 km. The character of fluxes and spectra are shown to compare favorably with those from a comprehensive Earth system model. In addition, as gridded altimeter data is affected by smoothening and interpolation, the qualitative robustness of the results is verified by sensitivity experiments using space and time filtered output from the Earth system model. Given the rotational character of the flow, this large‐scale inverse KE and smaller‐scale forward enstrophy transfer scenario is consistent with expectations from three‐dimensional rapidly‐rotating and strongly‐stratified turbulence studies as well as detailed analyses of spectra and fluxes in the upper‐level mid‐latitude troposphere. In further accord with results from the atmosphere, decomposing the currents into stationary and eddy components (demarcated here by variability greater and less than 100 days, respectively), it is seen that, in addition to the eddy‐eddy contribution, the stationary‐eddy and stationary‐stationary fluxes play a significant role in the inverse (forward) flux of KE (enstrophy). Thus, it is quite possible that, from about 200 km to 100 km, the altimeter is capturing the rotationally dominated portion of a surface oceanic counterpart of the upper tropospheric Nastrom‐Gage spectrum.

Climate model hindcasts are analyzed to reveal the impacts of the North Atlantic Oscillation (NAO) on the North Atlantic subpolar ocean, which exhibits variability on seasonal to decadal timescales. The ocean response to a single winter NAO event is separated into fast and slow responses. The fast response persists over winter–spring seasons, during which wind stress and heat flux anomalies associated with the NAO rapidly modify ocean temperatures via changes in Ekman transport and ocean‐atmosphere heat exchanges. The slow response persists for 3–4 years, during which overturning and gyre circulations redistribute opposing‐signed surface temperature anomalies created by the NAO. This redistribution modifies east‐west temperature contrasts altering the meridional heat transport associated with gyres and changing the strength of the overturning circulation. Hence, the fast and slow responses lead to opposing‐signed subpolar temperature anomalies in time from the competing effects of local forcing and horizontal heat convergence.

Here, we present an estimate for the ocean's global scale transfer of kinetic energy (KE), across scales from 10 to 40,000 km. Oceanic KE transfer between gyre scales and mesoscales is induced by the atmosphere’s Hadley, Ferrel, and polar cells, and the intertropical convergence zone induces an intense downscale KE transfer. Upscale transfer peaks at 300 gigawatts across mesoscales of 120 km in size, roughly one-third the energy input by winds into the oceanic general circulation. Nearly three quarters of this “cascade” occurs south of 15°S and penetrates almost the entire water column. The mesoscale cascade has a self-similar seasonal cycle with characteristic lag time of ≈27 days per octave of length scales; transfer across 50 km peaks in spring, while transfer across 500 km peaks in summer. KE of those mesoscales follows the same cycle but peaks ≈40 days after the peak cascade, suggesting that energy transferred across a scale is primarily deposited at a scale four times larger.

We expand on a recent determination of the first global energy spectrum of the ocean's surface geostrophic circulation (Storer et al., 2022, https://doi.org/10.1038/s41467-022-33031-3) using a coarse‐graining (CG) method. We compare spectra from CG to those from spherical harmonics by treating land in a manner consistent with the boundary conditions. While the two methods yield qualitatively consistent domain‐averaged results, spherical harmonics spectra are too noisy at gyre‐scales (>1,000 km). More importantly, spherical harmonics are inherently global and cannot provide local information connecting scales with currents geographically. CG shows that the extra‐tropics mesoscales (100–500 km) have a root‐mean‐square (rms) velocity of ∼15 cm/s, which increases to ∼30–40 cm/s locally in the Gulf Stream and Kuroshio and to ∼16–28 cm/s in the ACC. There is notable hemispheric asymmetry in mesoscale energy‐per‐area, which is higher in the north due to continental boundaries. We estimate that ≈25%–50% of total geostrophic energy is at scales smaller than 100 km, and is un(der)‐resolved by pre‐SWOT satellite products. Spectra of the time‐mean circulation show that most of its energy (up to 70%) resides in stationary eddies with characteristic scales smaller than (<500 km). This highlights the preponderance of “standing” small‐scale structures in the global ocean due to the temporally coherent forcing by boundaries. By coarse‐graining in space and time, we compute the first spatio‐temporal global spectrum of geostrophic circulation from AVISO and NEMO. These spectra show that every length‐scale evolves over a wide range of time‐scales with a consistent peak at ≈200 km and ≈2–3 weeks.

We describe an idealized primitive-equation model for studying mesoscale turbulence and leverage a hierarchy of grid resolutions to make eddy-resolving calculations on the finest grids more affordable. The model has intermediate complexity, incorporating basin-scale geometry with idealized Atlantic and Southern oceans and with non-uniform ocean depth to allow for mesoscale eddy interactions with topography. The model is perfectly adiabatic and spans the Equator and thus fills a gap between quasi-geostrophic models, which cannot span two hemispheres, and idealized general circulation models, which generally include diabatic processes and buoyancy forcing. We show that the model solution is approaching convergence in mean kinetic energy for the ocean mesoscale processes of interest and has a rich range of dynamics with circulation features that emerge only due to resolving mesoscale turbulence.