• Alford, M. H., J. A. MacKinnon, H. L. Simmons, and J. D. Nash, 2016: Near-inertial internal gravity waves in the ocean. Annu. Rev. Mar. Sci., 8, 95123, https://doi.org/10.1146/annurev-marine-010814-015746.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Buckingham, C. E., and Coauthors, 2016: Seasonality of submesoscale flows in the ocean surface boundary layer. Geophys. Res. Lett., 43, 21182126, https://doi.org/10.1002/2016GL068009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bühler, O., J. Callies, and R. Ferrari, 2014: Wave–vortex decomposition of one-dimensional ship-track data. J. Fluid Mech., 756, 10071026, https://doi.org/10.1017/jfm.2014.488.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., and R. Ferrari, 2013: Interpreting energy and tracer spectra of upper ocean turbulence in the submesoscale range (1–200 km). J. Phys. Oceanogr., 43, 24562474, https://doi.org/10.1175/JPO-D-13-063.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., R. Ferrari, J. M. Klymak, and J. Gula, 2015: Seasonality in submesoscale turbulence. Nat. Commun., 6, 6862, https://doi.org/10.1038/ncomms7862.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Capet, X., J. C. McWilliams, M. J. Molemaker, and A. J. Shchepetkin, 2008: Mesoscale to submesoscale transition in the California Current System. Part I: Flow structure, eddy flux, and observational tests. J. Phys. Oceanogr., 38, 2943, https://doi.org/10.1175/2007JPO3671.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chelton, D. B., M. G. Schlax, and R. M. Samelson, 2011: Global observations of nonlinear mesoscale eddies. Prog. Oceanogr., 91, 167216, https://doi.org/10.1016/j.pocean.2011.01.002.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chen, X., B. Qiu, Y. Dan, S. Chen, and Y. Qi, 2016: Interannual and interdecadal variability of the North Equatorial Countercurrent in the western Pacific. J. Geophys. Res. Oceans, 121, 77437758, https://doi.org/10.1002/2016JC012190.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • D’Asaro, E. A., 1978: Mixed layer velocities induced by internal wave. J. Geophys. Res., 83, 24372438, https://doi.org/10.1029/JC083iC05p02437.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • de Boyer Montégut, C., G. Madec, A. S. Fischer, A. Lazar, and D. Iudicone, 2004: Mixed layer depth over the global ocean: An examination of profile data and a profile-based climatology. J. Geophys. Res., 109, C12003, https://doi.org/10.1029/2004JC002378.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ducet, N., P.-Y. Le Traon, and G. Reverdin, 2000: Global high-resolution mapping of ocean circulation from TOPEX/Poseidon and ERS-1 and -2. J. Geophys. Res., 105, 19 47719 498, https://doi.org/10.1029/2000JC900063.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ferrari, R., and C. Wunsch, 2009: Ocean circulation kinetic energy: Reservoirs, sources, and sinks. Annu. Rev. Fluid Mech., 41, 253282, https://doi.org/10.1146/annurev.fluid.40.111406.102139.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fu, L.-L., and C. Ubelmann, 2014: On the transition from profile altimetry to swath altimeter for observing global ocean surface topography. J. Atmos. Oceanic Technol., 31, 560568, https://doi.org/10.1175/JTECH-D-13-00109.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fu, L.-L., D. B. Chelton, P.-Y. Le Traon, and R. Morrow, 2010: Eddy dynamics from satellite altimetry. Oceanography, 23 (4), 1425, https://doi.org/10.5670/oceanog.2010.02.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gill, A. E., 1982: Atmosphere–Ocean Dynamics, Academic Press, 662 pp.

  • Klein, P., B. L. Hua, G. Lapeyre, X. Capet, S. Le Gentil, and H. Sasaki, 2008: Upper ocean turbulence from high-resolution 3D simulations. J. Phys. Oceanogr., 38, 17481763, https://doi.org/10.1175/2007JPO3773.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lapeyre, G., and P. Klein, 2006: Dynamics of the upper oceanic layers in terms of surface quasigeostrophy theory. J. Phys. Oceanogr., 36, 165176, https://doi.org/10.1175/JPO2840.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lapeyre, G., P. Klein, and B. L. Hua, 2006: Oceanic restratification forced by surface frontogenesis. J. Phys. Oceanogr., 36, 15771590, https://doi.org/10.1175/JPO2923.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Locarnini, R. A., and Coauthors, 2013: Temperature. Vol. 1, World Ocean Atlas 2013, NOAA Atlas NESDIS 73, 40 pp.

  • Marshall, J. C., A. Adcroft, C. Hill, L. Perelman, and C. Heisey, 1997: A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers. J. Geophys. Res., 102, 57535766, https://doi.org/10.1029/96JC02775.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., 2016: Submesoscale currents in the ocean. Proc. Roy. Soc., A472, 20160177, https://doi.org/10.1098/rspa.2016.0117.

  • Menemenlis, D., J.-M. Campin, P. Heimbach, C. Hill, T. Lee, A. Nguyen, M. Schodlok, and H. Zhang, 2008: ECCO2: High-resolution global ocean and sea ice data synthesis. Mercator Ocean Quarterly Newsletter, No. 31, Mercator Ocean, Ramonville Saint-Agne, France, 13–21.

  • Mensa, J. A., Z. Garraffo, A. Griffa, T. M. Ozgokmen, A. Haza, and M. Veneziani, 2013: Seasonality of the submesoscale dynamics in the Gulf Stream region. Ocean Dyn., 63, 923941, https://doi.org/10.1007/s10236-013-0633-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Morrow, R., and P.-Y. Le Traon, 2012: Recent advances in observing mesoscale ocean dynamics with satellite altimetry. Adv. Space Res., 50, 10621076, https://doi.org/10.1016/j.asr.2011.09.033.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., and T. M. Joyce, 1992: Interannual variability in the mid- and low-latitude western North Pacific. J. Phys. Oceanogr., 22, 10621079, https://doi.org/10.1175/1520-0485(1992)022<1062:IVITMA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., S. Chen, P. Klein, H. Sasaki, and Y. Sasai, 2014: Seasonal mesoscale and submesoscale eddy variability along the North Pacific Subtropical Countercurrent. J. Phys. Oceanogr., 44, 30793098, https://doi.org/10.1175/JPO-D-14-0071.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., S. Chen, P. Klein, C. Ubelmann, L.-L. Fu, and H. Sasaki, 2016: Reconstructability of three-dimensional upper-ocean circulation from SWOT sea surface height measurements. J. Phys. Oceanogr., 46, 947963, https://doi.org/10.1175/JPO-D-15-0188.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., T. Nakano, S. Chen, and P. Klein, 2017: Submesoscale transition from geostrophic flows to internal waves in the northwestern Pacific upper ocean. Nat. Commun., 8, 14055, https://doi.org/10.1038/ncomms14055.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ray, R. D., and G. T. Mitchum, 1996: Surface manifestation of internal tides generated near Hawaii. Geophys. Res. Lett., 23, 21012104, https://doi.org/10.1029/96GL02050.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ray, R. D., and E. D. Zaron, 2011: Non-stationary internal tides observed with satellite altimetry. Geophys. Res. Lett., 38, L17609, https://doi.org/10.1029/2011GL048617.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ray, R. D., and E. D. Zaron, 2016: M2 internal tides and their observed wavenumber spectra from satellite altimetry. J. Phys. Oceanogr., 46, 322, https://doi.org/10.1175/JPO-D-15-0065.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Richman, J. G., B. K. Arbic, J. F. Shriver, E. J. Metzger, and A. J. Wallcraft, 2012: Inferring dynamics from the wavenumber spectra of an eddying global ocean model with embedded tides. J. Geophys. Res., 117, C12012, https://doi.org/10.1029/2012JC008364.

    • Search Google Scholar
    • Export Citation
  • Rio, M.-H., S. Guinehut, and G. Larnico, 2011: New CNES-CLS09 global mean dynamic topography computed from the combination of GRACE data, altimetry, and in situ measurements. J. Geophys. Res., 116, C07018, https://doi.org/10.1029/2010JC006505.

    • Search Google Scholar
    • Export Citation
  • Rocha, C. B., T. K. Chereskin, S. T. Gille, and D. Menemenlis, 2016a: Mesoscale to submesoscale wavenumber spectra in Drake Passage. J. Phys. Oceanogr., 46, 601620, https://doi.org/10.1175/JPO-D-15-0087.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rocha, C. B., S. T. Gille, T. K. Chereskin, and D. Menemenlis, 2016b: Seasonality of submesoscale dynamics in the Kuroshio Extension. Geophys. Res. Lett., 43, 11 30411 311, https://doi.org/10.1002/2016GL071349.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rodriguez, E., 2016: Surface Water and Ocean Topography Mission (SWOT) project: Science requirements document. NASA JPL Rep. D-61923, 28 pp., https://swot.oceansciences.org/docs/D-61923_SRD_Rev_A_20160318.pdf.

  • Sasaki, H., P. Klein, B. Qiu, and Y. Sasai, 2014: Impact of oceanic-scale interactions on the seasonal modulation of ocean dynamics by the atmosphere. Nat. Commun., 5, 5636, https://doi.org/10.1038/ncomms6636.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Savage, A. C., and Coauthors, 2017a: Frequency content of sea surface height variability from internal gravity waves to mesoscale eddies. J. Geophys. Res. Oceans, 122, 25192538, https://doi.org/10.1002/2016JC012331.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Savage, A. C., and Coauthors, 2017b: Spectral decomposition of internal gravity wave sea surface height in global models. J. Geophys. Res. Oceans, 122, 78037821, https://doi.org/10.1002/2017JC013009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Shriver, J. F., J. G. Richman, and B. K. Arbic, 2014: How stationary are the internal tides in a high-resolution global ocean circulation model? J. Geophys. Res. Oceans, 119, 27692787, https://doi.org/10.1002/2013JC009423.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stammer, D., 1997: Global characteristics of ocean variability estimated from regional TOPEX/Poseidon altimeter measurements. J. Phys. Oceanogr., 27, 17431769, https://doi.org/10.1175/1520-0485(1997)027<1743:GCOOVE>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., A. Tandon, and A. Mahadevan, 2008: Submesoscale processes and dynamics. Ocean Modeling in an Eddying Regime, Geophys. Monogr., Vol. 177, Amer. Geophys. Union, 17–38.

    • Crossref
    • Export Citation
  • Thompson, A. F., A. Lazar, C. Buckingham, A. C. Naveira Garabato, G. M. Damerell, and K. J. Heywood, 2016: Open-ocean submesoscale motions: A full seasonal cycle of mixed layer instabilities from gliders. J. Phys. Oceanogr., 46, 12851307, https://doi.org/10.1175/JPO-D-15-0170.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wang, J., L.-L. Fu, B. Qiu, D. Menemenlis, T. Farrar, Y. Chao, A. Thompson, and M. Flexas, 2018: An observing system simulation experiment for the calibration and validation of the Surface Water and Ocean Topography sea surface height measurement using in-situ platforms. J. Atmos. Oceanic Technol., 35, 281297, https://doi.org/10.1175/JTECH-D-17-0076.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zaron, E. D., 2017: Mapping the nonstationary internal tide with satellite altimetry. J. Geophys. Res. Oceans, 122, 539554, https://doi.org/10.1002/2016JC012487.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zhao, Z., M. H. Alford, J. B. Girton, L. Rainville, and H. L. Simmons, 2016: Global observations of open-ocean mode-1 M2 internal tides. J. Phys. Oceanogr., 46, 16571684, https://doi.org/10.1175/JPO-D-15-0105.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zweng, M. M., and Coauthors, 2013: Salinity. Vol. 2, World Ocean Atlas 2013, NOAA Atlas NESDIS 74, 39 pp.

    • Crossref
    • Export Citation
  • View in gallery
    Fig. 1.

    Surface EKE distributions from (a) MITgcm llc4320 simulation and (b) AVISO SSH product of November 2011–October 2012. For both maps, EKE is calculated geostrophically from the weekly, 1/4° × 1/4°–resolution SSH field.

  • View in gallery
    Fig. 2.

    Horizontal wavenumber–frequency distributions of surface kinetic energy in the 500 km × 500 km box centered on (a) 37°N, 164°E, (b) 21°N, 164°E, (c) 9°N, 164°E, and (d) 7°S, 164°E of the Pacific Ocean. Based on hourly MITgcm output of 1 Oct 2011–30 Sep 2012. Dashed white lines denote the inertial and tidal frequencies. Solid white lines denote the dispersion relation for inertia–gravity waves of the first 10 vertical modes. Solid black lines denote the lower-frequency boundary of either the IGW dispersion curve at the equatorward edge of the box or the local permissible tides. Dashed black lines denote the frequency of the 36-h filter. Tidal constituent S2, M2, K1, O1, and inertial frequencies are denoted by dashed white lines.

  • View in gallery
    Fig. 3.

    Horizontal wavenumber spectra of the surface eddy kinetic energy of balanced (thick black lines) and unbalanced (thin black lines) motions in the 500 km × 500 km box centered on (a) 37°N, 164°E, (b) 21°N, 164°E, (c) 9°N, 164°E, and (d) 7°S, 164°E of the Pacific Ocean. Black dashed lines denote the reference spectra and vertical gray dashed lines denote the transition scale Lt.

  • View in gallery
    Fig. 4.

    Distribution of the transition scale Lt indicating the dominance from balanced to unbalanced motions. Areas with Lt > 500 km are blanked out. White contours denote the 3000-m water depth.

  • View in gallery
    Fig. 5.

    Distributions of wavenumber-integrated EKE for (a) balanced and (b) unbalanced motions based on the dynamical delineation defined in Eqs. (2) and (3). Note that color scales for the EKE level are different for (a) and (b). White contours denote the 3000-m water depth.

  • View in gallery
    Fig. 6.

    Distributions of seasonal transition scale Lt in (a) February–April and (b) August–October. Areas with Lt > 500 km are blanked out.

  • View in gallery
    Fig. 7.

    Distribution of seasonal difference (ASO − FMA) in Lt or Fig. 6b minus Fig. 6a. Gray areas denote boxes in which either ASO or FMA Lt exceeds 500 km.

  • View in gallery
    Fig. 8.

    Distributions of seasonal difference (ASO − FMA) in wavenumber-integrated EKE for (a) balanced and (b) unbalanced motions with length scales shorter than 100 km.

  • View in gallery
    Fig. 9.

    Distributions of seasonal mixed layer depth in (a) FMA and (b) ASO based on the llc4320 simulation.

  • View in gallery
    Fig. 10.

    Distributions of seasonal difference (ASO − FMA) in wavenumber-integrated EKE at 500 m for (a) balanced and (b) unbalanced motions with length scales shorter than 100 km. Notice that the color scale is 10 times smaller in this figure than in Fig. 8.

  • View in gallery
    Fig. 11.

    Cumulative PDF distribution for wintertime Lt values (FMA for NH and ASO for SH) estimated in individual latitudinal bands of the World Ocean. Black solid lines denote the 80% threshold Lt values.

  • View in gallery
    Fig. 12.

    (a) Distribution of the transition scale based on equal SSH variance between the balanced and unbalanced motions. Areas with Lt > 500 km are blanked out. (b) Distribution of Lt. Gray areas denote boxes in which either or Lt exceeds 500 km.

  • View in gallery
    Fig. 13.

    Threshold lines of I = 1 for discrete vertical mode (n = 1 to 10) IGWs as a function of kh [see Eq. (14)]. Values are averaged across the latitudinal bands of the World Ocean. For a fixed vertical mode, I < 1 in the low-latitude and high-wavenumber space of the threshold line.

  • View in gallery
    Fig. 14.

    Distributions of seasonal difference (ASO − FMA). Gray areas denote boxes in which either ASO or FMA exceeds 500 km.

  • View in gallery
    Fig. 15.

    Seasonal wavenumber spectra for (a) surface KE and (b) SSH variance in the 500 km × 500 km box centered at 37°N, 164°E. Thick (thin) lines denote balanced (unbalanced) motion spectra and red (blue) lines denote FMA (ASO) spectra. Red (blue) circles denote the FMA (ASO) transition scales of Lt in (a) and in (b). The gray line in (b) shows the baseline requirement specified by the SWOT Science Team (Rodriguez 2016).

  • View in gallery
    Fig. 16.

    (a) Distribution of the transition scale Lt based on equal KE between the balanced and unbalanced motions after removing the stationary part of the diurnal and semidiurnal internal tidal signals. Areas with Lt > 500 km are blanked out. (b) Distribution of Lt difference between Fig. 4 and (a). Gray areas denote boxes in which Lt exceeds 500 km.

  • View in gallery
    Fig. 17.

    (a) As in Fig. 16a, but for the transition scale . (b) Distribution of difference between Fig. 12a and (a).

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1494 651 58
PDF Downloads 1290 442 51

