# Future changes to the upper ocean Western Boundary Currents across two generations of climate models

Alex Sen Gupta, Annette Stellema, Gabriel Pontes, Andréa Taschetto, Adriana Vergés, Vincent Rossi

► To cite this version:

Alex Sen Gupta, Annette Stellema, Gabriel Pontes, Andréa Taschetto, Adriana Vergés, et al. Future changes to the upper ocean Western Boundary Currents across two generations of climate models. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-88934-w. hal-03235034

HAL Id: hal-03235034

https://hal.science/hal-03235034

Submitted on 26 May 2021

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

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

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

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

OPEN

# Future changes to the upper ocean Western Boundary Currents across two generations of climate models

Alex Sen Gupta<sup>1,2,3✉</sup>, Annette Stellema<sup>1,2,3</sup>, Gabriel M. Pontes<sup>4</sup>, Andréa S. Taschetto<sup>1,2</sup>, Adriana Vergés<sup>3,5</sup> & Vincent Rossi<sup>6</sup>

Western Boundary Currents (WBCs) are important for the oceanic transport of heat, dissolved gases and nutrients. They can affect regional climate and strongly influence the dispersion and distribution of marine species. Using state-of-the-art climate models from the latest and previous *Climate Model Intercomparison Projects*, we evaluate upper ocean circulation and examine future projections, focusing on subtropical and low-latitude WBCs. Despite their coarse resolution, climate models successfully reproduce most large-scale circulation features with ensemble mean transports typically within the range of observational uncertainty, although there is often a large spread across the models and some currents are systematically too strong or weak. Despite considerable differences in model structure, resolution and parameterisations, many currents show highly consistent projected changes across the models. For example, the East Australian Current, Brazil Current and Agulhas Current extensions are projected to intensify, while the Gulf Stream, Indonesian Throughflow and Agulhas Current are projected to weaken. Intermodel differences in most future circulation changes can be explained in part by projected changes in the large-scale surface winds. In moving to the latest model generation, despite structural model advancements, we find little systematic improvement in the simulation of ocean transports nor major differences in the projected changes.

Anthropogenic climate change manifests as increases in surface temperature and sea level, rainfall distribution changes and increasing frequency and intensity of certain extreme events<fnref n="1" />. An often underappreciated aspect of climate change relates to upper ocean circulation. Surface ocean currents play a crucial role in the redistribution of heat², pollutants³, plastics⁴, biogeochemical tracers⁵ and species dispersal—particularly during passive early-life stages (e.g. eggs, larvae), influencing the distribution of marine species⁶.

Particularly important for heat transport, air-sea interactions and marine ecosystems are the energetic Western Boundary Currents (WBCs); narrow (~100–200 km), intense (~1 m/s) jets in the upper several hundred meters off the continental shelves in the western basins. Subtropical WBCs form the poleward flanks of the subtropical gyres. Their poleward extensions form high-variability regions where the WBCs often leave the coast. Equatorward flowing WBCs at low latitudes play an important role in tropical circulation, facilitating cross-equatorial and inter-basin connectivity.

The large-scale upper-ocean circulation is primarily driven by surface winds. In the subtropics, north-south surface wind gradients, associated with climatological high-pressure systems, drive an equatorward *Sverdrup* transport in the ocean interior⁷. To conserve potential vorticity this flow is largely compensated by intense, poleward WBCs forming partially closed subtropical gyres. In the Indian and Pacific basins, however, some water from the Agulhas Current (AC) and East Australian Current (EAC) leak from the subtropical gyres, escaping westward south of Africa (AC extension) and Australia (Tasman Leakage), respectively. Low-latitude WBCs (LLWBCs) also play a role in the leakage of water from the Pacific to Indian basins: the Mindanao Current

<footnote n="1">¹Climate Change Research Centre, University of New South Wales, Sydney, Australia. ²Australian Research Council Centre of Excellence for Climate Extremes, University of New South Wales, Sydney, Australia. ³Centre for Marine Science and Innovation, University of New South Wales, Sydney, Australia. ⁴Institute of Oceanography, University of São Paulo, São Paulo, Brazil. ⁵Centre for Marine Science & Innovation and Evolution & Ecology Research Centre, School of Biological, Earth and Environmental Sciences, UNSW Australia, Sydney, NSW, Australia. ⁶Mediterranean Institute of Oceanography (UM 110, UMR 7294), CNRS, Aix Marseille Univ., Univ. Toulon, IRD, 13288 Marseille, France. ⁷email: a.sengupta@unsw.edu.au</footnote>
(MC) and, to a lesser extent, the New Guinea Coastal Undercurrent (NGCU), feed the climatically important Indonesian Throughflow (ITF)⁸. The intense WBCs are typically associated with high levels of mesoscale activity that alter the characteristics of the WBCs; for example WBCs tend to be faster, narrower and better located in eddy-permitting/resolving models compared to models that cannot resolve mesoscale activity⁹.

Climate models project consistent future changes to the surface winds¹ including a poleward expansion of the Hadley Circulation¹⁰, a poleward intensification of the subtropical westerlies¹¹, and a slowdown of the Walker circulation and equatorial trade winds¹². Such changes would modify the upper ocean circulation. Using multiple lines of evidence, previous work¹³ has identified a centennial scale poleward intensifications of the subtropical WBCs (excluding the Gulf Stream, GS), although natural variability can mask these changes even on multi-decadal timescales¹⁴,¹⁵. These studies link the WBC changes to the poleward intensification of the subtropical westerlies and the expansion of the Hadley Cell in both hemispheres, and the slowdown of the Atlantic meridional overturning circulation in the case of the GS. The intensified heat transport and/or poleward shift of these WBCs have led to hotspots that have warmed 2–3 times faster than the global average¹⁶,¹⁷. This warming has increased the occurrence of marine heatwaves¹⁸ and affected larval dispersal, causing dramatic modifications to marine ecosystems. For example, the range expansion of some consumers has led to extensive overgrazing of kelp forests, with major cascading impacts on entire ecosystems¹⁹,²⁰.

The intensification of WBCs has been associated with the 'tropicalisation' of nearshore marine communities as warm-affinity species expand their distribution poleward²¹, impacting ecosystem functions like productivity²². Shifts in WBCs have already impacted fisheries in meaningful ways. For instance, a northward shift in the GS has been linked to the collapse of the Gulf of Maine cod fishery²³. Future climate-induced changes in ocean circulation are also projected to alter coastal retention and larval dispersal pathways²⁴.

