Physical Connectivity Between the NE Atlantic Seamounts

Within the Portuguese Exclusive Economic Zone, the Great Meteor and Madeira-Tore complexes are highly productive areas, which are likely to be classified as marine protected areas (MPAs) due to their ecological vulnerability. This was the main focus of the BIOMETORE project and, framed on it, the aim of the present study was to investigate the physical connectivity between both seamount complexes. Using the HYbrid Coordinate Ocean Model coupled with the Connectivity Modeling System (CMS) (a Lagrangian tool), a series of experiments was conducted in order to determine the influence of the main oceanographic phenomena governing the area in: (i) the origin of the particles that reach each complex, (ii) their capacity to capture and retain incoming particles, and (iii) the physical connectivity between them as well as the intra-connectivity within each seamount system. Due to the geographical location of both groups of seamounts, the Azores Current (AzC) and its associated eddies were identified as the main transport pathways, its influence being stronger at intermediate waters and decreasing with depth. Notwithstanding, the Great Meteor and the Madeira-Tore were mainly affected by the AzC southward and eastward branches, respectively, resulting in a non-significant connectivity between the two groups. Meanwhile, the inter-connectivity between seamounts slightly varied with depth at the Great Meteor complex while increasing at Madeira-Tore. In addition, the Plateau, Irving, and Cruiser (PIC) seamounts from the Great Meteor complex and Gorringe and Coral from the Madeira-Tore complex proved to incorporate the regional connectivity routes. Although containing the three smallest seamounts, Madeira-Tore showed the higher capturing capacity per square kilometer, highlighting the influence of the “sticky water effect.” In the Great Meteor complex, the “seamount effect” seems to be the main phenomenon responsible for the greater retention and self-recruitment abilities of these seamounts. The presented results provide valuable information for the design of a MPA to preserve these vulnerable habitats.


INTRODUCTION
Seamounts are wide and prominent underwater features of volcanic or tectonic origin that occur mainly around mid-ocean ridges, island-arc areas, and above mantle plumes (Etnoyer et al., 2010;Staudigel and Clague, 2010;Wessel et al., 2010). The size of the seamount is an important factor to consider and is usually defined by their height or relief (Kitchingman et al., 2007). According to Pitcher et al. (2007), elevations higher than 1,000 m are defined as large seamounts, whereas small seamounts are within the range of 100-1,000 m. Recent studies have estimated that the total number of large seamounts (>1,000 m) worldwide varies between 25,000 and 140,000, whereas small ones (100-1,000 m) are more abundant, approximately between 125,000 and 25 million (Wessel et al., 2010;Yesson et al., 2011). Nevertheless, seamount abundances can vary considerably, depending on the methodology applied and the definition of a seamount (Wessel et al., 2010). Despite this uncertainty, most authors highlight seamounts as one of the world's major (underwater) biomes (Etnoyer et al., 2010).
Seamounts are characterized by high productivity and food availability, being considered hotspots for several benthic and pelagic communities, although this is not true for all seamounts as their distinct morphological and oceanographic features may alter biodiversity and productivity patterns (Morato et al., 2016). The productivity of seamounts is closely related to their shallowness since the interaction between the topographical and the physical features (e.g., internal tides) enhances the vertical movement (upwelling) of nutrients toward the euphotic zone (Tuerena et al., 2019). If the local biomass production is sustained long enough, it will allow the transfer of energy to higher trophic levels. Otherwise, seamounts must rely on outer inputs of organic material to support resident populations (Rogers, 1994;White et al., 2008). Their communities are dominated by suspension feeders, dense aggregations of demersal and benthopelagic fishes, and elevated abundances of zooplankton and micronekton (Clark et al., 2010). Moreover, seamounts are important habitats for species with commercial value since they can create large aggregates that spawn or feed in the vicinity, thus reducing the fishing effort (Clark et al., 2016). However, previous studies showed that seamount-aggregating fishes have higher intrinsic vulnerability to fishing due to their life history attributes, such as slow development, longer lifetime, later sexual maturation, and lower natural mortality (Morato et al., 2006).
Due to the ecological vulnerability of these habitats, the Portuguese government plans to classify the Great Meteor and the Madeira-Tore complexes as marine protected areas (MPAs) under the European Marine Strategy Framework Directive (Jorge et al., 2016). The status of fisheries at Great Meteor seamounts is still unclear, even though evidences of exploratory fishing were found. The fish species most commonly caught in this area are Macroramphosus scolopax, Capros aper, Trachurus picturatus, and Anthias anthias (Fock et al., 2002). At the Madeira-Tore complex, the shallower summits of the Gorringe Bank are the main attraction for fishing fleets. The wreckfish Polyprion americanus and the European conger, Conger conger, are captured during spring and summer, while in autumn and winter the swordfish Xiphias gladius and the black scabbard fish Aphanopus carbo are the main targeted species (Hermida and Delgado, 2016;Campos et al., 2019). Nevertheless, it seems that these fisheries may have reached their limits of geographic expansion and are already overexploited (Hermida and Delgado, 2016). For this reason, the design of these two new MPAs should help balance the exploitation by fisheries and conserve these seamount habitats (Probert et al., 2007).
For the purpose of MPA design, it is necessary to understand population dynamics and the hydrodynamic processes related to seamounts. During the last years, several studies have been focused on describing seamounts biological communities and assessing patterns of benthic biodiversity and marine biogeography (Fock et al., 2002;Wienberg et al., 2013;Christiansen et al., 2015). Samadi et al. (2006) suggested that the genetic exchange between submarine mountains is more frequent than previously assumed and may be related to the dispersal transport induced by ocean currents. The extent to which populations are connected is also determined by the scale, intensity, direction, and frequency of the dispersal pathways (physical factors) and by the post-settlement mortality, growth, and condition from settlement to successful reproduction (biological factors; Pineda et al., 2007;Etter and Bower, 2015). However, the wide variety of interconnected mechanisms that promote or prevent the movement of individuals and their contribution to the gene pool are still poorly known (Shank, 2010). Recent methodological advances in larval tagging, parentage analysis, and otolith microchemistry shed some light on marine connectivity over relatively short spatial and temporal scales. Nonetheless, these methodologies are very expensive since they require intensive spatial and temporal sampling. On the other hand, biophysical numerical models offer the possibility to track virtual particles over longer spatial and temporal scales. This type of approach will contribute to identify the main physical processes that serve to connect or isolate seamount populations (Shank, 2010) as well as to determine their sourcesink relationships (Holstein et al., 2014).
In the northeast Atlantic Ocean, several studies focused on comprehending the oceanographic phenomena, as well as the transport and retention mechanisms, around several Portuguese small islands and archipelagos (e.g., Caldeira et al., 2002;Barbosa-Aguiar et al., 2011;Sala et al., 2013Sala et al., , 2016Caldeira and Reis, 2017). However, none of them focused on the connectivity between seamount complexes. To better understand how persistent oceanographic features affect the recruitment and the connectivity between Madeira-Tore and Great Meteor seamount populations, our study investigates the physical connectivity patterns in the region. For this purpose, a series of numerical experiments was performed using the Connectivity Modeling System (CMS), an off-line Lagrangian tool, attached to an ocean circulation model (as described in detail in section "Data and Methods"). The results of the different numerical experiments are presented in section "Results, " and the general discussion and main conclusions are in section "Discussion and conclusions."