Seasonality in Transition Scale from Balanced to Unbalanced Motions in the World Ocean

Bo QiuDepartment of Oceanography, University of Hawai‘i at Mānoa, Honolulu, Hawaii

Search for other papers by Bo Qiu in
Current site
Google Scholar
PubMed
Close
,
Shuiming ChenDepartment of Oceanography, University of Hawai‘i at Mānoa, Honolulu, Hawaii

Search for other papers by Shuiming Chen in
Current site
Google Scholar
PubMed
Close
,
Patrice KleinLOPS-CNRS, Ifremer, France, and California Institute of Technology, Pasadena, California

Search for other papers by Patrice Klein in
Current site
Google Scholar
PubMed
Close
,
Jinbo WangJet Propulsion Laboratory, California Institute of Technology, Pasadena, California

Search for other papers by Jinbo Wang in
Current site
Google Scholar
PubMed
Close
,
Hector TorresJet Propulsion Laboratory, California Institute of Technology, Pasadena, California

Search for other papers by Hector Torres in
Current site
Google Scholar
PubMed
Close
,
Lee-Lueng FuJet Propulsion Laboratory, California Institute of Technology, Pasadena, California

Search for other papers by Lee-Lueng Fu in
Current site
Google Scholar
PubMed
Close
, and
Dimitris MenemenlisJet Propulsion Laboratory, California Institute of Technology, Pasadena, California

