# Modeling Indian Ocean circulation to study marine debris dispersion: insights into high-resolution and wave forcing effects with Symphonie 3.6.6

Lisa Weiss, Marine Herrmann, Patrick Marsaleix, Matthieu Bompoil, Christophe Maes

► To cite this version:

Lisa Weiss, Marine Herrmann, Patrick Marsaleix, Matthieu Bompoil, Christophe Maes. Modeling Indian Ocean circulation to study marine debris dispersion: insights into high-resolution and wave forcing effects with Symphonie 3.6.6. Geoscientific Model Development, 2026, 19 (2), pp.827-865. 10.5194/gmd-19-827-2026. hal-05480635

HAL Id: hal-05480635

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

Submitted on 29 Jan 2026

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 CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

# Modeling Indian Ocean circulation to study marine debris dispersion: insights into high-resolution and wave forcing effects with Symphonie 3.6.6

Lisa Weiss<fnref n="1" />,<fnref n="2" />, Marine Herrmann<fnref n="1" />, Patrick Marsaleix<fnref n="1" />, Matthieu Bompoil<fnref n="2" />, and Christophe Maes<fnref n="2" />

**Correspondence:** Lisa Weiss (lisa.weiss@ird.fr)

Received: 22 April 2025 – Discussion started: 15 May 2025

Revised: 1 December 2025 – Accepted: 10 January 2026 – Published: 26 January 2026

**Abstract.** The Indian Ocean basin faces significant anthropogenic pressure due to its connection to over 2.2 billion people through river basins. Indian Ocean dynamics are characterized by strong regional and seasonal variability driven by the monsoon system and intense eddy activity. To support future studies on land-sea transfers and marine debris dispersion in this complex ocean, we developed a new circulation modeling configuration using the hydrodynamic model SYMPHONIE. Our configuration introduces a unique telescopic grid covering the entire basin, enabling the study of sub-basin connectivity while resolving meso and submesoscale processes in the coastal region, from the Mozambique Channel to the Bay of Bengal, at a resolution of 1 to 3 km. Additionally, we integrate the recently released high-resolution GloFAS river discharge dataset to force the physical simulations with daily freshwater inputs. Three annual experiments are conducted, exploring the effects of wave forcing and spatial resolution. The reference simulation (IndOc.HR) uses high resolution without wave forcing, IndOc.HR-Sto uses the same resolution but includes wave forcing through both the Stokes-Coriolis force and the Stokes drift, and IndOc.12 is a lower-resolution simulation without wave forcing. Comparisons of temperature, salinity and sea level with in situ and satellite data show the good performance of the simulations and the ability of the high resolution model to accurately capture the spatial and temporal variability of surface dynamics and water masses over the Indian Ocean. We further analyze energy budgets and perform Lagrangian experiments to illustrate the critical role of resolution and waves in shaping the circulation and the resulting

marine debris dispersion patterns. The effect of energy levels is particularly significant on trajectory statistics such as average travel distances and preferred spread direction. Notably, Stokes drift has a significant seasonal effect in the Arabian Sea during the southwest monsoon, while current field resolution strongly influences trajectories in the Mozambique Channel. Our results provide a robust modeling framework for studying Indian Ocean dynamics and exploring their effect on marine connectivity and the transport of matter, including pollutants, larvae or organic matter.

## 1 Introduction

The Indian Ocean is characterized by complex atmospheric and oceanic dynamics that influence marine debris dispersion and connectivity (Pattiaratchi et al., 2022). Its circulation is marked by a strong regional variability and subject to rapid warming and marine heatwaves (Phillips et al., 2021). In the northern hemisphere, the reversal of the monsoon-driven currents between seasons results in significant exchanges between the Arabian Sea and the Bay of Bengal, with potentially high stranding rates of marine debris along the coasts (van der Mheen et al., 2020). In the south, large quantities of debris are trapped in the subtropical gyre (Thibault et al., 2024) connected to the other southern hemisphere gyres of the Atlantic and Pacific Oceans through the Agulhas and Tasman leakages (Maes et al., 2018; Qu et al., 2019; Dobler et al., 2019). However, the general patterns of dispersion and accumulation in this basin remain debated and under-

<footnote n="1">¹Université de Toulouse, LEGOS (IRD, CNRS, CNES, UT3), OMP, F31400, Toulouse, France</footnote>

