Distribution and diel vertical migration of mesopelagic fishes in the Southern Sargasso Sea — observations through hydroacoustics and stratified catches

Vertical distribution patterns and relative abundance of mesopelagic fish species and other major taxonomic groups were investigated through vertically stratified trawl sampling and hydroacoustic analyses along the subtropical convergence zone from 52° W to 70° W in the oligotrophic Sargasso Sea. Persistent stationary layers and several migrating components of different scattering characteristics were detected. The results reveal varying vertical migration patterns, including different times of onset of diel vertical migration in different depths and a migrant pathway emerging daily from the lower deep scattering layer (DSL) at dusk and migrating through the upper DSL without affecting its composition. Fish species identification was made based on morphological characteristics and confirmed by genetic barcoding analyses of subsamples. In total, 5022 fish specimens from 27 families, 62 genera and 70 species were caught. In terms of relative abundance (A) and biomass (M), catches were dominated by species of the families Myctophidae (A=59.1%, M=47.4% of total fish catch) and Melamphaidae (A=22.5%, M=17.1%). Myctophidae and Stomiidae were the most species-rich families with 31 and 12 species, respectively. Catches at the two easternmost stations were dominated by Scopelogadus mizolepis and Nannobrachium cuprarium, while Bolinichthys photothorax and Ceratoscopelus warmingii were the most abundant species in catches from the two westernmost stations. This study provides insights into distribution and vertical migration behaviour of mesopelagic fish in the Sargasso Sea and adds to our understanding of the mesopelagic community in this large oceanic area.


