# First description of deep benthic habitats and communities of oceanic islands and seamounts of the Nazca Desventuradas Marine Park, Chile

Jan Tapia-Guerra, Ariadna Mecho, Erin Easton, María de Los Ángeles Gallardo, Matthias Gorny, Javier Sellanes

► To cite this version:

Jan Tapia-Guerra, Ariadna Mecho, Erin Easton, María de Los Ángeles Gallardo, Matthias Gorny, et al. First description of deep benthic habitats and communities of oceanic islands and seamounts of the Nazca Desventuradas Marine Park, Chile. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-85516-8. [https://doi.org/10.1038/s41598-021-85516-8](https://doi.org/10.1038/s41598-021-85516-8)

HAL Id: hal-03183901

https://hal.science/hal-03183901

Submitted on 4 Jun 2021

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

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

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

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

Check for updates

OPEN

# First description of deep benthic habitats and communities of oceanic islands and seamounts of the Nazca Desventuradas Marine Park, Chile

Jan M. Tapia-Guerra<sup>1,2,3</sup>, Ariadna Mecho<sup>3,4</sup>, Erin E. Easton<sup>3,5</sup>, María de los Ángeles Gallardo<sup>2,3</sup>, Matthias Gorny<sup>6</sup> & Javier Sellanes<sup>2,3✉</sup>

Seamounts and oceanic islands of the Chilean Exclusive Economic Zone at the intersection of the Nazca and Salas y Gómez ridges lie within one of the least explored areas in the world. The sparse information available, mainly for seamounts outside Chilean jurisdiction and shallow-water fauna of the Desventuradas Islands, suggests that the area is a hotspot of endemism. This apparent uniqueness of the fauna motivated the creation of the large Nazca-Desventuradas Marine Park (NDMP, ~300,000 km²) around the small islands San Felix and San Ambrosio in 2015. We report for the first time a detailed description of benthic microhabitats (i.e., centimeter to meter scale), macrohabitats (i.e., meter to kilometer-scale) and associated megafauna within the NDMP. Descriptions were based on analysis of fauna collected by trawling and ROV video observations from ~50 to 370 m depth. Rocky, coarse sand and silty sediment bottom habitats were observed at island slopes. In contrast, rocky and coarse sandy bottom habitats with a predominance of rhodoliths, thanatocoenosis, and other biogenic components were observed at seamounts. Mobile fauna and predators dominated the oceanic islands and nearby seamounts, whereas seamounts farther from the islands were dominated by sessile and hemisessile fauna that were mainly suspension and deposit feeders. Based on the register of 118 taxonomic units, our results provide an expanded and updated baseline for the benthic biodiversity of NDMP habitats, which seemed pristine, without evidence of trawling or anthropogenic debris.

Seamounts are topographic elevations that rise more than 1000 m above the seafloor<fnref n="1" />, and in some cases, they reach the euphotic zone or the surface (i.e., oceanic islands)². These ecosystems are considered vulnerable marine habitats because they are physically fragile or inherently rare³,⁴. Studied seamounts and oceanic islands are generally considered to play important roles by being highly productive locally, acting as regional centers of speciation or stepping-stones for dispersion. They comprise a diversity of substrates and thus habitats that support and provide refuge and feeding areas, not only for the benthic fauna but also for the associated pelagic and overly-ing surface water animals (e.g., birds, and marine mammals)⁵⁻⁷. Although seamounts are generally considered to play these important roles, each seamount experiences unique environmental conditions such as current characteristics and water-column environmental conditions (i.e., temperature, salinity and pressure) that differ with depth and current dynamics⁶. They also have unique characteristics such as height and morphology that result in unique physical dynamics that can result in increased nutrient availability relative to the surrounding

<footnote n="1">¹Programa de Magister en Ciencias del Mar Mención Recursos Costeros, Facultad de Ciencias del Mar, Universidad Católica del Norte, Larrondo 1281, Coquimbo, Chile. ²Sala de Colecciones Biológicas, Facultad de Ciencias del Mar, Universidad Católica del Norte, Larrondo 1281, Coquimbo, Chile. ³Millennium Nucleus for Ecology and Sustainable Management of Oceanic Islands (ESMOI), Departamento de Biología Marina, Facultad de Ciencias del Mar, Universidad Católica del Norte, Larrondo 1281, Coquimbo, Chile. ⁴UMR8212 Laboratoire des Sciences du Climat et de l'Environnement (LSCE), Paris, France. ⁵School of Earth, Environmental, and Marine Sciences, University of Texas Rio Grande Valley, 33363 Marine Lab Dr., South Padre Island, TX 78597, USA. ⁶Oceana Inc. Chile, Santiago, Chile. ✉email: sellanes@ucn.cl</footnote>
area (reviewed in$^6$). Such variability in environmental conditions results in variable ecological characteristics. Despite the above-mentioned importance of seamounts, only 4% of the ~40,000 reported have been sampled for scientific purposes$^{8-10}$. The available faunistic information for seamounts is limited (especially for depths > 300 m) and, for unstudied seamounts, usually comes from commercial fishing activities$^{11}$.

The Pacific Ocean basin is characterized by the existence of several relatively long quasi-linear chains of oceanic islands and seamounts$^{12}$. The oceanic islands and seamounts of the southeast Pacific lie within one of the most unexplored areas in the world. Most of the ~940 southeast Pacific seamounts are located along the Salas y Gómez (SGR) and Nazca (NR) ridges$^{13, 14}$, and only 24 (<3%) of the SGR and NR seamounts have been studied. Most of these studies were conducted by expeditions of the former USSR between 1973 and 1987 in international waters from ~80°W to 101°W at 162 to 1900 m depth$^{15}$. They reported unprecedented rates of endemism for benthic communities (e.g., ~41% for fishes and ~46% for invertebrates) that were similar to those rates observed for seamounts of the remote Norfolk Ridge, near New Caledonia$^{16}$. Despite studies being conducted on SGR and NR seamounts$^{15, 17-19}$, only one of the many seamounts within the Chilean EEZ adjacent to the Desventuradas Islands (San Felix and San Ambrosio islands) has been investigated. In addition to this seamount, known as Stockman Guyot$^{11, 15}$, the Desventuradas Islands were studied on two expeditions: Marine Research Cruises in Remote Areas (CIMAR 6) conducted by the National Oceanographic Committee of Chile (CONA) in 2001 and the "Pristine Seas Expedition" carried out by the National Geographic Society and Oceana in 2013. These studies, which were conducted at 40 to 2300 m$^{18, 20}$, revealed that the fauna of these islands are unique and irreplaceable, granting these areas a high conservation value$^{20, 21}$. These considerations prompted the creation, in 2015, of the large Nazca Desventuradas Marine Park (NDMP), protecting an area of ~300,000 km$^2$ including the Desventuradas Islands and seamounts located northwest of them$^{22}$. This marine park provides enormous advances in safeguarding the unique biodiversity of actual or past threats such as pelagic fishing of tuna and swordfish, as well as bottom trawling$^2$.

Of the available studies in the region, none has explored the importance of substrate type on species composition$^{23}$, nor the role of the habitat at different scales from centimeters to meters. Because distinct environmental conditions, such as internal tides and heterogeneity of substrates and oceanographic conditions (i.e., temperature, salinity and oxygen gradients), can differ considerably among and within seamounts$^{24}$, these and other variables can play a critical role in the distribution and behavior of megabenthos, operationally defined here as organisms with a body size > 1 cm. Differences in these environmental conditions can generate a mosaic of micro- and macrohabitats, with breaks in species distributions that can enhance biodiversity and drive changes in community composition within or among seamounts$^{23, 24}$.

The heterogeneity and topographic complexity of seamounts is affected by the presence of organisms that can change spatial conditions and, directly or indirectly, affect the availability of habitat and resources for other species$^{25}$. These organisms are called habitat-forming species and are characterized by being small-scale engineers of autogenous ecosystems. These organisms, and the biogenic structures they produce, contribute to increasing spatial complexity in soft-sediment benthic systems$^{25}$, promoting spatial heterogeneity on seamounts$^{26}$, but they are often easily overlooked because seamounts are sampled at much more general spatial scales$^{26}$. In recent years, the importance of sampling at smaller scales has been accentuated with the understanding that the characteristics of the bottom and the presence of habitat-forming species increases the heterogeneity of the habitat and fulfill the function of a structural habitat$^{27}$.

Several studies of continental shelf, slope and submarine canyon habitats have specified relationships of megafauna with certain microtopographic (microhabitat) characteristics (e.g., depressions, burrows and sessile fauna)$^{28-30}$. Because many of these studies lacked a standard habitat classification scheme, these observations cannot be compared efficiently among ecosystems and studies. To explore the importance of seabed habitat, Greene et al.$^{31}$ proposed a standardized classification, which divides the habitat according to depth, size and biogenic and abiotic variables, for in situ observations of the seabed. Within the size category are the subcategories of macrohabitat and microhabitat. Macrohabitats range in size from one to ten meters and include seafloor features such as boulders, crevices, cracks, caves, scarps, sinkholes and bedrock corals (solitary and reef-building). Meanwhile, microhabitats operate at centimeters and smaller scales and include seafloor features such as sand, silt, gravel, pebbles, small cracks, crevices and fractures$^{31, 32}$.

To explore relationships among micro- and macrohabitats and the fauna living on the summit of the seamounts and upper slope of oceanic islands of the NDMP, we describe the benthic habitats and associated fauna of the upper slope (50–370 m) of Desventuradas Islands and the summits and upper slopes (150–305 m) of seven nearby seamounts. We conducted remotely operated vehicle (ROV) surveys complemented with benthic fauna collections using an Agassiz trawl, aiming to achieve the following objectives: (1) to describe the benthic habitats of the oceanic islands and seamounts of the NDMP at different spatial scales (i.e. at centimeter to meter scale), (2) to describe the composition and the diversity of the benthic megafauna on these oceanic islands and seamounts, and (3) to evaluate the relationships between macrohabitat and megafauna species diversity and composition as well as feeding mode and movement type of the benthic communities. To our knowledge, this study is the first to describe benthic habitats and associated fauna at seamounts and oceanic islands of the NDMP.

## Materials and methods

**Study area.** The NR and SGR together form a sequential chain of seamounts of volcanic origin with an extension of ~2900 km$^{33}$ that is ~100 km wide along the SGR$^{34}$ and ~300 km wide along the NR$^{35}$. The Desventuradas Islands, constituted by San Felix and San Ambrosio islands (~26°S, 80°W), are located ~970 km off the coast of the Atacama Region, northern Chile. Seamounts of the Desventuradas Islands region are located northwest of these islands at the intersection of SGR and NR (Fig. 1)$^{36}$.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Location of the benthic stations sampled during CIMAR 22 cruise (13 October to 14 November 2016). Locations with Agassiz trawls and ROV deployments (filled upward triangle), only ROV deployments (filled square), and only Agassiz trawls (white upward triangle). Map generated using QGIS (version 3.10)⁷².</figcaption></figure>

**Data collection.** From 22 October to 1 November 2016, benthic surveys were performed at five stations on the upper slope of the Desventuradas Islands and the summits of seven nearby seamounts of the NR (Fig. 1). A Commander MK2 (Mariscope Meerestechnik, Kiel, Germany) ROV equipped with an HD Camcorder (Panasonic SD 909, resolution 1920×1080, 30 fps), a GoPro Hero camera (resolution 2704×1520, 60 fps) and laser pointers (10 cm apart) was used to survey benthic habitats and communities. A total of 11 ROV dives were conducted at ~40 to 370 m depth (Table S1); five on the slope of Desventuradas Islands (5 h of total bottom time) and six on the seamounts (6 h of total bottom time). The ROV observations were descriptive and did not aim to obtain quantitative data due to the potential overlapping of the paths of the surveys.

Benthic trawling was performed at stations where either multibeam or ROV observations indicated suitable soft bottoms. Trawls were conducted for ~4 to 15 min (bottom contact) with a constant speed of ~3 knots at each station (Table S1). The gear used was a modified Agassiz trawl, with a mouth of 1.5 m × 0.5 m (width × height) fitted with a net of 12 mm mesh at the cod end. A total of 10 trawls were conducted; four on the slope of the Desventuradas Islands and six on the seamount summits, with sampling depths ranging from 133 to 340 m (Table S1). Both trawl and ROV deployments were possible at all stations, except at ST18 and SF2 (only ROV data because the bottom was too rocky to safely conduct trawls) and at SF5 (only trawl data because strong currents prohibited ROV deployment) (Fig. 1).

Environmental data (i.e., temperature, salinity and oxygen) used in our analysis were obtained from the global three-dimensional CSIRO Atlas of Regional Seas (CARS) climatology (2009 version) as presented by Mecho et al.³⁷.

**Habitat description based on ROV videos.** The oceanic islands (OI) and seamounts (SM) were considered as two different environments (corroborated by statistical analyses presented below) and were referred to as subsystems, based on differences in physiography and depth. The characteristics of the seafloor habitats for each subsystem were described based on ROV data. Videos for each dive were analyzed in a time-lapse mode at half normal speed in VLC media Player 3.0.11³⁸. Habitats were classified at different spatial scales: (1) macrohabitats (i.e., meter-scale) and (2) microhabitats (i.e., centimeter-scale). The habitat types were described by geomorphology, including by slope (i.e., flat = 0°-5°, sloping = 5°-30°, steeply sloping = 30°-60°, vertical = 60°-90°), sediment type (i.e., rock, silty sediments, mixed substrate, coarse sand), texture (i.e., low rugosity, moderate rugosity, high rugosity) and "modifier" elements (i.e., biological communities, sedimentation and bio-perturbation) following Greene et al.³¹,³². The sections of the videos in which the ROV was steadily moving between 30 cm and 1 m above the ground along the bottom were used for the description of the macrohabitats (field of view per frame ~ 3 m²), and the sections in which the ROV remained sitting on the bottom (field of view per frame ~ 0.05 m²) were used
for the description of the microhabitats. The field of view per frame was estimated using the laser pointers of the ROV (positioned 10 cm apart).

**Benthic and demersal fauna description.** Specimens collected with the Agassiz trawl were preliminarily sorted, counted, and preserved onboard in 95% ethanol. Definitive counts and identification to the lowest possible taxonomic unit (hereafter operational taxonomic units, or OTUs) were performed at Sala de Colecciones Biológicas Universidad Católica de Norte (SCBUCN), where specimens were then assigned an ID number and cataloged. The observed fauna was classified according to relevant literature and previous reports for the area (e.g.,¹⁵, ³⁹⁻⁴⁵). Taxonomic assignments were further validated with resources such as the World Register of Marine Species (Worms, 2020. http://www.marinespecies.org) and the Ocean Biodiversity Information System (OBIS, 2020, Intergovernmental Oceanographic Commission of UNESCO. www.iobis.org).

Semi-quantitative data of abundance for each OTU were calculated based on the swept area for each trawl. The swept area (Abᵢ) was established as the product of the trawling speed (Vᵢ), effective trawl time (tᵢ) and width of the mouth of the trawl (Ahᵢ) and abundance was standardized to 10 m² following Barriga et al.⁴⁶:

$$Ab_i = t_i V_i Ah_i$$

The ROV recordings followed an exploratory methodology which restricted the data to the presence/absence of OTUs at each station. Videos for each transect were viewed at half their normal speed in VLC. The VLC tool "Interactive Zoom" was used to magnify images to observe the diagnostic characteristics of the taxa. Observations that were blurry or too distant were omitted because they did not provide sufficient details for identification.

All the sampling was performed under permission Res. Ext N°3685/2016 from SUBPESCA (National Fishing Authority of Chile) to Universidad Católica del Norte.

**Data analysis.** All analyses were performed using the R software, version 4.0.3⁴⁷ using the “vegan v2.5-6” package⁴⁸. Because both trawl and ROV data were not available for all stations (Table S1), the analysis approach considered either only trawl data or combined trawl and ROV data as follows. Only trawl data was used for the estimation of Shannon diversity (H') and Pielou's evenness (J) since abundance data is required. OTU richness (or total S) of each station was estimated based on the combined trawl and ROV data (except for ST18, SF2 and SF5 stations). In all cases, only "living" organisms were considered (i.e., empty shells or skeletons excluded) for the analyses. Species accumulation curves were constructed for each subsystem (OI and SM) and data gathering approaches (only trawl and combined trawl and ROV data) to estimate the rarified number of species and to assess sampling effort⁴⁸. The OTU richness and diversity indexes of the two subsystems (OI and SM) were compared using an unpaired two-sample Wilcoxon test.

To evaluate if there were differences in OTU composition among subsystems and stations a permutational analysis of variance (PERMANOVA), based on combined trawl and ROV data (ST18, SF2 and SF5 stations were excluded), was conducted. The routine “adonis2” on the Bray-Curtis index of similarity of the raw OTU presence/absence data (10,000 runs) with subsystems (two levels: OI and SM) and maximum depth as factors was used. Patterns in the structure of benthic communities among subsystems and stations were visualized using a cluster analysis UPGMA hierarchical clustering (using the functions “desvest” and “hclust”) and non-metric multidimensional scaling (nMDS; using the function “metaMD”) analyses based on Bray-Curtis dissimilarity elaborated from OTU presence/absence data. To analyze differences in functional diversity among sampling stations, OTUs were assigned to feeding modes (i.e., suspension, depositor, grazer, opportunistic or predator) and movement type (i.e., sessile, burrower/tube dweller, crawler or swimmer) based on functional traits and associated categories as proposed by Jones & Frid⁴⁹. A canonical correspondence analysis (CCA) was applied using the R function “cca” to relate the set of environmental parameters to: (1) OTUs presence/absence data, (2) feeding mode composition, and (3) movement type. Environmental variables included water hydrographic variables (temperature and oxygen), type of bottom (rock, mud and sand), depth, and geographic position (longitude and latitude). Salinity was not considered within the environmental variables in the nMDS and CCA because it co-varies with the oxygen. nMDS and CCA analysis were graphed using the library ggplot2 v3.3.2 in R⁵⁰.

## Results

**Description of macrohabitats.** Oceanic island macrohabitats (~43 to 370 m depth) were classified to one of three types (Table 1): (1) cobbles and bedrocks, predominated by sea urchins (i.e., *Centrostephanus sylviae*, Fig. 2A), (2) coarse sand, with anemones (mainly *Hormathia* sp.) (Fig. 2B), and (3) silty sediments with the presence of bioturbation (e.g., burrows) ripples with a maximum height of ~10 cm (Fig. 2C).

Seamount macrohabitats (~150 to 305 m depth), except at SF2, consisted of coarse sand and rhodolith beds (Fig. 2D,E, Table 1) with large patches of anemones and sea pens. Seamount SF2 substantially differed from the others by the predominance of hard substrates and a steep slope of compact cemented rock, with small patches of coarse sand and large rock formations. The sessile fauna of seamount SF2 consisted of a multitude of hydrozoans, sea pens, and small corals, among other sessile fauna associated with rocky substrates. Also, the presence of longitudinal trenches in the cemented rock gave shelter to numerous fishes (e.g., *Lotella* cf. *fernandeziana*, *Helicolenus lengerichi* and *Scorpaena thomsoni*) and large crustaceans (e.g., the Chilean jagged lobster *Projasus bahamondei* and the Juan Fernández carrier crab *Paromola rathbuni*) (Fig. 2F).

**Description of microhabitats.** Oceanic island microhabitats were classified into one of three types (Table 1): (1) joints, cracks, crevices, and overhangs (differentially eroded) covered by incrusting algae, sponges and stony corals (Fig. 3A) at 43–50 m; (2) coarse sand with sea pens and anemone colonies, which provide
<table><caption>**Table 1.** Description of benthic habitats of the oceanic islands (Desventuradas). Classification follows Greene et al.¹ and Greene et al.³². Subsystem OI = oceanic islands (upper slope) and SM = seamounts (top).</caption><thead><tr><td>Stations</td><td>Subsystem</td><td>Class (macrohabitat)</td><td>Subclass (microhabitat)</td><td>Environmental variables</td><td>Modifier</td></tr></thead><tbody><tr><td>ST17</td><td>OI</td><td>Flat bottom</td><td>Slope: flat (0-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.4<br>Temperature (C°): 13.9<br>Oxygen (mL/L): 5</td><td>Bottom morphology: regular-continuous<br>homogeneous bottom with little relief<br>with organic debris (coquina)<br>Biological processes: conspicuous<br>microhabitats/communities of anemone,<br>sea pens and kelp patch</td></tr><tr><td>ST18</td><td>OI</td><td>Sloping bottom, with vertical sections</td><td>Slope: steeply sloping (30°-60°)<br>Texture: moderate rugosity to high<br>rugosity surface<br>Sediment types: bedrock and boulder</td><td>Salinity: 34.7<br>Temperature (C°): 17.7<br>Oxygen (mL/L): 5.5</td><td>Bottom morphology: irregular bottom<br>with structures fractured, faulted and<br>folder<br>Bottom texture: occupied by clast or<br>rock<br>Biological processes: conspicuous<br>microhabitats/communities of coralline<br>incrusting algae, sponge and anemone</td></tr><tr><td>ST20</td><td>OI</td><td>Type I: steep slope bottom<br>Hard substrates and steep slope of<br>compact cemented rock<br>Type II: flat bottom</td><td>Type I<br>Slope: sloping (5°-30°)<br>Texture: moderate rugosity<br>Sediment types: bedrock and boulder<br>Type II<br>Slope: flat (0°-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.3<br>Temperature (C°): 12.9<br>Oxygen (mL/L): 4.7</td><td>Type I<br>Bottom morphology: irregular bottom<br>with structures fractured, faulted and<br>folder<br>Bottom texture: occupied by clast or<br>rock<br>Biological processes: conspicuous<br>microhabitats/communities of coralline<br>incrusting algae, sponge and anemone<br>Type II<br>Bottom morphology: regular-continuous<br>homogeneous bottom with little relief<br>Biological processes: conspicuous micro-<br>habitats/communities of anemones and<br>sea pens</td></tr><tr><td>ST21</td><td>OI</td><td>Flat bottom</td><td>Slope: flat (0°-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.3<br>Temperature (C°): 12.9<br>Oxygen (mL/L): 4.7</td><td>Bottom morphology: regular-continuous<br>homogeneous bottom with little relief<br>Biological processes: conspicuous micro-<br>habitats/communities of anemones and<br>sea pens</td></tr><tr><td>ST22</td><td>OI</td><td>Sloping bottom with bedform-sediment waves</td><td>Slope: sloping (5°-30°)<br>Texture: very low rugosity<br>Sediment types: silty sediments</td><td>Salinity: 34.5<br>Temperature (C°): 8.4<br>Oxygen (mL/L): 1.3</td><td>Bottom morphology: irregular bottom<br>(Undulated surface -ripples) with sedi-<br>ment waves (10 cm to amplitude)<br>Bioturbation: burrows and excavations</td></tr><tr><td>SF2</td><td>SM (top)</td><td>Steep slope bottom<br>hard substrates and steep slope of com-<br>pact cemented rock</td><td>Slope: steeply sloping (30°-60°)<br>Texture: moderate rugosity to high<br>rugosity surface, rock bottom</td><td>Salinity: 34.5<br>Temperature (C°): 10.4<br>Oxygen (mL/L): 1.6</td><td>Bottom morphology: rock bottom with<br>smalls accumulation of coarse sand<br>Biological processes: conspicuous<br>microhabitats/communities of sponges,<br>sea pens and white corals</td></tr><tr><td>SF6</td><td>SM (top)</td><td>Flat bottom</td><td>Slope: flat (0-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.3<br>Temperature (C°): 13.5<br>Oxygen (mL/L): 4.4</td><td>Bottom morphology: regular- continu-<br>ous homogeneous bottom with little<br>relief<br>Biological processes: conspicuous micro-<br>habitats/communities of anemones,<br>hydrozoan colonies, biocenosis (sea<br>urchin skeletons)</td></tr><tr><td>SF7</td><td>SM (top)</td><td>Flat bottom</td><td>Slope: flat (0-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.4<br>Temperature (C°): 13.5<br>Oxygen (mL/L): 4.4</td><td>Bottom morphology: regular- continu-<br>ous homogeneous bottom with little<br>relief<br>Biological processes: conspicuous micro-<br>habitats/communities of anemones,<br>rhodoliths, sponges, thanatocoenosis<br>(coquina- Pteropod shells)</td></tr><tr><td>SF8</td><td>SM (top)</td><td>Flat bottom</td><td>Slope: flat (0-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.4<br>Temperature (C°): 10.3<br>Oxygen (mL/L): 2.1</td><td>Bottom morphology: regular- continu-<br>ous homogeneous bottom with little<br>relief<br>Biological processes: thanatocoenosis<br>(coquina- bivalve shells); Bioturbation<br>(burrows and excavations)</td></tr><tr><td>SF9</td><td>SM (top)</td><td>Flat bottom</td><td>Slope: flat (0-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.2<br>Temperature (C°): 10.3<br>Oxygen (mL/L): 4.3</td><td>Bottom morphology: regular- continu-<br>ous homogeneous bottom with little<br>relief<br>Biological processes: conspicuous micro-<br>habitats/communities of anemones,<br>rhodoliths, sponges, thanatocoenosis<br>(coquina- Pteropod shells)</td></tr><tr><td>SFX</td><td>SM (top)</td><td>Flat bottom</td><td>Slope: flat (0-5°)<br>Texture: smooth surface<br>Sediment types: mixed sediment, coarse sand</td><td>Salinity: 34.4<br>Temperature (C°): 11.5<br>Oxygen (mL/L): 2.6</td><td>Bottom morphology: regular- continu-<br>ous homogeneous bottom with little<br>relief<br>Biological processes: conspicuous micro-<br>habitats/communities of anemones,<br>amphinomid polychaetas and detrital<br>patch</td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Benthic macrohabitats on the upper slope of Desventuradas islands (A–C) and summits of seamounts (D–F) at the Nazca-Desventuradas Marine Park. (A) Irregular rock bottom with fractures and faults, (B) regular continuous homogeneous bottom with coarse sand dominated by anemones (*Hormathia* sp), (C) silty sediments with the presence of bioturbation ripples, (D,E) regular continuous homogeneous bottom with little relief, coarse sand and rhodoliths, dominated by sponges (unidentified demosponges) and anemones (*Hormathia* sp. and cerianthids), (F) irregular rock bottom with structures fractured, faulted and folded, characterized by sea pens (*Scleroptilum* sp.), hydrozoan corals (*Stylaster* sp.) and some fishes (*Helicolenus lengerichi*). Scale bar: 25 cm. Image credits: ESMOI/OCEANA.</figcaption></figure>

harbor for fishes and invertebrates at 128–150 m (the most common habitat type observed on oceanic islands) (Fig. 3B), and (3) silty sediment bottom, absence of epibenthic cnidarians (e.g., anemones or sea pens). Besides, type 3 microhabitat was observed to have small holes on the surface of the sediment at 220–370 m, possibly caused by infaunal excavators. (Fig. 3C).

Cnidarians and sponges are the most conspicuous biological components of seamount microhabitats (~150 to 305 m depth), which included fields of sea pens (*Protoptilum* sp. and *Scleroptilum* sp.) and colonies of anemones (*Hormathia* sp. and Cerianthid), except at SF2. Besides, we also reported biogenic structures on coarse sand beds such as coquina of pteropod (Fig. 3D) and rhodoliths (Fig. 3E). Although the SF2 seamount is characterized by a hard and fractured substrate (faults and folds) covered with small corals (*Stylaster* cf. *marenzelleri*) and green
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Benthic microhabitats on the upper slope of Desventuradas islands (A–C) and summits of seamounts (D–F) at the Nazca-Desventuradas Marine Park. (A) Rocky bottom dominated by incrusting red algae, sea urchin *Centrostephanus sylviae*, the cracks and faults of the rocky bottom, are used as habitats for moray eels (*Gymnothorax porphyreus*) and small fishes (*Pseudolabrus cf. gayi*), (B) coarse sand flat bottom dominated by anemones (*Hormathia* sp.), (C) soft-sediment with burrows and excavations, (D) microhabitats of cerianthid anemones and hydrozoan colonies, (E) unattached nodules of crustose coralline red algae (Rhodoliths), demosponges, and coquina of pteropod shells, (F) irregular rock bottom with structures fractured, faulted and folded, characterized by sea pens (*Scleroptilum* sp.) and hydrozoan corals (*Stylaster* sp.). Scale bar: 10 cm. Image credits: ESMOI/OCEANA.</figcaption></figure>

sponges, small patches of coarse sand were observed between the rocks with sea pens of the genus *Scleroptilum* (Fig. 3F).

**Benthic community.** In total, 2414 individuals (only living specimens) were collected by trawling and assigned to 95 OTUs (OI: 44 OTUs and SM: 75 OTUs, 24 OTUs shared between subsystems), belonging to nine phyla (Table S2). In terms of total abundance, Cnidaria was the predominant major faunal group, comprising 29% of the total capture, followed by Porifera (26%), Arthropoda (13%), Echinodermata (13%), Mollusca (7%) Annelida (7%), Sipuncula (5%), Chordata (1%) and Bryozoa (0.2%). ROV observations resulted in a total of
61 OTUs (OI: 46, SM: 38, shared between subsystems: 23), of which 23 OTUs were observed only with this approach, 17 of these 23 OTUs were fish.

Combining collected specimens and ROV observations allowed allocation of a total of 118 OTUs (OI: 65, SM: 84, and shared between subsystems: 31). An additional 13 OTUs were identified just by shells or skeleton remains (e.g., pteropod or other mollusk shells), but were not included in the analyses. Invertebrates represented 75% of the OTUs, whereas vertebrates represented 25% (Fig. 4). Chordata (25%), Arthropoda (23%), Mollusca (13%) and Echinodermata (13%) were the most OTU-rich taxa in the total S data. Invertebrates were assigned to eight phyla: Arthropoda (23%), Mollusca (13%), Echinodermata (13%), Cnidaria (12%), Annelida (10%), Porifera (3%), Sipuncula (1%), Bryozoa (1%). Vertebrates (Chordata) were assigned to two classes of fishes (Table S2): Elasmobranchii (2 OTUs) and Actinopterygii (28 OTUs). In terms of overall OTU richness, fishes predominated at Desventuradas Islands (43%) and arthropods at seamounts (30%).

Combined OTU richness (total S data) for oceanic islands was highest at station ST20 (S=35), representing almost half of the total OTUs reported at the oceanic islands, followed by station ST21 (S=15). The lowest OTU richness at oceanic islands was observed for station ST18 (S=10) (Table 2). OTU richness at seamounts was highest at station SF9 (S=50), which represents ~40% of the total OTUs reported for seamounts, followed by station SF7 (S=35). The lowest OTU richness was observed at seamount stations SF8 (S=8) and SFX (S=13) (Table 2). However, no differences in total S (Wilcoxon test, W = 6.5, *p* = 0.284) or diversity (H'; Wilcoxon test, W = 8, *P* = 0.456) were observed between subsystems. Stations at oceanic islands, in general, were less diverse (*H'* = 1.7, mean) in comparison to the seamounts (*H'* = 1.9). The diversity for oceanic islands was highest at station ST20 (*H'* = 2.5) followed by station ST22 (*H'* = 1.6). The lowest diversity at oceanic island stations was observed at ST21 (*H'* = 1.3). The diversity at seamounts was highest at station SF7 (*H'* = 2.5), followed by station SF5 (*H'* = 2.2). The lowest diversity at seamounts was observed at station SFX (*H'* = 1.3). The evenness values were similar between subsystems (J: 0.7 for SM and 0.8 for OI, mean; Wilcoxon test, W = 15, *P* = 0.594). The evenness for oceanic islands was highest at station ST17 (*J* = 1.0) followed by station ST22 (*J* = 0.9). The lowest evenness at oceanic island stations was observed at ST21 (*J* = 0.5). The evenness at seamounts was highest at stations SF7 and SF8 (*J* = 0.9). The lowest evenness at seamounts was observed at station SF9 (*J* = 0.5) (Table 2). Species accumulation curves combining trawl and ROV data for each subsystem (seamount and oceanic island) showed that a sufficient level of sampling was not attained based on the species accumulation curves not reaching asymptotes (Fig. S1).

**Community comparisons.** The cluster (Fig. 5A) and the nMDS (Fig. 5B) analyses with combined data, revealed differences in the community structure between the stations, which clustered in two groups: (1) oceanic islands and (2) seamounts. This difference was corroborated by the PERMANOVA analysis based on combined data. Subsystem types (R² = 0.273, *F* = 3.234, *p* < 0.001) and depth (R² = 0.180, *F* = 2.132, *p* = 0.002) were significant variables in structuring the benthic communities of NDMP. The interaction between subsystem type and depth was not significant (R² = 0.126, *F* = 1.497, *p* = 0.097) (Table S3).

Besides, when performed for each type of sampling alone (i.e., only trawl or only ROV data), a similar grouping pattern was observed in the cluster and nMDS analyses. However, while oceanic islands and seamounts formed separate groups for only trawl data (Fig. S2), one seamount, SF2, grouped with oceanic islands when analysed by only ROV data (Fig. S3). In both cases, PERMANOVA analysis corroborated differences for subsystems (Trawl: R² = 0.155, *F* = 1.511, *p* = 0.002; ROV: R² = 0.205, *F* = 2.662, *p* = 0.001) (Table S3). Depth was a significant factor for ROV data (R² = 0.155, *F* = 2.018, *p* = 0.009) but not for trawl data (R² = 0.118, *F* = 1.146, *p* = 0.171). No significant interaction was observed between factors for either trawl or ROV data (Trawl: R² = 0.126, *F* = 1.497, *p* = 0.339; ROV: R² = 0.102, *F* = 1.333, *p* = 0.196) (Table S3).

The first two axes of the CCA represented ~38% of the total variance and two groups of stations were evident (Fig. 5C). The position and spacing of stations along CCA axis 1 shows that stations with negative scores were oceanic islands which positively correlated with increasing latitude, temperature, oxygen and rocky habitat. In contrast, the seamounts were grouped on the positive axis 1 and correlated with increasing depth, longitude and coarse sand habitats. Stations did not clearly distribute along CCA axis 2, for which oxygen, depth and substrate type were the primary contributors among the environmental data. Overlaying faunal contributions revealed Chordata were mainly associated with Desventuradas islands stations, whereas invertebrates such as Arthropoda, Echinodermata and Mollusca were associated with seamounts along the positive axis 1 of the CCA. A large part of the invertebrates (mollusks: 25 OTUs, crustaceans: 18 OTUs and echinoderms: 8 OTUs), in particular, the gastropods *Chryseofusus kazdailisi*, *Calliostoma* sp., *Atrimitra isolata*, *Cryptogemma praesignis*, the crustaceans *Munida diritas*, *Heteronida* sp., *Paragiopagurus* spp., *Ebalia sculpta*, and the irregular urchins *Brissopsis* sp. and *Spatangus* sp., were exclusively associated with seamounts F6, SF7 and SF9. Conversely, most of the fishes, including *Tetronarce* sp., *Callyonimus* sp., *Parapercis* cf. *dockinsi*, *Seriola lalandi*, *Nemadactylus gayi*, and *Gonorynchus greyi*, and a small portion of invertebrates (mollusks: 2 OTUs, crustaceans: 6 OTUs and echinoderms: 5 OTUs), including the bivalve *Arca cf. fernandezensis*, stomatopod *Hemisquilla ensigera*, snapping shrimp *Alpheus* cf. *romensky*, dollars *Clypeaster* spp. and the sea star *Astrostole platei*, were associated with stations around the islands (ST20 and ST21).

Regarding functional traits of benthic communities, the first two axes of the CCA analysis, based on movement type (Fig. 6A) and feeding modes (Fig. 6B), accounted for 92.2% and 81.9% of the total variance, respectively. This analysis showed that oceanic islands were mainly dominated by mobile fauna and predators (e.g., the fishes *Squalus mitsukurii*, *Lotella* cf. *fernandeziana*, *Helicolenus lengerichi* and *Scorpaena thomsoni*, and large crustaceans (e.g., *Projasus bahamondei* and *Paromola rathbuni*), whereas seamounts farther from the islands were dominated by sessile and hemisessile invertebrates (i.e. crawlers and burrowers/tube dwellers) that were mainly suspension and deposit feeders (e.g., *Hormathia* sp., *Protoptilum* sp., *Paragiopagurus* spp.).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Representative OTUs observed at the upper slope of Desventuradas islands and seamounts within the Nazca-Desventuradas Marine Park. (A) Sea pen *Protoptilum* sp. and anemone *Hormathia* sp., (B) anemone (*Hormathia* sp.) and hermit crab *Paragiopagurus boletifer*, (C) *Cryptogemma praesignis*, (D) *Anthenoides* sp., (E) *Clypeaster isolatus*, (F) *Notopogon fernandezianus*, (G) anemone colony *Hormathia* sp., (H) *Chloeia* sp., (I) *Paromola rathbuni*, (J) *Stereocidaris nascaensis* (left) and *Scrippsechinus fisheri* (right), (K) *Squalus mitsukurii*, (L) *Gonorhynchus greyi*, (M) tube anemone Ceriantharidae, (N) *Chryseofusus kazdailisi*, (O) *Projasus bahamondei*, (P) *Pseudarchaster* sp., (Q) *Tetronarce* sp., and (R) *Lotella fernandeziana*. Laser pointers, when visible, are 10 cm apart. Image credits: ESMOI/OCEANA.</figcaption></figure>
<table><caption>**Table 2.** OTU richness based on trawl data, total S (trawl+ROV) data and diversity and evenness indexes obtained from relative OTU abundances of only trawl data. Subsystem OI = oceanic islands and SM = seamounts. Only ROV (***).</caption><thead><tr><th>Site</th><th>Subsystem</th><th>OTUs richness (S)</th><th>Total S (trawl+ROV)</th><th>Diversity (H')</th><th>Evenness (J')</th></tr></thead><tbody><tr><td>ST17</td><td>OI</td><td>4</td><td>15</td><td>1.3</td><td>1.0</td></tr><tr><td>ST18</td><td>OI</td><td>-</td><td>10***</td><td>-</td><td>-</td></tr><tr><td>ST20</td><td>OI</td><td>27</td><td>35</td><td>2.5</td><td>0.8</td></tr><tr><td>ST21</td><td>OI</td><td>13</td><td>15</td><td>1.3</td><td>0.5</td></tr><tr><td>ST22</td><td>OI</td><td>6</td><td>11</td><td>1.6</td><td>0.9</td></tr><tr><td>Mean</td><td>OI</td><td>12</td><td>17</td><td>1.7</td><td>0.8</td></tr><tr><td>SF2</td><td>SM</td><td>-</td><td>14***</td><td>-</td><td>-</td></tr><tr><td>SF5</td><td>SM</td><td>32</td><td>-</td><td>2.2</td><td>0.6</td></tr><tr><td>SF6</td><td>SM</td><td>19</td><td>27</td><td>2.2</td><td>0.7</td></tr><tr><td>SF7</td><td>SM</td><td>19</td><td>35</td><td>2.5</td><td>0.9</td></tr><tr><td>SF8</td><td>SM</td><td>6</td><td>8</td><td>1.6</td><td>0.9</td></tr><tr><td>SF9</td><td>SM</td><td>45</td><td>50</td><td>1.9</td><td>0.5</td></tr><tr><td>SFX</td><td>SM</td><td>11</td><td>13</td><td>1.3</td><td>0.6</td></tr><tr><td>mean</td><td>SM</td><td>22</td><td>24</td><td>1.9</td><td>0.7</td></tr></tbody></table>

Discussion

This study constitutes the most complete and updated description of benthic habitats and fauna for the upper slope of Desventuradas Islands and summits of nearby seamounts within the NDMP. ROV videos allowed the first description of benthic habitats of the Desventuradas Islands and their surrounding seamounts at a fine scale (meters to centimeters). Although we are aware of the limitations of the present study (e.g., only one ROV dive and/or trawl per site), we consider that our information presents a fundamental base-line knowledge of the benthic fauna and environment in this relatively understudied region of the Pacific Ocean, with a high value for conservation.

Of the three macrohabitats (Fig. 2, Table 1) observed in this study, relatively homogeneous coarse sand with patches of habitat-forming sessile species and rhodoliths was the most commonly observed at oceanic islands and seamounts. These observations are consistent with the description of the macrohabitat made by “Pristine Seas Expedition” at depths of 20 to 350 m⁵¹ at the Desventuradas Islands and reports of substrates on other seamounts of the SGR¹⁵,¹⁹ as well as of the eastern Pacific Ocean⁵², northeast Atlantic⁵³ and equatorial Atlantic⁹. For example, Parin et al.¹⁵ reported, based on collections and observations made from a manned submersible, that the upper slopes of SGR seamounts located > 400 km NW of Desventuradas Islands had a soft bottom composed of biogenic sands derived from shells of local mollusks with some pelagic pteropods, foraminiferal material and fossil corals. We observed a change in the kind of sediment/ bottom type with predominantly coarse sands at the 150–180 m stations and predominantly fine sands at 220–370 m. This pattern of sediment change with depth has also been recorded in other deep-sea areas⁵⁴⁻⁵⁷.

The rocky macrohabitat observed on oceanic island ST18 (Fig. 2A) is similar to other islands located in the southwest Pacific. For example, Lord Howe Island (~33°S, 159°E, ~600 km E of Australia) is also an isolated, endemism hotspot with several similarities to Desventuradas Islands²⁰. The macrohabitat of Lord Howe Island at a depth of 50 m is characterized by having a rocky bottom covered with encrusting red algae and being dominated by urchins (*Tripneustes* and *Centrostephanus*)²⁰,³⁰.

The rocky substrate (i.e., steep slope of compact cemented rock) observed at seamount SF2 (Fig. 2F) is consistent with the macrohabitat of other seamounts of the NR (e.g., Professor Mesyatzev, Soldatov, Ikhtiar and Ekliptika; ~81–83°W) reported by Parin et al.¹⁵. This macrohabitat is characterized by the presence of a longitudinal trench in the cemented rock, which gives shelter to numerous species of fishes (e.g., *Helicolenus lengerichi* and *Scorpaena thomsoni*) and large crustaceans (e.g., the Chilean jagged lobster *Projasus bahamondei* and the Juan Fernández carrier crab *Paromola rathbuni*) (Fig. 2F).

Until recently, most seamount research related to benthic habitat heterogeneity and complexity has focused on larger-scale variations (macrohabitats)⁵⁸. Unlike hard-bottomed habitats, soft-bottom habitats have generally been considered unstructured, homogeneous environments; consequently, these habitats have been studied primarily from a landscape perspective (usually at a kilometer scale)⁵⁸. Hence, Clark et al.⁵⁹ suggested that the concept of seamounts as a single, relatively well-defined habitat type appears outdated, giving way to a growing recognition that within seamount variability can have different spatial scales. Although not statistically tested, our results suggest a higher faunal diversity associated with biogenic microhabitats compared to sites with few or no biogenic microhabitats (e.g., SF8 and SFX), further supporting the importance of such small-scale differences in understanding the mechanisms involved in maintaining species diversity and the linkages between habitat-forming species and users of these habitats (Table S4). The presence of multiple distinct morphologies of the deep seabed, such as, scars, channels, mud, rocks and sand, together with habitat-forming species (e.g., sponge, anemones and stylasterids corals) provide an environment that is home to numerous species, resulting in increases in diversity¹⁶,⁶⁰. Besides, habitat-forming species are known to promote local complexity and
<figure><img src="image_7.png" /><figcaption>**Figure 5.** (A) Cluster analysis (UPGMA method) based on Bray-Curtis dissimilarity, (B) non-metric multidimensional scaling (nMDS) ordination plot based on presence/absence data of the benthic megafauna, and (C) canonical correspondence (CCA). Analysis based on the community structure of benthic megafauna, grouped at major taxa level, from stations sampled at Desventuradas Islands (blue markers) and seamounts (red markers) of the Nazca Desventuradas Marine Park. Vectors in (B) and (C) represent contribution of environmental descriptors (salinity, oxygen, temperature, substrate type, depth, latitude and longitude), and ellipses in (B) represent the 95% confidence interval. Image generated using R software (version 4.0.3)⁴⁷.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Canonical correspondence analyses (CCA) analysis based on the community structure of benthic megafauna grouped by movement type (A) and feeding mode (B) from stations sampled at Desventuradas Islands (OI, blue markers) and seamounts (SM, red markers) of the Nazca Desventuradas Marine Park. Vectors represent the contribution of environmental descriptors (oxygen, temperature, substrate type, depth, latitude and longitude). Image generated using R software (version 4.0.3)⁴⁷.</figcaption></figure>

heterogeneity of the seafloor at small (microhabitat) to large (macrohabitat) scales and to provide suitable refuge and nursery zone for smaller mobile benthic organisms², 24, 60, 61. Consistent with this pattern, it appears that OTU richness may be higher where habitat-forming species, such as sea pens (*Protoptilum* and *Scleroptilum*), anemones (*Hormathia* sp.), tube-dwelling anemones and stylasterids corals were present (Table S4). Microhabitats within and among seamounts differed based on small modifications in the seafloor substrate and the presence of habitat-forming species such as anemones (*Hormathia* sp.), tube-dwelling anemones (cerianthids) and sea pens (*Protoptilum* sp. and *Scleroptilum* sp.). A marked pattern was observed between substrate type (soft vs. hard) and the type of habitat-forming species; anemones and sea pens of the genus *Protoptilum* were observed frequently on soft bottoms, whereas sea pen of the genus *Scleroptilum* and stylasterids corals were reported on hard bottoms. Similar patterns have been described for deeper habitats (> 400 m depth) of NR and SGR seamounts and oceanic islands¹⁵, ¹⁹ and also for New Caledonia, particularly on the Norfolk Ridge seamounts¹⁶.

Habitat-forming species were distributed as small patches that supported a higher concentration of fauna compared to the surrounding seafloor, providing a structure that can support more species than habitats without such complexity¹⁶. For example, at NDMP seamounts, squat lobsters (*Munida diritas* and *Phylladiorhynchus pusillus*), crabs (*Latreillia* sp.) and juveniles of *Caprodon longimanus* were observed hidden among anemones (Fig. 3D). This type of seamount microhabitat is mainly composed of coarse sand and rhodoliths with anemones, sea pens and tube-dwelling anemones, which likely provide nursery zones for some fishes, as suggested by ROV observations of the presence of juveniles of *C. longimanus* in these anemone patches (SF6 and SF7 seamounts). Another role of microhabitats is contributing to predator-prey dynamics; for example, the use of microhabitats not only plays a role in predator avoidance by providing shelter but can also play a role in ambush tactics by providing a focal site where the prey may not recognize the predator²⁸. In our study, squat lobsters (*Munida*
diritas) were observed hunting mysidaceans that foraged around anemone tentacles. In a recent description, Gallardo et al.⁶² mentioned that the white coloration of the squat lobster favors camouflage among coquina sediments. Similar observations have been reported by Auster et al.²⁸, who documented some fauna (e.g., red hake, ocean pout, longhorn sculpin, squid, conger eel, and black-bellied rosefish) using microtopographic features and resources of microhabitat to block visual and acoustic (i.e., proprioceptive) recognition of occupants in ambush tactics.

The NR and SGR are considered biologically-unique because of their remarkably high endemism¹¹, ²⁰. Previous studies have found that the Nazca-Desventuradas region has unique marine ecosystems that consist of a mixture of tropical, subtropical and temperate species with a strong affinity with the Indo-Pacific¹⁸, ²⁰. We report a total of 118 OTUs (131 if adding non-living specimens or shells), of which 50 invertebrates and 30 fishes were identifiable to species level (Table S2). The number of species reported was much higher than the review of Castilla⁶³ that reported 50 species (23 invertebrates and 27 fishes) and the Chilean expeditions (i.e., CIMAR 6 and Pristine Sea Expedition) that reported a total of 37 species. This increase in species number for the region is, most probably, due to the few previous scientific studies that focused on benthic megafauna⁶⁴. Although the sampling effort deployed in the area is still far from being satisfactory, as indicated by the species accumulation curves (Fig. S1), this study reports 43 genus-level OTUs not previously mentioned in the NDMP area (Table S2), of which only 14 were assigned to the level of species (new records for the area), with the remaining ones being potentially new species to science. New records from this study for the Nazca-Desventuradas Marine Park include: two genera of Pennatulaceans (*Protoptilum* and *Scleroptilum*); seven genera and five species of polychaetes (*Lanice* sp., *Chloelia* sp., *Eunice* decolorhami, *Phyllodoce* pseudopatagonica, *Trypanosyllis* cf. *zebra* and *Mesochaetopterus* minutus), bringing the total to 13 genera and five species¹⁵, ⁶⁵, ⁶⁶; one genus of Sipuncula (*Aspidosiphon* sp.); 14 genera and four species of gastropods (*Iniforis* cf. *limitaris*, *Atrimitra* isolata, *Chryseofusus* kazdailisi, *Cryptogemma* prae-signis, *Architectonica* karsteni), bringing the total to 40 genera and 31 species¹⁵, ²¹, ⁴³, ⁶⁷; 17 genera and 15 species of crustaceans (e.g., *Miersiella* haswelli, *Zarenkolambrus* minutus, *Z. epibranchialis*, *Heteronida* sp, *Munida* diritas and *Hispidolambrus* mironovi) (Table S2), bringing the total to 25 genera and 23 species⁴², ⁴⁴, ⁴⁵; and two genera of fishes (*Tetronarce* sp. and *Callionymus* sp.). Fishes is the best documented group in the NDMP area, including reports by Sepulveda⁶⁸, Pequeño & Lamilla⁴⁰, Dyer & Wesneat³⁹ and National Geographic & Oceana⁵¹ around the Desventuradas Islands from 5 to 2,215 m depth and the summit of the Stockman guyot (station SFX) at 200 to 370 m depth.

Therefore, it is probable that a greater sampling effort (i.e., higher areal coverage and including flanks and bases of seamounts) will continue revealing substantial numbers of novelties for the region. Although rare invertebrate species and possibly new species were observed in the present study, none of the new records appear to be endemic to a particular seamount (Table S2). However, finer-scale studies may help to uncover diversity patterns of species within and among seamounts in this region. Furthermore, because of limited sampling across the southeast Pacific, available information is still insufficient to assess the high endemism and regional affinity of the NDMP fauna.

Although out of the scope of the present study, we provide records that would help to better assess the hypothesis of a faunal transition zone or "break" proposed by Mironov & Detinova⁶⁹ and subsequently by Parin et al.¹⁵ and Mecho et al.³⁷ which divides the hard-bottom seamounts of south NR from the soft bottom seamounts of the east SGR. Parin et al.¹⁵ suggested a faunistic break at ~88°W. Mecho et al.⁷⁰ observed remarkable turnovers of echinoderms at ~101° and ~86°W of the SGR, where assemblages tended to differ more across seamounts, suggesting the effects of physical barriers to dispersion (e.g., currents) and habitat changes. Our data suggests a transition zone located at the intersection of SGR and NR (~82.5°W), as indicated by the separation among stations in the nMDS (Fig. 5B) and CCA (Fig. 5C) coinciding with the initial position (~82–84°W) proposed by Mironov & Detinova⁶⁹. Differences among the communities can also be seen in the functional groups and composition of the fauna of the seamounts of the NDMP separated by this approximate longitude (Fig. S2). Seamount SF2 is the only hard-bottomed seamount explored, with a mobile fauna with a type of feeding associated with predators such as the spiny lobster *Projasus bahamondei* and shark *Squalus mitsukurii*. These characteristics are shared with other seamounts such as Zvezda to Bolshaya further northeast in the NR¹⁵. In contrast, the other seamounts studied coincide more with the characteristics reported for the SGR because they have a soft bottom and sessile/hemisessile benthic megafauna with a generalist feeding mode¹⁵.

The CCAs (Fig. 6) suggest that differences in feeding mode and mobility of the megafauna appear to be linked to habitat type, depth, temperature and oxygen. Such patterns have been observed in other seamounts or deep habitats reported in literature¹⁵, ⁷⁰. The megafauna of the Desventuradas Islands was characterized by a comparable high frequency of occurrence and diversity of fishes (e.g., *Suezichthys rosenblatti*, *Paratrimma* sp., *Parapercis* cf. *dockinsi* and *Seriola* lalandi) and large predators such as the shark *Squalus mitsukurii*⁷. This predominance of fishes is likely the result of the islands having a biogenic structure that could support higher fish densities because of higher availability of nutrients, terrigenous contributions and eventual trophic subsidies by macroalgae (e.g., *Eisenia* cf. *cokeri*)⁶⁸, ⁷¹. In contrast, the seamount fauna was characterized by having hemisessile and sessile fauna that are opportunistic suspension feeders (e.g., *Hormathia* sp. *Protoptilum* sp. *Scleroptilum*, and *Stylaster* sp.) or deposit feeders (e.g., *Paragiopagurus* spp. and *Munida* diritas). This pattern suggests that the seamounts of NDMP might have short food webs and low guild complexity, i.e., the majority of the animals in seamounts consume similar resources at a low trophic level¹⁶, ²³. A similar pattern was described for the New Caledonian seamounts, where it has been hypothesized that the length of the food chain in seamounts is shorter than in other aquatic systems¹⁶. These authors also found that although the food web in seamounts is short, it is still relatively complex.
Conclusions

This study constitutes a first descriptive addition to the knowledge of benthic micro- and macrohabitats of the slope of Desventuradas Islands and summits of nearby seamounts within the NDMP, providing the most complete description of their benthic megafauna to date. The benthic habitats of oceanic island slopes and seamounts within the NDMP are home to a large diversity of habitats, including rock, rhodoliths, sandy bottoms, silty bottoms, vertical walls and caves. Our study highlights that faunal composition differed with seafloor habitat. For example, higher values of diversity and OTU richness were observed in areas with a predominance of habitat-forming species, which points to the need to study seamounts at finer scales (centimeters to meters) to better understand the links of habitat-forming species and faunal communities and the potential effects that would result from the loss of macro- and microhabitat-forming species. Currently, these habitats within the NDMP are in an apparent pristine state (i.e., without evident signs of human impacts, such as marine litter or trawled areas). Thus, these seamounts, now under protection, provide a regional unique opportunity to monitor effects of climate change without the confounding factor of human activities, such as fishing, deep seabed mining or bottom trawl fishing, and could also be useful for comparisons with other places in the world where such closures have been implemented for seamounts in the Azores, SW Indian Ocean, New Zealand, Hawaiian Islands and Palau, to mention a few.

Received: 5 June 2020; Accepted: 2 March 2021

Published online: 18 March 2021

References

1. Yesson, C., Clark, M. R., Taylor, M. L. & Rogers, A. D. The global distribution of seamounts based on 30 arc seconds bathymetry data. *Deep. Res. Part I Oceanogr. Res. Pap.* **58**, 442–453 (2011).

2. Preez, CDu., Curtis, J. M. R. & Clarke, M. E. The structure and distribution of benthic communities on a shallow seamount (Cobb Seamount, Northeast Pacific Ocean). *PLoS ONE* **11**, 1–29 (2016).

3. Auster, P. J. et al. Definition and detection of vulnerable marine ecosystems on the high seas: problems with the ‘move-on’ rule. *ICES J. Mar. Sci.* **68**, 254–264 (2011).

4. Watling, L. & Auster, P. J. Seamounts on the high seas should be managed as vulnerable marine ecosystems. *Front. Mar. Sci.* **4**, 1–4 (2017).

5. Cho, W. W. *Faunal Biogeography, Community Structure, and Genetic Connectivity of North Atlantic Seamounts* (Massachusetts Institute of Technology & Woods Hole Oceanographic Institution, 2008).

6. Rogers, A. D. *The Biology of Seamounts: 25 Years on*. *Advances in Marine Biology* vol. 79 (Elsevier, 2018).

7. Wagner, D. et al. The Salas y Gómez and Nazca ridges: a global diversity hotspot in need of protection. *28 (2020).*

8. Kvile, K. O., Taranto, G. H., Pitcher, T. J. & Morato, T. A global assessment of seamount ecosystems knowledge using an ecosystem evaluation framework. *Biol. Conserv.* **173**, 108–120 (2014).

9. Victorero, L., Robert, K., Robinson, L. F., Taylor, M. L. & Huvenne, V. A. I. Species replacement dominates megabenthos beta diversity in a remote seamount setting. *Sci. Rep.* **8**, 1–11 (2018).

10. Yesson, C. et al. Improved bathymetry leads to 4000 new seamount predictions in the global ocean. *UCL Open Environ. Preprint*, 1–12 (2020).

11. Gálvez Larach, M. Montes submarinos de Nazca y Salas y Gómez: una revisión para el manejo y conservación. *Lat. Am. J. Aquat. Res.* **37**, 479–500 (2009).

12. Jarrard, R. D. & Clague, D. A. Implications of Pacific Island and seamount ages for the origin of volcanic chains. *Rev. Geophys.* **15**, 57–76 (1977).

13. Chave, E. H. & Jones, A. T. Deep-water megafauna of the Kohala and Haleakala slopes, Alenuihaha Channel Hawaii. *Deep Sea Res. Part A Oceanogr. Res. Pap.* **38**, 781–803 (1991).

14. Kitchingman, A., Lai, S., Morato, T. & Pauly, D. How many seamounts are there and where are they located? In *Seamounts: Ecol- ogy, Fisheries & Conservation, Series 12* (eds Pitcher, T. J. et al.) 26–40 (Blackwell Publishing, 2008). https://doi.org/10.1002/97804 70691953.ch2.

15. Parin, N. V., Mironov, A. N. & Nesis, K. M. *Biology of the Nazca and Sala y Gómez submarine ridges, an outpost of the Indo-West Pacific fauna in the eastern Pacific ocean: composition and distribution of the fauna, its communities and history*. *Advances in Marine Biology* vol. 32 (1997).

16. Samadi, S., Schlacher, T. & Richer de Forges, B. Seamount benthos. In *Seamounts: Ecology, Fisheries and Conservation* (eds Pitcher, T. et al.) 119–140 (Wiley-Blackwell, 2007).

17. Mironov, A. N., Molodtsova, T. N. & Parin, N. V. Soviet and Russian studies on seamount biology. (2006).

18. Fernández, M., Pappalardo, P., Rodríguez-Ruiz, M. C. & Castilla, J. C. Síntesis del estado del conocimiento sobre la riqueza de especies de macroalgas, macroinvertebrados y peces en aguas costeras y oceánicas de Isla de Pascua e Isla Salas y Gómez. *Lat. Am. J. Aquat. Res.* **42**, 760–802 (2014).

19. Easton, E. E. et al. Chile and the Salas y Gómez Ridge. In *Mesophotic Coral Ecosystems* 477–490 (Springer, 2019). https://doi. org/10.1007/978-3-319-92735-0_27.

20. Friedlander, A. M. et al. Marine biodiversity in Juan Fernández and Desventuradas islands, Chile: global endemism hotspots. *PLoS ONE* **11**, e0145059 (2016).

21. Sellanes, J., Salisbury, R. A., Tapia, J. M. & Asorey, C. M. A new species of *Atrimitra* Dall, 1918 (Gastropoda: Mitridae) from seamounts of the recently created Nazca-Desventuradas Marine Park Chile. *PeerJ* **2019**, 1–16 (2019).

22. Gaymer, C. F. et al. *Plan General de Administración y su Valoración Económica. Informe final proyecto FIPA 2016-31 'Bases técnicas para la gestión del Parque Marino Nazca-Desventuradas y propuesta de Plan General de Administración' (2018).*

23. Clark, M. R. et al. The ecology of seamounts: structure, function, and human impacts. *Ann. Rev. Mar. Sci.* **2**, 253–278 (2010).

24. Henry, L. A. et al. Environmental variability and biodiversity of megabenthos on the Hebrides Terrace Seamount (Northeast Atlantic). *Sci. Rep.* **4**, 1–10 (2014).

25. Jones, C. G., Lawton, J. H. & Shachak, M. Organisms as ecosystem engineers. *Oikos* **69**, 373 (1994).

26. Morgan, N. B., Goode, S., Roark, E. B. & Baco, A. R. Fine scale assemblage structure of benthic invertebrate megafauna on the North Pacific Seamount Mokumanamana. *Front. Mar. Sci.* **6**, 1–21 (2019).

27. Davies, J. S. et al. Benthic assemblages of the Anton Dohrn Seamount (NE Atlantic): defining deep-sea biotopes to support habitat mapping and management efforts with a focus on vulnerable marine ecosystems. *PLoS ONE* **10**, 33 (2015).

28. Auster, P. J., Malatesta, R. J. & Larosa, S. C. Patterns of microhabitat utilization by mobile megafauna on the southern New England (USA) continental shelf and slope. *Mar. Ecol. Prog. Ser.* **127**, 77–85 (1995).
29. Uzmann, J. R., Cooper, R. A., Theroux, R. B. & Wigley, R. L. Synoptic comparison of three sampling techniques for estimating abundance and distribution of selected megafauna: submersible vs. camera sled vs. otter trawl. *Mar. Fish. Rev.* **39**, 11–19 (1977).

30. Valentine, J. P. & Edgar, G. J. Impacts of a population outbreak of the urchin *Tripneustes gratilla* amongst Lord Howe Island coral communities. *Coral Reefs* **29**, 399–410 (2010).

31. Greene, H. *et al*. A classification scheme for deep seafloor habitats. *Oceanol. Acta* **22**, 663–678 (1999).

32. Greene, H., O'Connell, V., Brylinsky, C. & Reynolds, J. Marine Benthic Habitat classification: What's Best for Alaska? In *Marine Habitat Mapping Technology for Alaska* (eds Reynolds, J. & Greene, H. G.) 169–184 (Alaska Sea Grant College Program University of Alaska Fairbanks, 2008). https://doi.org/10.4027/mhmta.2008.12.

33. Naar, D. F., Johnson, K. P., Wessel, D., Duncan, P. & Mahoney, J. Rapa Nui. 2001: Cruise report for Leg 6 of the Drift expedition aboard the R/V Revelle (2001).

34. Haase, K. M., Stoffers, P. & Garbe-Schönberg, C. D. The petrogenetic evolution of lavas from Easter Island and neighbouring seamounts, near-ridge hotspot volcanoes in the SE pacific. *J. Petrol.* **38**, 785–813 (1997).

35. Woods, M. T. & Okal, E. A. The structure of the Nazca Ridge and Sala y Gomez seamount chain from the dispersion of Rayleigh waves. *Geophys. J. Int.* **117**, 205–222 (1994).

36. Rodrigo, C., Foucher, N., Philippi, N. & Lara, L. E. Morfoestructuras volcánicas y sedimentarias de los montes submarinos de la región de las islas Desventuradas, basadas en el análisis de datos acústicos. 110–115 (2017).

37. Mecho, A. *et al*. Environmental drivers of mesophotic echinoderm assemblages of the Southeastern Pacific Ocean. *Front Mar. Sci.* **8**, 1–15 (2021).

38. VLC media player - Open Source Multimedia Framework and Player.

39. Dyer, B. S. & Westneat, M. W. Taxonomía y biogeografía de los peces costeros del Archipiélago de Juan Fernández y de las islas Desventuradas Chile. *Rev. Biol. Mar. Oceanogr.* **45**, 589–617 (2010).

40. Pequeño, G. & Lamilla, J. The Littoral Fish Assemblage of the Desventuradas Islands (Chile) Has Zoogeographical Affinities with the Western Pacific. *Glob. Ecol. Biogeogr.* **9**, 431–437 (2000).

41. Raines, B. & Huber, M. *Biodiversity Quadrupled-Revision of Easter Island and Salas y Gómez Bivalves*. *Zootaxa* **106** (2012).

42. Retamal, M. A. & Moyano, H. I. Zoogeografía de los crustáceos decápodos chilenos marinos y dulceacuícolas. *Lat. Am. J. Aquat. Res.* **38**, 302–328 (2010).

43. Sysoev, A. B. Gastropods of the family Turridae (Gastropoda:Toxoglosa) of the Nasca and Sala y Gómez underwater ridges. **124**, 245–260 (1990).

44. Zarenkov, N. A. Crabs of the family Leucosidae (subfamilies Ebalinae an Ilinae) collected in tropical water of Indian and Pacific oceans waters of Indian and Pacific oceans. *Bol. Nauk.* **10**, 16–26 (1969).

45. Zarenkov, N. A. Decapods (Stenopodidea, Brachyura, Anomura) of the underwater Nazca and Salas y Gómez Ridges. *Tr. Instituto Oceanol. AN USSR* **124**, 218–244 (1990).

46. Barriga, E., Salazar, C., Palacios, J., Romero, M. & Rodriguez, A. Distribución, abundancia y estructura poblacional del langostino rojo de profundidad Haliporoides diomedea (Crustacea: Decapoda: Solenoceridae). *Lat. Am. J. Aquat. Res.* **37**, 371–380 (2009).

47. R Core Team. R Core Team (2020). R: A language and environment for statistical computing. version 4.0.3. R Foundation for Statistical Computing, Vienna, Austria. https://www.r-project.org/ (2019).

48. Oksanen J *et al*. vegan: Community Ecology Package.R package version 2.5-7. https://cran.r-project.org/package=vegan (2020).

49. Jones, D. & Frid, C. L. J. Altering intertidal sediment topography: effects on biodiversity and ecosystem functioning. *Mar. Ecol.* **30**, 83–96 (2009).

50. Wickham, H. *ggplot2: Elegant Graphics for Data Analysis* (Springer, 2009).

51. National Geographic & Oceana. Islas Desventuradas. Biodiversidad marina y propuesta de conservación. 58 (2013).

52. Levin, L. A. & Nittrouer, C. A. Textural characteristics of sediment on deep seamounts in the eastern Pacific Ocean between 10°N and 30°N. In *Seamounts, Islands and Atolls*, 43 (eds Keating, B. *et al*. ) 187–203 (Geophysical Monograph, 1987).

53. Lourido, A., Parra, S. & Serrano, A. Preliminary Results on the Composition and Structure of Soft-Bottom Macrobenthic Communities of a Seamount: the Galicia Bank (NE Atlantic Ocean). *Thalassas* **35**, 1–9 (2019).

54. Flach, E., Muthumbi, A. & Heip, C. Meiofauna and macrofauna community structure in relation to sediment composition at the iberian margin compared to the goban spur (NE atlantic). *Prog. Oceanogr.* **52**, 433–457 (2002).

55. Levin, L. A. & Gooday, A. The deep Atlantic Ocean floor. In *Ecosystems of the Deep Oceans* (ed. Tyler, P.) 187–203 (Elsevier, 2003).

56. Thistle, D. The deep-sea floor: an overview. In *Ecosystems of the World, Ecosystems of the Deep Sea* (ed. Tyler, P. A.) 5–37 (Elsevier, 2003).

57. Louzao, M. *et al*. Historical macrobenthic community assemblages in the Avilés Canyon, N Iberian Shelf: Baseline biodiversity information for a marine protected area. *J. Mar. Syst.* **80**, 47–56 (2010).

58. Kon, K., Tsuchiya, Y., Sato, T., Shinagawa, H. & Yamada, Y. Role of microhabitat heterogeneity in benthic faunal communities in sandy bottom sediments of Oura Bay, Shimoda Japan. *Reg. Stud. Mar. Sci.* **2**, 71–76 (2015).

59. Clark, M. R., Schlacher, T. A., Rowden, A. A., Stocks, K. I. & Consalvey, M. Science priorities for Seamounts: research links to conservation and management. *PLoS ONE* **7**, e29232 (2012).

60. Zeppilli, D., Pusceddu, A., Trincardi, F. & Danovaro, R. Seafloor heterogeneity influences the biodiversity-ecosystem functioning relationships in the deep sea. *Sci. Rep.* **6**, 1–12 (2016).

61. de la Torriente, A. *et al*. Benthic habitat modelling and mapping as a conservation tool for marine protected areas: a seamount in the western Mediterranean. *Aquat. Conserv. Mar. Freshw. Ecosyst.* **29**, 732–750 (2019).

62. Gallardo, M., Macpherson, E., Tapia-Guerra, J. M., Asorey, C. M. & Sellanes, J. A new species of *Munida* Leach, 1820 (Crustacea: Decapoda: Anomura: Munididae) from seamounts of the Nazca-Desventuradas Marine Park. *PeerJ* https://doi.org/10.7717/peerj.10531 (2021).

63. Castilla, J. C. *Islas oceánicas chilenas: conocimiento científico y necesidades de investigación* (Ediciones Universidad Católica de Chile, 1987).

64. Bahamonde, N. San Félix y San Ambrosio, las islas llamadas Desventuradas 85–99 (1987).

65. Díaz-Díaz, O., Bone, D., Rodríguez, C. T. & Delgado-Blas, V. H. Poliquetos de Sudamérica. **Especial d**, 149 (2017).

66. Díaz-Díaz, O. F., Rozbaczylo, N., Sellanes, J. & Tapia-Guerra, J. M. A new species of *Eunice* Cuvier, 1817 (Polychaeta: Eunicidae) from the slope of the Desventuradas Islands and seamounts of the Nazca Ridge, southeastern Pacific Ocean. *A New Species Cuscus* **4860**, 211–226 (2020).

67. Kantor, Y. & Sysoev, A. Latiaxis (Babelomurex) naskensis, a new species of Coralliophilidae (Gastropoda) from South-Eastern Pacific. *Ruthenica* **2**, 163–167 (1992).

68. Sepulveda, J. I. Peces de las Islas Oceánicas Chilenas. In *Islas Oceánicas Chilenas: Conocimiento científico y necesidades de Investigaciones*. (ed. Castilla, J.) 225–246 (Ediciones Universidad Católica de Chile, 1987).

69. Mironov, A. & Detinova., N. Bottom fauna of the Nazca and Sala y Gomez ridges. *Plankton and benthos from the Nazca and Sala y Gomez Submarine Ridges* 269–278 (1990).

70. Lundsten, L. *et al*. Benthic invertebrate communities on three seamounts off southern and central California USA. *Mar. Ecol. Prog. Ser.* **374**, 23–32 (2009).

71. Rex, M. A. *et al*. Global bathymetric patterns of standing stock and body size in the deep-sea benthos. *Mar. Ecol. Prog. Ser.* **317**, 1–8 (2006).

72. QGIS.org. QGIS Geographic Information System. QGIS Association. Version 3.10. https://www.qgis.org (2020).
## Acknowledgements

We are indebted to the Captain and crew of R/V Cabo de Hornos of the Chilean Navy, and the scientific personnel participating in the CIMAR 22 cruise. Special thanks to Msc. Maria Valladares for helping us with environmental data and sampling onboard and Jorge Avilés from SCBUCN for his help during handling and curation of Agassiz trawl samples. We thank Dr. Nicolás Rozbaczylo, Dr. Oscar Díaz, Dr. Guillermo Guzmán, M. Sc. Juan Cañete, Dr. Cynthia Asorey and all the scientists who collaborated with the identification of the fauna of this study. Our thanks also go to OCEANA for providing the ROV used in the study and FONDEQUIP EQM 150109. Funding was provided by Grants CONA C22 16-09 to JS, AM, and EEE; FONDECYT #1181153 and 1180694 to JS and AM; ESMOI to JS, AM, EEE, MAG and JMT; UCN to MAG and JMT, Beca Magister Nacional CONICYT # 22190560 (now ANID); Conservation International to JMT and Beca de Postdoctorado Universidad Católica del Norte N°003 to MAG.

## Author contributions

J.M.T. processed and identified specimens, analyzed the data, and wrote the manuscript. J.S., A.M., E.E.E. and M.G. developed the proposal and secured funding for participation in the CIMAR22 cruise. A.M., E.E.E. and M.G. conducted fieldwork, assisted in specimen identifications and other lab work, and contributed to writing and editing the manuscript. J.S. and M.A.G. assisted in species identification, analyzed the data, contributed to writing and editing the manuscript. All authors reviewed the final version of the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

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

Correspondence and requests for materials should be addressed to J.S.

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

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

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

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

© The Author(s) 2021