Study Area and Oceanographic Context
Within the northeast (NE) Atlantic region (Figure 1a), the Great Meteor complex is located in the mid-Atlantic ridge, ∼900 km to the south of the Azores, and it is constituted by a series of extinct volcanoes rising from depths greater than 4,500 m to less than 300 m (Beckmann and Mohn, 2002;Geldmacher et al., 2006). The ancient volcanic cones are aligned in a north-south FIGURE 1 | (a) Representative map of the study area that includes the North Atlantic Subtropical Gyre and its associated currents: the Gulf Stream, the North Atlantic Current, the Portugal Current, the Azores Current, the Canary Current, and the North Equatorial Current. The solid line defines the HYCOM domain. The dotted lines represent the release boundaries for the particle retention experiment. (b) Zoomed-in view of the area that contains the seamounts under study [the dashed line on panel (a)]. The dashed lines represent the release lines for the reverse particle experiment. The black boxes represent the sink boxes for the particle retention experiment. The red boxes represent the connectivity polygons defined for the connectivity study.
orientation. This assembly includes Atlantis, Plateau, Irving, Cruiser, Tyro, Hyères, and the Great Meteor, as well as the smaller Small-Meteor and the Closs bank (Figure 1b). In turn, Madeira-Tore complex, located ∼700 km off the northwestern African coast, is a northeast-aligned submarine ridge in the central east Atlantic that rises from more than 4,000 m in water depth to ∼35 m below the sea surface (Figure 1b). Several hypotheses have been proposed for its origin (Tucholke and Ludwig, 1982;Morgan, 1983;Peirce and Barton, 1991;Royden, 1993;King and Anderson, 1998), but there is a lack of age and geochemical data (Geldmacher et al., 2006). The complex is composed of the Gorringe bank (with two mountain tops: Gettysburg and Ormonde), Seine, Josephine, Lion, and Unicorn seamounts. The present study is mainly focused on the beststudied groups in this area: Great Meteor, Plateau-Irving-Cruiser (PIC), Tyro, Lion-Josephine, Unicorn-Seine, and Gorringe bank. Additionally, Ampere and Coral seamounts were included due to their proximity to the Madeira-Tore complex.
These deep-sea areas exist in a highly variable oceanographic environment characterized by the prevailing ocean currents and the mesoscale phenomena that result from the interaction of passing ocean flows with the seamount bathymetry. Both complexes, Great Meteor and Madeira-Tore, are affected by the Azores Current (AzC), which is one of the branches of the Gulf Stream (GS) that composes the North Atlantic Subtropical Gyre (Figure 1a; Käse and Siedler, 1982;Gould, 1985;Onken, 1992). The AzC splits in one branch that goes into the Gulf of Cádiz (GoC), toward the Strait of Gibraltar, and another that deviates southward, connecting with the Canary Current (CaC; Johnson and Stevens, 2000). The CaC also receives a small contribution from the Portugal Current (PoC; Barton et al., 2001), which spreads in the north-south direction, connecting the North Atlantic Current, the northeastward branch of the GS, and the AzC (Dietrich et al., 1980). This current regime, along with other phenomena such as ocean eddies, internal tides, and variable meteorological forcing, contributes to mesoscale variability and, ultimately, the different connectivity levels occurring between the seamounts (Clark et al., 2010;Barbosa-Aguiar et al., 2011;Caldeira and Reis, 2017). There are two main processes associated to these features: the (i) "island mass effect" (or in this case "seamount effect"; Doty and Oguri, 1956), which is responsible for the enhanced production occurring around small islands/seamounts in comparison to the surrounding waters, and the (ii) "sticky water effect" (Wolanski, 1994), which generates slower mean currents in the surrounding regions of high-density topography such as a coral reef and, on a different scale, a group of seamounts, leading to an increase of larvae recruitment and trapping in the regions of weaker flow (Clark et al., 2010;Lavelle and Mohn, 2010;Sala et al., 2016).

