Description: Agulhas leakage, the warm and salty inflow of Indian Ocean water into the Atlantic Ocean, is of importance for the climate-relevant Atlantic Meridional Overturning Circulation. South of Africa, the eastward turning Agulhas Current sheds Agulhas rings, cyclones and filaments of order 100 km that carry the Indian Ocean water into the Cape Basin and further into the Atlantic. Here, we show that the resolution of submesoscale flows of order 10 km in an ocean model leads to 40 % more Agulhas leakage and more realistic Cape Basin water-masses compared to a parallel non-submesoscale resolving simulation. Moreover, we show that submesoscale flows strengthen shear-edge eddies and in consequence lee cyclones at the northern edge of the Agulhas Current, as well as the leakage pathway in the region of the filaments that takes place outside of mesoscale eddies. This indicates that the increase in leakage can be attributed to stronger Agulhas filaments, when submesoscale flows are resolved.

Description: Leakage of warm, salty waters from the Indian Ocean into the Atlantic increases by up to 40 % in high-resolution numerical ocean model simulations, suggesting that low-resolution models underestimate this key part of the global meridional overturning circulation.

Description: Agence Nationale de la Recherche (French National Research Agency) https://doi.org/10.13039/501100001665

Description: https://hdl.handle.net/20.500.12085/c572cde8-a82c-4c2d-9bd7-288dfc8f1939

Description: https://www.aoml.noaa.gov/phod/gdp/data.php

Description: https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=GLOBAL_REANALYSIS_PHY_001_030

Description: https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=SEALEVEL_GLO_PHY_L4_REP_OBSERVATIONS_008_047

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.

Description: With the increase in computational power, ocean models with kilometer-scale resolution have emerged over the last decade. These models have been used for quantifying the energetic exchanges between spatial scales, informing the design of eddy parametrizations, and preparing observing networks. The increase in resolution, however, has drastically increased the size of model outputs, making it difficult to transfer and analyze the data. It remains, nonetheless, of primary importance to assess more systematically the realism of these models. Here, we showcase a cloud-based analysis framework proposed by the Pangeo project that aims to tackle such distribution and analysis challenges. We analyze the output of eight submesoscale-permitting simulations, all on the cloud, for a crossover region of the upcoming Surface Water and Ocean Topography (SWOT) altimeter mission near the Gulf Stream separation. The cloud-based analysis framework (i) minimizes the cost of duplicating and storing ghost copies of data and (ii) allows for seamless sharing of analysis results amongst collaborators. We describe the framework and provide example analyses (e.g., sea-surface height variability, submesoscale vertical buoyancy fluxes, and comparison to predictions from the mixed-layer instability parametrization). Basin- to global-scale, submesoscale-permitting models are still at their early stage of development; their cost and carbon footprints are also rather large. It would, therefore, benefit the community to document the different model configurations for future best practices. We also argue that an emphasis on data analysis strategies would be crucial for improving the models themselves.

Description: Agulhas leakage, the warm and salty inflow of Indian Ocean water into the Atlantic Ocean, is of importance for the climate-relevant Atlantic Meridional Overturning Circulation. South of Africa, the eastward turning Agulhas Current sheds Agulhas rings, cyclones and filaments of order 100 km that carry the Indian Ocean water into the Cape Basin and further into the Atlantic. Here, we show that the resolution of submesoscale flows of order 10 km in an ocean model leads to 40 % more Agulhas leakage and more realistic Cape Basin water-masses compared to a parallel non-submesoscale resolving simulation. Moreover, we show that submesoscale flows strengthen shear-edge eddies and in consequence lee cyclones at the northern edge of the Agulhas Current, as well as the leakage pathway in the region of the filaments that takes place outside of mesoscale eddies. This indicates that the increase in leakage can be attributed to stronger Agulhas filaments, when submesoscale flows are resolved.

Description: © The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Nature Communications 6 (2015): 6862, doi:10.1038/ncomms7862.

Description: Although the strongest ocean surface currents occur at horizontal scales of order 100 km, recent numerical simulations suggest that flows smaller than these mesoscale eddies can achieve important vertical transports in the upper ocean. These submesoscale flows, 1–100 km in horizontal extent, take heat and atmospheric gases down into the interior ocean, accelerating air–sea fluxes, and bring deep nutrients up into the sunlit surface layer, fueling primary production. Here we present observational evidence that submesoscale flows undergo a seasonal cycle in the surface mixed layer: they are much stronger in winter than in summer. Submesoscale flows are energized by baroclinic instabilities that develop around geostrophic eddies in the deep winter mixed layer at a horizontal scale of order 1–10 km. Flows larger than this instability scale are energized by turbulent scale interactions. Enhanced submesoscale activity in the winter mixed layer is expected to achieve efficient exchanges with the permanent thermocline below.