Previous studies using ocean and coupled climate models from previous generations of the Coupled Model Intercomparison Project (CMIP) suggest robust future circulation changes in different basins¹⁷,²⁵⁻²⁷. Using the new CMIP6 ensemble, this study examines major WBCs and associated circulation features, including model representation of transport and seasonality in comparison to observations and reanalysis products, projected circulation changes and connections to surface wind changes. We also examine any changes compared to the CMIP5 model ensemble.

## Results

In the following, model transports and projected changes are expressed as ensemble interquartile ranges with individual model details provided in Tables S1-S10, reanalysis transport estimates are provided as the range across the three products examined (Table S11), and observational transports and associated references are provided in Table S12. Projections represent differences between the 1900–2000 *historical* mean and 2050–2100 means from the business-as-usual SSP5-8.5 (RCP8.5) scenarios for CMIP6 (CMIP5).

### Indian Ocean.
In the Indian Ocean, the South Equatorial Current (SEC) forms the northern limb of the subtropical gyre, carrying fresh ITF water to the western basin. The SEC bifurcates east of Madagascar, forming the Northeastern and Southeastern Madagascar Currents (NMC and SMC, Fig. 1). Along the African shelf, the NMC further splits southward through the Mozambique Channel (MZC) and northwards as the East African Coastal Current (EACC). Further south, the SMC and MZC transport combine into the Agulhas Current (AC). The AC extension continues westwards beyond the African cape where it retroflects, returning most water eastwards to the Indian basin³³, while a part of this water (~21Sv) escapes into the South Atlantic as Agulhas Leakage.

The EACC transport across the CMIP6 models (interquartile range: 16.5–19.9 Sv Fig. 2, Fig. 3a) is consistent with the observed 19 Sv peak near 5°S and lies within the broad range of reanalysis estimates (7.4–23 Sv). The simulated NMC (19.4–22.7 Sv) and SMC (–10.2 to –15.7 Sv) are generally weaker than the range of observations (27–48 Sv and 20–30 Sv, respectively), based on multiple short-term estimates (Table S12), but span similar ranges to the reanalysis (Fig. 3). Conversely, the simulated transports through the MZC (17–24.6 Sv) are slightly stronger than observations and reanalysis (15–19 Sv and 11.8–21 Sv, respectively). The simulated MZC transport seasonality, which is maximum around austral autumn, agrees well with observations and reanalysis (Figure S1). Further south, the CMIP6 AC transport increases to 50.8–61.6 Sv near Africa's southern tip, somewhat weaker than the observational (70–77 Sv) but overlapping the weaker reanalysis estimates (47.2–53.7 Sv). A recent 3-year campaign³⁴ found AC transport at ~27°E to be strongest in austral summer and weakest in winter, although large interannual variability was evident. This seasonality is qualitatively consistent with the models and reanalysis, although the observed seasonal range ~15 Sv is considerably larger than in the models ~3 Sv (Figure S1).

Previous work²⁸ showed a broad-scale projected slowdown of the south Indian Ocean circulation by the end of the twenty-first century in CMIP5 models. Their reported weakening of both the western boundary Agulhas system and eastern boundary Leeuwin Current system is consistent with our CMIP6 and CMIP5 results. There is near-unanimous agreement across CMIP6 for reduced transport for the MZC (3.3 to 5.3 Sv), SMC (0.9 to 1.8 Sv), NMC (–2.1 to –3.9 Sv) and AC (3.4–7.6 Sv) (Fig. 2, Fig. 3b). However, neither the CMIP5 nor CMIP6 models show a consistent change in the EACC. In contrast to the reduction in transport along much of the southern African coast, the westward flowing AC extension south of Africa intensifies in all models (–3.6 to –7.7 Sv at 25°E) – a ~15% strengthening.

### Atlantic Basin.
At the northern extent of the South Atlantic subtropical gyre, the westward SEC bifurcates with most of its water entering the equatorward North Brazil Current (NBC)—responsible for large upper-ocean cross-equatorial heat transport³⁵. The remainder flows southward from ~10°S forming the relatively weak Brazil Current (BC). In the North Atlantic, the poleward flow, partly fed by the NBC, follows the western boundary of
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Schematic showing projected changes in WBC transport. Background colours show the multi-model mean projected change in sea surface temperature divided by the global mean change, e.g. 150% implies a warming rate 1.5× the global average.</figcaption></figure>

Central America as the Caribbean and Yucatan Currents ultimately emerging via the Florida Straits to form the GS. The GS breaks away from the coast at ~40°N, feeding the north-eastward North Atlantic Current.

BC transport estimates from observations range from -19 and -23 Sv between 36 and 38°S. CMIP6 models generally simulate the maximum BC transport between about 35-40°S with values ranging from -13.9 to -25.8 Sv, which lies in the very broad range of reanalysis transports (-8.7 to -32.5 Sv). The observed NBC transport (23-26 Sv) is slightly underestimated by the ensemble (19-22 Sv, 5-10°S) with even weaker estimates from reanalysis (11.5-18 Sv). The models simulate maximum (minimum) transport in July-Aug (April-May, Fig. 4) (observed NBC seasonality estimates are not available at the latitudes examined). The BC and NBC forms at ~10°S³⁷ just north of the basin-averaged zero wind-stress curl latitude. This bifurcation typically sits about 10° too far south in the CMIP models, in part related to a systematic southward bias in the model Atlantic wind field (Fig. 2, Figure S2).

In the Northern Hemisphere, the complex circulation of the Caribbean Sea and Gulf of Mexico is represented very differently across the coarse resolution models. Compared to observations, most models (and reanalysis) underestimate the LLWBC transport of the Yucatan Current (30 Sv) with a model range of 13.5-23.3 Sv (reanalysis: 8.5-25.6 Sv). The GS transport intensifies moving northwards (>90 Sv) where it diverges from the coast, with the strongest transport occurring in boreal fall³⁸. This northward intensification is absent in the models and reanalysis: northward transports peaks at 38 to 42.6 Sv between about 28-33°N (reanalysis: 37.1-46.7 Sv). The simulated winter intensification of the GS is consistent from the western margin to the extension region (Figure S1). However, the models are generally out of phase with the observations that indicate maximum transports during summer at 26.5°N³⁹ and in the extension region³⁸. While the reanalysis seasonality matches the models along the coast, there is poor agreement in the extension region.

