# Ocean 2D eddy energy fluxes from small mesoscale processes with SWOT

Elisa Carli, Rosemary Morrow, Oscar Vergara, Robin Chevrier, Lionel Renault

► To cite this version:

Elisa Carli, Rosemary Morrow, Oscar Vergara, Robin Chevrier, Lionel Renault. Ocean 2D eddy energy fluxes from small mesoscale processes with SWOT. Ocean Science, 2023, 19 (5), pp.1413-1435. 10.5194/os-19-1413-2023. hal-04501353

HAL Id: hal-04501353

https://hal.science/hal-04501353v1

Submitted on 14 Mar 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Ocean 2D eddy energy fluxes from small mesoscale processes with SWOT

Elisa Carli<fnref n="1" />, Rosemary Morrow<fnref n="1" />, Oscar Vergara<fnref n="2" />, Robin Chevrier<fnref n="2" />, and Lionel Renault<fnref n="1" />

**Correspondence:** Elisa Carli (elisa.carli@univ-tlse3.fr)

Received: 27 May 2023 – Discussion started: 1 June 2023

Revised: 11 August 2023 – Accepted: 16 August 2023 – Published: 6 October 2023

**Abstract.** We investigate ocean dynamics at different scales in the Agulhas Current system, a region of important interocean exchange of heat and energy. While ocean observations and some of the most advanced climate models capture the larger mesoscale dynamics (> 100 km), the smaller-scale fronts and eddies are underrepresented. The recently launched NASA–CNES Surface Water and Ocean Topography (SWOT) wide-swath altimeter mission observes the smaller ocean geostrophic scales down to 15 km in wavelength globally. Here we will analyse different eddy diagnostics in the Agulhas Current region and quantify the contributions from the larger mesoscales observable today and the smaller scales to be observed with SWOT. Surface geostrophic diagnostics of eddy kinetic energy, strain, and energy cascades are estimated from modelled sea surface height (SSH) fields of the Massachusetts Institute of Technology general circulation model (MITgcm) latitude–longitude polar cap (LLC4320) simulation subsampled at 1/10°. In this region, the smaller scales (< 150 km) have a strong signature on the horizontal geostrophic strain rate and for all eddy diagnostics in the Western Boundary Current and along the meandering Agulhas Extension. We investigate the horizontal cascade of energy using a coarse-graining technique, and we observe that the wavelength range where the inverse cascade occurs is biased towards larger mesoscale wavelengths with today’s altimetric sampling. We also calculate the projected sampling of the eddy diagnostics under the SWOT swaths built with the NASA–CNES simulator to include the satellite position and realistic noise. For the swaths, a neural network noise mitigation method is implemented to reduce the residual SWOT random error before calculating eddy diagnostics. In terms of SSH, observable wavelengths of 15 to 20 km are

retrieved after neural network noise mitigation, as opposed to wavelengths larger than 40 km before the noise reduction.

## 1 Introduction

The Southern Ocean (SO) is a highly energetic region with an important role in our global climate, as it drives a large meridional overturning circulation and is a key region of the formation and modification of mode, intermediate, and deep waters (Bourgeois et al., 2022). The cold and oxygen-depleted surface water masses allow efficient storage of the anthropogenic heat and carbon, making the SO the greatest sink for these tracers (45%–62% of global ocean heat gain from 2005 to 2017) (IPCC, 2019). This SO circulation includes large-scale mean and highly turbulent transient flows due to eddies, fronts, and meanders (Morrow and Le Traon, 2012; Sallée, 2018), which is the most energetic component of the global ocean circulation (Wunsch, 2002; Chelton et al., 2007). The energetic eddy component of the SO circulation is key in transporting heat polewards across the zonally circulating Circumpolar Current (Volkov et al., 2010). Horizontal and vertical fluxes induced by eddies and meanders contribute to the redistribution of energy, momentum, and tracers such as heat, carbon, and nutrients (Thomas et al., 2008; Lévy et al., 2012; Sallée, 2018). SO eddy kinetic energy varies regionally in relation to the mean currents and orography but also evolves interannually in response to climate modes and wind forcing, with eddies rec>tifying the mean flow (Meredith and Hogg, 2006; Morrow et al., 2010; Sinha and Abernathey, 2016; Martínez-Moreno et al., 2019). Areas of strong eddy surface strain are also

<footnote n="1">¹Laboratoire d'Etudes en Géophysique et Océanographie Spatiales (CNES-CNRS-IRD-UPS), Toulouse, France</footnote>

<footnote n="2">²Collecte Localisation Satellites (CLS), Toulouse, France</footnote>
hotspots of smaller-scale (< 50 km) ageostrophic ocean vertical movements in the SO (Su et al., 2018; Siegelman et al., 2020), with geostrophic surface strain impacting the vertical distribution of heat, carbon, and nutrients (Zhang et al., 2019).

Within the SO, the Agulhas Current system, flowing south-
westwards along Africa's continental shelf break, plays a key
role in this global ocean circulation, feeding the upper branch
of the Atlantic Meridional Overturning Circulation (AMOC).
From the southern tip of the continental shelf (20° E), the cur-
rent presents an abrupt curvature back into the southern In-
dian Ocean known as the Agulhas Retroflection (Lutjeharms,
2006). Large anticyclonic eddies known as the Agulhas rings
are shed from the retroflection into the eastern boundary of
the South Atlantic and contribute to the inter-basin exchange
by leaking saline warm waters from the Indian Ocean to the
Atlantic (Olson and Evans, 1986; Gordon et al., 1992; Lutje-
harms, 2006; Beal et al., 2011). The Agulhas Current system
is one of the most energetic regions of the global oceans and
presents a series of large-scale (> 150 km) eddies and mean-
ders, thought to be the main driving mechanisms of the cur-
rent until a few decades ago (Boyd, 1994). This mesoscale
variability plays a strong role in the region's intense air-sea
interactions (Speich et al., 2007; Swart et al., 2008; Renault
et al., 2017), interocean exchange, and water mass modifi-
cation (Lutjeharms and Gordon, 1987; Gordon et al., 1992;
Boyd, 1994; Speich et al., 2007).

Most of these studies are based on eddy-resolving simu-
lations (1/12° or 1/10° spatial resolution) or global altimetry
maps that resolve the larger mesoscales > 100 km in diame-
ter but not the smaller mesoscales or submesoscales, yet the
ocean fine-scale dynamics (15 to 150 km), with their strong
gradients in ocean properties and temporal scales ranging
from days to weeks, are now understood to affect the ocean
physics and biomass up to climate scales (Ferrari and Wun-
sch, 2008). Over the last decade, high-resolution ocean mod-
els have made great advances in resolving finer-spatial- and
finer-temporal-scale dynamics, down to a few kilometres re-
gionally (Renault et al., 2018; Verron et al., 2020; Contr-
eras et al., 2023; Bendinger et al., 2023) and globally (Arbic
et al., 2022). These models provide a more detailed descrip-
tion of the dynamics of energetic Western Boundary Cur-
rent systems and in particular the role of meso- and subme-
soscale structures and their rectification on the mean circu-
lation (McWilliams, 2008; Gula et al., 2014; Renault et al.,
2016; Chassignet and Xu, 2017; Renault et al., 2017; Contr-
eras et al., 2023). Within the Agulhas Current system, a large
number of high-resolution modelling studies have addressed
the role of eddy-current interactions in the generation of the
large Natal Pulses (Krug and Penven, 2011), coastal Kelvin
waves (Sebille and Leeuwen, 2007), or the submesoscale fast
barotropic motions (Tedesco et al., 2019; Schubert et al.,
2020).

Although observational efforts such as the Argo pro-
gramme and remote sensing programmes provide global ob-

servations of the large-scale ocean circulation and larger
mesoscale ocean dynamics, a gap concerning global obser-
vations of the finer scales (< 150 km in wavelength) exists,
which could lead to misinterpretation or loss of key phys-
ical or biophysical mechanisms in ocean models (d'Ovidio
et al., 2019). Local in situ campaigns can target individ-
ual, small-scale features in 3D, but they only represent a
fraction of the ocean conditions, with limited time and spa-
tial coverage. Satellite sea surface temperature (SST) and
ocean colour observe these finer scales in non-cloudy con-
ditions, which limits their applicability over the SO. Along-
track satellite altimetry can resolve local spatial scales down
to 30 to 70 km when averaged through wavenumber spec-
tral analyses (Arbic et al., 2013; Dufau et al., 2016; Renault
et al., 2019; Vergara et al., 2023). Yet, a significant part of
the observed sea surface height (SSH) variability originates
at scales shorter than the 150 km in wavelength observed
by altimetric maps today (Sasaki et al., 2014; Callies et al.,
2015). Due to this observational gap, it is difficult to accu-
rately observe the ocean's energy cascades. We are missing
observations of how and where the small scales interact with
the large-scale field, in order to understand where the ener-
getic small scales compensate for or enhance the larger-scale
flow. Indeed, larger mesoscale eddies often appear and disap-
pear in the mid-ocean (Chelton et al., 2011) since we are not
correctly observing their key small-scale generation and dis-
sipation processes. Altimetric maps merging multi-mission
information also underestimate the eddy energy in energetic
boundary current regions due to the trade-off between the
spatial resolution and the coverage of the mapping technique
that smoothes out features smaller than 150 km. These are
all open questions but need more fine-scale observations to
close the global ocean energy budget (Ferrari and Wunsch,
2008). The recently launched NASA-CNES Surface Water
and Ocean Topography (SWOT) mission is expected to par-
tially fill this gap by increasing the current observational ca-
pabilities of satellite altimetry down to 15 km wavelength
(Desai, 2018; Morrow et al., 2019). The improved interfero-
metric synthetic aperture radar (SAR) technology will allow
us to observe the fine-scale sea surface height at 2 km res-
olution, a signature that responds to dynamics driven from
within and below the surface mixed layer. SWOT is ex-
pected to provide an unprecedented dataset to validate the
geostrophic part of eddy-resolving ocean models down to
scales of 15 km in wavelength (i.e. equivalent eddy diame-
ters of 8 km, Rossby radii of 4 km).

As for all satellite observations, the observing capabilities
of SWOT will be impacted by different sources of measure-
ment error due to instrument and atmospheric effects and
also the geophysical characteristics of the ocean circulation
(waves, tides, etc.) (Dibarboure et al., 2022). SWOT was
successfully launched in December 2022, and the first data
should be available in October 2023. In this paper, our ob-
jective is to diagnose SWOT observations in terms of sig-
nal and noise in order to better understand the mesoscale
dynamics observed today with multi-mission maps from the Data Unification and Altimeter Combination System and the Copernicus Marine Environment Monitoring Service (DU-ACS/CMEMS) in the Agulhas Current system and to investigate the new dynamical scales to be revealed by SWOT. To reproduce Surface Water and Ocean Topography (SWOT) we use a realistic global high-resolution model including tidal forcing, based on the Massachusetts Institute of Technology general circulation model (MITgcm) latitude-longitude polar cap (LLC4320) simulation (Rocha et al., 2015). Modelled SSH fields and their geostrophic approximation for surface currents and surface eddy diagnostics are analysed in the Agulhas Current region. We also perform equivalent nadir altimetric sampling of our model SSH fields and then reconstruct 2D “pseudo-DUACS” maps to quantify the scales observable with today’s mapping and interpolation schemes. Various eddy diagnostics are calculated on both datasets to compare the ocean dynamics resolved in today’s observations with those to be resolved by SWOT. These include the smaller scales < 150 km in wavelength and the important cross-terms within the energy cascade linking the smaller- and larger-scale dynamics. We then diagnose the impact of the instrumental and geophysical errors on the SSH observations by applying the SWOT ocean simulator (Gaultier et al., 2016). An innovative artificial intelligence (AI) processing (Tréboutte et al., 2023) has been developed and tested in the North Atlantic Ocean to reduce the impact of random instrumental noise in SWOT observations and is planned to be implemented in the global SWOT operational environment. We will test its performance in reducing the SWOT random noise in the Agulhas Current system and assess the adaptability of this AI method on a different model and different zone with respect to its training. Finally, we can estimate the scales of observability of the Agulhas Current system’s ocean dynamics in an ideal error-free scenario for noisy SWOT data and after noise reduction.

The paper is organized as follows. The data and methods are introduced in Sect. 2. Section 3 presents the results of the eddy diagnostics over the full Agulhas Current system model domain. In Sect. 4 we present the results of the eddy diagnostics extracted under the SWOT swaths and compare the observability of the ocean eddy diagnostics before and after the AI treatment of the residual KarIn random error. Finally, a discussion and a conclusion are provided in Sect. 5.

## 2 Methodology and data

This section describes the Massachusetts Institute of Technology general circulation model (MITgcm), which is the mesoscale resolving model used for the analysis, the reconstructed pseudo-DUACS data used as a reference for today’s altimetry measurements capability, and the SWOT simulator used to reconstruct SWOT orbits and simulate the observed position and errors.

### 2.1 LLC4320 simulation

The model data used in this study to represent the total ocean signal are the outputs of the mesoscale and submesoscale resolving MITgcm model (Marshall et al., 1997). The SWOT community has largely used the specific simulation LLC4320, where “LLC” refers to a latitude-longitude polar cap global model grid (Forget et al., 2015), which in its higher resolution has a nominal horizontal grid spacing of 1/48° globally, similar to the SWOT swath grid spacing of 2 km × 2 km, with an effective resolution of less than 20 km (Rocha et al., 2016). Outputs are hourly snapshots that span the period from 13 September 2011 to 15 November 2012 (Lin et al., 2020). There are 90 vertical levels ranging in thickness from 1 m at the surface to 480 m at a maximum model depth of 7000 m. The LLC4320 is forced at the surface with 6-hourly 0.14° atmospheric fields from the European Centre for Medium-Range Weather Forecasts (ECMWF) and by the 16 most significant components of the hourly tidal potential, applied as additional surface pressure, generating both barotropic and baroclinic tidal components (Chaudhuri et al., 2013). A few details on the validation of the model for the specific use of this study are given in Appendix A. Further information about the LLC4320, its construction, and its validation can be found in Torres et al. (2018).

