© 2005 International Council for the Exploration of the Sea
Acoustic detection of a scallop bed from a single-beam echosounder in the St. Lawrence
a Institut des Sciences de la Mer de Rimouski, Université du Québec à Rimouski 310 allée des Ursulines, Rimouski, Québec G5L 3A1, Canada
b Department of Fisheries and Oceans, Institut Maurice-Lamontagne 850 route de la Mer, Mont-Joli, Québec G5H 3Z4, Canada
*Correspondence to Y. Simard: tel: +1 418 723 1986, ext 1563; fax: +1 418 724 1842. e-mail: yvan_simard{at}uqar.qc.ca.
Single-beam seabed echoes combined with epi-macrobenthos photographs were used to remotely detect a scallop bed and characterize the specific acoustic signal of Iceland scallop (Chlamys islandica). A dense scallop bed was surveyed in 2002, with a QTC VIEW Series IV acoustic ground-discrimination system (AGDS) connected to a 38 kHz, 7° split-beam SIMRAD EK60 scientific echosounder. In 2003, a 50 kHz, 42° single-beam SUZUKI ES-2025 echosounder was connected to a QTC VIEW Series V AGDS. The QTC VIEW data were analysed with QTC IMPACT following the standard procedures and classified into acoustic classes. Several approaches were tested: unsupervised and supervised survey strategies directed to specific benthic communities. The SIMRAD EK60 seabed volume-backscattering strength (Sv) was submitted to a principal component analysis (PCA), before and after removal of a depth trend, and the scores on the first 10 principal components were classed by a K-means cluster analysis. The same seabed Sv data were submitted to stepwise discriminant analysis whose training data sets were defined with the ground-truth photographs using different groupings: biotope types, community types, and finally scallop-density classes. All the QTC AGDS approaches failed to reveal the scallop bed, community structures, or biotopes. The QTC classifications mimicked the bathymetry with a strong correlation of the acoustic classes with depth. The seabed Sv PCA + K-means approach presented similar depth-dependence, but, the PCA + K-means on the Sv residuals revealed the scallop bed. The discriminant analysis was the best solution for the scallop density with a general classification success rate of 75% and up to 91% for the highest density class. The Sv signature of the scallop bed is presented, and the most discriminant part of the acoustic signal is identified.
Keywords: acoustic ground discrimination, benthic biotopes, habitat mapping, QTC VIEW, remote sensing, scallop bed, seabed classification
Received 27 August 2004; accepted 10 March 2005.
| Introduction |
|---|
|
|
|---|
Seafloor mapping becomes more and more essential for the effective management of the marine environment and especially for fisheries as the exploitation of marine biological resources increases. Physical-, geological-, and biological-resource maps have proved to be essential aids in sustainable management by allowing the monitoring of environmental fluctuations and the estimation of anthropogenic influence on benthic communities and habitats (Siwabessy et al., 1999; Kostylev et al., 2001). However, accessing submarine landscapes has never been easy as sampling gears are remotely controlled and, consequently, often blind. Traditionally, many of the earlier mapping studies relied on physical sampling, for example with grabs, dredges, or both gears. This approach is not only time consuming and costly, but also highly disturbing for the benthic biotopes (i.e. habitats and their associated communities) and provides only separated, discrete data across the study area. Recent improvements in single-beam echosounders, sidescan sonar, swath-bathymetric systems, and signal processing now provide effective tools to explore the seabed as a complement to the physical sampling methods traditionally used to carry out benthic surveys (Kenny et al., 2003). More recently, "acoustic ground-discrimination systems" (AGDS) have been developed to detect the acoustic-reflectance properties of seabeds. Different reflectances are linked to differences in the physical and occasionally biological nature of these seabeds. Indeed, the nature of the bottom echoes is-influenced not only by basic sediment grain-size parameters, sediment sorting, microtopography, sediment density, and porosity, but also by the presence, concentration, and type of benthic fauna and flora (Tsemahman and Collins, 1997; Collins and Galloway, 1998; Bornhold et al., 1999; Hamilton et al., 1999; Kloser et al., 2001; Anderson et al., 2002).
RoxAnn (Chivers et al., 1990) and QTC VIEW (Collins et al., 1996) are two widely used AGDSs for surveying biotopes. They extract shape, energy, or both features contained in the bottom acoustic signals (Siwabessy et al., 2000). RoxAnn systems derive two parameters E1 and E2, respectively, the energy of the tail of the first echo, often related to the roughness of the sediment, and the total energy of the second echo, often related to hardness of the sediment (Burns et al., 1989). By associating the relationship between E1 and E2 with a substratum type, distinct bottom types can be mapped. It has been shown that RoxAnn performance is highly dependent on ship speed (Hamilton et al., 1999). In contrast, QTC VIEW systems examine the shape characteristics of the first echo and use a series of algorithms to translate it to an array of 166 descriptive variables (Collins et al., 1996), which are then reduced through principal component analysis (PCA) to three Q-values, Q1, Q2, and Q3 (Collins and McConnaughey, 1998). These three Q-values are plotted in three-dimensional Q-space and then run through a cluster analysis to distinguish acoustically distinct bottom types. The QTC VIEW classification accuracy has been shown to be greatly affected by bottom slopes exceeding approximately 58° (von Szalay and McConnaughey, 2002). According to Hamilton et al. (1999), QTC VIEW appears to be the more consistent and reliable of the two systems. Compared with traditional physical sampling, acoustic tools permit rapid, broad-scale, and non-intrusive sampling of the seabed. However, resulting AGDS classifications must be ground-truthed with physical sampling to assess their accuracy. In the so-called "supervised" approach, the classification process is trained a priori on a ground-truth data set.
Most of the acoustic classifications carried out in previous studies focused on identifying and mapping seafloor sediments, which were then associated with benthic habitats and ground-truthed to establish the distribution of organisms on the seabed (Greenstreet et al., 1997; Siwabessy et al., 2000; Kostylev et al., 2001; Anderson et al., 2002; Ellingsen et al., 2002; Freitas et al., 2003a). However, some studies have clearly identified influences from epi-benthic assemblages (e.g. urchin, cockle) in acoustic backscatter (Jumars et al., 1996; Self et al., 2001). Assuming that dense assemblages of benthic fauna and flora could modify the seabed acoustic backscatter, this study aims at detecting the specific signal of Iceland scallop and a known scallop bed located in the St. Lawrence Estuary, Canada. Several acoustic approaches are compared: QTC VIEW systems (Series IV and V) processing using QTC IMPACT software and custom processing of the SIMRAD EK60 scientific echosounder seabed volume-backscattering strength.
| Material and methods |
|---|
|
|
|---|
Study area
This study was conducted off Ile Rouge located in the St. Lawrence Estuary, Québec, Canada (Figure 1). The study area is approximately 22 km2 and the associated Iceland scallop bed is the most upstream one of the St. Lawrence Estuary (Arseneau et al., 2003). This scallop bed is located on Ile Rouge bank, at depths varying between 20 and 60 m. It is bounded by the 250300 m deep Laurentian Channel in the north and the 100 m deep South Channel in the south. The northern part of the scallop bed is located in the Saguenay-St. Lawrence Marine Park. It has one of the highest scallop densities in the Gulf of St. Lawrence. From dredge samples, the mean density of the scallop bed was estimated at 0.53 scallops m2 and the maximal density was 2.58 scallops m2 (MPO, 2000; Hartog et al., 2001). The bed has been exploited by commercial fishers since 1998, including the part located in the Marine Park (Arseneau et al., 2003).
|
Acoustic surveys
Two surveys with different acoustic equipment were run over the study area. The first survey occurred on July 15, 2002 aboard NGCC "F.G. Creed", an acoustic SWATH vessel particularly stable and relatively silent. Nineteen survey lines, orientated north/south, were run at an average speed of 16.7 km h1 (4.6 m s1) (Figure 2a). The acoustic equipment comprised a SIMRAD EK60 echosounder operating at 38, 120, and 200 kHz. Only data from the 38 kHz sounder are presented here. The 38 kHz sounder had a 7° split-beam transducer and was operated with a 1-ms pulse length triggered every second. The navigation system, connected to both the echosounder and the QTC VIEWTM system, was a Differential Global Positioning System (DGPS), providing positional accuracy of ±3 m. A vessel-attitude system (POSMV, Applanix, TX, USA) was connected to the EK60 system, which took the heave into account in bottom-depth measurement. A QTC VIEW Series IV (Collins, 1996) was connected to one quadrant of the split-beam transducer. This system detected, processed, and digitized the raw bottom echo trace read from the transducer cable, and extracted 166 echo-shape features from a five-pings, stacked echo with a series of algorithms for energy and shape characteristics in both frequency and time domains (Collins, 1996; Collins et al., 1996). Preston et al. (2004) defined the five algorithm families used to generate the 166 features as following: (i) "cumulative amplitude and ratios of samples of cumulative amplitude"; (ii) "amplitude quantiles"; (iii) "amplitude histogram"; (iv) "power spectrum", and (v) "wavelet packet transform". Furthermore, conductivitytemperaturedepth profiles were made with an SBE-19 CTD (Seabird Electronics Inc.) during the survey to compute the actual sound speed and absorption coefficient profiles for the analysis of the EK60 data.
|
The second acoustic survey was run in June 2003 aboard the trawler NGCC "Calanus II". Forty transects were surveyed at an average speed of 13 km h1 (3.6 m s1) (Figure 2b). Some tighter lines were also run in areas known as being representative of a particular bottom type or benthic faunal community (Figure 2c). A SUZUKI ES-2025 echosounder operating at the single frequency of 50 kHz with a 42° single-beam transducer and emitting 0.60-ms pulses every 2 s was used. The QTC VIEW system was the Series V, which achieves faster digitization of the signal with a better resolution and dynamic range (Freitas et al., 2003b), and permit classification over a depth range that includes very shallow waters, according to the manufacturer. In addition, Series V raw data are the full bipolar waveform of the echo trace, while the Series IV raw data are the pre-processed set of echo descriptors. The main differences between the acoustic equipments used for the two acoustic surveys are given in Table 1.
|
Ground-truth sampling
Biological data were collected in June 2002 aboard the RV "A.C. Horth". Seafloor photographs were taken with a high-resolution (4.2 megapixels) still camera fitted with a wide-angle lens and placed in a waterproof case and remotely controlled from the vessel. Two 250 W light sources allowed adjustment of the illumination according to the water turbidity, the position of the camera to the bottom, or both factors. The system was mounted on a tetrapod frame that included a reference ruler to evaluate the size of material of the seabed. Throughout the study area, 40 stations, located on a regular grid, but actually sampled slightly away from the plan because the ship drifted with the currents, were surveyed (Figure 2d) and, at each station, five 0.25 m2 images of the seafloor were collected when the camera reached the bottom. Furthermore, Smith-McIntyre grab (0.1 m2 area) and dredge samples were collected. Granulometry of the sediment was done at each station and used to compare grab samples and seafloor photographs (Arseneau et al., 2002). The grab samples underestimated large boulders and rocks, compared with photographs. Seafloor photographs were, therefore, used here since they are more appropriate to identify the relevant biotope characteristics.
Data analysis
The seafloor photographs were analysed through the image analysis software Image Pro Express v4.0 to quantify the surface covered by the sediment and the epi-macrobenthos. A grid of 100 uniformly distributed points was superimposed on the photographs, and what was under each point was identified to give an estimate of the percentages of the surface covered by each component. This is a modified method of that used by Archambault et al. (2001) directly in the field. For each station, these aerial percentages were averaged over the five photographs. For comparison with the acoustic classification, stations were first classed by their dominant organism, which represented more than 50% of the total aerial percentage represented by the macrobenthos. These biological dominance classes were: "scallop", "hydrozoa", "ophiuridae", or "encrusting algae", which were the four most common taxa. When there was no predominance, the station was then assigned to a "non-predominance" class. This allowed the generation of a single map showing the principal biological feature at each station. Furthermore, in the quantification of the epi-macrobenthos, the distinction between live and dead scallops was based on the presence of visible ocelli and a non-eroded shell margin. Second, the aerial percentages were submitted to multivariate statistical analyses to classify the stations by their similarity and to compare these results with the acoustic classification. The matrix of 40 stations by 49 variables, corresponding to the species- and sediment-type aerial percentages, was used to compute similarity matrices on selected variables using the BrayCurtis index (Clarke and Warwick, 1994). The similarity matrices were then submitted to average linkage hierarchical clustering to classify the stations (Legendre and Legendre, 1998). For each resulting classification, the mean relative aerial density and percentage of contribution to within groups similarity were computed for the main components of the group (for details, see Clarke and Warwick, 1994). Third, another classification for comparison with the acoustic classification was based on the aerial percentage of living scallops, which was objectively partitioned in homogenous classes by K-means clustering.
The QTC VIEW full-feature vectors data (FFV) were first analysed using an unsupervised classification method (Bornhold et al., 1999). The 166 variables, resulting from the QTC VIEW echo description, were reduced to three main dimensions (Q1, Q2, and Q3) using principal component analysis (PCA, cf. Legendre and Legendre, 1998) for sorting out the echoes (Collins and McConnaughey, 1998), which generally account for 90% or more of the total variance in a data set (Prager et al., 1995; Legendre et al., 2002). This reduced data matrix was submitted to K-means clustering based on a progressive-splitting process using QTC IMPACT v3.0 post-processing software. The detailed procedure of QTC IMPACT processing is described in Freitas et al. (2003b). The decision to split and merge clusters was assisted by provision of statistical information of each cluster, notably the total score and the Cluster Performance Index rate (CPI). The total score corresponds to the sum of the scores of the individual classes, while the individual score of each class is the product of the number of records and a "chi-square" value (
2) (Ellingsen et al., 2002). The total score should first decrease rapidly and then stabilize when the optimal split level is reached (QTC, 2002b), but see Legendre et al. (2002). The CPI is a measure of the ratio of the distance between the centres of clusters and the extent of the clusters in the Q-space. The CPI rate is defined as CPI(n) = (CPI(n) CPI(n 1))/CPI(n 1), where n is the split number. The CPI rate tends to be maximum at the optimal split level (QTC, 2002b). These two descriptors were taken into account, as recommended by QTC, to decide how to further subdivide the data set. At the end of the procedure, echoes with similar characters formed clusters that defined acoustic classes, which were mapped throughout the surveyed seabed. In addition, as the QTC classification accuracy can be affected by slopes (von Szalay and McConnaughey, 2002), data records corresponding to steep and moderate slopes at the margin of the study area were trimmed off before the analysis. The analysis was also done with all data and the resulting patterns (not shown here) pointed to similar conclusions. Only depths shallower than 70 m were kept and submitted to cluster analysis.
Following the supervised approach, a specific catalogue, which took into account a set of three, pre-determined seabed classes from the biological samples, was created. The FFV data within 150 m from a ground-truth station were assigned to that station and its dominant macro-benthic group ("scallop", "hydrozoa", and "ophiuridae"). These data were then used through QTC IMPACT to create a catalogue of known seabed classes, which served to classify the whole survey FFV data.
A separate acoustic classification was performed from the EK60 scientific echosounder seabed echo (only available for the 2002 acoustic survey). The volume-backscattering coefficient acoustic metric (Sv, m1) used in echointegration (MacLennan and Simmonds, 1992) was estimated for 40 intervals of 75 cm around the echo return from the seafloor to get the seabed echo. The echosounder was calibrated with the standard-sphere method (Foote et al., 1987). The EK60 raw data were first converted into standard HAC format through CH1 software (Simard et al., 1998). The Sv seabed-echo shape and amplitude were computed using CH2 (Simard et al., 2000), which took into account a time-varied gain (TVG) correction estimated from the measured sound-speed and absorption profiles. Thus, each ping was described by a multivariable composed of 40 Sv values. The entire study area represented 22 031 pings, which formed a data matrix of 22 031 rows and 43 columns including latitude, longitude, and depth. The 40 Sv variables were log transformed into volume-backscattering strength (Sv, dB re 1 m1). As for the QTC classification, only depths shallower than 70 m were retained for the analysis. The Sv echo shape was smoothed by a moving average on five pings. Then, this data matrix was submitted to two types of multivariate statistical analyses. First, a PCA was performed on the correlation matrix of the Sv variables. Similarly, these Sv variables were also regressed on depth and the residuals were submitted to a second PCA. Then, the scores on the PCA components representing more than 2.5% (i.e. 1/40) of the total variance were submitted to a K-means cluster analysis. Several acoustic classes were obtained this way. The optimal split level was intuitively determined by the observer as it best matched the biological pattern. The associated acoustic classes were then imported into a GIS and superimposed on the biological results. Second, discriminant analyses (Legendre and Legendre, 1998) were also applied to these seabed Sv data to determine an acoustic-classification solution for a set of pre-determined seabed classes from the biological samples. The pings within 150 m from a ground-truth station were assigned to the macro-benthic class of that station. A new acoustic-station matrix was then generated with the same above 40 Sv variables + 3 columns which included 2105 rows corresponding to the extracted pings. This Sv matrix was subjected to discriminant analysis using a series of pre-determined seabed classes obtained from several different hierarchical and K-means clustering on the ground-truth data set defined above. The classes' assignments were added to the acoustic-station matrix, which then formed the training set for the discriminant analysis. A stepwise analysis was used, in which only those variables which contribute the most to the discriminant function are included, one at a time. The discriminant solutions were then generalized to the whole survey area. The final classifications were mapped and superimposed with biological results.
| Results |
|---|
|
|
|---|
Ground-truth data
From analysis of the photographs, seafloor sediments on Ile Rouge bank are rather coarse with a dominance of gravel (>2 mm) and medium sands (0.250.5 mm). The gravel proportions in sediments decrease downstream from Ile Rouge, from 97% to less than 55%. Big boulders (>40 cm) were common (up to 20%) in the east/southeast of the bank. Scallop shells variably occurred over the study area with a maximum proportion of about 25% in the dredging area. The results of analysis of the photographs are well corroborated by the grab- and dredge-samples analysis, but these observations are not shown here. The sediment inter-sites variability was rather low. Consequently, the sediment pattern itself was not sufficiently variable to compare with an acoustic classification. Several benthic communities were identified from the analysis of the photographs and they tended to be distributed over distinct areas. The most common groups on the bank were ophiuridae, Iceland scallops, calcareous-encrusting algae, and hydrozoa. Two species of Ophiuroidea were observed, Ophiura robusta and Ophiopholis aculeate. They were combined for the analyses. These species were very abundant, especially in the northern and eastern parts of the bank and commonly in the crevices of rocks among scallops or under rocks and shells. They may form dense beds in high-energy, tidal area. They avoid light and prefer to hide beneath rocks and under shells. They were often observed attached to other animals through their adapted, long, flexible arms. Calcareous encrusting algae were observed on the shallower part of the bank relatively close to Ile Rouge, generally fixed on pebbles and boulders, which are rarely overturned. Hydrozoa, which populated the southern part of the bank, were observed on boulders (10 cm and more) and sometimes attached to living scallop shells. These organisms could be very abundant and reach a height >10 cm above the seafloor. These four groups of organisms were retained to map the macrobenthos dominance from the surveyed stations (Figure 3a).
|
Among the several similarity analyses and hierarchical clusterings realized with the ground-truth data matrix, only the ones which gave the best similarity with the discriminant analyses of the acoustic data were retained. First, the hierarchical clustering applied to the similarity matrix computed from the combination of sediments and macrobenthos aerial percentages identified five groups and one singleton that are detailed in Table 2 and mapped in Figure 3b. The five groups are actually mainly sorted by their sediment composition. Some of them can be compared with the previously described macrobenthos spatial distribution (Figure 3a) as follows: group 2 can be related to calcareous-encrusting algae habitat, group 4 to ophiuridae habitat, and group 6 to scallop habitat. Second, hierarchical clustering applied to the similarity matrix computed with only the aerial percentages of macrobenthos provides six groups and three singletons (Table 3 and Figure 3c). Those groups are characterized by one dominant taxon (groups 1 and 2) or macrofauna assemblages (groups 4, 6, 8, and 9) but, in this last case, one taxon is always predominant vs. the other taxa of the community. Thus, scallops are predominant in groups 6 and 9 but are also present in groups 4 and 8. In summary, these two cluster analyses gave consistent results that are both interesting to compare with the discriminant analyses applied to the acoustic data to compare the influence of sediments and macrobenthos combined on the one hand, and only macrobenthos on the other hand. In addition, K-means clustering was also applied to the living scallop aerial percentages only to identify scallop density for a third discriminant analysis with the acoustic data. Several density groups were produced, but only the results for three groups were kept because they best matched our prior knowledge of scallop spatial distribution throughout the survey area. These results are shown in Table 4 and mapped in Figure 3d.
|
|
|
Acoustic classifications
The results of the acoustic classifications of both QTC VIEW systems are given in Table 5. The optimal K-means split level of the QTC VIEW IV data unsupervised classification was already reached at the second split, as the total score began to stabilize and the CPI rate was maximal. Only Q1 and Q2 components contributed to the classification (Figure 4b). The optimal classification solution corresponded to three acoustic classes largely dominated by acoustic class 1. These three acoustic classes are superimposed with the dominance pattern at the biological stations in Figure 4a. This map did not indicate a close correspondence between the acoustic pattern and the main biological assemblages. It rather indicated a correlation of the Q-values with depth, which was significant for the Q1 and Q2 components responsible for the classification (Table 6). For the supervised classification, a specific catalogue was created with the selected training FFV data set for three biological dominances, i.e. "scallop", "ophiuridae", and "hydrozoa". The resulting QTC IMPACT specific catalogue was used to classify the whole survey. This supervised classification, also dominated by acoustic class 1 (Figure 4c), appeared as inconclusive as the above unsupervised classification. No correspondence between these QTC acoustic classification and the biological assemblages was found in this way. The results of the QTC VIEW V data classification up to the third split are presented in Table 7a for the 2003 data set. The total score and the CPI rate indicate that the optimal split level was reached at the second split. The QTC VIEW V system provided a classification pattern very similar to that of the QTC VIEW IV, with the same relationship with the bathymetry (Figure 5a). Here again, the Q-values were highly correlated with depth (Table 8). The results for the specific benthic-community area classification up to the third split are given in Table 7b and Figure 5c. The optimal split level was the second. The first class was dominant again, and the third class contained only a few cases. The resulting spatial distribution of this classification does not match the biological pattern either (Figure 5c).
|
|
|
|
|
|
The PCA results for the Sv depth-filtered data set are presented in Table 9. The first 10 principal components explained 76.4% of the total variance. The reduced matrix of the scores on these first 10 principal components was subdivided into groups by K-means clustering in order to produce acoustic classes. The optimal split level was subjectively determined for K = 4. The resulting acoustic classification did not show any correspondence with the biological pattern (Figure 6a). The map suggested a high correlation of the principal components with depth. This correlation was significant for the 10 principal components. Comparatively, when this depth relation was removed through a regression of Sv variables on depth and the residuals were submitted to PCA, 79.4% of the variance was explained by the first 10 principal components. The scores on these PCA were then separated in groups with K-means. The resulting acoustic classification matched the biological pattern better (Figure 6b). Only the partition with five identified acoustic classes is presented here because this is the one which best fitted the biological classification. The scallop-bed centre appears rather clearly (in red). The few stations with no biological predominance (blue crosses) that are included in this scallop-bed centre were not dominated by scallops but they still hold significant scallop concentrations.
|
|
The results of the three best discriminant analyses are summarized in Table 10, which shows the number of samples correctly assigned to each group (in the column assigned the same number as the row), the number of misclassifications (in other columns), and the good classification rates (in the last column). The attributes (i.e. the Sv bins variables) contributing the most to the discrimination between groups were also identified. The 0.75-m Sv bins were numbered from 1 to 40 according to their altitude, the last bins corresponding to the water column. With a mean sound speed of 1470 m s1 at acquisition, each bin corresponds to 5 x 104 s. First, Table 10a shows the results for the training data set obtained from the similarity analysis and hierarchical clustering on the sediments and macrobenthos percentages combined. Second, Table 10b shows the results for the training data set obtained from the similarity analysis and hierarchical clustering on the macrobenthos aerial percentages. Three of the first most discriminating variables in this analysis are the same as in the previous analysis (variables 40, 35, and 11), and the others are adjacent, which points towards the same part of the seabed signatures contributing to the discrimination. These two discriminant analyses gave comparable results with a rather good success rate, respectively, 69% and 66%. To classify the whole acoustic data set, the results from Table 10a and b were generalized to the whole survey area and the resulting classifications mapped (Figure 7a and b). Both maps highlighted rather well the centre of the scallop bed, from the acoustic classes 6 and 9, respectively, as well as the general spatial pattern of biotope or macrobenthos. Finally, Table 10c shows the discriminant classification results for the scallop-density, training data set. This is the best scallop classification so far. Altogether, 75% of the members from the three subgroups were correctly classified and 91% from the scallop were the richest group. The generalization to the whole survey area is mapped in Figure 7c. This classification probably gives a good idea of the scallop bed extent as defined by high concentration spots.
|
|
The scallop-specific acoustic signature was extracted from the ground-truth training data set comprising the scallop-density classes. Box-plots of the Sv values for each bin are presented in Figure 8 to highlight the seabed echo changes from the lowest to the highest density class to identify the specific "scallop signal" and the "no-scallop signal". The amplitude of the echo largely depends on the beam-pattern of the SIMRAD ES 38B 7° transducer. The peak on each plot (bin 33) corresponds to the backscattering from the surface of the seabed. Before the peak, i.e. just before reaching the seafloor (bin 35), the scallop-signal energy was weaker than the no-scallop signal. The seabed-echo tail is not discriminating between groups until a delay of 6.5 ms (bin 27). As the tail ended its steep decrease, between the time delays for seabed-surface echoe incidences of
2030°, where the most discriminant features are concentrated, the scallop signal held stronger energy content than the no-scallop signal. This situation reverses a little later, at a delay of 9 ms (bin 22), i.e. >30° incidence.
|
| Discussion |
|---|
|
|
|---|
The ground-truth analysis of photographs showed that the inter-site variability of the sediment types was too low to be discriminated acoustically contrary to the benthic-taxa distribution. To a certain extent, however, this substratum pattern determines the benthic-communities distribution as the sediment types (or substrata) are one of the main environmental descriptors that defines benthic habitats and, consequently, influences the composition of the communities. As the species encountered in our study live in preferred sediment types, we assumed that the sediment distribution is implicitly involved in the acoustic discrimination of the benthic communities, but is generally not controlling the observed patterns. Some ancillary factors are also used to characterize a habitat such as water depth, seafloor geomorphology, habitat complexity, current speed, food supply, temperature range, predation pressure, and disturbance by fishing activities (Kostylev et al., 2001, 2003). These environmental factors certainly influence the pattern of scallop distribution and should be taken into account in the interpretation of acoustic classifications. On Ile Rouge bank, scallops were found on areas with strong currents and gravel substratum two factors that are essential for larval dispersion and larval settlement (Bousfield, 1960). The region is also well known as an upwelling region that is very productive (Therriault and Levasseur, 1985) and where concentration processes are very active along topographic features (e.g. Cotté and Simard, in press). This could provide a large source of food for the scallops, and the hydrodynamic regime observed in this area (Simard et al., 2002) could favour the retention of this scallop population in this specific area. The hydrodynamic regime of this region is known to aggregate and retain large biomasses of swimming species such as krill (Lavoie et al., 2000).
The QTC AGDS classifications failed to highlight the benthos distribution. Actually, both QTC-unsupervised classifications were significantly related to bottom depth. But previously, as it was commonly observed that, in general, benthic communities and sediment structure changed with depth. In our case, however, no correlations between sediments or spatial patterns of communities and depth were found in the ground-truth data. Von Szalay and McConnaughey (2002) have shown that even modest bottom slopes could greatly affect the QTC performance. We removed all the data located on steep and moderate slopes from our initial QTC VIEW IV and QTC VIEW V data sets, but the resulting new QTC classifications were still dependent on depth. Similar depth-dependence with QTC AGDS classifications, have been encountered in previous studies. Anderson et al. (2002) demonstrated a high "depth-dependence" in the Q-values for some of their classifications, whereas Foster-Smith et al. (2004) encountered, in shallow waters, a marked depth trend of Q-values (especially Q1). This depth-dependence was also mentioned by Greenstreet et al. (1997) for the RoxAnn system, the performance of which might be affected by water depth, and by Bax et al. (1999), who reported the depth-dependence of the RoxAnn habitat indices E1 and E2. Reporting the case of a study whose first components from the QTC-extracted features were highly correlated with depth, Legendre (2003) interpreted this correlation with depth as being partly related to the QTC IMPACT processing that "only uses three axes in determining its partition [which] makes it especially sensitive to depth". To deal with this specific issue, Legendre (2003) recommended using all the principal components required to explain 95% or more of the variance in the data set. However, Preston (2003) asserted that the QTC IMPACT resampling process to standardize echo duration removes the dependence on depth. Though this may help to align echoes with angles of incidence, the problem of the footprint increasing with depth still remains when the transducer is fixed on the vessel for any acoustic systems, as noted by Morrison et al. (2001) and von Szalay and McConnaughey (2002). The QTC process averages the backscatter signal of five consecutive pings to generate a single-record to reduce ping-to-ping variability that may cloak sediment contribution to the echo (Preston et al., 2004). For the SIMRAD EK60 of the 2002 survey, given the vessel speed and the ping rate, this averaged QTC footprint (one quadrant of the 7° split-beam transducer) was equivalent to
7.3 m2 in 50-m depths with no overlapping between two individual pings. The equivalent for the SUZUKI ES-2025 echosounder of the 2003 survey was
1157 m2 in 50-m depths, i.e. 157 times more than in the 2002 survey, and with a very important overlapping percentage of
95% between two individual pings. Legendre (2003) stressed that a correlation with depth is "not avoidable because the size of the sonar footprint is a function of depth". The relationship between our QTC classifications and depth might be related to the increasing size of the averaged footprint with depth as the acoustic acquisition of this study was made in a relatively wide range of depths (4870 m for the depth-filtered data). Large-footprint sizes also hinder the resolution of small-scale patterns in benthos and bottom habitat, thus limiting the ability of the equipment to accurately discriminate, as mentioned by Kenny et al. (2003). Moreover, the supervised classification within the QTC VIEW IV data processing was not conclusive either, as the procedure failed to recognize the pre-determined subset of known seabed types in the rest of the survey area. The above dependence on depth of the QTC data is not likely to be one of the reasons behind this misclassification because even data in the same water depth were incorrectly classified. Other influences, such as the use of too small a fraction of total variance with only three principal components, or ineffective extraction of proper discriminant features, might be involved.
The bottom Sv data processing through PCA and K-means clustering first showed a similar dependence of the classification with depth in spite of the accurate TVG of the SIMRAD EK60 echosounder corrected for actual sound speed and absorption profiles, which properly compensated for spreading and absorption losses. Yet, a moving average was also performed on five pings to stabilize the signal (Lurton and Pouliquen, 1992; Sternlicht and de Moustiers, 2003). This finally provided an average footprint per ping of
29 m2 at 50 m. The comments made earlier on depth-dependence and footprint-size effects also apply here. Removing the depth trend by a regression was effective for the PCA + K-means clustering, and the resulting unsupervised acoustic classification matched rather well the biological classification. This partitioning of the acoustic variability was done with the first 10 principal components, which represented 79.4% of the variability in the data set. This is less than the 95% recommended by Legendre (2003), but all the significant components (eigenvalue > 1) were included. The statistical removal of the masking trend with depth through a regression produced a better classification, but the acoustic interpretation of this filter and its effect on distorting the seabed-discriminant information are unknown. The unsupervised approach is effective in partitioning the total variability in a study area, but it is not turned to highlight particular habitat or biotope, as our study hoped to do.
The discriminant analyses performed on the bottom Sv data allowed the percentage of classification success and the mean Sv signal of each group to be determined. They also permitted the identification of the variables responsible for the discrimination between the groups in order to help in understanding the relevant acoustic-scattering classification process. They provided the most interesting acoustic classifications. Indeed, the accuracy of this method strongly relies on the presence of discriminant features in the training data set. Several training data sets were tested without success. In descending order, the best correspondence between the resulting discriminant analysis, acoustic classification, and the associated biological classification was found for the three scallop-density classes, then for the six sediment-community classes, and finally for the nine community classes. The number of classes may influence the classification by increasing the overlap between the class signatures. Increasing the range of habitats or biotopes will surely enlarge the possibility of overlapping signatures and consequently reduce the discrimination. The classification accuracy may be improved by grouping habitats or biotopes in broader categories, but with the risk of missing small-scale but significant units (Foster-Smith and Sotheran, 2003). Nevertheless, the information extracted by the discriminant analysis for the three scallop-density classes, with a percentage of correct classification of 91% for the areas that hold the greatest scallop densities, is significant for the objectives of this study. The extent of the scallop bed was thus well defined, revealing patchiness in the distribution with some very high-density spots, essentially in the centre of the bank and extending to the south and west. Actually, commercial-size scallops were very abundant in the surveyed area with a mean dredge density of 0.75 and m2, higher values in the centre of the bed, whereas in the other parts of the St. Lawrence, the highest density was 0.4 m2 (MPO, 2000). This species is characteristic of cold-water regions (Eckman, 1953). It lives in areas with sediments characterized by a mixture of gravel and <10% sand (Fader et al., 1982). Gilkinson and Gagnon (1991) observed a high density of C. islandica associated with coarse sediment (dense gravelcobble). It is generally found attached by byssal threads to coarse substrata in areas with strong currents (Vahl and Clausen, 1980). However, the part of the discriminant information due to the presence of scallop and that attributable to the associated sediment composition remain to be determined. Verifications by additional ground truth at the discontinuities between classes would help validate the generalized classification for the whole survey area. Intense a priori and a posteriori ground-truthing is essential to verify the effectiveness and accuracy of all seabed acoustic-classification methods.
We can speculate on the causes of the specific scallop acoustic signature and the differences with the no-scallop acoustic signature and infer some characteristics of the scallop biotopes. First, there are only slight differences between the acoustic signatures of the seabed classes in this area, probably because the sediment across the whole area is fairly similar and relatively hard and reflective. The main particularity of the scallop signal vs. the no-scallop signal occurs just before the acoustic pulse encounters the seafloor. In areas with scallops, the slope effect on the signal may be less important than in areas without scallop, thus suggesting a flatter bottom for the scallop-rich areas, as the echo starts to return earlier and lasts longer over sloping bottoms. The roughness of the seabed at a small scale might be smaller in high-density scallop areas, or perhaps the interaction between the slope and the roughness of the seabed differs there. The second particularity occurred in the echo tail when the scallop signal appeared stronger than the no-scallop signal. This may be attributable to stronger interface backscattering detected from the first side-lobes of the SIMRAD ES 38B transducer, suggesting a higher roughness in scallop areas detected at incident angles between 20° and 30°, as the sampling surface becomes larger. Alternatively, this may be linked to the seabed volume backscattering in the presence of scallops that significantly differs from the volume backscattering in coarse sediments (e.g. sand) alone (cf. Sternlicht and de Moustiers, 2003). Finally, the third and last difference identified between the two signals occurred beyond 30° incidences, where higher backscatter was detected for areas with few or no scallops. One reason for this result could be that the sampling surface becomes larger, thus considering roughness variability at broader scales, which can then include boulder occurrence. The observation from ground-truth samples that boulders were more common outside the high-density scallop area is in agreement with this hypothesis. Once again, the potential interaction with volume backscattering at angles >30° cannot be overlooked, and may also contribute to the signal patterns.
In the light of these diverse results, we can draw some conclusions about the relative efficiency of each acoustic approach considered in our study. Both QTC systems (Series IV and V) failed to provide conclusive acoustic classifications compared with our prior knowledge of the benthic characteristics of the surveyed area, and this was the case for both sediment pattern, and for the community pattern. QTC results held a strong depth-dependence, which could not easily be filtered out because of the awkward access to the QTC 166 features. However, both QTC systems used in this study have been shown to be effective in mapping other marine benthic habitats (Collins and Galloway, 1998; Hamilton et al., 1999; Preston et al., 1999; Smith et al., 2001; Ellingsen et al., 2002; Freitas et al., 2003b). Furthermore, as suggested by Hamilton et al. (1999), who compared the performance of RoxAnn and QTC VIEW and found not only useful results but also misclassifications with the two systems, AGDSs remain empirical systems that may be efficient for some bottoms and inefficient for others. As an alternative to AGDS, the bottom Sv processing provided several different interesting acoustic classifications within the limits of the methods and equipment as mentioned above. The discriminant-analysis approach is attractive because it can be adapted to the specific objective of the classification looked for, such as sediments, communities, biotopes, or even a particular species, by organizing a posteriori the training data set accordingly. This is indeed not possible with the unsupervised approaches that are used to partition the total acoustic variability observed in a given area into different groups, or with the a priori-defined training data sets. Besides, the discriminant-analysis solution provides the level of uncertainty of the classification, which is not available for unsupervised approaches of the bottom Sv processing. It also identifies the variables that are largely contributing to the discrimination, which is helpful for understanding the related acoustic process. The use of simple acoustic signatures expressed in standard units, such as the Sv used here, is also essential for comparing results with other areas and for generalizing solutions.
In conclusion, the acoustic equipments and processing used in the present study can be of certain efficiency in the remote sensing of the nature of a given seabed. However, as extrapolation to non-surveyed areas is required because the coverage is not 100%, classification uncertainty is unavoidable. AGDS might be expected to detect differences between a limited number of very different biotopes that are well demarcated in space. However, large numbers of subtly different biotopes that merge into each other will be poorly discriminated by such systems (Brown et al., 2001). Given the fairly uniform sediments distribution over the whole study site, the limited resolution of the systems used in this study combined with the partial coverage of the study area might lead to underestimating the small-scale heterogeneity of the benthic habitats. To compensate for the shortcomings of the method using high-resolution acoustics, gears such as sidescan sonar (Brown et al., 2002) or multibeam systems (Kostylev et al., 2003) become more appropriate when looking for a complete coverage of the site. Those systems provide a more precise mapping of the seafloor topography, which could contribute to better estimates of scallop distribution and abundance, as notably shown by Kostylev et al. (2003) with a multibeam. The use of such a system is envisaged for future study of this area. However, the distance between the instrument and the seafloor affects spatial resolution for all systems. A constant-altitude platform is required to obtain a constant-sampling resolution. The distinction between the effect of the seabed slopes and the seafloor composition on the strength of the backscatter at different incident angles is another difficulty that is common to all gears.
| Acknowledgements |
|---|
We thank P. Goudreau for his great help and M. J. Arseneau, F. Hartog, M-F. Brisson, and M. Hardy for their help in the field and photo analyses. Special thanks also to the crews of the RV "Alcide C. Horth", NGCC "F.G. Creed", and NGCC "Calanus II". Comments from a reviewer helped to improve the text. This study was supported by Department of Fisheries and Oceans Canada (Fonds de Recherches Stratégiques en Science de l'Environnement) and the DFO chair in applied marine acoustics at ISMER-UQAR.
| References |
|---|
|
|
|---|
-
Anderson J.T., Gregory R.S., Collins W.T. (2002) Acoustic classification of marine habitats in coastal Newfoundland. ICES Journal of Marine Science 59:156167.
Archambault P., Banwell K., Underwood A.J. (2001) Temporal variation in the structure of intertidal assemblages following the removal of sewage. Marine Ecology Progress Series 222:5162.[Web of Science]
Arseneau M.J., Archambault P., Goudreau P. (2003) Effets de la pêche commerciale sur le gisement de pétoncles d'Islande (Chlamys islandica) de l'Ile Rouge dans l'estuaire du Saint-Laurent: évaluation des impacts sur le pétoncle et la communauté benthique associée. Rapport Technique Canadien des Sciences Halieutiques et Aquatiques, 2512. 43 pp.
Arseneau M.J., Hartog F., Archambault P., Pelletier L., Sévigny J-M. (2002) Évaluation du gisement de pétoncles d'Islande (Chlamys islandica) de l'Ile Rouge: caractérisation des communautés benthiques et évaluation des impacts de la pêche. Report available at Fisheries and Oceans Canada, Institut Maurice-Lamontagne, Mont-Joli, Québec. 41 pp.
Bax N., Kloser R., Williams A., Gowlett-Holmes K., Ryan T. (1999) Seafloor habitat definition for spatial management in fisheries: a case study on the continental shelf of southeast Australia. Oceanologica Acta 22:705719.[CrossRef][Web of Science]
Bornhold B.D., Collins B., Yamanaka L. (1999) Comparison of seabed characterization using sidescan-sonar and acoustic-classification techniques. Proceedings of the Canadian Coastal Conference pp. 893908.
Bousfield E.L. (1960) Canadian Atlantic Sea Shells(National Museum of Canada, Ottawa) 72 pp.
Brown C.J., Cooper K.M., Meadows W.J., Limpenny D.S., Rees H.L. (2002) Small-scale mapping of sea-bed assemblages in the eastern English Channel using sidescan sonar and remote-sampling techniques. Estuarine, Coastal and Shelf Science 54:263278.[CrossRef]
Brown C.J., Hewer W.J., Meadows D.S., Limpenny K.M., Cooper H.L., Rees A.J., Vivian C.M.G. (2001) Mapping of gravel biotopes and examination of the factors controlling the distribution, type and diversity of their biological communities. Science Series Technical Report 114. 43 pp.
Burns D.R., Queen C.B., Sisk H., Mullarkey W., Chivers R.C. (1989) Rapid and convenient acoustic seabed discrimination for fisheries applications. Proceedings of the Institute of Acoustics pp. 169178.
Chivers R.C., Emerson N., Burns D.R. (1990) New acoustic processing for underway surveying. The Hydrographic Journal 56:917.
Clarke K.R. and Warwick R.M. (1994) Change in marine communities: an approach to statistical analysis and interpretation. Natural Environment Research Council, UK. 144 pp.
Collins W.T. (1996) Echosounders used for seabed classification. International Dredging Review 15:1011.
Collins W.T. and Galloway J.L. (1998) Seabed classification and multibeam bathymetry: tools for multidisciplinary mapping. Sea Technology 39:4549.
Collins W.T., Gregory R., Anderson J. (1996) A digital approach to seabed classification. Sea Technology 37:8387.
Collins W.T. and McConnaughey R.A. (1998) Acoustic classification of the sea floor to address essential fish-habitat and marine-protected area requirements. Proceedings of the Canadian Hydrographic Conference '98 pp. 369377.
Cotté C. and Simard Y. (2005) The formation of dense krill patches under tidal forcing at whale-feeding "hot spots" in the St. Lawrence Estuary. Marine Ecology Progress Series 88:199210.
Eckman S. (1953) Zoogeography of the Sea(Sidgwick and Jackson Limited, London) 417 pp.
Ellingsen K.E., Gray J.S., Bjørnbom E. (2002) Acoustic classification of seabed habitats using the QTC VIEWTM system. ICES Journal of Marine Science 59:825835.
Fader G.B., King L.H., Josenhans H.W. (1982) Surficial geology of the Laurentian Channel and the western Grand Banks of Newfoundland. Paper 8122, Geological Survey of Canada. 37 pp.
Foote K.G., Knudsen H.P., Vestnes G., MacLennan D.N., Simmonds E.J. (1987) Calibration of acoustic instruments for fish-density estimation: a practical guide. ICES Cooperative Research Report, 144: 69 pp.
Foster-Smith R.L., Brown C.J., Meadows W.J., White W.H., Limpenny D.S. (2004) Mapping seabed biotopes at two spatial scales in the eastern English Channel. Part 2. Comparison of two acoustic ground-discrimination systems. Journal of the Marine Biological Association of the United Kingdom 84:489500.[CrossRef][Web of Science]
Foster-Smith R.L. and Sotheran I.S. (2003) Mapping marine benthic biotopes using acoustic ground-discrimination systems. International Journal of Remote Sensing 24:27612784.[Web of Science]
Freitas R., Rodrigues A.M., Quintino V. (2003) Benthic biotopes remote sensing using acoustics. Journal of Experimental Marine Biology and Ecology 285286:339353.
Freitas R., Silva S., Quintino V., Rodrigues A.M., Rhynas K., Collins W.T. (2003) Acoustic seabed classification of marine habitats: studies in the western coastal-shelf area of Portugal. ICES Journal of Marine Science 60:599608.
Gilkinson K.D. and Gagnon J-M. (1991) Substratum associations of natural populations of Iceland scallops, Chlamys islandica Müller 1776, on the northeastern Grand Bank of Newfoundland. American Malacological Bulletin 9:5967.[Web of Science]
Greenstreet S.P.R., Tuck I.D., Grewar G.N., Armstrong F., Reid D.G., Wright P.J. (1997) An assessment of the acoustic-survey technique, RoxAnn, as a means of mapping seabed habitat. ICES Journal of Marine Science 54:939959.
Hamilton L.J., Mulhearn P.J., Poeckert R. (1999) Comparison of Rox-Ann and QTC-View acoustic, bottom-classification system performance for the Cairns area, Great Barrier Reef, Australia. Continental Shelf Research 19:15771597.[CrossRef][Web of Science]
Hartog F., Brulotte S., Archambault P., Giguère M. (2001) Gisement de pétoncle de l'île Rouge: caractérisation et évaluation des impacts de la pêche. Saguenay-St. Lawrence Marine Park and Department of Fisheries and Ocean Canada joint project. Internal report, Institut Maurice-Lamontagne, Québec, Canada. 94 pp.
Jumars P.A., Jackson D.R., Gross T.F., Sherwood C. (1996) Acoustic remote sensing of benthic activity: a statistical approach. Limnology Oceanography 41:12201241.
Kenny A.J., Cato I., Desprez M., Fader G., Schüttenhelm R.T.E., Side J. (2003) An overview of seabed-mapping technologies in the context of marine-habitat classification. ICES Journal of Marine Science 60:411418.
Kloser R.J., Bax N.J., Ryan T., Williams A., Barker B.A. (2001) Remote sensing of seabed types in the Australian southeast fishery: development and application of normal incident acoustic techniques and associated "ground-truthing". Marine and Freshwater Research 52:475489.[CrossRef][Web of Science]
Kostylev V.E., Courtney R.C., Robert G., Todd B.J. (2003) Stock evaluation of the giant scallop (Placopecten magellanicus) using high-resolution acoustics for seabed mapping. Fisheries Research 60:479492.[CrossRef][Web of Science]
Kostylev V.E., Todd B.J., Fader G.B.J., Courtney R.C., Cameron G.D.M., Pickrill R.A. (2001) Benthic habitat mapping on the Scotian Shelf based on multibeam bathymetry, surficial geology and seafloor photographs. Marine Ecology Progress Series 219:121137.[CrossRef][Web of Science]
Lavoie D., Simard Y., Saucier F. (2000) Aggregation and dispersion of krill at channel heads and shelf edges: the dynamics in the Saguenay-St. Lawrence Marine Park. Canadian Journal of Fisheries and Aquatic Sciences 57:18531869.
Legendre P. (2003) Reply to the comment by Preston and Kirlin on "Acoustic seabed classification: improved statistical method". Canadian Journal of Fisheries and Aquatic Sciences 60:13011305.
Legendre P., Ellingsen K.E., Bjørnbom E., Casgrain P. (2002) Acoustic seabed classification: improved statistical method. Canadian Journal of Fisheries and Aquatic Sciences 59:10851089.
Legendre P. and Legendre L. (1998) Numerical Ecology 2nd English edn (Elsevier, Amsterdam) 853 pp.
Lurton X. and Pouliquen E. (1992) Automated seabed classification system for echosounders. Proceedings Oceans '92 pp. pp. 317321.
MacLennan D.N. and Simmonds J.E. (2001) Fisheries Acoustics. Fish and Fisheries Series 5(Chapman and Hall, London) 325 pp.
Morrison M.A., Thrush S.F., Budd R. (2001) Detection of acoustic-class boundaries in soft-sediment systems using the seafloor acoustic-discrimination system QTC VIEW. Journal of Sea Research 46:233243.[CrossRef]
MPO. (2000) Etat des stocks de pétoncles des eaux côtières du Québec. MPO Science Rapport sur l'état des stocks 2000/086: xi + 46 pp.
Prager B.T., Caughey D.A., Poekert R.H. (1995) Bottom classification: operational results from QTC VIEW. Proceedings of Oceans '95 pp. 18271835.
Preston J.M. (2003) Re-sampling sonar echo time-series primarily for sea-bed sediment classification. US Patent Application serial number 449914.
Preston J.M., Christney A.C., Beran L.S., Collins W.T. (2004) Statistical seabed segmentation from images and echoes to objective clustering. Proceedings of the Seventh European Conference on Underwater Acoustics.
Preston J.M., Collins W.T., Mosher D.C., Poeckert R.H., Kuwahara R.H. (1999) The strength of correlations between geotechnical variables and acoustic classifications. Proceedings of Oceans '99 pp. 11231128.
QTC. (2002a) QTC VIEW VTM. High-Speed, Data-Acquisition System. User Guide Version 2.1(Quester Tangent Corporation, Marine Technology Centre, Sidney, B.C., Canada) 58 pp.
QTC. (2002b) QTC IMPACTTM. Acoustic Seabed Classification. User Guide Version 3.00(Quester Tangent Corporation, Marine Technology Centre, Sidney, B.C., Canada) 110 pp.
Self R.F.L., A'Hearn P., Jumars P.A., Jackson D.R., Richardson M.D., Briggs K.B. (2001) Effects of macrofauna on acoustic backscatter from the seabed: field manipulations in West Sound, Orcas Island, Washington, U.S.A. Journal of Marine Research 59:9911020.[CrossRef][Web of Science]
Simard Y., Lavoie D., Saucier F.J. (2002) Channel head dynamics: capelin (Mallotus villosus) aggregation in the tidally driven upwelling system of the Saguenay-St. Lawrence Marine Park's whale-feeding ground. Canadian Journal of Fisheries and Aquatic Sciences 59:197210.
Simard Y., McQuinn I., Montminy M., Lang C., Stevens C., Goulet F., Lapierre J-P., Beaulieu J-L., Landry J., Samson Y., Gagné M. (2000) CH2, Canadian hydroacoustic data-analysis tool 2, user's manual (version 2.0). Canadian Technical Report of Fisheries and Aquatic Sciences, 2332.
Simard Y., McQuinn I., Montminy M., Samson Y., Lang C., Stevens C., Miller D. (1998) CH1, Canadian hydroacoustic data-analysis tool 1, user's manual (version 2.0). Canadian Technical Report of Fisheries and Aquatic Sciences, 2256.
Siwabessy P.J.W., Penrose J.D., Fox D.R., Kloser R.J. (2000) Bottom classification in the continental shelf: a case study for the Northwest and Southeast shelf of Australia. Proceedings of the Australian Acoustical Society Conference 6 pp.
Siwabessy P.J.W., Penrose J.D., Kloser R.J., Fox D.R. (1999) Seabed-habitat classification. International Conference on High-Resolution Surveys in Shallow Water 9 pp.
Smith G.F., Bruce D.G., Roach E.B. (2001) Remote acoustic habitat assessment techniques used to characterize the quality and extent of oyster bottom in the Chesapeake Bay. Marine Geodesy 24:171189.[CrossRef]
Sternlicht D.D. and de Moustiers C.P. (2003) Remote sensing of sediment characteristics by optimized echo-envelope matching. Journal of the Acoustical Society of America 114:27272743.[CrossRef][Web of Science][Medline]
Therriault J-C. and Levasseur M. (1985) Control of phytoplankton production in the lower St. Lawrence Estuary: light and freshwater runoff. Naturaliste Canadien 112:7796.
Tsemahman A.S. and Collins W.T. (1997) Acoustic seabed classification and correlation analysis of sediment properties by QTC VIEW. Proceedings of Oceans '97 pp. 921926.
Vahl O. and Clausen B. (1980) Frequency of swimming and energy cost of byssus production in Chlamys islandica (O. F. Müller). Journal du Conseil International pour l'Exploration de la Mer 39:101103.
von Szalay P.G. and McConnaughey R.A. (2002) The effect of slope and vessel speed on the performance of a single-beam, acoustic, seabed-classification system. Fisheries Research 56:99112.[CrossRef]
This article has been cited by other articles:
![]() |
V. Quintino, R. Freitas, R. Mamede, F. Ricardo, A. M. Rodrigues, J. Mota, A. Perez-Ruzafa, and C. Marcos Remote sensing of underwater vegetation using single-beam acoustics ICES J. Mar. Sci., November 13, 2009; (2009) fsp251v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. G. Fernandes Classification trees for species identification of fish-school echotraces ICES J. Mar. Sci., July 1, 2009; 66(6): 1073 - 1080. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



: predominance of ophiuridae;
: predominance of scallops;
: predominance of encrusting algae;
: predominance of hydrozoa;
: no predominance. (b) Spatial distribution of the affinity groups identified by hierarchical clustering on the similarity matrix of macrobenthos and sediments combined: five identified affinity groups and one singleton. (c) Spatial distribution of the affinity groups identified by hierarchical clustering on the similarity matrix of macrobenthos only: six identified affinity groups and three singletons. (d) Spatial distribution of the density groups identified by K-means clustering on living scallops aerial percentages: three identified density groups.