Ocean Model and Lagrangian Tool
In this study, we applied the HYbrid Coordinate Ocean Circulation Model (HYCOM), 1 which is capable of solving the hydrostatic primitive equations in a free surface mode (Bleck, 2002;Chassignet et al., 2007). The US Navy Research Laboratory global solutions (GLBa0.08, experiment 05.8) were used with a horizontal grid resolution of ∼9.25 km (1/12 • ) and 32 vertical layers. For the off-line Lagrangian simulations, data were extracted from a sub-domain encompassing 22 • to 47 • N in latitude and 55 • to 5 • W in longitude (solid line in Figure 1a). Frontiers in Marine Science | www.frontiersin.org A 9-year-solution was considered (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012). This model was previously validated by comparison with Argo (global array of temperature/salinity profiling floats) profiles by Sala et al. (2016), who concluded that the model ensures a good representation of the known water masses and the ocean currents in the region. The Connectivity Modeling System (CMS v2.0; Paris et al., 2013) was used to track passive particle movements. The CMS runs off-line, applying the 3D velocity fields (u, v, w) of the ocean circulation model to each particle using a 4 th order Runge-Kutta numerical discretization method, applied both in space and time (Van Sebille et al., 2012;Paris et al., 2013). To represent the subgrid-scale motion (physical turbulence effects) unresolved by the model, a value of 2 m 2 s −1 for horizontal diffusivity was added (following Sala et al., 2016). No vertical diffusivity FIGURE 3 | Particle density distribution for the reverse particle tracking study at 5 m depth for particles released at the different seamounts: (a) Great Meteor, (b) Lion-Josephine, (c) Unicorn-Seine, (d) Ampere-Coral, and (e) Gorringe Bank. The black bar corresponds to the final destination of the particles. The red boxes represent the nearby seamounts of interest (as identified in Figure 1b).
Frontiers in Marine Science | www.frontiersin.org was considered since these particles traveled in discrete layers while maintaining a constant depth. This Lagrangian tool was especially developed for larval dispersal modeling and to give probability estimates of population connectivity, with multiple modules available for the integration of biotic and abiotic data. Although some modules mimic larval behavior, none was activated in order to provide more accurate representations of the Lagrangian ocean circulation patterns around the Great Meteor and Madeira-Tore complexes. Therefore, the model was applied as a statistical representation of dispersal probabilities rather than as a deterministic study of individual particles and/or larval fates (Ross et al., 2016).

Reverse Particle Tracking Study
This case study was carried out to determine the origin of particles reaching these submarine archipelagos, Great Meteor and Madeira-Tore (also considering the Azorean seamounts), discerning the main mesoscale phenomena that influence the transport dynamics. For each seamount, a line following their longitudinal or meridional axis was considered (depending on their orientation, dashed lines in Figure 1b), with several release points defined at a regular spacing of 0.05 • (∼5.55 km): for Great Meteor 117 release points were considered, Lion-Josephine with 63, Unicorn-Seine with 42, and Ampere-Coral and Gorringe bank with 49 release points, adding up to a total of 320 release points. During this experiment, 10 particles were released every month at four different depth layers (5, 150, 500, and 1,000 m) at each release point, summing up to a total of 192,000 particles released at each layer. These release depths represent the common transport pathways previously identified for the northeast Atlantic Ocean by Sala et al. (2013). Trajectories of individual particles were tracked monthly ("backward" in time) from December 31, 2012 until January 1, 2008, considering that 5 years was enough time for the slowest particles to exit the model domain.

Particle Retention Study
This case study was conducted in order to evaluate the capacity of each seamount to capture and retain incoming particles. For these experiments, the particles were released from several release points regularly separated by 1 • (∼111.11 km) along the model domain boundaries: the North Boundary (NB) considered 39 release points, the West Boundary (WB) with 21, and the South Boundary (SB) with 32 (dotted lines in Figure 1a). At each release point, 10 particles were released every month at the same depth layers as the reverse particle case study. Therefore, in the NB, a total of 23,400 particles were released, in the WB 11,400 were released, and in the SB 19,200, were released. In this case, the trajectories were monitored ("forward" in time) also at a monthly interval from January 1, 2004 until December 31, 2008.
To assess the "capture capacity" and "retention time interval" of each group of seamounts relative to incoming particles, a set of five "sink boxes" was defined: Great Meteor, Lion-Josephine, Unicorn-Seine, Ampere-Coral, and Gorringe bank (Figure 1b). "Capture capacity" was calculated by counting the number of particles that reached each group of seamounts (i.e., each "sink box"), while the "retention time interval" (in days) was calculated as a frequency distribution, considering the number of days that each particle spent inside a "sink box." In this experiment, the journey of a particle did not end once this particle enters a "sink box." Each particle will continue its journey until the end of the 5-year simulation or until it leaves the model domain.