In the Southern Hemisphere, the BC is projected to intensify (4.2 to 6.0 Sv), especially south of 30°S, associated with an increased northward basin interior transport²⁷. This intensification is consistent with intensified westerlies across the Indian Ocean basin (Figure S3, Figure S4), which can increase northward Ekman transport and intensify the Indian Ocean input to the Atlantic via Agulhas Leakage³¹.

Conversely, WBC transports weaken northwards of ~15°S. The cross-equatorial NBC flow is projected to weaken (-1.7 to -4.7 Sv). Similarly, the GS reduces at all latitudes with a -4.9 to -10.8 Sv (~15%) decrease around the GS maximum. These changes are poorly explained by surface wind changes and are likely associated with a weakened Atlantic Meridional Overturning Circulation (see below).

**Pacific Basin.** In the South Pacific, the broad westward SEC bifurcates at the Australian margin forming the poleward EAC and equatorward Gulf of Papua Current (GPC). The EAC partially separates from the coast near 30°S forming the Tasman Front, which continues southward to the east of New Zealand as the East Auckland Current and East Cape Currents (ECC). The remaining EAC water feeds a series of eddies that move southwards, forming the EAC extension and Tasman Leakage that provides a high-latitude pathway of water to the Indian
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Historical meridional transport (left panels) and projected meridional transport change (right panels) by latitude along western boundaries shown in the map. Red/blue/green lines are multi-model median transport or transport change for CMIP6(SSP5-8.5)/CMIP5(RCP8.5)/CMIP6(SSP1-2.6) scenarios, associated shading indicates interquartile range (for high emission scenarios only). For projection panels lines are thickened where the multi-model median change is significant at the 95% level based on a two-sided Wilcoxon signed rank test. Black vertical lines and black polygons in the central map (along the WBC paths) show the location for the zonal and meridional transports presented in Fig. 3.</figcaption></figure>

Basin. The northward flowing GPC feeds the NGCU that exits the Northern Solomon Sea via multiple straits providing water to the subsurface Equatorial Undercurrent (EUC)¹⁷. In the Northern Hemisphere, the MC also feeds the EUC and forms the primary source of the ITF that transports warm water into the Indian Ocean¹⁷. Further north, the Kuroshio Current (KC) extends northwards from ~15°N along eastern Japan, where it eventually separates and continues eastward.

The observed EAC transport reaches about -22 Sv at 27°S¹⁰, with maximum/minimum transport in austral winter/summer⁴¹. CMIP6 transports are generally similar in strength (-20.3 to -23.4 Sv) and seasonality (Figure S1) to observations. While the seasonality is similar for the reanalysis products, they tend to underestimate the transport (-7 to -17 Sv). The EAC extension transport (~7 Sv) and Tasman Leakage (~8 Sv) are, however, systematically underestimated in the models (-1.4 to -7 Sv and -0.9 to -4 Sv, respectively), with some models simulating an EAC extension with northward mean flow, related to a poor representation of regional winds²⁵.

In the Northern Hemisphere, the KC intensifies from about 15 Sv at 18°N to over 20 Sv between 25 and 30°N. The CMIP6 models systematically overestimate the transport with an interquartile range of 30.1 to 44.2 Sv, which encompasses the reanalysis estimates of 38.6 to 39 Sv. Observations suggest that KC strength is weakest during winter, to the east of Taiwan⁴² while models and reanalysis display minimum transports earlier in autumn (Figure S1). Further north (28°N) observed transport is minimum in autumn⁴³, while the models show no distinct seasonality. In the extension region, surface transport is weakest in winter/spring and strongest in summer/autumn⁴⁴; while the model transports tend to peak in spring.
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Seasonal cycle of mean transport (upper panels) and projected change (lower panels) for selected currents, where the annual mean transports have been removed. Red line/shading indicate multi-model median/interquartile range for CMIP6 models; blue line/shading/dashed line indicate multi-model median/interquartile range/interdecile range for CMIP5 models. Grey shading in upper panels indicates the range of three ocean reanalysis.</figcaption></figure>
ocean projection<sup>29</sup> could not be explained by regional wind changes. Instead, these studies found that the changes are related to a slowdown in deep ocean waters entering the South Pacific.

For the majority of currents examined across the basins, there is no significant difference in ensemble mean historical transports between CMIP5 and CMIP6. Only the Tasman Leakage and SMC demonstrate significantly different MMM transports (Table S13). In the case of the Tasman Leakage the MMM flow reverses direction from weakly eastwards in CMIP5 to weakly westward in CMIP6. This constitutes an important regional improvement although 20% of models still have spurious eastward flow in CMIP6 (compared to 57% for CMIP5). Only the MMM projected change in the AC extension transport is significantly different between model ensembles, with the CMIP6 suggesting a 40% smaller intensification compared to CMIP5.

**Seasonal changes.** As described, many currents exhibit seasonal transport changes that are consistent across models (Figure S1). While seasonal timing is realistic for many currents, some simulated currents, for example the GS and Kuroshio system, poorly simulate the observed seasonal phase or amplitude. For some currents, comparison is hampered by uncertainties in the observed seasonality due to short observational records and large internal variability<sup>37,47,48</sup>. Both basin wide and local winds are important in setting transport seasonality, although the influence of remote winds may be lagged due to the slow propagation of ocean waves. The seasonal phase of wind stress curl in the models is broadly similar to the ERA5 reanalysis, although large discrepancies are evident, particularly near the boundaries of regions with strong seasonality differences (Figure S5a-c). We note that CMIP transports that show poor agreement with the observations or reanalysis seasonality (e.g. GS and KC extensions, NBC) are often associated with substantial biases in wind stress curl seasonality in the regions extending eastwards of the WBCs (Figure S5c).

A subset of currents also exhibits consistent projected changes in transport seasonality across both model generations (Fig. 4). In the South Indian Ocean, most models project a reduced seasonal cycle for the NMC and SMC (Figure S1). Likewise, in the South Atlantic, the models consistently simulate a substantial weakening of NBC seasonality. In the South Pacific, both the EAC extension and Tasman Leakage show an amplification in seasonality. Conversely, east of New Zealand, the seasonality of the ECC is projected to decrease. In the North Pacific, the seasonality in meridional transport where the KC separates from the coast shows consistent increases (decreases) in boreal winter (summer).

