Characterization of sound scattering layers in the Bay of Biscay using broadband acoustics, nets and video

Sound scattering layers (SSLs) are observed over a broad range of spatio-temporal scales and geographical areas. SSLs represent a large biomass, likely involved in the biological carbon pump and the structure of marine trophic webs. Yet, the taxonomic composition remains largely unknown for many SSLs. To investigate the challenges of SSL sampling, we performed a survey in a small study area in the Northern Bay of Biscay (France) by combining broadband and narrowband acoustics, net sampling, imagery and video recordings. In order to identify organisms contributing to the observed SSLs, we compared measured frequency spectra to forward predicted spectra derived from biological data. Furthermore, to assess the confidence in SSL characterization, we evaluated uncertainties in modeling, acoustical and biological samplings. Here, we demonstrate for the first time that SSL backscattering intensity in the Bay of Biscay can be dominated in springtime by resonant gas bearing organisms below 100 kHz, namely siphonophores and juvenile fishes and by pteropods at higher frequencies. Thus, we demonstrate the importance of broadband acoustics combined to nets, imagery and video to characterize resonant backscatterers and mixed mesozooplankton assemblages.


Introduction
Sound Scattering Layers (SSLs) are routinely observed with active acoustic devices in a great variety of ecosystems and over wide depth ranges in the global ocean [1][2][3][4].
Deep Scattering Layers [5] inhabiting the mesopelagic zone worldwide, are e.g. known to perform daily the largest migrations on earth [6] and their fish component might dominate the world total fishes biomass [4,7]. SSLs might hence play an important role in the biological carbon pump [8] and in the structure of marine trophic webs [9,10].
SSLs are generally composed of mixed species assemblages of great diversity in taxa, body size and acoustic backscattering. Mesozooplankton (0.2-20 mm, [11]), macrozooplankton (> 20 mm, mainly shrimp and gelatinous zooplankton, [11]) and micronekton   mainly small mesopelagic fish and other fish juveniles, [12]) have hence been reported to significantly contribute to SSLs backscatter. If SSLs are easy to observe with active acoustic devices, their taxonomic and size diversity makes them difficult to sample and identify [1]. Uncertainty of the species composition in the SSL's generally hinders the comprehension of their ecological functions. To overcome SSLs composition identification issues, we must to understand: i) the limits of each sampling method, and ii) how these methods can complement one another. Net sampling is a traditional way to identify organisms composing SSLs, but it is discrete in time and space and highly dependent on target catchability, as some organisms actively avoiding nets [13]. Furthermore, nets generally destroy fragile organisms such as jellyfish. Optical sampling is a good way to evaluate the presence of fragile organisms, albeit with low range and small sampling volumes. Still, like nets, optical sampling is affected by avoidance and behavior of organisms [14]. Acoustic sampling is continuous over various spatial scales with good resolution. It is however an indirect sampling method with low taxonomic resolution [15] and sensitivity to shadowing issues [16]. Additionally, the backscattering intensity of different types of scatterers is extremely variable [16] and the strongest backscatterers are not necessarily the most abundant. Until recently, most SSL acoustic studies were conducted using narrowband, multi-frequency acoustics [1,14,17]. However, the emergence of broadband acoustics [18] as standard monitoring tools has the potential to improve species identification based on their broadband frequency spectra [19][20][21].
Ubiquitous dense SSLs have been routinely observed during the PELGAS (PELagic GAScogne) multidisciplinary integrated survey which has taken place in springtime since 2000 [22]. Still to date, the composition of these layers remained largely unknown. Interestingly though, the presence of small gas-bearing targets in the Bay of Biscay water column was demonstrated by the detection of distinct acoustic resonance peaks [23,24].
Resonance occurs when some specific types of targets (as gas bubbles) increasingly vibrate upon insonification at a given frequency [25,26]. This natural frequency of the target depends on the target size and its physical properties as well as the one of the surrounding medium. This resonance can induce a dramatic increase in the backscattering intensity of small gasbearing organisms (e.g. physonect siphonophores or small swimbladdered fish), relatively to their size [4,14]. Previous studies suggested that Biscay dense SSLs were mainly composed of gaseous resonant target, as fish larvae, gelatinous gas-bearing siphonophores or phytoplankton [23,27]. But no biological sampling had supported these conclusions.
In this study, we investigated the composition of SSLs observed on the continental shelf of the northwestern Bay of Biscay. We combined broadband acoustic data with alternative sampling methods, including optical and net recordings. Our goals were to (i) evaluate the input of broadband acoustics in SSL studies, (ii) improve the identification of organisms dominating the SSL backscattering by using a forward approach [28], and (iii) investigate the various sources of uncertainty associated to each SSL characterization method.