Connectivity Study
This experiment analyzed the degree of connectivity between the different seamounts. For that, particles were released from the inside of 10 different 3D polygons (Figure 1b) at four different depth ranges following Sala et al. (2013): 1-10, 20-200, 300-500, and 600-2,000 m. Those 10 polygons represent the main seamounts that are currently the focus of this study (i.e., Great  Meteor, Plateau-Irving-Cruiser, Tyro, Lion, Josephine, Unicorn, Seine, Ampere, Coral, and Gorringe bank). At each polygon, 200 release points, randomly distributed (latitudinally, longitudinally, and in depth) and fixed in time, were considered. At every point, 10 particles were released at a monthly interval for 9 years (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012) and tracked "forward" in time. A particle's journey ends once it enters a polygon, but they can only settle 160 days after the release event. Connectivity matrices were calculated considering the number of particles released and captured by each polygon, reflecting the probability of a particle released from a given source box (i) to settle in a sink box (j), thus connecting the population i with the population j.
In order to quantify the physical connectivity within the network of seamounts, "source-sink indices" (SS) were calculated as the difference between the total number of particles being exported (K out ) and the total number of particles being imported (K in ) divided by the sum of the two (Holstein et al., 2014): where K represents the successful settlement "in" (i.e., particles imported) or "out" (i.e., particles exported) of a seamount.
Negative index values (SS < 0) indicate that the seamount imports more particles than it exports, acting as a net sink. On the other hand, positive index values (SS > 0) indicate that a seamount exports more particles than it imports, acting as a net source. Meanwhile, a value equal to 0 would indicate that the seamount imports as many particles as it exports. The "self-recruitment" (selfr) and "subsidy-recruitment" (subr) fractions of each seamount were calculated to identify which seamounts rely most on their own recruitment capacity versus the ones that depend more on the connections with its surroundings to maintain the local populations.
The selfr fraction was determined as the proportion of total recruited particles to a seamount (i) that originated from the seamount itself (Andrello et al., 2013): where c is the connection probability, i.e., the fraction of particles originating in the seamount i that ended up in the same location (c(i,i)) or in another seamount j (c(i,j)).
Meanwhile, the subr fraction was obtained as the proportion of total recruited particles that came from the surrounding seamounts (j) (Andrello et al., 2013), calculated as: Thus, if selfr is higher (lower) than 0.5, there are more (less) particles retained by self-recruitment than by subsidyrecruitment.

Sensitivity Tests
Several sensitivity tests were performed in order to examine the limitations of the model. A first test to determine the seasonal and inter-annual variability for the capture and retention capabilities of the seamounts showed a small diversion for the bulk transport pathways. Six commonly used indices were computed to compare the model results and assess its performance: the correlation coefficient, root mean square error, reliability index, average error, average absolute error, and modeling efficiency [for more detailed information, see Stow et al. (2009)]. Notwithstanding, the tests were limited to the most superficial layer where ocean dynamics is more complex, assuming that there are no major changes with increasing depth. For the reverse particle tracking study, all possible particle pathways were recalculated with half of the particles, obtaining very similar dispersal patterns to those obtained with all the particles. For the particle retention and connectivity studies, the capture capability of each group of seamounts was recalculated through the same method, and similar results were obtained ( Table 1). Although tens of thousands of particles were released per depth range, these tests proved that the model achieves statistical convergence using all or even half of the number of particles.
Additional sensitivity tests were performed for the connectivity study. The connectivity matrices with particles that are uniformly and randomly released (with respect to x and y coordinates) were compared, and the results were similar ( Table 1). With matrices in which the release events were set at and Gor -Gorringe. The y-axis represents seamounts acting as sources and the x-axis represents seamounts acting as sinks. The connectivity degree increases from blue to red. (b) Source-sink indices for each seamount. (c) Self-recruitment (selfr, blue color) and subsidy-recruitment (subr, red color) fractions for each seamount. The dashed line corresponds to the reference value of 0.5.
Frontiers in Marine Science | www.frontiersin.org a constant depth and at a random depth (z-coordinate) within the corresponding depth range (i.e., 1-10, 20-200, 300-500, and 600-2,000 m), the outcomes were also comparable ( Table 1). Nevertheless, as this study is not species specific, which are often confined (due to behavior) to a certain depth range, random releases at different depth ranges were considered in order to achieve a generic approach.
Finally, a comparison between the model and surface drifting buoys data was conducted. A total of 772 drifting buoys trajectories were downloaded from the CORIOLIS database, 2 taking into account the region (dashed line in Figure 1a) and the study period (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012). Considering the release point of each drifting buoy (i.e., longitude, latitude, and time) and its travel time (31 days maximum), a CMS experiment was performed. The resulting connectivity matrix was compared with the one obtained from the drifters (see Figure 2), revealing a high correlation between the ocean circulation model and the drifting buoys data (see Table 1).
The sensitivity studies, particularly the comparison with observations, gave the authors confidence to pursue a 2 http://www.coriolis.eu.org representative study of the physical connectivity between the NE Atlantic seamounts, of which the main findings are detailed in the section below.