<footnote n="2">²Université de Brest, LOPS (IRD, CNRS, Ifremer), IUEM, F29280, Plouzané, France</footnote>
studied (Pattiaratchi et al., 2022). In addition to this large-scale regional circulation, buoyant materials are concentrated by horizontal density fronts and mesoscale cyclonic activity (D'Asaro et al., 2018; Vic et al., 2022). Such frontal zones generate submesoscale structures and filaments, which actively converge drifting particles towards highly concentrated areas (Lodise et al., 2020; Weiss et al., 2024a). Submesoscale convergences are associated with higher vertical velocities than those found in mesoscale fronts, potentially enhancing mixing and accelerating spreading (D'Asaro et al., 2018; Lodise et al., 2020). But submesoscale features are ephemeral, and the clustering of marine debris in such structures may occur over a few to 10 d before slowly dissipating (D'Asaro et al., 2018). In a contamination monitoring context, understanding and interpreting convergence phenomena is challenging due to the short timescales involved. This is especially true in the case of plastic pollution, where the complexity of sampling and quantification makes it currently very difficult to verify the model-based hypotheses. Due to the scarcity of plastic debris concentration data in the ocean, it is currently impossible to capture the spatial and temporal variability of accumulation processes and concentration gradients associated with meso and submesoscales from observations.

Lagrangian dispersion modeling studies can help to fill this knowledge gap. They need in particular to be based on realistic current fields with grid resolutions permitting the representation of such structures and processes at different scales. In recent studies however, the finest horizontal resolutions achieved in ocean simulations over the entire Indian Basin available for Lagrangian scenarios range from 1/4 to 1/12° (~25 to 9 km), based on various configurations of global hydrodynamic models and products: the 1/12° CMEMS Global Ocean Physics Reanalysis GLORYS12V1 (Lellouche et al., 2021), the 1/12° Mercator Ocean NEMO-PSY4 analysis (Gasparin et al., 2018), the 1/12° HYCOM global product (0.08° zonal and 0.08 to 0.04° meridional resolution) (Cummings and Smedstad, 2013; Pottapinjara and Joseph, 2022; used by Chassignet et al., 2021 and van der Mheen et al., 2020 for Lagrangian plastic tracking), a 1/12° Indian Ocean child domain of ROMS product (designed and used for larval dispersal by Gamoyo et al., 2019), the 1/4° global ocean ECCO LLC270 product based on MITgcm (Kersalé et al., 2022) or the 1/4° ESA GlobCurrent-v3 product (Rio et al., 2014; used by van der Mheen et al., 2019 to study the Indian Ocean garbage patch). For shorter simulation periods, much higher resolution is available with the global LLC4320 configuration (1/48° ~ 2 km) for the year 2011 (Forget et al., 2015; Rocha et al., 2016). Limited coastal regions of the Southwest Indian Ocean are modeled with CROCO, with resolution reaching ~2 km over the northern part of Madagascar for a configuration developed to study coral reefs connectivity (Vogt-Vincent and Johnson, 2023) and 4.5–6 to 0.75 km for a configuration using a nesting approach

to study the Agulhas Current (Renault et al., 2017; Tedesco et al., 2019).

Several other factors are important to consider when modeling particles dispersion. River forcing is of course an essential factor in the study of marine debris dispersion in coastal environment (Weiss et al., 2021; Laverre et al., 2023; Weiss et al., 2024a). It indeed influences the amount of debris released as well as the submesoscale circulation that develops along the fronts in the coastal zone, which in turn may affect the Lagrangian dispersion statistics (Choi et al., 2017). Beside the grid resolution mentioned above, the atmospheric forcing influences the energetic production of submesoscale activity and thus the particles pathways and distribution (van Sebille et al., 2020; Weiss et al., 2024a). Wind and waves processes can significantly influence the dispersion and accumulation of buoyant material such as marine debris (van Sebille et al., 2020). Their impact occurs through multiple mechanisms (Onink et al., 2019): (i) the modification of Eulerian currents by waves, notably via the Stokes–Coriolis force, (ii) the direct advection of particles by Stokes drift when explicitly included in Lagrangian schemes, and (iii) the windage effect, i.e., the direct wind forcing on partially emerged objects. Different levels of complexity in representing wave effects have also been explored in the literature, from simplified parameterizations through one-way forcing to more comprehensive schemes that include multiple wave-induced processes and/or two-way coupling frameworks (Röhrs et al., 2012; Cunningham et al., 2022; Rühs et al., 2025; Bajon et al., 2023; Couvelard et al., 2020). However, only a few of the above-mentioned ocean model configurations for the Indian ocean account for all those parameters likely to influence the basin dynamics and Lagrangian dispersion patterns (van Sebille et al., 2020). None of them include wave forcing to advect drifting particle at the basin scale. While high resolution CROCO configurations (Vogt-Vincent and Johnson, 2023; Renault et al., 2017) integrate realistic atmospheric and tidal forcings, they don't include high-frequency freshwater discharges or wave effect, and their spatial coverage does not cover the entire basin. In particular, they do not account for the large number of daily freshwater inflows over the basin and do not use grid resolution and coverage allowing to address both small to mesoscale coastal dynamics and large scale circulation to study the basin-wide connectivity. One of the main challenge in modeling the dispersion of marine debris is indeed the necessity to study the continuity of dispersion patterns from small coastal scales to mesoscale structures to large scale ocean currents that govern connectivity and accumulation in the open ocean.

In this context, our aim is to build an ocean model framework to investigate scenarios of marine debris dispersion from the coastal zone – where the main sources of debris are located and where plastic stranding poses significant socio-environmental challenges – to the open basin. While the present study does not explicitly address plastic land-sea transfers, it provides the hydrodynamic framework necessary
for future analyses focusing on riverine inputs and coastal retention. The specific objectives of this paper are to (1) provide a state-of-the-art numerical representation of the Indian ocean dynamics, integrating updated atmospheric, tidal, riverine and wave forcings, and addressing the submesoscale horizontal resolution; (2) analyze the impact of resolution and wave forcing on seasonal energy budget; and (3) quantify and illustrate the variability in associated Lagrangian trajectory statistics. For this, we developed a new configuration of the SYMPHONIE hydrodynamic model covering the entire Indian Ocean with a stretched grid allowing a kilometric resolution in the coastal areas. We ran a reference and two sensitivity annual ocean simulations, investigating in particular the effect of wave forcing and resolution on dynamics and Lagrangian dispersion. The paper first presents the model configurations developed, the simulations performed and the datasets used for forcing and evaluation. Second we evaluate and compare the ability of those different simulations to reproduce the water masses and ocean dynamics characteristics over the studied area. Third, we evaluate the seasonal surface kinetic energy budget in each simulation to assess the sensitivity to grid resolution and wave forcing and to better understand the underlying dynamical driver. Finally, we present and analyze Lagrangian dispersion experiments to study buoyant particle trajectories and quantify the effect of the different factors considered here.

# 2 Methodology and data

## 2.1 Model setup

### 2.1.1 The SYMPHONIE model

To set up the new numerical configuration of ocean dynamics in the Indian Ocean, we use the 3D hydrodynamic model SYMPHONIE-v3.6.6 (developed by the SIROCCO service; Marsaleix et al., 2008, 2019). The model, based on Boussinesq approximation and hydrostatic equilibrium, solves the Navier-Stokes primitive equations for momentum, temperature, salinity and volume conservation. It is described and used in several recent regional applications: for the study of regional circulation and water, heat and salt budget (Trinh et al., 2024; Garinet et al., 2024) and upwelling (To Duy et al., 2022; Herrmann et al., 2023, 2024) in the South China Sea and Southeast Asian seas, and for the study of sediment dynamics (Estournel et al., 2023), organic carbon budget (Habib et al., 2023) and Lagrangian dispersion of larvae (Guizien et al., 2012) or microplastics (Weiss et al., 2024b, a) in the Mediterranean Sea.

### 2.1.2 Grids

The IndOc.HR configuration is based on a telescopic bipolar grid (Bentsen et al., 1999) (Fig. 1) that has been developed for the specific needs of our main objectives, with a view to

socio-ecological applications within the framework of multi-disciplinary regional collaborations. The grid is a curvilinear Arakawa C-grid (3000 × 2800 cells) covering the entire Indian Ocean from the Agulhas Current to the Southeast Asian seas and as far south as Antarctica, thereby limiting the number of open boundaries. It has a strong refinement toward the East African and Indian coasts, which are the two highly resolved poles with a horizontal resolution of 1 to 2 km. It allows to simulate the large-scale circulation throughout the whole basin while resolving the small-scale processes as eddies and filaments in the areas of interest. Vertically, the grid is defined by 60 vanishing quasi-sigma (VQS) vertical levels as described by Estournel et al. (2021) for the Mediterranean and Trinh et al. (2024) for the South China Sea. Bathymetry is built from an interpolation of the GEBCO-2021 product. Depth ranges from 1 to -7459 m in the domain.

The lower resolution configuration (IndOc.12) uses the same features as IndOc.HR, except for its horizontal discretization, which is regular at 1/12°.

Three regions of interest with contrasting ocean dynamics, corresponding to the best-resolved areas of the IndOc.HR coastal domain, are particularly analyzed in this study: the Arabian Sea (hereafter referred to as ArS), the Bay of Bengal (hereafter referred to as BoB) and the Mozambique Channel (hereafter referred to as MoC). These regions are used to calculate spatial averages of temperature, salinity, sea level and energy for the evaluation of the simulations (larger blue boxes on Fig. 1 used for Fig. 6 and smaller blue boxes on Fig. 1 used for Figs. 7 and 9).

### 2.1.3 Atmospheric and open ocean boundary conditions

Atmospheric forcing on the free surface boundary is calculated from the bulk formulae of Large and Yeager (2004), using wind stress, pressure, solar/thermal radiations, humidity and temperature provided by the ECMWF's ERA5 daily mean product at 0.25°. Open lateral boundary conditions and initial state for currents, temperature, salinity and sea surface elevation are provided by the daily mean CMEMS GLORYS12V1 product at 1/12° (Lellouche et al., 2021). A temporal nudging layer is configured in regions where the resolution of the telescopic grid is coarser than that of the 1/12° GLORYS forcing. This only occurs is the southeastern part of the domain, south of 50°S and around Australia, as shown by the grid resolution in Fig. 1. In these areas, temperature and salinity fields are relaxed toward their better-resolved GLORYS counterpart, over the full water column, using a relaxation time scale that decreases smoothly with the difference in resolution between the two grids. Relaxation is therefore greater at the southeastern edge of the domain, with a minimum time scale of 30 d. This nudging does not impact the three regions selected for model evaluation (ArS, BoB, MoC), ensuring that it does not artificially improve the agreement with observations in these areas. Tidal forcing is
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Model telescopic grid resolution (in km, gray scale), displaying every 50 gridlines (red lines). The south boundary reaches Antarc-
tica with a maximum mesh size of 23 km. Blue shaded color-bar and dots show the 336 daily freshwater discharges simulated in the config-
uration (the size of the dots is proportional to their mean annual discharge over the period 2016–2021). The names of the domain’s largest
rivers are displayed on the map. Blue boxes in the best-resolved regions (ArS: Arabian Sea, BoB: Bay of Bengal, MoC: Mozambique Chan-
nel) are used to calculate spatial averages of temperature, salinity, sea elevation, energy spectra and Lagrangian trajectories (larger boxes for
Fig. 6 and Table 2; smaller boxes for Figs. 2, 7, 9 and 10). The six yellow dots show the RAMA buoy locations for comparisons with in situ
profiles (Fig. A2). The three pink circles show the particle release locations for the Lagrangian sensitivity analysis.</figcaption></figure>

prescribed from the FES2014 global tidal model (Lyard et al.,
2021), taking into account 9 barotropic components in phase
and amplitude according to Pairaud et al. (2008, 2010).

**2.1.4 Rivers**

Our domain includes 336 rivers with a mean annual dis-
charge > 100 m³ s⁻¹ (Fig. 1), for a total annual discharge
into the domain around 9400 km³ yr⁻¹ (averaged over 2016–
2021). Those discharges are extracted from the gridded
CMEMS Global Flood Awareness System (GloFAS-v3.1;
Zsoter et al., 2021) forecasts characterized by unprecedented
temporal (daily) and spatial (0.1°) resolutions. It combines
satellite data, models (ECMWF and LISFLOOD) and in
situ measurements. We identify river mouths locations by
reconstructing the watercourses down to the coastal river
outlets where water flows are at their maximum. When a
river mouth corresponds to a delta, we ignore the differ-
ent branches of the delta and simulate a single inflow point
with a discharge corresponding to the sum of all the streams
(as for the Ganges-Brahmaputra Delta for example). Among
the 336 extracted rivers, 4 rivers have a mean annual dis-
charge > 10 000 m³ s⁻¹ (Ganges-Brahmaputra-Meghna, Ir-
rawaddy, Mekong, Zambezi), 11 rivers > 5000 m³ s⁻¹ (In-
dus, Salouen, Godavari, Kapuas, Chatt el Arab, Pulau, Ma-
hakam), 40 rivers > 1000 m³ s⁻¹ and 92 rivers > 500 m³ s⁻¹.

Each river daily inflow is introduced at the grid point clos-
est to the river mouth, using a hybrid Neumann/Dirichlet
boundary condition constrained to respect an imposed total
flow, described by Nguyen-Duy et al. (2021).

**2.1.5 Current-wave interaction**

In consideration of the recommendations put forth by van
Sebille et al. (2020) concerning the impact of waves on
transport, at the large scales and relatively low grid resolu-
tions considered here, the model incorporates a simplified
parametrization for the influence of waves on currents, re-
taining only the dominant terms in the right-hand sides of the
momentum Eqs. (2) and (3) below, i.e., the Stokes-Coriolis
force (van Sebille et al., 2020; Jordà et al., 2007; Xu and
Bowen, 1994). Besides, readers' attention is drawn to the fact
that a similar but more precise parameterization with more
terms is provided by Couvelard et al. (2020).

Here, the parametrization consists in an online one-way
forcing utilizing the three-component Stokes drift vector
(ûs, ûs, ûs) in m s⁻¹, extrapolated below the surface us-
ing the peak period, sourced in our study from the global
CMEMS WAVERYS 0.2° product (CMEMS, 2019). The
general wave-current equations of the SYMPHONIE model
distinguish between the components (u, v, w) of the Eulerian
current and the components (ûs, ûs, ûs) of the Stokes drift.
The Stokes drift enters the Eulerian equations through the
advection operator. The material derivative used in Eqs. (2) to (5) below includes the advection by the Lagrangian velocity ($u + \tilde{u}^s, v + \tilde{v}^s, w + \tilde{w}^s$), i.e. Eulerian current + Stokes drift, such as:

$$D_t = \partial_t + \partial_x [(u + \tilde{u}^s).] + \partial_y [(v + \tilde{v}^s).] + \partial_z [(w + \tilde{w}^s).] \quad (1)$$

The horizontal components of the Eulerian current are given by:

$$D_t u = +f(v + \tilde{v}^s) - \frac{1}{\rho_0} \frac{\partial p}{\partial x} + \frac{\partial}{\partial z} \left( K_m \frac{\partial u}{\partial z} \right) + \Delta^2 u \quad (2)$$

$$D_t v = -f(u + \tilde{u}^s) - \frac{1}{\rho_0} \frac{\partial p}{\partial y} + \frac{\partial}{\partial z} \left( K_m \frac{\partial v}{\partial z} \right) + \Delta^2 v \quad (3)$$

where $f$ is the Coriolis parameter, $\rho_0$ is the reference seawater density, $p$ is the pressure, and $\Delta^2$ is a horizontal bi-laplacian operator. Then, the equations for temperature $T$ and salinity $S$ are defined such as:

$$D_t T = \frac{\partial}{\partial z} K_h \frac{\partial T}{\partial z} + \Delta T + \frac{\partial}{\partial z} I \quad (4)$$

$$D_t S = \frac{\partial}{\partial z} K_h \frac{\partial S}{\partial z} + \Delta S \quad (5)$$

with $K_m$ and $K_h$ the vertical turbulent viscosity and diffusivity coefficients for momentum and tracers, calculated by the turbulent closure scheme ($k-\epsilon$), $I$ the radiative forcing term and $\Delta$ a horizontal diffusion operator.

The sea surface height $\eta$ is given by:

$$\partial_t \eta = -\partial_x \int_{-h}^{\eta} (u + \tilde{u}^s) dz - \partial_y \int_{-h}^{\eta} (v + \tilde{v}^s) dz + P - E \quad (6)$$

where $P-E$ is the precipitation and evaporation budget. The continuity equation is applied to the Lagrangian velocity, defined as the sum of the Eulerian current and the Stokes drift, to deduce its vertical component, as follow:

$$\partial_x (u + \tilde{u}^s) + \partial_y (v + \tilde{v}^s) + \partial_z (w + \tilde{w}^s) = 0 \quad (7)$$

Eulerian tracers $T$ and $S$ are therefore both displaced consistently by the Lagrangian velocity, ensuring coherence between Eulerian and Lagrangian transport formulations.

## 2.2 The simulations

### 2.2.1 Ocean simulations

Using the configuration and grid described above, we produced two twin simulations: a reference simulation, called IndOc.HR in the following and not including wave effects; and its twin simulation, called IndOc.HR-Sto in the following, including the effect of waves on Eulerian currents via the Stokes-Coriolis term in the momentum equations and via

the explicit advection by Stokes drift (see Sect. 2.1.5). To estimate the potential added-value of the telescopic grid, we also produced a simulation called IndOc.12 in the following, using the same forcings and settings as IndOc.HR, but on a regular horizontal grid covering the same domain at 1/12° resolution. Those three simulations are performed over the year 2017 after a two months spin up. Except for the wave effect, they are all subjected to the same forcings. The simulations are performed with a calculation time step of 200 s. Daily average outputs are stored, including current velocities, temperature, salinity, sea surface height and turbulent kinetic energy.

### 2.2.2 Lagrangian simulations

Three Lagrangian numerical experiments are performed offline using the SYMPHONIE Lagrangian tool (as described in Weiss et al., 2024b), each using the daily 3D velocity fields produced by the three simulations IndOc.HR, IndOc.HR-Sto and IndOc.12 (Sect. 2.2). The particle advection is based on a Runge-Kutta second order advection scheme with an integration time step of the particle trajectories of 10 min. The current velocities are interpolated linearly in space (3D) and time at each particle position. The time derivative of the 3D position vector $X(t) = (x, y, z)(t)$ is equal to the stored Lagrangian velocities provided by each hydrodynamic simulation.

- In the classical case (IndOc.HR and IndOc.12), the Lagrangian velocity corresponds to the resolved Eulerian currents $(u, v, w)$ only:

$$\frac{dX(t)}{dt} = \begin{cases} \frac{dx(t)}{dt} = u(t) \\ \frac{dy(t)}{dt} = v(t) \\ \frac{dz(t)}{dt} = w(t) + w_t(t) + w_s(t) \end{cases} \quad (8)$$

- In the wave-modified case (IndOc.HR-Sto), the Lagrangian velocity includes the Eulerian currents $(u, v, w)$ dynamically modified by the Stokes-Coriolis force and the Stokes drift $(\tilde{u}^s, \tilde{v}^s, \tilde{w}^s)$ as shown in Sect. 2.1.5, such as:

$$\frac{dX(t)}{dt} = \begin{cases} \frac{dx(t)}{dt} = u(t) + \tilde{u}^s(t) \\ \frac{dy(t)}{dt} = v(t) + \tilde{v}^s(t) \\ \frac{dz(t)}{dt} = w(t) + \tilde{w}^s(t) + w_t(t) + w_s(t) \end{cases} \quad (9)$$

In all cases (Eqs. 8 and 9), the vertical displacement of the particles considers (i) a vertical turbulent random-walk diffusion term $w_t$ based on the turbulent kinetic energy in the particle grid cell and a random number (Guizien et al., 2012), and (ii) a specific vertical velocity term $w_s$ reflecting the particle buoyancy (Weiss et al., 2024b). Only buoyant particles are considered here (illustrated Sect. 5) with a rising velocity $w_s = 1 \text{ mm s}^{-1}$, oriented toward the sea surface. This value is within the lower range reported in the literature
for specific plastic debris types such as light fibers, weath-
ered or biofouled items (Waldschläger et al., 2020; Fischer
et al., 2022; Kuizenga et al., 2022). It was chosen as an inter-
mediate order of magnitude of average vertical velocities of
surface currents, allowing such buoyant particles to be influ-
ences by vertical mixing processes in the surface layer rather
than remaining strictly at the surface.

Each Lagrangian experiment consists in releasing buoyant
particles within each of the three high-resolution zones an-
alyzed in this study, i.e. ArS, BoB and MoC. In each zone,
particles are randomly released at the sea surface, between
-1 and -2 m depth, within a circular emission area of 2°
radius centered on the predefined central point of the zone
(see Fig. 2). It results in three distinct emission regions, one
per zone, for each experiment. One particle is released per
minute on the first day of each month, followed by a 30 d drift
period, resulting in 1440 trajectories per month per zone.
Over a one-year period, a total of 51 840 particles are re-
leased in each Lagrangian experiment.

## 2.3 The observations

To evaluate the simulations, we use data-sets of Sea Surface
Temperature (SST), Sea Surface Salinity (SSS), Sea Level
Anomaly (SLA), temperature and salinity vertical profiles,
and river discharges.

2.3.1 Satellite datasets

The CMEMS Level 4 Operational Sea Surface Temperature and Ice Analysis (OSTIA) product (CMEMS, 2015), providing daily gap-free maps at resolution 0.05° built from in situ (drifting and moored buoys) and remote infrared (multiple spatial sensors) observations, is used to evaluate the SST. The ESA Climate Change Initiative (CCI) weekly product at 0.5° resolution (Boutin et al., 2021) combining observations from the Soil Moisture and Ocean Salinity (SMOS) and Aquarius and Soil Moisture Active Passive (SMAP) spatial sensors is used to evaluate the SSS. Global satellite altimetry gridded Level 4 data from the CMEMS platform (CMEMS, 2021) at resolution 0.25° is used to evaluate the SLA. For comparison purposes, we remove from the simulated and observed daily SLA timeseries their temporal average over the period of comparison (2017) (see Eq. A1).

2.3.2 In situ data

*In situ* data from six buoys from the NOAA Research Moored Array for African-Asian-Australian Monsoon Analysis and Prediction (RAMA) project (McPhaden et al., 2009) are used to evaluate the SST and SSS. The CMEMS GLORYS12V1 reanalysis product (CMEMS, 2018), assimilates Argo in situ vertical profiles, with a large number of Argo floats continuously active during the year 2017 in the studied domain (around 700, leading to more than 25 000 profiles). This

product is thus also used to evaluate the simulated temper-
ature and salinity vertical profiles.

2.3.3 Hydrology datasets

The evaluation of rivers discharge is performed thanks to two
global hydrology datasets. First, the HydroSHEDS atlas, at
15 arcsec resolution, i.e. 1/240° or ~0.5 km at the equator
(Linke et al., 2019, long-term average over 1971-2000), is
used to evaluate annual mean discharges and river catchment
areas. Second, the Dai & Trenberth (D&T) climatology (Dai,
2017, based on scattered measurements between 1900 and
2018), is used to evaluate annual and monthly climatological
mean discharges for the 47 individual rivers common to our
domain and their dataset.

2.4 Analysis of river discharge

In this section, we examine the representation of rivers, com-
paring daily and annual discharges from GloFASv3.1 with
HydroSHEDS, D&T climatology and the updated GloFASv4
datasets for some common rivers (Figs. 3, 4 and A4).

The mean annual discharge of the 336 rivers considered in
the IndOc domain (Indian Ocean and part of the South-East
Asian seas) is estimated at 297 464 and 216 295 m³ s⁻¹ by
GLOFASv3.1 and HydroSHEDS, respectively (i.e. +37% in GLOFASv3.1).
For the 47 rivers also listed in the D&T monthly climatology, the correlation coefficients between the
three datasets over the year 2017 are nevertheless highly
significant (> 0.96). The GLOFASv3.1 NRMSE calculated
from the mean discharges of the 47 rivers is 6% with a mean
bias of 946 m³ s⁻¹ compared to HydroSHEDS, and 24% with a mean bias of 1693 m³ s⁻¹ compared to D&T (Fig. 4).
Most rivers with a mean annual discharge bias higher than
100% are located on the East African coast, such as Jubba
in Somalia, Pangani, Ruvu and Rufiji in Tanzania, Zambezi
and Limpopo in Mozambique, Morondava in Madagascar,
and Tana in Kenya (Southern Hemisphere in Fig. 3). On the
contrary, rivers with a discharge bias smaller than 10% are
rather found in the BoB or in South-East Asia, such as the
Irawaddy in Myanmar, Subarnarekha, Brahmani, and Kali
in India, Mahaweli in Sri Lanka, Perak, Terengganu, and Pa-
hang in Malaysia (see Figs. 4 and A4).

These significant differences in mean discharge values
cannot be attributed to errors in river catchment identifica-
tion or boundary resolution. The total continental area of
the 336 river catchments connected to the domain is indeed
estimated to be equivalent, 13.247 and 13.168 million km²
in GloFASv3.1 (0.1° res.) and in HydroSHEDS (0.5 km
res.), respectively (i.e. a 0.6% difference only). The correla-
tion of river catchment areas between GloFASv3.1 and Hy-
droSHEDS is > 0.99, with an NRMSE of < 0.01% and an
average bias of 327 km² among the 47 common rivers of
Fig. 4. NRMSE and biases of river catchment areas compared
to D&T are higher, since this climatology is based on mea-
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Illustration of the emission and dispersion of the buoyant particle plumes in June 2017 in the three studied regions, ArS, BoB and MoC, in the IndOc.HR simulation. The initial positions are shown as pink dots and color shades show the particle count in 2 km grid cells, highlighting the 30 d trajectories. A regional view of particle release locations is shown in Fig. 1.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Daily discharges from GloFASv3.1 (model forcing, in black) and GloFASv4 (blue) for five major rivers in the Northern Hemisphere in Asia (Ganges-Brahmaputra-Meghna Delta, Irrawaddy, Indus, Salouen, Godavari) and in the Southern Hemisphere in East Africa (Zambezi, Tana, Rufiji, Betsiboka, Incomati) over the year 2017 (in m³ s⁻¹). The corresponding mean annual values calculated from GloFASv3.1, v4, and HydroSHEDS (green) are displayed with dashed lines and the D&T monthly climatology is in yellow.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Comparison of annual mean discharges (left) and corresponding upstream catchment areas (right) for 47 rivers recorded in three datasets: the assimilated GloFASv3.1 model (GF, 0.1°), HydroSHEDS GIS (HyS, 15 arcsec, 1/240°) and Dai & Trenberth climatology (D&T). Pearson correlation coefficients (R), NRMSE, and bias between GF and HyS and between GF and D&T are calculated from the 47 river mean values and shown in the headers.</figcaption></figure>

measurement stations sometimes far upstream of river mouths (NRMSE = 11% and bias = 32 500 km²). This also explains why GloFASv3.1 amplitudes in monthly mean discharges match poorly with D&T for some rivers like the Chatt-el-Arab, Manamblo or Betsiboka (Fig. A4).

This analysis thus reveals that the GloFASv3.1 version, used to prescribe river forcing in our simulations, tends to overestimate discharges, especially during peak flow periods (Fig. 3). The seasonality is however well represented. One of the explanations may be that GloFAS is forced by ERA5 atmospheric product, which tends to overestimate precipitation in the tropical zone (land and sea; Cucchi et al., 2020; Hasler and Lauer, 2021) with a tendency to show wet biases. Moreover, regional studies have observed systematic overestimation of river discharges, such as for the Ganges and Brahmaputra rivers (Hossain et al., 2023). This suggests that the structure of the hydrological model LISFLOOD on which the GloFAS-v3.1 product is based could also participate to the bias, by overestimating runoff from precipitation or underestimating evaporation from land and vegetation, which may vary between regions, depending on the disparity of the available observations used to calibrate the hydrological model (Alfieri et al., 2020). Those significant regional differences highlight the need for further investigations of hydrological processes in specific regions, such as East Africa

or South Asia. Despite these biases, using GloFAS is an improvement in terms of river discharge representation, as it allows many more rivers to be included at a higher temporal resolution than traditional climatology. However, due to the overestimation observed in v3.1, our subsequent work has adopted the more recent GloFASv4, which will be used for upcoming multi-annual simulations to support future realistic Lagrangian scenarios based on river sources.

## 3 Evaluation of the simulations and circulation analysis

In this section, we first assess in Sect. 3.1 the ability of the IndOc.HR reference simulation to reproduce the spatial and temporal variability of ocean dynamics, including large and mesoscale hydrodynamics, through an analysis of the temperature, salinity and SLA. We then examine the difference related to the spatial resolution with IndOc.12 and the wave-induced effect through Stokes-Coriolis force and Stokes drift with IndOc.HR-Sto in Sect. 3.2. For model-observation comparisons through bias and correlation calculations, the SYMPHONIE model outputs are interpolated onto the corresponding observational product grids, i.e. SST to OSTIA (0.05°), SSS to CCI (0.5°), and SLA to CMEMS (0.25°).
The year 2017 used for comparison is characterized by relatively weak positive Indian Ocean Dipole anomaly (Zhang et al., 2018), no significant ENSO pattern (transition year from El Niño to La Niña, Kersalé et al., 2022), and few cyclonic systems during the cyclone season in the Southwest Indian Ocean, with only one notable system being ENAWO hitting northern Madagascar in March. In the following, the description of the circulation and the naming of currents is based on Schott and McCreary (2001) and Schott et al. (2009) with a schematic representation on Fig. A1.

## 3.1 The IndOc.HR reference simulation

### 3.1.1 Seasonal surface spatial patterns

The mean simulated and observed SST, SSS, and SLA are mapped (Fig. 5, IndOc.HR simulation), along with their corresponding biases, for the boreal winter (December, January, February - DJF) and summer (June, July, August - JJA). These two periods correspond partly to the Northeastern (NE) and Southwestern (SW) monsoon seasons, respectively.

At the basin scale, the simulated SST from IndOc.HR exhibits a significant spatial correlation with OSTIA observations ($R = 0.99$ in DJF and JJA, Fig. 5 column 1). The simulation produces a cold SST bias on average, with a negative mean difference of $-0.17$ °C in DJF and $-0.27$ °C in JJA. Cold anomalies are mainly located to the east and north of Madagascar, encompassing the Northeast and Southeast Madagascar Current and the East African Coastal Current (see SST bias maps in Fig. 5-1c, -1f). During winter (DJF), the modeled SST fits perfectly with observations around the equator. This suggests an accurate representation of the cyclonic circulation between the Indonesian island of Sumatra to the east and the Somali Current to the west, that follows the Northeast Monsoon Current connecting the BoB to the ArS and the South Equatorial Counter-Current. Similarly, SST representation is good in the MoC and at the center of the BoB and ArS. On the other hand, warm (positive) biases are observed south of the MoC and north of the BoB (Fig. 5-1c). During summer (JJA), negative biases are more widespread, notably in the MoC, in the Andaman Sea ($-0.5$ to $-1$ °C) and north of the ArS. The South Java Current is also colder in summer in the simulation than in OSTIA. On the other hand, the Southern Gyre, the Great Whirl and associated upwelling wedges off the Somali coast tend to be warmer in the simulation, with a positive SST bias of around $0.4$ °C. The mesoscale instabilities in the southern frontal zone, corresponding to the Retroreflection of the Agulhas Current (around $40$° S) and within the Superygyre (below $20$° S), are simulated with temperature biases reaching $-2$ to $2$ °C, attributable to the chaotic nature of ocean eddies, which cause strong spatial variability in their location.

The simulated SSS from IndOc.HR also shows a strong spatial correlation with the CCI observations for both sea-

sons ($R = 0.93$ in DJF and $0.91$ in JJA, Fig. 5, column 2). There is a slight negative bias of $-0.035$ psu in DJF and $-0.032$ psu in JJA. Although the $0.5°$ resolution of the CCI product may not be sufficient to capture fine changes in coastal salinity, the model shows a tendency for river plumes to be fresher than the satellite observations. As suggested in the Sect. 2.4, these low salinity biases could be explained by the tendency of ERA5 to overestimate precipitation in the tropical zone (Cucchi et al., 2020; Hassler and Lauer, 2021, land and sea), which affects our simulations on two levels, firstly through the atmospheric forcing and secondly through the river forcing, as GloFAS is also forced by ERA5, overestimating river discharges. A second hypothesis is that the hydrological model used in GloFAS itself (LISFLOOD) could contribute to the biases by overestimating runoff from precipitation or underestimating evaporation from land and vegetation, which may vary between regions, depending on the disparity of the available observations used for calibration (Alfieri et al., 2020). This disparity is confirmed by the regional variability of the salinity bias. During both the NE monsoon (DJF) and the SW monsoon (JJA), the Gulf of Thailand and the Java Sea are fresher in the model compared to the CCI product ($-0.5$ to $-2$ psu). This is likely due to the influence of the Mekong freshwater plume to the north and to the flow of the Barito River to the south, the second main river of the Indonesian island of Borneo. Although the water in the Java Sea is fresher, this characteristic is not reflected in the Indian Ocean south of Java (with anomalies of $0$ to $+0.5$ psu), strongly supplied by the Indonesian Throughflow (Zhang et al., 2022). Fresher biases are also observed off the Indus estuary NE of the ArS (especially in summer during the SW monsoon with $-1$ to $-2$ psu) and off the Zambezi River in the MoC (particularly during the austral summer (DJF) with the rainy season, i.e. $-0.5$ to $-1$ psu). On the contrary, the modeled BoB SSS is found to be saltier than the observations, with biases of up to $+2$ psu in the north, where the SSS is typically low around $31-32$ psu. These significant differences could be linked to a poor representation of the circulation of the freshwater released by the Ganges-Brahmaputra-Meghna (GBM) delta in the extreme north of the BoB.

The spatial correlation between the simulated and observed SLA is lower than for SST and SSS, but remains statistically significant (at more than 99% with $p$-values < 0.01), with $R = 0.58$ in DJF and $0.47$ in JJA (Fig. 5, column 3). No significant mean bias is observed ($+0.001$ to $+0.005$ m in both seasons). These lower SLA correlation values are mainly due to discrepancies in the mesoscale representation, particularly in the Southern Hemisphere, such as along the chaotic frontal zone of the Agulhas Current Retroreflection, in the MoC and along the eddies of the Southeast Madagascar Current.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Spatial distribution of average SST (°C, column 1), SSS (psu, column 2) and SLA (m, column 3) over the boreal winter (DJF) and summer (JJA) of 2017: (a, d) in the IndOc.HR reference simulation and (b, e) in satellite observations, followed by (c, f) the average seasonal bias (model minus observations). The spatial Pearson correlation coefficient (R) and the average bias (B) are displayed over the maps.</figcaption></figure>

### 3.1.2 Regional surface annual cycles

Figure 6 shows the modeled and observed annual cycle of daily SST, SSS, and SLA in each regional subdomain (BoB, ArS and MoC). Table 2 evaluate the scores of the three simulations compared with the satellite observations, through the calculation of the correlation coefficient, the bias and the NRMSE.

The correlation of IndOc.HR with OSTIA SST exceeds 0.99 in all regions, with NRMSE below 12%. SST annual biases remain moderate, from -0.34 °C in BoB to -0.46 °C in MoC and -0.42 °C in ArS (Table 2). The fit is also consistent with in situ RAMA buoy data (Fig. A2), where similar cold biases are observed (ranging from -0.71 to -0.15 °C). These negative anomalies are known and reported in other
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Time series of daily SST (°C, first column), SSS (psu, second column), SLA (m, third column) anomalies over 2017, spatially averaged over the three sub-regions (ArS, BoB and Moc, see Fig. 1). The observations are in black and the three modeled configuration in blue (IndOc.HR), red (IndOc.HR-Sto) and yellow (IndOc.12).</figcaption></figure>

regional studies (Trinh et al., 2024). It is related to cold biases in ECMWF-ERA5 atmospheric forcing, estimated by Yang et al. (2021) at -0.2 °C on average compared with OS-TIA over a 15-year trend over the Indian Ocean (based on Fig. 10d of the reference study).

For SSS, the performance of IndOc.HR varies between basins. The temporal correlation with the CCI product reaches 0.98 in BoB and 0.95 in MoC, but is lower in ArS ($R = 0.68$, Table 2). The highest NRMSE is found in the MoC (73%), reflecting the strong variability in salinity along the African coasts, especially in the vicinity of the Zambezi River. Biases are generally negative compared to observations, with fresh biases down to -0.30 psu in BoB, -0.18 psu in ArS, and -0.09 psu in MoC. These fresh anomalies (down to -0.5 psu, Fig. 6) occur in SON for the ArS and BoB, and in MAM for the MoC, corresponding to periods following strong freshwater river discharges (starting in June in the Northern Hemisphere, and in December in the Southern Hemisphere, see Fig. 3). Spatially, the low salinity biases were mostly localized in the coastal zone near the mouths of the main rivers (i.e. Ganges-Brahmaputra-Meghna plume in the BoB and Indus in the ArS in JJA, Zambezi plume in the MoC in DJF, Fig. A3). It should be noted that CCI satellite data at 0.5° may not reliably represent the salinity of the coastal zone, it has indeed been shown that SSS from SMOS and SMAP satellites composing the CCI product have a good

consistency in representing seasonal SSS variations, but tend to underestimate these variations substantially (Fournier and Lee, 2021), as it is the case in comparison with our simulation. Conversely, during the drier seasons, saltier anomalies (up to +0.3 psu) tend to appear in the model, notably in MAM for the ArS and BoB, and in SON for the MoC. Positive anomalies mainly occur offshore like in the BoB (Fig. 5).

The SLA temporal correlation for IndOc.HR is strongest in ArS ($R = 0.93$), moderate in BoB ($R = 0.62$), and low in MoC ($R = 0.34$). The associated NRMSE values are 13%, 25%, and 35% respectively (Table 2). SLA biases remain weak overall (ranging from -1.2 mm in MoC to +1.1 mm in BoB), although the relatively lower values in MoC may be due to the dominance of mesoscale and submesoscale dynamics in this region with high intrinsic variability and no seasonality observed in SLA variation (Fig. 6).

### 3.1.3 Regional vertical profiles

The IndOc.HR simulation accurately reproduces the properties of the water masses, as demonstrated by the good fit between the simulated vertical temperature and salinity profiles and the GLORYS assimilated products (Fig. 7). The stratification varies significantly between sub-regions.

In the ArS, the vertical profiles confirm the modeled slight cold (-0.36 °C) and fresh (-0.085 psu) bias in the surface layer (down to 70 m depth). Between 100 and 250 m depth, a
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Temperature (°C) and Salinity (psu) vertical profiles averaged over the three sub-regions (ArS, BoB and MoC shown in Fig. 1) over 2017. All daily and mean profiles and biases (Eq. A4) spatially averaged, between the model outputs (in blue for IndOc.HR, red for IndOc.HR-Sto and yellow for IndOc.12) and the observations (in black) are represented. The Mean Absolute Error (MAE, Eq. A5) and the Standard Deviation (SD, Eq. A6) of the biases are shown.</figcaption></figure>

warm bias appears between 0.5 and 1 °C, likely related to an overly deep thermocline. Despite these biases, the mean absolute error (MAE) and the standard deviation (SD) in temperature remain low (0.057 and 0.13 °C, respectively). The salinity of the intermediate layer is better represented and the errors are minimal (MAE = 0.007 psu and SD = 0.017 psu).

In the BoB, the temperature profiles are similar, typical of tropical zones, but slightly colder at depth than in the ArS (with an average of 6.5 °C at 1000 m in the BoB instead of 8 °C in the ArS): the surface and intermediate waters of the ArS are denser than those of the BoB. The IndOc.HR simulation reproduces well the temperature and salinity characteristics in the surface layer (down to 30 m), although it has a slight cold (−0.26 °C) and salty bias (0.076 psu). Then, at the thermocline depths (between 30 and 150 m), the IndOc.HR simulation goes from a −0.4 cold to a 0.4 °C warm bias. The shallower halocline (between 20 and 100 m depth) goes from a salty to a fresh bias around 55 m depth. The MAE in temperature and salinity are equal to 0.063 °C and 0.006 psu, respectively,

and the SD equal to 0.098 °C and 0.013 psu, respectively.

In the MoC, IndOc.HR exhibits weaker vertical gradients with less important river discharges (Fig. 1) than in the BoB. On average, surface water temperatures are colder (26.4 °C) than in the northern zones (28.0 °C in the ArS and 28.6 °C in the BoB), with a less pronounced thermocline. The IndOc.HR temperature profiles are quite close to the observations in the surface layer. Contrary to other regions, the temperature biases at depth (up to more than 1000 m) remain significant with ±0.2 °C (MAE = 0.15 °C and SD = 0.18 °C). The surface and intermediate waters of the IndOc.HR show a fresh bias compared to the observations (MAE = 0.011 psu and SD = 0.023 psu).

## 3.2 Sensitivity to spatial resolution and wave forcing

We now evaluate how the two sensitivity experiments, IndOc.HR-Sto (including wave forcing) and IndOc.12 (lower
spatial resolution), differ from the reference simulation IndOc.HR.

### 3.2.1 Seasonal surface spatial patterns

The SST spatial correlation remains equally high in all three simulations ($R = 0.99$ for both DJF and JJA), indicating robust large-scale pattern agreement with OSTIA in all configurations. The cold bias is slightly enhanced in the sensitivity simulations, with mean values of $-0.22$ to $-0.35$ °C (IndOc.HR-Sto) and $-0.20$ to $-0.37$ °C (IndOc.12), compared to $-0.17$ to $-0.27$ °C in IndOc.HR (Table 1, column 1). This indicates that both the inclusion of wave forcing and the coarser resolution tend to increase the negative SST bias, especially in JJA.

The SSS patterns are well captured in all simulations, with correlation coefficients ranging from $R = 0.91$ to $0.93$. Seasonal biases are slightly reduced in the sensitivity experiments, with negative mean biases of $-0.027$ to $-0.016$ psu in IndOc.HR-Sto and $-0.006$ to $+0.001$ psu in IndOc.12 (Table 1, column 2). These results indicate that both sensitivity configurations tend to simulate slightly saltier surface waters than IndOc.HR, especially in summer. The general patterns and freshwater plume locations remain consistent.

The SLA correlations are slightly lower than those for SST and SSS, due to the complexity of mesoscale dynamics and intrinsic variability, as well as the higher sensitivity of SLA to the location and structure of mesoscale features. The SLA discrepancy might be further amplified by the coarser resolution of the satellite-derived SLA product (0.25°) compared to the model resolution (2–9 km), which can lead to higher decorrelation when comparing fine scale dynamics. SST and SSS are generally less affected by such discrepancies due to their smoother spatial gradients. However, all correlations remain statistically significant ($p < 0.01$). The spatial correlation coefficients vary very little between the three simulations, ranging from $R = 0.47$ in IndOc.HR-Sto in JJA to $R = 0.58$ in IndOc.HR-Sto and IndOc.12 in DJF (Table 1, column 3). The mean SLA biases remain small in all simulations between $+0.001$ m in JJA and $+0.005$ m in DJF. Both sensitivity tests don't change the SLA correlations much, as the observations have a coarser resolution (0.25°) than our three simulations, and in particular don't capture the small scales represented by IndOc.HR and IndOc.HR-Sto.

### 3.2.2 Regional surface annual cycles

Compared to IndOc.HR, the IndOc.HR-Sto simulation shows similar performance in SST, with equally high correlation coefficients ($R = 0.99$) and slightly increased cold biases: $-0.49$ °C in ArS, $-0.39$ °C in BoB, and $-0.47$ °C in MoC. NRMSE values remain close (10%–13%). For SSS, IndOc.HR-Sto shows a slight improvement in correlation in ArS ($R = 0.77$ vs. 0.68) and BoB ($R = 0.98$ unchanged), while maintaining a strong score in MoC ($R = 0.95$). The

fresh biases are slightly reduced in MoC ($-0.08$ psu vs. $-0.09$ psu), but unchanged in ArS ($-0.18$ psu). NRMSE decreases slightly in all regions (e.g., from 73% to 71% in MoC). SLA results improve slightly with the inclusion of wave forcing in ArS ($R = 0.95$ vs. $0.93$), and in MoC ($R = 0.41$ vs. $0.34$), and remain similar in BoB ($R = 0.60$). The NRMSE also decreases slightly in MoC (28% vs. 35%). This suggests a small positive effect of wave forcing on SLA realism, particularly in regions dominated by mesoscale activity. This might be attributed to the Stokes-Coriolis force, which may modify the surface dynamics and current shear (McWilliams and Restrepo, 1999).

The IndOc.12 simulation shows slightly more different results. SST correlation remains high ($R \ge 0.98$ in all regions), but biases are slightly stronger, especially in MoC ($-0.62$ °C) and BoB ($-0.43$ °C). NRMSE also increases in MoC (14%), while remaining comparable elsewhere. For SSS, IndOc.12 shows lower NRMSE and bias than the HR configurations, particularly in BoB (NRMSE = 6%, bias = $-0.04$ psu) and MoC (NRMSE = 33%, bias < 0.01 psu). Temporal correlations also remain high ($R > 0.94$ everywhere). This result, although counterintuitive, may reflect a better consistency between the coarser CCI resolution (0.5°) and the smoother representation of coastal salinity gradients in our low resolution IndOc.12 simulation, rather than an improved physical realism. The SLA temporal fit in IndOc.12 decreases in BoB ($R = 0.36$ vs. 0.60–0.62 in HR runs), and is close to zero in MoC ($R = -0.09$). The SLA NRMSE is also the highest in all regions (33% in BoB, 43% in MoC), indicating that the coarser resolution degrades the representation of SLA variability, particularly in eddy-rich and dynamically complex zones.

### 3.2.3 Regional vertical profiles

In the ArS, both IndOc.12 and IndOc.HR-Sto simulate temperature profiles in line with IndOc.HR, with MAE values in the same range (0.054 °C) and low bias variability ($SD = 0.12$ to 0.14 °C). However, IndOc.12 shows a stronger cold bias ($-0.47$ °C) in the surface layer compared to HR ($-0.36$ °C). Salinity errors, already low in IndOc.HR (MAE = 0.007 psu), are slightly reduced in IndOc.12 (0.005 psu) with SD going from 0.017 in IndOc.HR to 0.015 in IndOc.HR-Sto and 0.012 psu in IndOc.12. The performance in all three simulations remains very close in this region.

In the BoB, where vertical stratification is more pronounced, more differences appear in temperature and salinity profiles. In the surface layer (down to 30 m), the cold bias is enhanced in IndOc.12 ($-0.42$ °C vs. $-0.26$ °C in HR) as well as the the salty bias ($+0.24$ psu vs. $+0.076$ psu in HR). Then, the thermocline bias transitions (between 30 and 150 m) observed in IndOc.HR (from cold to warm) is very attenuated in IndOc.HR-Sto, while IndOc.12 shows only a
<table><caption>**Table 1.** Spatial Pearson correlation coefficient (Eq. A2) and seasonal spatial bias (Eq. A4) over the whole domain between the three annual simulations and the observations for SST, SSS and SLA. Associated *p*-values are all < 0.01.</caption><thead><tr><th rowspan="3"></th><th colspan="4">SST - OSTIA</th><th colspan="4">SSS - CCI</th><th colspan="4">SLA - ALTI</th></tr><tr><th colspan="2">R</th><th colspan="2">Bias</th><th colspan="2">R</th><th colspan="2">Bias</th><th colspan="2">R</th><th colspan="2">Bias</th></tr><tr><th>DJF</th><th>JJA</th><th>DJF</th><th>JJA</th><th>DJF</th><th>JJA</th><th>DJF</th><th>JJA</th><th>DJF</th><th>JJA</th><th>DJF</th><th>JJA</th></tr></thead><tbody><tr><td>IndOc.HR</td><td>0.99</td><td>0.99</td><td>-0.17</td><td>-0.27</td><td>0.93</td><td>0.91</td><td>-0.035</td><td>-0.032</td><td>0.56</td><td>0.50</td><td>0.004</td><td>0.001</td></tr><tr><td>IndOc.HR-Sto</td><td>0.99</td><td>0.99</td><td>-0.22</td><td>-0.35</td><td>0.93</td><td>0.91</td><td>-0.027</td><td>-0.016</td><td>0.58</td><td>0.47</td><td>0.005</td><td>0.001</td></tr><tr><td>IndOc.12</td><td>0.99</td><td>0.99</td><td>-0.20</td><td>-0.37</td><td>0.93</td><td>0.92</td><td>-0.006</td><td>0.001</td><td>0.58</td><td>0.52</td><td>0.001</td><td>0.001</td></tr></tbody></table>

<table><caption>**Table 2.** Temporal Pearson correlation coefficient (Eq. A2), bias (Eq. A4) and NRMSE (Eq. A3 in %) between the three annual simulations and the observation datasets for SST, SSS and SLA averaged over the three sub-regions. All *p*-values are < 0.01.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="3">SST - OSTIA</th><th colspan="3">SSS - CCI</th><th colspan="3">SLA - ALTI</th></tr><tr><th>R</th><th>Bias (°C)</th><th>NRMSE</th><th>R</th><th>Bias (psu)</th><th>NRMSE</th><th>R</th><th>Bias (mm)</th><th>NRMSE</th></tr></thead><tbody><tr><td rowspan="3">ArS</td><td>IndOc.HR</td><td>0.99</td><td>-0.42</td><td>12</td><td>0.68</td><td>-0.18</td><td>41</td><td>0.93</td><td>-0.4</td><td>13</td></tr><tr><td>IndOc.HR-Sto</td><td>0.99</td><td>-0.49</td><td>13</td><td>0.77</td><td>-0.18</td><td>39</td><td>0.95</td><td>0.4</td><td>10</td></tr><tr><td>IndOc.12</td><td>0.99</td><td>-0.45</td><td>12</td><td>0.82</td><td>-0.09</td><td>25</td><td>0.93</td><td>0.7</td><td>11</td></tr><tr><td rowspan="3">BoB</td><td>IndOc.HR</td><td>0.99</td><td>-0.34</td><td>10</td><td>0.98</td><td>-0.30</td><td>19</td><td>0.62</td><td>1.1</td><td>25</td></tr><tr><td>IndOc.HR-Sto</td><td>0.99</td><td>-0.39</td><td>12</td><td>0.98</td><td>-0.22</td><td>15</td><td>0.60</td><td>0.4</td><td>25</td></tr><tr><td>IndOc.12</td><td>0.98</td><td>-0.43</td><td>13</td><td>0.99</td><td>-0.04</td><td>6</td><td>0.36</td><td>1.4</td><td>33</td></tr><tr><td rowspan="3">MoC</td><td>IndOc.HR</td><td>&gt; 0.99</td><td>-0.46</td><td>11</td><td>0.95</td><td>-0.09</td><td>73</td><td>0.34</td><td>-1.2</td><td>35</td></tr><tr><td>IndOc.HR-Sto</td><td>&gt; 0.99</td><td>-0.47</td><td>11</td><td>0.95</td><td>-0.08</td><td>71</td><td>0.41</td><td>-1.0</td><td>28</td></tr><tr><td>IndOc.12</td><td>&gt; 0.99</td><td>-0.62</td><td>14</td><td>0.94</td><td>&lt; 0.01</td><td>33</td><td>-0.09</td><td>-0.5</td><td>43</td></tr></tbody></table>

cold bias of -1 °C. Temperature MAE is slightly better in IndOc.HR-Sto (0.058 °C), and variability increases in IndOc.12 (SD = 0.15 °C), compared to IndOc.HR (0.098 °C) and HR-Sto (0.082 °C). For salinity, IndOc.12 is also slightly worse (MAE = 0.008 psu and SD = 0.025 psu) compared to IndOc.HR (MAE = 0.006 psu; SD = 0.013 psu).

In the MoC, where vertical gradients are weaker, the differences between IndOc.12 and both HR simulations are accentuated. IndOc.12 shows a -1.4 °C cold bias (at 55 m depth) while both HR simulations are quite close to the observations. The MAE in temperature are similar for the three simulation, while the SD in IndOc.12 is higher (0.23 °C vs. 0.18–0.19 °C in HR). For salinity, all simulations show small MAE (0.011–0.013 psu), but IndOc.12 again exhibits more variable biases (SD = 0.025 psu), compared to 0.022–0.023 psu in HR.

### 3.2.4 Synthesis of simulations performance

The comparison of spatial patterns and regional annual cycles reveals that no clear hierarchy can be established regarding the quality of the three simulations performed: comparable correlation coefficients and error values (bias, NRMSE) for SSS and SST, in space and time (Tables 1 and 2) for IndOc.HR, HR-Sto and IndOc.12. However, there are differences when considering the SLA and vertical profiles, reveal-

ing significant better performances for IndOc.HR and HR-Sto.

The high-resolution configurations (IndOc.HR and IndOc.HR-Sto) perform better in representing SLA variability, particularly in mesoscale-active regions such as BoB and MoC. This is likely due to their finer representation of ocean dynamics, which is not captured by the lower resolution IndOc.12. However, SLA differences remain low due to the relatively coarse resolution of the compared observations (0.25°). Interestingly, IndOc.12 shows similar or even lower SSS errors in some regions, possibly benefiting from a resolution closer to that of the CCI product. However, IndOc.12 exhibits stronger biases and variability in vertical temperature and salinity profiles, particularly in BoB and MoC, reflecting a less accurate representation of stratification and subsurface structure. All three configurations reproduce key observed features, but HR simulations offer the most balanced performance, combining a good agreement with surface observations and vertical structures.

### 3.3 Seasonal variability of surface tracers and inferred dynamics

This section interprets the seasonal variability of SST and SSS, and the associated regional circulation patterns, which
aligns with key dynamical features described by Schott et al. (2009), using the same current naming (see Fig. A1).

In the Northern Hemisphere, the SST peak in the ArS and BoB occurs in mid-May, with observed regional averages of 30.2 and 30.8 °C respectively (Fig. 6). Temperatures then decrease due to the summer monsoon-driven circulation reversal: cold waters from the Southern Hemisphere flow up into the ArS following the South Equatorial Current, the Northeast Madagascar Current and the East African Current. Significant upwelling wedges along the coasts of Somalia and Oman, associated with major eddies, i.e. the Southern Gyre, Great Whirl and Socotra Gyre, accentuate the surface cooling (well visible in the ArS in May–October, Fig. A3). The Southwest Monsoon Current connects the cooled waters from the ArS to the southern BoB and the Andaman Sea. During this season, strong precipitation and freshwater discharges also contribute to the cooling of surface waters in the BoB (Fig. 3), together with the strong vertical mixing caused by Southwest monsoon winds. The SST in those northern regions reaches an observed average summer minimum around 27.5 to 28.5 °C in August.

Starting in September, the second inter-monsoon season is characterized by the second circulation reversal to the south of India and a strong decrease in the intensity of the East African Current, as well as the Somalia and Oman upwellings (see the maps of mean kinetic energy (MKE, Eq. A7) in DJF in Fig. A5). As a result, the SST increases by 1.5 °C in the ArS and 0.5 °C in the BoB until November (Fig. 6). This period marks the sea surface winter cooling by the strong Northeast monsoon winds blowing to the North of the sub-basins (as shown in Figs. 5 and A3 during DJF and November–April). The regional SST reaches its minimum in February, around 26.3° in the ArS and 27.1° in the BoB.

In the Southern Hemisphere, the annual cycle of SST is reversed, and only two seasons are observed (unlike the four seasons in the North, Fig. 6). The hot and rainy season in the MoC (austral summer in November–April, Fig. A3) causes SST peaks in February–March, with observed mean regional SST at 29 °C. The cold and dry season (austral winter in May–October, Fig. A3) brings an abrupt drop in SST, with a mean regional minimum of around 24.5 °C in August.

In the Northern Hemisphere, the annual cycle of the SSS differs greatly between sub-basins: the mean SSS in the ArS is 3.7 psu higher than in the BoB (35.9 and 32.2 psu respectively) due to much lower freshwater inflows (Fig. 1) and higher evaporation. Moreover, the SSS variability in BoB has a higher amplitude. The summer monsoon, which starts in June (when the regional mean SSS is at its maximum of 32.9 psu), brings heavy rainfall to the continent, resulting in a rapid increase in river discharges to the BoB. The maximum discharge of around 160 000 m³ s⁻¹ in the Ganges-Brahmaputra-Meghna delta in Bangladesh occurs in mid-August (Fig. 3) and reduces the SSS by almost 2 psu in the sub-basin, with a minimum of 31 psu in mid-October. The decrease in river discharges to the north from September

(Fig. 3), combined with the inversion of circulation from the SW summer monsoon to the NE winter monsoon, leads to a further increase in SSS in the BoB from December. This is due to the significant export of cold, low-salty water to the ArS in the south of India (Fig. A3), which induces its mean SSS drop of almost 0.3 psu (Fig. 6).

To the South, the annual variability of SSS in the MoC is low (0.4 psu) around an annual regional average of 35.1 psu. The SSS starts to decrease in December, which marks the beginning of the austral rainy season, reaching its minimum in April. This season is characterized by high river discharges in East Africa and Madagascar, the major river in the area being the Zambezi, with a GloFASv3.1 maximum discharge of around 50 000 m³ s⁻¹ (see Fig. 3). The austral dry winter (starting in May), is then characterized by a significant discharge reduction and thus a SSS increase in the coastal zone (Fig. 6).

# 4 Energy budgets analysis

The energy associated with ocean currents plays a critical role in the transport of marine debris, influencing their dispersion velocity, spread direction and resulting distribution. A detailed analysis of energy levels in circulation simulations is thus essential to understand and discuss the effect of simulated current fields on Lagrangian trajectories. Although Lagrangian kinetic energy (KE) computed along particle trajectories, would provide a more direct estimate of the energy experienced by marine debris, we focus here on Eulerian KE derived from the model gridded velocity field. This allows for a consistent comparison between simulations. In this section, we compare the energy budgets between our three simulations: IndOc.HR, .HR-Sto and .12.

## 4.1 Kinetic Energy budget

Figure 8 depicts the daily temporal evolution over 2017 (in addition to the two months of spin-up) of the KE averaged over the same IndOc domain for the three simulations conducted in this study (IndOc.HR, .HR-Sto and .12), for the surface (Fig. 8a) and over the whole depth (Fig. 8b). The energy in the 3D domain stabilizes over a two-months spin-up period (November to December 2016). The three simulations exhibit similar temporal variability at different time scales (daily to seasonal), but with varying levels of energy. The IndOc.HR-Sto configuration with wave forcing produces the most energetic simulation at the surface with an averaged KE of 62.88 m⁻² s⁻² (Fig. 8a), 8% higher than IndOc.HR (Fig. 8c). In contrast, the IndOc.12 simulation has a significantly lower surface energy level than the two others with an averaged KE of 48.22 m⁻² s⁻² (Fig. 8a), 17% smaller than in IndOc.HR, due to different scale and energy cascade representation (Fig. 8c). Although the differences are attenuated when integrating vertically, resulting
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Time series of (a) the mean surface kinetic energy (KE in $m^{-2}s^{-2}$), (b) the mean 3D KE over the whole depth, and (c) the ratio of the mean KE of IndOc.HR-Sto (red) and IndOc.12 (yellow) simulations over IndOc.HR (blue) for the surface (solid lines) and the 3D (dashed lines). The three simulations are averaged over the same entire IndOc domain. The mean values over 2017 and the two-month spin-up are shown in the colored boxes (corresponding to the dotted lines).</figcaption></figure>

in negligible differences on average over the depth and period (Fig. 8b, 6.59, 6.58 and $6.64 \times 10^{-2}$ for IndOc.HR, IndOc.12 and IndOc.HR-Sto respectively), the low-resolution IndOc.12 simulation produces less energy during the most energetic periods (August to December), as expected.

## 4.2 Energy spectral density

The observed difference in energy amplitude between the HR simulations (1–3 km over the NW half domain) and the 1/12° (≈9 km) simulation can be attributed to the representation of mesoscale (50–300 km) and submesoscale (<50 km) activity (Herrmann et al., 2008; McWilliams, 2016; Waldman et al., 2017). To investigate dynamics spatial and temporal scales and eddy intensity across each region and grid resolutions, we examine their wavenumber-frequency power spectral density of KE as shown in Figs. 9, A6 and A7. For comparison purposes, the IndOc.HR, IndOc.HR-Sto and IndOc.12 daily average surface outputs were interpolated onto a 1/50° (~2 km) regular grid prior to spectral analysis. The

wavenumber-frequency spectra (Figs. A6 and A7) are cut at the frequency of 2 d due to the daily time sampling of the outputs. Calculations are performed over the DJF and JJA seasons for the three studied regions (ArS, BoB, MoC). To complement, surface vorticity snapshots (Fig. 10) illustrate the distinct spatial dynamics resolved by the IndOc.HR and IndOc.12 simulations at meso and submesoscales.

The most energetic region is MoC, followed by BoB and then ArS (an order of magnitude separates the energy levels of MoC and ArS). This is supported by the vorticity snapshots in Fig. 10, where the MoC shows the largest and most intense mesoscale coherent eddies whatever the season. In general, the KE density is systematically higher in both HR simulations than in IndOc.12, for all regions and in every season, except for BoB during DJF where KE levels associated with large scale motions are higher in the IndOc.12 simulation (Fig. 9). Moreover, in every region and every season, there is a change in the KE spectral slopes as we move towards the small scales (below 50 km for ArS and MoC, below 40 km for BoB): the energy dissipates faster in the IndOc.12 simulation than in the IndOc.HR with more turbulence, and the energy gap between both resolutions increases. Vorticity snapshots on Fig. 10 show that sub-mesoscale structures such as small eddies and filaments are much more numerous and better represented in IndOc.HR than IndOc.12. The KE density spectra for IndOc.HR and IndOc.HR-Sto are quite similar in the three regions, but the IndOc.HR-Sto develops more energy than IndOc.HR during the JJA season in the ArS (at all scales) and during DJF in the MoC (at large scales), and less than IndOc.HR during JJA in the MoC (at large scales).

The disparities in terms of wavenumber-frequency spectra (Figs. A6 and A7) are highlighted in the high-frequency (<10 d) and high-wavenumber (<50 km) values that correspond to short-lived submesoscale balanced motions, but also in mesoscale motions (>50 km). In the majority of 2D spectra, the energy peak is centered on structures with wavelengths of 100 to 200 km and periods of 10 to 30 d. These spectra show that a better representation of fine scales not only results in an increase in energy in the fine scales, but that it also leads to an inverse cascade of energy from the small to the large scales, also raising the energy levels of the large-scale structures with short and long periods (Gula et al., 2016; Klein et al., 2019). The exception of the BoB region, where the IndOc.12 simulation shows energy levels higher than IndOc.HR in winter (DJF), could be explained by the energy of the atmospheric forcings not dissipating towards the underdeveloped fine scales and therefore being concentrated in the large scales.

## 4.3 Scale representation and seasonality

For all regions, wavenumber spectra in Fig. 9 show that the energy differences between IndOc.12 and both HR simulations are greatest in winter (corresponding to DJF for ArS
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Wavenumber spectra of surface KE (m² s⁻²) for the IndOc.HR (blue), IndOc.12 (yellow) and Indoc.HR-Sto (red) simulations, over the DJF (solid lines) and JJA (dashed lines) seasons, for the three studied regions (ArS, BoB and MoC in Fig. 1).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 10.** Surface vorticity snapshots (s⁻¹) for IndOc.HR and IndOc.12 simulations on 1 January and 1 August 2017 over the three sub-regions (ArS, BoB and MoC, see blue boxes in Fig. 1).</figcaption></figure>
and BoB and JJA for MoC) and particularly in the motions scales below 100 km (through dissipation).

Figure 9 shows that the ArS exhibits the greatest seasonal differences in all simulations, with the summer season (JJA), corresponding to the southwest monsoon, generally more energetic than the winter season (DJF). In the western part of the sub-basin, baroclinic and barotropic instabilities are enhanced during this JJA period (Zhan et al., 2020), forming climatological mesoscale eddies such as the Great Whirl and the Socotra Eddy in the Somali Current region (see Fig. A1 and MKE, EKE maps on Fig. A5). Eddies are advected from the western formation zone towards the eastern part of the ArS (Trott et al., 2018), as shown by the EKE maps (Fig. A5). This west-to-east transport of energy is enhanced in the IndOc.HR simulation compared to IndOc.12. The JJA season in the ArS is also marked by maximum wind stress in the region (Fischer et al., 2002; Varna et al., 2023), leading to a stronger Stokes drift (Fig. A8), which in turn enhances the KE levels in IndOc.HR-Sto at this period.

In the BoB, seasonality is more visible for scale motions smaller than 100 km (Fig. 9), with more mesoscale and sub-mesoscale structures developed in winter (DJF), as shown on the vorticity snapshots (Fig. 10). However, we can see from the EKE (Fig. A5) that although the energy linked to eddies and sub-mesoscale is more homogeneous throughout the BoB in winter, it is in summer that the Southwest Monsoon Current, which bypasses Sri Lanka and runs up the west coast of the Gulf, is the strongest, associated with more intense mesoscale eddies. In the Northern Hemisphere, the temporal variability of the eddy field is thus partly related to instabilities resulting from changes in monsoon regimes and intraseasonal equatorial zonal winds (Cheng et al., 2018).

In the MoC, the differences between the seasons are the smallest, with large scales (> 50 km) slightly more energetic during DJF (warm rainy and cyclone season), whereas fine scales are slightly more energetic during JJA (cold dry season) in the IndOc.HR simulation. The MoC is part of the larger Agulhas Current system and is characterized by intense mesoscale eddy activity, dominated by large anti-cyclones (~300 km) propagating southward, and occurring regularly over the year, as reported by Halo et al. (2014). These mesoscale structures are accurately reproduced by both IndOc.HR and IndOc.12 simulations, as illustrated by the surface vorticity snapshots (Fig. 10). However, only the IndOc.HR simulation represents the submesoscale structures and filaments along the coast and between the large eddies, which is more energetic during JJA (the austral winter).

## 4.4 Synthesis of energy levels and scale dynamics

The energy analysis highlights that the high-resolution simulation (IndOc.HR) more accurately captures both mesoscale and submesoscale dynamics compared to the lower resolution of IndOc.12 simulation, as expected. IndOc.HR exhibits higher KE levels at nearly all scales and in all regions,

particularly in the 10–100 km range, where submesoscale processes are active. These fine-scale structures are associated with enhanced vorticity and filaments, especially along boundary currents and in energetic regions like the MoC. The IndOc.12 simulation dissipates energy more rapidly at small scales, failing to represent the full spectrum of eddy activity. This limitation leads to an underestimation of KE levels in both mesoscale and submesoscale ranges, with exceptions such as the BoB during DJF, where forcing energy accumulates at larger scales due to insufficient dissipation. The inclusion of the wave forcing in IndOc.HR-Sto through the Stokes-Coriolis force and Stokes drift leads to an increase in surface KE, particularly in regions and seasons of high wind stress (ArS during JJA, MoC during DJF), highlighting the role of wave-current interactions in energizing upper-ocean dynamics. Spectral results also indicate signs of an inverse cascade, an upscale transfer of energy from submesoscale to larger mesoscale structures, especially visible in the HR configurations. This finding aligns with previous theoretical and numerical work showing that submesoscale turbulence can energize larger scales through nonlinear interactions (Gula et al., 2016; Klein et al., 2019).

These differences emphasize the importance of resolving fine-scale dynamics for realistic ocean current simulations, with implications for tracer dispersion, vertical transport, and Lagrangian trajectory predictions.

# 5 Lagrangian trajectories

## 5.1 Trajectory statistics

The IndOc simulations are designed with the main objective of modeling the dispersion of marine debris in the Indian Ocean. In this final section, we present the effect of the three simulated current fields (IndOc.HR, IndOc.HR-Sto and IndOc.12) evaluated above on the representation of Lagrangian trajectories in the ocean surface layer. We simulate trajectories in our three regions with contrasting dynamics (ArS, BoB, MoC), and we take into account the temporal variability of the circulation through monthly releases of particles during a year (experiment described in Sect. 2.2.2 and illustrated in Fig. 2). Figure 11 shows the trajectory statistics obtained for each of the three simulations (IndOc.HR, IndOc.HR-Sto and IndOc.12), considering the plumes from the three regions (ArS, BoB, MoC) combined. Figure 12 shows these same trajectory statistics for each region separately. Finally, Fig. 13 shows trajectory roses with the main trajectory directions and distances traveled by the particles throughout the year for each region.

Figure 11a shows that the PDF of trajectory lengths reaches a maximum 100–150 km smaller for particles released in IndOc.12 than in both HR simulations. The mean trajectory length in IndOc.12 (737 km) is 168 km shorter than in IndOc.HR (905 km) and 157 km shorter than in IndOc.HR-
<figure><img src="image_13.png" /><figcaption>**Figure 11.** Statistical analysis of the Lagrangian trajectories obtained from the three ocean current simulations, i.e. IndOc.HR (blue), IndOc.HR-Sto (red) and IndOc.12 (yellow), considering the plumes from the three studied regions. (a) Histogram of monthly trajectory lengths, i.e. actual distances traveled. (b) Time evolution of the mean displacement of the plumes, i.e. the mean straight-line distance from the release position during the 30 d drift of each particle over 2017. (c) Monthly variability of the mean monthly trajectory length and (d) median monthly trajectory length with [min, max] distances traveled by the particles, and, (e) mean monthly displacement between the release position and the position after 30 d of drift (mean straight-line distance).</figcaption></figure>

Sto (894 km). Similarly, the median trajectory length obtained from IndOc.12 is 122 km shorter than in IndOc.HR and 143 km shorter than in IndOc.HR-Sto. The spread of the trajectory length distribution is greater in the case of the HR simulations compared with IndOc.12, as shown by the standard deviation values (455 km in IndOc.12 vs. 578 and 537 km in IndOc.HR and IndOc.HR-Sto respectively). Figure 11b shows that the mean displacement of the particle plume is accentuated not only by the HR current fields but also by the Stokes drift. The particles end their 30 d trajectories on average 6% closer to the release position in the IndOc.12 simulation than in IndOc.HR, and 8% further away in the IndOc.HR-Sto simulation. The deviation between both HR simulations and IndOc.12 is visible from the second day

of drift. The deviation between IndOc.HR and IndOc.HR-Sto is visible from the ninth day of drift. These deviations in trajectory lengths and plume displacements reflect the differences in KE levels generated at the ocean surface by the three different simulations (Fig. 8c): the higher KE in IndOc.HR and IndOc.HR-Sto enables stronger and more variable currents that enhance particle dispersion, while the lower KE in IndOc.12 results in slower and less dispersed trajectories. Looking at the monthly variability over the year (Fig. 11c, d, e), the statistical differences between the three simulations can vary over time, but the trajectories produced by both HR simulations are on average always longer (except in February but still very close) and more spread out (except in January,
<figure><img src="image_14.png" /><figcaption>**Figure 12.** Statistical analysis of the Lagrangian trajectories for the three studied regions separately, i.e. ArS (column 1), BoB (column 2) and MoC (column 3), obtained from the three ocean current simulations of the study, i.e. IndoC.HR (blue), IndOc.HR-Sto (red) and IndOc.12 (yellow), (a) Histogram of the monthly trajectory lengths of the particles. (b) Monthly evolution of the mean ± SD, and (c) median trajectory lengths with [min, max] values, during 30 d over the year 2017.</figcaption></figure>

October and December but still very close) than those produced by the IndOc.12 simulation.

## 5.2 Regional analysis

In order to evaluate the variability of these statistics, it is preferable to analyze the three regions separately in order to link the trajectories generated to the regional dynamic processes at play. Whatever the region, the average trajectory length of the particles in 30 d is always lower in IndOc.12, and the maximum of the PDF is reached on average for shorter trajectory lengths (Fig. 12). The distances traveled by the particles are much longer and the distribution of trajectories is wider in the MoC, followed by trajectories in the BoB and then in the ArS. This is consistent with the KE ratios reported in the spectra of Fig. 9, which show a MoC region that is much more energetic than BoB and ArS, with the presence of coherent mesoscale eddies of high intensity (Fig. 10). The most pronounced difference between the trajectories produced by both HR simulations and

IndOc.12 also corresponds to the MoC region (Fig. 12a). This suggests that in highly dynamic regions with a strong geostrophic component, the resolution of the current fields is a critical factor for the study of Lagrangian dispersion in terms of spread velocities and trajectory lengths.

In the ArS region, the trajectories from HR simulations are clearly longer than in IndOc.12 from March to August, then tend to be equivalent on average from September to February (Fig. 12b, c). The particle dispersion is accentuated in the IndOc.HR-Sto simulation that exhibits the most significant seasonal variability of trajectory statistics (which is consistent with the seasonal variation in KE spectra of Fig. 9, where JJA exhibits increased energy levels particularly at meso and submesoscales, enhancing particle dispersion during the monsoon season): mean and median trajectory lengths are almost 300 km greater in IndOc.HR-Sto current fields from March to August compared to IndOc.12 and from May to August compared to IndOc.HR. This annual variability corresponds to the monsoon cycle (observed for
<figure><img src="image_15.png" /><figcaption>**Figure 13.** Annual trajectory roses for the three regions (ArS, BoB, MoC) for the three simulations (IndOc.HR, IndOC.HR-Sto, IndOC.12). The size of the arrows indicates the frequency of particles that follow the corresponding direction (%). The colors indicate the length of the particle 30 d trajectories in the corresponding direction.</figcaption></figure>

example in SLA curves in Fig. 6). From March to August (positive SLA in the ArS), the reversal of the circulation in the Northern Hemisphere is followed by the south-west summer monsoon (JJA). This is a period of maximum wind stress in the region (Fischer et al., 2002; Varna et al., 2023), producing the highest Stokes drift, as shown in Fig. A8, with a higher KE particularly in large and meso-scales and a more intense vorticity field than the rest of the year, as illustrated in Fig. 9 (JJA curves) and Fig. 10. This demonstrates the link between dispersion velocity and submesoscale representation associated with higher energy levels. It is interesting to note that during the south-west monsoon season, the prevailing direction of Lagrangian trajectories in the ArS is East-Southeast, following the total current direction (Fig. A9, JJA), while the Stokes drift points Northeast-East (Fig. A8, JJA). The trajectory roses demonstrate that the wave forcing can play a key role in shaping the buoyant particle drift,

centuating significantly the particle velocities and thus the trajectory lengths during strong wind periods. Throughout the year, particles move most frequently in an East-Southeast direction (Fig. 13, ArS), demonstrating the dominance of the southwest monsoon over the northeast monsoon in dispersion patterns.

In the BoB region, the closer energy levels between the three simulations (Fig. 9) is reflected by closer trajectory statistics (Fig. 12a). While KE associated with large to meso scales (> 60 km) in the IndOc.12 simulation is higher than both HR simulations in DJF, there is no clear differences in the average trajectory length. However, the fact that HR simulations still have more energy associated with fine scales (< 40 km) is sufficient to produce trajectory statistics on average longer than in IndOc.12 (Fig. 12a). The effect of the monsoonal circulation on trajectory statistics is not as marked as in the ArS (Fig. 12b, c). This is associated with
much greater variability in particle trajectories in terms of direction and length (as shown by Fig. 13, BoB). No seasonality is indeed observed in energy levels associated with large scale (> 100 km) in the HR simulations, while accentuated energy levels are associated with small scales (< 100 km) during DJF compared to JJA (Fig. 9, BoB). During the summer monsoon (JJA), the prevailing direction of particle trajectories is Northeast-East (Fig. A10). The presence of Stokes drift in the IndOc.HR-Sto simulation does not accentuate the trajectory lengths, nor does it accentuate the energy densities to any great extent (Fig. 9, BoB). It is interesting to note that during this JJA period, there is more variability in trajectory directions in the IndOc.12 than in both HR simulations. This probably reflects the lack of coherent mesoscale structures in the lower resolution IndOc.12, resulting in more dispersed particle directions. In contrast, the HR simulations better resolve mesoscale features that guide particles along more consistent pathways, reducing directional variability.

In the MoC, a slight decrease in the trajectory lengths is visible from May to October, corresponding to the cold and dry season. The KE levels associated with large to mesoscales (> 40 km) are indeed lower in the JJA period than in the DJF period (Figs. 9, A6, A7). However, the dissipation through submesoscale (< 40 km) is enhanced during JJA compared with DJF. The difference in statistics between IndOc.HR-Sto and IndOc.HR is smaller in the MoC compared to the other two regions, indicating that the effect of waves is less significant here, both on the distances traveled (Fig. 12) and on the direction of the trajectories (Fig. A11). The Stokes drift in this region is southeast-tonorth on average throughout the year with lower intensities than in the Northern Hemisphere (Fig. A8), contrasting with the predominant energy transfers of the oceanic circulation from northeast to southwest (Fig. A5). As a result, particles move most frequently in an East-Southeast direction (Fig. 13, MoC). Generally, in the IndOc.HR-Sto simulation, the particles tend to be channeled more in privileged directions, whereas in the IndOc.12 simulation the variability of the trajectories is accentuated in all directions (as in the other two regions, Fig. 13). The wave forcing effect may be negligible compared with the strong geostrophic component of ocean dynamics in this region.

## 5.3 Synthesis of resolution and wave forcing effects on particle dispersion

The IndOc.HR simulation, used as reference, produces longer and more dispersed Lagrangian trajectories than the lower resolution of IndOc.12, highlighting the crucial role of high-resolution current fields in capturing fine-scale structures that enhance dispersion. The differences in trajectory lengths and plume displacement are visible from the early stages of drift and are consistent with higher kinetic energy levels at the surface in IndOc.HR. Compared to IndOc.HR, the IndOc.HR-Sto simulation, which includes the wave ef-

fect through the Stokes-Coriolis force and Stokes drift, generates even more extended and anisotropic trajectories, especially during the monsoon season in the ArS, where wave-induced transport becomes significant. However, in regions dominated by geostrophic dynamics, such as the MoC, the Stokes drift contribution remains marginal. IndOc.HR offers a better representation of mesoscale and submesoscale variability than IndOc.12, while the implementation of wave forcing in IndOc.HR-Sto further improves Lagrangian transport under strong wind periods. These results confirm that both resolution and wave effects are key factors to accurately simulate surface drift in the Indian Ocean.

## 6 Conclusions

The goal of this study is to address current limitations in modeling the Indian Ocean circulation and its implications for marine debris dispersion in the region, by developing a new regional configuration with enhanced realism. Existing large-scale simulations of the Indian Ocean generally rely on horizontal resolutions of 1/4 to 1/12°, and rarely combine the key forcings that influence Lagrangian dispersion, such as tides, wind, waves, and freshwater inputs. A more accurate representation of debris transfer, transport and accumulation in the ocean was therefore crucial, requesting to develop new appropriate modeling approaches. In this study, we integrate those key forcings and address the sub-mesoscale dynamics with a new telescopic grid, allowing both a kilometric resolution over the areas of interest and a coverage of the whole basin circulation. We compare our annual reference simulation (IndOc.HR) with two annual sensitivity simulations to respectively assess the impact of resolution (IndOc.12 at 1/12° ≈ 9 km) and wave effect (IndOc.HR-Sto including a parameterization of the Stokes-Coriolis force and the explicit Stokes drift) on ocean dynamics and marine debris dispersion.

Our reference simulation IndOc.HR successfully reproduces the spatial and temporal variability of sea surface (temperature, salinity, height) and vertical (temperature, salinity) properties and the main characteristics of ocean dynamics, particularly in the three contrasting regions studied: the Arabian Sea (ArS), the Bay of Bengal (BoB) and the Mozambique Channel (MoC). Although a systematic cold SST bias (particularly in the western subtropical region north and east of Madagascar) and low-salinity bias in coastal regions are observed, these biases are consistent with previous studies that used ERA5 or GloFASv3.1 forcings. The reference simulation realistically reproduces the regional circulation over the Indian Ocean, in good agreement with observations and literature, in particular the seasonal reversal of the monsoon-driven circulation in the Northern Hemisphere and the intense eddy activity characteristic of the Southern Hemisphere. The HR and 1/12° simulations show similar agreement with observations in terms of basin-scale surface cir-
culation, making it difficult to clearly establish a hierarchy between the two at that scale, however the HR simulation shows notable improvements in the representation of vertical profiles of temperature and salinity.

The high-resolution simulations provide higher energy budget and a more detailed representation of fine-scale dynamics, especially at the sea surface, as compared to IndOc.12. The high resolution logically resolve smaller eddies and filaments, increasing the energy associated with fine scales but also meso- to large scales, depending on the region and the season. The differences between the high-resolution and 1/12° simulations are largest in winter (corresponding to DJF for ArS and BoB, and JJA for MoC). The wave forcing has a more significant effect on energy budgets during summer periods (southwest monsoon in ArS and rainy season in MoC), increasing the energy levels especially associated with large to meso- scales (> 50 km). The most energetic region is the MoC, with the largest and most intense mesoscale coherent eddies and enhanced vertical mixing, followed by the BoB and last by the ArS, which is more stratified. The ArS exhibits the greatest seasonal variations in all simulations.

In terms of Lagrangian statistics, the three simulations exhibit interesting differences. Both HR simulations consistently generate longer 30 d trajectory lengths compared to the 1/12° simulation (about 160 km longer on average), and larger mean displacement of the plumes (+6% to +14% with wave effect). The distances traveled by the particles are much longer and the distribution of trajectories is wider in the MoC, followed by trajectories in the BoB and then in the ArS, consistent with the energy ratios. Regardless of the region, the resolution of the current fields is a critical factor for the study of Lagrangian dispersion in terms of spread velocities and trajectory lengths. This is accentuated in highly dynamic regions with a strong geostrophic component, like in the MoC. Seasonal variability is enhanced in more stratified regions with large-scale seasonal circulation (such as ArS). The inclusion of wave forcing in IndOc.HR-Sto further amplifies this effect, particularly during strong wind periods, aligning the Stokes drift with the monsoon-driven currents. The Stokes drift increases the energy levels at meso- to large scales, which is most pronounced during the southwest monsoon and rainy season, and significantly influences the direction of particle dispersion. Stokes drift has a key role in accentuating particle velocities and travel distances, especially when it reinforces the main surface circulation directions, leading to more coherent directional dispersion and less chaotic spreading compared to the lower-resolution simulation.

In conclusion, our evaluation and energy budget analysis highlights and quantifies the added value of high spatial resolution and wave forcing accounting for numerical studies on buoyant marine debris dispersion, emphasizing their role in the representation of the wide range of scales involved in buoyant particle trajectories. These results provide essential insights for future dispersion scenarios from

the coastal zone, where the main sources of marine debris are located and where plastic stranding poses significant socio-environmental challenges, to the open basin. While the present release scenarios are idealized, the implications for realistic debris dispersion will only be fully assessed once coastal and riverine sources will be included. Our regional configuration and outputs will be useful for future studies on pollution spreading scenarios but also marine connectivity studies in the Indian Ocean, and are now freely available to the interested scientific community.

## Appendix A: Formulas for the described variables and additional figures

Sea Level Anomaly calculated from the Sea Surface Height:

$$SLA(t) = SSH(t) - \overline{SSH} \quad (A1)$$

Pearson correlation coefficient:

$$R = \frac{\sum_{i=1}^{n} (x_{\text{mod},i} - \bar{x}_{\text{mod}})(x_{\text{obs},i} - \bar{x}_{\text{obs}})}{\sqrt{\sum_{i=1}^{n} (x_{\text{mod},i} - \bar{x}_{\text{mod}})^2 \sum_{i=1}^{n} (x_{\text{obs},i} - \bar{x}_{\text{obs}})^2}} \quad (A2)$$

Normalized Root Mean Square Error:

$$NRMSE = \frac{\sqrt{\frac{1}{n} \sum_{i=1}^{n} (x_{\text{mod},i} - x_{\text{obs},i})^2}}{x_{\text{obs,max}} - x_{\text{obs,min}}} \quad (A3)$$

Bias between model and observations:

$$\text{Bias} = \frac{1}{n} \sum_{i=1}^{n} (x_{\text{mod},i} - x_{\text{obs},i}) \quad (A4)$$

Mean Absolute Error of the bias vertically weighted:

$$\text{MAE} = \frac{\sum_{i=1}^{n} |\text{Bias}_i| \Delta z_i}{\sum_{i=1}^{n} \Delta z_i} \quad (A5)$$

Standard Deviation of the bias vertically weighted:

$$SD = \sqrt{\frac{\sum_{i=1}^{n} (\text{Bias}_i)^2 \Delta z_i}{\sum_{i=1}^{n} \Delta z_i}} \quad (A6)$$

with:

- n is the number of temporal or spatial samples used in the calculation.

- $x_{\text{mod},i}$, $x_{\text{obs},i}$ are the modeled and observed parameter values at index i in time or space. $\bar{x}_{\text{mod}}$, $\bar{x}_{\text{obs}}$ denote temporal means over the 2017 period or over the whole grid for spatial correlation.

- $x_{\text{obs,max}}$, $x_{\text{obs,min}}$ are the maximum and minimum observed parameter values over the 2017 period.

- $\Delta z_i$ is the thickness of the i'th vertical layer [m].
Mean Kinetic Energy:

$$MKE = \frac{1}{2} (\bar{u}^2 + \bar{v}^2) \qquad (A7)$$

Eddy Kinetic Energy:

$$EKE = \frac{1}{2} (u - \bar{u})^2 + (v - \bar{v})^2 \qquad (A8)$$

Current vorticity:

$$\zeta = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} \qquad (A9)$$