Introduction
Mesopelagic fish, living in depths of 200 to 1000 m, are dominating the fish biomass in large parts of the world's oceans. A recent estimate of their global biomass was even a magnitude higher than the former estimate of approximately 1,000 million tons (Lam and Pauly 2005;Irigoien et al. 2014). Through extended diel vertical migrations (DVM) into shallower layers predominantly for feeding, the sheer biomass of mesopelagic fish has large influences on trophic connectivity and the carbon flux in the oceans (Saba et al. 2021). While the importance of pelagic fish communities for the nutrient flow in the oceans is widely acknowledged (e.g. Davison et al. 2013;Anderson et al. 2019;Saba et al. 2021), our understanding and knowledge of the vertical ecology of mesopelagic organisms in relation to their physico-chemical environment is still limited (St. John et al. 2016;Olivar et al. 2017;Proud et al. 2019;Romero-Romero et al. 2019) and studies on globally available data sets have demonstrated profound differences in the depth of day-time sound scattering layers between productive and oligotrophic regions (Bianchi and Mislan 2016;Klevjer et al. 2016).
From an oceanographic perspective, the Sargasso Sea is a particularly relevant transitional area in the western part of the North Atlantic. It is situated within the North Atlantic Subtropical Gyre and surrounded by ocean currents, which can trap the water in the core of the Sargasso Sea for decades (Maximenko et al. 2012). This large marine ecoregion is influenced by various hydrographical features, including complex patterns of thermal fronts, mesoscale eddies, advective transport of water masses and seasonal convective overturn, which all have an influence on the biota and lead to a pronounced spatio-temporal variability in planktonic productivity (McGillicuddy et al. 1998;Hansell and Carlson 2001;Palter et al. 2005;Eden et al. 2009;Hellenbrecht et al. 2019). Based on the prevailing low nutrient levels, the Sargasso Sea is generally classified as oligotrophic. However, in spite of this fact, it has a high net primary production rate per unit area (Steinberg et al. 2001;Laffoley et al. 2011).
Besides environmental factors like primary production, oxygen concentration and light intensity (Irigoien et al. 2014;Klevjer et al. 2016;Aksnes et al. 2017), the distribution of pelagic fish in the Sargasso Sea may also be influenced by distinct temperature fronts in the upper 150 m of the ocean. These fronts may form zoogeographic boundaries for certain mesopelagic fish (Backus et al. 1969) and cephalopod species (Lischka et al. 2017) caused by temperature-related changes in e.g. productivity or stratification. Changes in mesopelagic fish species composition across fronts have been reported, for example at frontal systems in the southern California Current region, demonstrating incongruent patterns in the spatial distribution of migrators and non-migrators (Netburn and Koslow 2018). In the Sargasso Sea, these fronts are caused by the confluence of warm and cold-water masses in the subtropical convergence zone (STCZ) (Katz 1969). It has also been demonstrated that distribution and transport of midwater fish species can be driven by eddies (Olson and Backus 1985;Godø et al. 2012) including also the transport into or out of the Sargasso Sea (Craddock et al. 1992).
A large number of studies describe vertical distribution patterns of mesopelagic fish (e.g. Sutton 2013;Klevjer et al. 2016;Olivar et al. 2016;Proud et al. 2017;Sutton et al. 2017). The distribution of mesopelagic fish species in deep scattering layers (DSL) and the species composition of those layers have been identified on a regional scale in the temperate northeastern Atlantic and in the subtropical and tropical western Atlantic and adjacent waters (Roe and Badcock 1984;Wienerroither et al. 2009;Peña et al. 2014;Ariza et al. 2016;D'Elia et al. 2016). Abundance and distribution of mesopelagic fish have also been investigated in the Sargasso Sea (Backus et al. 1969;Backus et al. 1970;Jahn and Backus 1976;Backus and Craddock 1977;Gartner Jr et al. 1989;Sutton et al. 2010;Ayala et al. 2016). However, as for most ocean areas, current data on mesopelagics from this remote region are scarce.
To further increase our knowledge about the distribution and abundance of mesopelagic fish species in the Sargasso Sea, a hydroacoustic analysis of scattering layers and a series of depth stratified mesopelagic trawl samplings were conducted during an extensive multipurpose pelagic fishery survey along the subtropical convergence zone from 52°W to 70°W . Stratified catches were compared with hydroacoustic data, in order to assess the vertical migration behaviour of mesopelagic fish and to characterize the composition of deep scattering layer, providing insights into the horizontal and vertical distribution of fish species over this large oceanic area.

Fishing gear and sampling methodology
Sampling took place from March 20 to April 6, 2014, during an extensive multipurpose fishery survey (WH373) on the German R/V Walther Herwig III in the Sargasso Sea with a pelagic trawl (Engel Netze, Bremerhaven, Germany). The net had an opening width of 30 m, a height of 20 m, a length of 145 m, and mesh sizes (knot to knot) from 90 cm decreasing stepwise to 40, 20, 10, 5, 4, 3, 2 cm, with a 1.5 cm mesh in the 27 m long codend. It was equipped with a multi-closing system that enabled fishing in three defined depth strata by a time-controlled closing mechanism (multisampler).
Trawl duration at all stations was 225 min, with each of the 3 nets open for 75 min, apart from station 233, where net 1 was opened for 90 min and nets 2 and 3 for 105 min each (Table 1). All hauls were towed obliquely from shallower to deeper layers. Towing speed through water was 3 knots at all stations. The multisampler codends were cleaned from remaining fish after each haul.

Hydroacoustic data
Hydroacoustic data were collected continuously with a Simrad EK60 scientific echosounder operated at 18, 38, 120 and 200 kHz with hull-mounted transducers at a ship speed of approximately 12 knots during transit and 3 knots during fishing operations. Only data from the 18 kHz and 38 kHz transducer were used for further analyses. Beam widths of the transducers were 11°(18 kHz) and 7°(38 kHz), respectively, and both were operated at 2000 W. The pulse length was 1.024 ms with the ping rate set to maximum. To account for surface turbulence and the transducer near-field, acoustic data were collected from 15 m below the surface down to 1000 m. The echosounder could not be calibrated prior, during or after survey operations, but had been calibrated on a dedicated hydroacoustic survey a few months prior with the standard sphere method (Foote et al. 1987;Demer et al. 2015). These  Ryan et al. 2015).
High backscattering values originate from targets with a density that differs strongly from the surrounding seawater that resonate when their dimensions are near the wavelength of a given frequency. Organisms with a density that is very similar to the surrounding medium produce a much weaker echo (Simmonds and MacLennan 2006). The former echoes mostly originate from organisms that bear a gas-filled structure, i.e. swimbladdered fishes, but also physonect siphonophores (Stanton et al. 1998;Korneliussen and Ona 2003;Proud et al. 2019), the latter from so-called fluid-like scatterers like cephalopods and crustaceans (Korneliussen and Ona 2003). Based on these backscattering characteristics, organisms that produce resonance at different incident frequencies can be classified.
A corresponding classification based on 18 and 38 kHz acoustic data was conducted. Acoustic backscatter was translated into volume backscattering (S v, dB re 1 m −1 ) and binned into cells of 2 min x 1 m depth for the duration of trawl sampling on each station and for 10 min x 10 m depth for a representative 72-h time-series of hydroacoustic data recorded on a transect section at 70°W illustrating several cycles of diel vertical migration. For each resulting cell, mean S v was calculated at an integration threshold of −80 dB, and synthetic variable Δ Sv was created by calculating S v 18-S v 38. This variable was used to identify functional groups of mesopelagic scatterers following a classification tree described by D'Elia et al. (2016). This classification is based on the size and scattering properties of different organism groups and allocates four categories to bins according to the following properties: small swimbladdered fishes (including small nonswimbladdered fishes and crustaceans; sm.(N)SB/Crust.) are −14 dB < Δ Sv18kHz−38kHz < −3 dB; large non-swimbladdered fishes (Lrg.NSB) are −3 dB < Δ Sv18kHz−38kHz < 0 dB; gelatinous zooplankton, cephalopods and pteropods, i.e. fluid-like scatterers (FL), are 0 dB < Δ Sv18kHz−38kHz < 3 dB, and large swimbladdered fishes (Lrg.SB) are 3 dB < Δ Sv18kHz−38kHz < 12 dB. Cells in which S v was below the threshold (−80 dB, D'Elia et al. 2016) in either the 38 kHz or the 18 kHz data were classified as dominant Dom18 and Dom38, respectively. The corresponding acoustic scattering measured in both such cells can most likely be attributed to swimbladder-bearing fishes (Love 1978).
To evaluate the contribution of different taxonomic (fish) groups to the post-processed echograms, a classification of fishes according to the presence of a gas-filled (i.e. resonant) swimbladder was conducted and the corresponding families were categorized accordingly. A corresponding categorization was achieved using (identification) literature and swimbladder catalogues (Marshall 1960;Whitehead et al. 1986;Saenger 1989). Accordingly, the families Anoplogastridae, Bregmacerotidae, Chiasmodontidae, Gonostomatidae, Howellidae, Melamphaidae, Myctophidae, Phosichthyidae, Scombrolabracidae and Sternoptychidae were classified as fishes with swimbladder, whereas Evermannelidae, Notosudidae, Paralepididae and Stomiidae were classified as fishes without swimbladder. Species for which no information could be gathered were categorized into "other fish". Ontogenetic changes in swimbladder-structure (i.e. gas-filled in juveniles and lipid-filled in adults) as occurring in some families/genera were not taken into account.