Reverse Particle Tracking Study
With regard to the Great Meteor complex located at ∼1,600 km from the nearest continental coast, at 5 m depth, the majority of the particles came from the west/northwest side of the domain (Figure 3a), whereas for the most coastal mountains (i.e., Lion-Josephine, Unicorn-Seine, Ampere-Coral, and Gorringe bank), the particles came mainly from the Iberian Peninsula (IP), GoC, and the northwest African coast (Figures 3b-e). A particularly high-density pathway was observed between the south IP and the Gorringe bank (Figure 3e). The trajectories for particles released at 150 and 500 m exhibited similar results (data not shown).
At 1,000 m depth, the particles came mostly from the west and the northeast toward the Great Meteor (Figure 4a). At the remaining seamounts, the particles arrived from the nearby  (Figures 4b,c,e), with a major input from the southwest side (Figures 4b,d).

Particle Retention Capability Study
At Great Meteor, Lion-Josephine, and Unicorn-Seine seamounts, a higher percentage of particles was retained in the most superficial layer (5 m; Table 2). For almost all depths, the main source of particles was the NB, with the exception of the deepest layer where the SB outstood as the main contributor ( Table 2). The captured particles spent a minimum of 30 days and a maximum of 250 days inside a box ( Table 3). On the other hand, at Ampere-Coral and Gorringe bank, more particles were retained at intermediate waters (500 m; Table 2). The main origin of the total trapped particles for the first two depth layers (5 and 150 m) was the NB, whereas for the remaining two (500 and 1,000 m), the southern contribution exceeded the northern contribution ( Table 2). The majority of the retained particles spent a minimum of 40 days and a maximum of 100 days inside one of these boxes ( Table 3).

Connectivity Study
The connectivity matrices calculated for each depth range are presented in Figures 5-8. In order to be dimensionless, each element was normalized by the highest number of recruited particles. Therefore, the scale ranges from 0 to 1, with 0 being the lack of connectivity. The X-axis and the Y-axis represent the different seamounts considered in this study, geographically arranged from southwest to northeast, with the main diagonal providing self-recruitment information. It is also important to note that each row represents a release seamount (source), while each column corresponds to a settlement seamount (sink). In addition, the SS index will determine if a seamount is acting as a sink or as a source of particles. Meanwhile, the selfr and subr fractions will distinguish if a seamount relies more on self-recruitment or on subsidy-recruitment, respectively.

Particles Released Between 1 and 10 m
At this depth range, the intra-seamounts connectivity was stronger in the Great Meteor complex (i.e., Great Meteor, Plateau, Irving, Cruiser, and Tyro) than in the Madeira-Tore complex (i.e., Lion, Josephine, Unicorn, Seine, Ampere, Coral, and Gorringe). A higher self-recruitment was found for the western group of seamounts (i.e., Great Meteor complex; main diagonal, Figure 5a). The Great Meteor and PIC acted as net sinks, with the highest self-recruitment (∼0.37 for both) in comparison to the other seamounts (Figure 5a). These results coincided with the computed SS index < 0 for both mountains, which is indicative of a sink behavior (Figure 5b). For the PIC mountains, the main sources of recruited particles were Tyro, Josephine, Lion, and Great Meteor seamounts. With the exception of Great Meteor, those mountains acted mainly as net sources (SS > 0) and were characterized by low selfr values (Figures 5b,c). Although the PIC seamounts acted as net sinks (SS < 0; Figure 5b), they also presented the highest selfr value, retaining slightly more of their own particles rather than from other seamounts (selfr > 0.5; Figure 5c). Tyro presented the lower value of self-recruitment of the Great Meteor complex (0.31; main diagonal, Figure 5a), which coincides with a positive value of SS index (i.e., acted as source; Figure 5b). Additionally, it revealed a relatively low value of selfr (Figure 5c). On the eastern group of seamounts, where the interaction intra-seamounts was very weak, a stronger connection occurred between the Gorringe bank and Coral patch. For the Coral patch, the highest input of particles came from Gorringe bank, followed by a high self-recruitment pattern, with a low contribution from the remaining sites (Figure 5a). Thus, Gorringe acted as a net source (SS > 0), with a low selfr fraction (Figures 5b,c). As expected, a higher recruitment was observed for seamounts geographically closer to the source (Figure 5a). Moreover, most of the recruited particles were above the main diagonal, evidencing a westward transport.