This work focuses on sea surface height and geostrophic currents that would be observed with satellite altimetry in the Agulhas region, represented in Fig. 1, and our diagnostics are based on the modelled SSH. To reproduce an equivalent altimetry-like SSH field, the model’s barotropic tides and the high-frequency barotropic response to the atmospheric forcing (the so-called dynamical atmospheric correction: DAC) are estimated and removed. The 1/48° simulation is known to have minor forcing errors in the tides and an offset on the atmosphere forcing fields (by 6 h) (Arbic et al., 2022). This results in barotropic and baroclinic tidal artefacts and wind-forcing offsets on the SSH signal. Unrealistic grid-type features at the spatial scale of the atmospheric forcing are also observed in the 1/48° SSH field, creating discontinuities when deriving eddy diagnostics from the geostrophic components of the velocity field.

The data used in this work are the outputs of the original high-resolution run of the MITgcm LLC4320, smoothed on a regular global grid at 1/20° and subsampled at 1/10°. This allows us to greatly reduce the discontinuities previously mentioned. Since we do not aim to analyse the ocean geostrophic dynamics at scales smaller than the Surface Water and Ocean Topography (SWOT) resolution (around 15 km at best), we will rely on the 1/10° subsampled modelled simulations for our study. This 1/10° version was also used to calculate the pseudo-DUACS dataset and for simulating the global SWOT cross-calibration. This dataset is referred to as LLC10 in the following.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Superposition of 2D geographic map and swath total geostrophic current on 1st January 2012 in the Agulhas Current region. The dashed black line defines the zone of this study. The yellow and purple contours define the 1500 and 3000 m topography respectively. Data on the swath include instrumental and geophysical measurement errors. In the fast-sampling phase of the mission, four tracks will sample the region: 20, 5, 3, and 18.</figcaption></figure>

## 2.2 DUACS data and pseudo-DUACS reconstructed dataset from the LLC10 simulation

Today's multi-mission Data Unification and Altimeter Combination System (DUACS) altimetry products deliver global and regional sea level and geostrophic current maps for oceanographic applications (distributed by CMEMS). The latest version, the level-4 DUACS-DT (Delayed Time), is constructed by optimal interpolation of level-3 altimeter observations (i.e. with along-track corrections, editing, and processing applied) and is produced on a regular 1/4° × 1/4° grid with daily sampling (Taburet et al., 2019). However, Chelton et al. (2011) found that the capability of retrieving small mesoscale dynamics in level-4 data is limited by the optimal interpolation framework used in DUACS. Ballarotta et al. (2019) analysed the effective spatial and temporal resolution of level-4 sea level anomaly (SLA) globally. For the time span of the LLC4320 simulation, the satellites in use for the generation of the DUACS product are Jason-2, ENVISAT, HY-2A, and CryoSat-2 (Taburet et al., 2019). Further details on this observational product and its benchmarking against its surrogate derived from LLC10 can be found in Appendix A.

If we assume that our modelled SSH represents reality, our first step is to verify how this SSH would be observed by today's along-track altimetric sampling and mapping, so we built a pseudo-DUACS product. This product consists of the same DUACS mapping technique described above; however, the data are subsampled from the LLC10 fields along the altimeter mission tracks. These altimeter-like model data are then processed using the DUACS optimal interpolation

mapping algorithm, applying the same mapping parameters used for the real DUACS maps, thereby providing pseudo-DUACS daily maps on a regular 1/10° grid. Since the original model has no data assimilation, its eddy field may be shifted compared to the real eddy position, and the advantage of the pseudo-DUACS is to have better colocation of the modelled large- and small-scale features. So in this study, the pseudo-DUACS product is used as a proxy for the current observational capability of satellite altimetry in our region, in terms of spatial and temporal resolution. In the region of the Agulhas Current, the effective resolution of the pseudo-DUACS gridded maps is larger than 150 km and represents dynamics with temporal variability greater than 10 d. SWOT will potentially observe all scales included in the LLC10 SSH fields. In the computation, the smaller-scale processes observable with SWOT have been defined as the residual between the total LLC10 fields (corrected for the high-frequency motions) and the pseudo-DUACS fields. This derived product, referred to as "SSH residuals", is analysed to investigate the new dynamics that SWOT will observe.

Figure 2 illustrates the SSH spectral power of the different products derived from the LLC10 in terms of time and spatial frequencies, calculated over the annual model period. The top panels show the LLC10 original SSH (a) and the same field (b) after removing the high-frequency correction for the barotropic tide and DAC. As expected, most of the higher-frequency tidal energy (6 h, 12 h, 1 d) disappears in the corrected version of the spectrum. There are weaker residual energy signals at tidal frequencies that may be due to errors in the barotropic tide correction or at smaller scales, internal tides. Although the MITgcm LLC4320 is known to overestimate the tide forcing by a factor of 1.1121 (Arbic et al., 2022), in our Agulhas region, the residual internal tide signal is very low compared to the energetic mesoscale SSH field shown in Fig. 2b. The DAC correction removes much of the large-scale, high-frequency energy at timescales < 10 d; since this has a large-scale structure, it does not impact the eddy diagnostics we derive that are shaped by the smaller-scale SSH gradients. The SSH spectral power for the reconstructed pseudo-DUACS product is shown in Fig. 2c and for the SSH residual high frequencies in Fig. 2d. The time frequencies of the pseudo-DUACS product and the SSH residuals are cut at the Nyquist frequency of 2 d because of their daily time sampling. We note that even if the pseudo-DUACS mapping decorrelation scales are around 150 km in the Agulhas region (Ballarotta et al., 2019), its frequency-wavenumber spectra have some residual power in scales smaller than the effective resolution of 150 km in Fig. 2c and smaller lobes of energy around 30 km in wavelength, whereas the SSH residuals (d) represent the residual SSH power (b–c) at scales < 150 km and between 2–10 d well.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** SSH wavenumber–frequency energy spectra. (a) Total SSH from the LLC10 simulation, (b) SSH from LLC10 corrected for tide and DAC, (c) large-scale pseudo-DUACS data in the LLC10 simulation, and (d) residual from (b)–(c) retaining the high-frequency and smaller-scale SSH.</figcaption></figure>

## 2.3 The SWOT mission and the SWOT ocean simulator