Sampling
The PELGAS survey take place every year with the French Navy (Marine Nationale) permission. The 2016 permit number is: IF 14 CECLANT/OPS/OPSCOT. Sampling took place on the R/V Thalassa at the end of the PELGAS2016 survey [29] in the Northern Bay of Biscay, between the 27 th and the 28 th of May. SSLs were sampled over a diel cycle in one area off south Brittany (Fig 1). Here, we present the daytime biological sampling and the corresponding acoustic data that were collected at speeds between 2 and 3 knots. [31]), following the 2016 USA-Norway EK80 Workshop recommendations [18]. These code routine are implemented in MOVIES3D freeware [32]. An important step in broadband data post-processing is the choice of the Fast Fourier Transform (FFT) window width (i.e. the amount of data points used to compute frequency spectra) [21]. Reducing the FFT window length increases the spatial resolution, but decreases the frequency resolution. We used a frequency resolution of 2 kHz for the spectrum analysis, corresponding to an effective temporal resolution of 0.34 m. A 0.6-Tukey shading window was applied to the data before applying the FFT to reduce the leakage of spectral analysis. Individual S v spectra were obtained by an echo-integration of 1.5 m vertically and 5 pings horizontally, which covered a horizontal distance of approximately 15-23 m. This integration volume contained approximately 20 points, and was a trade-off between a stable frequency response and spatial resolution (ideally around 100 points, [33]).
Biological sampling. Biological sampling was performed with a Methot Isaac Kidd (MIK) net for specimens over 2 mm and with a Multinet (Hydro-Bios, Germany) for individuals under 2 mm. Net characteristics are presented in Table 2, using [34] tow type nomenclature. Net tow depths and durations are presented in S1 Table in the Supporting Informations. All nets were equipped with a depth sensor. The multinet was equipped with 5 nets and an opening/closing system triggered by a pressure sensor. The MIK net was equipped with a GOPRO Hero 4 video camera encased in a pressure container fitted with two 2000 lumens white LED spotlights. The video camera system was used to detect the presence of fragile centimetric organisms that would have been destroyed by the net. The camera was oriented towards the Additionally, we performed Conductivity-Temperature-Density (CTD, Seabird, USA) vertical profiles. The CTD probe was equipped with a fluorometer (Seabird, USA), a turbidimeter (Campbell Scientist, USA) and an oxygen sensor (Seabird, USA). The CTD profiles are presented in S3  We sampled two different depth layers, corresponding to the depth ranges of the densest and most structured SSLs observed in real time on pulse compressed echograms (Fig 2): i) a surface layer (from 10 to 24m depth, located at the thermocline/halocline level) and ii) a deep layer (from 70 to 96m depth). Each layer was sampled by two MIK tows and two multinet net tows.
Micronekton and gelatinous zooplankton caught by the MIK nets were identified and measured immediately after collection, or later in the laboratory after formalin fixation (4%). Multinet mesozooplankton samples were also imaged immediately after collection using the ZooCAM in-flow imaging system [35]. Mesozooplankton between 1 and 3 mm from the MIK samples were imaged using a ZooScan [36] after formalin fixation (4%) in the laboratory. The Ecotaxa web application [37] was used to identify, count and archive images originating from both instruments. The video footage recorded during MIK tows were visually analysed to detect the presence of large gelatinous organisms, such as jellyfish or siphonophores.
Density and size of fish and euphausiids larger than 2 mm were estimated from MIK samples. Density and size of mesozooplankton (copepods, pteropods) and euphausiids smaller than 2 mm were estimated from Multinet samples.
Large physonect siphonophore bodies were visually counted on GOPRO videos. As they are always composed of colonial bodies with one single gaseous inclusion (the pneumatophore), densities of pneumatophores were estimated by dividing organism counts by the net sampling volume. We assumed that pneumatophores larger than 0.3 mm found in the catches came from large siphonophore bodies observed on the videos recorded during the same tow. Siphonophores bearing pneumatophores smaller than 0.3 mm were too small to be observed on videos. Their densities were estimated based on the number of pneumatophores found in the Multinet samples which we processed using imagery. A small number of other weak acoustic scatterers (small jellyfish, comb jelly and swimming crabs) were also caught but were not included in the forward approach.

