Interaction of the Antarctic Circumpolar Current With Seamounts Fuels Moderate Blooms but Vast Foraging Grounds for Multiple Marine Predators

In the Antarctic Circumpolar Current region of the Southern Ocean, the massive phytoplankton blooms stemming from islands support large trophic chains. Contrary to islands, open ocean seamounts appear to sustain blooms of lesser intensity and, consequently, are expected to play a negligible role in the productivity of this area. Here we revisit this assumption by focusing on a region of the Antarctic Circumpolar Current zone which is massively targeted by marine predators, even if no island fertilizes this area. By combining high resolution bathymetric data, Lagrangian analyses of altimetry-derived velocities and chlorophyll a observations derived from BGC-Argo floats and ocean color images, we reveal that the oligotrophic nature of the study region considered in low chlorophyll a climatological maps hides in reality a much more complex environment. Significant (chlorophyll a in excess of 0.6 mg/m3) phytoplankton blooms spread over thousands of kilometers and have bio-optical signatures similar to the ones stemming from island systems. By adopting a Lagrangian approach, we demonstrate that these moderate blooms (i) originate at specific sites where the Antarctic Circumpolar Current interacts with seamounts, and (ii) coincide with foraging areas of five megafauna species. These findings underline the ecological importance of the open ocean subantarctic waters and advocate for a connected vision of future conservation actions along the Antarctic Circumpolar Current.


INTRODUCTION
Marine phytoplankton production drives ecosystem trophic structures (Iverson, 1990;Frederiksen et al., 2006;Boersma et al., 2009;Block et al., 2011) and ecosystem services, ranging from carbon export (Michaels and Silver, 1988;Legendre and Rassoulzadegan, 1995;Heinze et al., 2015) to fisheries production (Chassot et al., 2010). Consequently, understanding the factors underpinning phytoplankton productivity has major implications for process studies. These include the parameterization of nutrient input adopted in biogeochemical models (Aumont and Bopp, 2006;Lancelot et al., 2009;Tagliabue et al., 2014) and the identification of biological hotspots for the management and conservation of marine living resources (Gove et al., 2016;Santora et al., 2017).
In the Atlantic and Indian sectors of the Antarctic Circumpolar Current (ACC) region (Figure 1A), iron limitation generates intense contrasts in the biological productivity (e.g., Martin, 1990;Sullivan et al., 1993;De Baar et al., 1995;Boyd et al., 2000) affecting the distribution of pelagic ecosystems, from primary producers to upper predators (De Broyer et al., 2014). The major source of iron for the shallow ocean is provided by the interaction of the eastward ACC flow with the shallow bathymetry around islands (Sullivan et al., 1993;Thomalla et al., 2011;Ardyna et al., 2017). The island mass effect stimulates dramatic phytoplankton blooms, such as the ones occurring downstream of South Georgia, Crozet, and the Kerguelen Plateau (Sohrin et al., 2000;Blain et al., 2001Blain et al., , 2007Moore and Abbott, 2002;Lancelot et al., 2009;Pollard et al., 2009;Borrione et al., 2014;d'Ovidio et al., 2015; see also Figure 1A). This large phytoplankton biomass is at the basis of the pelagic food web which sustains numerous colonies of predators breeding in the subantarctic islands (De Broyer et al., 2014). Pinnipeds and seabirds intensively forage in the phytoplankton plumes extending in the Antarctic Polar Front (PF) region downstream of South Georgia (Barlow and Croxall, 2002;Staniland et al., 2004;Scheffer et al., 2010;Rogers et al., 2015) and the Kerguelen Plateau (Guinet et al., 2001;Lea et al., 2006;Cotté et al., 2015;O'Toole et al., 2017).
Away from the islands, the interaction of the ACC with large and deep (>2500 m) topographic features, as mid-ocean ridges, is recognized to generate intense upward movements (Sokolov and Rintoul, 2007;Rintoul, 2018). These vertical motions regularly sustain some phytoplankton productivity through the injection of iron into the euphotic layer (Sokolov and Rintoul, 2007). However, with the notable exception of isolated contributions from hydrothermal vents (Ardyna et al., 2019), these vast open ocean regions show a yearly averaged phytoplankton biomass much lower than the island systems. Thereby, their contribution to the primary production of the ACC region is considered comparatively negligible (Graham et al., 2015;Ardyna et al., 2017). Furthermore, these regionally contrasted quantitative assessments of primary production are accompanied by differences in phytoplankton community composition. Islands blooms are dominated by diatoms, while the open ocean low productive (low phytoplankton) regions are generally described by other lower size (nano-, pico-plankton) phytoplankton communities, like for instance Phaeocystis (Detmer and Bathmann, 1997;Gervais et al., 2002;Smetacek et al., 2004;Queguiner, 2013). These High Nutrient Low Chlorophyll (HNLC) communities are less likely to support pelagic ecosystems that include large predators as marine Abbreviations: ACC, Antarctic Circumpolar Current; BGC-Argo, Biogeochemical Argo; CCAMLR, Commission for the Conservation of Antarctic Marine Living Resources; chl a, chlorophyll a; PF, Antarctic Polar Front; SRS, Shona Ridge System; SWIR, South-West Indian R. mammals or seabirds (Cushing, 1989;Moline et al., 2004;Pakhomov et al., 2019).
It is becoming therefore tempting to see the ACC as an (iron limited) HNLC phytoplankton-desert, punctuated by oases in the wake of islands or hydrothermal vents. However, this view is not consistent with megafauna observations. One of the supposedly low productive regions has been recently associated to intense foraging activity of five species of swimming predators breeding at the Prince Edward Islands: Antarctic Fur Seals (Arctocephalus gazella), Southern Elephants Seals (Mirounga leonina), King Penguins (Aptenodytes patagonicus), Macaroni Penguins (Eudyptes chrysolophus), and Southern Rockhopper Penguins (Eudyptes chrysocome chrysocome) (Arthur et al., 2017;Pistorius et al., 2017;Reisinger et al., 2018). Differently from the colonies of the South Georgia and Kerguelen Islands, penguins and seals, after moving to the PF region, travel upstream and explore an apparently oligotrophic area ( Figure 1A). Although some of the animals focus on an area associated to intense phytoplankton blooms sustained by upwelled water masses enriched by hydrothermal iron (downstream of 30 • E; Ardyna et al., 2019), some others travel and forage further upstream, in the supposedly HNLC biologically-desert area. These long trips are recurrent and energy-consuming and require to be balanced by profitable prey captures. But what is the origin of the phytoplankton production sustaining this prey field?
This paper addresses this question, reassessing from an ecological viewpoint the productivity and nutrient delivery of vast open ocean regions in the ACC area. In particular, our work proposes to extend the biogeochemical island mass effect to the numerous seamounts located upstream of the foraging grounds of marine high-level predators ( Figure 1C). In respect to previous studies which assessed the seamounts effect as negligible compared with the island one (e.g., Graham et al., 2015;Ardyna et al., 2017), we propose to investigate it further through two new main datasets. First of all, we use the high resolution seamount inventory provided by Yesson et al. (2011) (an illustrative example is depicted in the bottom panels of Figure 1). Previous studies analyzed coarser bathymetric dataset, thus not allowing to identify most of the seamounts. Secondly, we assume a regional ecological perspective, while previous works had a large-scale biogeochemical focus (Graham et al., 2015;Ardyna et al., 2017). In this aim, we investigate the seamount effect over seasonal time scales and we complement the satellite images with profiling floats' observations, which allow us to infer the total phytoplankton biomass contained in the water column. Indeed, from a circumpolar biogeochemical viewpoint, moderate phytoplankton biomass may appear almost invisible on annual mean climatology when compared to the dramatic production occurring in the wake of the subantarctic islands. However, from an ecological perspective, seasonal secondary blooms may represent a critical input of biomass for local trophic chains. This is especially true when considering the possible effect of these blooms in sustaining foraging grounds of predators who are tied to a specific region for energetic or physiological purposes.
We stress that understanding the factors behind the foraging activity of land-breeding predators is of urgent concern, given the increasing attention that the subantarctic open ocean is receiving within the Commission for the Conservation of Antarctic Marine Living Resources (e.g., Sc-Camlr-XXXVII (2018), Annex 7, para. 3.74; Ccamlr-XXXVII (2018), para. 6.60). A key step to that end is to improve our comprehension of the biophysical processes which sustain biological productivity up to megafauna species. Using a Lagrangian approach (Lehahn et al., 2018), we propose here to test how seamounts complement the island mass effect in fueling pelagic ecosystems within the ACC domain and how they affect the biological production inferred from landbased top predators.