Search for other papers by Dimitris Menemenlis in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

The transition scale Lt from balanced geostrophic motions to unbalanced wave motions, including near-inertial flows, internal tides, and inertia–gravity wave continuum, is explored using the output from a global 1/48° horizontal resolution Massachusetts Institute of Technology general circulation model (MITgcm) simulation. Defined as the wavelength with equal balanced and unbalanced motion kinetic energy (KE) spectral density, Lt is detected to be geographically highly inhomogeneous: it falls below 40 km in the western boundary current and Antarctic Circumpolar Current regions, increases to 40–100 km in the interior subtropical and subpolar gyres, and exceeds, in general, 200 km in the tropical oceans. With the exception of the Pacific and Indian sectors of the Southern Ocean, the seasonal KE fluctuations of the surface balanced and unbalanced motions are out of phase because of the occurrence of mixed layer instability in winter and trapping of unbalanced motion KE in shallow mixed layer in summer. The combined effect of these seasonal changes renders Lt to be 20 km during winter in 80% of the Northern Hemisphere oceans between 25° and 45°N and all of the Southern Hemisphere oceans south of 25°S. The transition scale’s geographical and seasonal changes are highly relevant to the forthcoming Surface Water and Ocean Topography (SWOT) mission. To improve the detection of balanced submesoscale signals from SWOT, especially in the tropical oceans, efforts to remove stationary internal tidal signals are called for.

Denotes content that is immediately available upon publication as open access.

© 2018 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Bo Qiu, bo@soest.hawaii.edu

Abstract

The transition scale Lt from balanced geostrophic motions to unbalanced wave motions, including near-inertial flows, internal tides, and inertia–gravity wave continuum, is explored using the output from a global 1/48° horizontal resolution Massachusetts Institute of Technology general circulation model (MITgcm) simulation. Defined as the wavelength with equal balanced and unbalanced motion kinetic energy (KE) spectral density, Lt is detected to be geographically highly inhomogeneous: it falls below 40 km in the western boundary current and Antarctic Circumpolar Current regions, increases to 40–100 km in the interior subtropical and subpolar gyres, and exceeds, in general, 200 km in the tropical oceans. With the exception of the Pacific and Indian sectors of the Southern Ocean, the seasonal KE fluctuations of the surface balanced and unbalanced motions are out of phase because of the occurrence of mixed layer instability in winter and trapping of unbalanced motion KE in shallow mixed layer in summer. The combined effect of these seasonal changes renders Lt to be 20 km during winter in 80% of the Northern Hemisphere oceans between 25° and 45°N and all of the Southern Hemisphere oceans south of 25°S. The transition scale’s geographical and seasonal changes are highly relevant to the forthcoming Surface Water and Ocean Topography (SWOT) mission. To improve the detection of balanced submesoscale signals from SWOT, especially in the tropical oceans, efforts to remove stationary internal tidal signals are called for.

Denotes content that is immediately available upon publication as open access.

© 2018 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Bo Qiu, bo@soest.hawaii.edu

1. Introduction

Advent of nadir-looking satellite altimetry in the 1990s has revolutionized our ability to measure with high precision the global sea surface height (SSH) field and to explore, through geostrophy, the upper-ocean circulation dynamics. This exploration has significantly advanced our understanding of mesoscale oceanic variability; see Fu et al. (2010) and Morrow and Le Traon (2012) for comprehensive reviews. A critical limitation of the nadir-looking altimeters is their 100–300-km spacing between the satellite ground tracks. Even with combined data merged from multiple altimeters, the spatial resolution in a two-dimensional SSH map is typically on the order of 200 km in wavelength (Ducet et al. 2000; Chelton et al. 2011). This length scale is inadequate to fully capture the mesoscale oceanic signals that contain 90% of the kinetic energy of the ocean (Ferrari and Wunsch 2009) and misses completely the submesoscales that have length scales shorter than the Rossby radius of deformation.

With the wide-swath radar interferometry, the next-generation Surface Water and Ocean Topography (SWOT) satellite will improve the measured SSH resolution potentially down to the spectral wavelength of 15 km, allowing us to investigate for the first time the global upper-ocean circulation variability at the short mesoscale and submesoscale ranges (Fu and Ubelmann 2014). Dynamically, upper-ocean processes with length scales of O(15–200) km are important because they determine the equilibrium state of the upper ocean through the turbulent kinetic energy cascade and energy dissipation (e.g., Ferrari and Wunsch 2009; McWilliams 2016). They are also crucial to how the surface ocean communicates with the subsurface interior ocean, affecting the mixed layer evolution and upper-ocean thermal anomalies (e.g., Lapeyre et al. 2006; Capet et al. 2008; Thomas et al. 2008; Klein et al. 2008). Because of the correlation existing between surface and interior ocean potential vorticity anomalies, meso- and submesoscale SSH information in the 15–200-km range can potentially be used to reconstruct the three-dimensional, upper-ocean circulation structures, including the balanced vertical velocity field (Lapeyre and Klein 2006). A rich literature is available relating to this research topic and an up-to-date review of the topic and its application to the planned SWOT SSH measurements can be found in Qiu et al. (2016).

With the decrease in horizontal scale, the kinetic energy level of the geostrophically balanced motion has been observed to weaken proportional to , where kh is the horizontal wavenumber (Callies and Ferrari 2013; Bühler et al. 2014; Rocha et al. 2016a). For the 15–200-km range pertinent to the SWOT mission, it is possible that the balanced motion can lose its dominance and be overtaken by the unbalanced wave motions, including near-inertial flows, internal tides, and inertia–gravity waves. In such a case, the altimeter-measured SSH information can no longer be used readily to infer the time-varying surface geostrophic flows.

The transition scale separating the balanced and unbalanced motions was explored in our recent analyses of repeat shipboard ADCP measurements along 137°E across the wind-driven tropical and subtropical gyres in the northwestern Pacific (Qiu et al. 2017). By analyzing the repeat ADCP surveys from 2004 to 2016, we found that the observed transition scale depends sensitively on the energy level of local mesoscale eddy variability. In the eddy-abundant western boundary current band of Kuroshio, the transition scale is short, <15 km, and it increases to ~50 km in the moderately unstable Subtropical Countercurrent (STCC) band in the central subtropical gyre. Along the path of relatively stable North Equatorial Current, the transition scale was observed to exceed 200 km. In addition to these geographical changes, the repeat ADCP measurements along 137°E further revealed that the transition scale is decreased in winter and lengthened in summer in the Kuroshio and STCC bands because of the occurrence of mixed layer instability in winter (Sasaki et al. 2014; Qiu et al. 2014). No apparent seasonal difference in transition scale was observed in the tropical gyre.

Information about the seasonally and geographically varying transition scale is important for the SWOT mission because it provides an important length scale measure below which the SSH data may no longer be used to accurately infer the surface geostrophic velocity. To expand our analyses on transition scale to the global ocean, we adopt in this study the state-of-the-art global ocean simulation output from the Massachusetts Institute of Technology general circulation model (MITgcm) llc4320. With its horizontal resolution at 1/48° and with the realistic tidal forcings, the output of llc4320 has been used recently to explore the submesoscale dynamics in different parts of the World Ocean (Rocha et al. 2016a,b), to compare internal gravity wave SSH spectra against globally available moored measurements (Savage et al. 2017a,b), and to design calibration and validation (CalVal) approaches for SWOT-measured SSH signals (Wang et al. 2018).

This paper is organized as follows: In section 2, we provide a brief description about the MITgcm llc4320; its performance in resolving the global mesoscale eddy variability is compared to that observed by available nadir-looking satellite altimeters. In section 3, we introduce the transition scale based on kinetic energy (KE) equality of the balanced and unbalanced motions and examine its geographical variations. Seasonal variations in transition scale are examined in section 4 where opposing seasonal contributions from the balanced and unbalanced motions are emphasized. Section 5 explores the transition scale that is based on the SSH variance equality of the balanced and unbalanced motions. Its differences from the KE-based transition scale will be contrasted. In section 6, we summarize the analysis results in the context of the planned SWOT mission and discuss the need to reduce the unbalanced motion KE by removal of stationary internal tidal signals.

2. MITgcm llc4320 simulation

The data product analyzed in this study is based on the state-of-the-art global ocean simulation of MITgcm (Marshall et al. 1997) in a latitude–longitude–cap (llc) configuration with a polar cap that has 4320 × 4320 grid cells. Referred to as llc4320, the model has a 1/48° horizontal resolution and 90 vertical levels. To better resolve the submesoscale upper-ocean variability, the model’s vertical levels have a ~1-m vertical resolution at the surface and increases to ~30 m near the 500-m depth. The llc4320 simulation is initialized from the output of the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2), project (Menemenlis et al. 2008), with the subsequent model resolution increased to llc1080 (1/12°), llc2160 (1/24°), and the present llc4320 (1/48°). The llc4320 model is forced by the 6-hourly ERA-Interim atmospheric reanalysis as well as by a synthetic surface pressure field consisting of the 16 most dominant tidal constituents. For our analysis, we use the hourly output of the modeled surface height η and horizontal velocity (u, υ) data from 1 November 2011 to 31 October 2012 (366 days).