These currents with modified seasonality generally occur at latitudes where the phase of the projected wind stress curl seasonality also show large projected changes (except for the MC where the changes are just upstream of the MC latitudes examined; Figure S5e). These projected changes in wind stress curl seasonality are zonally oriented, occurring at transition zones where the historical wind stress curl seasonality changes rapidly with latitude (Figure S5b), suggesting that the changes are associated with a poleward expansion of the wind fields, a well-established consequence of anthropogenic climate change<sup>10</sup>, and their associated seasonality (Figure S5f.).

**Emergent constraints.** For a subset of currents, there appears to be a significant inter-model relationship between historical and projected transports (Figure S6). These relationships may provide emergent constraints to narrow the uncertainty associated with the large spread in projections, although such constraints may be biased by common structural errors<sup>49</sup>. For example, for the EAC extension and Tasman Leakage, models that underestimate mean transport or that have flow in the wrong direction tend to project the largest increases in southward or westward flow, respectively. Given observed EAC extension transports (~7 Sv, Table S12), a more moderate future change (~5 Sv) may therefore be more credible than the more extreme changes projected by some models. Similarly, given an observed mean transport of ~15 Sv, it is likely that the ITF decrease would be 3-4 Sv rather than the more extreme model estimates.

**Connections to surface wind changes.** Neglecting friction, non-linear processes and interactions with deep ocean circulation, the depth-integrated meridional transport away from the western boundaries can be related to gradients in the surface wind field via Sverdrup dynamics<sup>7</sup>. In particular, a positive (negative) wind stress curl drives northward (southward) flow in the ocean interior. WBCs provide a return flow for much of this meridional transport and a significant part of the inter-model differences in the historical mean WBC transport can be related to differences in interior transport (Table 1, Figure S7). The offset between the WBC-interior regression lines and the one-to-one line in Figure S7 for certain currents relate to inter-basin leakage or flow compensation in the deep ocean as part of the overturning circulation. For example, the ~10 Sv offsets for the EAC, EAC extension, GPC and NGCU result from a leakage of water via the ITF. The offset is ~5 Sv less than the ITF transport as there is also a net upwelling into the upper Pacific from below 1000 m. The offsets for the Atlantic basin currents, including the GS and NBC, result from the deep return flow below 1000 m.

For many currents intermodel difference in interior transports can be explained to some degree by differences in the surface wind field via Sverdrup dynamics (Table 1). As a result, up to 50% of the intermodel variance in WBC transports can be related to differences in the surface winds. Other factors, including different overturning rates, different inter-basin transports and non-linear dynamics must be invoked to explain the wide range of mean WBC transports.

Similarly, a significant fraction of intermodel projected WBC differences can be related to changes in surface wind stress curl, for most WBCs investigated (Table 1). In general, WBC whose mean differences are well explained by differences in their surface winds tend to be those whose projected transport differences are also well explained by differences in surface wind changes. The particularly poor relationship noted for the BC probably relates to the fact that the Sverdrup calculation becomes poorly defined as the eastern boundary lies at the southern tip of Africa. Other weak relationships in the Atlantic likely stem from large projected changes in the
<table><caption>**Table 1.** Correlation between: interior (to the east of the WBC) and WBC transport (column 2), interior and derived Sverdrup transport (column 3), WBC and Sverdrup transport (column 4). Associated correlations for projected changes shown in columns 4, 5 and 6. Outliers (values exceeding 3 × scaled median deviations) are removed prior to the calculation of correlations. <sup>+1</sup>EAC extension includes transport to the east of New Zealand. <sup>+2</sup>MAD includes the WBC to the east and west of Madagascar. Scatter plots of Interior vs WBC and interior vs Sverdrup transport for the combined CMIP5 & 6 ensemble shown in Figure S7. Bold correlations indicate significant correlations at 95% level, based on Spearman Rank correlation.</caption><thead><tr><th rowspan="2"></th><th colspan="3">Historical</th><th colspan="3">Projected change</th></tr><tr><th>Interior vs WBC</th><th>Interior vs Sverdrup</th><th>WBC vs Sverdrup</th><th>Interior vs WBC</th><th>Interior vs Sverdrup</th><th>WBC vs Sverdrup</th></tr></thead><tbody><tr><td>EAC core</td><td>-0.91</td><td>0.76</td><td>-0.68</td><td>-0.81</td><td>0.49</td><td>-0.58</td></tr><tr><td>EAC extension<sup>+1</sup></td><td>-0.87</td><td>0.85</td><td>-0.72</td><td>-0.96</td><td>0.8</td><td>-0.7</td></tr><tr><td>MC</td><td>-0.82</td><td>0.53</td><td>-0.66</td><td>-0.91</td><td>0.6</td><td>-0.66</td></tr><tr><td>KC South</td><td>-0.93</td><td>0.47</td><td>-0.34</td><td>-0.93</td><td>0.76</td><td>-0.79</td></tr><tr><td>KC Central</td><td>-0.98</td><td>0.2</td><td>-0.15</td><td>-0.94</td><td>0.41</td><td>-0.35</td></tr><tr><td>KC North</td><td>-0.77</td><td>0.62</td><td>-0.46</td><td>-0.98</td><td>0.76</td><td>-0.77</td></tr><tr><td>GPC</td><td>-0.95</td><td>0.49</td><td>-0.45</td><td>-0.86</td><td>0.83</td><td>-0.6</td></tr><tr><td>NGCU</td><td>-0.76</td><td>0.75</td><td>-0.75</td><td>-0.92</td><td>0.89</td><td>-0.76</td></tr><tr><td>BC</td><td>-0.95</td><td>0.34</td><td>-0.32</td><td>-0.78</td><td>0.22</td><td>-0.05</td></tr><tr><td>NBC</td><td>-0.42</td><td>0.84</td><td>-0.38</td><td>-0.2</td><td>0.86</td><td>-0.35</td></tr><tr><td>GS</td><td>-0.54</td><td>0.57</td><td>-0.04</td><td>-0.78</td><td>0.48</td><td>-0.35</td></tr><tr><td>AC</td><td>-0.8</td><td>0.55</td><td>-0.26</td><td>-0.85</td><td>0.42</td><td>-0.29</td></tr><tr><td>MAD<sup>+2</sup></td><td>-0.58</td><td>0.61</td><td>-0.45</td><td>-0.58</td><td>0.57</td><td>-0.31</td></tr><tr><td>EACC</td><td>-0.75</td><td>0.53</td><td>-0.48</td><td>-0.8</td><td>0.72</td><td>-0.76</td></tr></tbody></table>