Forward approach
Using the forward approach [28,38], a predicted Volume backscattering Strength (S v , dB re 1m −1 , [31]) at a given frequency (f) is estimated with individual Target Strength (TS, dB, [31]) models, parameterized by the biological samples (size and shape). Predicted S v (f ) were then compared to our in situ S v (f ) averaged in the area where biological samplings were performed.
The predicted S v (f ) (S v predicted (f )) were calculated for each depth layer as the sum in natural scale of the frequency spectrum S vj (f ) of each sampled taxa j (Eq 1): S vj (f ) was calculated as: Where TS i represents the modeled Target Strength of length class i in dB re 1m 2 and D i the sampled density of the organisms (in individual by m 3 ) belonging to the i th length class (Eq 2). Models and parameters used for each organism are presented in the next section. The measured S v (f ) (S v measured (f )) of each layer were calculated as the median of the S v (f ) in echo-integration cells along the sampling net trajectory. These cells were defined as 2 echointegration cells located below the top of the net. Uncertainty in the predicted and measured S v (f ) curves are considered below.

Scattering models
Fluid-like organisms. A scattering model based on the distorted-wave Born approximation (DWBA) [39] was used to model the backscattering of Fluid-like (FL) organisms [16], for all frequencies and averaged over a normal distribution of orientation angles. A prolate spheroid model was used for copepods and an uniformly-bent and tapered cylinder model for shrimp-like organisms and fish lacking swimbladders. Model parameters are summarized in Table 3. Euphausiids model parameters come from [2], copepods model parameters from [40] and fish without swimbladders model parameters from [41].
In the case of net samples processed with imagery, the length of FL organisms was estimated as the major axis of their silhouette area observed on imagery pictures. The radius a of an organism's cross section was estimated as the semi-minor axis of the planar projection of the model shape (rectangle for a cylinder, ellipse for a spheroid) whose area was equal to the silhouette area. The mean ratio L/a of body Length (L) over cross section radius (a) was calculated on a subset (n > 30 by scatterer types) of imagery pictures without flexion or cut body part. In the case of organisms processed manually, the mean ratio L/a was estimated based on a subset of a dozen individuals.
Elastic shell organisms. Thecosome gastropods (common name pteropods) produce relatively high backscattering in geometric scattering compared to FL organisms of similar size [16]. Pteropods are caracterized by a hard aragonite shells, with a large discontinuity: the opercular opening [2]. Based on the review of Elastic Shell (ES) scattering models by [2], we selected a high-pass dense-fluid sphere model with an empirically derived reflection coefficient (R = 0.5) to predict the backscattering intensity of ES organisms [16].
Pteropods sampled in the nets were measured by plankton imagery. Their cross section radius a was estimated as the radius of a disc of the same area as the organism's silhouette image.
Gas-Bearing organisms. We used a gaseous prolate spheroid model [42], supplemented with [26]'s damping equation (see details of the model in S2 Appendix) to predict the backscattering of Gas Bearing (GB) organisms [16]. This model is valid for Rayleigh scattering, i.e. for ka < kb < 0.1; where a and b are respectively the semi-minor and semi-major axes of the spheroid in meter and k is the wave number (m −1 ).
Swimbladdered fish viscosity and surface tension parameters come from [43] and [44]. Swimbladder semi-minor axis a and length-to-width ratios (L/2a) of transparent fish (juveniles or larvae of Crystallogobius linearis and Carapus acus) were visually measured in a sample of thirty specimens, all of which were in good condition. A mean ratio relating fish body length to swimbladder length was derived based on these measurements. This ratio was used to predict the swimbladder cross-section of damaged individuals.
When direct measurements were impossible, we approximated the fish body volume as a prolate spheroid, with the fish length as the major axis and 1/6 of the fish length as the minor axis (based on the mean of all visual measurements). Following the estimations and measurements of [45], [46] and [47], swimbladder volume was then estimated as 2.5% of the body volume, using an empirical swimbladder length-to-width ratio of 1.5. These values were in line with the swimbladder to body volume estimations made on transparent fishes.
Pneumatophores of physonect siphonophores displayed a pneumatophore mean length-towidth ratio (L/2a) of 2.35, regardless of their size. As no specific value of viscosity and surface  [44]) and water (0.0013 Pa.s −1 and 0.075 N.m −1 ) and to provide the best fit between the measured and predicted S v (f ). Parameters used in GB organisms scattering models are presented in Table 4.

