Deep-Sea Habitats and Megafauna on the Slopes of the São Paulo Ridge, SW Atlantic

The São Paulo Ridge (SPR) is a 350 km-long linear geological feature located in the Continental Margin off Brazil (Latitude 28–29°S, Longitude 40–45°W). In 2013, the region was mapped during the SW Atlantic “Iata-Piuná” expedition and explored by a series of deep-sea dives of the manned submersible Shinkai 6500. A digital bathymetric model analyzed for seafloor morphology, delimited four major bathymetric sectors namely plateau, ridge crest, ridge escarpment and ridge foot. These sectors further enclosed 12 morphological features at smaller spatial scales (structural classes) including plains, valleys, peaks, terraces, and troughs. Video profiles across the depth gradient (4,219–2,644 m depths) revealed that the slopes of the SPR southern flank were gentle and terraced, mostly covered by biogenic sediments and interrupted by rocky cliffs/crests, dispersed outcrops and loose particles. The North Atlantic Deep Water (NADW) and Antarctic Bottom Water (AABW) overlaid at the escarpment along which they established colder (0.4–1.0°C; 4,200–3,400 m) and warmer (2.0–3.0°C; 3,400–2,600 m) habitats, respectively. Physical components were used to define seven seascape units in the ridge foot (2), escarpment (3), and plateau-ridge crest (2), where a total of 914 organisms of the epibenthic and benthopelagic megafauna were recorded. Over 70% of these records were sessile suspension feeders, including sponges (61.5%) and anthozoans (11.4%). Most taxonomic groups concentrated above 3,800 m, under the influence of NADW, where densities reached maximum values (mean 0.26 organisms.m–2; 0.024–0.027 organisms.m2 95% CI). Also, nearly half of megafauna records concentrated in patches delimited by the 3,800–3,300 m and 2,900–2,700 m isobaths. The deepest patch (3800–3300 m) coincided with the interface zone between AABW and NADW, where mixing processes create a density gradient. Evidences suggested that topography-related deep-water flow dynamics, and not substrate availability, drives benthic megafauna distribution at meso-habitat scale.


INTRODUCTION
In recent decades, marine sciences, industry and conservation initiatives have turned their attention toward the Southwest Atlantic. The region has a critical role in the Atlantic Meridian Overturning Circulation (AMOC), whose dynamics influence poleward heat flux and global climate (Garzoli and Matano, 2011;Frajka-Williams et al., 2019). Deep components of AMOC are strongly affected by distinctive seafloor topographies, whose origin and morphology derive from geological events established throughout the history of South Atlantic Ocean expansion (Bassetto et al., 2000;Ussami et al., 2012). In association with topographic features, valuable mineral deposits have been mapped and explored (Hein et al., 2013) and deep ecosystems and biodiversity have been increasingly described (e.g., Perez et al., 2012;Kitazato et al., 2017;Jovane et al., 2019).
Particularly relevant in this context is the east-west trending system of topographic features that include the Rio Grande Rise, Vema Channel, São Paulo Ridge and São Paulo Plateau (Figure 1). With depths spanning 600-5,000 m, these rises and troughs interpose and provide channels for the flow of the deepest water masses of the Atlantic Ocean: the North Atlantic Deep Water (NADW) and the Antarctic Bottom Water (AABW). The former flows southwards along the West Atlantic, compensating the northward circulation of surface, central and intermediate waters, and maintaining mass balance in the Atlantic (Garzoli and Matano, 2011;Frajka-Williams et al., 2019). Transported by the Deep Western Boundary Current between 1,500 and 3,000 m depths, NADW circulation is constrained westward by the topography of Brazil's continental slope and transversely oriented seamounts and ridges, most notably the Vitoria-Trindade Chain, Rio Grande Rise and São Paulo Ridge (Stramma and England, 1999;McDonagh et al., 2002). Below 3,000 m, the AABW flows northwards, along the Southwest Atlantic basin and into the North Atlantic, often via abyssal conduits, most notably the 5,000 m-deep Vema Channel (McDonagh et al., 2002;Morozov et al., 2010). These deep-water masses, with their distinctive physical and chemical properties, interact with the seafloor contributing to the establishment of variable sedimentation regimes, habitats and biological communities. Descriptions of these relationships in the region are generally scarce but have progressively increased often driven by initiatives addressing the needs for conservation and future sustainable mineral exploitation (e.g., Sumida et al., 2016;Hajdu et al., 2017;Perez et al., 2018;Montserrat et al., 2019).
Bathyal communities tend to change continuously along depth gradients essentially because of depth-correlated conditions such as pressure, temperature and dissolved oxygen (Carney, 2005). Drastic changes, however, may be driven by spatial discontinuities in physical and chemical conditions, food supply, sedimentary regime, bottom currents, and major topographic features (Rex and Etter, 2010). These effects may be conspicuous in the São Paulo Ridge (SPR) whose southern flank forms a steep 2,200 m depth gradient that extends transversely to the predominant flow of the SW Atlantic deep circulation (Alberoni et al., 2019). The region was subject to geological studies in the 1970s and 1980s (e.g., Gamboa and Kumar, 1977;Gamboa and Rabinowitz, 1981) but remained mostly unexplored in terms of habitat configuration and biodiversity. In 2013 the SPR was targeted by a global expedition "Quelle 2013 -Quest for the limits of life" led by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC), which searched for extreme deep-sea environments with the RV Yokosuka and the manned submersible Shinkai 6500  1 . An important finding of this exploration in the SPR was a deep whale fall with a well-established chemosynthetic community (Sumida et al., 2016;Cavalett et al., 2017;Shimabukuro et al., 2017;Shimabukuro and Sumida, 2019). In addition, the geological setting of the explored areas and a report of a fossil whale skull was provided by Ichisima et al. (2017). In this study we describe nonchemosynthetic benthic meso-and macro-habitats (sensu Greene et al., 2007) and megafauna along the SPR depth gradient and explore the effect of associated abiotic factors, including substrate, seafloor morphology and deep-water mass stratification. We further estimate large-scale distribution of benthic habitats in the SPR by analyzing bathymetry-derived terrain variables and classification of seafloor features. Benthic megafauna variability across different spatial scales will be discussed as a baseline and hypotheses for more detailed ecological studies in the future.