Particles Released Between 20 and 200 m
At this depth range, the intra-seamounts connectivity from both complexes was intensified, especially in the eastern complex (i.e., Madeira-Tore complex). The exchange of particles between both groups was still present, although most of the recruited particles were below the main diagonal, evidencing an eastward transport (Figure 6a). On the western group, at Great Meteor, the input of particles from PIC slightly exceeded its self-recruitment (Figure 6a). However, Great Meteor acted as a net source (SS > 0), exporting more than it imports (Figure 6b). Moreover, its low selfr fraction proved that more particles from other mountains were retained than its own (Figure 6c). PIC showed the highest self-recruitment value of all mountains (Figure 6a). Besides that, it also received a large contribution from Tyro and Great Meteor (Figure 6b). Thus, PIC showed a selfr fraction value close to 0.5, confirming that the number of particles selfrecruited was similar to the number of particles imported from other seamounts (Figure 6c). Meanwhile, Tyro acted as a source, showing a SS index close to 1, which means that the proportion of exported particles was much larger than the imported ones (Figure 6b). At the eastern group of seamounts (Madeira Tore complex), the self-recruitment of almost all mountains was higher at these depths when compared to the most superficial layer. The main source of particles for Lion was Josephine, with a contribution slightly superior than its own self-recruitment (Figure 6a), whereas for Josephine, the input from Gorringe was FIGURE 9 | Bar plot representing the captured particles per square kilometer, per sink box, by depth. The red stars represent the total number of captured particles per sink box. The sink boxes are arranged in descending order of sizes (i.e., Great Meteor, Lion-Josephine, Gorringe bank, Ampere-Coral, and Unicorn-Seine). Both datasets were normalized so as to be comparable. higher than its own self-recruitment (Figure 6a). Meanwhile, for Unicorn, the greatest input came from Josephine and Lion (0.35 and 0.29, respectively; Figure 6a). For Seine seamount, self-recruitment exceeded inputs from other sources (Figure 6a). Lastly, Gorringe bank acted as the main sink (Figure 6b) for particles released by the surrounding submarine mountains (Figure 6a). Thus, all the eastern group seamounts, with the exception of Gorringe Bank and Coral patch, acted as net sources, albeit the relatively low values (Figure 6b), although at <0.5, Gorringe Bank showed the highest selfr fraction, importing more particles from the surrounding mountains than its own via selfrecruitment (Figure 6c).

Particles Released Between 300 and 500 m
The results obtained at this depth-range were very similar to those presented previously. Again, a strong eastern transport was observed. Higher values of connectivity were found intraseamounts of the western group (i.e., Great Meteor, PIC, and Tyro), while for the eastern group, the strongest connections were established mainly with Gorringe bank (Figure 7a). High self-recruitment was observed (Figure 7a) at the western group, especially in PIC, which continued to act as the main sink (SS < 0) in the Great Meteor complex (Figures 7a,b). PIC also presented similar values of selfr and subr fractions, with the number of self-recruited particles concurrent to the number of imported particles (Figure 7c). At the eastern group of seamounts, Gorringe stood out once again as the main sink for particles released in the surrounding submarine mountains, therefore with a selfr value lower than subr (Figure 7).

Particles Released Between 600 and 2,000 m
Once more the results obtained for this depth range were similar to the previous ranges. However, the intra-seamounts connectivity in both complexes and the eastward connectivity between them were less intense (Figure 8a). All mountains at the western group, particularly PIC stood out, showing a high selfrecruitment (Figure 8a), and acting as the main sink (SS < 0) (Figures 8a,b). Nevertheless, the selfr fraction was higher than 0.5, with the number of self-recruited particles greater than the number of imported particles (Figure 8c). At the eastern group, a higher connectivity was observed for Unicorn, Seine, Ampere, Coral, and Gorringe Bank, with the last two seamounts acting mostly as sinks (Figures 8a,b).