The Surface Water Ocean Topography (SWOT) mission, launched on 16 December 2022, will provide the first global SAR interferometry observations providing collocated 2D surface topography and 250 m resolution SAR images (Fu et al., 2012; Morrow et al., 2019). SWOT observations are made over two 50 km wide swaths, with a Jason-class nadir altimeter in the centre. An onboard processor allows for a 250 m × 250 m expert product and a precise 2 km × 2 km grid product over all oceans. The measurement is very precise, with an expected root mean square (rms) SSH noise at 2 km of 1.37 cm, an order of magnitude smaller than conventional Jason-class altimeters (Chelton, 2019). SWOT will have two orbit phases: a daily revisit over a limited number of tracks during the first 6 months (January to June 2023) of the mission’s commissioning and a calibration and validation phase (the CalVal orbit). From July 2023, the satellite had moved into the science orbit with a 21 d revisit and global coverage. Both orbits cover latitudes of up to 78° (Dibaboure et al., 2022). SWOT is the first altimeter mission whose science requirements are specified in terms of wavenumber spectra and observable wavelengths; these requirements are that the SSH signal should exceed the noise (for 2 m wave height conditions) at 15 km wavelength globally (Desai, 2018). During the daily CalVal phase of the SWOT mission, a series of international in situ validation campaigns are proposed (https://www.swot-adac.org/, last access: 16 August 2023). One of these campaigns, QUICCHE (Quantifying Interocean fluxes in the Cape Cauldron Hotspot of

Eddy kinetic energy), targets the Agulhas Cape Cauldron region (https://www.swot-adac.org/campaigns/quicche/, last access: 16 August 2023) where a SWOT cross-over exists (SWOT tracks 5 and 20; Fig. 1). This region is used as a test site in our study to evaluate the projected SWOT signal-to-noise ratio.

### 2.3.1 SWOT simulator

The SWOT simulator provides the swath position and instrument errors along SWOT tracks, for both the daily fast sampling and the 21 d repeat science orbits as described by Gaultier et al. (2016). Here the main characteristics are given; for more in-depth information on the simulator refer to https://swot-simulator.readthedocs.io/en/latest/ (last access: 3 May 2023). The SWOT simulator allows the extraction of SSH data from an ocean general circulation model and the interpolation of the latter onto the SWOT swath grid. For our study, the signal is derived from the LLC10. The simulator interpolates the original model grid onto the along-track-cross-track SWOT 2 km grid, the SSH interpolation is linear in space, and no interpolation is performed in time. The results presented in this study are from the CalVal phase, so the outputs are daily cycles on the fast-sampling orbit.

The SWOT simulator then allows us to generate most sources of error for the SWOT measurement, as foreseen in the SWOT error budget documentation (Esteban-Fernandez, 2013) and described in Gaultier et al. (2016) and Rodriguez et al. (2017). The main sources will be instrumental: random noise (which varies with the sea state and other inho-
mogeneities) generated by KarIn, the Ka-band radar interferometer, the main instrument on board SWOT, and systematic errors (satellite roll, phase errors, baseline dilation, timing errors, orbital errors). There are also contributions from geophysical errors at mesoscales (mainly wet troposphere). The SWOT simulator uses a global spectral wavenumber estimate of each error and projects that back locally onto the swaths in space and time. Figure 1 is an example of the four daily fast-sampling orbits that pass in the Agulhas Current region, showing a snapshot of the total geostrophic current computed with the SWOT simulator data (with the full spectrum of errors added to the data) over the LLC10 simulation. It is clear that the full spectrum of the SSH error and their gradients makes it impossible to interpret the underlying geostrophic current field and that the finer noise structures could be misinterpreted as mesoscale structures, as already observed by Fu and Ubelmann (2014) and Chelton et al. (2022). Dibarboure et al. (2022) estimate that systematic errors alone contribute tens of centimetres in SSH. However the SWOT project’s cross-calibration techniques are efficient in reducing these simulated systematic errors, and since their wavelengths are larger than the mesoscales observed in this paper, the impact on our eddy statistics is minimal.

This study focuses on the influence of the KarIn random error on the SSH data and observable wavelengths in sea surface height (SSH), eddy kinetic energy (EKE), and strain rate. The KarIn random error is defined as a Gaussian noise with a zero-centred distribution and is strongly dependent on the sea state. As shown in Esteban-Fernandez (2013), its standard deviation depends on the distance to the nadir and on the SWH. Since a realistic significant wave height (SWH) is not available in the LLC4320 simulation, the KarIn random error for our analysis has been generated from the spectral SWOT requirement specification with a constant SWH of 2 m.

## 2.3.2 Mitigation of the KarIn random error

The step of de-noising along-track or swath altimetry data is fundamental for exploiting the real altimeter signal. In this work, we only focus on the reduction in the random error associated with the KarIn instrument. The SWOT project will perform cross-calibration to reduce the larger-scale errors (Dibarboure et al., 2022). Wang et al. (2019) compute the wavelengths observable with SWOT when the residual random error is added to the SSH signal, based on a realistic wave field and no noise reduction strategy. At higher latitudes with higher wave fields such as in the Southern Ocean, they find that without de-noising, wavelengths smaller than around 40 km will be hidden by the random noise. The observable wavelength increases as higher-order derivatives are computed on the SSH, such as geostrophic velocities (first-order derivative) or diagnostics like the strain rate (second-order derivative). To address this problem of the random noise in a 2D field such as SWOT, several methods have

been explored, such as the spatial-based median and Lanczos filters (Fan et al., 2019) or the boxcar, Gaussian, or Laplacian filters. Gómez-Navarro et al. (2020) and Gómez-Navarro et al. (2018) show that these spatial filters do not give satisfying results. In Gómez-Navarro et al. (2018), a variational filter was specifically created to treat SWOT data. The objective was to minimize a cost function by optimizing a parameter called $\lambda_2$ which depends on the study area and the season and maintains smooth second-order derivatives (e.g. vorticity). Tréboutte et al. (2023) developed an alternative method based on a convolutional neural network (CNN), called U-Net, which is more efficient than previous filtering techniques and preserves small-scale ocean features, including those near coasts and islands. In this work, we evaluate the performance of the U-Net noise mitigation technique to reduce the random noise in our pseudo-SWOT observations.

The de-noising model is based on a U-Net architecture, which is trained and tested on the 5 years of the eNATL60 simulation (https://zenodo.org/record/4032732, last access: 15 May 2023) in the North Atlantic region. The basic training of U-Net is performed by comparing the de-noised SSH to the true reference SSH using a loss function (Tréboutte et al., 2023). If the residual is not improved for a chosen number of steps N, the training stops to avoid overfitting. For information on the noise reduction technique, the reader is referred to Tréboutte et al. (2023). The robustness of the technique has been tested for different scenarios, including a de-noising of the global 1/12° reanalysis from Mercator Ocean International, GLORYS, based on training only in the North Atlantic with the eNATL60 model.

Our objective is to test the U-Net de-noising algorithm in the region of the Agulhas Current with the LLC10 model outputs along SWOT tracks generated using the SWOT simulator. This dataset represents different dynamics from the North Atlantic training zone, being in a different and very energetic zone, with a different model. The U-Net algorithm has not been retrained, and here we test its capabilities of adaptability to a different, energetic region and a different model in order to estimate if this method will be applicable to the early real SWOT data.

The scores put in place to assess the robustness of the noise reduction technique are the following: the root mean square error (RMSE), the variance of the SSH residuals, and the wavelength observable over the swaths. We compute the observability following Wang et al. (2019) as the ratio between the spectral content of the noisy signal ($h_{noisy}$) and the spectral content of the reference SSH ($h_{true}$), simulated via the SWOT simulator without added noise. In this work, the noisy signal is either the LLC10 SSH with the total added KarIn random error or the LLC10 SSH with the residual noise after treatment with U-Net. The observable wavelength is the wavelength for which the NSR, defined in Eq. (1), reaches 1. The two results are compared to understand the influence of
the U-Net de-noising on SWOT observability.

$$NSR = \frac{\text{PSD}(h_{\text{noisy}} - h_{\text{true}})}{\text{PSD}(h_{\text{true}})} = 1, \qquad (1)$$

where the PSD is the power spectral density of the signal.

### 3 Diagnostics on the 2D maps

This section aims at showing the potential of the SWOT mission by analysing eddy diagnostics based on the model-simulated SSH alone, with no SWOT errors applied. We compare the results for the larger mesoscales (> 150 km) observed today, represented by the pseudo-DUACS product, and the smaller and faster dynamics that should be newly observed by SWOT, represented by the residuals of the LLC10 dataset after removing the pseudo-DUACS. Although SWOT aims to observe the full spectra of larger and smaller geostrophic mesoscale structures and their interactions, we focus on the residuals and the smaller scales.

The main geostrophic eddy diagnostics investigated in this work are described below.

- Relative vorticity. This diagnoses the SSH field variability in terms of its small-scale and turbulence content. We compute the normalized form as

$$\zeta = \left( \frac{\partial v_g}{\partial x} - \frac{\partial u_g}{\partial y} \right) / f, \qquad (2)$$

where $u_g$ and $v_g$ are the surface geostrophic velocities derived from the altimetric SSH ($h$):

$$\begin{aligned} u_g &= -\frac{g}{f} \frac{\partial h}{\partial y} \\ v_g &= \frac{g}{f} \frac{\partial h}{\partial x}, \end{aligned} \qquad (3)$$

where $g$ is gravitational acceleration and $f = 2\Omega\sin\theta$ is the Coriolis frequency, proportional to $\Omega$, the Earth's rotation rate, and to $\theta$, the latitude.

- EKE. Eddy variability is identified by EKE, and altimetric EKE is often used to validate the realism of ocean models. This component of the energy is directly related to the temporal evolution of the fluid parcels. It is computed from the anomalies of the zonal and meridional components of the geostrophic velocities. The anomalies are defined as the velocity minus its temporal mean at each grid point.

$$\text{EKE} = \frac{1}{2} (u_g^{\prime 2} + v_g^{\prime 2}) \qquad (4)$$

We note that when separating our EKE analysis into larger scales represented by CMEMS mapping and residual smaller scales, this EKE equation needs to be

expanded, and important cross-terms that mix the large and small scales result. The result is in terms of the geostrophic velocity anomaly of the total LLC10 simulation $((u, v)_{gLLC10})$ and that of the pseudo-DUACS data $((u, v)_{gPD})$ where the subscript $g$ for geostrophy and the ' for anomaly are omitted for simplicity.

We first calculate the EKE from the corrected full model ($EKE_{LLC10} = 1/2(u_{LLC10}^2 + v_{LLC10}^2)$) and from the pseudo-DUACS reconstruction ($EKE_{PD} = 1/2(u_{PD}^2 + v_{PD}^2)$). Their difference contains the contribution of the small scales only ($EKE_{SS} = EKE_{LLC10} - EKE_{PD}$). The EKE in Eq. (5), computed from the residual small-scale geostrophic velocities $((u, v)_{\text{res}} = (u, v)_{LLC10} - (u, v)_{\text{PD}})$, contains both the contribution of the small scales and the cross-terms:

$$\begin{aligned} EKE_{\text{res}} &= \frac{1}{2}(u_{\text{LLC10}} - u_{\text{PD}})^2 + \frac{1}{2}(v_{\text{LLC10}} - v_{\text{PD}})^2 \\ &= \frac{1}{2} \left( u_{\text{LLC10}}^2 + u_{\text{PD}}^2 - 2u_{\text{LLC10}}u_{\text{PD}} \right. \\ & \qquad \left. + v_{\text{LLC10}}^2 + v_{\text{PD}}^2 - 2v_{\text{LLC10}}v_{\text{PD}} \right) \\ &= \frac{1}{2} \left( u_{\text{LLC10}}^2 + v_{\text{LLC10}}^2 \right) + \frac{1}{2} \left( u_{\text{PD}}^2 + v_{\text{PD}}^2 \right) \\ & \qquad - u_{\text{LLC10}}u_{\text{PD}} - v_{\text{LLC10}}v_{\text{PD}}. \end{aligned} \qquad (5)$$

The cross-term EKE in Eq. (6) can be found by subtracting the small-scale EKE from the residual EKE:

$$\begin{aligned} EKE_{\text{ct}} = EKE_{\text{res}} - EKE_{\text{SS}} &= u_{\text{PD}}^2 + v_{\text{PD}}^2 \\ & \quad - u_{\text{LLC10}}u_{\text{PD}} - v_{\text{LLC10}}v_{\text{PD}}. \end{aligned} \qquad (6)$$

- Strain rate. This is an estimate of the surface deformation of the current. At the surface, regions of high strain are associated with a secondary circulation characterized by high vertical velocities. These areas can be key for the vertical exchange of heat, carbon, and nutrients and are biologically more productive (Zhang et al., 2019). Its formulation is defined from surface geostrophic velocities. As with the EKE analysis, an expanded version of the strain rate is needed when calculating the difference between large CMEMS scales and residual small scales, which also include cross-terms. This component is analysed in the Results section.

$$S_g = \sqrt{\left(\frac{\partial u_g}{\partial x} - \frac{\partial v_g}{\partial y}\right)^2 + \left(\frac{\partial v_g}{\partial x} + \frac{\partial u_g}{\partial y}\right)^2} \qquad (7)$$

- Energy cascades between different scales. Traditionally the fluxes of energy between different scales have been computed with a spectral method or more recently with the coarse-graining method (Leonard, 1975; Germano, 1992), used for the first time in oceanography by Aluie et al. (2017). The coarse-graining method is
more efficient than the spectral method and allows the generation of 2D maps of energy cascades at a fixed scale L (Schubert et al., 2020; Contreras et al., 2023). It does not require windowing nor the hypothesis of an isotropic field, which avoids a large amount of data loss (and energy loss). Here, the coarse-graining method has been used to study the geostrophic eddy fluxes in the Agulhas Current with the LLC10 simulation. Consequently, the wavelength limit considered is 10 km at 45° S, and smaller scales would not be interpretable from the model and are below the scales observable with SWOT.

The term $\Pi$ representing the scale transfer of kinetic energy can be derived by convoluting the equation of motion:

$$ \Pi = -\rho_0 \left[ (\overline{u^2} - \overline{u^2})\overline{u}_x + (\overline{uv} - \overline{u\overline{v}})(\overline{u}_y + \overline{v}_x) + (\overline{v^2} - \overline{v^2})\overline{v}_y \right], \quad (8) $$

where the overbar represents the convolution of the variable with a top-hat kernel with a diameter size equal to the chosen wavelength to analyse. Refer to Schubert et al. (2020) for specific details about the mathematical derivation of the method. A difficulty in this approach is related to the choice to be made close to continental boundaries. Here, we treated land as water with zero velocity, as in Aluie et al. (2017).

## 3.1 Geographical distribution of EKE and strain rate

Our first objective is to analyse the different eddy diagnostics based on the larger and smaller mesoscale dynamics observable with today's altimetric mapping and with the future SWOT observations, using the LLC10 as our ocean reality in the Agulhas Current region. Figure 3 shows the geographical distribution of the average EKE and the strain standard deviation (SD) over the 1-year period of the LLC10 simulation for three separate cases: the top row is calculated from the total model (LLC10) geostrophic currents, computed from the SSH after correction for barotropic tide and DAC, and represents the total field to be observed by SWOT; the middle row represents the results computed with the pseudo-DUACS data, a proxy of the capability of products derived from current-generation altimeters to measure the larger mesoscale structures; and the bottom data represent their difference: the new scales to be observed with SWOT. We use it to tease out the circulation variability currently not included in the operational products.

In terms of mean surface geostrophic EKE over this 1-year period, the maximum for the total corrected LLC10 current occurs in the Agulhas Retroflection (0.256 m² s⁻²), which is close to the reported observations of the geostrophic EKE computed from DUACS DT2021 data on CMEMS over a longer period (from 1 January 1993 to 3 June 2020, not

shown). The larger scales in the pseudo-DUACS (Fig. 3c) are responsible for most of the energy in the Agulhas Retroflection (0.171 m² s⁻²) compared to the smaller-scale residuals (0.085 m² s⁻²) (Fig. 3e). The larger-scale energy follows the meanders of the main current in the Agulhas Extension. The LLC10 simulation represents physics down to 20 km (Fig. 2b), whereas the pseudo-DUACS is limited by the interpolation scheme. The residual and smaller scales, representing the missing dynamics, have a weaker signature in EKE except along the Agulhas Shelf, especially near Gqeberha (formerly Port Elizabeth) (around 26° E) where the LLC10 current dynamics present greater smaller-scale turbulence. In the CMEMS observations in Fig. A1a, the larger mesoscale develops further south, whereas the model develops weaker large-scale energy offshore, all along the Agulhas Current, with smaller scales dominating on the shelf. This is characteristic of the DUACS interpolation that forces the reconstructed dynamics close to zero near the coasts and thus underestimates the boundary currents and their variability. The strain root mean square variations show a different behaviour. The mean annual strain structure is dominated by the larger scales (not shown). However, Fig. 3 shows that the small-scale, high-frequency dynamics (f) are predominant in terms of variance over the full Agulhas Current system and that as an average over the year, the larger scales (d) are quite stable. The majority of the small-scale strain variability comes from the coastal Agulhas Current and the Agulhas Retroflection and loses its intensity along the Agulhas Extension.

## 3.2 EKE and strain rate variability in four specific regions

As illustrated in Fig. 2, the SSH variance in the Agulhas region is essentially dominated by the rich variability occurring at synoptic to intra-seasonal timescales. In order to analyse the spatial coherency of such variability, we compute the Hovmöller diagrams of the LLC10 geostrophic EKE (Fig. 4), focusing on the region delimited by the time-average current position (grey contours in Fig. 3) along the Agulhas coast and extending our focus to the Agulhas Retroflection (orange line in the geographic plot of Fig. 4). Four zones can be distinguished in the area: zone 1 is the northern part of the coastal region, from where the Agulhas Current starts flowing southwards in the study area. Zone 2 is the coastal extension south of Gqeberha. Zone 3 is the retroflection of the current, and zone 4 is the Cape Cauldron region to the west, surveyed by the QUICCHE CalVal campaign. The first three zones are located along the mean current and are marked on the Hovmöller diagram, from which two sets of distinctive physical processes are distinguishable.

The first is associated with two strong, large mesoscale events called Natal Pulses, represented by the dashed blue lines in the plots. Natal Pulses are solitary meanders that develop near the Natal Bight (around 31° E); they can reach 50 to 200 km in amplitude and have largely been studied with
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Average EKE (left) and strain rate standard deviation (right) over the full LLC10 simulation period (September 2011 to November 2012) (a, b), the pseudo-DUACS product (c, d), and the residuals between the simulation and the pseudo-DUACS product (e, f). The grey contours are the yearly mean current, and the purple contours represent the 3000 m bathymetry. Note that the residual small scales are computed as the difference between the LLC10 full model and the pseudo-DUACS dominated by the larger scales (> 150 km).</figcaption></figure>

satellite altimetry (Schouten et al., 2002; Lutjeharms, 2006). They originate south of the Mozambique Channel and travel downstream, westwards towards the tip of the retroreflection. Up until the moment they reach the Agulhas Bank (around 24° E), continued processes of dissipation and merging occur (Krug and Penven, 2011; Krug and Tournadre, 2012; Krug et al., 2014), and only a fraction of the Natal Pulses reaches this point. The Agulhas Current becomes increasingly unstable downstream, and in addition to the Natal Pulses, shear-edge eddies develop (Lutjeharms and Gordon, 1987), with diameters between 50 and 100 km, which are difficult to observe with altimetry maps. During 2012 in the LLC10 model simulation, two Natal Pulses were reproduced in the coastal region of the Agulhas Current, propagating at an average speed of 6 km d⁻¹. At the northern limit of the domain, they have a size of around 150 km, consistent with the observations and literature. Around mid-March, when the first Natal Pulse reaches Gqeberha, its structure is stretched along the mean current direction and can reach 300 km. The DUACS-reconstructed data, representing the potential of current altimetry, is able to resolve these slow, large-scale processes (in a similar Hovmöller diagram, not shown) in terms of diameter and phase, but their amplitude is underestimated by about 40%. This is a well-documented characteristic of the grid-ded products as a consequence of the mapping algorithms (Ballarotta et al., 2019). Figure 4 highlights that the first Natal Pulse passed box 1 in December 2011, propagated at 6.3 km d⁻¹, and reached box 2 in early March 2012 and then the retroreflection box 3 in early May 2012. An EKE anima-

tion (not shown) shows that after passing the unstable zone of Gqeberha, part of the Natal Pulse is split into one Agulhas ring and propagates in the region of Cape Cauldron. The rest merges in the Agulhas Current and continues down through the Agulhas Extension.

The Hovmöller diagram highlights a second type of much smaller and faster dynamics that have time-varying amplitudes along the Agulhas Shelf from zones 1 to 2 and are almost always present in the zone of the Agulhas Bank, except in the periods when the large Natal Pulses pass. Similar dynamics are modelled in the Gulf Stream (Gula et al., 2016). We find that these submesoscale frontal eddies propagate at a speed of around 20 km d⁻¹ and have a diameter of a tens of kilometres, too small to be observed with gridded altimetry maps. Tedesco et al. (2019) and Krug et al. (2017) demonstrate that southwards from Gqeberha, barotropic instabilities are the main cause of smaller mesoscale eddy generation in a period when no Natal Pulses are occurring. Our LLC10 fields are in accordance, and in the proximity of Gqeberha the Natal Pulse loses most of its energy, and the flow becomes highly unstable, generating faster smaller-scale eddies.

We can quantify how well the different datasets represent these dominant large- and small-scale processes by focusing on the temporal evolution of the EKE for each of the four boxes defined in Fig. 4. Figure 5 shows the EKE on the full period of the LLC4320 simulation from the LLC10 data as well as from the pseudo-DUACS data and their residual. The mean EKE and strain value for each box are shown in Table 1. In the highly energetic Agulhas region, the larger
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Hovmöller diagram of the Agulhas Current coastal region EKE, computed from the *corrected SSH* of the LLC10 simulation. Zones 1, 2, and 3 define specific boxes, shown in the geographical section of the plot. The y axis represents the time, and the x axis represents the distance from the first northern point, following the mean current shown in orange in the geographical plot of Fig. 4. The blue lines show the slower and larger-scale dynamics, so-called Natal Pulses generated in the Mozambique Channel, whereas the red lines represent the faster and smaller-scale ones generated in the Agulhas Bank after the pass of the Natal Pulse.</figcaption></figure>

mesoscale represented by the pseudo-DUACS product (or-ange line in Fig. 5) structures the flow in the coastal box (a), but the large and small scales have similar mean values. In the retroflection (c) and Cape Cauldron (d) boxes the large scales dominate the EKE, with a mean value close to dou-ble the small scales in both cases. Yet, in comparison with the total LLC10 data (blue line), the pseudo-DUACS under-estimates the magnitude of the larger eddies in strong events by 30%–40%. The biggest features here are the two Natal Pulses propagating downstream from coastal box 1, along the shelf break (box 2) to the retroflection area (box 3), as shown in Fig. 4.

In these three boxes, the amplitude of the smaller-scale dy-namics generally increases in periods when the larger-scale eddies are active, suggesting that energy transfer processes are activated during this period. Certain smaller-scale energetic processes are smoothed out by the pseudo-DUACS re-construction, such as in November 2011 in box 3 and the large peak in late July 2012 in the Cape Cauldron region. Box 2 is different, being a shelf region dominated by small-scale, rapidly propagating dynamics that are stronger on av-erage than the larger scales in the same region, except during the periods when the two Natal Pulses propagate downstream from box 1. Figure 4 shows that these small-scale, rapid dy-namics have speeds of 20 km d⁻¹ on average, similar to the submesoscale frontal eddies generated by barotropic insta-bilities on the Agulhas Plateau as described by Ruijter et al. (1999), Krug et al. (2017), and Tedesco et al. (2019). They are weakly present in the early period of the simulation, but their amplitude increases during the 3 months in autumn fol-lowing the passage of the first Natal Pulse. During the fol-

lowing winter-spring months (July–October), the dynamics return to a less energetic state in box 2, as this barotropic instability occurs when the Agulhas Current is not in a me-andering state (Tedesco et al., 2019). This critical small-scale adjustment process is undetected in today’s altimetric maps.

To complete the analysis we focus on the flow’s deforma-tion with the potential of generating secondary vertical over-turning circulation at different scales (Zhang et al., 2019). Figure 6 shows the geostrophic strain rate variability over two of the boxes (box 1: a–b, box 2: c–d). The left column shows the spatial average of strain over the box, whereas the right column shows the spatial SD. This figure, together with Table 1 which quantifies the average strain rate values over the boxes and over 1 year, confirms that the larger-scale pseudo-DUACS geostrophic strain rate contributes around 60% of the total average LLC10 strain rate amplitude in box 1 (and in boxes 3 and 4, not shown). The residual rapidly evolving small-scale geostrophic strain contributes 40% of the mean strain rate but most of its variability as shown in Figs. 3d, 3f, 5b, 5d.

The Natal Pulses, the strongest events in terms of EKE, only have a small signature on the strain, especially on the pseudo-DUACS strain. They propagate through the region (a, b, c, d) with very high SSH and EKE values but with a less strong impact on the deformation of the flow: the very large EKE peaks in December 2011 and August 2012 in box 1 (Fig. 5a) have moderate DUACS strain peaks (Fig. 6a, b), but other periods with weaker EKE also have moderate pseudo-DUACS strain peaks in box 1, of the same order of magni-tude. Indeed, the passage of the Natal Pulses tends to damp out the smaller-scale strain (evident in the SD, Fig. 6b), but
<table><caption>**Table 1.** EKE and strain mean on the four zones identified in Fig. 4.</caption><thead><tr><th rowspan="2"></th><th colspan="2">LLC10</th><th colspan="2">Pseudo-DUACS</th><th colspan="2">Residuals</th></tr><tr><th>eke<br>[m<sup>2</sup>s<sup>-2</sup>]</th><th>strain<br>×10<sup>-5</sup> [1 s<sup>-1</sup>]</th><th>eke<br>[m<sup>2</sup>s<sup>-2</sup>]</th><th>strain<br>×10<sup>-5</sup> [1 s<sup>-1</sup>]</th><th>eke<br>[m<sup>2</sup>s<sup>-2</sup>]</th><th>strain<br>×10<sup>-5</sup> [1 s<sup>-1</sup>]</th></tr></thead><tbody><tr><td>Coast</td><td>0.111</td><td>2.3</td><td>0.06</td><td>1.4</td><td>0.051</td><td>0.9</td></tr><tr><td>Coastal extension</td><td>0.255</td><td>2.9</td><td>0.104</td><td>1.4</td><td>0.152</td><td>1.5</td></tr><tr><td>Retroflection</td><td>0.256</td><td>2.2</td><td>0.171</td><td>1.3</td><td>0.085</td><td>0.9</td></tr><tr><td>Cape Cauldron</td><td>0.072</td><td>1.3</td><td>0.047</td><td>0.7</td><td>0.025</td><td>0.6</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** EKE (m² s⁻²) in the four regions defined in Fig. 4: the upstream coastal zone (a), coastal extension (b), Agulhas Retroflection (c), and Cape Cauldron zone (d). The four superposed lines correspond to the full LLC10 simulation (blue), the pseudo-DUACS data (orange), the residual high-frequency smaller scale (green), and the large-/small-scale cross-terms (red).</figcaption></figure>

the total strain in box 1 remains fairly constant throughout the year, with a mix of larger- and smaller-scale contributions. In contrast, the effect of the Natal Pulse is evident on the cross-terms (red line in Fig. 6a), whose effect increases during the pass of the Natal Pulses and causes the dampening of the smaller scales.

An exception is box 2 (Fig. 6b, c), where the small-scale strain is the most active all year, except when the Natal Pulse passes. Small-scale strain remains present even during the low-EKE period (August to October 2012), with implications for vertical overturning and mixing at small scales in this region. In contrast with Sasaki et al. (2014), there is no evidence during this 1-year simulation of small mesoscale EKE in winter and spring feeding energy to larger mesoscale EKE in late spring and summer. Indeed, the seasonal averages in this 1-year simulation are dominated by more individual eddy events. The seasonal mixed-layer depth remains shallow all year round, in contrast to the strong seasonal mixed-layer depth changes in the North Pacific (see climatology in Fig. 4 of Johnson and Lyman, 2022). So seasonal dynamics in the Agulhas Retroflection are weak compared to

mesoscale events and do not strictly follow the paradigm of mixed-layer instabilities, in contrast to Sasaki et al. (2014).

### 3.3 Geographical distribution of the energy cascade

One of the main aims of SWOT will be to understand and monitor a wider range of mesoscale structures, whether they act to reinforce or compensate for the larger-scale variability, and to observe the interaction and transfers of energy between the different spatial scales. A direct (or downwards) cascade is defined as the transfer of kinetic energy towards smaller scales. An inverse (or upwards) cascade is the transport of kinetic energy to progressively larger scales. We analyse the geostrophic energy fluxes in our region using the coarse-graining method, as displayed in Fig. 7. Note that although the SSH observed by SWOT includes both geostrophic and ageostrophic components that dominate sub-mesoscale features, only the geostrophic ones are accessible and can be used for computing energy transfers.

In this analysis, our subsampled model at 1/10° around 45° in latitude has a Nyquist sampling of around 15 km, in line with the SWOT observability (not far off from 15 km) from Sect. 4. Since the objective is to perform a compar-
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Time series of geostrophic strain rate (1 s⁻¹) in the two coastal boxes for the full LLC10 simulation (blue), the pseudo-DUACS data (orange), the residual high-frequency smaller scale (green), and the large-/small-scale cross-terms (red). The panels (a, c) represent the time series of the spatial average over box 1 (a) and box 2 (c). The panels (b, d) are the time series of the spatial SD over box 1 (b) and box 2 (d).</figcaption></figure>

ison with the pseudo-DUACS products, Fig. 7 shows the geostrophic energy flux at 60 km (top) and 150 km (bottom) for the LLC10 data (a, c) and for the pseudo-DUACS product (b, d). Shorter wavelengths would not be representative of the pseudo-DUACS data. In line with Schubert et al. (2020), the inverse cascade’s (blue) strength and spatial distribution increase at the larger 150 km scale. The full “SWOT-like” LLC10 presents a stronger inverse cascade, and the change of sign between upscale and downscale fluxes occurs at shorter wavelengths. This has already been observed by Renault et al. (2019) in the Agulhas Current and the Gulf Stream and indicates that current altimetric maps are unable to accurately estimate the spatial scales at which the energy cascade changes from an inverse to a direct cascade. The “patched” shape of the flux may be due to the limited duration of our model: only 1 year of data may not be long enough in this active region to provide stabilized mean fluxes, especially given the very energetic and diverse events shown in Fig. 4. Both Schubert et al. (2020) and Contreras et al. (2023) use a multi-year simulation for their flux computation which shows an inverse cascade dominated by the balanced flow, becoming stronger as larger wavelengths are analysed. At 150 km, the zone of the direct cascade is larger than at 60 km, and it is concentrated east of Gqeberha (around 26° E). It corresponds to the zone and wavelengths where Krug and Penven (2011) found the generation of smaller eddies after the passage of the Natal Pulses, in relation with the downstream decrease in the number of Natal Pulses south of Port Edward (31° S, 30° E). In Fig. 5a the Natal Pulse propagating downstream from the Natal Bight breaks into many smaller and faster dynamics, shown in Fig. 5b, creating a direct geostrophic en-

ergy cascade. As seen in the time series, this very large and energetic event is the only one that is fairly well represented by the pseudo-DUACS product, although the inverse cascade it is largely underestimated.

This coarse-graining diagnostic will not be reproduced at the maximum resolution along the SWOT tracks due to the narrow width of the two 50 km swaths. The diagnostics could be revisited once SWOT is included in future multi-mission CMEMS maps.

## 4 Eddy diagnostics observed by SWOT

The EKE and strain rate diagnostics have been computed on the swaths generated with the SWOT simulator by interpolating LLC10 onto the 2 km regular SWOT grid. The aim is to understand how these diagnostics will be observed with SWOT data, how much we will lose in terms of swath width, and how the residual instrumental error will impact the observability of the SSH and the eddy diagnostics derived from it.

The strain rate and the EKE are both invariant with respect to the system coordinates (see Appendix B). Therefore, we can compare these diagnostics from the 2D LLC10 maps computed with the geostrophic velocities in the zonal–meridional framework with similar diagnostics derived from SSH on the swaths computed in the along-track–across-track framework.
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Geographical distribution of the yearly average (September 2011 to November 2012) direct (red) and inverse (blue) energy cascade in the Agulhas region computed with the coarse-graining method for 60 km (top) and 150 km (bottom) scales. Panels (a, c) are computed with the LLC10 data and (b, d) with the pseudo-DUACS data. The grey lines are the mean SSH contours, and the purple and yellow lines are the 3000 and 1000 m bathymetry respectively.</figcaption></figure>

## 4.1 Impact of the SWOT instrumental noise

The objective of this section is to analyse the impact of KarIn noise on the SSH measurements and therefore on the derived eddy diagnostics. We also infer the capability of the U-Net noise mitigation technique, trained on eNATL60, to treat a different model (LLC10) in a different and very energetic zone. After de-noising, we expect to be able to see the smallest possible ocean structures. Gómez-Navarro et al. (2018), Chelton et al. (2022), and Tréboutte et al. (2023) showed that KarIn's random error adds noise to the signal features that could lead to a misinterpretation of the observed variable or completely mask the signal, making it impossible to interpret the data. The impact of the small-scale random noise increases when making first- or second-order spatial derivatives, such as the EKE or the strain. Figure 8 demonstrates this point, showing the SSH (left), EKE (centre), and strain (right) in the ideal LLC10 model scenario (top), when random noise is added (centre) and after U-Net is applied to remove the noise (bottom). All plots refer to pass number 5 (geographically shown in Fig. 1). The noisy SSH (b) presents the same large-scale features with respect to the simulated true field (a), but smaller-scale features are introduced or modified, leading to a misrepresentation of the SSH variability at small spatial scales. The noisy EKE (e) only shows the strongest feature in the middle of the swath, and the strain (h) is completely covered by the noise, preventing any interpretation of this variable. The U-Net noise mitigation technique, with parameters derived from a different model in the North Atlantic, suitably restores the input signal in the three cases. For the SSH, for wavelengths between 15 and 50 km, the noise is reduced by 1 order of magnitude. However, U-Net also removes part of the signal by a few percent points over this interval. Thus, the ocean SSH dynamics observed after noise reduction will have the correct positioning and

phase, but their amplitude is slightly underestimated. The interval in which the signal is slightly underestimated is larger for higher-order derivatives, being between 10 and 100 km for the strain.

The scores used to assess the impact of U-Net are defined in Sect. 2.3.2. The RMSE is shown in Fig. 9. A synthetic overview of the overall noise mitigation efficiency in the Agulhas Current zone with the MITgcm model is given in Table 2, comparing the scores and observability of this study to the original U-Net simulation based on the eNATL60 North Atlantic study (Tréboutte et al., 2023) and the use of the filter implemented by Gómez-Navarro et al. (2020). The results are promising, with a homogeneous RMSE for the four daily-sampled tracks in the Agulhas region, ranging between 0.20 and 0.45 cm, which is in line with the findings from Tréboutte et al. (2023) for the North Atlantic. The RMSE across-track shape reflects the simulated random error that increases towards the inner and outer swath borders (Esteban-Fernandez, 2013).

## 4.2 SWOT observability

Ocean observability is defined as the wavelength for which the NSR equals 1, meaning that the noise and the signal have equivalent energy. Following Tréboutte et al. (2023), for noisy fields where the error is always more energetic than the signal (NSR always higher than 1), we consider the observability to be larger than 1000 km. Previous findings by Wang et al. (2019) show that SWOT observability is expected to be degraded and to resolve wavelengths up to 35–45 km at high latitudes due to SWH-induced instrumental noise. This is especially true in the Antarctic Circumpolar Current (ACC) where the wave-induced random error has a large seasonality. In their simulation, Wang et al. (2019) used the realistic WAVEWATCH III (WW3) model as input data for waves, and they showed that spatial filtering of the data to smooth
<figure><img src="image_10.png" /></figure>

**Figure 8.** SSH (a, b, c), EKE (d, e, f), and strain (g, h, i) in the ideal case of the non-noisy field (a, d, g), the noisy field (b, e, h), and after noise mitigation with the U-Net method (c, f, i). All plots refer to pass number 5 in Fig. 1. The swath is represented horizontally with its distances along-track (a-track) and across-track (x-track) in kilometres. The SSH is a snapshot of cycle 112 on 1 January 2012. EKE and strain are averaged over 3 months, simulating the CalVal scenario (January–March 2012). The left side of the swath corresponds to the southern section of the track, and the right side reaches the South African coast.

<figure><img src="image_11.png" /></figure>

**Figure 9.** SSH root mean square error for each of the SWOT 1d CalVal tracks in our zone, averaged over 3 months (January–March 2012).

<figure><img src="image_12.png" /></figure>

**Figure 10.** Estimating the NSR as the ratio of the power spectral density of the SSH error compared to the full LLC10 SSH fields for two cases: in green, is the spectrum of the noisy SSH, and in orange is the spectrum of the U-Net-treated SSH. The vertical dashed lines indicate the wavelength at which the signal reached the noise (NSR = 1). The PSD of the error and the noise are calculated for along-track PSD over the 1600 km segment and averaged for each across-track position, and for each track and cycle.

the noisy smaller wavelengths destroys much of the information contained in the small-scale SSH field. Figure 10 shows the LLC10 SSH noise-to-signal ratio, before and after applying the U-Net method over our four SWOT tracks. Even though we are using a noise simulated for a 2 m wave field, as in Wang et al. (2019), the wavelength SWOT should resolve in this region before de-noising is 40 km, whereas when U-Net is applied to the noisy data, an observable wavelength of 17.5 km can be restored. The resolved observable scale is also calculated on the EKE and on the strain rate. For noisy data, the observability decreases by a factor of 3 for first-

order derivatives such as EKE, resolving scales of more than 120 km with noisy data. For second-order derivatives such as geostrophic strain, the noise dominates the signal for all wavelengths, since the NSR spectrum is always > 1. After U-Net application, the observability is restored to wavelengths of the order of 15 km in both cases.

Considering the results for these de-noising scores and the comparison with other noise reduction techniques used in the literature, the U-Net method is very promising for the Agulhas Current when analysing a realistic simulation. In terms
of SSH mean RMSE, without any noise treatment, the values in the Agulhas region and in the North Atlantic are comparable, 1.17 cm and 1.27 respectively. Using the U-Net method in the Agulhas region, the final RMSE is 0.24 cm, which is the same found by Gómez-Navarro et al. (2020) with specific tuning of their algorithms for the North Atlantic. The U-Net performance in the North Atlantic has an even better error reduction of 0.19 cm. In terms of the variance of the SSH residuals, the results are similar, with the best result found with U-Net in the North Atlantic and U-Net in the Agulhas region and the Gomez filter in the North Atlantic having comparable values. This result was expected because the Gomez filter was parameterized specifically on their study zone, and the U-Net training was performed in the North Atlantic with a different model. What is remarkable about the U-Net performance is that even with a different model and in a different zone, the mean RMSE and variance of SSH residuals are comparable with the original parametrization of the other methods. The SSH observable wavelengths give the most surprising result in this sense. The two zones have comparable observability before noise reduction (40.5 km in the Agulhas region and 42 km in the North Atlantic). U-Net in the Agulhas region manages to retrieve wavelengths down to 17 km. This result in the Agulhas region could be improved in the future by performing new training in this specific zone. This is encouraging for the early SWOT processing, since the U-Net parameters, trained in the North Atlantic, could be applied directly to the early SWOT data and then be improved by retraining U-Net once enough SWOT data have been retrieved.

Finally, this study was conducted in preparation for the real SWOT data that will represent a new challenge, since in the early months of the mission's data analysis, no training dataset will be available for the U-Net algorithm. Note that in all the noise mitigation techniques presented in Table 2, the simulated random noise is estimated in the same way, based on the SWOT project's best spectral estimates of the random noise for 2 m SWH (Gaultier et al., 2016). In reality, the SWOT random noise will vary in amplitude and could be impacted by small-scale anomalies that will need extra editing (e.g. presence of isolated ships, icebergs, platforms, extreme waves). However, the U-Net performance in a different region and using a different model from the original training is very encouraging, not just for the restored SSH, but also for the key high-order eddy diagnostics such as EKE and strain.

## 5 Conclusions and discussion

Our study addresses the Agulhas Current mesoscale dynamics and its observation with current altimetry maps and with the new 2D interferometric SAR swath sampling brought by Surface Water and Ocean Topography (SWOT). We investigated how eddy diagnostics can be observed with a synthetic pseudo-DUACS dataset reconstructed from the LLC10 simulation

data and how SWOT will change the current paradigm by observing new space and time frequencies. Finally, we studied the effect of instrumental random error on SWOT observations and quantified the potential of a new de-noising algorithm based on a neural network approach.

To reproduce altimetry-like data from our model, we corrected the LLC10 SSH for the barotropic tide and the DAC as a standard altimetric data processing step and then calculated the geostrophic currents from the corrected SSH. Some high-frequency signals remain in the corrected SSH fields due to residual barotropic tides, internal tides, and internal gravity waves (see Fig. 2). These are minimized in our analysis since we are dealing with daily-averaged model data for the pseudo-DUACS and therefore with the daily residuals of small-scale fields. In reality, SWOT will fly across this region at $7 \text{ km s}^{-1}$ (passing over the $15^\circ$ in latitude in 4 min), and so during the CalVal phase, SWOT will measure daily 2D snapshots including these residual high-frequency signals. To verify the impact of this on our eddy diagnostics, we compared the daily-averaged EKE derived from the full LLC10 fields versus a daily snapshot from the hourly LLC10 model data. In the high-energy Agulhas region, the differences were minimal (see Appendix C). Even in the area of higher internal tides radiating westwards from the Benguela Current, the impact on the spatial derivatives of SSH such as EKE was small. The Agulhas Current is one of the most energetic regions of the global oceans, and other less energetic regions may need to use specific internal tide corrections (Zaron and Ray, 2017) or more sophisticated techniques to separate the high-frequency internal gravity wave field from the rotation-dominated eddy fields (Le Guillou et al., 2021).

The separation between large scales ($> 150 \text{ km}$) from the DUACS-like reconstruction and the residual small scales ($< 150 \text{ km}$) with the total LLC10 simulation data has allowed us to quantify the EKE that is missing in today's observations during strong events like the Natal Pulses. In the 2011–2012 period, LLC10 shows two Natal Pulses flowing southwards from the Natal Bight to Gqeberha. These large structures are well known and have previously been documented with altimetry and satellite sea surface temperature (SST) data (Krug and Penven, 2011; Krug et al., 2017). Our analysis shows that their location and phase are correctly estimated by the DUACS-like reconstructed observations, but we miss 30%–40% of the magnitude due to the interpolation used in the mapping process. SWOT should also observe snapshots of the smaller and faster dynamics southwards of Gqeberha, which are mainly created by barotropic instabilities (Krug et al., 2017; Tedesco et al., 2019) and are totally missing in current altimetry products due to their small diameters. Our EKE analysis also highlighted that the cross-terms between the small- and larger-scale dynamics are also important contributors to the total EKE field that are not observable today.

The Agulhas region is also characterized by strong mesoscale to submesoscale strain, estimated at $6–8 \times$
<table><caption>**Table 2.** Summary of the main noise reduction and observability scores for the North Atlantic and the Agulhas regions. North Atlantic de-noising compares the U-Net method and the filter from Gómez-Navarro et al. (2020) to the case with no filtering of the random noise. The Agulhas Current has only been treated with the U-Net method.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Agulhas region</th><th colspan="3">North Atlantic</th></tr><tr><th>No filter</th><th>U-Net</th><th>No filter</th><th>U-Net</th><th>Gomez</th></tr></thead><tbody><tr><td>SSH RMSE [cm]</td><td>1.17</td><td>0.24</td><td>1.27</td><td>0.19</td><td>0.24</td></tr><tr><td>Variance in SSH res [cm<sup>2</sup>]</td><td>1.5</td><td>0.06</td><td>1.63</td><td>0.04</td><td>0.07</td></tr><tr><td>SSH wavelength [km]</td><td>40.5</td><td>17.5</td><td>42</td><td>10</td><td>27</td></tr></tbody></table>

10⁻⁶ s⁻¹ by Zhang et al. (2019) based on gridded altime-try estimates (their figure may be saturated at its highest strain scales). Our modelled strain rate is a factor of 2–4 times larger, with mean values of 1.5 to 3 × 10⁻⁵ s⁻¹ for all scales. We showed that although the modelled average strain rate magnitude is dominated by the large mesoscales in most regions, except in box 2 where the smaller-scale barotropic instabilities emerge, most of the strain variability occurs at smaller scales. Currently, available observations completely miss this information as the pseudo-DUACS data reveal a fairly constant and low strain variability in all regions and times of the year, even during the strongest changes in the EKE carried by the Natal Pulses.

Correctly assessing these dynamics is crucial because the Agulhas Current nearshore region has a rich upwelling system (Blanke et al., 2009; Goschen et al., 2015; Jacobs et al., 2022). Small-scale and rapid changes in the SSH, EKE, and strain as described here are key for the stability and generation of local upwelling pulses and the associated rich biomass (Largier et al., 1992; Zhang et al., 2019).

The resolution of these smaller dynamics will be essential to correctly observe the ocean's energy cascade and in particular the inverse cascade which today is underestimated in magnitude and shifted to larger wavelengths (Renault et al., 2019). In future SWOT measurements, the smaller-scale ageostrophic currents will not be accessible, and even the geostrophic eddy cascade will have a spatial resolution of around 20 km at best. Thus, the interaction between balanced and unbalanced motions, which at small scales Contreras et al. (2023) showed to be a key component for the knowledge of the forward cascade in the Gulf Stream, will not be accessible. However, geostrophic balanced motions dominate for the SWOT-observable scales in both the Gulf Stream (Contreras et al., 2023) and the Agulhas region (Schubert et al., 2020). Thus, within the scales of interest for SWOT, we will not capture most of the inverse cascade and the geostrophic contributions of the direct cascade at scales larger than 15 km such as the barotropic instabilities observed near Gqeberha. The coarse-graining method will not be directly implemented on the SWOT swaths as the two 50 km wide swaths are not wide enough to correctly implement the 2D methodology used here. However, in 2024, a new gridded high-resolution 2D DUACS product will be

implemented that will also include SWOT data. We expect this new product to be able to represent a wider spectrum of spatial scales, making it possible to interpret the geostrophic energy fluxes between scales.

The results presented in this study are subject to how well the LLC10 simulation and the SWOT simulator represent reality. A few caveats are known for the model simulation, as discussed in Sect. 2. The simulator also includes estimates of SWOT's systematic errors (satellite roll, phase errors, baseline dilation, timing errors, orbital errors), and the SWOT project is developing cross-calibration techniques to estimate and remove these additional errors (Dibarboure et al., 2022). Due to their nature, these platform errors are expected to be significant only for wavelengths larger than 1000 km. However, high-frequency residuals could still be present at spatial scales shorter than 1000 km and therefore would translate into the SSH measurements and the higher-order derivatives. These residual errors, if present, would be subject to higher-level mission calibration against the available altime-try constellation (Dibarboure et al., 2022). Even after this cross-calibration step is applied, some residual errors may remain, and their impact on the derivation of eddy diagnostics needs to be addressed in future work. The main differences with the real data are related to the generation and treatment of KarIn's random noise. Here a Gaussian distribution and a constant SWH of 2 m have been used to estimate the noise. The real SWH, however, can reach up to 7 m in the Southern Ocean, which would potentially produce a higher random noise and lower overall observability, and the noise may not have the same statistical distribution. However, early estimates by the SWOT project suggest that the SWOT random noise may be smaller than in our simulations, unexpected and good news for future eddy observations.

Finally, the U-Net noise mitigation is very promising. The technique slightly reduces the signal (by a few percent) but does not over-smooth the gradients and retains the main structures and anomalies up to the coast. This is a strong benefit for eddy diagnostics and provides better noise mitigation with respect to the other filters analysed by Tréboutte et al. (2023). Our study shows that the U-Net trained for the eNAtl60 model (North Atlantic zone) is still providing good results in the Agulhas region during a different period and with data from another simulation that accounts for different
ocean SSH dynamics. One of the reasons for the U-Net effi-
ciency is that the U-Net has already been trained in the North
Atlantic for different wave heights and seasons trying to be
representative of ocean dynamics in other regions. Moreover,
the simulated SWOT random errors are based on global spec-
tral estimates and are not regionally varying. This may not be
the case for the early SWOT data in 2023, potentially with a
geographically varying random error. In applying the U-Net
technique, we also take care that the input SSH is not too dif-
ferent from the one used for the training to maintain a good
performance. For this, the training SSH and the SSH are both
normalized with respect to their variance and mean before the
U-Net technique is applied, and then the inverse normaliza-
tion is applied to recover the correct input SSH. One of the
disadvantages of the U-Net method is that it needs a solid
training dataset that will not be readily available for the early
SWOT data. However, the promising results of this Agulhas
study, based on a different model and region, are very en-
couraging. The parameters derived from the U-Net technique
trained on simulated data and noise may be used for the first
random noise correction for SWOT, and it will be interesting
to see the outcome and correct for potential differences when
applying it to the real data.

Appendix A: Validation of the LLC4320 simulation

The LLC4320 has been validated by numerous studies: Rocha et al. (2016, 2015) compared the model kinetic energy spectra in the Drake Passage against long series of in situ acoustic Doppler current profiler (ADCP) data, finding good agreement for both the rotational and the divergent components of the 1D Helmholtz decomposition. Drushka et al. (2018) analysed the internal tide component of the 1/48° MITgcm data and found good agreement with in situ data from a glider. Wang et al. (2018) demonstrated that the simulation reproduces the location and amplitude of kinetic energy peaks well, comparing the simulation to 25 repeat ADCP surveys in the northwestern Pacific Ocean. In the North Pacific Ocean, Savage et al. (2017) tested the performance of two global ocean simulations (the HYbrid Coordinate Ocean Model at 1/24° resolution and the LLC4320) to reproduce the diurnal, semidiurnal, and supertidal variance in the SSH against nine McLane profilers. Both models agree on the diurnal and semidiurnal tidal ranges, but discrepancies were found at supertidal frequencies. The LLC4320 variance in the SSH was found to be closer to the profiler variance because its higher spatial resolution allows it to better reproduce the energy transfer out of the inertial and semidiurnal bands.

We estimate in Fig. A1 the geostrophic EKE temporal
mean for the LLC10, the observed DUACS all-sat product
from AVISO, and the pseudo-DUACS product. In order to
benchmark the realism of the model physics, we want to
highlight the large-scale circulation variability rather than

<figure><img src="image_13.png" /></figure>

**Figure A1.** Geographical maps of EKE averaged over the period of September 2011 to November 2012. (a) Observed altimetry DUACS all-sat maps. (b) Pseudo-DUACS product. (c) LLC10. The grey contour lines represent the mean current. The purple lines represent the bathymetry at 3000 m.

the mean state by removing the annual cycle from the full
data for the three cases (not shown). For the derivation of the
EKE, see Sect. 3. The time average over the full simulation
period (September 2011 to November 2012) shows that due
to their smoothness, the real DUACS product and the pseudo-
DUACS underestimate the EKE levels of energy with respect
to the LLC10 simulation, while retaining the spatial features
of the current system, in agreement with other studies (Chel-
ton et al., 2019), in particular along the coastal Agulhas Cur-
rent. The pseudo-DUACS underestimates the observed DU-
ACS/CMEMS reconstruction in the retroflection region (37–
42° S, 15–20° E) (see Fig. A1a and b), the mean EKE in the
retroflection being 0.157 m² s⁻² for the pseudo-DUACS and
0.178 m² s⁻² for the real DUACS, whereas the LLC10 has a
mean EKE of 0.256 m² s⁻². The mean axis of the coastal Ag-
ulhas Current, its retroflection around 20° E, and the mean-
dering extension are more clearly delineated with the LLC10
and much more diffuse in the real DUACS maps. The LLC10
field has more energy at the retroflection south of 40° S and
along the coastal Agulhas Current. The mean DUACS EKE
in the extension is slightly shifted to the north with respect
to the modelled LLC10 and pseudo-DUACS mean current.
The coastal Agulhas Current region has larger differences: at
<figure><img src="image_14.png" /><figcaption>**Figure A2.** Snapshot of normalized surface relative vorticity on 1 January 2012. The grey contour lines represent the mean current. The purple and yellow lines represent the bathymetry respectively at 3000 and 1500 m. (a) Pseudo-DUACS-reconstructed data from LLC10. (b) LLC10.</figcaption></figure>

33° S, the mean EKEs in the DUACS maps and in the LLC10 are 0.039 and 0.111 m² s⁻² respectively.

Figure A2 represents a snapshot of the geostrophic surface relative vorticity (see Sect. 3) for the LLC10 and for the pseudo-DUACS product. This qualitative comparison demonstrates that the LLC10 model has energetic velocity gradients and gives a good representation of small-scale vorticity structures, whereas the reconstructed DUACS represents the amplitude and the location of the main, large-scale vorticity structures but, as expected, does not represent the small-scale vorticity in terms of rotational and strained flow. Since the LLC10 provides a good representation of the EKE and the small-scale vorticity and since the objective of this study is to understand how SWOT will observe the variability compared to what is achieved today, the LLC10 is a good fit and has been used in our analysis. Furthermore, it has the advantage of consistent sampling with the additional fields used in this study: the pseudo-DUACS data (see Sect. 2.2) and the SWOT cross-calibration data mentioned in the “Conclusions and discussion”.

## Appendix B: Well-defined variables invariant to coordinate system rotation

In order to calculate the Eulerian diagnostics on a rotated along-track-cross-track field, we need to use well-defined variables. The full mathematical derivation can be found in Wirth (2015). In the following, we have a 2D flow field composed by the velocity vector $(u(x, y, t), v(x, y, t))$. Some scalar quantities, referred to as well defined, are invariant with respect to the coordinate system, which might be moved or rotated without any impact on the quantity's values. As an example, the speed is a well-defined variable as opposed to

the velocity vector, which changes when the coordinate system is rotated. The linear deformation of a fluid volume is characterized by the strain tensor. If $u$ and $v$ are components of the velocity vector in a first reference frame in a 2D flow, the strain tensor is defined as

$$ (\nabla \mathbf{u}')^t = \begin{pmatrix} \partial_x u & \partial_y u \\ \partial_x v & \partial_y v \end{pmatrix}. \qquad (B1) $$

For a vector, the rotation to a second coordinate system is given by the rotation of the first system by an angle $\alpha$, through the left multiplication by the rotation matrix $\mathbf{A}$:

$$ \begin{pmatrix} u' \\ v' \end{pmatrix} = \begin{pmatrix} \cos\alpha & \sin\alpha \\ -\sin\alpha & \cos\alpha \end{pmatrix} \begin{pmatrix} u \\ v \end{pmatrix} = \mathbf{A} \mathbf{u}. \qquad (B2) $$

Its gradient transforms with the right multiplication by the transpose of the matrix $\mathbf{A}$. Thus, the strain tensor transforms as follows:

$$ \begin{pmatrix} \partial_{x'} u' & \partial_{y'} u' \\ \partial_{x'} v' & \partial_{y'} v' \end{pmatrix} = \mathbf{A} \begin{pmatrix} \partial_x u & \partial_y u \\ \partial_x v & \partial_y v \end{pmatrix} \mathbf{A}^t, \qquad (B3) $$

where $\mathbf{A}' = \mathbf{A}^{-1}$. Note that several well-defined quantities can be derived from the strain tensor, and they depend linearly on the strain tensor:

- trace of tensor $d = \partial_x u + \partial_y v$,

- vorticity $\zeta = \partial_x v - \partial_y u$,

- determinant $D = \partial_x u \partial_y v - \partial_y u \partial_x v$,

- square of the components of the strain matrices $H = d^2 + \zeta^2 - 2D$,

- strain rate $s^2 = d^2 + \zeta^2 - 4D = H - 2D = (\partial_x u - \partial_y v)^2 + (\partial_x v + \partial_y u)^2$,

- Okubo-Weiss parameter $OW = s^2 - \zeta^2 = d^2 - 4D$.

The strain rate is hence invariant with respect to the system coordinates and so is the EKE, derived from the speed. This justifies the comparison of the diagnostics on the 2D maps computed with the geostrophic velocities in the zonal-meridional framework with the diagnostics on the swaths computed in the along-track-across-track framework.

## Appendix C: Impact of internal gravity waves in the Agulhas region

The Agulhas region has very strong mesoscale energy and cascades that largely dominate any internal tide signals at these relatively high latitudes. We made a careful analysis of the impact of the high-frequency residuals in the studied zone.

We compared modelled geostrophic EKE and strain rate statistics computed from (a) the full hourly time series
<figure><img src="image_15.png" /><figcaption>**Figure C1.** Percentage of variance explained by the low-frequency dynamics relative to the full LLC10 hourly EKE. The larger the percentage, the more the full dynamics are due to slow (> 24 h) dynamics. Low percentages indicate that high-frequency (< 24 h) processes contribute to the total dynamics, but these do not have the typical structure of internal tides.</figcaption></figure>

and (b) the daily average of the full time series. Each of these datasets represents slightly different high-frequency dynamics. The full hourly dataset represents the corrected LLC10 SSH dynamics, including the hourly evolution of the mesoscale and submesoscale structures, internal gravity waves, and internal tides. The daily averaged data minimize the variations < 24 h, and we use them as a proxy for the geostrophic motion.

Figure C1 quantifies the percentage of variance explained by the low-frequency (LF) dynamics EKE (averaged over 1 d, dynamics > 24 h) relative to the full hourly EKE, as per Eq. (C1):

$$ \left( 1 - \frac{\text{var}(\text{LLC10}_{\text{EKE}} - \text{LF}_{\text{EKE}})}{\text{var}(\text{LLC10}_{\text{EKE}})} \right) \times 100. \quad (\text{C1}) $$

This tells us how much the low- and high-frequency dynamics contribute to the total SD. We see that the percentage of variance explained by the low-frequency (> 24 h) dynamics is large everywhere, meaning that the full dynamics are mainly represented by the lower-frequency dynamics, except on the shallow plateau and off the Benguela coast where higher-frequency dynamics dominate. The white contours represent a mean EKE of 0.04 and 0.16 m² s⁻² (see Fig. 3a). In the region of stronger EKE (> 0.04 m² s⁻²), 7% of the total SD is due to some high-frequency (HF) dynamics that are difficult to clearly separate from the mean dynamics but that do not have the typical structure of internal gravity waves or internal tides.

*Code availability.* The code developed for the SSH-based diagnostics was created using Python 3.9 and is available upon request. The SWOT simulator has been developed by the Centre National d'Études Spatiales (CNES) and is publicly available on the CNES GitHub (https://github.com/CNES/swot_simulator, Ubelmann et al., 2021). The U-Net noise mitigation algorithm has been developed by CLS and is under a proprietary licence that pre-

vents its diffusion. Its use has been granted in the frame of this PhD thesis.

*Data availability.* The LLC10 simulation data are a re-gridded version of the LLC4320 simulation made by Collecte Localisation Satellites (CLS) in Toulouse, France, and made available upon request (last access: 25 January 2022). The global LLC4320 SSH surface fields are publicly available on https://sextant.ifremer.fr/geonetwork/srv/api/records/b2bcb9af-f335-45b6-a2a9-e460e4132879 (Ponte, 2020). The specific pseudo-DUACS product on the Agulhas region was built by CLS with a DUACS processing (as in Taburet et al., 2019) and is available upon request by contacting Elisa Carli (elisa.carli@univ-tlse3.fr).

*Author contributions.* Conceptualization and methodology: EC, RM, OV; software: EC and RC; writing: EC, RM, OV; supervision: RM, OV, LR. All the authors have read and have agreed to the published version of the paper.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

*Acknowledgements.* Elisa Carli was supported by the Centre National d'Études Spatiales (CNES) and Collecte Localisation Satellites (CLS) with a co-funded research grant to conduct her PhD. This work was financed through the CNES SWOT TOSCA project.

*Financial support.* This research has been supported by the Centre National d'Études Spatiales (CNES/CLS PhD co-funding grant).
*Review statement.* This paper was edited by Meric Srokosz and reviewed by two anonymous referees.

## References

Aluie, H., Hecht, M., and Vallis, G.: Mapping the Energy Cascade in the North Atlantic Ocean: The Coarse-Graining Approach, J. Phys. Oceanogr., 48, 225–244, https://doi.org/10.1175/JPO-D-17-0100.1, 2017.

Arbic, B. K., Polzin, K. L., Scott, R. B., Richman, J. G., and Shriver, J. F.: On Eddy Viscosity, Energy Cascades, and the Horizontal Resolution of Gridded Satellite Altimeter Products, J. Phys. Oceanogr., 43, 283–300, https://doi.org/10.1175/JPO-D-11-0240.1, 2013.

Arbic, B. K., Elipot, S., Brasch, J. M., Menemenlis, D., Ponte, A. L., Shriver, J. F., Yu, X., Zaron, E. D., Alford, M. H., Buijsman, M. C., Abernathey, R., Garcia, D., Guan, L., Martin, P. E., and Nelson, A. D.: Near-Surface Oceanic Kinetic Energy Distributions From Drifter Observations and Numerical Models, J. Geophys. Res., 127, 10, https://doi.org/10.1029/2022JC018551, 2022.

Ballarotta, M., Ubelmann, C., Pujol, M.-I., Taburet, G., Fournier, F., Legeais, J.-F., Faugère, Y., Delepoulle, A., Chelton, D., Dibarboure, G., and Picot, N.: On the resolutions of ocean altimetry maps, Ocean Sci., 15, 1091–1109, https://doi.org/10.5194/os-15-1091-2019, 2019.

Beal, L. M., De Ruijter, W. P. M., Biastoch, A., Zahn, R., SCOR/WCRP/IAPSO Working Group 136, Cronin, M., Hermes, J., Lutjeharms, J., Quartly, G., Tozuka, T., Baker-Yeboah, S., Bornman, T., Cipollini, P., Dijkstra, H., Hall, I., Park, W., Peeters, F., Penven, P., Ridderinkhof, H., and Zinke, J.: On the role of the Agulhas system in ocean circulation and climate, Nature, 472, 429–436, https://doi.org/10.1038/nature09983, 2011.

Bendinger, A., Cravatte, S., Gourdeau, L., Brodeau, L., Albert, A., Tchilibou, M., Lyard, F., and Vic, C.: Regional modeling of internal-tide dynamics around New Caledonia – Part 1: Coherent internal-tide characteristics and sea surface height signature, Ocean Sci., 19, 1315–1338, https://doi.org/10.5194/os-19-1315-2023, 2023.

Blanke, B., Penven, P., Roy, C., Chang, N., and Florian, K.: Ocean variability over the Agulhas Bank and its dynamical connection with the southern Benguela upwelling system, J. Geophys. Res., 114, C12028, https://doi.org/10.1029/2009JC005358, 2009.

Bourgeois, T., Goris, N., Schwinger, J., and Tjiputra, J. F.: Stratification constrains future heat and carbon uptake in the Southern Ocean between 30° S and 55° S, Nat. Commun., 13, 340, https://doi.org/10.1038/s41467-022-27979-5, 2022.

Boyd, A. F.: Physical forcing and circulation patterns on the Agulhas Bank, S. Afr. J. Sci., 90, 143–154, https://doi.org/10.10520/AJA00382353_4624, 1994.

Callies, J., Ferrari, R., Klymak, J. M., and Gula, J.: Seasonality in submesoscale turbulence, Nat. Commun., 6, 6862, https://doi.org/10.1038/ncomms7862, 2015.

Chassignet, E. P. and Xu, X.: Impact of Horizontal Resolution (1/12° to 1/50°) on Gulf Stream Separation, Penetration, and Variability, J. Phys. Oceanogr., 47, 1999–2021, https://doi.org/10.1175/JPO-D-17-0031.1, 2017.

Chaudhuri, A. H., Ponte, R. M., Forget, G., and Heimbach, P.: A comparison of atmospheric reanalysis surface products over the ocean and implications for uncertainties in air–sea boundary forcing, J. Climate, 26, 153–170, 2013.

Chelton, D.: The Wavenumber Spectra and Standard Deviations of Uncorrelated Errors in SWOT Measurements of Sea-Surface Height for Various Footprint Sizes, Tech. rep., Oregon State University, Corvallis, Oregon, https://swot.jpl.nasa.gov/system/documents/files/2253_2253-Chelton_2019_SWOT_Measurement_Noise_190523.pdf (last access: 18 April 2023), 2019.

Chelton, D., Schlax, M. G., and Samelson, R.: Global observations of nonlinear mesoscale eddies, Prog. Oceanogr., 91, 167–216, https://doi.org/10.1016/j.pocean.2011.01.002, 2011.

Chelton, D., Samelson, R., and Farrar, J.: The Effects of Uncorrelated Measurement Noise on SWOT Estimates of Sea-Surface Height, Velocity and Vorticity, J. Atmos. Ocean. Technol., 39, 1053–1083, https://doi.org/10.1175/JTECH-D-21-0167.1, 2022.

Chelton, D. B., Schlax, M. G., Samelson, R. M., and de Szoeke, R. A.: Global observations of large oceanic eddies, Geophys. Res. Lett., 34, 15, https://doi.org/10.1029/2007GL030812, 2007.

Chelton, D. B., Schlax, M. G., Samelson, R. M., Farrar, J. T., Molemaker, M. J., McWilliams, J. C., and Gula, J.: Prospects for future satellite estimation of small-scale variability of ocean surface velocity and vorticity, Prog. Oceanogr., 173, 256–350, https://doi.org/10.1016/j.pocean.2018.10.012, 2019.

Contreras, M., Renault, L., and Marchesiello, P.: Understanding Energy Pathways in the Gulf Stream, J. Phys. Oceanogr., 53, 719–736, https://doi.org/10.1175/JPO-D-22-0146.1, 2023.

Desai, S.: SWOT science requirements document, Tech. rep., Jet Propulsion Laboratory: Pasadena, CA, USA, https://swot.jpl.nasa.gov/system/documents/files/2176_2176_D-61923_SRD_Rev_B_20181113.pdf (last access 15 May 2023), 2018.

Dibarboure, G., Ubelmann, C., Flamant, B., Briol, F., Peral, E., Bracher, G., Vergara, O., Faugère, Y., Soulat, F., and Picot, N.: Data-Driven Calibration Algorithm and Pre-Launch Performance Simulations for the SWOT Mission, Remote Sens., 14, 6070, https://doi.org/10.3390/rs14236070, 2022.

d'Ovidio, F., Pascual, A., Wang, J., Doglioli, A. M., Jing, Z., Moreau, S., Grégori, G., Swart, S., Speich, S., Cyr, F., Legresy, B., Chao, Y., Fu, L., and Morrow, R. A.: Frontiers in Fine-Scale in situ Studies: Opportunities During the SWOT Fast Sampling Phase, Front. Mar. Sci., 6, 168, https://doi.org/10.3389/fmars.2019.00168, 2019.

Drushka, K., Rainville, L., and Menemenlis, D.: Internal waves and eddies from gliders and the MITgcm, https://www.aviso.altimetry.fr/fileadmin/documents/user_corner/SWOTST/SWOTST2018/Day2O_1015_Drushka_montreal_v0.pptx.pdf (last access: 29 March 2023), 2018.

Dufau, C., Orsztynowicz, M., Dibarboure, G., Morrow, R., and Le Traon, P.-Y.: Mesoscale resolution capability of altimetry: Present and future, J. Geophys. Res., 121, 4910–4927, https://doi.org/10.1002/2015JC010904, 2016.

Esteban-Fernandez, D.: SWOT Mission Performance and Error Budget; NASA/JPL Document (Reference: JPL D-79084), https://swot.jpl.nasa.gov/system/documents/files/2178_2178_SWOT_D-79084_v10Y_FINAL_REVA_06082017.pdf (last access: 4 May 2023), 2013.
Fan, L., Zhang, F., Fan, H., and Zhang, C.: Brief review of image de-noising techniques, Visual Computing for Industry, Biomedicine, and Art, 2, 7, https://doi.org/10.1186/s42492-019-0016-7, 2019.

Ferrari, R. and Wunsch, C.: Ocean Circulation Kinetic Energy: Reservoirs, Sources, and Sinks, Annu. Rev. Fluid Mech., 41, 253–282, https://doi.org/10.1146/annurev.fluid.40.111406.102139, 2008.

Forget, G., Campin, J.-M., Heimbach, P., Hill, C. N., Ponte, R. M., and Wunsch, C.: ECCO version 4: an integrated framework for non-linear inverse modeling and global ocean state estimation, Geosci. Model Dev., 8, 3071–3104, https://doi.org/10.5194/gmd-8-3071-2015, 2015.

Fu, L.-L. and Ubelmann, C.: On the Transition from Profile Altimeter to Swath Altimeter for Observing Global Ocean Surface Topography, J. Atmos. Ocean. Technol., 31, 560–568, https://doi.org/10.1175/JTECH-D-13-00109.1, 2014.

Fu, L.-L., Alsdorf, D., Morrow, R., Rodriguez, E., and Mognard, N.: SWOT: The Surface Water and Ocean Topography Mission, Tech. rep., Jet Propulsion Laboratory, California Institute of Technology Pasadena, California, https://swot.jpl.nasa.gov/system/documents/files/2179_SWOT_MSD_final-3-26-12.pdf (last access: 25 April 2023), 2012.

Gaultier, L., Ubelmann, C., and Fu, L.-L.: The Challenge of Using Future SWOT Data for Oceanic Field Reconstruction, J. Atmos. Ocean. Technol., 33, 119–126, https://doi.org/10.1175/JTECH-D-15-0160.1, 2016.

Germano, M.: Turbulence - The filtering approach, J. Fluid Mech., 238, 325–336, https://doi.org/10.1017/S0022112092001733, 1992.

Gordon, A., Weiss, R., Smethie, J., and Warner, M.: Thermocline and Intermediate Water Communication Between the South Atlantic and Indian Oceans, J. Geophys. Res., 97, 7223–7240, https://doi.org/10.1029/92JC00485, 1992.

Goschen, W. S., Bornman, T. G., Deyzel, S. H. P., and Schumann, E. H.: Coastal upwelling on the far eastern Agulhas Bank associated with large meanders in the Agulhas Current, Cont. Shelf Res., 101, 34–46, https://doi.org/10.1016/j.csr.2015.04.004, 2015.

Gula, J., Molemaker, M. J., and McWilliams, J. C.: Submesoscale Cold Filaments in the Gulf Stream, J. Phys. Oceanogr., 44, 2617–2643, https://doi.org/10.1175/JPO-D-14-0029.1, 2014.

Gula, J., Molemaker, M. J., and McWilliams, J. C.: Submesoscale Dynamics of a Gulf Stream Frontal Eddy in the South Atlantic Bight, J. Phys. Oceanogr., 46, 305–325, https://doi.org/10.1175/JPO-D-14-0258.1, 2016.

Gómez-Navarro, L., Fablet, R., Mason, E., Pascual, A., Mourre, B., Cosme, E., and Le Sommer, J.: SWOT Spatial Scales in the Western Mediterranean Sea Derived from Pseudo-Observations and an Ad Hoc Filtering, Remote Sens., 10, 4, https://doi.org/10.3390/rs10040599, 2018.

Gómez-Navarro, L., Cosme, E., Sommer, J. L., Papadakis, N., and Pascual, A.: Development of an Image De-Noising Method in Preparation for the Surface Water and Ocean Topography Satellite Mission, Remote Sens., 12, 4, https://doi.org/10.3390/rs12040734, 2020.

IPCC: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, edited by: Pörtner, H.-O., Roberts, D.C., Masson-Delmotte, V., Zhai, P., Tignor, M., Poloczanska, E., Mintenbeck, K., Alegría, A., Nicolai, M., Okem, A., Petzold, J.,

Rama, B., and Weyer, N.M.: Tech. rep., Intergovernmental Panel on Climate Change, https://www.ipcc.ch/site/assets/uploads/sites/3/2019/12/SROCC_FullReport_FINAL.pdf (last access 22 May 2023), 2019.

Jacobs, Z., Roberts, M., Jebri, F., Srokosz, M., Kelly, S., Sauer, W., Bruggeman, J., and Popova, E.: Drivers of productivity on the Agulhas Bank and the importance for marine ecosystems, Deep-Sea Res., 199, 105080, https://doi.org/10.1016/j.dsr2.2022.105080, 2022.

Johnson, G. C. and Lyman, J. M.: GOSML: A Global Ocean Surface Mixed Layer Statistical Monthly Climatology: Means, Percentiles, Skewness, and Kurtosis, J. Geophys. Res., 127, 1, https://doi.org/10.1029/2021JC018219, 2022.

Krug, M. and Penven, P.: New perspectives on Natal Pulses from satellite observations, J. Geophys. Res., 116, C7, https://doi.org/10.1029/2010JC006866, 2011.

Krug, M. and Tournadre, J.: Satellite observations of an annual cycle in the Agulhas Current, Geophys. Res. Lett., 39, 15, https://doi.org/10.1029/2012GL052335, 2012.

Krug, M., Tournadre, J., and Dufois, F.: Interactions between the Agulhas Current and the eastern margin of the Agulhas Bank, Cont. Shelf Res., 81, 67–79, https://doi.org/10.1016/j.csr.2014.02.020, 2014.

Krug, M., Swart, S., and Gula, J.: Submesoscale cyclones in the Agulhas current, Geophys. Res. Lett., 44, 346–354, https://doi.org/10.1002/2016GL071006, 2017.

Largier, J. L., Chapman, P., Peterson, W. T., and Swart, V. P.: The western Agulhas Bank: circulation, stratification and ecology, S. Afr. J. Marine Sci., 12, 319–339, https://doi.org/10.2989/02577619209504709, publisher: Taylor & Francis, 1992.

Le Guillou, F., Lahaye, N., Ubelmann, C., Metref, S., Cosme, E., Ponte, A., Le Sommer, J., Blayo, E., and Vidard, A.: Joint Estimation of Balanced Motions and Internal Tides From Future Wide-Swath Altimetry, Journal of Advances in Modeling Earth Systems, 13, e2021MS002613, https://doi.org/10.1029/2021MS002613, publisher: John Wiley & Sons, Ltd, 2021.

Leonard, A.: Energy Cascade in Large-Eddy Simulations of Turbulent Fluid Flows, in: Turbulent Diffusion in Environmental Pollution, edited by: Frenkiel, F. N. and Munn, R. E., vol. 18 of Advances in Geophysics, 237–248, Elsevier, https://doi.org/10.1016/S0065-2687(08)60464-1, 1975.

Lin, H., Liu, Z., Hu, J., Menemenlis, D., and Huang, Y.: Characterizing meso- to submesoscale features in the South China Sea, Prog. Oceanogr., 188, 102420, https://doi.org/10.1016/j-pocean.2020.102420, 2020.

Lutjeharms, J. R. E.: The Agulhas Current, Springer Berlin Heidelberg, https://doi.org/10.1007/3-540-37212-1, 2006.

Lutjeharms, J. R. E. and Gordon, A. L.: Shedding of an Agulhas ring observed at sea, Nature, 325, 138–140, https://doi.org/10.1038/325138a0, 1987.

Lévy, M., Iovino, D., Resplandy, L., Klein, P., Madec, G., Tréguier, A. M., Masson, S., and Takahashi, K.: Large-scale impacts of submesoscale dynamics on phytoplankton: Local and remote effects, Ocean Modell., 43-44, 77–93, https://doi.org/10.1016/j.ocemod.2011.12.003, 2012.

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

Martínez-Moreno, J., Hogg, A. M., Kiss, A. E., Constantinou, N. C., and Morrison, A. K.: Kinetic Energy of Eddy-Like Features From Sea Surface Altimetry, J. Adv. Model. Earth Syst., 11, 3090–3105, https://doi.org/10.1029/2019MS001769, 2019.

McWilliams, J. C.: The nature and consequences of oceanic eddies, in: Geophysical Monograph Series, edited by: Hecht, M. W. and Hasumi, H., vol. 177, 5–15, American Geophysical Union, Washington, D. C., https://doi.org/10.1029/177GM03, 2008.

Meredith, M. P. and Hogg, A. M.: Circumpolar response of Southern Ocean eddy activity to a change in the Southern Annular Mode, Geophys. Res. Lett., 33, 16, https://doi.org/10.1029/2006GL026499, 2006.

Morrow, R. and Le Traon, P.-Y.: Recent advances in observing mesoscale ocean dynamics with satellite altimetry, Adv. Space Res., 50, 1062–1076, https://doi.org/10.1016/j.asr.2011.09.033, 2012.

Morrow, R., Ward, M. L., Hogg, A. M., and Pasquet, S.: Eddy response to Southern Ocean climate modes, J. Geophys. Res., 115, C10, https://doi.org/10.1029/2009JC005894, 2010.

Morrow, R., Fu, L.-L., Ardhuin, F., Benkiran, M., Chapron, B., Cosme, E., d'Ovidio, F., Farrar, J. T., Gille, S. T., Lapeyre, G., Le Traon, P.-Y., Pascual, A., Ponte, A., Qiu, B., Rascle, N., Ubelmann, C., Wang, J., and Zaron, E. D.: Global Observations of Fine-Scale Ocean Surface Topography With the Surface Water and Ocean Topography (SWOT) Mission, Front. Mar. Sci., 6, 232, https://doi.org/10.3389/fmars.2019.00232, 2019.

Olson, D. B. and Evans, R. H.: Rings of the Agulhas current, Deep-Sea Res., 33, 27–42, https://doi.org/10.1016/0198-0149(86)90106-8, 1986.

Ponte, A.: LLC4320 surface fields [data set], https://sextant.ifremer.fr/geonetwork/srv/api/records/b2bcb9af-f335-45b6-a2a9-e460e4132879 (last access: 25 September 2023), 2020.

Renault, L., Molemaker, M. J., Gula, J., Masson, S., and McWilliams, J. C.: Control and Stabilization of the Gulf Stream by Oceanic Current Interaction with the Atmosphere, J. Phys. Oceanogr., 46, 3439–3453, https://doi.org/10.1175/JPO-D-16-0115.1, 2016.

Renault, L., McWilliams, J. C., and Penven, P.: Modulation of the Agulhas Current Retroflection and Leakage by Oceanic Current Interaction with the Atmosphere in Coupled Simulations, J. Phys. Oceanogr., 47, 2077–2100, https://doi.org/10.1175/JPO-D-16-0168.1, 2017.

Renault, L., McWilliams, J. C., and Gula, J.: Dampening of Submesoscale Currents by Air-Sea Stress Coupling in the Californian Upwelling System, Sci. Rep., 8, 13388, https://doi.org/10.1038/s41598-018-31602-3, 2018.

Renault, L., Marchesiello, P., Masson, S., and McWilliams, J. C.: Remarkable Control of Western Boundary Currents by Eddy Killing, a Mechanical Air-Sea Coupling Process, Geophys. Res. Lett., 46, 2743–2751, https://doi.org/10.1029/2018GL081211, 2019.

Rocha, C., Chereskin, T., Gille, S., and Menemenlis, D.: Mesoscale to Submesoscale Wavenumber Spectra in Drake Passage, J. Phys. Oceanogr., 46, 151222135934003, https://doi.org/10.1175/JPO-D-15-0087.1, 2015.

Rocha, C., Gille, S., Chereskin, T., and Menemenlis, D.: Seasonality of submesoscale dynamics in the Kuroshio Extension, Geophys. Res. Lett., 43, 601–620, https://doi.org/10.1002/2016GL071349, 2016.

Rodriguez, E., Fernadez, D., Peral, E., Chen, C., Bleser, J.-W., and Williams, B.: Wide-Swath Altimetry: A Review, 71–112, CRC Press, https://doi.org/10.1201/9781315151779-2, 2017.

Ruijter, W. P. M. d., Leeuwen, P. J. v., and Lutjeharms, J. R. E.: Generation and Evolution of Natal Pulses: Solitary Meanders in the Agulhas Current, J. Phys. Oceanogr., 29, 3043–3055, https://doi.org/10.1175/1520-0485(1999)029<3043:GAEONP>2.0.CO;2, 1999.

Sallée, J.-B.: Southern Ocean Warming, Oceanography, 31, 52–62, https://www.jstor.org/stable/26542651 (last access: 31 March 2023), 2018.

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

Savage, A. C., Arbic, B. K., Alford, M. H., Ansong, J. K., Farrar, J. T., Menemenlis, D., O'Rourke, A. K., Richman, J. G., Shriver, J. F., Voet, G., Wallcraft, A. J., and Zamudio, L.: Spectral decomposition of internal gravity wave sea surface height in global models, J. Geophys. Res., 122, 7803–7821, https://doi.org/10.1002/2017JC013009, 2017.

Schouten, M. W., de Ruijter, W. P. M., and van Leeuwen, P. J.: Upstream control of Agulhas Ring shedding, J. Geophys. Res., 107, 23–1–23–11, https://doi.org/10.1029/2001JC000804, 2002.

Schubert, R., Gula, J., Greatbatch, R., Baschek, B., and Biastoch, A.: The Submesoscale Kinetic Energy Cascade: Mesoscale Absorption of Submesoscale Mixed-Layer Eddies and Frontal Downscale Fluxes, J. Phys. Oceanogr., 50, 2573–2589, https://doi.org/10.1175/JPO-D-19-0311.1, 2020.

Sebille, E. V. and Leeuwen, P. J. V.: Fast Northward Energy Transfer in the Atlantic due to Agulhas Rings, J. Phys. Oceanogr., 37, 2305–2315, https://doi.org/10.1175/JPO3108.1, 2007.

Siegelman, L., Klein, P., Rivière, P., Thompson, A. F., Torres, H. S., Flexas, M., and Menemenlis, D.: Enhanced upward heat transport at deep submesoscale ocean fronts, Nat. Geosci., 13, 50–55, https://doi.org/10.1038/s41561-019-0489-1, 2020.

Sinha, A. and Abernathey, R. P.: Time Scales of Southern Ocean Eddy Equilibration, J. Phys. Oceanogr., 46, 2785–2805, https://doi.org/10.1175/JPO-D-16-0041.1, 2016.

Speich, S., Arhan, M., Ansorge, I., Boebel, O., Sokov, A., Gladyshev, S., Farbach, E., Byrne, D., Klepikov, A., Garzoli, S., and Rodriguez, M. A.: Good-Hope/Southern Ocean: A study and monitoring of the Indo-Atlantic connections, Tech. Rep. 27, Mercator Newsletter, https://www.coriolis.eu.org/content/download/729/4975/file/2005_Speich_OceanAustral_NewsMerc07.pdf (last access 11 June 2023), 2007.

Su, Z., Wang, J., Klein, P., Thompson, A. F., and Menemenlis, D.: Ocean submesoscales as a key component of the global heat budget, Nat. Commun., 9, 775, https://doi.org/10.1038/s41467-018-02983-w, 2018.

Swart, S., Speich, S., Ansorge, I. J., Goni, G. J., Gladyshev, S., and Lutjeharms, J. R. E.: Transport and variability of the Antarctic Circumpolar Current south of Africa, J. Geophys. Res., 113, C9, https://doi.org/10.1029/2007JC004223, 2008.
Taburet, G., Sanchez-Roman, A., Ballarotta, M., Pujol, M.-I., Leg-
eais, J.-F., Fournier, F., Faugere, Y., and Dibarboure, G.: DUACS
DT2018: 25 years of reprocessed sea level altimetry products,
Ocean Sci., 15, 1207–1224, https://doi.org/10.5194/os-15-1207-
2019, 2019.

Tedesco, P., Gula, J., Ménesguen, C., Penven, P., and Krug,
M.: Generation of Submesoscale Frontal Eddies in the
Agulhas Current, J. Geophys. Res., 124, 7606–7625,
https://doi.org/10.1029/2019JC015229, 2019.

Thomas, L. N., Tandon, A., and Mahadevan, A.: Submesoscale processes and dynamics, in: Geophysical Monograph Series, edited by: Hecht, M. W. and Hasumi, H., Vol. 177, 17–38, American Geophysical Union, Washington, D. C., https://doi.org/10.1029/177GM04, 2008.

Torres, H. S., Klein, P., Menemenlis, D., Qiu, B., Su, Z., Wang, J., Chen, S., and Fu, L.-L.: Partitioning Ocean Motions Into Balanced Motions and Internal Gravity Waves: A Modeling Study in Anticipation of Future Space Missions, J. Geophys. Res., 123, 8084–8105, https://doi.org/10.1029/2018JC014438, 2018.

Tréboutte, A., Carli, E., Ballarotta, M., Carpentier, B., Faugère, Y., and Dibarboure, G.: KaRIn Noise Reduction Using a Convolutional Neural Network for the SWOT Ocean Products, Remote Sens., 15, 2183, https://doi.org/10.3390/rs15082183, 2023.

Ubelmann, C., Gaultier, L., Fu, L.-L., and Briol, F.: SWOT Sim-
ulator for Ocean Science, Jet Propulsion Laboratory, Califor-
nia Institute of Technology, CNES [data set], https://github.com/
CNES/swot_simulator (last access: 27 February 2023), 2021.

Vergara, O., Morrow, R., Pujol, M.-I., Dibarboure, G., and Ubelmann, C.: Global submesoscale diagnosis using along-track satellite altimetry, Ocean Sci., 19, 363–379, https://doi.org/10.5194/os-19-363-2023, 2023.

Verron, J., Bonnefond, P., Andersen, O., Ardhuin, F., Bergé-
Nguyen, M., Bhowmick, S., Blumstein, D., Boy, F., Brodeau,
L., Cretaux, J., Dabat, M., Gerald, D., Fleury, S., Garnier, F.,
Gourdeau, L., Marks, K., Queruel, N., Sandwell, D., Smith,
W., and Zaron, E.: The SARAL/AltiKa mission: A step for-
ward to the future of altimetry, Adv. Space Res., 68, 808–828,
https://doi.org/10.1016/j.asr.2020.01.030, 2020.

Volkov, D. L., Fu, L.-L., and Lee, T.: Mechanisms of the meridional heat transport in the Southern Ocean, Ocean Dynam., 60, 791–801, https://doi.org/10.1007/s10236-010-0288-0, 2010.

Wang, J., Fu, L.-L., Qiu, B., Menemenlis, D., Farrar, J. T., Chao, Y., Thompson, A. F., and Flexas, M. M.: An Observing System Simulation Experiment for the Calibration and Validation of the Surface Water Ocean Topography Sea Surface Height Measurement Using In Situ Platforms, J. Atmos. Ocean. Technol., 35, 281–297, https://doi.org/10.1175/JTECH-D-17-0076.1, 2018.

Wang, J., Fu, L.-L., Torres, H. S., Chen, S., Qiu, B., and Menemen-
lis, D.: On the Spatial Scales to be Resolved by the Surface Water
and Ocean Topography Ka-Band Radar Interferometer, J. Atmos.
Ocean. Technol., 36, 87–99, https://doi.org/10.1175/JTECH-D-
18-0119.1, 2019.

Wirth, A.: A Guided Tour Through Physical Oceanography. Master. Physical Oceanography, France, France, cel-01134110v4, HAL Id: cel-01134110, https://hal.archives-ouvertes.fr/cel-01134110v4 (last access: April 2023), 2015.

Wunsch, C.: What Is the Thermohaline Circulation?, Science, 298, 1179–1181, https://doi.org/10.1126/science.1079329, 2002.

Zaron, E. D. and Ray, R. D.: Using an altimeter-derived internal tide model to remove tides from in situ data, Geophys. Res. Lett., 44, 4241–4245, https://doi.org/10.1002/2017GL072950, 2017.

Zhang, Z., Qiu, B., Klein, P., and Travis, S.: The influ-
ence of geostrophic strain on oceanic ageostrophic mo-
tion and surface chlorophyll, Nat. Commun., 10, 2838,
https://doi.org/10.1038/s41467-019-10883-w, 2019.