Predators Distribution
Predators distribution is estimated from tracking data of seabirds and seals from the Prince Edward Islands between 2003 and 2014 and collected from research programs supported by the National Research Foundation and the South African Department of Environmental Affairs  1 . The animals' trajectories we use in this work have been analyzed in detail in Reisinger et al. (2018): here we report only the information relevant for our study, which is the largescale geographical extension of their habitat utilization, i.e., the whole smoothed spatial range of predators. In this study we consider only the tracking data of seals and penguins, because they display a more selective habitat than flying seabirds (Reisinger et al., 2018). The Subantarctic Fur Seals (Arctocephalus tropicalis) are not considered because the atsea distribution of this colony is mostly located north of the Subantarctic Front which bound our study region [see Figure 2 in Reisinger et al. (2018)]. Kernel densities of Antarctic Fur Seals (Arctocephalus gazella), Southern Elephants Seals (Mirounga leonina), King Penguins (Aptenodytes patagonicus), Macaroni Penguins (Eudyptes chrysolophus) and Southern Rockhopper Penguins (Eudyptes chrysocome chrysocome) from the Prince Edward Islands were calculated using the "kde2d" function from the R-package "MASS" on the observed tracks (Venables and Ripley, 2002). The habitat utilization for each colony is identified with a threshold of 95% of each predator density (complete kernel utilization distributions are available for each species in Reisinger et al., 2018; Figure 2). This allows us to identify the whole spatial range used by multiple subantarctic predators including different phases and behavior of their trajectories (usually traveling and foraging).

Seamounts
Seamounts are defined as geomorphological features having a vertical elevation of more than 1000 m from the seafloor with limited extent across the top summit region (Menard, 1964;Staudigel et al., 2010). The seamounts' summit location is extracted from Yesson et al. (2011). We retain here only the shallower structures having a maximal water depth of 1000 m, in order to compare them with the ocean upper layer horizontal velocity fields derived from altimetry. In the following, a seamount is defined by a disk centered in the seamount summit and with a 70 km radius. In the same way, the region under the influence of Bouvet Island is defined by a disk centered in the island (3.36 • E; 54.42 • S) with a 70 km radius.