Atlantic overturning circulation<sup>50</sup>. Indeed, projected NBC decreases in CMIP5 are largely compensated by a weakening of North Atlantic Deep Water transport<sup>27</sup>.

**Near-surface transport.** WBCs affect the distribution of marine species via the dispersal of early-life stages and modulation of local thermal regimes<sup>24,51</sup>. However, ecosystem impacts will be most sensitive to near-surface circulation changes within the euphotic zone where most marine life thrives. As such, we also examine WBC transport changes in the top 100 m of the water column.

For most currents examined, the change in the near-surface flow is of the same sign as the 1000 m integrated transport. An exception is the KC system, where the full-depth WBC is projected to weaken slightly along most of its length, while the near-surface flow is projected to intensify weakly north of 25°N (Figure S8, Figure S9). Previous work suggested that this intensification is associated with differences in warming rates across the KC, leading to an enhanced baroclinic flow<sup>52</sup>. As a tight connection between the state of the KC and the regional marine food webs has been documented<sup>53</sup>, this surface intensification may have consequences for the ecosystem. In contrast, the full-depth MC, which is projected to weaken, typically intensifies near the surface south of 7°S.

In general, when the direction of a WBC is aligned with the warming signal (e.g. in the subtropics), a poleward intensified WBC will assist species dispersal at poleward range edges. In contrast, when the WBC flow opposes climate change velocities (e.g. in the tropics), strengthening would hinder dispersal at the poleward edges with greater propagule dispersal at the warming, equatorward edges<sup>24</sup>. Weakened/strengthened WBCs are likely to directly modify larval transport and thermal regimes, affecting rates of poleward range shifts<sup>51</sup>. In addition, other more subtle changes such as WBC broadening or modified coastal retention or dispersal pathways may also impact marine life<sup>24</sup>.

**Summary and conclusions**

We have examined projected changes to the upper ocean western boundary circulation and associated currents. Many circulation features examined exhibit consistent projected changes across climate models (Fig. 1) despite large differences in model structure, e.g. bathymetry (Figure S10), resolution, parameterisation of unresolved processes. This suggests common large-scale drivers of change. Comparing CMIP5 and CMIP6, we found little robust differences in either mean current strengths (including their seasonality) or their projected changes between the model generations. A significant fraction of the inter-model differences in the projected changes can be related to differences in projected wind stress curl changes through linear dynamics. As such, the strong similarity between regional surface wind changes from CMIP5 and CMIP6 (Figure S3, Figure S4) would suggest little inter-generational difference in the circulation projections. However, other model differences including inter-basin exchanges (e.g. via the ITF) and overturning circulation changes (e.g. the slowdown of the Atlantic overturning circulation) will also be important.

Basin-scale changes in the wind field can drive changes to the WBCs on the multi-decadal timescales required for planetary waves to propagate across basins and establish a steady state. However, in a few regions, we also found consistent seasonality changes. These changes likely result from regional rather than basin-scale surface
winds change, particularly at higher latitudes, given the months to years for oceanic planetary waves triggered by remote winds to propagate to the western boundaries or from changes to the overturning circulation. While we did not perform a detailed analysis, currents with projected seasonality changes were often located at or to the west of regions with substantial projected wind stress curl seasonal phase changes. These changes likely relate to seasonal wind regimes shifting poleward

Such seasonal changes may be particularly relevant for understanding how future climate will impact marine life as many key biological processes (e.g. spawning, migrations) exhibit marked phenology. Circulation driven changes in temperature seasonality are also likely to be important. For example, winter, rather than annually averaged, temperatures in the Gulf Stream best explains the inter-annual variation in species' distribution and biomass⁵⁴.

It is interesting to note that the projected changes along the WBCs do not typically align with the historical mean currents: it would be wrong to conceptualise these circulation changes as increases or decreases in the speed of a conveyor belt. For example, the core EAC exhibits little future change, while the EAC extension shows large consistent increases in transport. The main Agulhas system shows a similar projected decrease in transport along its length, while south of Africa, the AC extension intensifies.

A major criticism of coarse resolution climate models is that WBCs are poorly simulated e.g. flow speeds are too weak and too wide. Furthermore, WBCs are usually associated with strong mesoscale eddy activity not resolved in CMIP models. Despite this, simulations carried out in eddy-resolving or eddy-permitting ocean models find projected transport changes in the Agulhas system³¹, EAC system³⁰,⁵⁵, ITF²⁹ and globally²⁴ that are broadly consistent with the CMIP models. This suggests that mesoscale processes are of secondary importance and that changes in basin-scale winds³⁰,³¹, or deep circulation in some cases like the ITF²⁵,²⁹,⁵⁶, explain a large component of the projected changes in transport. Indeed, while mesoscale processes affect certain characteristics of WBCs, transport projections for these narrow WBCs may be more reliable than projections for other variables that depend critically on local processes. This is not to say that other aspects of WBC changes may be unrealistic, including changes in current speeds or changes in eddy activity. Indeed, previous eddy-resolving ocean-biogeochemical model simulation⁵⁵ found projected increases in Tasman Sea eddy activity, nutrient pumping and productivity associated with the enhanced EAC extension, contrary to projected changes in a coarse resolution control model. It will be important to examine how the projections from coarse resolutions models discussed here differ from the high-resolution climate projections from HighResMIP⁵⁷ in progress (at the time of writing no output is available for the second half of the twenty-first century).

These future circulation changes could have important implications for ecosystems and the fisheries that rely on them. The intensification of the EAC extension has already caused profound shifts in the distribution of fish⁵⁸ and almost complete giant kelp forests losses, with major impacts on valuable fisheries⁵⁹. Conversely, new fisheries are emerging as a result of range expansions of high-value species⁶⁰. Future intensified warming and circulation increases in this region are likely to exacerbate these changes, with implications for conservation and management⁶.