with:

- $u, v$ are the zonal and meridional components of the total current velocity [m s⁻¹].

- $\bar{u}, \bar{v}$ are the mean velocities averaged over the specified period (DJF or JJA).

<figure><img src="image_16.png" /><figcaption>**Figure A1.** Spatial distribution of mean annual MiP specific fluxes (number of particles 10⁶ yr⁻¹ km⁻²) discharged by model river catchments into the IndOc domain. The estimations are based on the empirical model of Weiss et al. (2021) and HydroSHEDS river catchment delineation (Linke et al., 2019). The schematic near-surface circulation during the southwest (July–August in orange) and northeast monsoon (in green) is taken from Schott et al. (2009). Black arrows indicate annual mean flows with no seasonal reversals in the Southern Hemisphere. Current branches indicate the South Equatorial Current (SEC), South Equatorial Counter-Current (SECC), Northeast and Southeast Madagascar Current (NEMC and SEMC), East African Coastal Current (EACC), Southern Gyre (SG) and Great Whirl (GW) and associated upwelling wedges during the southwest monsoon only (pink shades), Southwest and Northeast Monsoon Currents (SMC and NMC), Indonesian Through-flow (ITF).</figcaption></figure>
<figure><img src="image_17.png" /></figure>

**Figure A2.** Daily SST and SSS time series for the year 2017 comparing the IndOc.HR, IndOc.HR-Sto and IndOc.12 simulations with in situ measurements from six RAMA buoys. The time bias for each simulation is shown above the figures. Buoy coordinates are plotted in Fig. 1.
<figure><img src="image_18.png" /><figcaption>**Figure A3.** Spatial distribution of the averaged SST (°C, columns 1–2), SSS (psu, columns 3–4), SLA (m, columns 5–6) during the Winter Monsoon in the North, i.e. warm rainy season in the South (November–April) and Summer Monsoon in the North, i.e. cold dry season in the South (May–October), in model outputs and satellite observations for the three sub-regions: the Arabian Sea (ArS), the Bay of Bengal (BoB) and the Mozambique Channel (MoC).</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Figure A4.** Time series of climatological monthly mean discharge from GloFASv3.1 forcing (over 2016–2021) and from Dai & Trenberth climatology (scattered data between 1900 and 2018 by river).</figcaption></figure>
<figure><img src="image_20.png" /></figure>