Uncertainty analysis
Uncertainty analysis of scattering models. Model parameters used in this study to estimate organisms backscatter were not specific to the Bay of Biscay. We performed an uncertainty analysis for each scattering model to assess the reliability of model predictions, when accounting for uncertainty in model parameters derived from the literature. Reference values for parameters other than length (that was measured on all organisms) are presented in Tables  3 and 4. The standard deviations of the parameter distributions were based on the range of the literature. Alternative values for the parameters were drawn from statistical distributions centered on reference values. The mean, standard deviation and quantiles of these model parameter distributions are presented in Table 5. The uncertainty on the FL organism's orientation was represented by using a mean angle distribution in the DWBA model (Table 3) and was therefore not included in the uncertainty analysis.
The uncertainty analysis consisted of 1000 model simulations run in each 0.1 μm scatterer size class. The parameter values were drawn from the distributions presented in Table 5, using Latin Hyper Cube sampling [48]. This near-random sampling method ensured that each parameter was sampled according to its distribution, and that the parameter space was homogeneously explored. Confidence intervals around scattering model predictions were defined as the 5% and 95% quantiles of the simulation results.
Variability of in situ frequency responses. The SSLs spatial and frequential heterogeneity along net trajectories was estimated to assess the representativity of biological sampling. Clusters of echo-integration cells were defined on echograms based on S v (f ) and net tracks were superimposed on the segmented echograms. To assess frequency spectra heterogeneity along each net track, we compared median frequency spectra computed along net tracks to the median S v (f ) of each cluster. Echogram segmentation was performed using Expectation Maximization (EM) clustering [49], following [50] recommendations. Data used in EM clustering were standardized by the mean of the S v (f ) spectrum weighted by each transducer bandwidth, as follows:  [43] and [44] b No specific value for viscosity/surface tension of the pneumatophores, they were fixed arbitrarily https://doi.org/10.1371/journal.pone.0223618.t004 Characterization of sound scattering layers in the Bay of Biscay using broadband acoustics, nets and video with s vn the s v mean value (in natural scale) over the whole bandwidth of the transducer n: with Sv b the backscattering at the b frequency of the transducer n bandwidth, f bandwidth the number of frequency bins within the bandwidth, fmin and fmax the bandwidth boundaries of transducer n, presented in Table 1.
We hypothesized that layers were stratified horizontally (coherent with the echograms on Fig 2), thus each "k" cluster was initialized by the mean (in natural scale) S v (f ) of "k" equal horizontal bands stratifying the layer. The number of clusters used in the final segmentation was set to maximize the diversity of S v (f ) spectra shapes while minimizing the number of clusters.