To verify the surface circulation changes modeled in the llc4320 simulation, we use the global SSH anomaly dataset generated by Ssalto/Duacs and distributed by AVISO with support from CNES (http://marine.copernicus.eu/). This dataset merges along-track SSH measurements from all satellite altimeter missions after October 1992 and has a 7-day temporal resolution and a 1/4° longitude Mercator spatial resolution (Ducet et al. 2000). The data period available for this study extends from January 1993 to December 2016.

Figure 1 compares the surface eddy kinetic energy (EKE) distribution simulated by llc4320 and that derived from the AVISO SSH product for the period of November 2011–October 2012. For comparison with the AVISO-derived EKE values, we average the llc4320 modeled SSH data into a weekly time series on the 1/4° longitude Mercator grid and compute the EKE assuming geostrophy:
e1
where g is gravitational acceleration, f is the Coriolis parameter, and η′ is the SSH anomalies from the mean SSH of November 2011–October 2012. By and large, llc4320 simulates faithfully the mesoscale eddy variability captured by the satellite altimeter missions that have spectral scales larger than O(150) km (e.g., Chelton et al. 2011). One noticeable discrepancy between the two EKE maps is within the equatorial band of ±10°, where the AVISO-derived EKE values in Fig. 1b appear consistently larger than those based on llc4320 (Fig. 1a). One likely reason for this discrepancy can be related to the measurement errors in the AVISO SSH data that are amplified when converting to EKE because of the small f values within the equatorial band.
Fig. 1.
Fig. 1.

Surface EKE distributions from (a) MITgcm llc4320 simulation and (b) AVISO SSH product of November 2011–October 2012. For both maps, EKE is calculated geostrophically from the weekly, 1/4° × 1/4°–resolution SSH field.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

Evaluating the llc4320's skill in simulating the unbalanced wave motions is more challenging because of the lack of in situ high-frequency SSH measurements in the World Ocean. By utilizing the hourly dynamic height time series from nine McLane profiler moorings available in the Pacific, North Atlantic, and south Indian Oceans, Savage et al. (2017b) found recently that the wave frequency spectra of SSH variance simulated by llc4320 agreed reasonably well with those observed by the profiler moorings. Given these high-frequency SSH comparisons and those for the mesoscale eddy variability shown in Fig. 1, we believe the llc4320 results are adequate for quantifying the relative importance between the balanced versus unbalanced motions.

3. Balanced–unbalanced motion decomposition and Lt

Delineation of balanced and unbalanced motions can be conducted in several different ways. One common method is to select a temporal filter with the low-pass (high pass) filtered signals regarded as the balanced (unbalanced) motions. For example, Richman et al. (2012) used a 48-h filter in their study on KE spectral slopes relevant to the mesoscale variability. While straightforward in its implementation, such a single-frequency filter may not be ideal for decomposing the balanced and unbalanced motions on the global scale. Figure 2 compares the typical horizontal wavenumber–frequency (khω) spectra of surface KE in four representative areas in the Pacific Ocean from the hourly output of llc4320:1 (i) the highly turbulent Kuroshio Extension, (ii) moderately unstable STCC, (iii) stable North Equatorial Current, and (iv) the tropical southwest Pacific. In the plots, elevated unbalanced motion KE is discernible in the near-inertial, diurnal, and semidiurnal tidal frequency bands as well as along the dispersion curves of inertia–gravity waves (IGWs) of discrete vertical modes: (see white solid lines), where cn (n = 1 to 10) are evaluated based on the World Ocean Atlas 2013 (WOA2013) climatological temperature–salinity profiles (Locarnini et al. 2013; Zweng et al. 2013). From Fig. 2, it is clear that the cutoff period delineating the balanced and unbalanced motions is highly sensitive to the local inertial period: near the local inertial frequency band, balanced and unbalanced motions can coexist with large and small wavenumbers, respectively. As a result, this renders the adoption of a single-frequency filter, such as the dashed black line showing the 36-h filter, undesirable.

Fig. 2.
Fig. 2.

Horizontal wavenumber–frequency distributions of surface kinetic energy in the 500 km × 500 km box centered on (a) 37°N, 164°E, (b) 21°N, 164°E, (c) 9°N, 164°E, and (d) 7°S, 164°E of the Pacific Ocean. Based on hourly MITgcm output of 1 Oct 2011–30 Sep 2012. Dashed white lines denote the inertial and tidal frequencies. Solid white lines denote the dispersion relation for inertia–gravity waves of the first 10 vertical modes. Solid black lines denote the lower-frequency boundary of either the IGW dispersion curve at the equatorward edge of the box or the local permissible tides. Dashed black lines denote the frequency of the 36-h filter. Tidal constituent S2, M2, K1, O1, and inertial frequencies are denoted by dashed white lines.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

In the present study, we adopt a more dynamically based method by defining the balanced (unbalanced) motions below (above) the lower frequency of either the local tenth vertical-mode IGW dispersion curve or the permissible tides (see thick black lines in Fig. 2). To evaluate the transition scale Lt from balanced to unbalanced motion, we first select the hourly surface (u, υ) data from llc4320 in each fixed 500 km × 500 km box. After removing the best-fit two-dimensional spatial slopes from the (u, υ) data and applying a two-dimensional Hanning window, we obtain the horizontal wavenumber–frequency spectrum of surface KE, , such as those shown in Fig. 2. The transition scale Lt is defined as the wavelength at which KE integrated over the balanced motion frequencies;
e2
is equal to that integrated over the unbalanced motion frequencies
e3
where ωyr denotes the annual frequency, ωcr(kh) is the dynamically delineating frequency defined above, and ωNy is the Nyquist frequency at 0.5 cph. As an illustration, we plot in Fig. 3 KEbal and KEunb as a function of kh in the same four representative boxes as those selected in Fig. 2. By our dynamical definition, the transition scale Lt is given by the vertical dashed lines where KEbal = KEunb.
Fig. 3.
Fig. 3.

Horizontal wavenumber spectra of the surface eddy kinetic energy of balanced (thick black lines) and unbalanced (thin black lines) motions in the 500 km × 500 km box centered on (a) 37°N, 164°E, (b) 21°N, 164°E, (c) 9°N, 164°E, and (d) 7°S, 164°E of the Pacific Ocean. Black dashed lines denote the reference spectra and vertical gray dashed lines denote the transition scale Lt.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

Figure 4 shows the global Lt distribution using the dynamical delineation method described above. Areas where Lt exceeds 500 km are blanked out. In accordance with the shipboard ADCP data results observed along 137°E in the northwestern Pacific (Qiu et al. 2017), Lt is short, <20 km, in the western boundary current Kuroshio region south of Japan. It increases to ~50 km in the weakly baroclinically unstable STCC band of 14°–28°N. The Lt value exceeds 250 km along the dynamically stable North Equatorial Current band of 9°–14°N. In the North Equatorial Countercurrent (NECC) band of 2°–9°N, Lt decreases to ~200 km. Notice the Lt value in the NECC band is larger than that inferred from the shipboard ADCP measurements (~80 km). This difference could be because the 1-yr llc4320 simulation is too short to account for the dominant interannual variability of the NECC, which was captured in the long-term shipboard ADCP data (Qiu and Joyce 1992; Chen et al. 2016).

Fig. 4.
Fig. 4.

Distribution of the transition scale Lt indicating the dominance from balanced to unbalanced motions. Areas with Lt > 500 km are blanked out. White contours denote the 3000-m water depth.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

These results in the northwestern Pacific are quite representative for the rest of the global ocean. The Lt values smaller than 20 km can be seen in all other western boundary current regions: the Kuroshio Extension, Gulf Stream, Agulhas Current, and Brazil–Malvinas Confluence. One exception is in the East Australian Current (EAC) region, which will be commented on below. Another region with small Lt is along the Antarctic Circumpolar Current (ACC) in the Southern Ocean, with the exceptions where prominent bottom topographic features exist, for example, the Kerguelen Plateau near 70°E, the East Scotia Basin east of the Drake Passage, and the Spiess Seamount near the southern tip of the Mid-Atlantic Ridge. In the Northern Hemisphere’s temperate latitudes of 15°–30°N, Lt falls generally in the 50–100-km range. The corresponding latitudes in the Southern Hemisphere tend to have Lt in the larger O(150) km range. In the tropics equatorward of 15°, the Lt values often exceed 200 km with the unbalanced motion KE overwhelming the balanced motion KE. In the high latitudes, one area where Lt exceeds 200 km is in the Alaskan Gyre, where the balanced oceanic mesoscale eddy variability is extremely weak (recall Fig. 1).

That Lt is small in the western boundary current and ACC regions is not surprising because these are the regions in the World Ocean known for their large mesoscale eddy variability (e.g., Fu et al. 2010; Morrow and Le Traon 2012). Indeed, Fig. 5a shows the global map of the balanced motion KE integrated over the wavenumber space, and it exhibits a spatial distribution very similar to that of the mesoscale EKE map in Fig. 1a. It is, however, important to note that the level of mesoscale eddy variability is not the sole determinant for Lt. As shown in Fig. 5b, the wavenumber-integrated unbalanced motion KE level is similarly geographically nonuniform. It is much higher in regions with prominent bottom topographic features, like the Philippine Sea in the northwestern Pacific; the Aleutian Islands in the subarctic North Pacific; the Solomon, Coral, and Tasman Seas in the southwestern Pacific; the Tuamoto Archipelago in the central South Pacific; northeast of Madagascar in the Indian Ocean; and the continental slopes off the northeast coast of United States. Many of these locations, not coincidently, correspond to the major generation sites of internal tides in the World Ocean (Ray and Zaron 2016; Zhao et al. 2016; and references therein). In fact, it is this enhanced KE level associated with the unbalanced motions that is responsible for the elevated Lt values at ~150 km in the Coral and Tasman Seas in Fig. 4, albeit the presence of the western boundary current EAC in the region. Likewise, the abyssal seamount and midocean ridge–enhanced unbalanced motions are responsible for the spatially dependent Lt modulations along the ACC path in the Southern Ocean.

Fig. 5.
Fig. 5.

Distributions of wavenumber-integrated EKE for (a) balanced and (b) unbalanced motions based on the dynamical delineation defined in Eqs. (2) and (3). Note that color scales for the EKE level are different for (a) and (b). White contours denote the 3000-m water depth.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

4. Seasonal variations in Lt

There exists increasing evidence in recent years that near-surface meso- and submesoscale variability in the 10–100-km range has a distinct annual cycle relating to the occurrence of mixed layer instability. Such evidence emerges from both in situ measurements in different parts of the World Ocean (e.g., Callies et al. 2015; Buckingham et al. 2016; Thompson et al. 2016; Qiu et al. 2017) and high-resolution ocean general circulation model simulations (e.g., Mensa et al. 2013; Sasaki et al. 2014; Qiu et al. 2014). Since this 10–100-km range overlaps extensively with the Lt values obtained in the preceding section, it is natural to inquire if our above-defined transition scale Lt is influenced by the seasonally modulating meso- and submesoscale variability. Our interest in the seasonal Lt changes is further motivated by a recent study by Rocha et al. (2016b), who found that the unbalanced motion kinetic energy in the Kuroshio Extension region in the MITgcm llc4320 exhibited a seasonal maximum in summer, opposite to that of the balanced motion kinetic energy.

To examine the seasonal changes in Lt, we define winter (summer) as February–April (August–October) and evaluate Lt according to Eqs. (2) and (3) with the year-long hourly (u, υ) data replaced by the seasonal ones (Figs. 6a,b). While the overall geographical dependence of the seasonal Lt values is similar to the annual-mean pattern shown in Fig. 4, consistent changes with larger (smaller) Lt values in the summer (winter) hemisphere are visually identifiable in Fig. 6. To better quantify the seasonal difference, we plot in Fig. 7 the Lt difference of summer [August–October (ASO)] minus winter [February–April (FMA)]. In agreement with the visual changes seen between Figs. 6a and 6b, the ASO − FMA Lt is by and large positive in the Northern Hemisphere (NH) and negative in the Southern Hemisphere (SH). Notice that the differential Lt values are in the 50–150-km range in much of the off-equatorial ocean, and this range is on par with the magnitude of the annual-mean Lt values presented in Fig. 4.

Fig. 6.
Fig. 6.

Distributions of seasonal transition scale Lt in (a) February–April and (b) August–October. Areas with Lt > 500 km are blanked out.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

Fig. 7.
Fig. 7.

Distribution of seasonal difference (ASO − FMA) in Lt or Fig. 6b minus Fig. 6a. Gray areas denote boxes in which either ASO or FMA Lt exceeds 500 km.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

An interesting exception to the SH’s general negative ASO − FMA Lt occurs in the Pacific and Indian sectors of the Southern Ocean where the SH winter-minus-summer Lt is small but remains positive. To clarify the cause responsible for this regional exception, we evaluate the ASO − FMA balanced and unbalanced motion KE levels in the wavelength range <100 km of our interest (Figs. 8a,b). Because of the prevalence of winter mixed layer instability, the balanced motion KE level in ASO is clearly lower (higher) than in FMA in the NH (SH) oceans. This is no exception in the Pacific and Indian sectors of the Southern Ocean. For the unbalanced motion ASO − FMA KE level, Fig. 8b reveals that while it is overall positive in the NH and negative in the SH, the differential KE level remains weakly positive in the Pacific and Indian sectors of the Southern Ocean. In other words, it is the lack of summer enhancement by the unbalanced motion KE that causes ASO − FMA Lt to be weakly positive in the Pacific and Indian sectors of the Southern Ocean.

Fig. 8.
Fig. 8.

Distributions of seasonal difference (ASO − FMA) in wavenumber-integrated EKE for (a) balanced and (b) unbalanced motions with length scales shorter than 100 km.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

In their study of submesoscale seasonality in the Kuroshio Extension region, Rocha et al. (2016b) argued that the enhancement in the unbalanced motion KE in summer is due to the shoaling of mixed layer and strengthening of upper-ocean stratification. Both of these effects work to trap the inertia–gravity waves to the shallow surface ocean, enhancing the unbalanced motion’s projection onto the mixed layer (D’Asaro 1978). An inspection of the seasonal mixed layer depth changes in llc4320 reveals that the summer (FMA) mixed layer along the path of ACC (Fig. 9a) extends to below 100 m and this is particularly true in the Pacific and Indian Ocean sectors where the mixed layer can reach to deeper than 200 m.2 This 100–200-m mixed layer depth in the summer ACC should be contrasted to the 20-m mixed layer depth value observed typically in the summer NH oceans (see Fig. 9b) or the rest of the summer SH oceans. It is this deep summer mixed layer that is likely to be responsible for the ASO − FMA unbalanced motion KE to remain positive along the ACC shown in Fig. 8b.

Fig. 9.
Fig. 9.

Distributions of seasonal mixed layer depth in (a) FMA and (b) ASO based on the llc4320 simulation.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

It is important to emphasize that the seasonal difference in the unbalanced motion KE relating to the summer trapping to within the surface mixed layer does not penetrate deep into the subsurface ocean. At the main thermocline depth of 500 m, for example, the ASO − FMA unbalanced motion KE (Fig. 10b) shows an oppositely signed distribution to that of the surface map (cf. Fig. 8b; note that the color scales are different between Figs. 8 and 10 because of the shallower depth scales associated with the submesoscale flows). Below the surface mixed layer, the wintertime unbalanced motion KE exceeds that in summer because of the enhanced surface wind forcing that elevates the KE level of the wind-induced superinertial wave motions (e.g., Alford et al. 2016). For the ASO − FMA balanced motion KE with scales < 100 km at the 500-m depth, Fig. 10a reveals that the hemispheric contrast that exists in the surface layer (recall Fig. 8a) becomes less distinctive. The ASO − FMA difference shown in Fig. 10a is spatially patchy, although the general tendency for the wintertime balanced motion KE to surpass that in summer is still discernible.

Fig. 10.
Fig. 10.

Distributions of seasonal difference (ASO − FMA) in wavenumber-integrated EKE at 500 m for (a) balanced and (b) unbalanced motions with length scales shorter than 100 km. Notice that the color scale is 10 times smaller in this figure than in Fig. 8.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

In concluding this section, we plot in Fig. 11 the cumulative probability density function (PDF) distribution for Lt in the NH and SH winter as a function of global ocean latitude. Taking the thick black lines as an example, it is clear from Fig. 11 that 80% of Lt falls below 20 km in most of the NH oceans between 25° and 45°N and in all of the SH oceans south of 25°S. That the wintertime Lt is shorter than 20 km in most of the extratropical NH and SH oceans indicates that the entanglement between balanced and unbalanced motions will likely pose less of a challenge for the SWOT mission in these locations to detect the balanced submesoscale features in winter.

Fig. 11.
Fig. 11.

Cumulative PDF distribution for wintertime Lt values (FMA for NH and ASO for SH) estimated in individual latitudinal bands of the World Ocean. Black solid lines denote the 80% threshold Lt values.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

5. SSH variance-based

In the preceding sections, we have explored the global distribution of Lt based on kinetic energy equality between the balanced and unbalanced motions. From a practical point of view, Lt thus derived could be regarded as the length scale below which surface circulation inferred geostrophically from the SSH gradient becomes unreliable. With the high-resolution SSH information from the planned SWOT mission, it is instructive to also evaluate the transition scale that is based on the SSH variance equality between the balanced and unbalanced motions. To differentiate this transition scale based on the SSH variance equality from the Lt defined in section 3, we denote hereafter the SSH variance-based Lt as .

To evaluate from the hourly output of llc4320, we adopt the same dynamical delineation method detailed in section 3 and replace the KE spectra by the SSH variance spectra in Eqs. (2)(3). Figure 12a shows the distribution of calculated from the year-long llc4320 output. Geographically, exhibits a spatial pattern very similar to that of the KE-based Lt shown in Fig. 4. Like for the balanced and unbalanced motion kinetic energy levels, the wavenumber-integrated SSH variance associated with the balanced and unbalanced motions are similarly spatially inhomogeneous (figures not shown). By and large, has a larger value than Lt as compared in Fig. 12b and the general trend is such that their difference, Lt, decreases progressively from the equator toward higher latitudes.

Fig. 12.
Fig. 12.

(a) Distribution of the transition scale based on equal SSH variance between the balanced and unbalanced motions. Areas with Lt > 500 km are blanked out. (b) Distribution of Lt. Gray areas denote boxes in which either or Lt exceeds 500 km.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

It is of interest to consider here why exceeds in general Lt and what contributes to the diminishing difference in Lt from the equator toward higher latitudes. For simplicity, let us consider a linearized shallow-water ocean in which surface velocity u and υ and SSH η are governed by (Gill 1982)
e4
e5
e6
where , and H denotes the water depth (or equivalent depth for a baroclinic mode in the stratified ocean). Rearranging Eqs. (4)(6), we obtain
e7
e8
and the wavenumber–frequency spectra for KE and SSH variance and are in this case related through
e9
where . For the balanced motion, ωf, and its KE spectrum defined in Eq. (2) can be approximated by
e10
where HVbal(kh) denotes the SSH variance spectrum associated with the balanced motion. That the KE spectral slope in Eq. (10) is shallower than the SSH variance spectral slope by has been extensively used in previous studies of mesoscale eddy variability based on satellite altimeter data (e.g., Stammer 1997).
For the unbalanced motion, combining Eqs. (3) and (9) leads to
e11
As demonstrated in Fig. 2, the unbalanced motion KE resides largely in the near-inertial, diurnal/semidiurnal tide and discrete IGW bands. Since the near-inertial motions have weak projections onto SSH (e.g., Alford et al. 2016), contributions to the RHS integration in Eq. (11) come mostly from the tidal and IGW motions. If is assumed to be dominated by the tidal or IGW motion with a particular frequency ω0, we may approximate , where δ is the Dirac delta function. In this case, Eq. (11) becomes
e12
and combining Eqs. (10) and (12) results in
e13
where . From Eq. (13), it is clear that the difference between Lt and depends on the magnitude of I. If I > 1, Lt will surpass and, if I < 1, will exceed Lt. For the semidiurnal tides, I is smaller (larger) than unity equatorward (poleward) of 35.1° latitudes; for the diurnal tides, this threshold latitude shifts to 16.7°. Both diurnal and semidiurnal tidal motions, as a consequence, act to lengthen with respect to Lt in the tropics and shorten it in higher latitudes. These tidal effects are consistent qualitatively with the Lt pattern shown in Fig. 12b.
For the discrete IGWs that follow the dispersion relationship , I becomes
e14
where cn denotes the phase speed of the mode-n internal gravity wave. In Eq. (14), I < 1 if and vice versa. For a fixed mode cn and wavenumber kh, the tendency is again for to exceed Lt in lower latitudes. Figure 13 shows the threshold lines of I = 1 for different vertical mode IGWs averaged across the latitudinal bands of the World Ocean. As in Fig. 2, cn (n = 1 to 10) are calculated based on the WOA2013 climatological temperature–salinity profiles (Locarnini et al. 2013; Zweng et al. 2013). Unlike the tidal effect noted above, the I = 1 threshold for IGWs not only depends on latitude but also on the wavenumber. Specifically, I can remain less than unity in high-latitude oceans when kh is large in the 10–100-km wavelength range of our interest. It is in fact the contributions from the low-mode IGWs that keep > Lt in some of the high-latitude areas seen in Fig. 12b.
Fig. 13.
Fig. 13.

Threshold lines of I = 1 for discrete vertical mode (n = 1 to 10) IGWs as a function of kh [see Eq. (14)]. Values are averaged across the latitudinal bands of the World Ocean. For a fixed vertical mode, I < 1 in the low-latitude and high-wavenumber space of the threshold line.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

Like the seasonal changes in Lt, the SSH variance-based is also subject to distinct seasonal modulations. Figure 14 shows the difference between ASO and FMA, and its global pattern is spatially very similar to that for Lt shown in Fig. 7. The seasonal amplitude in , on the other hand, is about half that obtained in Fig. 7. One reason for this is that the balanced motion SSH variance spectral slope is steeper than the KE spectral slope; recall Eq. (10) and compare the spectra slopes for KE and SSH variance in Fig. 15 at a representative site of 37°N and 164°E in the North Pacific. In contrast, the unbalanced motion KE and SSH variance spectra have similar kh slopes (cf. dashed lines in Fig. 15). This renders the seasonally varying crossings by balanced and unbalanced motion SSH variance spectra more restricted in the kh space than those by balanced and unbalanced motion KE spectra.

Fig. 14.
Fig. 14.

Distributions of seasonal difference (ASO − FMA). Gray areas denote boxes in which either ASO or FMA exceeds 500 km.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

Fig. 15.
Fig. 15.

Seasonal wavenumber spectra for (a) surface KE and (b) SSH variance in the 500 km × 500 km box centered at 37°N, 164°E. Thick (thin) lines denote balanced (unbalanced) motion spectra and red (blue) lines denote FMA (ASO) spectra. Red (blue) circles denote the FMA (ASO) transition scales of Lt in (a) and in (b). The gray line in (b) shows the baseline requirement specified by the SWOT Science Team (Rodriguez 2016).

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

6. Discussion and summary

We have in this study explored geographical and seasonal variations in the transition scale at which the geostrophically balanced motion loses its dominance over the unbalanced wave motions. As the metric for the transition scale, we adopted Lt defined by the wavelength that has equal balanced and unbalanced motion KE spectral densities. By analyzing the hourly output from the 1/48° horizontal resolution MITgcm llc4320, we found that Lt is geographically highly inhomogeneous. Broadly speaking, Lt falls below 40 km in the western boundary current and ACC regions and increases to 40–100 km in the vast subtropical and subpolar world oceans. An exception to the latter case is in the Alaskan Gyre where Lt can exceed 200 km. In the tropical oceans equatorward of 20° latitudes, Lt is found, in general, to be larger than 150 km. These llc4320-derived, latitude-dependent Lt values are, by and large, consistent with our recent estimates based on the repeat shipboard ADCP measurements along 137°E in the northwestern Pacific Ocean (Qiu et al. 2017).

On a regional scale, the Lt value depends not only on the local level of mesoscale eddy variability, but also on the energy level set by the unbalanced motions consisting of near-inertial flows, internal tides, and inertia–gravity waves. The unbalanced motion energy level is spatially nonuniform, and high values are commonly detected where prominent bathymetric features exist. Examples of effects on Lt by enhanced unbalanced motions include the East Australian Current region where Lt reaches 100 km (rather than 20 km as in other western boundary current regions) and along the ACC path where increased Lt values are found to be collocated with the major topographic features in the Southern Ocean.

The energy levels of both the balanced and unbalanced motions modulate seasonally and, as a result, they can alter the Lt value constructively. In winter, small-scale balanced motions are energized because of the ubiquitous occurrence of mixed layer instability. Shoaling of the mixed layer and increasing of density jump at the mixed layer base in summer, on the other hand, works to trap and elevate the surface unbalanced motion kinetic energy (Rocha et al. 2016b). The combined effect of these seasonal kinetic energy changes results in a seasonal contrast in Lt that can typically reach >100 km in significant parts of the global ocean. One important exception to this Lt change caused constructively by the seasonally varying balanced and unbalanced motions is seen in the Pacific and Indian sectors of the Southern Ocean. In these sectors, the regional mixed layer in summer remains deep, failing as a result to trap effectively the unbalanced motion kinetic energy to within the shallow surface layer.

In addition to the metric Lt based on KE equality, we have in this study also examined defined by the wavelength at which balanced and unbalanced motions have equal SSH variance spectral densities. The shows similar geographical and seasonal patterns as those of Lt, highlighting the common dynamics controlling the spatial and temporal variations between these two transition-scale metrics. While exhibiting similar patterns, the values are in general larger than Lt, especially in low-latitude oceans. The difference Lt tends to diminish toward the high-latitude oceans, where may become shorter than Lt. The reason for this tendency is that in the wavenumber space, the KE and SSH variance spectra are functionally related differently for the balanced and unbalanced motions. While the balanced motions between the two are proportional to , the unbalanced motions are connected via convolutions that are sensitive functions of latitude.

Since Lt signifies the transition scale below which the SSH data may no longer be used to accurately infer the geostrophically balanced flows, our findings about its geographical and seasonal variations are highly relevant for the forthcoming SWOT mission. With the 15-km spectral resolution expected for the SWOT mission (Fu and Ubelmann 2014), unbalanced motion signals are likely to pose little problems for detecting the balanced submesoscale flows in the mesoscale-rich regions like the ACC and the western boundary currents in the North Pacific, North Atlantic, South Atlantic, and south Indian Oceans. For the East Australian Current and in the regions of moderate mesoscale activities (e.g., the Subtropical Countercurrent bands in temperate latitudes), caution is likely warranted in diagnosing the surface geostrophic flows in the 50–100-km range.

It is worth emphasizing that the annual-mean Lt values noted above can decrease significantly in winter because the seasonal balanced and unbalanced motion KE levels fluctuate out of phase, shortening Lt constructively in winter. Indeed, our analysis revealed that 80% of the NH oceans between 25° and 45°N and all of the SH oceans south of 25°S will have Lt shorter than 20 km during the winter season, including the EAC and the Subtropical Countercurrent bands cautioned above. Since winter is the season when the balanced submesoscale motions prevail, this implies that the SWOT-measured SSH signals can be used potentially to explore the wintertime submesoscale variability down to the 15–20-km range in a significant part of the extratropical world oceans.

The serious challenge for the SWOT mission to capture the small-scale geostrophically balanced motion will be in the tropical oceans within ±20° latitudes and in the Alaskan Gyre, where Lt commonly exceeds 150 km in all seasons. This challenge stems from the fact that these regions have a high unbalanced motion KE level as compared to the balanced motion. One important source behind the high unbalanced motion KE is the internal tides generated by the interaction of barotropical tides with bottom topography. It is important to note that a significant portion of the low-mode internal tides have been known to exhibit phase stability, or coherence, in their SSH expressions (e.g., Ray and Mitchum 1996). By analyzing the internal tidal signals from the 17-yr combined record of TOPEX/Poseidon and Jason altimeters, Ray and Zaron (2011) estimated that the stationary variance is generally 75% or more of the average variance of mode-1 internal tides. More recent numerical modeling and data analysis studies by Shriver et al. (2014) and Zaron (2017) have revealed that the fraction explained by the stationary internal tides can be highly spatially inhomogeneous. In the context for Lt determination, it is obvious that any removal of internal tidal signals that are stationary and predictable will result in a reduction in the Lt value.

A judicious separation between the stationary and nonstationary internal tides requires multiyear in situ measurements or numerical model output. Rather than pursuing this separation based on the year-long llc4320 output, we attempt to quantify the reduction in Lt by assuming that a fraction of the llc4320-modeled diurnal and semidiurnal internal tidal signals are stationary and could be, as a result, removed. Specifically, we follow Zaron (2017, his Fig. 9) and assume that the variance fraction explained by the stationary internal tides has a simple latitude-dependent form:
e15
where y is the latitude in degrees. Using this form, we subjectively reduced the wavenumber KE spectra (recall Fig. 2) by 1 − α(y) in the diurnal and semidiurnal frequency bands over the background continuum and reevaluate KEunb(kh) in Eq. (3). Here, we define the diurnal (semidiurnal) frequency band as bounded by ω = 1/27.0–1/23.0 (1/13.5–1/11.5) cph and the background continuum by linearly interpolating the frequency spectra within these chosen bounds by neighboring spectral values. Figure 16a shows the annual-mean Lt distribution based on the newly evaluated KEunb(kh), while its difference from the Lt values shown in Fig. 4 is presented in Fig. 16b. From Fig. 16b, it is clear that removal of the stationary internal tidal signals can reduce Lt by 50–100 km in many regions of the World Ocean; significant reduction in Lt is seen in the Alaskan Gyre, the northwestern tropical Pacific, and the broad tropical and subtropical SH oceans. Because of the smaller fraction of the stationary internal tides, the reduction in Lt is modest, generally less than 50 km, within the ±10° equatorial band. For completeness, we present the same analysis results in Fig. 17 with respect to . Although having smaller amplitude when compared to Lt, the reduction in caused by the removal of stationary internal tidal signals exhibits a similarly geographical pattern.
Fig. 16.
Fig. 16.

(a) Distribution of the transition scale Lt based on equal KE between the balanced and unbalanced motions after removing the stationary part of the diurnal and semidiurnal internal tidal signals. Areas with Lt > 500 km are blanked out. (b) Distribution of Lt difference between Fig. 4 and (a). Gray areas denote boxes in which Lt exceeds 500 km.

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

Fig. 17.
Fig. 17.

(a) As in Fig. 16a, but for the transition scale . (b) Distribution of difference between Fig. 12a and (a).

Citation: Journal of Physical Oceanography 48, 3; 10.1175/JPO-D-17-0169.1

The results shown in Figs. 16 and 17 are based on our current knowledge about the distributions of the stationary internal tides. Because of the background changes in ocean stratification, currents, and mesoscale eddies, isolating and removing the stationary internal tidal signals from the SWOT SSH measurements will be a formidable endeavor. Nevertheless, we believe the results of Figs. 16 and 17 suggest that this endeavor is critical and should be pursued in order to improve detection of the balanced submesoscale signals and maximize the science returns by the SWOT mission.

Acknowledgments

We thank Rosemary Morrow, Clement Ubelmann, and Ed Zaron for insightful discussions. Constructive comments made by two anonymous reviewers helped improve an early version of the manuscript. B. Q. and S. C. acknowledge support from NASA SWOT and OSTST missions (NNX16AH66G and NNX17AH33G). P. K. acknowledges the support of the Agence Nationale pour la Recherche [ANR-10-LABX-19-01 (LabexMER)] and of the SWOT project (JPL, NASA). J. W., H. T., L. F., and D. M.’s research presented in this paper was carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under contract with NASA. They acknowledge support from the SWOT projects.

REFERENCES

  • Alford, M. H., J. A. MacKinnon, H. L. Simmons, and J. D. Nash, 2016: Near-inertial internal gravity waves in the ocean. Annu. Rev. Mar. Sci., 8, 95123, https://doi.org/10.1146/annurev-marine-010814-015746.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Buckingham, C. E., and Coauthors, 2016: Seasonality of submesoscale flows in the ocean surface boundary layer. Geophys. Res. Lett., 43, 21182126, https://doi.org/10.1002/2016GL068009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bühler, O., J. Callies, and R. Ferrari, 2014: Wave–vortex decomposition of one-dimensional ship-track data. J. Fluid Mech., 756, 10071026, https://doi.org/10.1017/jfm.2014.488.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., and R. Ferrari, 2013: Interpreting energy and tracer spectra of upper ocean turbulence in the submesoscale range (1–200 km). J. Phys. Oceanogr., 43, 24562474, https://doi.org/10.1175/JPO-D-13-063.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Callies, J., R. Ferrari, J. M. Klymak, and J. Gula, 2015: Seasonality in submesoscale turbulence. Nat. Commun., 6, 6862, https://doi.org/10.1038/ncomms7862.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Capet, X., J. C. McWilliams, M. J. Molemaker, and A. J. Shchepetkin, 2008: Mesoscale to submesoscale transition in the California Current System. Part I: Flow structure, eddy flux, and observational tests. J. Phys. Oceanogr., 38, 2943, https://doi.org/10.1175/2007JPO3671.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chelton, D. B., M. G. Schlax, and R. M. Samelson, 2011: Global observations of nonlinear mesoscale eddies. Prog. Oceanogr., 91, 167216, https://doi.org/10.1016/j.pocean.2011.01.002.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chen, X., B. Qiu, Y. Dan, S. Chen, and Y. Qi, 2016: Interannual and interdecadal variability of the North Equatorial Countercurrent in the western Pacific. J. Geophys. Res. Oceans, 121, 77437758, https://doi.org/10.1002/2016JC012190.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • D’Asaro, E. A., 1978: Mixed layer velocities induced by internal wave. J. Geophys. Res., 83, 24372438, https://doi.org/10.1029/JC083iC05p02437.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • de Boyer Montégut, C., G. Madec, A. S. Fischer, A. Lazar, and D. Iudicone, 2004: Mixed layer depth over the global ocean: An examination of profile data and a profile-based climatology. J. Geophys. Res., 109, C12003, https://doi.org/10.1029/2004JC002378.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ducet, N., P.-Y. Le Traon, and G. Reverdin, 2000: Global high-resolution mapping of ocean circulation from TOPEX/Poseidon and ERS-1 and -2. J. Geophys. Res., 105, 19 47719 498, https://doi.org/10.1029/2000JC900063.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ferrari, R., and C. Wunsch, 2009: Ocean circulation kinetic energy: Reservoirs, sources, and sinks. Annu. Rev. Fluid Mech., 41, 253282, https://doi.org/10.1146/annurev.fluid.40.111406.102139.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fu, L.-L., and C. Ubelmann, 2014: On the transition from profile altimetry to swath altimeter for observing global ocean surface topography. J. Atmos. Oceanic Technol., 31, 560568, https://doi.org/10.1175/JTECH-D-13-00109.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fu, L.-L., D. B. Chelton, P.-Y. Le Traon, and R. Morrow, 2010: Eddy dynamics from satellite altimetry. Oceanography, 23 (4), 1425, https://doi.org/10.5670/oceanog.2010.02.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gill, A. E., 1982: Atmosphere–Ocean Dynamics, Academic Press, 662 pp.

  • Klein, P., B. L. Hua, G. Lapeyre, X. Capet, S. Le Gentil, and H. Sasaki, 2008: Upper ocean turbulence from high-resolution 3D simulations. J. Phys. Oceanogr., 38, 17481763, https://doi.org/10.1175/2007JPO3773.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lapeyre, G., and P. Klein, 2006: Dynamics of the upper oceanic layers in terms of surface quasigeostrophy theory. J. Phys. Oceanogr., 36, 165176, https://doi.org/10.1175/JPO2840.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lapeyre, G., P. Klein, and B. L. Hua, 2006: Oceanic restratification forced by surface frontogenesis. J. Phys. Oceanogr., 36, 15771590, https://doi.org/10.1175/JPO2923.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Locarnini, R. A., and Coauthors, 2013: Temperature. Vol. 1, World Ocean Atlas 2013, NOAA Atlas NESDIS 73, 40 pp.

  • Marshall, J. C., A. Adcroft, C. Hill, L. Perelman, and C. Heisey, 1997: A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers. J. Geophys. Res., 102, 57535766, https://doi.org/10.1029/96JC02775.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McWilliams, J. C., 2016: Submesoscale currents in the ocean. Proc. Roy. Soc., A472, 20160177, https://doi.org/10.1098/rspa.2016.0117.

  • Menemenlis, D., J.-M. Campin, P. Heimbach, C. Hill, T. Lee, A. Nguyen, M. Schodlok, and H. Zhang, 2008: ECCO2: High-resolution global ocean and sea ice data synthesis. Mercator Ocean Quarterly Newsletter, No. 31, Mercator Ocean, Ramonville Saint-Agne, France, 13–21.

  • Mensa, J. A., Z. Garraffo, A. Griffa, T. M. Ozgokmen, A. Haza, and M. Veneziani, 2013: Seasonality of the submesoscale dynamics in the Gulf Stream region. Ocean Dyn., 63, 923941, https://doi.org/10.1007/s10236-013-0633-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Morrow, R., and P.-Y. Le Traon, 2012: Recent advances in observing mesoscale ocean dynamics with satellite altimetry. Adv. Space Res., 50, 10621076, https://doi.org/10.1016/j.asr.2011.09.033.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., and T. M. Joyce, 1992: Interannual variability in the mid- and low-latitude western North Pacific. J. Phys. Oceanogr., 22, 10621079, https://doi.org/10.1175/1520-0485(1992)022<1062:IVITMA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., S. Chen, P. Klein, H. Sasaki, and Y. Sasai, 2014: Seasonal mesoscale and submesoscale eddy variability along the North Pacific Subtropical Countercurrent. J. Phys. Oceanogr., 44, 30793098, https://doi.org/10.1175/JPO-D-14-0071.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., S. Chen, P. Klein, C. Ubelmann, L.-L. Fu, and H. Sasaki, 2016: Reconstructability of three-dimensional upper-ocean circulation from SWOT sea surface height measurements. J. Phys. Oceanogr., 46, 947963, https://doi.org/10.1175/JPO-D-15-0188.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Qiu, B., T. Nakano, S. Chen, and P. Klein, 2017: Submesoscale transition from geostrophic flows to internal waves in the northwestern Pacific upper ocean. Nat. Commun., 8, 14055, https://doi.org/10.1038/ncomms14055.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ray, R. D., and G. T. Mitchum, 1996: Surface manifestation of internal tides generated near Hawaii. Geophys. Res. Lett., 23, 21012104, https://doi.org/10.1029/96GL02050.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ray, R. D., and E. D. Zaron, 2011: Non-stationary internal tides observed with satellite altimetry. Geophys. Res. Lett., 38, L17609, https://doi.org/10.1029/2011GL048617.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ray, R. D., and E. D. Zaron, 2016: M2 internal tides and their observed wavenumber spectra from satellite altimetry. J. Phys. Oceanogr., 46, 322, https://doi.org/10.1175/JPO-D-15-0065.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Richman, J. G., B. K. Arbic, J. F. Shriver, E. J. Metzger, and A. J. Wallcraft, 2012: Inferring dynamics from the wavenumber spectra of an eddying global ocean model with embedded tides. J. Geophys. Res., 117, C12012, https://doi.org/10.1029/2012JC008364.

    • Search Google Scholar
    • Export Citation
  • Rio, M.-H., S. Guinehut, and G. Larnico, 2011: New CNES-CLS09 global mean dynamic topography computed from the combination of GRACE data, altimetry, and in situ measurements. J. Geophys. Res., 116, C07018, https://doi.org/10.1029/2010JC006505.

    • Search Google Scholar
    • Export Citation
  • Rocha, C. B., T. K. Chereskin, S. T. Gille, and D. Menemenlis, 2016a: Mesoscale to submesoscale wavenumber spectra in Drake Passage. J. Phys. Oceanogr., 46, 601620, https://doi.org/10.1175/JPO-D-15-0087.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rocha, C. B., S. T. Gille, T. K. Chereskin, and D. Menemenlis, 2016b: Seasonality of submesoscale dynamics in the Kuroshio Extension. Geophys. Res. Lett., 43, 11 30411 311, https://doi.org/10.1002/2016GL071349.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rodriguez, E., 2016: Surface Water and Ocean Topography Mission (SWOT) project: Science requirements document. NASA JPL Rep. D-61923, 28 pp., https://swot.oceansciences.org/docs/D-61923_SRD_Rev_A_20160318.pdf.

  • Sasaki, H., P. Klein, B. Qiu, and Y. Sasai, 2014: Impact of oceanic-scale interactions on the seasonal modulation of ocean dynamics by the atmosphere. Nat. Commun., 5, 5636, https://doi.org/10.1038/ncomms6636.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Savage, A. C., and Coauthors, 2017a: Frequency content of sea surface height variability from internal gravity waves to mesoscale eddies. J. Geophys. Res. Oceans, 122, 25192538, https://doi.org/10.1002/2016JC012331.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Savage, A. C., and Coauthors, 2017b: Spectral decomposition of internal gravity wave sea surface height in global models. J. Geophys. Res. Oceans, 122, 78037821, https://doi.org/10.1002/2017JC013009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Shriver, J. F., J. G. Richman, and B. K. Arbic, 2014: How stationary are the internal tides in a high-resolution global ocean circulation model? J. Geophys. Res. Oceans, 119, 27692787, https://doi.org/10.1002/2013JC009423.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stammer, D., 1997: Global characteristics of ocean variability estimated from regional TOPEX/Poseidon altimeter measurements. J. Phys. Oceanogr., 27, 17431769, https://doi.org/10.1175/1520-0485(1997)027<1743:GCOOVE>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thomas, L. N., A. Tandon, and A. Mahadevan, 2008: Submesoscale processes and dynamics. Ocean Modeling in an Eddying Regime, Geophys. Monogr., Vol. 177, Amer. Geophys. Union, 17–38.

    • Crossref
    • Export Citation
  • Thompson, A. F., A. Lazar, C. Buckingham, A. C. Naveira Garabato, G. M. Damerell, and K. J. Heywood, 2016: Open-ocean submesoscale motions: A full seasonal cycle of mixed layer instabilities from gliders. J. Phys. Oceanogr., 46, 12851307, https://doi.org/10.1175/JPO-D-15-0170.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wang, J., L.-L. Fu, B. Qiu, D. Menemenlis, T. Farrar, Y. Chao, A. Thompson, and M. Flexas, 2018: An observing system simulation experiment for the calibration and validation of the Surface Water and Ocean Topography sea surface height measurement using in-situ platforms. J. Atmos. Oceanic Technol., 35, 281297, https://doi.org/10.1175/JTECH-D-17-0076.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zaron, E. D., 2017: Mapping the nonstationary internal tide with satellite altimetry. J. Geophys. Res. Oceans, 122, 539554, https://doi.org/10.1002/2016JC012487.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zhao, Z., M. H. Alford, J. B. Girton, L. Rainville, and H. L. Simmons, 2016: Global observations of open-ocean mode-1 M2 internal tides. J. Phys. Oceanogr., 46, 16571684, https://doi.org/10.1175/JPO-D-15-0105.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zweng, M. M., and Coauthors, 2013: Salinity. Vol. 2, World Ocean Atlas 2013, NOAA Atlas NESDIS 74, 39 pp.

    • Crossref
    • Export Citation
1

In light of SSH information from the SWOT mission, barotropical tidal signals are removed throughout this study. Specifically, this is done by eliminating the best-fit two-dimensional spatial slopes from the (u, υ) data in individual 500 km × 500 km boxes prior to the spectral wavenumber–frequency decompositions.

2

Similar seasonally varying spatial patterns of the mixed layer depth can be found in the mixed layer depth climatology based on in situ measurements; see, for example, Fig. 5 in de Boyer Montégut et al. (2004).

Save