DISCUSSION AND CONCLUSIONS
Recent studies raised concerns on the vulnerability of seamount communities to human impact due to the intrinsic biological factors characteristic of seamount species (e.g., slow growth rates or late maturation; Rogers, 2004;Clark et al., 2010). Of special distress are the consequences of the development of large-scale bottom trawl fishing for sessile fauna (Clark et al., 2010), with coral reefs being particularly vulnerable to damage . In the NE Atlantic seamounts, the level of endemism in coral species is very low (<3%), which means that populations rely on external contributions. In addition, the fauna associated to them is significantly different from that recorded at the same depths on the continental slope (Hall-Spencer et al., 2007). Thus, it is necessary to advance in the protection of these (less abundant) native species. In order to further perceive how these species are sustained, understanding the influence of the oceanographic phenomena in the connectivity between the different seamounts located in this region is of special interest. This study proposes a first approximation to pursue this goal using a numerical modeling approach, focusing on the Great Meteor and Madeira-Tore complexes.
Due to the geographical location of Great Meteor and Madeira-Tore, our results highlighted the main role of the AzC and their associated eddies in the transport, dispersion, and retention of particles in this region, as well as their influence in the connectivity among the different seamounts that comprise those complexes (Figure 1). The Great Meteor complex, located offshore ∼1,500 km away from the African coast, is influenced by both branches of the AzC: the one that travels eastward toward the GoC and the one that veers southward near Madeira. The Madeira-Tore complex, located ∼550 km away from the IP coast, is mainly influenced by the eastward branch. When this eastward branch arrives at the GoC, it produces a front characterized by high mesoscale variability, probably being responsible for the transport of particles from the coastal regions of the GoC toward this complex (Figures 3, 4). Moreover, although on a small scale, the PoC also influences the Madeira-Tore dynamics, transporting particles from the western coast of the IP to the open ocean.
The AzC and its associated eddies are the main regional features injecting particles in our study area; thus, the highest input from the north and west boundaries was expected. Our results showed that the NB was the main source of captured particles by all seamounts at the upper layers (i.e., 5 and 150 m; Table 2). However, the contribution from the WB was only evident at the Great Meteor complex ( Table 2) since it is located closer to this boundary (∼2,000 km) than Madeira-Tore (∼3500 km). An intensified AzC between 150 and 500 m in depth and 32 • and 35 • N (Klein and Siedler, 1989;Comas-Rodríguez et al., 2011;Sala et al., 2016) resulted in the increase of captured particles from the WB by all seamounts. Due to the decrease in intensity of the AzC with depth, a lower influence was observed at 1,000 m depth. Moreover, the released particles had a maximum travel time of 5 years that proved to be insufficient to travel a greater distance (Figure 4d and Table 2). While the NB contribution decreased with depth, the contribution from the SB increased due to the intensification of the north-westward and westward transport mediated by the eddy corridors also associated with the AzC (Sangrà et al., 2009;Barbosa-Aguiar et al., 2011;Comas-Rodríguez et al., 2011;Sala et al., 2016).
The capture capacity (i.e., number of particles that reached each sink box) and the retention time interval (i.e., the number of days that each particle spent inside a sink box) were only partially related to seamount size as they depend also on the proximity of each seamount to the different source boundaries and on the predominant oceanic processes determining their oceanic pathways. As expected, the percentage of captured particles decreased with depth ( Table 2), while the retention time increased for all boxes except at the Gorringe bank, which did not show a clear pattern of variability (Table 3). In the upper mixed layer (i.e., 5 m), Great Meteor (the largest box) withheld 25% of captured particles, followed by Lion-Josephine (twice as small as the Great Meteor) with 11%. Besides the higher capture capabilities, they also had higher mean retention times. Meanwhile, Unicorn-Seine, the smallest sink box (seven times smaller than Great Meteor), showed the third highest capture capacity (6%; Table 2) but with a shorter retention time interval (0-30 days; Table 3). At the other release depths, the variability between the different sink boxes was lower (Tables 2, 3). Notwithstanding, when the capture capacity per square kilometer is calculated for all the sink boxes, Unicorn-Seine stands out with the highest values at all depths (Figure 9). Great Meteor (the largest box) only exceeds in the most superficial layer. From 150 to 1,000 m deep, the smallest sink boxes (Gorringe bank, Ampere-Coral, and Unicorn-Seine, from large to smaller-sized seamounts) showed a capture capacity (per km 2 ) greater than the larger sink boxes (i.e., Great Meteor and Lion-Josephine). This higher capture capacity of the smallest sink boxes could be related to their location and orientation with respect to the incoming AzC (Figure 1). In addition, these aggregated seamount chains are influenced by the "sticky water effect, " which increases their ability to capture and retain incoming particles, including coral eggs and fish larvae (Andutta et al., 2012). However, as their delimited areas are smaller, the particles stay for a shorter time. So, the particles are likely to be retained longer on seamounts occupying larger areas. While the Madeira-Tore complex is mainly influenced by the "sticky water effect, " the Great Meteor complex, due to its orientation to the incoming AzC, must be stirring the incoming current and causing a "seamount effect, " a process analogous to the "island mass effect", which enhances productivity in the lee of oceanic islands (see e.g. Dower and Mackas, 1996). Other physical processes, such as lower rates of flow during neap tides, can still induce particle retention (Andutta et al., 2012). Thus, future studies could consider variable flow rates (e.g., tidal currents).
Considering the results, the connectivity between both complexes and the intra-seamount interaction should be stronger following the direction of the predominant stream/flow. Nevertheless, it should be noted that the complexes are ∼1,000 km apart, and assuming a constant speed for the AzC of 50 cm m −1 (Pingree, 1997), a particle would need at least ∼23 days to travel from one complex to another. At the upper mixed layer (5 m depth), a westward connectivity was observed between Great Meteor and Madeira Tore, with PIC seamounts being the main sink ( Figure 5). This can be related to the westward flow previously reported by Barbosa-Aguiar et al. (2011) between 18 • and 24 • W along 32.5 • N. This counterflow feeds on a southward recirculation east of Madeira, travels westward, and veers to the south at about 24 • W. However, when AzC is intensified (i.e., 20 and 500 m depth), the connectivity seemed to follow an eastward pattern from Great Meteor complex toward Madeira-Tore (Figures 6, 7), whereas in deeper layers (600-2,000 m), where the currents are slower, there was almost no interaction between the two submerged archipelagos (Figure 8).
Regarding the intra-seamount connectivity, this interaction may depend on the submesoscale processes associated with AzC, given the almost perpendicular orientation of the Great Meteor complex (Great Meteor, PIC, and Tyro) relative to the current. In fact, previous observations have shown that the AzC recirculates in westward counterflows located both north and south of the main current, which might enable north-south as well as south-north interactions (Pingree, 1997). This intra-seamount connectivity did not vary substantially with depth. The PIC mountains always stood out with the strongest self-recruitment values, concurrently acting as the main sink for Tyro and Great Meteor, probably because of its larger volume compared with the others (see Figure 1b). The PIC and Great Meteor polygons are closer than PIC and Tyro, allowing a stronger interaction.
Meanwhile, the connectivity within the Madeira-Tore complex should be higher, considering that the seamounts are directly exposed to the AzC flow (Figure 1). At the most superficial layer, the intra-seamount connectivity was very weak. Furthermore, at this depth, all of them acted as net sources. This pattern probably relates with the fact that the particles could only settle 160 days after the release event ( Figure 5). However, intra-seamount connectivity increased with depth, being strongest when the AzC is more intense (i.e., 20 and 500 m; Figures 6, 7). Between 20 and 2000 m depth, Gorringe and Coral seamounts, the two polygons closest to the GoC, acted as prime sinks of this complex. Meanwhile, the other mountains (sited further west) continued to act as sources.
The design of MPAs seems an appropriate way to protect the Great Meteor and Madeira-Tore complexes as vulnerable habitats, helping to balance exploitation and conservation (Probert et al., 2007). According to previous studies, their design should (i) assess the necessary degree of conservation, (ii) define the goals that justify the establishment of new MPAs, (iii) integrate the biological information of threatened species, such as dispersal distances and overall distribution, and (iv) select suitable areas to create MPA networks for biodiversity conservation and fishery management (Jones et al., 2007;Clark et al., 2011). Nonetheless, to maximize the success and benefits of these ocean sanctuaries, knowledge about larval retention and connectivity should also be considered (Jones et al., 2007). In this context, our results provide valuable information for this goal. Only physical connectivity was considered in our study, and the particles behaved as passive tracers without any ecological behavior. Therefore, our results will only be relevant for species with at least one planktonic pelagic phase (i.e., unable to swim against a current). Previous studies highlighted the roles of larval behavior and duration in dispersal and connectivity control (Fox et al., 2016;Lett et al., 2019;Taninaka et al., 2019). For instance, vertical larval migration may allow deep populations from nearby seamounts to repopulate shallow areas impacted by overfishing, which might increase the intraconnectivity within each seamount system (Fox et al., 2016). On the other hand, these swimming abilities might disperse larvae over long distances (Majoris et al., 2019), possibly increasing the inter-connectivity between the two systems. Longer-living larvae are expected to travel greater distances, thus enhancing the existing connectivity pathways (Fox et al., 2016). For this reason, future studies including biological factors (i.e., diel vertical migration, reproduction season, pelagic larval duration, and mortality rates) should target a selection of key species with known and/or well-documented behavior. Apart from focusing on a generic first study, the scientific literature on specific species behavior for Great Meteor and Madeira-Tore seamounts is virtually inexistent. On the other hand, considering that our results showed a very low connectivity between the Great Meteor and Madeira-Tore complexes, the connectivity patterns with the surrounding landscape should continue to be further investigated, for instance, in relation to the Macaronesian Islands and the western Iberian and African coasts. This will allow us to better understand the regional dynamics, including the proposed stepping-stone theory whereby the North Atlantic populations expand to the different islands through seamount chains (e.g., Cho and Shank, 2010).