Description: J.C. and R.F. acknowledge financial support under grants ONR-N00014-09-1-0458 and NSF-OCE-1233832; J.M.K. under grants ONR-N00014-11-1-0165 and NSERC-327920-2006; J.G. under grant ONR-N00014-12-1-0105.

Description: With the increase in computational power, ocean models with kilometer-scale resolution have emerged over the last decade. These models have been used for quantifying the energetic exchanges between spatial scales, informing the design of eddy parametrizations, and preparing observing networks. The increase in resolution, however, has drastically increased the size of model outputs, making it difficult to transfer and analyze the data. It remains, nonetheless, of primary importance to assess more systematically the realism of these models. Here, we showcase a cloud-based analysis framework proposed by the Pangeo project that aims to tackle such distribution and analysis challenges. We analyze the output of eight submesoscale-permitting simulations, all on the cloud, for a crossover region of the upcoming Surface Water and Ocean Topography (SWOT) altimeter mission near the Gulf Stream separation. The cloud-based analysis framework (i) minimizes the cost of duplicating and storing ghost copies of data and (ii) allows for seamless sharing of analysis results amongst collaborators. We describe the framework and provide example analyses (e.g., sea-surface height variability, submesoscale vertical buoyancy fluxes, and comparison to predictions from the mixed-layer instability parametrization). Basin- to global-scale, submesoscale-permitting models are still at their early stage of development; their cost and carbon footprints are also rather large. It would, therefore, benefit the community to document the different model configurations for future best practices. We also argue that an emphasis on data analysis strategies would be crucial for improving the models themselves.

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Chelton, D. B., Schlax, M. G., Samelson, R. M., Farrar, J. T., Molemaker, M. J., McWilliams, J. C., & Gula, J. Prospects for future satellite estimation of small-scale variability of ocean surface velocity and vorticity. Progress in Oceanography, 173, (2019):256-350, doi:0.1016/j.pocean.2018.10.012.

Description: Recent technological developments have resulted in two techniques for estimating surface velocity with higher resolution than can be achieved from presently available nadir altimeter data: (1) Geostrophically computed estimates from high-resolution sea surface height (SSH) measured interferometrically by the wide-swath altimeter on the Surface Water and Ocean Topography (SWOT) Mission with a planned launch in 2021; and (2) Measurements of ocean surface velocity from a Doppler scatterometer mission that is in the early planning stages, referred to here as a Winds and Currents Mission (WaCM). In this study, we conduct an analysis of the effects of uncorrelated measurement errors and sampling errors on the errors of the measured and derived variables of interest (SSH and geostrophically computed velocity and vorticity for SWOT, and surface velocity and vorticity for WaCM). Our analysis includes derivations of analytical expressions for the variances and wavenumber spectra of the errors of the derived variables, which will be useful to other studies based on simulated SWOT and WaCM estimates of velocity and vorticity. We also discuss limitations of the geostrophic approximation that must be used for SWOT estimates of velocity. The errors of SWOT and WaCM estimates of velocity and vorticity at the full resolutions of the measured variables are too large for the unsmoothed estimates to be scientifically useful. It will be necessary to smooth the data to reduce the noise variance. We assess the resolution capabilities of smoothed estimates of velocity and vorticity from simulated noisy SWOT and WaCM data based on a high-resolution model of the California Current System (CCS). By our suggested minimum threshold signal-to-noise (S/N) variance ratio of 10 (a standard deviation ratio of 3.16), we conclude that the wavelength resolution capabilities of maps of velocity and vorticity constructed from WaCM data with a swath width of 1200 km are, respectively, about 60 km and 90 km in 4-day averages. For context, the radii of resolvable features are about four times smaller than these mesoscale wavelength resolutions. If the swath width can be increased to 1800 km, the wavelength resolution capabilities of 4-day average maps of surface velocity and vorticity would improve to about 45 km and 70 km, respectively. Reducing the standard deviation of the uncorrelated measurement errors from the baseline value of m s−1 to a value of 0.25 m s−1 would further improve these resolution capabilities to about 20 km and 45 km. SWOT data will allow mapping of the SSH field with far greater accuracy and space–time resolution than are presently achieved by merging the data from multiple nadir altimeter missions. However, because of its much narrower 120-km measurement swath compared with WaCM and the nature of the space–time evolution of the sampling pattern during each 21-day repeat of the SWOT orbit, maps of geostrophically computed velocity and vorticity averaged over the 14-day period that is required for SWOT to observe the full CCS model domain are contaminated by sampling errors that are too large for the estimates to be useful for any amount of smoothing considered here. Reducing the SSH measurement errors would do little to improve SWOT maps of velocity and vorticity. SWOT estimates of these variables are likely to be useful only within individual measurement swaths or with the help of dynamic interpolation from a data assimilation model. By our criterion, in-swath SWOT estimates of velocity and vorticity have wavelength resolution capabilities of about 30 km and 55 km, respectively. In comparison, in-swath estimates of velocity and vorticity from WaCM data with m s−1 have a wavelength resolution capability of about 130 km for both variables. Reducing the WaCM measurement errors to m s−1 would improve the resolution capabilities to about 50 km and 75 km for velocity and vorticity, respectively. These resolutions are somewhat coarser than the in-swath estimates from SWOT data, but the swath width is more than an order of magnitude wider for WaCM. Instantaneous maps of velocity and vorticity constructed in-swath from WaCM data will therefore be much less prone to edge effect problems in the spatially smoothed fields. Depending on the precise value of the threshold adopted for the minimum S/N ratio and on the details of the filter used to smooth the SWOT and WaCM data, the resolution capabilities summarized above may be somewhat pessimistic. On the other hand, aspects of measurement errors and sampling errors that have not been accounted for in this study will worsen the resolution capabilities presented here. Another caveat to keep in mind is that the resolution capabilities deduced here from simulations of the CCS region during summertime may differ somewhat at other times of year and in other geographical regions where the signal variances and wavenumber spectra of the variables of interest differ from the CCS model used in this study. Our analysis nonetheless provides useful guidelines for the resolutions that can be expected from SWOT and WaCM.