**Figure A5.** Seasonal spatial distribution of the MKE (Mean Kinetic Energy, Eq. A7) and EKE (Eddy Kinetic Energy, Eq. A8) from the IndOc.HR and IndOc.12 simulations, during the boreal winter (DJF) and summer (JJA) of 2017.
<figure><img src="image_21.png" /><figcaption>**Figure A6.** Frequency-Wavenumber spectra of KE (m² s⁻²) for the high (IndOc.HR, 1–3 km, blue) and low (IndOc.12, 9 km, yellow) resolution simulations over the DJF season for the three studied regions (ArS, BoB and MoC, smaller blue boxes in Fig. 1). MBM stands for mesoscale balanced, SBM for submesoscale balanced.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>**Figure A7.** Frequency-Wavenumber spectra of KE (m² s⁻²) for the high (IndOc.HR, 1–3 km, blue) and low (IndOc.12, 9 km, yellow) resolution simulations over the JJA season for the three studied regions (ArS, BoB and MoC, smaller blue boxes in Fig. 1). MBM stands for mesoscale balanced, SBM for submesoscale balanced.</figcaption></figure>
<figure><img src="image_23.png" /><figcaption>**Figure A8.** Monthly average Stokes drift (m s⁻¹) from the global CMEMS WAVERYS 0.2° product (CMEMS, 2019) used to force the IndOc.HR-Sto simulation.</figcaption></figure>
ArS

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

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

