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

Abstract

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.

About the Authors:

Arthur Blanluet

Roles Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Writing – original draft

* E-mail: [email protected]

Affiliation: Unité Écologie et Modèles pour l’Halieutique, Ifremer, Nantes, France

ORCID logo http://orcid.org/0000-0002-1372-2171

Mathieu Doray

Roles Conceptualization, Project administration, Resources, Supervision, Validation, Writing – review & editing

Affiliation: Unité Écologie et Modèles pour l’Halieutique, Ifremer, Nantes, France

Laurent Berger

Roles Data curation, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing

Affiliation: Service Acoustique Sous-marine et Traitement de l’Information, Ifremer, Brest, France

Jean-Baptiste Romagnan

Roles Data curation, Validation, Visualization, Writing – review & editing

Affiliation: Unité Écologie et Modèles pour l’Halieutique, Ifremer, Nantes, France

Naig Le Bouffant

Roles Data curation, Software, Validation

Affiliation: Service Acoustique Sous-marine et Traitement de l’Information, Ifremer, Brest, France

Sigrid Lehuta

Roles Methodology, Validation

Affiliation: Unité Écologie et Modèles pour l’Halieutique, Ifremer, Nantes, France

Pierre Petitgas

Roles Funding acquisition, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing

Affiliation: Unité Écologie et Modèles pour l’Halieutique, Ifremer, Nantes, France

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–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 (20-50 mm, 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–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 gas-bearing 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.

Materials and methods

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 27th and the 28th 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.

[Figure omitted. See PDF.]

Fig 1. Study area in the north-western Bay of Biscay (average depth: 110 m).

https://doi.org/10.1371/journal.pone.0223618.g001

Acoustic sampling.

R/V Thalassa was equipped with two hull-mounted Simrad EK60 split-beam narrowband (Continuous Wave pulse, CW) echosounders transmitting at 18 kHz and 38 kHz and four hull-mounted Simrad EK80 split-beam echo-sounders operated in broadband mode (Frequency Modulated pulse, FM) transmitting from 47 kHz to 420 kHz in plural upward sweeps (Table 1). The 6 echosounders were set to ping sequentially at 2 pings per second (PPS) (i.e. one ping every 3 seconds in each band) to avoid cross-channel interferences. The transducers were located at 6m below sea surface. Acoustic data collected above 10 m deep were discarded to exclude the 18 kHz echosounder post transmit ringing and to ensure that acoustical data were collected within the far field of each echosounders.

[Figure omitted. See PDF.]

Table 1. Echosounders settings.

https://doi.org/10.1371/journal.pone.0223618.t001

Calibration was performed using the SIMRAD EK80 software (version 1.11.1) with default calibration settings [30] (see details in S1 Appendix). Frequency regions with noise spikes or low backscattering sensibility were removed from any further analysis (e.g. 240-260 kHz, part of 200 kHz and 333 kHz bandwidth when noise peak at depth interfered with the signal, S1 Fig in the Supporting Informations). Due to its limited range, the 333 kHz band was not used in the following cluster analysis.

EK80 broadband raw data were pre-processed using customized MATLAB codes developed at Ifremer to obtain measurements of frequency-dependent volume backscattering strength (Sv, dB re 1m−1, [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 Sv 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 MIK codend to ensure that all organisms entering the net would pass in the camera field of view.

[Figure omitted. See PDF.]

Table 2. Characteristics of the nets used for biological sampling.

https://doi.org/10.1371/journal.pone.0223618.t002

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 Fig in the Supporting Informations.

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.

[Figure omitted. See PDF.]

Fig 2. Pulse compressed Sv echograms recorded during net tows.

The vertical blackline separate data collected on May, 27th and 28th. Blue and red horizontal lines represent the boundaries of surface and deep sampled layers, respectively.

https://doi.org/10.1371/journal.pone.0223618.g002

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 (Sv, 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 Sv(f) were then compared to our in situ Sv(f) averaged in the area where biological samplings were performed.

The predicted Sv(f) (Sv predicted(f)) were calculated for each depth layer as the sum in natural scale of the frequency spectrum Svj(f) of each sampled taxa j (Eq 1):(1) Svj(f) was calculated as:(2)Where TSi represents the modeled Target Strength of length class i in dB re 1m2 and Di the sampled density of the organisms (in individual by m3) belonging to the ith length class (Eq 2). Models and parameters used for each organism are presented in the next section.

The measured Sv(f) (Sv measured(f)) of each layer were calculated as the median of the Sv(f) in echo-integration cells along the sampling net trajectory. These cells were defined as 2 echo-integration cells located below the top of the net. Uncertainty in the predicted and measured Sv(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].

[Figure omitted. See PDF.]

Table 3. Scattering model parameters used for Fluid-like organisms.

SD: Standard Deviation.

https://doi.org/10.1371/journal.pone.0223618.t003

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-to-width ratio (L/2a) of 2.35, regardless of their size. As no specific value of viscosity and surface tension could be found in the literature for siphonophores pneumatophores, those parameters were arbitrarily set to 0.1 Pa.s−1 and 15 N.m−1, respectively. Those values were chosen to lay between values reported for fish flesh (1 Pa.s−1 and 200 N.m−1, from [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 Sv(f). Parameters used in GB organisms scattering models are presented in Table 4.

[Figure omitted. See PDF.]

Table 4. Scattering model parameters used for Gas Bearing organisms.

https://doi.org/10.1371/journal.pone.0223618.t004

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.

[Figure omitted. See PDF.]

Table 5. Statistical distributions for the model parameters used in the uncertainty analysis.

μ: mean, SD: standard deviation, 5% and 95% quantiles of the distribution: Confidence Interval (0.05-0.95) the, L/2a: length/width ratio, L/a length/half-width ratio.

https://doi.org/10.1371/journal.pone.0223618.t005

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 Sv(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 Sv(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 Sv(f) spectrum weighted by each transducer bandwidth, as follows:(3)

with the sv mean value (in natural scale) over the whole bandwidth of the transducer n:(4)

with Svb the backscattering at the b frequency of the transducer n bandwidth, fbandwidth 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) Sv(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 Sv(f) spectra shapes while minimizing the number of clusters.

Results

Biological sampling

The density, mean length and Standard Deviation (SD) of sampled organisms in each zone and layer are presented in Table 6. Illustrations of the sampled main scatterers can be find in the S2 Fig in the Supporting Informations. CTD profiles are presented in the S3 Fig in the Supporting Informations.

[Figure omitted. See PDF.]

Table 6. Composition and length of each layer main scatterers.

SD: Standard Deviation.

https://doi.org/10.1371/journal.pone.0223618.t006

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.

[Figure omitted. See PDF.]

Fig 3. EM clustering results for (a) the surface layer and (b) the deep layer.

Upper panels: echo-integration cells belongings to cluster; left of vertical black line: data collected on May, 27th, right of black line: data collected on May, 28th. Middle panel: net trajectories (black lines). Lower panel: median volume backscattering strength spectra at frequency f (Sv(f)) in the whole depth layer (red and blue curve with 25% and 75% quantiles confidence intervals) and in each cluster.

https://doi.org/10.1371/journal.pone.0223618.g003

The surface layer was split into two vertically stratified clusters (Fig 3(a)). The Sv(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 Sv(f) spectrum of each cluster 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 Sv(f) curves of each layer are presented with their confidence intervals in Fig 4. We also represented the narrowband equivalent Sv spectra (small crosses on the spectrum at 18, 38, 70, 120, 200 and 333 kHz).

[Figure omitted. See PDF.]

Fig 4.

Forward approach results: modeled frequency spectra (Sv(f)) are represented in the left column (Fig (a) for surface layer, (c) for deep layer). Coloured areas: confidence intervals (5th and 95th percentiles of uncertainty analysis results). Measured Sv(f) (red or blue lines) and predicted Sv(f) (black line) frequency spectra are represented in the right column, (Fig (b) for surface layer,(d) for deep layer). Coloured areas around measured Sv(f): confidence intervals (25th and 75th percentiles). Volume backscattering strength Sv level at 18, 38, 70, 120, 200 and 333 kHz and their 25th and 75th percentiles confidence intervals are represented by crosses.

https://doi.org/10.1371/journal.pone.0223618.g004

The predicted Sv(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 Sv(f) of 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 Sv(f) shapes although 3-4 dB lower in level. However, the peak at 38 kHz detected in the measured Sv(f) was not observed in the modeled Sv(f). Futhermore, the predicted Sv(f) was 5 dB higher than the measured Sv 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 Sv(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 Sv(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 Sv(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 (18-50 kHz). 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 Sv(f) spectrum that could be attributed to gaseous scatterers, and the other showed a bump in its Sv(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 Sv(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 Sv(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 Sv(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 un-sampled 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 Sv(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 Sv(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 Sv(f) curves to assess the catchability bias of the nets, when measured Sv(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 resonance-like 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.

Supporting information

[Figure omitted. See PDF.]

S1 Appendix. Echosounder calibrations.

https://doi.org/10.1371/journal.pone.0223618.s001

(PDF)

S2 Appendix. Gaseous target backscattering model.

https://doi.org/10.1371/journal.pone.0223618.s002

(PDF)

S1 Fig. Noise spectrum.

Noise measurements (volume backscattering in decibel, color key) performed with the echosounders in passive mode during the PELGAS2017 survey. x-axis: acoustic frequency (kHz), y-axis: range from transducer (m).

https://doi.org/10.1371/journal.pone.0223618.s003

(TIF)

S2 Fig. Images of main scatterers.

Images of main scatterers in the samples taken with (a), (b) and (c): zooCAM; (d) and (e): Zooscan; (f), (g) and (h): binocular microscope.

https://doi.org/10.1371/journal.pone.0223618.s004

(TIF)

S3 Fig. Hydrological context.

Temperature (a), salinity (b), fluorescence (c) and turbidity (d) profiles. Each dashed line represent one profile, the solid lines represent the mean profile for daytime (red) and night time (blue).

https://doi.org/10.1371/journal.pone.0223618.s005

(TIF)

S1 Table. Description of used mesozooplankton and micronekton nets.

Depth and time of the net tows. Each gear was associated with a station number. The multinet open each of its nets at a specified depth, presented in net depth interval. Only the multinet nets sampling in a sampled layer are presented in this table.

https://doi.org/10.1371/journal.pone.0223618.s006

(PDF)

Acknowledgments

The authors would like to thank the R/V Thalassa Captain Loic Provost, his crew and all the scientists involved in the 2016 PELGAS surveys. We would also like to thank SIMRAD and Lars N. Andersen, for lending the EK80 echosounders. Special acknowledgements to Fabien Lombard for his priceless help on siphonophore identification, to Cyrill Gallut and Nalani Schnell for micronekton identification and Julien Simon for the use of cameras on nets. We want additionally to thank Sven Gastauer, Dezhang Chu and all the anonymous reviewers for all their pertinent comments and remarks over the manuscript.

Citation: Blanluet A, Doray M, Berger L, Romagnan J-B, Le Bouffant N, Lehuta S, et al. (2019) Characterization of sound scattering layers in the Bay of Biscay using broadband acoustics, nets and video. PLoS ONE 14(10): e0223618. https://doi.org/10.1371/journal.pone.0223618

References

1. Mair AM, Fernandes PG, Lebourges-Dhaussy A, Brierley AS. An investigation into the zooplankton composition of a prominent 38-kHz scattering layer in the North Sea. Journal of Plankton Research. 2005;27(7):623–633.

2. Lavery AC, Wiebe PH, Stanton TK, Lawson GL, Benfield MC, Copley N. Determining dominant scatterers of sound in mixed zooplankton populations. The Journal of the Acoustical Society of America. 2007;122(6):3304. pmid:18247742

3. Peña M, Olivar MP, Balbín R, López-Jurado JL, Iglesias M, Miquel J, et al. Acoustic detection of mesopelagic fishes in scattering layers of the Balearic Sea (western Mediterranean). Canadian Journal of Fisheries and Aquatic Sciences. 2014;71(8):1186–1197.

4. Proud R, Handegard NO, Kloser RJ, Cox MJ, Brierley AS, Handling editor: Demer David. From siphonophores to deep scattering layers: uncertainty ranges for the estimation of global mesopelagic fish biomass. ICES Journal of Marine Science. 2018.

5. Tont SA. Deep scattering layers: patterns in the Pacific. Calif Coop Ocean Fish Investig Rep. 1976;18:112–117.

6. Dietz R. Deep scattering layer in the Pacific and Antarctic Oceans. J mar Res. 1948;7(3):430–442.

7. Irigoien X, Klevjer TA, Røstad A, Martinez U, Boyra G, Acuña JL, et al. Large mesopelagic fishes biomass and trophic efficiency in the open ocean. Nature Communications. 2014;5(1). pmid:24509953

8. Davison PC, Checkley DM, Koslow JA, Barlow J. Carbon export mediated by mesopelagic fishes in the northeast Pacific Ocean. Progress in Oceanography. 2013;116:14–30.

9. Banse K. Zooplankton: Pivotal role in the control of ocean production. ICES Journal of Marine Science. 1995;52(3-4):265–277.

10. Beaugrand G, Brander KM, Lindley JA, Souissi S, Reid PC. Plankton effect on cod recruitment in the North Sea. Nature. 2003;426(6967):661–664. pmid:14668864

11. Sieburth JM, Smetacek V, Lenz J. Pelagic ecosystem structure: Heterotrophic compartments of the plankton and their relationship to plankton size fractions 1. Limnology and Oceanography. 1978;23(6):1256–1263.

12. Brodeur RD, Seki MP, Pakhomov EA, Suntsov AV. Micronekton—What are they and why are they important? PICES Press. 2005;13:7–11.

13. Kaartvedt S, Staby A, Aksnes D. Efficient trawl avoidance by mesopelagic fishes causes large underestimation of their biomass. Marine Ecology Progress Series. 2012;456:1–6.

14. Kloser RJ, Ryan TE, Keith G, Gershwin L. Deep-scattering layer, gas-bladder density, and size estimates using a two-frequency acoustic and optical probe. ICES Journal of Marine Science: Journal du Conseil. 2016;73(8):2037–2048.

15. Benoit-Bird KJ, Lawson GL. Ecological Insights from Pelagic Habitats Acquired Using Active Acoustic Techniques. Annual Review of Marine Science. 2016;8(1):463–490. pmid:26515810

16. Stanton T, Wiebe PH, Chu D, Benfield MC, Scalon L, Martin LV, et al. On acoustic estimates of zooplankton biomass. ICES Journal of Marine Science. 1994;51(4):505–512.

17. Davison PC, Koslow JA, Kloser RJ. Acoustic biomass estimation of mesopelagic fish: backscattering from individuals, populations, and communities. ICES Journal of Marine Science. 2015;72(5):1413–1424.

18. Demer DA, Andersen LN, Bassett C, Berger L, Chu D, Condiotty J, et al. USA–Norway EK80 Workshop Report: Evaluation of a wideband echosounder for fisheries and marine ecosystem science. NOAA’s Southwest Fisheries Science Center in San Diego, California, USA: ICES; 2017. 336.

19. Stanton TK, Sellers CJ, Jech JM. Resonance classification of mixed assemblages of fish with swimbladders using a modified commercial broadband acoustic echosounder at 1–6 kHz. Canadian Journal of Fisheries and Aquatic Sciences. 2012;69(5):854–868.

20. Jech JM, Lawson GL, Lavery AC. Wideband (15–260 kHz) acoustic volume backscattering spectra of Northern krill (Meganyctiphanes norvegica) and butterfish (Peprilus triacanthus). ICES Journal of Marine Science. 2017;74(8):2249–2261.

21. Bassett C, De Robertis A, Wilson CD, Ratilal HeP. Broadband echosounder measurements of the frequency response of fishes and euphausiids in the Gulf of Alaska. ICES Journal of Marine Science. 2018;75(3):1131–1142.

22. Doray M, Petitgas P, Romagnan JB, Huret M, Duhamel E, Dupuy C, et al. The PELGAS survey: Ship-based integrated monitoring of the Bay of Biscay pelagic ecosystem. Progress in Oceanography. 2018;166:15–29.

23. Lezama-Ochoa A, Ballón M, Woillez M, Grados D, Irigoien X, Bertrand A. Spatial patterns and scale-dependent relationships between macrozooplankton and fish in the Bay of Biscay: an acoustic study. Marine Ecology Progress Series. 2011;439:151–168.

24. Remond B. Les couches diffusantes du golfe de Gascogne: caractérisation acoustique, composition spécifique et distribution spatiale [PhD Thesis]. Université Pierre et Marie Curie. Institut Français pour l’Exploitation de la MER (IFREMER), Unité Ecologie et Modèles pour l’Halieutiques (EMH); 2015. Available from: http://archimer.ifremer.fr/doc/00267/37784/.

25. Stanton TK, Chu D, Jech JM, Irish JD. New broadband methods for resonance classification and high-resolution imagery of fish with swimbladders using a modified commercial broadband echosounder. ICES Journal of Marine Science. 2010;67(2):365–378.

26. Ainslie MA, Leighton TG. Review of scattering and extinction cross-sections, damping factors, and resonance frequencies of a spherical gas bubble. The Journal of the Acoustical Society of America. 2011;130(5):3184–3208. pmid:22087992

27. Trenkel VM, Berger L. A fisheries acoustic multi-frequency indicator to inform on large scale spatial patterns of aquatic pelagic ecosystems. Ecological Indicators. 2013;30:72–79.

28. Love RH. Predictions of volume scattering strengths from biological trawl data. The Journal of the Acoustical Society of America. 1975;57(2):300.

29. Doray M, Duhamel E, Huret M, Petitgas P. PELGAS 2016 cruise, Thalassa R/V; 2016.

30. Demer DA, Berger L, Bernasconi M, Eckhard B, Boswell K, Chu D, et al. Calibration of acoustic instruments. Denmark: ICES; 2015. 326. Available from: http://ices.dk/sites/pub/Publication%20Reports/Cooperative%20Research%20Report%20(CRR)/crr326/CRR326.pdf.

31. Maclennan ND, Fernandes PG, Dalen J. A consistent approach to definitions and symbols in fisheries acoustics. ICES Journal of Marine Science. 2002;59(2):365–369.

32. Trenkel VM, Berger L, Bourguignon S, Doray M, Fablet R, Massé J, et al. Overview of recent progress in fisheries acoustics made by Ifremer with examples from the Bay of Biscay. Aquatic Living Resources. 2009;22(4):433–445.

33. De Robertis A, McKelvey DR, Ressler PH. Development and application of an empirical multifrequency method for backscatter classification. Canadian Journal of Fisheries and Aquatic Sciences. 2010;67(9):1459–1474.

34. Wiebe PH, Allison D, Kennedy M, Moncoiffe G. A vocabulary for the configuration of net tows for collecting plankton and micronekton. Journal of Plankton Research. 2015;37(1):21–27.

35. Colas F, Tardivel M, Perchoc J, Lunven M, Forest B, Guyader G, et al. The ZooCAM, a new in-flow imaging system for fast onboard counting, sizing and classification of fish eggs and metazooplankton. Progress in Oceanography. 2018;166:54–65.

36. Gorsky G, Ohman MD, Picheral M, Gasparini S, Stemmann L, Romagnan JB, et al. Digital zooplankton image analysis using the ZooScan integrated system. Journal of Plankton Research. 2010;32(3):285–303.

37. Pilcheral M, Colin S, Irisson JO. EcoTaxa, a tool for the taxonomic classification of images; 2017. Available from: http://ecotaxa.obs-vlfr.fr.

38. Benoit-Bird KJ. The effects of scattering-layer composition, animal size, and numerical density on the frequency response of volume backscatter. ICES Journal of Marine Science. 2009;66(3):582–593.

39. Chu D, Foote KG, Stanton TK. Further analysis of target strength measurements of Antarctic krill at 38 and 120 kHz: Comparison with deformed cylinder model and inference of orientation distribution. The Journal of the Acoustical Society of America. 1993;93(5):2985.

40. Chu D, Wiebe P. Measurements of sound-speed and density contrasts of zooplankton in Antarctic waters. ICES Journal of Marine Science. 2005;62(4):818–831.

41. Gorska N, Ona E, Korneliussen R. Acoustic backscattering by Atlantic mackerel as being representative of fish that lack a swimbladder. Backscattering by individual fish. ICES Journal of Marine Science. 2005;62(5):984–995.

42. Ye Z. Low-frequency acoustic scattering by gas-filled prolate spheroids in liquids. The Journal of the Acoustical Society of America. 1997;101(4):1945–1952.

43. Scoulding B, Chu D, Ona E, Fernandes PG. Target strengths of two abundant mesopelagic fish species. The Journal of the Acoustical Society of America. 2015;137(2):989–1000. pmid:25698030

44. Love RH. Resonant acoustic scattering by swimbladder-bearing fish. The Journal of the Acoustical Society of America. 1978;64(2):571.

45. Alexander R. Physical aspects of swimbladder function. Biological Reviews. 1966;41(1):141–176. pmid:5323464

46. Horn MH. Swimbladder state and structure in relation to behavior and mode of life in stromateoid fishes. Fishery bulletin. 1975;73:95–109.

47. Blaxter J, Batty R. Swimbladder “behaviour” and target strength. Rapports et Proces-verbaux des Réunions du Conseil International pour l’Exploration de la Mer. 1990;189:233–244.

48. Faivre R, Iooss B, Mahévas S, Makowski D, Monod H. Analyse de sensibilité et exploration de modèles: application aux sciences de la nature et de l’environnement. Editions Quae; 2013.

49. McLachlan G, Krishnan T. The EM algorithm and extensions. vol. 382. John Wiley & Sons; 2007.

50. Peña M. Robust clustering methodology for multi-frequency acoustic data: A review of standardization, initialization and cluster geometry. Fisheries Research. 2018;200:49–60.

51. Schmidt J. On the larval and post-larval development of the Argentines (Argentina silus (Ascan.) and Argentina sphyraena (Linné)) with some notes on Mallotus villosus (O. F. Müller). In: MEDDLELSER FRA KOMMISSIONEN FOR HAVUNDERSOGELSER. vol. 2 of FISKERI; 1906. p. 20.

52. Gorska N. Evaluation of sound extinction and echo interference in densely aggregated zooplankton. OCEANOLOGICA. 2000;42(3):315–334.

53. Stranne C, Mayer L, Weber TC, Ruddick BR, Jakobsson M, Jerram K, et al. Acoustic Mapping of Thermohaline Staircases in the Arctic Ocean. Scientific Reports. 2017;7(1). pmid:29123176

54. Fromant G. Mesure de Matières En Suspension (MES) dans la colonne d’eau par combinaison de méthodes acoustiques et optiques. Université de Bretagne Occidentale. Institut Universitaire Européen de la Mer (IUM), au Laboratoire Domaines Océaniques (LDO—UMR 6538); 2015.

55. Marshall NB. Swimbladder Structure of Deep-sea Fishes in Relation to Their Systematics and Biology. vol. 31 of Discovery reports. National institute of oceanography ed. University Press, 1960; 1960.

56. Doray M, Berger L, Le Bouffant N, Coail JY, Vacherot JP, de La Bernardie X, et al. A method for controlled target strength measurements of pelagic fish, with application to European anchovy (Engraulis encrasicolus). ICES Journal of Marine Science. 2016;73(8):1987–1997.

57. Knutsen T, Hosia A, Falkenhaug T, Skern-Mauritzen R, Wiebe PH, Larsen RB, et al. Coincident Mass Occurrence of Gelatinous Zooplankton in Northern Norway. Frontiers in Marine Science. 2018;5(158).

58. Benfield MC, Lavery AC, Wiebe PH, Greene CH, Stanton TK, Copley NJ. Distributions of physonect siphonulae in the Gulf of Maine and their potential as important sources of acoustic scattering. Canadian Journal of Fisheries and Aquatic Sciences. 2003;60(7):759–772.

59. Barham EG. Deep Scattering layer migration and composition: observations from a diving saucer. Science. 1966;151(3716):1399–1403. pmid:17817303

60. Pugh PR. The distribution of siphonophores in a transect across the North Atlantic Ocean at 32 N. Journal of Experimental Marine Biology and Ecology. 1975;20(1):77–97.