Description: We thank Ralph Milliff, Bo Qiu, Ernesto Rodríguez and Lee-Lueng Fu for many helpful editorial comments and suggestions that improved the manuscript. This research was funded by NASA Grants NNX13AD78G, NNX14AM72G, NNX13AE32G, NNX14AM66G, NNX16AH76G,NNX14AM71G and NNX17AH54G. The two North Atlantic Ocean simulations in this study were performed using HPC resources from GENCI-TGCC with support from Grant 2017-A0010107638 for Jonathan Gula.

Description: Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 46(5), (2019):2704-2714, doi:10.1029/2019GL081919.

Description: Seismic images and glider sections of the Gulf Stream front along the U.S. eastern seaboard capture deep, lens‐shaped submesoscale features. These features have radii of 5–20 km, thicknesses of 150–300 m, and are located at depths greater than 500 m. These are typical signatures of anticyclonic submesoscale coherent vortices. A submesoscale‐resolving realistic simulation, which reproduces submesoscale coherent vortices with the same characteristics, is used to analyze their generation mechanism. Submesoscale coherent vortices are primarily generated where the Gulf Stream meets the Charleston Bump, a deep topographic feature, due to the frictional effects and intense mixing in the wake of the topography. These submesoscale coherent vortices can transport waters from the Charleston Bump's thick bottom mixed layer over long distances and spread them within the subtropical gyre. Their net effect on heat and salt distribution remains to be quantified.

Description: J. G. gratefully acknowledges support from the French government, managed by the French National Agency for Research (ANR), through programs ISblue (ANR‐17‐EURE‐0015) and LabexMER (ANR‐10‐LABX‐19) and from LEFE/IMAGO through the project AO2017‐994457‐RADII. Simulations were performed using HPC resources from GENCI‐TGCC (grant 2017‐A0010107638). Simulations output is available upon request. Seismic data were processed using the Echos software package by Paradigm, Matlab, and Generic Mapping Tools. The Eastern North America Margin Community Seismic Experiment was funded by the National Science Foundation under grant OCE‐1347498 and UNOLS; cruise data are freely available via the Marine Geoscience Data System Academic Seismic Portal at Lamont‐Doherty Earth Observatory (http://www.marine-geo.org/portals/seismic/). Spray glider observations in the Gulf Stream are available from http://spraydata.ucsd.edu and should be cited using the following DOI (10.21238/S8SPRAY2675; Todd & Owens, 2016). Spray glider operations were funded by the National Science Foundation (OCE‐1633911) and the Office of Naval Research (N000141713040).

Description: 2019-08-27