**Figure A9.** Monthly trajectory roses in the ArS region for the three simulations (IndOc.HR, IndOC.HR-Sto, IndOC.12). The size of the arrows indicates the frequency of particles that follow the corresponding direction. The colors indicate the length of the 30 d trajectories in the corresponding direction (pink: < 300, orange: 300–600, green: 600–900, blue: > 900 km).
BoB

<figure><img src="image_26.png" /><figcaption>**Figure A10.** Monthly trajectory roses in the BoB region for the three simulations (IndOc.HR, IndOC.HR-Sto, IndOC.12). The size of the arrows indicates the frequency of particles that follow the corresponding direction. The colors indicate the length of the particle 30 d trajectories in the corresponding direction (pink: < 400, orange: 400–800, green: 800–1200, blue: > 1200 km).</figcaption></figure>
MoC

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

**Figure A11.** Monthly trajectory roses in the MoC region for the three simulations (IndOc.HR, IndOC.HR-Sto, IndOC.12). The size of the arrows indicates the frequency of particles that follow the corresponding direction. The colors indicate the length of the particle 30 d trajectories in the corresponding direction (pink: < 700, orange: 700–1400, green: 1400–2100, blue: > 2100 km).
*Code and data availability.* The current version 3.6.6 of the SYMPHONIE source code (hydrodynamics and Lagrangian modules), together with both Indian Ocean configuration grid files used to produce the results presented in this paper, are archived under the following DOI: https://doi.org/10.6084/m9 figshare.28903067 (Weiss, 2025). Last version of the model can also be obtained from the INSU-certified service SIROCCO at https://sirocco obs.mip.fr (last access: 23 January 2026). Considering their size, raw three-dimensional simulation outputs are not stored publicly but are available upon request to the corresponding author. Forcing data are described in the section along with references and upload addresses. New simulations based on this configuration can be produced for any specified years of interest and for any ocean parameters to be studied.