CTD and sea surface temperature satellite data
In situ hydrographic measurements were conducted with a SBE 9/11 conductivity-temperature-depth probe (CTD) (Sea-Bird Electronic, Bellevue, WA, USA). CTD casts were made at all stations (depth 300-1000 m), recording conductivity, temperature, salinity and pressure (depth).
Prominent temperature fronts along the STCZ in the study area were observed via sea surface temperature (SST) data. SST data were derived from the Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) system (https:// opendap.jpl.nasa.gov/opendap/OceanTemperature/ghrsst/ data/L4/GLOB/UKMO/OSTIA/2014/028/contents.html). Figure 1 depicts SST on March 28, 2014, as this date represents the mid-point of the sampling period.

Catch analysis and species identification
Total catches were sorted and divided into major taxonomic groups immediately after each haul and frozen at −20°C for transport and further utilization. After thawing, all fish specimens were identified to the lowest possible taxonomic level by the use of region-specific identification keys (Whitehead et al. 1986;FWNA 1989;Carpenter 2002;Richards 2005) and standard length and weight were assessed individually to the lowest 1.0 mm and rounded to the nearest 0.1 g, respectively. Regarding invertebrates, the analysis was limited to the assessment of total weight per group (mollusca, crustacea, gelatinous plankton) and haul.
Genetic analyses of subsamples of fish species were used to verify the morphological identification (for analytical methods see below). In case the results did not match with morphological species identification, morphological examination and identification were repeated and species names were revised or assigned to higher taxonomic ranks.
The influence of sea surface temperature on the abundance of the most dominant species was tested by a linear regression model (R Core Team 2019).

Genetic analysis
Tissue samples of representative individuals from each species or group were preserved in ethanol (96% abs.) for subsequent genetic analysis. DNA was extracted using Chelex100 (Walsh et al. 1991) and stored at +4°C or at −80°C for longterm storage. For DNA barcoding, the mitochondrial markers Cytochrome c oxidase I (COI) and Cytochrome b (Cytb) and the nuclear marker Myosin heavy chain 6 (myh6) were amplified by polymerase chain reactions (PCR). All PCRs were carried out using High-Fidelity PCR Master Mix with HF Buffer (Phusion, New England Biolabs), 0.25 μM of each primer (Table 2), 3 μl template (using the supernatant of the Chelex extracted samples). Finally, nuclease-free water was added until a final volume of 20 μl was reached. Amplification was accomplished with a T3 Thermocycler (Biometra/Analytik Jena, Germany) with the following cycling conditions: initial denaturation at 98°C for 30 s, followed by 32 cycles for COI and 35 cycles for Cytb, consisting of 98°C for 10 s, primer annealing for 30 s (COI at 53°C, Cytb at 51°C), extension at 72°C (COI for 20 s, Cytb for 40 s) and a final extension step at 72°C for 8 min. For myh6 amplification, a special cycling protocol was used with the following conditions: initial denaturation at 98°C for 30 s, followed by a primer annealing touchdown step with 8 cycles, consisting of 98°C for 10 s, primer annealing decreasing with 1°C from 62°C to 54°C in each cycle for 30 s, extension at 72°C for 30 s, followed by 30 cycles with 98°C for 10 s, 54°C for 30 s, 72°C for 30 s and final extension at 72°C for 7 min. As a quality check, 5 μl of each PCR product was analysed with standard agarose gel electrophoresis (1% agarose, stained with ethidium bromide). PCR products showing strong sharp and clear bands with the right amplicon size were diluted with ultra-pure water. Products showing weak bands were purified and concentrated using PCR and DNA Cleanup Kits (Monarch, New England Biolabs, T1030) prior to sequencing. Sequence raw data were checked and edited with CodonCode Aligner Software (Centerville MA, USA) by cutting off primer sites and generating consensus sequences. To verify the results of the morphological identification, all obtained sequences were compared with NCBIs Nucleotide Sequence Database by BLAST (Altschul et al. 1990).

Hydrography
Depth profiles of temperature and salinity varied among sampling stations (Fig. 2). A well-defined thermocline was present at all stations between 100 and 200 m depth. Sea surface temperatures varied between 22.5 and 24.9°C, with the lowest temperatures being measured at northern-, western-and easternmost stations (stations 284, 228 and 316) (Fig. 1). Within the upper 200 m, temperatures dropped at about 4-5°C to 18.5-20.7°C and at 500 m depth temperatures between 15.5 and 17.2°C were measured. Trawls at regular stations took place in depths with temperatures between 18 and 22°C, while deep stations reached temperatures below 7°C at greater depths. Sea surface temperature data showed a distinct temperature front between 24°N and 30°N in the western part of the study area, with temperatures decreasing northwards from above 25°C to below 22°C (Fig. 1). In the eastern part of the study area, the front was less pronounced.
Water salinity decreased vertically from around 36.5 near the surface to 35 at 1000 m depth. Within the upper 500 m, salinity values ranged from 36.1 to 37.2.

Hydroacoustics: general mesopelagic habitat structure
From hydroacoustic data, a common pattern of scattering layers was evident throughout the survey area and along the sampled transects, with persistent stationary and several migrating components of different scattering characteristics (Figs. 3, 4 and 5).
In the epipelagic zone, dense stationary echoes were visible especially in the upper 150 m of the water column (epipelagic layer, hereafter referred to as layer 1). In the mesopelagic zone, two separate deep scattering layers (DSLs) were evident: one dense layer was most prominent in the 18 kHz echodata between 450 and 600 m (Figs. 3, 4 and 5b) and between 400 and 700 m in the 38 kHz data (upper mesopelagic scattering layer, hereafter referred to as principal DSL and layer 2). A weaker scattering layer around 800-900 m depth was most prominent in the 18 kHz data (Figs. 3, 4 and 5a) (lower mesopelagic scattering layer, hereafter referred to as secondary DSL and layer 3). Based on the classification of acoustic data through Δ Sv , layer 1 contained echoes originating from a variety of organisms, namely swimbladdered as well as non-swimbladdered fishes, fluidlike scatterers and crustaceanswith varying contributions during day-and night-time. The upper part of the principal DSL (layer 2) (400-500 m) consisted mostly of (large) swimbladdered and non-swimbladdered fishes with clear contributions of fluid-like scatterers, while the deeper parts of the principal DSL (layer 2) (500-700 m) were dominated by small swimbladdered and non-swimbladdered fishes and crustaceans (Figs. 3, 4 and 5c). From the echograms recorded at either frequency, a clear and regular diel pattern was evident with different components emerging from both DSLs at different times and undertaking diel vertical migration (Fig. 3).
Beginning from ca. 16:00 h local time (ca. 2 h before sunset), an upward migrating layer emerged from the principal DSL (layer 2) and by 19:00 h merged with layer 1. A second cohort of organisms undertaking DVM emerged from the deeper secondary DSL (layer 3) around 17:00 h. This group ascended through the principal DSL (layer 2) and merged during its ascent with the previous cohort from this layer shortly before merging with layer 1. The migrating components appeared to mostly consist of swimbladdered and non-swimbladdered fishes as well as fluid-like scatterers. For the next ca. 9 h, layer 1 appeared both denser (in terms of acoustic backscatter) and of different composition (through the contribution of the migrating components). While during day-time the epipelagic layers seemed to be mostly dominated by small swimbladdered and non-swimbladdered fishes and crustaceans, distinct additional contributions of fluid-like scatterers as well as large swimbladdered and non-swimbladdered fishes were evident during night (Fig. 3c). Starting from ca. 04:00 h, several descending layers became evident leaving layer 1 and the epipelagic zone, with the components from the secondary DSL (layer 3) showing a faster descent into that depth than the components of the principal DSL (layer 2). By ca. 07:00 h, the faster descending migrating component of layer 3 merged into that layer again, while the migrating organisms from layer 2 merged into the stationary component ca. 1 h later.  conducted with uniform trawl paths and depth profiles. There, the upper net 1 (ca. 150-220 m) partly covered the dense epipelagic layer (layer 1) that appeared during night-time after migrating organisms had completed their DVM. Layer 1 was classified as most likely consisting of a mixture of components, with contributions of swimbladdered, non-swimbladdered and fluid-like scatterers. Nets 2 (ca. 220-280 m) and 3 (ca. 280-360 m) covered depths below layer 1, in which echoes were less dense and most likely originated from resonant, swimbladdered fish and, to a lesser extent, of fluid-like scatterers (Fig. 4c). In all depth layers covered by the multisampler tows, swimbladdered fish species dominated the catches by total catch weight, and fishes without a swimbladder only marginally contributed to the bulk catch across all depth layers (Fig. 6). In these nighttime catches, Myctophidae dominated especially the shallowest samples taken from layer 1 (net 1), while their number decreased in deeper samples (nets 2 and 3) ( Table 3). The second most dominant family was Melamphaidae, with highest numbers in net 2 and lowest numbers in the deepest samples. Other important families that contributed to the catch in all sampled layers, albeit in distinctly lower numbers, were Gonostomatidae (all Sigmops elongatus), Stomiidae, Scombrolabracidae (all Scombrolabrax heterolepis) and Evermannellidae (mostly Coccorella atlantica). In all but the Gonostomatidae (i.e. S. elongatus), numbers were highest in shallower samples (nets 1 and 2). The latter occurred in higher numbers in the deeper tows (nets 2 and 3). Sternoptychidae were only sampled in the deeper layers (nets 2 and 3) with highest numbers in the deepest tows. Other fish families were caught in distinctly lower numbers.

Depth distribution of mesopelagic fish and taxonomic composition of scattering layers
No clear trend was evident in the fractions of other organisms across hauls. The second most important organism group were molluscs (mostly cephalopods). In general, molluscs contributed stronger to the total catch weight of the shallowest net than of the deeper nets and they are assumed to contribute to the fluid-like backscatter measured in these depths. Crustaceans were also present in all hauls with similar fractions across all sampled depth layers (Fig. 6). Gelatinous zooplankton only marginally contributed to the catches. . Cells in which S v was below the threshold in either the 38 kHz or the 18 kHz data were classified as dominant Dom18 and Dom38, respectively. The corresponding acoustic scattering measured in both such cells can most likely be attributed to swimbladder-bearing fishes (Love 1978). White, dashed vertical lines indicate sunset and sunrise, respectively

Deep stations
On the deep day-time station (233), located at the western part of the study area, an increase in total catch (abundance and biomass) was evident with increasing depth and a dominance of different organism groups in terms of catch weight became evident (Table 4). Samples collected at that station were taken in depths covering the layers 1 (net 1, 108-362 m) and 2 (net 2, 362-450 m and net 3, 450-698 m). Net 1 was dominated by molluscs (ca. 60% of catch weight), followed by crustaceans and swimbladdered fishes (both ca. 20% of catch weight) (Fig. 6). According to the scattering properties described above, this depth contained swimbladdered fishes as well as fluid-like scatterers and contributions of small swimbladdered and non-swimbladdered fishes and crustanceans (Fig. 5c). Net 2, i.e. samples from the upper part of layer 2, was dominated by gelatinous zooplankton (ca. 50% of catch weight), followed by swimbladdered fishes (ca. 30% of catch weight) and crustaceans (ca. 20% of catch weight). Corresponding echo signals were classified as consisting mainly of swimbladdered fishes and fluid-like scatterers (Fig. 5c). Net 3, the main part of layer 2, showed a similar catch composition, with gelatinous zooplankton contributing ca. 75% to the catch weight, swimbladdered fishes accounting for ca. 20% and crustaceans for ca. 5%. The depth range covered by net 3 was dominated by the acoustic classes of small swimbladdered and nonswimbladdered fishes and crustaceans (Fig. 5c).
Myctophidae, followed by Gonostomatidae, was the most important fish family caught at that station. The bulk of specimens was caught in the deepest net 3 sampling layer 2 (ca. 450-700 m), with distinctly lower numbers in net 2 and lowest numbers in net 1. Gonostomatidae were absent from the shallower nets 1 and 2.
The deep night-time station (316) was located approximately 1,790 km further east from the deep day-time station (233) and covered the water column down to 965 m. At this station, layer 1 (net 1, 32-152 m) and the upper part of layer 2 (net 2, 482-606 m) as well as layer 3 (net 3, 775-965 m) were sampled (Fig. 5c) and a decrease in the number of individuals became evident with increasing depth. Catches in terms of weight were dominated by swimbladdered fishes in the shallowest net (net 1, ca. 85%) with small contributions of gelatinous zooplankton and molluscs (ca. 6% each) and . Cells in which S v was below the threshold in either the 38 kHz or the 18 kHz data were classified as dominant Dom18 and Dom38, respectively. The corresponding acoustic scattering measured in both such cells can most likely be attributed to swimbladder-bearing fishes (Love 1978). For each station, trawl paths, i.e. regions and layers covered by the three nets of the multisampler, are indicated by red rectangles in panel c negligible fractions of crustaceans. A mixture of the corresponding categories had also been allocated to the hydroacoustic data collected from this layer. Fishes collected with that net were almost exclusively Myctophidae, with Phosichthyidae and Melamphaidae only marginally contributing to the catch (Table 4). Catches from net 2, sampling layer 2, were also dominated by fishes, with fishes without a swimmbladder constituting the bulk (50%) and fishes with a swimbladder and uncategorized fishes contributing smaller fractions (17 and 11%, respectively) (Fig. 6). Crustaceans contributed with 15%, while gelatinous zooplankton and molluscs each contributed with less than 5% to the total catch weight. Among fishes, Melamphaidae, Myctophidae and Gonostomatidae were dominant, albeit in far lower numbers than in the upper 150 m (Table 4). Other families, such as Evermannellidae, Stomiidae, Chiasmodontidae and Sternoptychidae, also contributed to the catch with several individuals of each family. In the deepest net (net 3), gelatinous zooplankton was the dominant fraction (40%), followed by crustaceans (20%) and fish (swimbladdered, no swimbladder and others with 19, 14 and 4%, respectively). Fish families sampled at that depth were mostly Stomiidae, Myctophidae, Gonostomatidae and Sternoptychidae (Table 4). Species that were caught exclusively at depth below 360 m were Bonapartia pedaliota (N=67), Cyclothone spec. (N=9), Serrivomer beanii (N=4) and Poromitra crassiceps (N=1).

Taxonomic and quantitative fish catch composition
The total catch consisted of 5022 fish specimens from 27 families, 62 genera and 70 species of which 4050 individuals (80.7%) were identified to species level, 813 (16.2%) to genus level, and 137 (2.7%) to family level. Twenty-two specimens (0.4%) remained entirely unidentified due to mechanical damages in the net during trawling. An overview of the total catch at stations and depths is presented in Table 5, detailed information about catch composition is given in Table 6, and length and weight data as well as results from genetic analyses are presented in Table 7.
In terms of abundance (N), relative abundance (A) and biomass (M), the total fish catch was dominated by Myctophidae (N=2970, A=59.1%, M=47.4% of total catch), (including small non-swimbladdered fishes and crustaceans). Cells in which S v was below the threshold in either the 38 kHz or the 18 kHz data were classified as dominant Dom18 and Dom38, respectively. The corresponding acoustic scattering measured in both such cells can most likely be attributed to swimbladder-bearing fishes (Love 1978  The most species-rich families were Myctophidae (14 genera, 31 species) and Stomiidae (12 genera, 12 species) followed by Sternoptychidae and Phosichthyidae (3 genera, 4 species), and Melamphaidae, Evermannellidae, Gonostomatidae and Paralepididae (3 genera, 3 species). Genera with more than two species were Diaphus (10 species) and Hygophum (4) (both myctophids).
Genetic analysis of samples from 77 different species and groups of higher taxa were performed. In 58 cases, the morphologic identification was confirmed by genetic results. In 11 of these cases, specimens were assigned to higher taxonomic ranks, because morphologic and genetic identification did not match at species level (i.e. results did not match with sufficient    Chauliodus spec. 2 2 quality to any species in identification keys and NCBI Nucleotide Sequence Database). In 19 cases, a genetic confirmation was not possible because the identified taxa were not included in the Database, and in 40 cases, genetic analyses were not performed. The relative abundance of Myctophidae was high (42.9-67.7% of total catch) in catches from all regular stations (stations 228, 240, 259, 284, 300), whereas Melamphaidae showed lower abundances at western stations (stations 228 and 240; 17.4-19.1%) compared to the central and eastern sampling area (stations 259, 284, 300; 31.0-36.7%) (Fig.  7). The dominance of single species in catches was more pronounced in the eastern part of the survey area compared to the west. At the easternmost station (300), the two most abundant species accounted for 51.8% of total catch, while at intermediate stations three species and at the westernmost station (228)  Among regular stations, species richness was highest in catches at the northernmost station (station 284), where 50 different species have been caught. Evermannella indica, Macroparalepis brevis, Melanostomias tentaculatus, Nemichthys scolopaceus and Odontostomops normalops were exclusively caught at that station. Catches at regular stations further south contained 44 to 46 different species.

Discussion
The catch composition of mesopelagic fish species during this survey is similar to catches reported from previous investigations in the Sargasso Sea and other North-Atlantic areas (Backus et al. 1969;Ross et al. 2010;Olivar et al. 2017). Stratified night-time catches in layers 1 (epipelagic layer) and 2 (principal DSL) at regular stations generally reflected the taxonomic and quantitative composition that was encountered in previous surveys in the eastern part of the Atlantic in the same ecoregion "Central North Atlantic" (Ariza et al. 2016;     Olivar et al. 2017;Sutton et al. 2017). It is known that organisms from epi-and upper mesopelagic layers exhibit strong DVM and generally show the strongest differences in day-and night-time distribution with aggregations in the epipelagic zone during night (Roe and Badcock 1984;Olivar et al. 2017  and 800 m during the day and a more widespread distribution during night-time (Barlow and Sutton 2008;Sutton et al. 2008;Sutton 2013;Olivar et al. 2017). Some gonostomatid species are not known to display extensive DVM and usually are found in layers below the epipelagial both during day-and night-time (McClain et al. 2001;Olivar et al. 2017). While in this study similar observations were also made from the deep hauls, where Gonostomatidae were sampled in the corresponding layers, the species sampled in epipelagic nighttime catches are characteristic of the upper mesopelagic zone and are known to undertake DVM (Sutton 2013).
The here presented catches of stomiids reflected their known distribution and migration pattern as previously reported by Kenaley (2008) and Olivar et al. (2017), with night catches of individual specimens in layer 1, while during day-time this family was only sampled below 400 m depth. The comparatively large number of species of the subfamily Melanostomiinae caught in this study during night-time in layers 1 and 2 is a further indication of vertical migration behaviour in these species. This corroborates accumulating evidence from oligotrophic, but also productive regions, that these species, despite remaining likely invisible in hydroacoustic recordings, contribute considerably to vertical energy fluxes (Cook et al. 2013;Eduardo et al. 2020a;Czudaj et al. submitted). Sternoptychidae in general are considered limited or only partial migrants (Kinzer and Schulz 1988) and usually occur preferentially in deeper layers, although some shallower occurrences can be observed during night hauls, especially in oligotrophic regions as shown for the south western Atlantic (Olivar et al. 2017;Eduardo et al. 2020b). This is well in line with the observations in this study, as night-time catches of sternoptychids, mainly dominated by Argyopelecus aculeatus, were present in the deeper nets of the epipelagic zone (below layer 1) and also occurred in layers 2 and 3 both during day and night ( Table 6). The general observation of characteristic scattering layers in different depths in this study that show rather unique backscattering characteristics and regular DVM of organisms originating from different depth zones, was consistent with observations from adjacent areas and also from across the pelagic zones of the global ocean (Sutton 2013;Peña et al. 2014;Ariza et al. 2016;D'Elia et al. 2016;Klevjer et al. 2016). The observed presence of an epipelagic layer (layer1), a strong upper mesopelagic deep scattering layer (layer 2, principal DSL) as well as a weaker lower deep scattering layer (layer 3, secondary DSL) was in line with general descriptions of the ocean's deep mesopelagic zone (Proud et al. 2017). The detected layers with different scattering properties highly resembled layers characterized at the Canary Islands in the eastern Atlantic and the Northeastern Atlantic in general: an epipelagic layer with strong diel differences in scattering characteristics; a stationary upper layer of the upper principal DSL (upper layer 2) dominated by 18 kHz (450-550 m); a stationary lower layer of the principal DSL (layer 2) with a dominance of 38 kHz (ca. 600-700 m) and a mixed layer between; weak or no signals between 700 and 800 m; a permanent weak layer dominant at 18 kHz (secondary DSL, layer 3) between ca. 800 and 1000 m (Ariza et al. 2016;Klevjer et al. 2016). Additionally, similar patterns of a different but recurrent temporal onset of DVM of groups with differing backscattering characteristics at nightfall and the later descent into deeper layers with sunrise, as shown in this study, have been observed in nearby regions to the southwest (D'Elia et al. 2016) and east (Ariza et al. 2016) of the survey area. However, in this study, the sampling regime did not allow targeting and resolving the different upward and downward migrating fractions from different layers or the composition of the different layers itself. Nonetheless, based on the backscattering characteristics, and in comparison with previous studies, it seems likely that the upper part of the principal DSL (layer 2) mostly consisted of small, swimbladdered fishes resonant at 18 kHz like myctophids that also contributed the bulk of migrating fish emerging from that layer at nightfall (Peña et al. 2014(Peña et al. , 2020Ariza et al. 2016). While the lower part of the principal DSL (layer 2) remained rather stationary, migrant fauna also contributed to this layer, albeit to a lower extent. Although it was only caught in low numbers in this study, the gonostomatid Cyclothone spp. seem to be the dominant fish in the principal DSL (layer 2) (compare Ariza et al. (2016)). In contrast to Ariza et al. (2016) and Sutton (2013), who described the lower mesopelagic zone/secondary DSL (800-1000 m) as permanent and stationary, but identified a weak migrant signal between 700 and 800 m that disappeared at night, in this study a migrant pathway clearly became visible emerging from the secondary DSL (layer 3) and migrating upward through the principal DSL (layer 2), with no clear changes in the signal of the stationary fraction of that zone. Based on the dominance of the 18 kHz echoes, it can be assumed that the migrating and the stationary signal also originate from swimbladdered fishes and gas-bearing organisms. From that layer not only a peak in non-migrant fishes like Cyclothone spp. was reported earlier, but also of migrating fishes like the myctophid Notoscopelus spp. (Badcock and Merrett 1976;Roe and Badcock 1984;Sutton 2013;Ariza et al. 2016). The latter haveamong otherscontributed to epipelagic night-time catches made during this study. It has to be mentioned that a certain amount of backscattering in the lower frequency range also is highly likely to originate from physonect siphonophores carrying a gas-filled pneumatophore for buoyancy (Kloser et al. 2016;Proud et al. 2019). While comparatively little is known on actual vertical distribution or migration of these organisms (Pugh 1975;Pugh 1984;Lüskow et al. 2019), siphonophores are known to inhabit a broad depth range in the epi-as well as the upper mesopelagial, and some species also undertake DVM. Despite gelatinous zooplankton not being identified to lower taxonomic levels due to the fragility of the organisms and their condition in the codend of the multisampler nets, the contribution of gelatinous zooplankton to the total catch weight especially in the deep tows (net 3 of stations 233 (450-698 m) and 316 (774-965 m)) indicates that a certain fraction of the scattering layers with a dominance at 18 and 38 kHz may be assigned to a contribution from siphonophores. The presence and condition of swimbladders changes among and within taxonomic groups of mesopelagic fishand even may vary within genera and species based on developmental stage and length (Marshall 1960). It has been shown that e.g. biomass estimates of mesopelagic fish based on acoustic data collected at 38 kHz can be complicated because of the small physical size of mesopelagic fauna, ontogenetic changes in swimbladder morphology, inflation and regression (Davison et al. 2015). Additionally, the backscatter is depthsensitive and non-linear with respect to size; at the same time, the size structure of mesopelagic fish is skewed with abundance driven by the smallest and biomass driven by the largest fishes. Echograms accordingly rather reflect the distribution of the strongest scattererse.g. (small) fishes with swimbladdersthan the actual distribution of biomass (Davison et al. 2015). While some of the individual species identified in the multinet hauls may have different swimbladder characteristics than other genera/species from the same family, it is assumed that the main drivers of backscatter characteristics of scattering layers are still correctly identified through the classification used in this study.
A classification of hydroacoustic data into categories based on their scattering properties (e.g. D´Elia et al. 2016) and in combination with trawl net hauls is feasible, but a comprehensive interpretation of the aggregated results is challenging and corresponding sampling biases have been reported ). This also affects contributions of gelatinous zooplankton and other invertebrates that may significantly contribute to backscatter through resonance at low frequencies but will be virtually absent from trawl haul catches due to the small and/or fragile nature of these organisms. It can safely be assumed that the corresponding layers identified using the method described by D' Elia et al. (2016) are not exclusively inhabited by the dominant taxonomic groups triggering the classification, but also by a magnitude of other species whose acoustic signal is masked by the dominant scatterers as well as the range limitation inherent in hydroacoustic data from higher frequencies.
An indication for such a masking and/or a "missed classification" is evident from the catch composition of some tows that were conducted within layers with typical characteristics and dominance of gas-bearing organisms (i.e. resonant at 18 and 38 kHz), but also contained fluid-like zooplankton organisms like crustaceans and molluscs.
As the scientific echosounder could not be calibrated directly prior to, during or after the survey, but had been calibrated with good results on a previous survey and was again calibrated a few weeks after the survey, the utilization of transducer parameters from a preceding calibration updated with ambient physical measurements from the current survey is considered sufficiently precise for the classification attempted here. The classification approach followed in this study is based on Sv intervals that are used to differentiate biological groups. D' Elia et al. (2016) derived the corresponding parameters from length measurements of representative organisms from concurrent net samples and on theoretical models relating these lengths to target strengths and subsequently to the Δ Sv intervals. Since this study was conducted in a different, albeit adjacent, ecoregion with potentially different species and length compositions (e.g. Sutton et al. 2017), a certain degree of misclassification cannot be ruled out. Cells that did not match Δ Sv intervals specified in D' Elia et al. (2016) could nevertheless be allocated to e.g. swimbladder-bearing fishes based on their scattering properties (Love 1978). Accordingly, the overall classification is considered robust enough for the analyses conducted.
According to Backus et al. (1970), who defined oceanic areas in the Atlantic on the basis of characteristic water masses, and supported by the hydrographic data obtained during this survey, the here investigated area is part of the Southern Sargasso Sea. This zone is characterized by distinct temperature fronts (Fig. 1), which are present from fall to spring (Halliwell et al. 1991). In this study, species richness was slightly higher at the northernmost station (lowest SST) compared to warmer stations further south with five species exclusively caught at the northernmost station. However, since only single individuals of these species were caught, the observed higher species richness north of the front might also be an artefact of low sampling effort. Nonetheless, similar results were reported by Backus et al. (1969), who caught a number of species north of the temperature front that were absent further south. In this study, northern and southern fish communities could not be compared in detail, since the sampling grid was not designed accordingly, but the total survey area comprised the Central Atlantic Mesopelagic Ecoregion that in general has a distinct faunal composition, albeit with spatial differences that are based primarily on abundance and rank order rather than presence or absence of species (Sutton et al. 2017).
Variations were observed between eastern and western stations with regard to the predominant species in catches (i.e. S. mizolepis and N. cuprarium vs. B. photothorax and C. warmingii). These differences were apparently not driven by water temperature, since no significant correlation was evident between SST and the abundance of the four most abundant fish species. However, in case of C. warmingii (p=0.06, R 2 =0.74), the absence of a significant correlation might be the result of the limited number of sampling stations. The relevance of temperature as a major driver for fish distribution might be lower in species that experience strong temperature gradients on a daily basis caused by DVM. Other environmental conditions like primary production, oxygen concentration and light attenuation where also shown to affect the distribution of mesopelagic fish species (Irigoien et al. 2014;Klevjer et al. 2016;Aksnes et al. 2017), but were not analysed in this study. Hence, the SST data collected in this study do not provide a full picture of environmental influences on the mesopelagic fish community, but contain information about the potential effect of horizontal temperature structures like temperature fronts on the distribution of mesopelagic fish. Nonetheless, the absence of a temperature effect on changes in the distribution of the dominant species suggests that the here observed changes in catch composition reflect differing regional influences from adjacent water masses.
Due to methodological constraints in sampling (large mesh size, relatively small number of stations, sampling not only in peak scattering layers), the data obtained and presented here do not provide a fully comprehensive and representative picture of the mesopelagic species composition of the investigated area. Stations, catch depths and trawl paths of this study were defined by the needs of the original purpose of the survey and not for the investigation of the mesopelagic fish community, which was only an additional benefit of the survey. As a consequence, no day-time samples were collected from the epipelagic scattering layer (layer 1) and the deep scattering layers (layers 2 and 3) could only be sampled at two stations at the far ends of the survey area. This study design hampers the discernibility of permanent residents and migrators in the different layers and complicates a more detailed interpretation of catches with regard to DVM. Another limitation of this study is the relatively large mesh size of the mesopelagic trawl. This might have led to the under-representation of certain small and thin species in catches, despite a general capture efficiency of the gear also for small-sized species < 30 mm (Fock and Czudaj 2018). For example, bristlemouths of the genus Cyclothone are usually found in high abundance in the investigated area and have been reported as resonant scatterers at depth (Peña et al. 2014). In this study, Cyclothone represents only a minor fraction of the trawl samples, which most likely is due to the large mesh sizes in the codend of the trawl net employed (1.5 cm), as it already was the case in the eastern tropical North Atlantic (Fock and Czudaj 2018;Czudaj et al. submitted). Large mesh sizes may also be the reason why Ranzania laevis, a small species of the sunfish family Molidae, was not present in any of the hauls, even though spawning activity and increased larval and post-larval abundance during the March and April in the area were recently described by Hellenbrecht et al. (2019). It is noteworthy though that nine post-larval (Molacanthus stage) specimens of Masturus lanceolatus, another fairly unexplored species of Molidae, were caught in this study.
Despite aforementioned limitations, the current study provides valuable insights into the distribution and vertical migration behaviour of mesopelagic fishes in the Sargasso Sea and adds to our understanding of the mesopelagic community in this large oceanic area. To assess how and to what extent the fish community in this ecoregion is affected by the influence of adjacent areas and by (changing) hydrographic conditions requires additional effort and extensive further investigations.
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/.