What’s all that racket! Soundscapes, phenology, and biodiversity in estuaries

S9 Fig. Estimated mean low frequency sound pressure level with Posterior Probability (PP) of change at station 14M. Probability of change during springtime (light blue line) with corresponding estimated mean low (50–1200 Hz) frequency SPL (dark blue line) and corre- sponding water temperature (red line) in years (A) 2013, (B) 2014, (C) 2015, (D) 2016, (E) 2017, and (F) 2018. Stars indicate first positive (i.e. detected change in estimated mean SPL due to an increase not a decrease in SPL values) PP � 0.5. (TIF) S10 Fig. Estimated mean broadband frequency sound pressure level with Posterior Proba- bility (PP) of change at station 14M. Probability of change during springtime (light blue line) with corresponding estimated mean (1–40,000 Hz) broadband frequency SPL (dark blue line) and corresponding water temperature (red line) in years (A) 2013, (B) 2014, (C) 2015, (D) 2016, (E) 2017, and (F) 2018. Stars indicate first positive (i.e. detected change in estimated mean SPL due to an increase not a decrease in SPL values) PP � 0.5. (TIF) S11 Fig. Estimated mean high frequency sound pressure level with Posterior Probability (PP) of change at station 37M. Probability of change during springtime (light blue line) with corresponding estimated mean high (7000–40,000 Hz) frequency SPL (dark blue line) and cor- responding water temperature (red line) in years (A) 2013, (B) 2014, (C) 2015, (D) 2016, (E) 2017, and (F) 2018. Stars indicate first positive (i.e. detected change in estimated mean SPL due to an increase not a decrease in SPL values) PP � 0.5. Gray box = no data. (TIF)

S12 Fig. Estimated mean low frequency sound pressure level with Posterior Probability (PP) of change at station 37M. Probability of change during springtime (light blue line) with corresponding estimated mean low (50–1200 Hz) frequency SPL (dark blue line) and corre- sponding water temperature (red line) in years (A) 2013, (B) 2014, (C) 2015, (D) 2016, (E) 2017, and (F) 2018. Stars indicate first positive (i.e. detected change in estimated mean SPL due to an increase not a decrease in SPL values) PP � 0.5. Gray box = no data. (TIF) S13 Fig. Estimated mean broadband frequency sound pressure level with Posterior Proba- bility (PP) of change at station 37M. Probability of change during springtime (light blue line) with corresponding estimated mean broadband (1–40,000 Hz) frequency SPL (dark blue line) and corresponding water temperature (red line) in years (A) 2013, (B) 2014, (C) 2015, (D) 2016, (E) 2017, and (F) 2018. Stars indicate first positive (i.e. detected change in estimated mean SPL due to an increase not a decrease in SPL values) PP � 0.5. Gray box = no data. (TIF) S14 Fig. Time series of broadband frequency Sound Pressure Levels (SPLs) (i.e. with physi- cal sounds and anthropogenic noise files removed), species richness, and abundance from 2016 until 2018. Heat maps represent temporal and spatial patterns of broadband (1–40,000 Hz) frequency SPLs with corresponding species richness (black line), species abundance (blue dotted line), and temperature (red line) at stations (A) 9M, (B) 14M, and (C) 37M. Files with

physical sounds and anthropogenic noise removed = white and no data = gray box. (TIF) S15 Fig. Correlation and regression analysis of species richness and Shannon-Wiener diversity index, sound pressure levels, and temperature. Pearson’s correlation between monthly averages of species richness and monthly averages of sound pressure levels (SPLs) in (A) high (7000–40,000 Hz), (B) low (50–1200 Hz), and (C) broadband (1–40,000 Hz) PLOS ONE Estuarine soundscapes and biodiversity PLOS ONE | https://doi.org/10.1371/journal.pone.0236874 September 3, 2020 14 / 18