*Author contributions.* LW, MH and CM contributed to the study conceptualization and methodological strategy. LW developed the model configuration, performed the hydrodynamic and Lagrangian simulations as well as data analysis and figures. PM supported the model configuration development and numerical parametrization, in particular the wave forcing. MB contributed to evaluate the simulation performance compared with observation databases. All authors contributed to the interpretation of the results. LW wrote the original draft and all authors reviewed and edited the manuscript. All authors have read and agreed to the final 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 made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.

*Acknowledgements.* This project was provided with computing HPC and storage resources by GENCI at TGCC thanks to the grants A0120110098, A0140110098 and A0160110098 on the supercomputer Joliot Curie's SKL partition, and by CALMIP on projects P20055 and P1325. We would like to kindly thank Nadia Ayoub for leading the scientific project submitted to the TGCC and sharing the computing time required to produce the simulations studied, as well as Thomas Duhaut for providing tutorials on the SYMPHONIE grid tool called Bathy Maker. The authors would like to thank the IRD OCEANS department for the postdoctoral fellowship of Lisa Weiss. Finally, we thank Bruno Sainte-Rose and the anonymous referee for their valuable comments and helpful suggestions that improved the manuscript.

*Financial support.* This research was supported by the postdoctoral fellowship of Lisa Weiss funded by the French National Research Institute for Sustainable Development (IRD, 2022–2023).