Climate-associated shifts in stocks of commercial species, such as those related to WBC changes, pose major governance challenges as species move beyond national and international boundaries, which can lead to conflict over shared resources⁶¹. Improving our understanding of future WBC changes will be key to credible forecasts of the distribution of commercially-important species, enabling the development of governance solutions⁶². While CMIP-class models may provide credible forecasts of transport changes, the associated changes in current speeds, eddy activity and cross shelf processes will require the use of global or regional models that can represent mesoscale and smaller processes.

## Materials and methods

We use zonal and meridional ocean velocity, wind stress and sea surface temperature from 28 CMIP5 (Table S1) and 25 CMIP6 (Table S2) models. Meridional western boundary transports are integrated from the surface to 1000 m (or to 100 m) and from the coast to 7° offshore (15° offshore for the NGCU) except where the boundary currents do not run close to the coast. In such cases, we manually identified the western edge of WBCs by an examination of the historical mean meridional transport for each model. Unlike in the real ocean, WBCs in coarse resolution climate models can have widths of several hundreds of kilometres. Transports are relatively insensitive to the boundary current width used as WBC flow is generally strongest towards the shore (Figure S11). Zonal transports are calculated at selected longitudes to capture the flow associated with WBC extensions (e.g. AC extension, Tasman Leakage and the ITF). The northern and southern extents of the currents were manually identified based on examination of the historical mean meridional transport (for approximate locations for all sections, see Fig. 2).

We compare transports averaged over 1900 to 2000 (historical scenario) with averaged transports for 2050 to 2100 (RCP8.5 scenario for CMIP5; SSP5-8.5 for CMIP6). These long timeframes ensure that differences are primarily associated with the imposed anthropogenic forcing, as much of the intrinsic variability is removed when considering the long-term average. Given the large spread in the number of available ensemble members, our results are based on a single ensemble member from each model. Testing with selected models indicates negligible differences in climatological transports based on different ensemble members. We use the 'business as usual' scenario as it provides the highest signal to noise ratio. Results from lower emission scenarios produce commensurately smaller projected changes; multi-model median results for SSP1-2.6, which is commensurate with the policy relevant 2°C warming by the end of the century, are presented for comparison for a subset of analysis (Fig. 2, Fig. 3, Figure S8, Figure S9). In the manuscript text, we quote interquartile transport ranges across the CMIP6 ensemble, individual model transports for the upper 1000 m (and 100 m) and interquartile ranges for CMIP5 and CMIP6 are shown in Fig. 2, Fig. 3 and Table S1-10.
The significance of projected changes was tested using a non-parametric two-sided Wilcoxon signed rank test. Inferences based on these statistical tests likely overestimate the level of significance since models share common components and are so not truly independent of each other.

Where possible model transports are compared to estimates from three ocean reanalysis products (Table S11): Ocean Reanalysis System 5 (ORAS5)⁶³, Global Ocean Data Assimilation System (GODAS)⁶⁴ and CMCC Global Ocean Reanalysis System (C-GLORS)⁶⁵; and observations cited in the literature (Table S12). While the latter is not a like-for-like comparison as locations, integration depths and time periods vary, this provides a useful order of magnitude evaluation of the models. Model winds are compared against the ERA5 reanalysis⁶⁶.

All figures generated in MATLAB 2020b update 3. Labels, legends and schematics were done using Adobe Illustrator 25.2. Coastlines generated using m_map toolbox⁶⁷.

Received: 22 December 2020; Accepted: 19 April 2021

Published online: 05 May 2021

## References

1. IPCC. *Climate Change 2007 - The Physical Science Basis: Working Group I Contribution to the Fourth Assessment Report of the IPCC*. (Cambridge University Press, 2007).

2. Ganachaud, A. & Wunsch, C. Large-scale ocean heat and freshwater transports during the world ocean circulation experiment. *J. Clim.* **16**, 696–705 (2003).

3. Rossi, V., Van Sebille, E., Sen Gupta, A., Garon, V. & England, M. H. Multi-decadal projections of surface and interior pathways of the Fukushima Cesium-137 radioactive plume. *Deep-Sea Res. Part Oceanogr. Res. Pap.* **80**, 37–46 (2013).

4. Eriksen, M. *et al*. Plastic pollution in the South Pacific subtropical gyre. *Mar. Pollut. Bull.* **68**, 71–76 (2013).

5. Williams, R. G. & Follows, M. J. Physical Transport of Nutrients and the Maintenance of Biological Production. in *Ocean Biogeochemistry: The Role of the Ocean Carbon Cycle in Global Change* (ed. Fasham, M. J. R.) 19–51 (Springer, 2003). doi: https://doi.org/10.1007/978-3-642-55844-3_3.

6. Dubois, M. *et al*. Linking basin-scale connectivity, oceanography and population dynamics for the conservation and management of marine ecosystems. *Glob. Ecol. Biogeogr.* **25**, 503–515 (2016).

7. Sverdrup, H. U. Wind-driven currents in a Baroclinic Ocean; with application to the equatorial currents of the Eastern Pacific. *Proc. Natl. Acad. Sci. USA* **33**, 318–326 (1947).

8. Godfrey, J. S. The effect of the Indonesian throughflow on ocean circulation and heat exchange with the atmosphere: A review. *J. Geophys. Res. Oceans* **101**, 12217–12237 (1996).

9. Hewitt, H. T. *et al*. Resolving and parameterising the ocean mesoscale in earth system models. *Curr. Clim. Change Rep.* **6**, 137–152 (2020).

10. Hu, Y., Tao, L. & Liu, J. Poleward expansion of the hadley circulation in CMIP5 simulations. *Adv. Atmos. Sci.* **30**, 790–795 (2013).

11. Swart, N. C. & Fyfe, J. C. Observed and simulated changes in the Southern Hemisphere surface westerly wind-stress. *Geophys. Res. Lett.* **39**, L16711 (2012).

12. Vecchi, G. A. & Soden, B. J. Global warming and the weakening of the tropical circulation. *J. Clim.* **20**, 4316–4340 (2007).

13. Yang, H. *et al*. Intensification and poleward shift of subtropical western boundary currents in a warming climate. *J. Geophys. Res. Oceans* **121**, 4928–4945 (2016).

14. Curry, R. G. & McCartney, M. S. Ocean gyre circulation changes associated with the North Atlantic Oscillation. *J. Phys. Oceanogr.* **31**, 3374–3400 (2001).

