The Use of Soundscapes to Monitor Fish Communities: Meaningful Graphical Representations Differ with Acoustic Environment

: Many marine animals produce sounds in several phases of their life cycles, either actively or as a byproduct of their activities, such as during mate attraction or when moving. Recent studies of underwater soundscapes have proved passive acoustic monitoring to be a cost-effective, non-invasive tool to understand ecological processes, especially when sampling in adverse conditions or at great depth. Four days of sound recordings at three seamounts from the Azorean archipelago were examined to assess the suitability of different sound graphical representations to characterize different acoustic environments that contrast in the contribution of vocal fish communities. Long-term spectrograms, sound pressure level, spectral probability densities and the Acoustic Complexity Index (ACI) were computed for two shallow seamounts (Formigas and Princesa Alice, c. 35 m) and one deep seamount (Condor, 190 m) using graphics with different time spans. Only in Formigas, which presented the highest occurrence of fish sounds, was it possible to observe temporal patterns of fish vocal activity in the graphical representations. We highlight that habitats with a higher diversity and abundance of sounds are the most suitable targets for these methods, while in locations with a low prevalence of fish sounds a combination of several methods would be recommended.


Introduction
In the marine environment, sound is a more efficient source of information than vision, touch, taste or smell. Sound provides fast (around five times faster than in air), long-distance information to marine animals during day and night [1][2][3]. Indeed, many marine species, including fish, use sound to navigate, select habitats, detect predators and preys, and to communicate [4][5][6][7][8][9]. The acoustic scene (or soundscape) that encompasses biotic, abiotic and manmade sounds is explored by animals as a source of information and may reveal biodiversity, ecosystem dynamics and ecological processes, and habitat health [10][11][12][13]. The urgent need to assess and monitor marine biodiversity has led humans to explore the information within soundscapes, which has driven significant advances in passive acoustic monitoring (PAM) and soundscape analysis [13]. PAM allows cost-effective, longterm studies with a negligible impact in the biota while revealing information on acoustic sources, their amplitudes, the interaction between the different sources and how marine species may interpret and react to sounds [14].
The biogenic component of soundscapes is mostly related to cetacean, fish and invertebrate sounds, and its study can bring insights about the community composition of vocal/acoustic species [10,15,16]. PAM can be particularly useful to monitor populations over long periods of time, detect reproductive aggregations, or follow spatiotemporal activities of soniferous species [17][18][19]. For example, Caiger et al. [19] used PAM to detect spawning activity trends of Atlantic cod over ten consecutive winter spawning seasons in Massachusetts Bay (USA), providing valuable insights into its spawning ecology, conservation and fishery managers. The Atlantic cod, like many other fish species, form dense breeding vocal aggregations often showing site-fidelity, accounting for a major component of soundscapes, especially at the lower frequency range (<2 kHz) explored by fishes [19][20][21][22][23].
Due to the presence and prevalence of different communities, different habitats have distinct soundscapes that are dynamic in space and time [10,[24][25][26]. To extract information from long-term soundscape recordings and capture this dynamic, a wide variety of data reduction (e.g. acoustic indices) and data visualization techniques (e.g. long-term sound pressure level representations, longterm spectrograms, power spectral density percentiles) have been used [12,13,25]. Data visualization methods may be particularly useful to target particular communities or species as they enable rapidly scanning for patterns in the soundscape, and have been increasingly applied to aquatic environments, mainly focusing on marine mammals but also on fishes [23,27]. In turn, data reduction, namely acoustic indices, can provide important information on biodiversity, e.g., the Acoustic Complexity Index (ACI), the Acoustic Entropy Index, or Acoustic Richness [12,28,29]. Nevertheless, soundscape data reduction and visualization remain a significant research area [12,30,31], and whether the same methods apply to, and allow for, comparisons of different ecosystems remain to be ascertained. For example, a given metric or resolution for data visualization may be suitable for a community that is rich in frequency-modulated biological sound events (e.g., cetacean sounds) but inadequate for one dominated by short broadband sounds (such as the clicks of invertebrates), or with few relatively narrow-band low frequency sounds, such as those of fishes [12]. Finding the best approaches for typical communities is needed and especially important when the standardization of marine soundscape measurements is warranted [13].
Despite recent advances in marine soundscape analysis [6,11,14,24,32,33], significant gaps still need to be addressed, including the identification of vocal fish species and their periods of vocal activity, as acoustic indices characterize mostly vocal species during key phases of their lives [10,19,20]. In addition, there is a need to standardize marine soundscape measurements and refine soundscape exploration and visualization techniques [12][13][14].
Here, we aim to characterize the acoustic environment of shallow and deeper seamounts from the Azorean archipelago, focusing on the contribution of vocal fish to the overall soundscape, and to investigate the best visualization methods and acoustic metrics to detect fish communities that may vary in their acoustic presence across different locations and habitats. We highlight the importance of the vocal fish communities to marine soundscapes and put into evidence that long-term visualization methods may only be informative in habitats with a high acoustic abundance and diversity, such as shallow water seamounts. In contrast, at locations where biological sounds are less frequent, visualization methods should have a higher temporal resolution and be complemented with other methods, such as automatic detection or a combination of acoustic indices. We propose that soundscape analysis can be an effective tool for monitoring biodiversity and for conservation efforts but needs to be adjusted to the particularities of the sampled ecosystems.

Study Site
This study was conducted in the Azorean archipelago. The Azores are located in the Northeast Atlantic, about 1600 km from Portugal. The archipelago comprises nine islands and several seamounts of volcanic origin. Seamounts are important biodiversity hotspots for marine life harboring several fish species [34]. Sound recordings from three seamounts were investigated: Condor (depth, 190 m), located at 17 km from Faial Island, was established as a protected area in June of 2010, when most fisheries were forbidden [35]. Princesa Alice (depth, 36 m) is located at 83 km from Faial Island, and until 2016 it had no protection status. Formigas (depth, 35 m) is located at 55 km from Santa Maria Island and has been classified as a natural reserve since 1988 [36] (Figure 1).

Acoustic Recordings
Ecological Acoustic Recorders (EARs) [37,38] were deployed in three seamounts, Condor, Princesa Alice and Formigas, to obtain long-term recordings of the soundscape of these marine ecosystems ( Figure 1) [39]. In Condor, the hydrophone was suspended at about 8-10 m from the sea bed with a floating device, while in Princesa Alice and Formigas both EARs were deployed on the sea bed by scuba-diving.
The EAR is an autonomous system that consists of a Sensor Technology SQ26-01 hydrophone with a response sensitivity of −193.14/−194.17 dB re 1 V/μPa (with some inter-deployments variation), with a frequency response within ±1.5 dB in the range 1 Hz to 28 kHz. A Burr-Brown ADS8344 A/D converter was used with a zero-to-peak voltage of 1.25 V, and a total system gain of 47.5 dB re 1μPa. Sounds were recorded in Formigas for 30 s every 600 s; in Princesa Alice and Condor for 90 s every 900 s. The recordings were down-sampled to 8 kHz to improve processing speed and search for fish sounds in spectrograms up to 4 kHz.

Data Analysis
For the analysis, we considered recordings made during June 2010, one of the months where greater vocal activity would be expected [40,41]. Since we could not examine all the available recordings, to maximize the probability of selecting recording days with more animal sounds the Acoustic Complexity Index was calculated on the overall data set. A sub-sample of four June days exhibiting higher index values was then selected to be closely inspected. Notice that several authors reported a positive relation between the ACI and the number and diversity of fish sounds [11,12,15], but the index is greatly dependent on the selected settings [12]. Therefore, from the available circa 36 h, 15 h and 72 h of recordings, we selected about 5 h, 10 h and 10 h at Formigas, Princesa Alice and Condor, respectively. Additionally, we also analyzed one day in November of 2008 (2.5 h) with the highest abundance of fish sounds at Condor in the period of 2008 to 2012 [42] to highlight intra-site variability.
The ACI algorithm computes the difference between adjacent spectra for each frequency bin and sums these differences within each time clump-for details on the ACI computation see [28,43]. Because biological sounds are usually transient and/or with variable spectral features, higher values of ACI are expected to be associated with higher sound abundance and diversity [44]. In contrast, more continuous geophony and anthrophony sounds are usually represented by lower acoustic complexity values, although random geophonic noises or intermitting anthropic sounds may cause an increase in ACI values [12,28]. The ACI calculations were made for the frequency ranges most used by fish,  Hz and 2000-4000 Hz, using the WaveSurfer software and the plug-in SoundscapeMeter [43] (FFT window size: 512, frequency resolution: 8 kHz, clump of the ACI algorithm: 0.5 s, noise filter: 0).
Acoustical analyses to identify and quantify the number of occurrences of fish sound sequences (a set of several isolated sounds that arises in a repetitive and stereotyped pattern, including from 1 to several sounds, [23]), were conducted by manually scrolling through the recordings and inspecting visually and aurally the spectrograms with Adobe Audition 3.0 software, following Carriço et al. [23]. The sound sequences used are described in Carriço et al. [23]. A 4 day and a 24 h stacked longduration spectrograms spanning 0-2000 Hz were created for each location with the software WaveSurfer (sampling rate: 8 kHz; FFT size: 512; window type: Hamming; 50% overlap). The spectrogram size was set to 256 pixels in the vertical (frequency) axis, and 1800 or 7200 pixels (24 h and 4-day spectrograms, respectively) in the horizontal (time) axis. The spectrogram was computed with a traditional Sparse Fast Fourier Transform (SFFT) method, and with this temporal representation each pixel depicted the mean value of circa 18-38 time bins (depending on the duty cycle: Formigas with 72 min per day, and Condor and Princesa Alice with 144 min per day). The longduration spectrograms were accompanied by detailed plots with the associated occurrences of the sound sequence types. The frequency range used for each sound sequence was in accordance with Carriço et al. [23]. The 24 h spectrogram was accompanied by plots of the variation of sound pressure level (SPL), computed for the same frequency range throughout the day and the variation of the Acoustic Complexity Index. The ACI values represent the average of all clumps for each track for the frequencies 15-2000 Hz. The expanded spectrograms (sampling rate: 8 kHz; FFT size: 512; window size: 64 points; window type: Hamming, 50% overlap) of the most important sound sequence types were also made for each location. Spectral Probability Density plots (SPD) were computed for the same four days to examine the energy distribution along the frequency. These were compared with the power spectra of some of the most abundant fish sounds at each location. The sound pressure level (SPL) and the Power Spectral Density (PSD) were computed and calibrated in the computing environment R [45], adapting the code available by Merchant et al. [46] (SPL: FFT 512, Hann window, 50% overlap, average for each file; PSD: FFT 4096, Hann window, 50% overlap, 1 s average). For each deployment, the correction factor for calibration was calculated using the hydrophone sensitivity, system gain and zero-to-peak voltage of the analog-to-digital converter.

Results
The presence and diversity of fish sounds were tallied for the three seamounts by manual scrolling of a subsample of four days from June 2010. Formigas had the highest abundance of fish sounds, with 93% of the files (5 h of 30 s recordings every 600 s) containing fish sounds, and presenting 15 sound sequence types in total. This was the only location exhibiting choruses. In Princesa Alice, the abundance of sounds was considerably lower, with sounds present in only 37% of the files (10 h of 90 s recordings every 900 s) but with a slightly higher total number of sound sequence types (17). Condor presented the lowest abundance and diversity, with fish sounds in only 9% of the analyzed files (10 h of 90 s recordings every 900 s) and five sound sequence types. The most prevalent sound sequences were #1, #5, #10, #12, #15, and #50 for Formigas; #1, #5, #10, #12, and #15 for Princesa Alice and #1, #4, #5; and #10 for Condor. Less common sounds were categorized as "others". These sound sequences were identified and characterized in Carriço et al. [23].
A In the long-term spectrograms from Formigas (4 days and 24 h, Figure 2a,b), a higher presence of fish sounds, namely #12 and #50, was noticeable in the night period (9 p.m.-6 a.m.), whereas other sound types seem to have a wider occurrence during day time (#15, #10, others). In the expanded spectrograms (6 h, 90 s and 2 s, Figure 2c-e), there is more energy in the frequency band 50-1500 Hz, which is related to the occurrence of fish sounds. In Figure 2e , some of the sound sequence types (#12, #15 and #50) previously described in Carriço et al. [23] are depicted. Sound pressure level follows a similar trend to the spectrogram (Figure 2b, orange line), presenting higher values associated with the occurrence of fish sounds, but the same trend is not obvious for ACI levels ( Figure  2b, grey line). Interestingly, in Formigas seamount, the sound sequence types distribution depicted below the spectrograms (Figure 2a-d) highlights a frequency partitioning of the acoustic environment where sound sequence #50 seems to occupy the lower frequencies and with a distribution mostly at dusk and dawn. Sound sequence #15 seems to be ubiquitous during the day but occupying intermediate frequencies.
In contrast, sound sequence #12 was mostly present at night and occupied the higher frequencies. Sound sequences #1, #5 and #10 were much less abundant, and were recorded at different times of the day, overlapping their frequencies with other sound sequences present at lower frequencies (Figure 2a,b).
The long-term spectrograms for Princesa Alice presented a high between-day variability ( Figure  3a,b) and did not show an obvious daily pattern, likely related with the lower number of fish sound sequences. Although there was an increase in the SPL values associated to the presence of some fish sounds, the general patterns were not so obvious. The ACI levels were higher than in Formigas and did not present any noticeable pattern. Only at the most expanded spectrograms (Figure 3c,d) fish sounds could be detected. The sound sequence occurrence presented a similar pattern to the one observed in Formigas, where sound sequence #12 was mostly present at night, and sounds #1, #5, #10, #15 and "others" had a more ubiquitous distribution throughout the day (Figure 3a,b).   In the Condor seamount, where the EAR was deployed much deeper (190 m, and at 8-10 m above the seafloor), and the total number of sound sequences was much lower than those recorded in shallow waters, the information collected from the long-term spectrograms was much less representative, being ineffective to demonstrate the presence of the fish vocalizations (Figure 4a,b). As observed for the other locations, sounds #1 and #5 were also recorded at different times of the day despite of their low representation and abundance (Figure 4). The SPL at Condor (Figure 4b) was more constant throughout the day than in the shallow seamounts, but it was nevertheless close to or even higher at some times of the day than at Princesa Alice. However, it was not possible to observe the expected SPL increase associated with the presence of fish sounds noticeable in the spectrograms. As in Princesa Alice, the ACI levels did not present any clear pattern. To highlight the intra-site variability, Figure 4e illustrates the diel variation of one day registered at a different year in Condor (2008), depicting a much higher abundance of sounds: 193 sound sequences recorded in the 2 November 2008 versus 31 sound sequences (daily mean, SD = 27.3) recorded in the 4 days of June 2010. On this day (2008), fish sounds, dominated by the fish sound sequence #1, were visible in the spectrogram in the frequencies mostly below 500 Hz. There was also a marked presence of wind/waves, anthropogenic noise and some cetaceans.
The Spectral Probability Density plots in Figure 5 illustrate the contribution of the different sound sources to soundscapes at the three locations for the selected 4 days in June 2010, depicting the different percentiles (1%, 5%, median, 95%, 99%) and the Root Mean Square level (RMS). The SPD from Formigas exhibits three peaks, especially for the lower percentiles (1%, 5%) and the median, which seem to be highly related with the power spectra of the sound sequences #12 (higher amplitude at about 1300 Hz), #15 (higher amplitude at about 430 Hz) and #50 (higher amplitude at about 250 Hz) (Figures 5 and 6). The strongest SPD peak (400-500 Hz) is observed in all percentiles, which may be related with the presence of the sound sequences #15 near and far from the hydrophone.  The SPD from Princesa Alice presents the two most pronounced peaks at about 400 Hz and 1000-1300 Hz, especially at the higher percentiles (95%, 99%), which should have a strong influence of anthropogenic noise, such as boat noise. The power spectrum of the sound sequence #10 (with higher amplitudes at about 200-400 Hz) may slightly contribute for the 400 Hz peak, while the sound sequence #12 (with higher amplitude at about 1200-1300 Hz) may have a smaller contribution for the peak at around 1200 Hz (Figures 5 and 6). Higher SPD percentiles (95%, 99%) highlight the influence of anthropogenic noise at higher frequencies; those anthropogenic sounds can be clearly recognized in the morning of 30 June (see Figure 3a).
The SPD from Condor present peaks at 127 Hz (at percentiles: 99% and RMS), at 200-500 Hz (percentiles: 99% and RMS) and at higher frequencies (1500 Hz and 2000 Hz). Fish sound sequences #1 (with higher amplitude at about 100-400 Hz), #4 (with higher amplitude at 100-200 Hz) and cetacean sounds may help to explain the 127 Hz and 200-500 Hz peaks. The other peaks may be influenced by the sounds from cetaceans and anthropogenic noise, whose presence was observed in the various frequency ranges.

Discussion
Fishes are prolific contributors to marine underwater sound around the world and understanding the relation between fish sound patterns and ecosystem health would be a major advance towards using PAM as a major conservation tool [47]. A first step to understand these patterns is to use meaningful data visualization and data reduction methods to relay adequate information on local ecosystems that may vary in the contribution of biological sound sources. Here, we addressed how temporal resolution affects the information gained from common data visualization methods and acoustic metrics like long-term spectrograms, SPD plots, SPL or the ACI when applied to soundscapes with different contributions from vocal fish communities. We show that, despite focusing only on marine seamounts, regarded as biodiversity hotspots, the chosen temporal frame significantly affects the retrievable information on vocal fish communities from the applied methods, depending on the abundance of fish sounds.
The selected locations are three contrasting habitats with distinct characteristics: Formigas is a shallow seamount, and has been a protected area since 1988 [36]. Princesa Alice, with a similar depth, is an oceanic seamount, and a famous fishing and recreational diving spot without any protection status at the time of the recordings. Condor is the deepest seamount (190 m) and due to its difficult access is a location that could greatly benefit from being monitored with PAM. Condor has also been covered by a protection status since June 2010 to enable research and monitoring activities [35]. All these particularities should be reflected on the local fish communities and associated biophony. According to Carriço et al. [23,42], acoustic signals produced by fishes on these seamounts can inform us about local biodiversity, presenting similar trends to those obtained with other surveying methods like fishing [42] or visual censuses [48]. Here, we investigated whether graphical representations of soundscapes may inform us about vocal fish communities at these locations.
Formigas presented a notable example of a location where PAM can provide valuable information on the acoustic activity patterns of fish communities, particularly using long-term spectrograms, and SPL and SPD plots. In addition, any of the chosen temporal frames were suitable to illustrate the contribution of fishes to the soundscape. Spectrograms, were useful to reveal a partitioning of fish sounds in frequency and time (Figure 2), clearly indicating that different putative species (possibly associated with different sound sequence types) may vocalize at different periods and exploit different frequency bands. The variation in SPL values was associated with the occurrence of fish sounds observed in the spectrograms. Also, the higher SPL values (higher than those recorded in the other studied locations) are likely associated with the overall stronger biophony at this seamount. This suggests that passive acoustic monitoring associated with long-term soundscape graphical representations are invaluable tools for the soundscape assessment in these habitat types (with a high abundance of vocal fish species likely resulting from a 32-year protection status [36]). Note that higher SPL values may also be influenced by louder sources or reflecting animals being closer to the hydrophone.
In contrast, the Princesa Alice seamount (with no protection status at the time of the recordings) had a lower number of fish sound sequences, and did not exhibit robust patterns with the visualization methods used to inspect the vocal fish community, suggesting that monitoring this habitat with long-term spectrograms, SPL or SPD plots is ineffective. Although some of the sound sequences were similar to Formigas, their lower abundance resulted in less obvious patterns in the spectrograms and lower SPL levels. Only in the expanded spectrograms (270 s and 2 s) could fish sounds be visible. Also, fish sounds were only discreetly present in the lowest percentiles (1%, 5%, median) of the SPD.
Condor presented the lowest fish biophony and SPL levels close to those recorded in Princesa Alice. This was observed even in November 2008 (a month with a known high fish vocal activity [41]), where there was a pronounced presence of two fish sounds (#1 and #10), with a marked daily pattern (mostly present at dusk, night and dawn) (Figure 4e). The recordings from these two months showed a temporal variability for the same location [42]. However, although the abundance of sounds was much higher in November 2008, it was still below the values found for Formigas. Similar to Princesa Alice, the different visualization methods did not offer robust tools to monitor fish communities. Even if there may be a smaller number of sounds produced by vocal fish species, which may be related with the fact that greater depth accommodates less diversity, an important note to consider is that the hydrophone deployed at Condor was placed at about 8-10 m above the bottom, likely hindering the detection of low-intensity sounds from vocal benthonic fishes [49]. The propagation distance of the sounds produced by some fish species may vary from a few centimeters to several meters [50][51][52], which, depending on the species, makes the distance from the emitter to the hydrophone a crucial factor in the efficiency to monitor vocal fish. On the other hand, in Formigas and Princesa Alice, the hydrophone was deployed on the seafloor, allowing better recordings of benthonic species which could have caused additional differences beyond those inherent to the distinct community and habitat type. Consequently, the abundance of fish sounds and the distance of the hydrophone to the source seem to be two determinant factors to allow the detection and evaluation of the presence of vocal species in long-term recordings.
The different calling activity among different locations may have been influenced by different factors, such as water depth, food availability, spawning opportunities and protection status, which are determinant for the spatial distribution of the vocal fish species. Although it was not possible to identify any of the source species responsible for the sound sequences found in the seamounts, all the sequences exhibited general characteristics common to fish sounds such as presenting low frequency (< 2 kHz), little or no frequency modulation, short duration, and temporal patterns within the range of other reported fish sounds [23,41,53]. For example, the sound sequence #12 presented many similarities with the sounds recorded in the Mediterranean belonging to the Scorpaena genus [23,54], which has some species present in Azores.
Of the different visualization techniques, long-term spectrograms seem particularly informative to study ecological processes as they allow the observation of a temporal-spectral acoustic structure that may contain relevant ecological information [55]. For example, Parsons et al. [56] note that the separation of the acoustic signals in frequency and time observed in fish choruses [53], following the acoustic niche hypotheses [25], may help to increase the detectability of sounds in the presence of other competing biological or anthropogenic sounds [57,58] that are increasingly present nowadays in the marine soundscapes. This may provide an adaptive advantage to these species, since, by being more effective in their communication, they should increase their reproductive success [59]. Additionally, Kaplan et al. [60] state that the assessment of coral reefs sounds through long-term acoustic recordings enabled the observation of a relationship between the biological sound production and fish density, and coral cover, which reinforces the idea that the acoustic data can be used as an approximation for the presence and abundance of fish species, as well as to inform and bring insights about habitat features [10,15].
The utility of a commonly used acoustic index, the ACI, was also assessed. The ACI proved to be ineffective in reporting the inherent complexity of fish biophony. Although in Formigas the values were higher, in accordance with its high fish biophony, it showed no clear pattern except a slight increase in the index values associated with the presence of some fish sounds at the end of 4 June (Figure 2b), but it was not representative of fish sound patterns in the remaining studied period. In the same way, for the other two places, it was not possible to see any obvious relationship between the ACI and fish biophony. Bohnenstiehl and colleagues only observed an increase of the ACI associated to an increase of fish calling rates when using a finer spectral resolution in the analysis but obtained opposite results when using coarse spectral resolutions [61]. Likewise, Lindseth and Lobel [62] also state that the index may overlook finer details in the presence of one dominant vocalization. Although, the ACI had been successfully applied in several studies in the marine environment [11,63,64], care is still needed when making comparisons between different locations or with other studies, since they are highly dependent on the chosen settings [12,61] and may be driven by differences in biological sources, weather conditions and anthropogenic activities [65].
Is PAM only useful in locations with a very strong fish biophony? Parsons et al. [47] raise the enquiry whether the metrics and techniques used should be environmental or species-specific. We additionally suggest that it is mostly important to adapt the soundscape analysis methods to the level of fish biophony to extract meaningful information. Therefore, the best approach to be used at sites with low fish biophony is to associate short-term spectrograms (e.g., Figures 2d,e, 3d,e and 4c,d) of selected hours of the day and combine them with other techniques like manual annotation (not possible with large data sets) or automatic recognition of fish sounds [66,67], and the use of a set of several acoustic indices and/or other acoustic data reduction methods-on this, see the work of [30]. These additional techniques can be especially useful for data reduction applied to long-term acoustic recordings. Automatic recognition methods enable fast detection of underwater bioacoustic signals produced by marine animals. These methods can use machine learning to find a certain acoustic pattern similar to the sound of interest from a set of various acoustic recordings [68]. Using a different approach, acoustic indices characterize the distribution of the acoustic energy in the recordings in a single value and can be used as approximate measures of species abundance and diversity [30,38].

Conclusions
This study contributes to highlight the usefulness of some currently used graphical representations of acoustic data to study different ecosystems, and how much information they can relay according to the biological sound contribution to the soundscape type and characteristics of the habitat. Here we show that, for locations with a high abundance and diversity of fish sounds, PAM associated with visualization methods for soundscape analysis is a useful tool to assess and monitor fish communities. However, at other locations where sound abundance is lower, the analysis of longterm spectrograms, SPD, SPL or ACI may not be by themselves the best approach. Thus, visualization techniques must be complemented with other methods, such as automatic recognition of fish sounds, [66][67][68][69] or with the combined application of several acoustic indices [29,31,47] to avoid missing or underrepresenting some important acoustic sources.