*Review statement.* This paper was edited by Vassilios Vervatis and reviewed by Bruno Sainte-Rose and one anonymous referee.

## References

Alfieri, L., Lorini, V., Hirpa, F. A., Harrigan, S., Zsoter, E., Prudhomme, C., and Salamon, P.: A global streamflow re-analysis for 1980–2018, Journal of Hydrology X, 6, 100049, https://doi.org/10.1016/j.hydroa.2019.100049, 2020.

Bajon, R., Huck, T., Grima, N., Maes, C., Blanke, B., Richon, C., and Couvelard, X.: Influence of waves on the three-dimensional distribution of plastic in the ocean, Marine Pollution Bulletin, 187, 114533, https://doi.org/10.1016/j.marpolbul.2022.114533, 2023.

Bentsen, M., Evensen, G., Drange, H., and Jenkins, A. D.: Coordinate Transformation on a Sphere Using Conformal Mapping, Monthly Weather Review, 127, 2733–2740, https://doi.org/10.1175/1520-0493(1999)127<2733:ctoasu>2.0.co;2, 1999.

Boutin, J., Vergely, J.-L., Reul, N., Catany, R., Koehler, J., Martin, A., Rouffi, F., Arias, M., Chakroun, M., Corato, G., Estella-Perez, V., Guimbard, S., Hasson, A., Josey, S., Khvorostyanov, D., Kolodziejczyk, N., Mignot, J., Olivier, L., Reverdin, G., Stammer, D., Supply, A., Thouvenin-Masson, C., Turiel, A., Vialard, J., Cipollini, P., and Donlon, C.: ESA Sea Surface Salinity Climate Change Initiative (Sea_Surface_Salinity_cci): weekly and monthly sea surface salinity products, v03.21, for 2010 to 2020, NERC EDS Centre for Environmental Data Analysis, https://doi.org/10.5285/5920A2C77E3C45339477ACD31CE62C3C, 2021.

Chassignet, E. P., Xu, X., and Zavala-Romero, O.: Tracking Marine Litter With a Global Ocean Model: Where Does It Go? Where Does It Come From?, Frontiers in Marine Science, 8, https://doi.org/10.3389/fmars.2021.667591, 2021.

Cheng, X., McCreary, J. P., Qiu, B., Qi, Y., Du, Y., and Chen, X.: Dynamics of Eddy Generation in the Central Bay of Bengal, J. Geophys. Res.-Oceans, 123, 6861–6875, https://doi.org/10.1029/2018jc014100, 2018.

Choi, J., Bracco, A., Barkan, R., Shchepetkin, A. F., McWilliams, J. C., and Molemaker, J. M.: Submesoscale Dynamics in the Northern Gulf of Mexico. Part III: Lagrangian Implications, Journal of Physical Oceanography, 47, 2361–2376, https://doi.org/10.1175/jpo-d-17-0036.1, 2017.

Couvelard, X., Lemarié, F., Samson, G., Redelsperger, J.-L., Ardhuin, F., Benshila, R., and Madec, G.: Development of a two-way-coupled ocean-wave model: assessment on a global NEMO(v3.6)–WW3(v6.02) coupled configuration, Geosci. Model Dev., 13, 3067–3090, https://doi.org/10.5194/gmd-13-3067-2020, 2020.

Cucchi, M., Weedon, G. P., Amici, A., Bellouin, N., Lange, S., Müller Schmied, H., Hersbach, H., and Buontempo, C.: WFDE5: bias-adjusted ERA5 reanalysis data for impact studies, Earth Syst. Sci. Data, 12, 2097–2120, https://doi.org/10.5194/essd-12-2097-2020, 2020.

Cummings, J. A. and Smedstad, O. M.: Variational Data Assimilation for the Global Ocean, Springer Berlin Heidelberg, 303–343, ISBN 9783642350887, https://doi.org/10.1007/978-3-642-35088-7_13, 2013.
Cunningham, H. J., Higgins, C., and van den Bremer, T. S.: The Role of the Unsteady Surface Wave-Driven Ekman–Stokes Flow in the Accumulation of Floating Marine Litter, J. Geophys. Res.-Oceans, 127, https://doi.org/10.1029/2021jc018106, 2022.

Dai, A.: Dai and Trenberth Global River Flow and Continental Discharge Dataset, NSF National Center for Atmospheric Research [data set], https://doi.org/10.5065/D6V69H1T, 2017.

D'Asaro, E. A., Shcherbina, A. Y., Klymak, J. M., Molemaker, J., Novelli, G., Guigand, C. M., Haza, A. C., Haus, B. K., Ryan, E. H., Jacobs, G. A., Huntley, H. S., Laxague, N. J. M., Chen, S., Judt, F., McWilliams, J. C., Barkan, R., Kirwan, A. D., Poje, A. C., and Özgökmen, T. M.: Ocean convergence and the dispersion of flotsam, Proceedings of the National Academy of Sciences, 115, 1162–1167, https://doi.org/10.1073/pnas.1718453115, 2018.

Dobler, D., Huck, T., Maes, C., Grima, N., Blanke, B., Martinez, E., and Ardhuin, F.: Large impact of Stokes drift on the fate of surface floating debris in the South Indian Basin, Marine Pollution Bulletin, 148, 202–209, https://doi.org/10.1016/j.marpolbul.2019.07.057, 2019.

Estournel, C., Marsaleix, P., and Ulses, C.: A new assessment of the circulation of Atlantic and Intermediate Waters in the Eastern Mediterranean, Progress in Oceanography, 198, 1026–1049, https://doi.org/10.1016/j.pcean.2021.102673, 2021.

Estournel, C., Mikolajczak, G., Ulses, C., Bourrin, F., Canals, M., Charmasson, S., Doxaran, D., Duhaut, T., de Madron, X. D., Marsaleix, P., Palanques, A., Puig, P., Radakovitch, O., Sanchez-Vidal, A., and Verney, R.: Sediment dynamics in the Gulf of Lion (NW Mediterranean Sea) during two autumn–winter periods with contrasting meteorological conditions, Progress in Oceanography, 210, 10294–102942, https://doi.org/10.1016/j.pcean.2022.102942, 2023.

E.U. Copernicus Marine Service Information (CMEMS): Global Ocean OSTIA Sea Surface Temperature and Sea Ice Reprocessed, Marine Data Store (MDS) [data set], https://doi.org/10.48670/MOI-00168, 2015.

E.U. Copernicus Marine Service Information (CMEMS): Global Ocean Physics Reanalysis, Marine Data Store (MDS) [data set], https://doi.org/10.48670/MOI-00021, 2018.

E.U. Copernicus Marine Service Information (CMEMS): Global Ocean Waves Reanalysis, Marine Data Store (MDS) [data set], https://doi.org/10.48670/MOI-00022, 2019.

E.U. Copernicus Marine Service Information (CMEMS): Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Reprocessed (1993–Ongoing), Marine Data Store (MDS) [data set], https://doi.org/10.48670/MOI-00148, 2021.

Fischer, A. S., Weller, R. A., Rudnick, D. L., Eriksen, C. C., Lee, C. M., Brink, K. H., Fox, C. A., and Leben, R. R.: Mesoscale eddies, coastal upwelling, and the upper-ocean heat budget in the Arabian Sea, Deep Sea Research Part II: Topical Studies in Oceanography, 49, 2231–2264, https://doi.org/10.1016/s0967-0645(02)00036-x, 2002.

Fischer, R., Lobelle, D., Kooi, M., Koelmans, A., Onink, V., Laufkötter, C., Amaral-Zettler, L., Yool, A., and van Sebille, E.: Modelling submerged biofouled microplastics and their vertical trajectories, Biogeosciences, 19, 2211–2234, https://doi.org/10.5194/bg-19-2211-2022, 2022.

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.

Fournier, S. and Lee, T.: Seasonal and Interannual Variability of Sea Surface Salinity Near Major River Mouths of the World Ocean Inferred from Gridded Satellite and In-Situ Salinity Products, Remote Sensing, 13, 728, https://doi.org/10.3390/rs13040728, 2021.

Gamoyo, M., Obura, D., and Reason, C. J.: Estimating Connectivity Through Larval Dispersal in the Western Indian Ocean, J. Geophys. Res.-Biogeo., 124, 2446–2459, https://doi.org/10.1029/2019JG005128, 2019.

Garinet, A., Herrmann, M., Marsaleix, P., and Pénicaud, J.: Spurious numerical mixing under strong tidal forcing: a case study in the south-east Asian seas using the Symphonie model (v3.1.2), Geosci. Model Dev., 17, 6967–6986, https://doi.org/10.5194/gmd-17-6967-2024, 2024.

Gasparin, F., Greiner, E., Lellouche, J. M., Legalloudec, O., Garric, G., Drillet, Y., Bourdallé-Badie, R., Traon, P. Y. L., Rémy, E., and Drévillon, M.: A large-scale view of oceanic variability from 2007 to 2015 in the global high resolution monitoring and forecasting system at Mercator Océan, Journal of Marine Systems, 187, 260–276, https://doi.org/10.1016/j.jmarsys.2018.06.015, 2018.

Guizien, K., Belharet, M., Marsaleix, P., and Guarini, J. M.: Using larval dispersal simulations for marine protected area design: Application to the Gulf of Lions (north-west Mediterranean), Limnol. Oceanogr., 57, 1099–1112, https://doi.org/10.4319/lo.2012.57.4.1099, 2012.

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

Habib, J., Ulses, C., Estournel, C., Fakhri, M., Marsaleix, P., Pujo-Pay, M., Fourrier, M., Coppola, L., Mignot, A., Mortier, L., and Conan, P.: Seasonal and interannual variability of the pelagic ecosystem and of the organic carbon budget in the Rhodes Gyre (eastern Mediterranean): influence of winter mixing, Biogeosciences, 20, 3203–3228, https://doi.org/10.5194/bg-20-3203-2023, 2023.

Halo, I., Backeberg, B., Penven, P., Ansorge, I., Reason, C., and Ullgren, J.: Eddy properties in the Mozambique Channel: A comparison between observations and two numerical ocean circulation models, Deep Sea Research Part II: Topical Studies in Oceanography, 100, 38–53, https://doi.org/10.1016/j.dsr2.2013.10.015, 2014.

Hassler, B. and Lauer, A.: Comparison of Reanalysis and Observational Precipitation Datasets Including ERA5 and WFDE5, Atmosphere, 12, 1462, https://doi.org/10.3390/atmos12111462, 2021.

Herrmann, M., Somot, S., Sevault, F., Estournel, C., and Déqué, M.: Modeling the deep convection in the northwestern Mediterranean Sea using an eddy-permitting and an eddy-resolving model: Case study of winter 1986–1987, J. Geophys. Res.-Oceans, 113, https://doi.org/10.1029/2006jc003991, 2008.

Herrmann, M., To Duy, T., and Estournel, C.: Intraseasonal variability of the South Vietnam upwelling, South China Sea: influence of atmospheric forcing and ocean intrinsic variability in the northwestern Mediterranean, Deep Sea Research Part II: Topical Studies in Oceanography, 100, 38–53, https://doi.org/10.1016/j.dsr2.2013.10.015, 2014.
ity, Ocean Sci., 19, 453–467, https://doi.org/10.5194/os-19-453-2023, 2023.

Herrmann, M., To Duy, T., and Marsaleix, P.: Mechanisms and in-terseasonal variability in the South Vietnam Upwelling, South China Sea: the role of circulation, tides, and rivers, Ocean Sci., 20, 1013–1033, https://doi.org/10.5194/os-20-1013-2024, 2024.

Hossain, S., Cloke, H. L., Ficchì, A., Gupta, H., Speight, L., Has-san, A., and Stephens, E. M.: A decision-led evaluation approach for flood forecasting system developments: An application to the Global Flood Awareness System in Bangladesh, Journal of Flood Risk Management, https://doi.org/10.1111/jfr3.12959, 2023.

Jordà, G., Bolaños, R., Espino, M., and Sánchez-Arcilla, A.: Assessment of the importance of the current-wave cou-pling in the shelf ocean forecasts, Ocean Sci., 3, 345–362, https://doi.org/10.5194/os-3-345-2007, 2007.

Kersalé, M., Volkov, D. L., Pujiana, K., and Zhang, H.: Interan-nual variability of sea level in the southern Indian Ocean: lo-cal vs. remote forcing mechanisms, Ocean Sci., 18, 193–212, https://doi.org/10.5194/os-18-193-2022, 2022.

Klein, P., Lapeyre, G., Siegelman, L., Qiu, B., Fu, L., Torres, H., Su, Z., Menemenlis, D., and Le Gentil, S.: Ocean-Scale In-teractions From Space, Earth and Space Science, 6, 795–817, https://doi.org/10.1029/2018ea000492, 2019.

Kuizenga, B., van Emmerik, T., Waldschläger, K., and Kooi, M.: Will it Float? Rising and Settling Velocities of Com-mon Macroplastic Foils, ACS ES&T Water, 2, 975–981, https://doi.org/10.1021/acsestwater.1c00467, 2022.

Large, W. and Yeager, S.: Diurnal to decadal global forcing for ocean and sea-ice models: The data sets and flux climatologies, Tech. rep., University Corporation for Atmospheric Research, https://doi.org/10.5065/D6KK98Q6, 2004.

Laverre, M., Kerhervé, P., Constant, M., Weiss, L., Charrière, B., Stetzler, M., González-Fernández, D., and Ludwig, W.: Heavy rains control the floating macroplastic inputs into the sea from coastal Mediterranean rivers: A case study on the Tét River (NW Mediterranean Sea), Science of The Total Environment, 877, 162733, https://doi.org/10.1016/j.scitotenv.2023.162733, 2023.