SUMMARY
• The AzC and its associated eddies were identified as the main transport pathways for particles on both groups of seamounts, especially at intermediate waters.
• The Great Meteor complex was mainly affected by the AzC southward branch, while the Madeira-Tore was mainly influenced by its eastward branch. • Low connectivity was found between the two seamount complexes. • The inter-connectivity between seamounts increased with depth at Madeira-Tore, while minor variations were found at Great Meteor. • The PIC seamounts from the Great Meteor complex and Gorringe and Coral from the Madeira-Tore complex led the main connectivity routes. • Despite their size, the Madeira-Tore seamounts outstood the others with its higher capturing capacity per square kilometer, suggesting a greater influence of the "sticky water effect." • In the Great Meteor complex, the "seamount effect" seems to be responsible for the greater retention and selfrecruitment capabilities. • These results provide a valuable first approach, from the (generic) physical point of view, to support MPA design efforts. However, since no coherent biological traits were found and/or considered, strong conclusions cannot be drawn at this stage. • Future species-specific research including biological factors such as diel vertical migration, reproduction season, pelagic larval duration, and mortality rates are imperative.

DATA AVAILABILITY STATEMENT
The dataset analyzed in this study (HYCOM GLBa0.08, experiment 05.8), from 2004 to 2012, cannot be provided since only the latest data are publicly available. The dataset used for the sensitivity analysis was downloaded from CORIOLIS database (http://www.coriolis.eu.org/Data-Products/Data-Delivery/Dataselection), adapted to our time and spatial domains.

AUTHOR CONTRIBUTIONS
RC proposed the original ideas which lead to the design of the study. ML and IS were involved in the modeling, data collection, and analysis. All authors were involved in the writing and review process and all authors approved it for publication. An earlier version of this work was part of ML bachelor final project.

FUNDING
The BIOMETORE Project (PT02-0018) has partially funded this work. ML and RC were financially supported by the Oceanic Observatory of Madeira Project (M1420-01-0145-FEDER-000001-Observatório Oceânico da Madeira-OOM). IS was supported by a grant of the FPI fellowship program.