15. Seager, R., Kushnir, Y., Naik, N. H., Cane, M. A. & Miller, J. Wind-driven shifts in the latitude of the Kuroshio-Oyashio extension and generation of SST anomalies on decadal timescales*. *J. Clim.* **14**, 4249–4265 (2001).

16. Wu, L. *et al*. Enhanced warming over the global subtropical western boundary currents. *Nat. Clim. Change* **2**, 161–166 (2012).

17. Hu, D. *et al*. Pacific western boundary currents and their roles in climate. *Nature* **522**, 299–308 (2015).

18. Oliver, E. C. *J. et al*. Longer and more frequent marine heatwaves over the past century. *Nat. Commun.* **9**, 1324 (2018).

19. Ling, S. D., Johnson, C. R., Ridgway, K., Hobday, A. J. & Haddon, M. Climate-driven range extension of a sea urchin: Inferring future trends by analysis of recent population dynamics. *Glob. Change Biol.* **15**, 719–731 (2009).

20. Wernberg, T. *et al*. Climate-driven regime shift of a temperate marine ecosystem. *Science* **353**, 169–172 (2016).

21. Verges, A. *et al*. The tropicalization of temperate marine ecosystems: Climate-mediated changes in herbivory and community phase shifts. *Proc. R. Soc. B-Biol. Sci.* **281**, (2014).

22. Vergés, A. *et al*. Tropicalisation of temperate reefs: Implications for ecosystem functions and management actions. *Funct. Ecol.* **33**, 1000–1013 (2019).

23. Pershing, A. *J. et al*. Slow adaptation in the face of rapid warming leads to collapse of the Gulf of Maine cod fishery. *Science* **350**, 809–812 (2015).

24. van Gennip, S. *J. et al*. Going with the flow: the role of ocean circulation in global marine ecosystems under a changing climate. *Glob. Change Biol.* **23**, 2602–2617 (2017).

25. Sen Gupta, A. *et al*. Future changes to the Indonesian Throughflow and Pacific circulation: The differing role of wind and deep circulation changes. *Geophys. Res. Lett.* https://doi.org/10.1002/2016GL067757 (2016).

26. Sen Gupta, A., Ganachaud, A., McGregor, S., Brown, J. N. & Muir, L. Drivers of the projected changes to the Pacific Ocean equatorial circulation. *Geophys. Res. Lett.* **39**, L09605 (2012).

27. Pontes, G. M., Gupta, A. S. & Taschetto, A. S. Projected changes to South Atlantic boundary currents and confluence region in the CMIP5 models: the role of wind and deep ocean changes. *Environ. Res. Lett.* **11**, 094013 (2016).

28. Stellema, A., Sen Gupta, A. & Taschetto, A. S. Projected slow down of South Indian Ocean circulation. *Sci. Rep.* **9**, 17705 (2019).

29. Feng, M., Zhang, X., Sloyan, B. & Chamberlain, M. Contribution of the deep ocean to the centennial changes of the Indonesian Throughflow. *Geophys. Res. Lett.* **44**, 2859–2867 (2017).

30. Oliver, E. C. *J. & Holbrook, N. J*. Extending our understanding of South Pacific gyre “spin-up”: Modeling the East Australian Current in a future climate. *J. Geophys. Res. Oceans* **119**, 2788–2805 (2014).

31. Biastoch, A. & Böning, C. W. Anthropogenic impact on Agulhas leakage. *Geophys. Res. Lett.* **40**, 1138–1143 (2013).

32. Duan, J., Chen, Z. & Wu, L. Projected changes of the low-latitude north-western Pacific wind-driven circulation under global warming. *Geophys. Res. Lett.* **44**, 4976–4984 (2017).

33. Lutjeharms, J. R. E. Three decades of research on the greater Agulhas Current. *Ocean Sci.* **3**, 129–147 (2007).

34. Daher, H., Beal, L. M. & Schwarzkopf, F. U. A new improved estimation of agulhas leakage using observations and simulations of lagrangian floats and drifters. *J. Geophys. Res. Oceans*. https://doi.org/10.1029/2019JC015753 (2020).

35. Beal, L. M., Elipot, S., Houk, A. & Leber, G. M. Capturing the transport variability of a western boundary jet: Results from the agulhas current time-series experiment (ACT). *J. Phys. Oceanogr.* **45**, 1302–1324 (2015).
36. da Silveira, I. C. A., de Miranda, L. B. & Brown, W. S. On the origins of the North Brazil Current. *J. Geophys. Res. Oceans* **99**, 22501–22512 (1994).

37. Johns, W. E. et al. Annual Cycle and Variability of the North Brazil Current. *J. Phys. Oceanogr.* **28**, 103–128 (1998).

38. Sato, O. T. & Rossby, T. Seasonal and low frequency variations in dynamic height anomaly and transport of the Gulf Stream. *Deep Sea Res. Part Oceanogr. Res. Pap.* **42**, 149–164 (1995).

39. Zhao, J. & Johns, W. Wind-driven seasonal cycle of the atlantic meridional overturning circulation. *J. Phys. Oceanogr.* **44**, 1541–1562 (2014).

40. Sloyan, B., Ridgway, K. & Cowley, R. The East Australian Current and Property Transport at 27 S from 2012–2013. *J. Phys. Oceanogr.* **46**, 160108151222003 (2016).

41. Ridgway, K. R. & Godfrey, J. S. Seasonal cycle of the East Australian Current. *J. Geophys. Res. Oceans* **102**, 22921–22936 (1997).

42. Hsin, Y.-C., Qiu, B., Chiang, T.-L. & Wu, C.-R. Seasonal to interannual variations in the intensity and central position of the surface Kuroshio east of Taiwan: VARIATIONS OF KUROSHIO EAST OF TAIWAN. *J. Geophys. Res. Oceans* **118**, 4305–4316 (2013).

43. Wei, Y., Pei, Y. & Zhang, R.-H. Seasonal variability of the Kuroshio Current at the PN Section in the East China Sea based on in-situ observation from 1987 to 2010. *Acta Oceanol. Sin.* **34**, 12–21 (2015).

44. Yang, Y. & San Liang, X. On the Seasonal Eddy Variability in the Kuroshio Extension. *J. Phys. Oceanogr.* **48**, 1675–1689 (2018).

45. Sprintall, W., Jijffels, S. E., Molcard, R. & Jaya, I. Direct estimates of the Indonesian Throughflow entering the Indian Ocean: 2004–2006. *J. Geophys. Res. Oceans* **114**, n/a-n/a (2009).