Study Area
The São Paulo Ridge (SPR) is a 350 km-long linear geological feature located in the Continental Margin off Brazil, between 28-29 • S and 40-45 • W ( Figure 1A). It is a component of the east-west trending Rio Grande Fracture Zone alignment that delineates the southern boundary of the São Paulo Plateau (Bassetto et al., 2000). As other aseismic ridges, the SPR is asymmetric with (a) a flat 2,000 m-deep northern flank, buried by sediments of the São Paulo Plateau, and (b) a steep southern flank (known as the São Paulo Escarpment), diving from 2,500 m-deep crests to a 4,200 m-deep foot (Gamboa and Rabinowitz, 1981;Alberoni et al., 2019). It is believed that SPR current morphology is largely derived from irregular submarine vulcanism that started in the Aptian period (∼120 m.a.), during the early South Atlantic Ocean expansion, when the SPR acted as a barrier obstructing the marine circulation toward the north and promoted the formation of a shallow environment where the deposition of a thick layer of evaporites resulted in the formation of the São Paulo Plateau (Gamboa and Rabinowitz, 1981;Bassetto et al., 2000;Ichisima et al., 2017). The southern flank interposes the northern flow of the Antarctic Bottom Water (AABW) generating countercurrents and eddies that mobilize sediments and form a trough delineating the deep contour of the ridge, also known as the São Paulo Channel (Figure 1B; Gamboa and Kumar, 1977;Alberoni et al., 2019). The NADW flows southward over the São Paulo Plateau, across the SPR crests and down over the escarpment where it overlays the AABW.

Analyzed Data
Data on SPR seafloor benthic habitats and megafauna were acquired during a research cruise conducted by the RV "Yokosuka" in 2013, under the "Iata -Piuná" consortium established between the Japan Agency for Marine-Earth Science and Technology (JAMSTEC), the Oceanographic Institute of University of São Paulo (IOUSP) and the Geological Survey of Brazil (CPRM) . The SPR area was explored between April 23 and 27 and included swath bathymetry and deep-sea dives of the manned submersible Shinkai 6500.
Bathymetric data were acquired during along-ridge transects with a 12 kHz hull-mounted Multibeam Echosounder System (MBES), integrated Differential Global Positioning System (DGPS) and compensated by the Inertial Measurements System (IMU). The data were filtered to remove spurious depth records and submitted to methods of linear interpolation to produce a Digital Bathymetric Model (DBM) with a cell size of 123 m. Video transects were produced during four deep-sea dives along a complete SPR southern flank depth profile, from 4,219 to 2,644 m depths (Figure 1 and Table 1).
Dives 6K1333 and 6K1336 explored the abyssal region at the foot of the ridge. Part of these dives were dedicated at the study of a whale-fall carcass environment, encountered at 4,204 m depth, previously described by Sumida et al. (2016) and excluded from this study. Dives 6K1334 and 6K1335 described the seafloor along the depth gradient of the São Paulo Escarpment. All dives involved approximately 4 h of activities near the seafloor, including photo/video recording, geological and biological sampling. Video data was acquired by two HD-TV color video cameras, both positioned at the bow, 1.7 m above the vehicle's bottom (Nakajima et al., 2014). Camera 1 angled obliquely 40 • toward the seafloor and recorded continuously a fixed area ahead the bow of the submersible (horizontal acceptance = 90 • , vertical acceptance = 57 • ). Camera 2 was mobile (pan -tilt) and was used for detailed observations of habitat features and megafauna species. Continuous information of date/time, depth (and altitude in meters) and the vehicle's heading (in degrees) were overlaid in the videos. Horizontal position (latitude, longitude) was estimated by SSBL (Super Short Base Line) method which required a transponder mounted on the submersible and an where α is the altitude of the camera, θ and ω are the camera's horizontal and vertical acceptance angles, respectively, and δ is the angle of the camera from vertical (Nakajima et al., 2014). Measurable uncertainty in area estimations derived from variability in α along video transects. Bootstrap 95% confidence intervals of MLW were calculated for video transects or segments of it, using the quantile method. Megabenthos specimens were collected by the submersible's manipulators and slurp-gun, photographed on board and stored in ethanol (75%) and formalin for posterior identification by specialists in the National Museum (Federal University of Rio de Janeiro) where they were cataloged.
Videos, sample records and oceanographic data and metadata are available in "DARWIN -Data and Sample Research System for Whole Cruise Information in JAMSTEC 2 ."