Biological sampling
The density, mean length and Standard Deviation (SD) of sampled organisms in each zone and layer are presented in Table 6 Table 4 https://doi.org/10.1371/journal.pone.0223618.t005 Characterization of sound scattering layers in the Bay of Biscay using broadband acoustics, nets and video We found Pteropods (Limacina sp.) and Copepods in high density near the surface (5300 ind.m −3 and 1403 ind.m −3 respectively). High densities of small euphasiids (4.26 mm of mean length, 14 ind.m −3 ) were also found in the surface layer, while large euphausiids were sampled in the deep layer (14 mm of mean length, 1.5 ind.m −3 ).
We sampled two sizes of physonect siphonophore pneumatophores: small siphonophores (mean pneumatophore length 0.3 mm) present in the surface layer (8.3 ind.m −3 ) and larger siphonophores (mean pneumatophore length 2 mm) caught at a lower density in the deep layer (0.5 ind.m −3 ). The body length of large siphonophore was estimated at about 1 meter, based on video recordings, using the known MIK net opening as scale. The body size of small siphonophores was unknown as they were not observed by the video camera.
We only caught one fish in the surface layer, a Ciliata mustela larvae. In the deep layer we mostly caught non-swimbladdered juveniles of Argentina sphyranea [51] which were more abundant in the catch than swimbladdered fish species (0.012 ind.m −3 and 0.0065 ind.m −3 respectively). The swimbladdered fish were mostly composed of Crystallogobius linearis juveniles (50% in number) and Carapus acus larvae (20%). The remaining 30% was composed of juvenile gobiidae, blenniidae, triglidae and gadidae.
Lastly, we found a small amount of jellyfish (Pandeiidae spp.), Polybius henslowii swimming crabs and Pleurobrachia pileus comb jellies in the samples (Table 6).

Sound scattering layers internal variability
The results of the EM clustering of echo-integration cells within the layers sampled by the nets are presented in Fig 3. The surface layer was split into two vertically stratified clusters (Fig 3(a)). The S v (f ) spectra of both cluster displayed a similar shape above 100 kHz with a decrease in the overall backscattering intensity with depth. In contrast, both spectra differed between 38 and 95 kHz, with a single broad peak around 80 kHz for cluster 2, and two peaks at 38 and 80 kHz for cluster 1. Multinet U0354 and MIK net tows each sampled exclusively one cluster, cluster 1 and 2, respectively. However, the multinet U0341 sampled clusters 1 and 2 equally.
The deep layer was separated into three vertically stratified clusters, whose mean backscattering intensity increased with depth (Fig 3(b)). The median S v (f ) spectrum of each cluster Characterization of sound scattering layers in the Bay of Biscay using broadband acoustics, nets and video displayed specific shapes: a peak at 38 kHz then a rise for cluster 3; a flat spectrum for cluster 2 and a bump between 18 and 150 kHz for the cluster 1, followed by a flat curve beyond 150 kHz. Most net tows sampled the deepest cluster (1), with the exception of the W-shaped MIK U0352 tow that sampled the whole layer and the multinet U0341 that sampled the interface between cluster 1 and cluster 2.