46. Hill, K. L., Rintoul, S. R., Ridgway, K. R. & Oke, P. R. Decadal changes in the South Pacific western boundary current system revealed in observations and ocean state estimates. *J. Geophys. Res. Oceans* **116**, 5926 (2011).

47. Beal, L. M. et al. Five years of Florida Current structure and transport from the Royal Caribbean Cruise Ship Explorer of the Seas. *J. Geophys. Res. Oceans* **113**, 4154 (2008).

48. Garzoli, S. L. North Brazil Current retroflection and transports. *J. Geophys. Res.* **109**, C01013 (2004).

49. Bracegirdle, T. J. & Stephenson, D. B. On the Robustness of Emergent Constraints Used in Multimodel Climate Change Projections of Arctic Warming. *J. Clim.* **26**, 669–678 (2012).

50. Weijer, W., Cheng, W., Garuba, O. A., Hu, A. & Nadiga, B. T. CMIP6 Models Predict Significant 21st Century Decline of the Atlantic Meridional Overturning Circulation. *Geophys. Res. Lett.* **47**, e2019GL086075 (2020).

51. Wilson, L. J. et al. Climate-driven changes to ocean circulation and their inferred impacts on marine dispersal patterns: Ocean circulation dispersal patterns. *Glob. Ecol. Biogeogr.* **25**, 923–939 (2016).

52. Chen, C., Wang, G., Xie, S.-P. & Liu, W. Why Does Global Warming Weaken the Gulf Stream but Intensify the Kuroshio?. *J. Clim.* **32**, 7437–7451 (2019).

53. Yatsu, A. et al. Climate forcing and the Kuroshio/Oyashio ecosystem. *ICES J. Mar. Sci.* **70**, 922–933 (2013).

54. Morley, J. W., Batt, R. D. & Pinsky, M. L. Marine assemblages respond rapidly to winter climate variability. *Glob. Change Biol.* **23**, 2590–2601 (2017).

55. Matear, R. J., Chamberlain, M. A., Sun, C. & Feng, M. Climate change projection of the Tasman Sea from an Eddy-resolving Ocean Model. *J. Geophys. Res. Oceans* **118**, 2961–2976 (2013).

56. Sun, S. & Thompson, A. F. Centennial Changes in the Indonesian Throughflow Connected to the Atlantic Meridional Overturning Circulation: The Ocean's Transient Conveyor Belt. *Geophys. Res. Lett.* **47**, (2020).

57. Haarsma, R. J. et al. High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6. *Geosci. Model Dev.* **9**, 4185–4208 (2016).

58. Last, P. R. et al. Long-term shifts in abundance and distribution of a temperate fish fauna: a response to climate change and fishing practices. *Glob. Ecol. Biogeogr.* **20**, 58–72 (2011).

59. Barange, M. et al. *Impacts of climate change on fisheries and aquaculture: synthesis of current knowledge, adaptation and mitigation options* (2018).

60. Champion, C., Hobday, A. J., Tracey, S. R. & Pecl, G. T. Rapid shifts in distribution and high-latitude persistence of oceanographic habitat revealed using citizen science data from a climate change hotspot. *Glob. Change Biol.* **24**, 5440–5453 (2018).

61. Pinsky, M. L. et al. Preparing ocean governance for species on the move. *Science* **360**, 1189–1191 (2018).

62. Rogers, L. A. et al. Shifting habitats expose fishing communities to risk under climate change. *Nat. Clim. Change* **9**, 512–516 (2019).

63. Zuo, H., Balmaseda, M. A., Tietsche, S., Mogensen, K. & Mayer, M. The ECMWF operational ensemble reanalysis–analysis system for ocean and sea ice: A description of the system and assessment. *Ocean Sci.* **15**, 779–808 (2019).

64. Behringer, D. & Xue, Y. Evaluation of the global ocean data assimilation system at NCEP: The Pacific Ocean. /paper/EVALUATION-OF-THE-GLOBAL-OCEAN-DATA-ASSIMILATION-AT-Behringer-Xue/2824da5d850b9c12d20ec304919c311b508d909e (2003).

65. Storto, A., Dobricic, S., Masina, S. & Di Pietro, P. Assimilating along-track altimetric observations through local hydrostatic adjustment in a global ocean variational assimilation system. *Mon. Weather Rev.* **139**, 738–754 (2011).

66. Hersbach, H. et al. The ERA5 Global Reanalysis. *Q. J. R. Meteorol. Soc.* **146**(730), 1999–2049. https://doi.org/10.1002/qj.3803 (2020).

67. Pawlowicz, R. M_Map: A mapping package for MATLAB", version 1.4m, [Computer software]. www.eoas.ubc.ca/~rich/map.html (2020)

## Acknowledgements

We acknowledge the World Climate Research Programme, which, through its Working Group on Coupled Modelling, coordinated and promoted CMIP. We thank the climate modelling groups for producing and making available their model output, the Earth System Grid Federation (ESGF) for archiving the data and providing access, and the multiple funding agencies who support CMIP and ESGF. GODAS data provided by the NOAA/ OAR/ESRL PSL, Boulder, Colorado, USA, from their Web site at https://psl.noaa.gov/. This study has been conducted using E.U. Copernicus Marine Service Information including the C-GLORS and ORAS5 datasets downloaded from https://resources.marine.copernicus.eu. This work is supported by Australian Research Council DP180101251. AST is supported by the Australian Research Council FT160100495. Code for downloading, processing, analysing and plotting model and reanalysis data is available at https://doi.org/10.5281/zenodo.4616306.

## Author contributions

A.S.G. performed the analysis and the initial paper draft. All authors reviewed the manuscript, provided portions of text and helped in the interpretation of the results.

## Funding

V.R. acknowledge funding support from (i) the SEAMoBB project (Solutions for sEmi-Automated Monitoring of Benthic Biodiversity), funded by European ERA-Net Mar-TERA (id. 145) and managed by ANR (grant no.
ANR-17-MART-0001-02) and (ii) the international OFC project (Ocean Front CHANGE), funded by the Bel-
mont Forum and managed by the ANR (grant no. ANR-20-BFOC-0006-04).

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-88934-w.

**Correspondence** and requests for materials should be addressed to A.S.G.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021