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

and spatial patterns of low (i.e. 50–1200 Hz) frequency SPLs reflecting fish and lower portion of snapping shrimp acoustic activity at stations (A) 9M, (B) 14M, and (C) 37M in the May River. Time is shown between noon and noon of the next day. Gaps in data = gray, files with physical sounds and anthropogenic noise removed = white, temperature = black line, and day- light hours = dotted line. Green dots indicate first posterior probability (PP) � 0.5 detected during springtime. At station (C) 37M first PP was not calculated for spring 2017 due to miss- ing acoustic data. (TIF) S4 Fig. Time series of broadband frequency Sound Pressure Levels (SPLs) from 2013 until 2018 with physical sounds and anthropogenic noise files removed. Heat maps represent temporal and spatial patterns of broadband (i.e. 1–40,000 Hz) frequency SPLs reflecting all biological activity at stations (A) 9M, (B) 14M, and (C) 37M in the May River. Time is shown between noon and noon of the next day. Gaps in data = gray, files with physical sounds and anthropogenic noise removed = white, temperature = black line, and daylight hours = dotted line. Green dots indicate first posterior probability (PP) � 0.5 detected during springtime. At station (C) 37M first PP was not calculated for spring 2017 due to missing acoustic data. (TIF) S5 Fig. Estimated mean high frequency Sound Pressure Level (SPL) with Posterior Proba- bility (PP) of change at station 9M. Posterior probability of change during springtime (light blue line) with corresponding estimated mean high (7000–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 esti- mated mean SPL due to an increase not a decrease in SPL values) PP � 0.5. (TIF) S6 Fig. Estimated mean low frequency sound pressure level with Posterior Probability (PP) of change at station 9M. 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) S7 Fig. Estimated mean broadband frequency sound pressure level with Posterior Proba- bility (PP) of change at station 9M. Probability of change during springtime (dark blue line) with corresponding estimated mean broadband (1–40,000 Hz) frequency SPL (light 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) S8 Fig. Estimated mean high 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 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. (TIF) PLOS ONE Estuarine soundscapes and biodiversity PLOS ONE | https://doi.org/10.1371/journal.pone.0236874 September 3, 2020 13 / 18