Seafloor Morphometric Analysis
The DBM was analyzed for seafloor morphology and segmentation (Brown et al., 2011). This process included the transformation of bathymetry data into secondary-derived layers (terrain variables) using the algorithm package Benthic Terrain Modeler (BTM) (Walbridge et al., 2018) contained in ArcGIS Desktop R 10.2.2. Variables with the greatest potential for the description of benthic habitats at the corresponding spatial scales were chosen, including surface gradients (slope, aspect), relative depth (Bathymetric Position Index -BPI) and surface rugosity (Vector Roughness Measurement -VRM) (Wilson et al., 2007;Walbridge et al., 2018). BPI is a neighborhood analysis function of the mean depth around each cell in the DBM. Positive values correspond to features and regions that are higher than the surrounding area which characterize ridges; negatives values would represent depressions on the seafloor. 2 http://www.godac.jamstec.go.jp/darwin/ Values equal or near zero are either flat or constant slope area (Walbridge et al., 2018).
The procedures for seafloor segmentation (Figure 2) started with overlaying 23 bathymetric profiles on the DBM, transversal to the bathymetric gradient of the São Paulo Escarpment (Supplementary Figure 1) for visual interpretation. Along each profile, the variables BPI-broad (annulus = 300 pixels ∼ 36,900 m), BPI-fine (annulus = 60 pixels ∼ 7,380 m) and slope were extracted and analyzed individually in order to identify Zonal and Structural classes (sensu Erdey-Heydorn, 2008). Zonal classification described broad scale surficial characteristic of the seafloor. Classes representing the main bathymetric components were visually defined, and had their boundaries delimited by depth and values of BPI-broad and slope (i.e., averages across the 23 transversal bathymetric profiles). Within each Zonal class, surficial characteristics of the seafloor were then described to a greater detail (Structural classification) using terrain variables BPI-fine and slope. Structural classes were again visually defined and delimited by maximum and minimum values of depth, BPIfine and slope values (i.e., averages across the 23 bathymetric profiles). A classification table was then built, summarizing lower and upper limits of depth, slope, BPI-broad and BPI-fine for all visually defined Zonal and Structural classes (Supplementary Table 1). The classification table was then used as input data for semi-automated classification of new adjusted classes by BTM's Classify Benthic Terrain tool. Such classification procedure (followed by spatial representation of the resulting classes) was performed several times until the spatial segmentation was considered satisfactory (Figure 2). Each model run was conducted after manual adjustments of the terrain variable values delimiting structural classes, and/or the creation of new intermediate classes.

Seascape Classification
Each video produced by camera 1 was initially observed to record changes in depth, altitude, submersible activities and image visibility. Also, seafloor features relevant for habitat classification (following Greene et al., 1999) were selected, namely: substrate texture, particle sizes, relief and others. These features were organized in a scale of "substrate types, " represented by capital letters, where: R = rocky ridge with rugged surface; F = flat rocky pavement with plain or granulated surface; B = loose boulders (>25.5 cm); C = cobbles (>6.5 cm and <25.5 cm); P = pebbles (>2 cm and <6.5 cm); G = gravel (>4 mm and <2 cm); U = unconsolidated substrate varying from fine mud to coarse sand.
A second analysis of these videos included only segments when the submersible was moving ahead. During these segments, substrate types were attributed to 1-min video intervals using a two capital letters system; the first referring to the substrate type covering more than 50% of the visible seafloor, and the second to the substrate covering between 30 and 50%. For example, a segment where seafloor was mostly unconsolidated (e.g., covered by biogenic sediment) with variable amounts of scattered pebbles was codded as UP. Also, when only one type of substrate was visible during the observed interval the correspondent letter was duplicated (e.g., UU = seafloor completely covered by sediments). These combinations of two substrate types were defined as "bottom types" (Greene et al., 2007;Tissot et al., 2007;Tissot, 2008).
Seascapes units were mainly defined by one or more dominant bottom types regularly combined along a continuous segment of the seafloor. Adjacent seascapes were delimited by an abrupt change in the dominant bottom type (or bottom type combination) observed for more than 10 s along the video track, indicating the end of a seascape unit and the beginning of a new one. Bottom types that contrasted with the dominant ones but did not persist long enough along the video track (e.g., a patch), were not considered a new seascape unit, but as part of the current seascape substrate variability. Descriptions of seascape units were complemented by "modifying" elements (e.g., currents, bioturbation), relief and slope (Greene et al., 1999), as well as by the influence of AABW and NADW on the seafloor. This influence was estimated by temperature and salinity data recorded continuously during the dive track by the submersible's CTD. Mixing percentages of these water masses during the dives were calculated according to Mamayev (1975).