Lellouche, J.-M., Greiner, E., Bourdallé Badie, R., Garric, G., Melet, A., Drévillon, M., Bricaud, C., Hamon, M., Le Galloudec, O., Regnier, C., Candela, T., Testut, C.-E., Gasparin, F., Ruggiero, G., Benkiran, M., Drillet, Y., and Le Traon, P.-Y.: The Copernicus Global 1/12° Oceanic and Sea Ice GLORYS12 Reanalysis, Frontiers in Earth Science, 9, https://doi.org/10.3389/feart.2021.698876, 2021.

Linke, S., Lehner, B., Ouellet Dallaire, C., Ariwi, J., Grill, G., Anand, M., Beames, P., Burchard-Levine, V., Maxwell, S., Moidu, H., Tan, F., and Thieme, M.: Global hydro-environmental sub-basin and river reach characteristics at high spatial res-olution, Scientific Data, 6, https://doi.org/10.1038/s41597-019-0300-6, 2019.

Lodise, J., Özgökmen, T., Gonçalves, R. C., Iskandarani, M., Lund, B., Horstmann, J., Poulain, P.-M., Klymak, J., Ryan, E. H., and Guigand, C.: Investigating the Formation of Sub-mesoscale Structures along Mesoscale Fronts and Estimating Kinematic Quantities Using Lagrangian Drifters, Fluids, 5, 159, https://doi.org/10.3390/fluids5030159, 2020.

Lyard, F. H., Allain, D. J., Cancet, M., Carrère, L., and Pi-cot, N.: FES2014 global ocean tide atlas: design and perfor-mance, Ocean Sci., 17, 615–649, https://doi.org/10.5194/os-17-615-2021, 2021.

Maes, C., Grima, N., Blanke, B., Martinez, E., Paviet-Salomon, T., and Huck, T.: A Surface “Superconvergence” Pathway Connecting the South Indian Ocean to the Subtropical South Pacific Gyre, Geophysical Research Letters, 45, 1915–1922, https://doi.org/10.1002/2017gl076366, 2018.

Marsaleix, P., Auclair, F., Floor, J. W., Herrmann, M. J., Estour-nel, C., Pairaud, I., and Ulses, C.: Energy conservation issues in sigma-coordinate free-surface ocean models, Ocean Modelling, 20, 61–89, https://doi.org/10.1016/j.ocemod.2007.07.005, 2008.

Marsaleix, P., Michaud, H., and Estournel, C.: 3D phase-resolved wave modelling with a non-hydrostatic ocean circulation model, Ocean Modelling, 136, 28–50, https://doi.org/10.1016/j.ocemod.2019.02.002, 2019.

McPhaden, M. J., Meyers, G., Ando, K., Masumoto, Y., Murty, V. S. N., Ravichandran, M., Syamsudin, F., Vialard, J., Yu, L., and Yu, W.: RAMA: The Research Moored Array for African–Asian–Australian Monsoon Analysis and Prediction, Bulletin of the American Meteorological Society, 90, 459–480, https://doi.org/10.1175/2008bams2608.1, 2009.

McWilliams, J. C.: Submesoscale currents in the ocean, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 472, 20160117, https://doi.org/10.1098/rspa.2016.0117, 2016.

McWilliams, J. C. and Restrepo, J. M.: The Wave-Driven Ocean Circulation, Journal of Physical Oceanog-raphy, 29, 2523–2540, https://doi.org/10.1175/1520-0485(1999)029<2523:twdoc>2.0.co;2, 1999.

Nguyen-Duy, T., Ayoub, N. K., Marsaleix, P., Toublanc, F., De Mey-Frémaux, P., Piton, V., Herrmann, M., Duhaut, T., Tran, M. C., and Ngo-Duc, T.: Variability of the Red River Plume in the Gulf of Tonkin as Revealed by Numerical Modeling and Clustering Analysis, Frontiers in Marine Science, 8, 1–25, https://doi.org/10.3389/fmars.2021.772139, 2021.

Onink, V., Wichmann, D., Delandmeter, P., and van Sebille, E.: The Role of Ekman Currents, Geostrophy, and Stokes Drift in the Ac-cumulation of Floating Microplastic, J. Geophys. Res.-Oceans, 124, 1474–1490, https://doi.org/10.1029/2018jc014547, 2019.

Pairaud, I., Lyard, F., Auclair, F., Letellier, T., and Marsaleix, P.: Dynamics of the semi-diurnal and quarter-diurnal in-ternal tides in the Bay of Biscay. Part 1: Barotropic tides, Continental Shelf Research, 28, 1294–1315, https://doi.org/10.1016/j.csr.2008.03.004, 2008.

Pairaud, I. L., Auclair, F., Marsaleix, P., Lyard, F., and Pi-chon, A.: Dynamics of the semi-diurnal and quarter-diurnal internal tides in the Bay of Biscay. Part 2: Baro-clinic tides, Continental Shelf Research, 30, 253–269, https://doi.org/10.1016/j.csr.2009.10.008, 2010.

Pattiaratchi, C., van der Mheen, M., Schlundt, C., Narayanaswamy, B. E., Sura, A., Hajbane, S., White, R., Kumar, N., Fernan-des, M., and Wijeratne, S.: Plastics in the Indian Ocean – sources, transport, distribution, and impacts, Ocean Sci., 18, 1–28, https://doi.org/10.5194/os-18-1-2022, 2022.

Phillips, H. E., Tandon, A., Furue, R., Hood, R., Ummenhofer, C. C., Benthuysen, J. A., Menezes, V., Hu, S., Webber, B., Sanchez-Franks, A., Cherian, D., Shroyer, E., Feng, M., Wijesekera, H., Chatterjee, A., Yu, L., Hermes, J., Murtugudde, R., Tozuka, T., Su, D., Singh, A., Centurioni, L., Prakash, S., and Wiggert, J.:
Progress in understanding of Indian Ocean circulation, variability, air-sea exchange, and impacts on biogeochemistry, Ocean Sci., 17, 1677–1751, https://doi.org/10.5194/os-17-1677-2021, 2021.

Pottapinjara, V. and Joseph, S.: Evaluation of mixing schemes in the HYbrid Coordinate Ocean Model (HYCOM) in the tropical Indian Ocean, Ocean Dynamics, 72, 341–359, https://doi.org/10.1007/s10236-022-01510-2, 2022.

Qu, T., Fukumori, I., and Fine, R. A.: Spin-Up of the Southern Hemisphere Super Gyre, J. Geophys. Res.-Oceans, 124, 154–170, https://doi.org/10.1029/2018JC014391, 2019.

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, Journal of Physical Oceanography, 47, 2077–2100, https://doi.org/10.1175/JPO-D-16-0168.1, 2017.

Rio, M., Mulet, S., and Picot, N.: Beyond GOCE for the ocean circulation estimate: Synergetic use of altimetry, gravimetry, and in situ data provides new insight into geostrophic and Ekman currents, Geophysical Research Letters, 41, 8918–8925, https://doi.org/10.1002/2014gl061773, 2014.

Rocha, C. B., Chereskin, T. K., Gille, S. T., and Menemenlis, D.: Mesoscale to Submesoscale Wavenumber Spectra in Drake Passage, Journal of Physical Oceanography, 46, 601–620, https://doi.org/10.1175/jpo-d-15-0087.1, 2016.

Röhrs, J., Christensen, K. H., Hole, L. R., Broström, G., Drivdal, M., and Sundby, S.: Observation-based evaluation of surface wave effects on currents and trajectory forecasts, Ocean Dynamics, 62, 1519–1533, https://doi.org/10.1007/s10236-012-0576-y, 2012.

Rühs, S., van den Bremer, T., Clementi, E., Denes, M. C., Moulin, A., and van Sebille, E.: Non-negligible impact of Stokes drift and wave-driven Eulerian currents on simulated surface particle dispersal in the Mediterranean Sea, Ocean Sci., 21, 217–240, https://doi.org/10.5194/os-21-217-2025, 2025.

Schott, F. A. and McCreary, J. P.: The monsoon circulation of the Indian Ocean, Progress in Oceanography, 51, 1–123, https://doi.org/10.1016/s0079-6611(01)00083-0, 2001.

Schott, F. A., Xie, S., and McCreary, J. P.: Indian Ocean circulation and climate variability, Reviews of Geophysics, 47, https://doi.org/10.1029/2007rg000245, 2009.

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

Thibault, M., Weiss, L., Fernandez, R., Avargues, N., Jaquemet, S., Lebreton, L., Garnier, J., Jaeger, A., Royer, S.-J., Cartraud, A., ter Halle, A., Marsaleix, P., Chevillon, L., Tourmetzj, J., and Le Corre, M.: Barau’s petrel, Pterodroma baraui, as a bioindicator of plastic pollution in the South-West Indian Ocean: A multifaceted approach, Marine Environmental Research, 202, 106709, https://doi.org/10.1016/j.marenvres.2024.106709, 2024.

To Duy, T., Herrmann, M., Estournel, C., Marsaleix, P., Duhaut, T., Bui Hong, L., and Trinh Bich, N.: The role of wind, mesoscale dynamics, and coastal circulation in the interannual variability of the South Vietnam Upwelling, South China Sea – answers from a high-resolution ocean model, Ocean Sci., 18, 1131–1161, https://doi.org/10.5194/os-18-1131-2022, 2022.

Trinh, N. B., Herrmann, M., Ulses, C., Marsaleix, P., Duhaut, T., To Duy, T., Estournel, C., and Shearman, R. K.: New insights into the South China Sea throughflow and water budget seasonal cycle: evaluation and analysis of a high-resolution configuration of the ocean model SYMPHONIE version 2.4, Geosci. Model Dev., 17, 1831–1867, https://doi.org/10.5194/gmd-17-1831-2024, 2024.

Trott, C. B., Subrahmanyam, B., Chaigneau, A., and Delcroix, T.: Eddy Tracking in the Northwestern Indian Ocean During Southwest Monsoon Regimes, Geophysical Research Letters, 45, 6594–6603, https://doi.org/10.1029/2018gl078381, 2018.

van der Mheen, M., Pattiaratchi, C., and van Sebille, E.: Role of Indian Ocean Dynamics on Accumulation of Buoyant Debris, J. Geophys. Res.-Oceans, 124, 2571–2590, https://doi.org/10.1029/2018JC014806, 2019.

van der Mheen, M., van Sebille, E., and Pattiaratchi, C.: Beaching patterns of plastic debris along the Indian Ocean rim, Ocean Sci., 16, 1317–1336, https://doi.org/10.5194/os-16-1317-2020, 2020.

van Sebille, E., Aliani, S., Law, K. L., Maximenko, N., Alsina, J. M., Bagaev, A., Bergmann, M., Chapron, B., Chubarenko, I., Cózar, A., Delandmeter, P., Egger, M., Fox-Kemper, B., Garaba, S. P., Goddijn-Murphy, L., Hardesty, B. D., Hoffman, M. J., Isobe, A., Jongedijk, C. E., Kaandorp, M. L. A., Khatmullina, L., Koelmans, A. A., Kukulka, T., Laufkötter, C., Lebreton, L., Lobelle, D., Maes, C., Martinez-Vicente, V., Morales Maqueda, M. A., Poulain-Zarcos, M., Rodríguez, E., Ryan, P. G., Shanks, A. L., Shim, W. J., Suaria, G., Thiel, M., van den Bremer, T. S., and Wichmann, D.: The physical oceanography of the transport of floating marine debris, Environmental Research Letters, 15, 023003, https://doi.org/10.1088/1748-9326/ab6d7d, 2020.

Varna, M., Jithin, A., and Francis, P.: Characteristics and dynamics of mesoscale eddies in the eastern Arabian sea, Deep Sea Research Part II: Topical Studies in Oceanography, 207, 105218, https://doi.org/10.1016/j.dsr2.2022.105218, 2023.

Vic, C., Hascoët, S., Gula, J., Huck, T., and Maes, C.: Oceanic Mesoscale Cyclones Cluster Surface Lagrangian Material, Geophysical Research Letters, 49, https://doi.org/10.1029/2021gl097488, 2022.

Vogt-Vincent, N. S. and Johnson, H. L.: Multidecadal and climatological surface current simulations for the southwestern Indian Ocean at 1/50° resolution, Geosci. Model Dev., 16, 1163–1178, https://doi.org/10.5194/gmd-16-1163-2023, 2023.

Waldman, R., Herrmann, M., Somot, S., Arsouze, T., Benshila, R., Bosse, A.,nut, J., Giordani, H., Sevault, F., and Testor, P.: Impact of the Mesoscale Dynamics on Ocean Deep Convection: The 2012–2013 Case Study in the Northwestern Mediterranean Sea, J. Geophys. Res.-Oceans, 122, 8813–8840, https://doi.org/10.1002/2016jc012587, 2017.

Waldschläger, K., Born, M., Cowger, W., Gray, A., and Schüttrumpf, H.: Settling and rising velocities of environmentally weathered micro- and macroplastic particles, Environmental Research, 191, 110192, https://doi.org/10.1016/j.envres.2020.110192, 2020.

Weiss, L.: Model code and regional configuration for 'Modeling Indian Ocean circulation to study marine debris dispersion: insights into high-resolution and Stokes drift effects', figshare [code], https://doi.org/10.6084/m9.figshare.28903067, 2025.

Weiss, L., Ludwig, W., Heussner, S., Canals, M., Ghiglione, J. F., Estournel, C., Constant, M., and Kerhervé, P.: The missing
ocean plastic sink: Gone with the rivers, Science, 373, 107–111,
https://doi.org/10.1126/science.abe0290, 2021.

Weiss, L., Estournel, C., Marsaleix, P., Mikolajczak, G., Con-
stant, M., Ghiglione, J.-F., and Ludwig, W.: From source to
sink: part 2 – seasonal dispersion of microplastics discharged
in the NW Mediterranean Sea by the Rhone River in south-
ern France, Environmental Science and Pollution Research,
https://doi.org/10.1007/s11356-024-35364-6, 2024a.

Weiss, L., Estournel, C., Marsaleix, P., Mikolajczak, G., Constant, M., and Ludwig, W.: From source to sink: part 1 – characterization and Lagrangian tracking of riverine microplastics in the Mediterranean Basin, Environmental Science and Pollution Research, https://doi.org/10.1007/s11356-024-34635-6, 2024b.

Xu, Z. and Bowen, A. J.: Wave- and Wind-Driven Flow in Water of Finite Depth, Journal of Physical Oceanography, 24, 1850–1866, https://doi.org/10.1175/1520-0485(1994)024<1850:wavdfi>2.0.co;2, 1994.

Yang, C., Leonelli, F. E., Marullo, S., Artale, V., Beggs, H., Nardelli, B. B., Chin, T. M., De Toma, V., Good, S., Huang, B., Merchant, C. J., Sakurai, T., Santoleri, R., Vazquez-Cuervo, J., Zhang, H.-M., and Pisano, A.: Sea Surface Temperature Intercomparison in the Framework of the Copernicus Climate Change Service (C3S), Journal of Climate, 34, 5257–5283, https://doi.org/10.1175/jcli-d-20-0793.1, 2021.

Zhan, P., Guo, D., and Hoteit, I.: Eddy-Induced Transport and Ki-
netic Energy Budget in the Arabian Sea, Geophysical Research
Letters, 47, https://doi.org/10.1029/2020gl090490, 2020.

Zhang, L., Du, Y., and Cai, W.: A spurious positive Indian Ocean Dipole in 2017, Science Bulletin, 63, 1170–1172, https://doi.org/10.1016/j.scib.2018.08.001, 2018.

Zhang, Z., Wang, J., and Yuan, D.: Mixed Layer Salinity Balance in the Eastern Tropical Indian Ocean, J. Geophys. Res.-Oceans, 127, https://doi.org/10.1029/2021jc018229, 2022.

Zsoter, E., Harrigan, S., Barnard, C., Wetterhall, F., Ferrario, I., Mazzetti, C., Alfieri, L., Salamon, P., and Prudhomme, C.: River discharge and related historical data from the Global Flood Awareness System. v3.1. Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], https://doi.org/10.24381/CDS.A4FDD6B9, 2021.