Forward approach
Predicted and measured S v (f ) curves of each layer are presented with their confidence intervals in Fig 4. We also represented the narrowband equivalent S v spectra (small crosses on the spectrum at 18, 38, 70, 120, 200 and 333 kHz). The predicted S v (f ) spectrum in the surface layer (Fig 4(a)) was dominated by small siphonophores (mean pneumatophore size 0.27 mm) and limacina pteropods. The modeled S v (f ) of  (Fig (a) for surface layer, (c) for deep layer). Coloured areas: confidence intervals (5 th and 95 th percentiles of uncertainty analysis results). Measured S v (f ) (red or blue lines) and predicted S v (f ) (black line) frequency spectra are represented in the right column, (Fig (b)  Characterization of sound scattering layers in the Bay of Biscay using broadband acoustics, nets and video small siphonophores displayed a widely spread resonance peak between 38 and 150 kHz. At low frequencies (18-160 kHz) the modeled siphonophore contribution roughly followed the measured S v (f ) shapes although 3-4 dB lower in level. However, the peak at 38 kHz detected in the measured S v (f ) was not observed in the modeled S v (f ). Futhermore, the predicted S v (f ) was 5 dB higher than the measured S v for frequencies greater than 225 kHz due to the high contribution of pteropods (Fig 4(b)). Of note, narrowband frequency responses did not capture peaks at 38 and 80 kHz, which would have potentially lead to misinterpretation, had narrowband acoustics been used alone.
The predicted S v (f ) of the deep layer (Fig 4(c)) was dominated by pteropod and large siphonophore (mean pneumatophore size 1.23 mm) backscatter. The siphonophore and fish resonant peaks occurred between the 18 and 38 kHz narrowband point. Importantly, predicted and measured S v (f ) were close at frequencies below 38 kHz. At higher frequencies, we observed large discrepancies between both curves, with non-overlapping uncertainty. Indeed, the frequency content between 38 and 250 kHz on the measured frequency spectrum, similar to those produced by the small siphonophores in the S v (f ) surface layer (Fig 4(a)), could not be reproduced based on our biological sampling data. Notice that in this case, narrowband and broadband frequency spectra displayed similar shapes.

Discussion
In this paper, the composition of two SSLs detected in the Bay of Biscay in springtime was characterized by jointly analyzing acoustic and groundtruthing data. In doing so, fine scale heterogeneity within SSLs was evidenced and its importance for sampling was highlighted. The SSLs were generated by a variety of scatterers belonging to several taxa. Their frequency responses from 18 to 100 kHz were dominated by GB siphonophores. As reported in previous studies [1,2], the resonance frequency of small GB organisms fall within the range of frequencies commonly used in fisheries acoustics, namely at lower frequencies . Our study confirms that small GB organisms can dominate the total backscattering at these frequencies, even at low density. Pteropods were the main contributors at higher frequencies (> 100 kHz), with very high densities recorded in the surface SSL (4000-5000 ind.m −3 ). To our knowledge, this is the first mention of such high pteropod density.
The unsupervised clustering of broadband echograms revealed a high frequential heterogeneity within the SSLs. SSLs were stratified vertically, with mild horizontal heterogeneity at our sampling scale. Interestingly, the surface layer was composed of two sub-layers indicating a potential shift in the siphonophore size composition. These two sub-layers were sampled differently by the nets. Two sub-layers out of three were sampled in the deep layer: one had a flat S v (f ) spectrum that could be attributed to gaseous scatterers, and the other showed a bump in its S v (f ) curve which not be assigned to any known scatterer. Unfortunately, these SSL frequency heterogeneities could not easily be assessed during the survey by simply looking at the pulse compressed echograms displayed in real time. Hence, net sampling targeted at best the larger scale features of SSLs. The lack of real-time information on these SSL fine scale frequency heterogeneities thus hindered the characterization of their fine scale composition.
Despite some uncertainty in the net positioning in the surface layer, the measured and predicted S v (f ) matched reasonably well in the forward approach, when considering the uncertainties around the two median curves. The two dominant scatterers producing these SSLs were likely siphonophores and pteropods. However the size diversity of siphonophores is suspected to be greater in the water than in the nets, as a peak at 38 kHz in the measured S v (f ) remained unexplained by the catches. Additionally, some of the mismatch observed at higher frequency (> 200 kHz) might be explained by complex scattering [52] due to the high density of pteropods, or by a non-random distribution of pteropods in their position or angle (e.g. because of aggregation due to mucus or opercular opening in the beam direction).
The physical vertical gradients were weak in the study area, and the surface backscattering at 18 kHz was low (-77 dB) relatively to higher frequencies. We therefore neglected the potential backscattering caused by physical phenomena such as microstructures [2] in the surface layer, as microstructures normally present a high backscattering at 18 kHz [53], and a relatively stable frequency response at higher frequencies [2].
In the deep layer, discrepancies between predicted and measured S v (f ) curves could not be explained by the uncertainty in model parameters. Those differences likely resulted from the under-sampling of important scatterers. The resonance of larger siphonophores and fish juveniles was positioned between 18 and 38 kHz, while the spectrum bump at higher frequency could not be explained by the biological sampling. This could be caused by the presence of unsampled small siphonophores, or by a hypothetical accumulation of centimetric flocks of suspended matter [54] in relation with the high turbidity ( S3 Fig in the Supporting Informations).
Moreover, Argentina sphyranea juveniles were modeled as FL organisms. The timing of Argentina sphyranea juveniles swimbladder development is however unclear [51,55]. Therefore, the modeled S v (f ) of Argentina sphyranea juveniles might have been largely underestimated at lower frequencies if a functional swimbladder was in fact present.
The forward approach required accurate acoustic backscattering models of dominant scatterers and a good parameterization of these models. The uncertainty analysis enabled to represent the potential errors on model parameters and better assess the discrepancies between predicted and measured S v (f ) curves. However, [42]'s model, which was used here, is designed for ka < 0.1, and is less reliable for a higher ka. This decrease in GB model precision at higher frequencies has low impact on the overall results, as GB backscattering was generally low at higher frequencies relatively to other scatterers. Nonetheless, it could be necessary to supplement [42]'s model with a gaseous cylinder model (e.g. [56]) when GB organism backscattering dominates the entire frequency spectrum.
The largest source of bias and uncertainty however likely resided in the biological sampling, due to avoidance [13] and fragile organism destruction by the nets issues [57]. Biological sampling errors were however not included in our uncertainty analysis, due to difficulties to quantify them. One way forward might be to use the discrepancies between the measured and predicted S v (f ) curves to assess the catchability bias of the nets, when measured S v (f )'s are higher than predicted ones [13].
Our results illustrate the potential of broadband acoustics to (i) fill-in the gaps between narrowbands, (ii) investigate the SSLs fine frequential heterogeneity and (iii) ascertain the presence or absence of resonance peaks [21] and derive their actual maximum backscattering intensity (Fig 4(b)). Narrow-band acoustics can lead to serious miss-interpretation when resonant peaks are present far away from the sampled frequency (Fig 4(d)).
Siphonophores, especially smaller ones, were difficult to sample, as they are partially destroyed by classical plankton nets [58]. This difficulty (here, one net out of two sampled them) can lead to an underestimation of pneumatophore size diversity and density. The use of ZooCAM imaging [35] was essential for detecting their presence and counting the small siphonophore pneumatophores, too small to be seen on in situ video recordings. The in situ video recordings were essential to ascertain the presence and count large siphonophores [59].
Importantly, the presence of physonect siphonophores was confirmed in the Bay of Biscay [24,60], as well as their potential strong backscattering intensity from 18 to 120 kHz. This finding could explain the frequent occurrence of widely spread SSLs displaying resonancelike peaks in the Bay of Biscay during springtime ( [23,24,27], whose origin remained unexplained.
In conclusion, the combination of broadband acoustics, nets, imagery and video presented in this paper paves the way for the characterization of composite SSLs produced by complex assemblages of scatterers, of various taxa and size class. This opens up a perspective for acoustic monitoring of the spatio-temporal dynamics of well identified SSL, which might be an essential step in a general comprehension of pelagic ecosystems in the context of global climate change.