Megafauna Diversity
Videos produced by cameras 1 and 2 were analyzed for visible megafauna. Life forms smaller than approximately 5 cm and/or visible in images taken above 2-m altitude were not included in the analysis. Otherwise, videos were stopped at each sighting and recorded the organism type ("morphotype") along with associated information, including date, time, depth (m), altitude (m), and heading (in degrees). Morphotypes were classified in higher taxa (Phylum, Class, Orders) and their consistency was double checked by repeating the analyses of videos produced by both camera 1 and 2. When no morphotypes could be confidently assigned to a given observed organism after these repeated analyses, it accounted for higher taxa quantification only. Some morphotype identifications to family, genus and species level were possible through collaboration with deep-sea fish and invertebrate taxonomists, and with the aid of deep-sea fauna image databases (e.g., OER's Benthic Deepwater Animal Identification Guide 3 and others).
Megafauna spatial distribution and abundance was analyzed by representing density of total recorded megafauna and/or taxonomic groups as a function of the depth strata and seascape units. Density was calculated by dividing organism numbers by the mean estimated areas covered by video transects (expressed as individuals. m −2 ), and their 95%CI limits.

Seafloor Segmentation
The SPR is depicted as the southern margin of the São Paulo Plateau, which extends to the north as 3,000 m-deep plane areas (Figure 1). Along the border of the plateau, prominent ridge crests rise to 2,200 m depths, over 1,000 m above the plateau level. These crests are separated by plane areas and channels that extend to the edge of the plateau, possibly characterizing paths of deepwater flow. The margin of the plateau is continuously bordered by the São Paulo Escarpment that forms a south -southeast facing wall, which is (a) higher/steeper in sectors adjacent to the plateau crests (Supplementary Figure 2, e.g., profiles 3, 11, and 19) and (b) lower/gentler in sectors adjacent to plane areas in between the ridge crests (Supplementary Figure 1, e.g., profiles 8, 15, and 22). The deep edge of the escarpment is connected to the 4,100-4,300 m-deep São Paulo Channel, a trough extending continuously along the SPR. To the south of this channel, the seafloor rises into a 4,000 m-deep plain area that is part of the Santa Catarina Plateau (Figure 1).
The spatial representation of BPI-broad delineated the areas comprised by broad bathymetric components, including the São Paulo Plateau, plateau crests, São Paulo Escarpment, São Paulo Channel and Santa Catarina Plateau (Figure 3). Secondary seafloor structures were evidenced from the spatial analysis of BPI-fine and included valleys and channels on the São Paulo Plateau, as well as discontinuities of the São Paulo Channel. Slope, aspect (easterness and northerness) and roughness (VRM) provided refined spatial representations of (a) channels and valleys bordering the plateau crests, and crossing the margin of the São Paulo Plateau, and (b) the escarpment, depicted as a steep, southern-faced roughed terrain (Figure 3). These were seafloor topography subject to classification by the morphometric analysis (see below).
The seafloor classification procedure resulted in 12 Structural Classes (Figure 4) enclosed within four Zonal Classes. Descriptors (upper and lower limits of depth, terrain variables and slope) and names attributed to each on of them are presented in Supplementary Table 1. The plateau (Zonal Class II) comprises plains (Class II.5), valleys (Class II.6), and gentle slopes (Class II.7) which were shallower than 3,700 m and with less than 5 • slope. The SPR crest (Class I) is composed of seamount-like summits of high (Class I.2), moderate (Class I.3), and low (Class I.4) altitudes, some of the former topped by 2,240 m-deep peaks (Class I.1). Bordering the edge of the plateau,

Seascapes
The deep-sea dives crossed four zonal classes: ridge foot, ridge escarpment, plateau and ridge crest ( Figure 4C). Seafloor texture was dominated by biogenic sediments usually mixed with scattered rocky particles (UP, UC), outcrops and crests (UF, UR) (Supplementary Table 2). These bottom types were recorded for over 56% of the observation time, followed by bedrock outcrops (RU + RC = 33.5%). Areas entirely covered by rocky substrata (e.g., RC, RR, RP, FP, PC) were recorded in 7.1% of the observed time. Most of the area observed (>90%) corresponded to gentle slopes (5-30 • ). Steep slopes (>30 • ) occurred in 2% of the observed area and only in the ridge escarpment.
Texture classification, slope, water masses and modifying elements allowed the differentiation of seven seascape units ( Table 2). The "Abyssal Mud Field" (AMF), was characterized by an undulated sediment surface covering the São Paulo Channel seafloor, with isolated patches of pebbles and bedrock outcrops (UP, UR). This seascape dominated images produced by dive 6K1333 that explored a linear track along the foot of the ridge, in eastward direction over the 4,000-4,100 m isobaths. The effect of AABW flux was noticeable in this seascape, usually producing regular sand waves ( Table 2). Adjacent do AMF, the "Debris Field" (DF) extended along the deep border of the escarpment generally characterized by a mixture of cobbles, pebbles, and boulders, interspersed with biogenic sediments (UC and UP). These bottom types, recorded in 64.5% of the total observation time of this seascape (Supplementary Table 2), were often associated with landslide signs and coated by a thin sediment cover ( Figure 5A).
Moving into the lower section of the ridge escarpment (above 4,100 m depths) the seascape was dominated by a steep and roughed bedrock surface, the "Bedrock Cliff " (BC) ( Table 2). Rocky ridges sometimes covered by piles of loose cobbles (RU and RC) dominated the seafloor (64.5% of the observation time, Supplementary Table 2) ( Figure 5B). Small amount of sediments accumulated in cracks and crevices, and thinly covered bedrock particles. This seascape unit was interrupted at approximately 4,070-3,998 m depths by a narrow terrace formation, the "Granular Flat" (GF) ( Table 2 and Figure 5C), where flat bedrock pavements were usually covered by pebbles and variable amounts of sediments (FP + FU = 71.4% of the observation time, Supplementary Table 2). A positive (convex) relief characterized the upper section of the escarpment (above 3,366 m depths), where the seascape "Bedrock Crest" (BCR) was defined by a mixed substrate formed by bedrock outcrops and sediment ponds (RU) (Table 2 and Figure 5D). This seascape unit was also observed in the ridge crest, between 2,578 and 2,657 m depths (Figure 5).  At the ridge plateau and the adjacent ridge crest, flat sedimented areas (UU), occasionally with sparse patches of pebbles and cobbles (UP, UC), formed a seascape unit called "Mud Terrain" (MT) ( Table 2). Sediment surface was generally smooth, but often modified by animal tracks (Figure 5E). The ridge plateau also included a slightly convex "Granular Crest" (GCR) (Table 2 and Figure 5F) covered by pebbles, cobbles and soft sediments (RU + UP = 82.2% of observation time Supplementary Table 2).

Megafauna
A total of 914 organisms of the epibenthic and benthopelagic megafauna were recorded during 700 min of SPR seafloor observation. Over 70% of the records were sessile suspension feeders, including sponges (62.8%) and anthozoans (11.6%). These were followed by shrimp-like crustaceans (11.3%) and fish (Actinoperygii, 8.7%) ( Table 3). A list of identified taxa is presented in Table 4.
Megafauna records were usually sparse (<0.01 individuals.m −2 ) along the SPR depth gradient, but two patches of moderate concentrations (0.07-0.30 individuals.m −2 ) mostly of suspension feeders, occurred at 3,800-3,300 m, and 2,900-2,700 m depth intervals (Figure 6). The former patch occurred in the ridge escarpment and the BC and BCR seascapes. The latter patch occurred at the ridge crest and within the reach of the GCR, MT and BCR seascapes. Both depth zones were 2 | Classification of seafloor explored during deep-sea dives in the São Paulo Ridge, SW Atlantic into seascape units according with Greene et al. (1999).

Seascape unit
Class ( covered by hard substrata (e.g., outcrops or loose particles), but also included a variable coverage of sediments. Most taxonomic groups were observed above 2,900 m and were associated with mixed substratum in the BCR seascape (Table 3). Nearly 60% of all records of Porifera occurred in this seascape being largely represented by two species: a pedunculated sponge, Family Dendoricellidae (cf Pyloderma sp., Demospongiae), and Poliopogon amadou (Hexactinellida) (Figure 7). Corals (Anthipataria and Alcyonacea) also tended to occur deeper, between 3,300 and 3,800 m, on the ridge escarpment (Seascape BC). 74.3 and 74.0% of crustaceans and echinoderms, respectively, were recorded below 2,900 m ( Table 3). The latter was dominated by Holothuroidea, Order Elasipodida, including Psychropotes semperiana, Benthodites sp., and Enypniastes sp. These were observed in BCR usually occurring on sediment ponds. Nearly 44% of fish records concentrated in the 2,700-2,900 m depth stratum (Table 3). One   species, Acanthonus armatus (Ophidiiformes, Ophidiidae) was particularly abundant in association with the GCR seascape ( Figure 7D). Epibenthic and benthopelagic megafauna was recorded within a total area estimated in 35609.4 m 2 (34076.6-37236.8 m 2 95% CI). Mean density was 0.026 individuals.m −2 (0.024-0.027 individuals.m 2 95% CI) decreasing progressively with depth (0.130-0.005 individuals.m −2 , Table 5). The sectors above 3,300 m, including GCR, MT and, BCR seascapes, exhibited more elevated concentrations of megafauna, all of them under a prevailing influence of NADW (usually > 80%, Table 5). Upper sectors of BCR (<2,900 m) exhibited the highest fauna concentrations (0.328 individuals.m −2 -0.300-0.360 95%CI) largely dominated by sponges cf Pyloderma sp. and P. amadou (Table 5 and Figure 7C). Estimated densities of taxonomic groups are presented in Supplementary Table 3.

DISCUSSION
Deep-sea habitats and megafauna biodiversity were described in a limited area across the SPR depth gradient and related to variability of physical factors. Comparable to other ridges in the Atlantic, the slopes of the SPR southern flank were gentle and terraced, mostly covered by biogenic sediments and interrupted by rocky cliffs/crests, dispersed outcrops and loose particles (e.g., the Mid-Atlantic ridge: Priede et al., 2013;Niedzielski et al., 2013;Bell et al., 2016;Alt et al., 2019). Benthic and benthopelagic fauna were generally scarce but spatial variations were noticeable, and potentially driven by depth (and depth covariates), water masses and seascapes zonal distribution. Insufficient sampling precluded attributing causal factors to faunal distribution observed at a meso-and macrohabitat scales. It is possible, however, to assume that observed patterns result from the hierarchical effects of structuring factors operating at different spatial scales (Levin  , 2001;Williams et al., 2010). Allied to bathymetry-derived seafloor morphology segmentation (Brown et al., 2011), reported habitat and megafauna distribution data provided elements to identify scales of spatial variation and to draw hypotheses of potential driving factors. Biodiversity patterns in the deep SW Atlantic are partly determined by regional (>1,000 km) biogeochemical processes at surface waters. Most of the region down to 30 • S lays under the influence the South Atlantic Gyral biogeochemical province (Longhurst, 1995) where surface primary production is generally low due to a stable water column structure and rapid vertical remineralization, resulting in minimal POC export fluxes to the seafloor (Mouw et al., 2016). Schlitzer et al. (2003) estimated POC flux values as low as 0.05-0.1 mol C m −2 yr −1 in the central South Atlantic, 10-30 times lower than values estimated along the continental margins of West Africa and South America off Brazil. This reduced POC flux reaching the seafloor would tend to support reduced megafauna densities in the SPR (Sibuet et al., 1989;Smith et al., 2008;Wei et al., 2010), as reported in this study (mean 0.026 individuals.m −2 ; range 0.004-0.350 individuals.m −2 ). When compared to megafauna densities reported from bathyal areas (2,000-3,500 m) of the North Atlantic (compiled in Levin and Gooday, 2003) (Bell et al., 2016;Alt et al., 2019). Densities recorded in comparable depths at the SPR (<2,900 m depths; 0.001-0.33 individuals.m −2 ) were approximately 100 times lower than those recorded to the north of the CGFZ but approximated those recorded to the south.
The regional effect of POC flux on megafauna distribution and abundance is altered by prominent topographic features, such as the SPR and Rio the Grande Rise (RGR), which generate variability in a provincial-scale (∼100-1,000 km). These features disrupt the SW Atlantic basin general morphology (a) producing abrupt discontinuities in depth and slope, (b) modifying the countercurrent flow of the NADW and AABW, and (c) exposing areas of rocky seafloor. The SPR forms a long and linear slope that extends transversely to the South American continental margin north -south orientation and separates a 2,000 m-deep sedimentary plateau (São Paulo), that extends to the north, from a 4,000 m-deep sedimentary-tectonic plateau (Santa Catarina) that extends to the south. This morphology was determined by the tectonic evolution of the South American continental margin and subsequent sedimentary processes, mostly associated with along-slope action of bottom currents and the downslope mass-transport and sediment fluxes (Alberoni et al., 2019). These are processes associated with the flow of NADW and AABW that overlay each other at the slopes of the SPR escarpment (∼3,400 m depth), exposing benthic habitats to downslope changes in physical and chemical conditions. In the explored area, temperature increased by nearly 3.0 • C (∼0.4-3.2 C) over a 1,500 m depth range from the ridge foot (4,200 m) to the ridge crest (2,700 m). Where both water masses mix, at the 3,500-3,300 m depth interval, temperature increased by 1.2 • C delimiting lower colder (0.4-1.0 • C) and upper warmer (2.0-3.0 • C) habitats (Supplementary Figure 2). In addition, habitats below this depth interval, under the influence of AABW, also tend to be less oxygenated and less saturated with CaCO 3 than those in the upper slope where the influence of NADW predominates (Chung et al., 2003;Rijkenberg et al., 2014). These are contrasting conditions that may partly explain the nearly 10-fold difference in the mean megafauna densities estimated in the SPR depth gradient, below and above 3,400 m depths (Table 5). Furthermore, the SPR escarpment interposes the SW Atlantic depth horizons for CaCO 3 Aragonite saturation ( Arg = 1), Aragonite compensation (ACH), and Calcite saturation ( Cal = 1) at ∼2,600, ∼3,400, and ∼4,000 m depths, respectively (Melguen and Thiede, 1974;Thunell, 1982;Chung et al., 2003). These conditions would limit growth of scleractinian cold-water corals in the SPR but not alcyonaceans (Octocorals), which build calcitic skeletons (Roberts et al., 2009). In the explored area such conditions are consistent with the general absence of scleractinian corals and the existing records of bamboo corals (Alcyonacea, Isididae), whose distribution patterns tend to be driven by Calcite saturation levels (Yesson et al., 2012). Finally, topography, depth, temperature, salinity, dissolved O 2 and POC fluxes were all environmental proxies used to define bathyal and abyssal biogeographic provinces (Watling et al., 2013). The SPR is the boundary between a lower bathyal (2,000-3,500 m depths) South Atlantic province and two adjacent abyssal provinces (3,500-6,500 m), the Argentine Basin and the Brazil Basin provinces (Watling et al., 2013). The area explored across the SPR southern flank included the South Atlantic and the Argentine Basin provinces. The drastic changes on megafauna abundance and composition above 3,400 m depths locally support these provinces and the proposed depth boundary.
A variety of megahabitats (∼1-100 km) were inferred by seafloor morphology segmentation in the ridge crest, escarpment, plateau and ridge foot zones. Some seafloor structural classes were ground truthed along the depth profile explored by video cameras, despite the different spatial resolution of both methods. In the ridge escarpment, the structural classes "escarp" (Class III.8) and "terrace" (Class III.9) differentiated a lower steeper zone from an upper gently sloping zone. This differentiation was evident in the video seafloor analysis, as a steep bedrock cliff changed into a gently sloping bedrock crest seascapes (Figure 5), also with a significant increase on the occurrence of sessile suspension feeders (e.g., sponges) and deposit feeders (e.g., holothurians) ( Table 3). Similar correspondences were also noted at the ridge foot, plateau and ridge crest, suggesting that bathymetry-derived seafloor morphological units were biologically relevant and could express habitat heterogeneity along the SPR. Allied to observations on fauna composition and distribution (see below), biotopes could be predicted and mapped (Brown et al., 2011). However, considerable additional sampling effort would be needed for that purpose (e.g., Robert et al., 2015;Anderson et al., 2016), which was generally beyond the scope of this study .
Nearly half of megafauna records (48%) along the explored depth gradient concentrated in 700-900 m-long patches delimited by the 3,800-3,300 m and 2,900-2,700 m isobaths. The lower zone was covered mostly by rough rocky surfaces whereas the upper zone was characterized by mixed substrates with a predominance of sediments. Despite such differences in substrate composition, megafauna observed in both patches were dominated by sessile suspension feeders (cnidarians and sponges, 73-74% of recorded organisms), followed by benthopelagic organisms (swimming shrimps, fish and cephalopods, 19%) and the less frequent soft bottom dwellers (mostly echinoderms, 6-8%). Suspension feeders did not seem to be limited by hard bottom availability and were commonly recorded even when only a few loose particles (cobbles, pebbles, small outcrops) were available interspersed with dominant sediment substrate. Soft bottom dwellers, on the other hand, were generally scarce despite largely available sedimented areas. Whereas these organisms may be limited by low recruitment rates in the area, records of lebensspuren on the sediment surface could indicate that part of this fauna could be buried in the sediments and not visible in the images. These observations suggest that megafauna spatial distribution at a mesohabitat scale (∼10 s of meters to km) in the explored area of the SPR was (a) less affected by substrate availability, as also reported by studies on the flanks of the Mid-Atlantic Ridge (Bell et al., 2016;Alt et al., 2019), and (b) driven by the availability of suspended food particles, as generated by topography-related current flow patterns over the seafloor.
Both patches of megafauna occurred in the vicinity of crests formed at the edge of terraces (Figure 6), where flow dynamics of the NADW may be particularly favorable for sessile suspension feeders (Genin et al., 1986). As a dense water mass flowing over the São Paulo Plateau, the NADW develops a bottom friction transport perpendicular to the depth contour, which, when overlaying a topographic depression such as the SPR escarpment, would drive the NADW flow inside this depression (Wahlin, 2002). This downslope flow, altered by abrupt changes in slope as observed in the edge of the terraces, could generate areas of vorticity and suspended particle concentrations. In the upper depth zone, evidence of such current action included sand waves observed in sediment cover and movements of the highly dense and flexible sponges Pyloderma. Downsloping NADW will encounter AABW opposing flow at 3,400 m depth, within the depth zone of the lower patch of sessile suspension feeders. This water mass interface formed a temperature and salinity gradient as well as a density variation which may characterize a physical boundary condition often associated with POM enrichment (Dullo et al., 2008). Additionally, this density surface (Supplementary Figure 2) may also be associated with turbulent mixing (Zhao and Thrunherr, 2017) and along-ridge currents as derived from changes in the AABW transport vorticity as it collides with the SPR southern flank (topographic steering, White and Dorschel, 2010). If these are physical processes occurring to some extent between 3,800 and 3,300 m depths, resuspension or acceleration of advected food particles could favor local settlement and growth of suspension feeders at this depth range. The NADW and AABW interactions extend along the SPR potentially acting at a provincial scale. Yet, because topography is determinant in such physical processes, its effect on megafauna distribution may vary along the ridge, driving variability at smaller spatial scales. In fact, channels and valleys were features modifying the ridge crest, plateau and escarpment overall morphology, potentially altering water flow and sustaining a variety of mesohabitats and patches of suspension feeding fauna.
An Antarctic minke whale skeleton fall was a major driver of megafauna distribution at macrohabitat scale in the debris field seascape (4,204 m depth, Sumida et al., 2016). A series of dispersed vertebrae and intervertebral disks spread over approximately 2 m over the seafloor, comprised over 40 species some of them occurring at high densities in the bone and sediment area (2-70 individuals.m −2 ). Except from one crustacean (gen. Munidopsis), species recorded at and around the whale fall were adapted to organic fall environments and not recorded anywhere else in the SPR explored area.
Most taxa occurrences were recorded from single or very few observations. The exceptions were the sponges P. amadou and cf. Pyloderma sp., and the benthopelagic fish A. armatus, relatively abundant in particular segments of the SPR depth gradient. Despite its uncertain identification, the sponge cf. Pyloderma formed the densest epifauna concentrations in the explored area of the SPR. P. amadou is the only species of this genus to occur in the Atlantic Ocean, reported in dense patches (up to 5 individuals.m −2 ) in the Great Meteor seamount (29 • 30 N; 28 • 17 W) between 2,675 and 2,765 m depths (Xavier et al., 2015); and in the Tropic Seamount (23 • 55 N; 20 • 45 W) where it occurred from 1,960 to 3,660 m depths (Ramiro-Sánchez et al., 2019). In the SPR the species was observed in much lower concentrations (∼0.03-0.06 individuals.m −2 ) around 3,000 m depth, attached to outcrops and loose rocky particles. A. armatus occurs at bathyal and abyssal depths (1,500-4,415 m) in all tropical and subtropical oceans, being particularly abundant in the western Atlantic (Nielsen et al., 1999). The species has been recorded in Brazil's continental margin off Bahia, from 1,171 to 1,929 m (Mincarone et al., 2008) and in the Caribbean between 1,500 and 4,150 m (Polanco et al., 2019). A comparable exploration in the Rio Grande Rise, revealed nearly 3× more fish morphotypes (30) during half the observation time (Perez et al., 2018), potentially because these dives explored much shallower area (1,233-600 m depth). Bathypterois, Spectrunculus, and Aldrovandia were the only genera found in both the SPR and the Rio Grande Rise. In the latter, they were observed only in deepest sectors (1,200-900 m).
The SPR has been characterized as an important area of environmental transitions in the SW Atlantic basin, mostly driven by the depth gradient of its southern flank, where the two main deep-water masses of the Atlantic dynamically interact. Important along-slope chemical and physical gradients are established along the ridge extension, but seafloor morphology, substrate type and topography-driven current flow processes may generate habitat heterogeneity at varying spatial scales, all relevant do deep fauna distribution. Unprecedented observations along the SPR depth profile have generally indicated that megafauna, although scarce, may respond to such drivers and vary considerably in the mesohabitats established along the ridge, justifying future explorations and studies designed to test the effect of the hypothesized drivers (e.g., Bell et al., 2016;Alt et al., 2019).

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article including videos, sample records and oceanographic data and metadata are available in 'DARWIN -Data and Sample Research System for Whole Cruise Information in JAMSTEC' (http://www.godac. jamstec.go.jp/darwin/).

AUTHOR CONTRIBUTIONS
JP oversaw all aspects of this research, including specimen, sample, data collection, and analysis. All authors conducted the research, analyzed the data, and contributed to the manuscript.

FUNDING
Funding of Brazilian scientists in the "Iata-Piuná" cruise was provided by a grant from CAPES (Program CAPES -JSPS, AUXPE-JSPS-0059-2013, Ministry of Education, Brazil). The senior author was supported by a CPNq productivity fellowship (Process 307992/2019-5). This study was within the umbrella of the National Institute of Science and Technology -Integrated Oceanography Centre (INCT -Mar COI, CNPq).

ACKNOWLEDGMENTS
We thank Brazilian and Japanese governments, and members of JAMSTEC, IOUSP, and CPRM, whose efforts allowed this unprecedented study in the SW Atlantic, as part of JAMSTEC's "Quelle" expedition. We owe the crews of the RV Yokosuka and the submersible Shinkai 6500 the acquisition of all analyzed data and samples. Eduardo Hajdu, Renato Ventura (Museu Nacional -UFRJ), Thayse Fonseca, and Richard Schwarz (UNIVALI) provided invaluable help with the process of identification of fauna from biological samples and video images. Angelica Maffini Mastella, Carine Eccel, and Anna Caroline Silva de Andrade contributed with video analysis procedures.
Katz Fujikura (JAMSTEC) kindly made available high-resolution bathymetry data.