Satellite Observations and Reanalysis Data
Geostrophic velocities, surface chlorophyll a (Chl a) and temperature reanalysis are evaluated between 2000 and 2015 from different products distributed by the European Copernicus Marine Environment Monitoring Service 2 .
Satellite surface Chl a concentration is analyzed with the Global Ocean Color product (product id: OCEANCOLOUR_ GLO_CHL_L4_REP_OBSERVATIONS_009_082) produced by ACRI-ST. The product evaluated has a 1/24 • spatial resolution. Here, we consider monthly-mean composites in order to reduce the cloud coverage and at the same time preserve the mesoscale variability of the Chl a blooms. Similarly to previous studies, the term "bloom" refers to regions of relative elevated Chl a compared to the background concentration (Sokolov and Rintoul, 2007). The statistical significance of the Chl a enhancement is evaluated with the bootstrap method (a detailed description of the method used for the identification of Chl a blooms is provided in the section "Seamount Mass Effect: Detection of the Lagrangian Plume Stemming From the Seamounts"). Chl a values are investigated in spring and summer (from October to February) ahead of the maximum feeding activities of higher trophic levels (Laws, 1985;Moline et al., 2004;De Broyer et al., 2014).
Phytoplankton dominance is evaluated from the PHYSAT observations (Alvain et al., 2005(Alvain et al., , 2008 between 2000 and 2011 from the month of December and January, when we expect a larger diatoms dominance (Smetacek et al., 1990;Alvain et al., 2008). The PHYSAT method has already been validated at the global and regional scale in the Southern Ocean (Alvain et al., 2005(Alvain et al., , 2008(Alvain et al., , 2013; see in particular d 'Ovidio et al., 2010 Supplementary Material for a validation at the latitudes used here). The method identifies probable specific dominant phytoplankton groups based on the analysis of spectral measurements computed by ocean-color sensors. These are nano-eukaryotes, Prochlorococcus, Synechococcus, diatoms, Phaeocystis-like and coccolithophores. Note that a phytoplankton group is defined as dominant when it is a major contributor to the total pigment (Alvain et al., 2005). As demonstrated in Della Penna et al. (2018), a large number of pixels of PHYSAT daily maps are flagged as not classified in this oceanic region, so that composite maps must span very large periods in order to provide meaningful spatial patterns. For this reason, and because we study large scale features, in our analysis we have used multi-year PHYSAT composites, without resolving interannual variability.
Vertical distribution of temperatures is evaluated with the global ocean reanalysis GLORYS2V4 produced by Mercator Ocean for the Global Ocean and Sea Ice Physics (product id: GLOBAL_REANALYSIS_PHY_001_025). This product is performed with NEMOv3.1 ocean model in configuration ORCA025_LIM and has horizontal spatial resolution of 1/4 • with 75 vertical levels. For the analyses we used the product version averaged by month. The water masses associated to the position of the Antarctic Polar Front (PF) are identified as the northern extent of the presence of the 2 • C isotherm minimum in the subsurface (0-500 m depth) (Botnikov, 1963).

Biogeochemical Argo Floats
The vertical profiles of Chl a concentration and backscattering measured from the Biogeochemical Argo (thereafter BGC-Argo) floats have been extracted from the CORIOLIS Global Data Assembly Center 3 . Quality control tests have been achieved on the Chl a and backscattering data (Schmechtig et al., 2014(Schmechtig et al., , 2018. Non-photochemical quenching on daytime profiles are also corrected on the Chl a data following (Xing et al., 2012). In this aim, the maximal Chl a value above the mixed layer depth is extrapolated until the surface and corrected with a specific correction factor. In order to adapt the correction procedure to the study region, the non-photochemical quenching is corrected with the calibration factor of 3.3 used by Ardyna et al. (2019) in the same region. The mixed layer depth is evaluated in each profile as the depth of a density difference of 0.03 kg/m 3 in respect to the density at 10 m depth. The density is computed based on International Thermodynamic Equation of Seawater (Sérazin, 2011).

Seamount Mass Effect: Detection of the Lagrangian Plume Stemming From the Seamounts
In order to identify the effects of seamounts on phytoplanktonic blooms we analyze whether possible water parcels displaying enhanced Chl a have, or have not, been in contact with the seamounts. To achieve this, we compute a backward-in-time, altimetry-based Lagrangian calculation from the bloom periods. This calculation allows us to identify the water parcels that transited over a seamounts area, defined here as a disk of 70 km radius. This value is chosen considering that altimetry data are expected to resolve features on the order of 70 km or larger (Pascual et al., 2006). The aim is to compare the water dispersion pathways computed downstream of the seamounts, in the following referred to as the Lagrangian water plume, to the Chl a plume. This approach has already been validated for the island mass effect, to reproduce the Kerguelen bloom (d 'Ovidio et al., 2015;Sanial et al., 2015) and the Crozet bloom (Sanial et al., 2014) at mesoscale precision.
For each studied period the analysis is conducted in two steps: (1) Identification of the seamounts that are located within deep iron-enriched hydrographic regions. For each seamount we verify if, in the 4 months before the bloom, it was located within the PF. This current core of the ACC is associated to enhanced dissolved iron concentrations (De Baar et al., 1995;Löscher et al., 1997;Measures and Vink, 2001;Mengelt et al., 2001;Viljoen et al., 2018) presumably upwelled from the deep layer (Sokolov and Rintoul, 2007). The PF is identified for the whole period in the vertical temperature reanalysis as described in section "Satellite Observations and Reanalysis Data." We applied the same approach to the water parcels that have been in contact with Bouvet Island, which is also a potential source of iron to the environment (Perissinotto et al., 1992;Croot et al., 2004). The latter contact is estimated by detecting the water masses which transited over the island in the 4 months before the bloom.
(2) Identification of the Lagrangian water plume from the selected seamounts. Each water parcel of the study region is advected backward-in-time from the blooming period for 4 months. The parcels which are found to come from an area where a seamount is present are then given a timestamp with the time, in days, since they left that topographic features. The 4-month time window chosen for the integration assumes that the extension springtime bloom is pre-conditioned by the wintertime advection of iron rich waters (Mongin et al., 2009).
The enhancement of the Chl a concentration within the 100day Lagrangian water plume is evaluated with the bootstrap method in order to obtain its statistical significance. The bootstrap is conducted on sub-sampled Chl a maps, selecting one pixel each five both in longitude and latitude directions. This is done to remove spatial correlations of the order of the typical size of a mesoscale eddy of this region (about 100 km). For the region under the influence of Bouvet, the Lagrangian water plume spreading downstream of the island is removed from the Lagrangian water plume reconstructed with the seamounts. This allows us to distinguish the plume from the seamounts in respect to the one induced by the island.

Chlorophyll a Enhancement Observed From BGC-Argo Downstream of Seamounts
Two BGC-Argo floats (WMO 6901585 and 6901688) sampled the Shona Ridge region. They crossed the seamount chain around 49 • S between the end of December 2014 and January 2015 with a time lag of a few days (Figure 2). A local increase of the depth-integrated biomass over the mixed layer downstream of seamounts (8-10 • E) is recorded by the profiling floats: the Chl a content doubled in the float WMO 6901585 (Figure 2A) and increased by approximately one third in the float WMO 6801688 ( Figure 2B). These observations contrast with the homogeneous low Chl a concentration (about 0.3 mg/m 3 ) observed by satellite in January 2015 in the near surface (background images in Figures 2A,B). When looking at the Chl a content in the water column, neglecting the variability of the mixed layer depth which is deeper than 100 m upstream of the seamounts and shallower than 40 m over them, the intensification of the Chl a signature downstream of the seamounts is evident in both floats (Figures 2C,D). The vertical profiles depict an increase in Chl a and backscattering values in the first 100 m depth downstream of the seamounts (respectively, about 1 mg/m 3 and 3.0 × 10 −3 m −1 ) compared to the water masses upstream (respectively, about 0.6 mg/m 3 and 2.0 × 10 −3 m −1 ). Note that the yellow triangles in Figures 2C,E (respectively, Figures 2D,F) locate the time of the closest position of the BGC Argos trajectories with the yellow seamount in Figure 2A (respectively, Figure 2B). They symbolize the crossing of the seamount chain. The phytoplankton enhancement downstream of seamounts is further investigate with satellite observations which have larger spatial and temporal covers.

Altimetry-Based Lagrangian Reconstruction of the Seamount Mass Effect
We assume here the hypothesis of a seamount mass effect on phytoplankton biomass, in analogy to the well known island mass effect. This hypothesis implies the formation of phytoplanktonic blooms induced by natural iron fertilization spreading downstream of the seamounts. We discuss the possible mechanisms of iron upwelling in section "Discussion".
Two seamount chains are susceptible to lift iron-laden deep waters to the surface layer and consequently enhance the phytoplankton biomass in the studied region: the Shona Ridge System (SRS), which is composed by the Shona Ridge and the Meteor Rise, and the South-West Indian Ridge (SWIR) (see Figures 1B,C for their localization). Yesson et al. (2011) identified 36 seamounts shallower than 1000 meters depth associated with the SRS, with one third of them shallower than 100 m depth, and 26 seamounts shallower than 1000 m associated to the SWIR in the study region (between 0 and 30 • E). When looking at Chl a images not yearly averaged, but over time periods of 1 month during spring or summer, Chl a signatures significantly above the background (around 2-fold) do appear downstream of both the seamount chains (as illustrated in Figures 3A,B, 4A).
In order to further sustain the hypothesis of a seamount mass effect, we study the morphology of the satellite-observed Chl a plumes detected in the vicinity of seamounts. In particular, we analyze whether the shape of these plumes is consistent with a dispersal pattern stemming from a seamount. In order to achieve this, we reproduce the Lagrangian plumes associated with the seamounts during the blooms event using altimetry-derived currents. The analysis consists in detecting the water parcels which transited above the seamounts in a certain time period with a Lagrangian approach. As explained in section (Seamount Mass Effect: Detection of the Lagrangian Plume Stemming From the Seamounts), for each bloom period the Lagrangian plumes are derived only from the subset of seamounts which were in contact with potentially deep iron-enriched hydrographic provinces. Here we consider two regions that have been identified in the literature for displaying higher iron concentration, the water masses downstream of Bouvet Island (Perissinotto et al., 1992;Croot et al., 2004) and the water masses located in the ACC along the jet associated with the PF (De Baar et al., 1995;Löscher et al., 1997;Measures and Vink, 2001;Mengelt et al., 2001;Viljoen et al., 2018). This hypothesis will be discussed in section "Drivers of the Seamount Mass Effect." We illustrate the analysis region by region.

Shona Ridge System
The ACC crosses longitudinally the SRS between 0 and 10 • E. In both February and December 2004, a distinct Chl a plume topping at 0.6 mg/m 3 spreads downstream of the SRS aligned with the PF current core (Figures 3A,B). The monthly-mean composite images of Chl a are compared to the altimetry-derived Lagrangian reconstruction of the water plumes emanating from the seamounts of the SRS which have been in contact with ironenriched water masses (Figures 3C,D). These plumes represent the water transported by the ACC which has transited over shallow seamounts in the 4 months before the blooms. The color of Figures 3C,D represents the number of days since last contact, i.e., a pixel color coded with "30 days" was on the top of a seamount 30 days before the date of the image. The 100-day Lagrangian water plumes (black contours in Figures 3C,D) are overlapped to the Chl a plumes (Figures 3E,F). The Lagrangian analysis reveals that, in both studied periods, the Chl a plume extending downstream of the seamounts matches up to fine details the reconstructed water masses that transited over the seamounts few months previous the bloom.
The good prediction of the extension and shape of Chl a gives us confidence about the role of the selected seamounts. Remarkable examples of the match between ocean-color and altimetry-derived features are the correspondence of meanders and mesoscale features between the two signatures, such as the distinct Chl a patch trapped by a mesoscale eddy around 17 • E 48.5 • S in February 2004 ( Figure 3A). Note that even if the Chl a plume is aligned with the PF current core, particularly in December 2004 (Figure 3D), the complexity of the Chl a plume is captured just with the advection of water masses that surrounded the seamounts. This is the case for instance for the latitudinal extension of the plume close to the SRS, between 3 and 9 • E, in December 2004 (Figures 3D,F) or the extension of the plume between 15 and 30 • E in February 2004 (Figures 3A,C). We remark also that the selected seamounts are different in February and December 2004 and that similarly, the same seamounts may be located in different productivity regimes at the two periods. For instance, in December 2004 the shallow seamounts at 48 • S were not in contact with the current core of the ACC (they are not identified in Figure 3D) and seems to not display enhanced phytoplankton biomass ( Figure 3F). Conversely, the   (Figures 3C,E). This different behavior suggests a temporal variability of the biological response associated to submerged topography features.
During the study periods, February 2004 and December 2004, the mean Chl a concentration in the water masses that transited in the last 100 days over seamounts in contact with the current core is, respectively, 49 and 40% higher than elsewhere in the region depicted in Figure 3. In order to test whether the observed differences in the Chl a content inside and outside of the Lagrangian plumes are significant, a bootstrap sampling is carried out for each studied period over the whole image size. The null hypothesis states that the observed difference in Chl a content inside and outside the Lagrangian plume is not significant. Standard deviation of the bootstrap distributions inside and outside the plumes are, respectively, 1.2 × 10 −3 and 5.0 × 10 −4 mg/m 3 for the year 2003/2004 and 2.1 × 10 −3 and 7.1 × 10 −4 mg/m 3 for the year 2004/2005. In all the study cases the null hypothesis is rejected with >99% confidence. These analyses provide observational evidence of a seamount mass effect associated to the SRS.

South-West Indian Ridge
Bouvet Island (3.36 • E and 54.42 • S, indicated by a square in Figure 4) lies on a volcanic system about 50 km west to the SWIR (Georgen et al., 2001) and represents a known iron input in the south-west of the study region (Croot et al., 2004). The island has been associated with enhanced Chl a concentrations and primary production compared to the surrounding areas. However, its effect seems to be locally restricted (e.g., Perissinotto et al., 1992;Croot et al., 2004, refers also to Figure 1A). Differently to what was found in previous studies, in October 2002 a Chl a plume spreads over almost 2000 km downstream of Bouvet Island ( Figure 4A). However, the altimetry-derived Lagrangian plume stemming from Bouvet Island (green area in Figure 4A) overlaps only partially the Chl a plume observed from the satellite in this period. Indeed, the Chl a plume extended over a much larger area nearby 26 seamounts shallower than 1000 m (Figure 4A).
Similarly, to the analysis of the SRS, here we consider whether the seamounts of the SWIR may explain the Chl a enhancement which is not explained by the island alone. Specifically, we reproduce the altimetry-derived Lagrangian plume from the seamounts that have been in contact with the water masses enriched in iron at Bouvet Island in the few months before the bloom. Indeed, the few seamounts in the north edge of the SWIR that were not in contact with the iron-enriched water masses stemming from the island did not display enhanced Chl a concentrations. The Lagrangian water plume emanating from the selected seamounts allows us to reproduce the Chl a plume integrally (black contour in Figure 4C), differently to the island Lagrangian plume (green contour). As for the SRS, the Chl a plume and the Lagrangian plume stemming from the seamounts are qualitatively very similar. In addition, the mean Chl a concentration inside the Lagrangian plume is two thirds higher than the mean concentration outside (respectively, 0.26 ± 0.098 mg/m 3 inside and 0.16 ± 0.071 mg/m 3 outside).
Therefore, we infer that the region with highest Chl a values corresponds to water masses that transited, in the 100 days before the bloom, over the seamounts (black contours in Figures 4B,C). A bootstrap sampling carried out on the part of the SWIR plume which is not superposed with the island's plume rejects the null hypothesis with >99% confidence. Standard deviation of the bootstrap distributions inside and outside the plumes are, respectively, 2.4 × 10 −3 and 4.8 × 10 −4 mg/m 3.

Comparison of Predators' Distribution, Phytoplankton Biomass and Seamounts Locations in the Southern Agulhas Basin
A 15-year climatology of the December Chl a concentration in the southern Agulhas Basin is shown in Figure 5, which corresponds to the time of Chl a peak in the region (Sokolov and Rintoul, 2007). A 15-year climatology of the Lagrangian water plume issued from the seamounts from the SHS and the SWIR (black triangles), computed as in section "Altimetry-Based Lagrangian Reconstruction of the Seamount Mass Effect, " is overlapped to the Chl a map (black contour). The region covered by the Lagrangian water plume is upstream of the pinnipeds and penguins ranges (white contours). This indicates that the water masses targeted by the predators' colonies have been under the influence of the seamounts in the previous months (the Lagrangian water plume is represented with a 100-day threshold, as in the analyses of section "Altimetry-Based Lagrangian Reconstruction of the Seamount Mass Effect"). The similarity between the climatological Chl a and Lagrangian plumes is less evident compared to the monthly analyses previously considered. This is primarily due to the smoothed signal in the Chl a climatology. However, where the mean pigment distribution is more contrasted, as in the northern edge of the Lagrangian plume, the two signals appear to have qualitatively similar patterns.
In order to investigate the signal of the seamounts in the climatological phytoplankton biomass we compute the mean Chl a concentration in a 2 • latitude narrow band centered on the current core of the ACC which we identify by the PF, between 2000 and 2015 ( Figure 5B). Figure 5B shows that the highest values (mean Chl a about 0.7 mg/m 3 ) are observed upstream of the study region, between 20 and 5 • W and decrease eastward. Intense Chl a concentrations extend downstream of South Georgia (36 • W, 54 • S) over wide distances due to the intense transport (e.g., Korb et al., 2004;Sokolov and Rintoul, 2007;De Jong et al., 2012;Borrione et al., 2014;Graham et al., 2015) and to the regular iron inputs linked to flow-topography interactions (Sokolov and Rintoul, 2007). The averaged Chl a level do not significantly increase in correspondence of the interaction between the ACC and the shallow seamounts of the SRS (Figure 5B). However, we note the high standard deviation of the signal between 0 and 15 • E compared to the downstream signal (respectively, σ ≈ 0.2 mg/m 3 and σ ≈ and 0.1 mg/m 3 ) indicating an important interannual variability. Compared to the SHS, the Ob and Lena seamounts close to 50 • E display a lower mean Chl a concentration ( µ ≈ 0.1 mg/m 3 ) and a lower standard deviation ( σ ≈ 0.1 mg/m 3 ), indicating a lower phytoplankton biomass and a low temporal variability.

Phytoplankton Composition: Comparison of the Seamount Region, Kerguelen Bloom and Open Ocean Region
The possibility for primary production of sustaining the large trophic chain up to megafauna, is known to largely depend not only on its abundance but also on the phytoplankton composition (e.g., Moline et al., 2004). We thus analyze our study region with maps produced with the PHYSAT algorithm (Alvain et al., 2005(Alvain et al., , 2008. The latter cluster together the ocean color spectral signal into classes associated to typical phytoplanktonic types which are statistically dominant for those bio-optical anomalies. PHYSAT patterns downstream of the SRS and SWIR seamounts (0-30 • E, 50-53 • S in Figure 6) mainly contain the same classes found in island-induced blooms (70-100 • E, 50-53 • S in Figure 6) and differ from the type classes found in the Australian-Antarctic Basin (i.e., downstream of the Kerguelen bloom about 100-130 • E) ( Figure 6A). According to the PHYSAT observations, during December and January between 2000 and 2011, the seamount region is dominated more than 55% of the time by diatoms (σ = 30%) and less than 1% of the time by Phaeocystis (σ = 5%). By comparison, the Kerguelen bloom (and respectively, the offshore region downstream of this bloom) is dominated 70% (respectively 16%) of the time by diatoms and 3% (respectively 80%) of the time by the small Phaeocystis (Figures 6B-D).

DISCUSSION
Enhanced primary production has been episodically reported in the past over some deep (>2000 m depth) submerged topographic features of the Southern Ocean, such as the Northwest Georgia Rise (Meredith et al., 2003) and the Maud Rise in the Weddell Gyre (Van Bennekom et al., 1988). At a circumpolar scale, Sokolov and Rintoul (2007) observed that blooms are initiated where the ACC crosses most of the deep bathymetric features and they persist for long distances. However, the biogeochemical focus of subsequent works has led to the conclusion that Chl a enhancements downstream of submerged features are negligible because on average they are much less intense than the ones in the wake of the subantarctic islands (e.g., Graham et al., 2015;Ardyna et al., 2017). These ideas have reinforced the view of the region in between the subantarctic islands as a HNLC phytoplankton-desert area. Here we revisit this apparently low chlorophyll region investigating its ecological implications, in particular for land breeding marine predators for which the possibility of accessing a large islands' bloom far away from their colony may be associated to excessive energetic locomotion costs. This is the case for instance with several species of penguins and pinnipeds breeding at the Prince Edward Islands. A very limited island mass effect is observed around and downstream Prince Edward Islands because of its very narrow plateau (Perissinotto et al., 2000). Competition for food is thus expected to be high in this small productive area. Animals colonizing this island appear indeed to prefer the lower productive but extended inter-island waters than other targets, such as the small productive area around Prince Edward Islands (where competition may be high) or the larger blooms further away, costly to reach.
Our study targets two topographic systems of the southern Agulhas Basin: the Shona Ridge System and part of the South-West Indian Ridge, between 0 • and 30 • E. By focusing on shorter time scales and on higher bathymetric resolution than former studies, our results identify new phytoplankton blooms associated to the water masses that transited over these shallow topographic highs and provide a detailed description of them. Firstly, an increase in phytoplankton biomass is observed by BGC-Argo floats in the wake of the seamounts during the summer period, between the middle of January and the middle of February. This late Chl a enhancement suggests a possible positive biological effect linked to the presence of the topographic highs which is not simply linked to the seasonal bloom phenology (Sallée et al., 2015;Ardyna et al., 2017). Due to the important limitation of the primary production by iron in this region (e.g., Martin, 1990;Sullivan et al., 1993;De Baar et al., 1995), a relative enhancement of the Chl a values should be related to iron inputs in the surface waters. Secondly, the seamount effect is further investigated with satellite observations which have larger spatial and temporal covers. By identifying isolated Chl a plumes on monthly-mean maps and comparing them to the Lagrangian water plumes stemming from the topographic features we demonstrate the widespread influence of seamounts on the surface seasonal phytoplankton distribution. Therefore, the advection of water parcels from shallow topography can be considered as a powerful tool for identifying the effect of isolated submerged topography, similarly, to the analyses computed on islands (e.g., Mongin et al., 2009;Sanial et al., 2014Sanial et al., , 2015d'Ovidio et al., 2015;Ardyna et al., 2017). By comparison, Graham et al. (2015) investigated the spatial variability of the annual mean Chl a at the regional scale in the study area and concluded that, at this time scale, the phytoplankton biomass is not visibly enhanced over or downstream of the seamounts. However, we note that most of seamounts identified in our high resolution bathymetric data were not identified in this previous analysis (see for instance their Figure 5A), which could also influence their conclusions. Similarly, the Lagrangian pathways computed at the circumpolar scale using a bathymetric threshold also do not identify Chl a plumes in the southern Agulhas Basin, because no shallow topography was detected there with these data (Ardyna et al., 2017). These different outcomes highlight the need to consider both high resolution bathymetric data and short time scales for investigating a potential biological effect of seamounts.
In the following we firstly discuss the possible mechanisms at the origin of the seamount mass effect observed in the present study, the limits of the methods employed, and their applicability to other study cases. Secondly, we explore the ecological role of seamounts in the southern Agulhas Basin, from primary production to predators' distribution, as well as some implications of these new insights for conservation activities.

Drivers of the Seamount Mass Effect
What is the origin of the blooms observed in the wake of seamounts in the southern Agulhas Basin? In most of the previous studies, the Chl a anomalies observed over the seamounts and in the surrounding waters were linked to the shallow seamount depths. However, the spatial and temporal variability of biological effects were rarely understood due to the complex nature of the physical processes involved there (e.g., Venrick, 1991;Pitcher et al., 2007;Pauly et al., 2017). The process presented here is framed by the limiting effect of iron on primary production (Martin, 1990) and the preponderant structuring effect of the horizontal transport on biochemical properties and ecosystems in the ACC region. The two necessary conditions identified for the fertilization of water masses in the open ocean southern Agulhas Basin are thus: (i) the presence of hydrographic regions associated, according to the literature, to potential enhanced dissolved iron content at intermediate depth; and (ii) upwelling dynamics due to the interaction of seamounts shallower than 1000 m depth with an energetic branch of the ACC.
In terms of iron content, two possibly enriched regions are identified in the southern Agulhas Basin. The first one is composed by the water masses that have transited over Bouvet Island, which is a known iron source (Perissinotto et al., 1992;Croot et al., 2004). The second one, is composed by the water masses located in the ACC, along the jet associated with the PF. Numerous in situ observations have recorded enhanced iron concentrations in the surface waters close to this ACC front, both in the Atlantic sector (De Baar et al., 1995;Löscher et al., 1997;Viljoen et al., 2018) and elsewhere Mengelt et al., 2001). In the deeper ocean, iron measurements are less frequent (e.g., Tagliabue et al., 2012). Nevertheless, enhanced iron concentrations appear to occur in the Circumpolar Deep Water close to the PF region in correspondence of the Mid-Atlantic Ridge (Löscher et al., 1997). Furthermore, the analysis of the global GEOTRACES database indicates a sharp ferricline in the Atlantic ACC region reflecting the signature of subsurface lateral transfer of dissolved iron from the subantarctic islands (Tagliabue et al., 2014). Indeed numerous studies have concluded that the intense velocities of the PF current core result in a strong transport of nutrients downstream of South Georgia (Korb et al., 2004;Venables and Meredith, 2009;Venables et al., 2012;Borrione et al., 2014) and in a regular iron supply to the surface water due to the ACC-bathymetry interactions (Sokolov and Rintoul, 2007).
In terms of upward dynamics, the interaction of the ACC with sharp topographic features is well known to generate upwelling hotspots carrying deep waters to the surface (Viglione and Thompson, 2016;Tamsitt et al., 2017) and enhancing surface Chl a concentrations (Sokolov and Rintoul, 2007). Three different mechanisms have been identified in the past for explaining upwelling in the presence of seamounts. The first mechanism is eddy fluxes which are intensified when the ACC encounters a rough topography due to an increase in baroclinicity. The eddy activity enhanced in this way at seamounts can vigorously drive upward northern-sourced deep water along steeply sloping isopycnals in the PF region (Tamsitt et al., 2017;Rintoul, 2018). The enhanced eddy activity at seamounts is also associated to intense flow shear, which in turn preconditions an ageostrophic submesoscale circulation characterized by deep reaching vertical motion (Lévy et al., 2018;Siegelman et al., 2019). A second possible mechanism for upwelling is the bottom pressure torque created when the ACC impinges over bathymetric obstacles. This process induces upwelling upstream of a seamount and downwelling downstream of it (Hughes, 2005;Sokolov and Rintoul, 2007;Thompson and Naveira Garabato, 2014;Rintoul, 2018). This situation has been described in the Atlantic sector of the ACC region, at the South and North Scotia Ridges, the South Orkney Islands and the South-West Indian Ridge (Sokolov and Rintoul, 2007). Finally, a third mechanism which is able to connect deep water to the surface layer is the generation of internal lee waves at seamounts. These waves propagate upward and break, enhancing diapycnal mixing (Watson et al., 2013) and fluorescence signature over the seamounts (Brink, 1995;Kunze and Sanford, 1997;Eriksen, 1998). This has been observed for instance at the seamounts of the SWIR, at the northeast of the southern Agulhas Basin (Read and Pollard, 2017).
The interaction of a strong current with a topographic sharp relief is also known to produce Taylor columns, that is, semipermanent anticyclones on the top of the topographic obstacle which can trap water for periods of months or longer (Taylor, 1923;Meredith et al., 2015). These features can have a special ecological interest because they may retain over the seamounts aging trophic webs, thus providing predictable preying hotspots (e.g., Pitcher et al., 2007;Morato et al., 2010;Cascão et al., 2017). Such circulations, more intense at high latitudes due to the enhanced Coriolis force (Owens and Hogg, 1980;Chapman and Haidvogel, 1992), have been observed in the Southern Ocean at the Maud Rise (Van Bennekom et al., 1988) and conditions conducive to their formation have been found at the Georgia Rise, South Scotia Ridge and at the Prince Edward plateau (Perissinotto and Rae, 1990;Meredith et al., 2003Meredith et al., , 2015. All together, the mechanisms here proposed for the Atlantic ACC region are schematized in Figure 7. The iron content increases in the surface layer at the South Georgia Island and sinks to the subsurface along the flow direction, with regular upward movements due to the interactions between the ACC and the topography. The Chl a blooms identified in the study cases (Figures 3, 4) and the important temporal variability of the climatological surface Chl a content close to the seamounts (Figure 5) are the surface expression of this topography-ACC interactions, usually hidden in the mean climatological Chl a signal when compared with subantarctic island-blooms (e.g., Graham et al., 2015;Ardyna et al., 2017). By highlighting the role of iron content of the upstream water masses, the mechanism here suggested could explain why seamounts seem to display reenhanced phytoplankton biomass when they are close to rich systems, such as coastal margins (Van Bennekom et al., 1988;Perissinotto and Rae, 1990;Meredith et al., 2003;Graham et al., 2015;Oliveira et al., 2016;Lemos et al., 2018), or whether they are located in productive frontal zones (Djurhuus et al., 2017). On the other hand (and following Sokolov and Rintoul, 2007) our results identify processes extending at much larger spatial scales compared with the local phytoplankton enhancements described in previous studies which have focused on seamounts in other regions (e.g., Oliveira et al., 2016;Lemos et al., 2018). This larger extent is due to intense advection by the ACC, in analogy to the spread of non-consumed iron from the subantarctic islands in wintertime which preconditions the size of the islands' blooming area (Mongin et al., 2009;d'Ovidio et al., 2015;Graham et al., 2015).
Our approach does not explain all the patterns of primary production in the ACC region. The ACC core also interacts with a group of shallow seamounts (<1000 m depth) between 10 • W and 6 • E ( Figure 5C) which corresponds to climatological Chl a concentrations topping at 0.7 mg/m 3 (Figure 5B). Differently from the analysis of the SHS and SWIR (black contour in Figure 5A) the Lagrangian plume spreading from the seamounts close to 10 • W only partially explain the Chl a plume in this region (result not shown). Most of Chl a seems to be aggregated close to the current core, which, constrained by the Mid-Atlantic Ridge, shifts northward close to 10 • W (Boehme et al., 2008). The specific regional circulation catalyzed by the topography may explain the Chl a distribution. On the other hand, our analysis leaves open the question concerning why seamounts of other regions do not seem to be associated to intense Chl a plumes, such as the shallow Ob and Lena seamounts (∼50 • E in Figure 5). This could be explained by the mechanism schematized in Figure 7. The deep iron, primarily upwelled to the intermediate layer at the South Georgia Island (Korb et al., 2004;Borrione et al., 2014), is assumed to slowly decay along the flow direction from west to east. Following this hypothesis, it seems plausible that the seamount effect may be less intense for topography features located further downstream of the South Georgia Island (∼35 • W), as the Ob and Lena seamounts (∼50 • E), compared to the closer ones (as the SRS and SWIR). We note that the mechanism here proposed consider primarily external iron inputs, i.e., the upstream iron transport. However, some local processes have also been identified as sources of surface nutrient supply at seamounts. For instance, the interaction of the porewater with sandy sediment releases reduced metal species to the water column (Huettel et al., 1998) and leads to an enhanced erosion of the sediments due to the physical and chemical actions of potential benthic communities (Rhoads and Boyer, 1982). For all these reasons, more detailed studies should be done in order to unambiguously identify the mechanisms at the origin of the Chl a enhancements observed in this study. FIGURE 7 | Schematic summarizing of phytoplankton productivity in the Atlantic ACC region. The interaction of the ACC with topographic highs generates upward movements and facilitates the injection of iron in the euphotic layer (symbolized by the burgundy-colored jets). Lightening of the phytoplankton blooms (green) is an indication of the decreasing of the Chl a concentration. Subantarctic pinnipeds and penguins from Prince Edward Islands (symbolized by the gray arrow) are foraging on this complex biological system.

Ecological Significance of the Seamount Mass Effect in the Southern Agulhas Basin and Implications for Conservation
The surface Chl a concentrations observed in satellite images downstream of seamounts are minor compared to the ones occurring downstream of islands. For instance, in the Crozet bloom Chl a concentrations exceed 2 mg/m 3 (Bakker et al., 2007), four times the ones observed in our analysis. However, we stress that, to assess the total phytoplankton biomass of a given region, one must integrate Chl a concentration over the mixed layer (ML) (e.g., Ardyna et al., 2019). The ML in the seamount region is at least between 10 and 30% deeper of the ML in the Kerguelen bloom region during January and February and between 30 and 50% deeper during December (Supplementary Figure S1). When compared to the Crozet and the South Georgia blooming regions the seamount ML increases up to 60-100% (Supplementary Figure S1). Therefore, the phytoplankton biomass integrated over the ML may present less intense contrasts between the seamount open ocean and the island bloom regions than the surface signals. Furthermore, we note that a deeper ML may also induce a lower Chl a signal in the satellite images, intensifying the observed surface differences between the two regions.
The analysis of the phytoplankton composition of the seamount plumes provides further insights about the relevance of this biological production from an ecological perspective. The phytoplankton community identified in the Atlantic ACC shows a bio-optical signature similar to the ones found in diatom-dominated island blooms, such as the region downstream of Kerguelen, and different from the nearby oligotrophic waters, characterized by Phaeocystis-types (Figure 6). This nanophytoplankton is known to be rapidly consumed by the microbial food web (Detmer and Bathmann, 1997;Pakhomov et al., 2019) and therefore its occurrence is not associated with profitable foraging grounds for top predators (Cushing, 1989;Moline et al., 2004;Pakhomov et al., 2019). Conversely, diatomdominated blooms are more likely to support large trophic webs which lead to fish and crustaceans and higher predators (Cushing, 1989;Moline et al., 2004;Pakhomov et al., 2019). Presence of a large and rich (in terms of diversity) megafauna in our study region is confirmed by biologging data (Arthur et al., 2017;Pistorius et al., 2017;Reisinger et al., 2018). The similarity between the seamounts region to the island blooms is supported also by in situ measurements showing diatom enhanced concentrations at the Atlantic PF region in the vicinity of the Islas Orcadas Rises and Mid-Atlantic Ridge (Laubscher et al., 1993;Froneman et al., 1995;Bathmann et al., 1997;Detmer and Bathmann, 1997;Smetacek et al., 1997). On the other hand, the contrasted phytoplankton community between the seamount region and the classic HNLC areas (as the one between 110 and 130 • E in Figure 6) would emphasize the interest of regions downstream of topographic highs. However, the lack of in situ phytoplankton observations in this low chlorophyll region prevent us to conclude on the potential zonal differences in the circumpolar phytoplankton communities of the PF region. Observation of mid trophic levels would also be valuable for linking the primary producers to the higher predators and verifying these hypotheses. Flynn and Williams (2012) showed that the topography-flow interaction over the Macquarie Ridge (160 • E) seems to induce an increased biomass of the lanternfishes, a key prey for pinnipeds and seals (e.g., De Broyer et al., 2014). Further analyses of intermediate trophic levels should be done for explaining the biological productivity inferred from land-based top predators and assessing how each colony specifically relies on seamounts production. For instance, the distribution of surface and low trophic level foragers, such as Southern Rockhopper Penguin and Macaroni Penguin (Brown and Klages, 1987;Adams and Brown, 1989;Wilson et al., 1997;Crawford et al., 2003;Whitehead et al., 2016), might be more closely related to the surface productivity than the distribution of mesopelagic foragers.
Our analyses advocate for a connected vision of the biological productivity of the ACC region. Following the ACC, the high Chl a concentration decreases from west to east downstream of the South Georgia Island around the ACC core, but it is reenhanced after the ACC interacts with submerged topographic features, in particular the SRS (as schematized in Figure 7). A more comprehensive view may have consequences especially in conservation policies. Currently, conservation planning in the ACC region is dominated by national programs and consequently displays important discontinuities constrained by the jurisdictional framework (e.g., exclusive economic zones). National Marine Protected Areas have been established around Heard and Mcdonald Islands, Kerguelen Islands, Crozet Islands, and Prince Edward Islands in order to conserve the rich and diverse ecosystems hosted in these oceanic waters. By comparison, the high seas of the southern Agulhas Basin, which are already managed for the fishing activities (Supplementary Figure S2), are currently not considered in conservation plans. However, our observations, along with the numerous foraging hotspots of numerous predators occurring in this latitudinal band, both from Bouvet Island and other subantarctic colonies (Biuw et al., 2010;Hindell et al., 2016;Arthur et al., 2017), highlight the need to consider the subantarctic range as a continuum. Current attempts made under the Commission for the Conservation of Antarctic Marine Living Resources (CCAMLR) aim at characterizing the ecosystems structure and function in the subantarctic region and eventually identifying priority areas for conservation. Therefore, efforts shall address the continuity of this oceanic region and the biophysical processes across the high seas and the multiple economic exclusive zones.

SUMMARY AND CONCLUSION
In the Antarctic Circumpolar Current (ACC) region, island blooms stretch for hundreds of kilometers along the current. Island related production plays a primary role in the biogeochemical cycles of the Southern Ocean and sustains large trophic chains up to megafauna species. Conversely, submerged topography is also known to stimulate primary production in its interaction with the ACC, but this productivity is comparatively much less intense (e.g., Sokolov and Rintoul, 2007;Graham et al., 2015;Ardyna et al., 2017). As a consequence, its contribution has been considered negligible at the basin scale, and the ACC waters are often described as a high nutrient low phytoplankton desert, punctuated by massive blooms in the wake of the islands. Conservation efforts are in line with this description, with subantarctic marine protected areas only occurring in the economic exclusive zones around the islands.
Here we revisited the role of seamounts in this framework, focusing on a region of the Indo-Atlantic ACC zone. This area is massively targeted by marine predators, even if no islands fertilize it. By combining high resolution bathymetric data, multi-satellite data, and BGC-Argo observations, we found that the oligotrophic nature of the study region apparent in Chl a climatological maps hides a much more complex situation. Important production (in excess of 0.6 mg/m 3 , Figures 2, 3) occurs there and spreads over huge distances, comparable to the extension of the islands' blooms.
By applying a Lagrangian analysis to these lesser blooms we connected unambiguously these blooms to specific seamount systems. Our Lagrangian scheme was then depicted to track at mesoscale precision the thousands of kilometers long water plume influenced by each individual seamount system (Figures 3, 4).
In order to estimate the type of phytoplankton present in these blooms, we used the PHYSAT product, which provides maps of putative planktonic functional types. The PHYSAT algorithm revealed that the signatures of seamounts' bloom are similar to the ones found in the islands' wake, and not nearby less productive waters (Figure 6). This is consistent with the presence in both systems of foraging grounds for megafauna species.
The ecological role of seamounts in sustaining large trophic chains was finally confirmed by comparing climatologies of Chl a to the predators' distributions ( Figure 5). Our results allow us to conclude that, although temporally intermittent, less predictable, and less effective than the ones associated with the island mass effect, the blooms downstream of seamounts cannot be neglected. Our findings advocate for a connected vision of the ecosystems along the ACC and may support future conservation actions by providing the tools to track the areas under the seamount effect.

DATA AVAILABILITY STATEMENT
All the data used in this research are freely available and may be downloaded through the links detailed in the section "Materials and Methods."

ETHICS STATEMENT
Ethical review and approval was not required for the animal study because this study analyzed observations of vertebrate animals previously published by another study reviewed and approved by an animal ethics committee (doi: 10.1111/ddi.12702).