Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Photosynthesis by marine algae produces sound, contributing to the daytime soundscape on coral reefs

  • Simon E. Freeman ,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    simon.freeman@navy.mil, simon.freeman@gmail.com

    Affiliations Naval Undersea Warfare Center, Newport, RI, United States of America, Underwater Acoustics and Signal Processing Division, U.S. Naval Research Laboratory, Washington D.C., United States of America

  • Lauren A. Freeman,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Writing – original draft, Writing – review & editing

    Affiliation Naval Undersea Warfare Center, Newport, RI, United States of America

  • Giacomo Giorli,

    Roles Data curation, Investigation, Methodology, Resources, Validation, Writing – review & editing

    Affiliation National Institute of Water and Atmospheric research, Greta Point, Wellington, New Zealand

  • Andreas F. Haas

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – review & editing

    Affiliation Department of Marine Microbiology and Biogeochemistry, NIOZ Royal Netherlands Institute for Sea Research and Utrecht University, Texel, Netherlands

Abstract

We have observed that marine macroalgae produce sound during photosynthesis. The resultant soundscapes correlate with benthic macroalgal cover across shallow Hawaiian coral reefs during the day, despite the presence of other biological noise. Likely ubiquitous but previously overlooked, this source of ambient biological noise in the coastal ocean is driven by local supersaturation of oxygen near the surface of macroalgal filaments, and the resultant formation and release of oxygen-containing bubbles into the water column. During release, relaxation of the bubble to a spherical shape creates a monopole sound source that ‘rings’ at the Minnaert frequency. Many such bubbles create a large, distributed sound source over the sea floor. Reef soundscapes contain vast quantities of biological information, making passive acoustic ecosystem evaluation a tantalizing prospect if the sources are known. Our observations introduce the possibility of a general, volumetrically integrative, noninvasive, rapid and remote technique for evaluating algal abundance and rates of primary productivity in littoral aquatic communities. Increased algal cover is one of the strongest indicators for coral reef ecosystem stress. Visually determining variations in algal abundance is a time-consuming and expensive process. This technique could therefore provide a valuable tool for ecosystem management but also for industrial monitoring of primary production, such as in algae-based biofuel synthesis.

Introduction

Coral reef degradation and algal smothering are caused by human impacts such as overfishing, pollution through nutrient runoff and climate change effects. Algal dominance diminishes the value of reefs [1], is likely permanent [2, 3], and serves as a clear indicator of ecosystem degradation [4]. Timely monitoring of reef state is crucial to quantifying impacts and mitigating ecosystem shifts [5], but present-day assessment methods are resource intensive and infrequent [6]. Because coral reefs are highly soniferous ecosystems where most sound comes from biological processes, ecologists have sought to apply the integrative acoustic monitoring approaches used successfully in terrestrial soundscape ecology [7, 8]. These techniques can yield spatially integrative and long-term continuous observations that can be collected autonomously with inexpensive recorders, but the sources of many reef sounds are challenging to identify [9]. Previously we established positive correlations (Pearsons ρ) between benthic macroalgal cover and aspects of coral reef soundscapes (acoustic pressure spectral density levels in the 2 to 20 kHz band) at 17 sites over a 2570 km transect from Kure Atoll to the Island of Hawaii [10]. These correlations were only observed during the day when overall sound levels were typically lower (Table 1), and correlations of this nature did not exist between sound and any other non-acoustic environmental metric. To further deconvolute the soundscape sources, we performed tank-based passive acoustic experiments with the invasive Hawaiian algae Salicornia gracilaria. Ambient sound was recorded in combination with dissolved oxygen and time-lapse photography of algae under daylight and night conditions to investigate the mechanism behind the observed associations of sound with algal presence.

thumbnail
Table 1. Correlation between benthic algal cover and coral reef soundscapes.

Pearson’s correlation coefficients between the percentage of benthic macroalgal cover at 17 reef sites throughout the Hawaiian Islands and ambient acoustic pressure spectral density in three bands between 2 and 20 kHz, obtained simultaneously at the same locations. Intensity-filtered spectra were averaged over one hour beginning at the indicated times.

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

Short-time Fourier analysis on intensity filtered soundscape recordings [10] provide a more spectrally detailed view of individual transient soundscape components (i.e., most reef-based biological sound). Typical Fourier analysis approaches in underwater acoustics integrate over time periods that are much longer than each individual biological sound. While longer integration results in increased frequency resolution, ideal for detecting tonal and/or narrow-band sources, the approach can also spectrally smear multiple transient sounds together and/or reduce peak level estimates through giving equal weighting to quiet periods between transient arrivals. If a persistent environmental noise is present, it can mask some transient biological sounds sampled through this approach. Our approach involved selecting each transient using an intensity filter, ensuring the transform length encompassed only the transient, and then assessing the spectral qualities of each transient individually. As a consequence, reef spectra shown here cannot be directly compared to spectra obtained without intensity filtering. Depending on ambient and preamplifier noise levels, lower-level sounds from individual biological sources may not be discernable at significant distances from the source. However, many of these events occurring over a distributed area create an extended sound source for which propagation models show decreased attenuation with distance [11].

Results

With the onset of photosynthesis during light periods, bubbles could readily be observed with the naked eye on the surface of macroalgae (Fig 1A and 1B). As the bubbles detached from the plant, they created a short ‘ping’ sound. Acoustic recordings of bubble release separated by quiet periods consequently appear as an irregular pulse-train-like time series (Fig 1C). The sound created by perturbation of the bubble during release is naturally transient and decays exponentially (Fig 1D). Oversampled recordings (100 kS·s-1) permit more detailed spectral analysis of each sound (Fig 1E) and fundamental frequency estimates for each bubble can be made from the resultant data (Fig 1F). Tank resonance modes [12] were restricted through the presence of algae and the drafted, acoustically absorptive nature of the tank walls and did not appear to substantially influence the recorded spectra (Fig 1G).

thumbnail
Fig 1. Characterization of bubble sound from algae.

(A) Gracilaria salicornia actively creating gas bubbles during photosynthesis. White arrows indicate the locations of larger bubbles. The scale bar is 5 mm in length (B) A closer view of one gas bubble about to detach from the algae surface. The scale bar is 2 mm in length. (C) A 9 s time series of ambient sound from inside the aquarium. (D) A high-resolution view of a typical transient waveform as shown in (C). (E) Spectrogram showing the time-varying spectral content of the received waveform in (D). (F) Pressure spectral density (PSD) estimates showing mean, median and maximum pressure spectral densities of received level from the 2 ms period as shown in (D) and (E), indicating the spectral peak of the waveform as 13.07 kHz (vertical dashed line). (G) An averaged spectrum from transient sounds recorded over one hour overlaid with standing-wave resonant tank modes. Mode sums indicate the sum of mode numbers in the horizontal (length, width) and vertical (depth) directions.

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

Experiments with S. gracilaria in controlled laboratory settings demonstrate that sounds produced by algal photosynthesis are similar in nature to components of the soundscape recorded in shallow water regions where S. gracilaria and other benthic macroalgae species are common. With the intermittent application of Photosynthetically Active Radiation (PAR, μ = 381, σ = 67 μmol·m-2s-1, roughly equivalent to the available level at water depths of 10 m on a cloudy day in the equatorial pacific) the link between photosynthesis, dissolved oxygen and acoustic emissions through bubble formation becomes apparent (Fig 2 and S1 Fig). Over the duration of PAR application both the number of bubbles (R2 = 0.76, S2A Fig) and their mean size (R2 = 0.47, S2B Fig) increased with dissolved oxygen levels. The process was reversed with the removal of light. These increases in the size and rate of bubble formation during the illuminated period lower the frequency distribution of sounds produced by the algae (Fig 2A) and cause an increase in the Sound Exposure Level (SEL, R2 = 0.82, Fig 2B and S2C Fig), a measure analogous to acoustic work. The onset of bubbling rate and size increase from the application of light was delayed due to the time required for photosynthesizing algae to build supersaturation conditions in the surrounding water, and for bubbles to reach a size large enough to create sufficient buoyancy force to detach from the algae. After removal of the light source there was a comparable delay in the reduction of bubble production and sound, due to bubbles remaining on the algae and saturation state remaining above 100 percent.

thumbnail
Fig 2. The acoustic response of algae to light over time.

(A) Ten-minute averaged time-frequency histograms showing the distribution of Minnaert frequencies from bubbles produced by algae with the application and removal of PAR. The color scale indicates the number of bubble events acoustically characterized over a 10-min period, per frequency bin (195 Hz bin width). (B) Sound Exposure Level (SEL, 10 min. averages, ± 1 S.E., left axis) and dissolved oxygen concentration (right axis) showing an increase, decrease, and increase of dissolved oxygen and SEL with the application, removal and application of PAR, respectively. The grey regions indicate the period when the light source was removed.

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

A comparison between acoustically derived bubble size distributions and photographically obtained measurements of bubble size (Fig 3) optically validates the acoustic estimation of algae-driven bubble radii (99 percent significance: no evidence the distributions were unequal, Wilcoxon 2-sided test, p = 0.9087).

thumbnail
Fig 3. Size distribution comparisons between bubble radii distributions simultaneously derived through photography and passive acoustic recording.

(A) Bubble size distributions obtained from photographically imaged measurements (histogram, left axis) and Minnaert radii (solid line, right axis) every 150 s over a one-hour period of active bubble formation. (B) Boxplots indicating 25th, 50th, 75th percentiles and 5th and 95th percentile outliers (crosses) of optically imaged and acoustically derived bubble radii distributions over the same period.

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

Discussion

Algae release oxygen as a byproduct of photosynthesis. While these waste molecules are formed intracellularly in solution, the nucleation of oxygen gas bubbles on the surface of macroalgal tissue takes place when localized supersaturation of dissolved oxygen occurs at a nucleation site. Depending on the timescale of bubble formation and the total gas tension, the diffusion of nitrogen will also contribute to total bubble volume. Previous work has shown the relationship between oxygen bubble production by algae and oxygen supersaturation at their surface microenvironment [13]. The bubbles grow with the addition of more waste oxygen produced by photosynthesis and nitrogen diffusion, ultimately separating from the algae through a combination of buoyancy and surface tension forces. As a spherical bubble is perturbed through release, it oscillates in volume and pressure with exponential decay at a frequency inversely proportional to its radius R0, the roots of the specific heat ratio of the gas γ and ambient fluid pressure Pfl, and the inverse root of ambient fluid density, ρfl, a relationship first derived by Minnaert [14] and now referred to as the Minnaert frequency [15] ωMinn: (1)

This approximation was first developed with the assumption of negligible heat flow (adiabatic conditions) and negligible surface tension, but has since been shown to be a good approximation for bubbles of radii between 30 nm to 300 μm. Consequently, passive acoustic estimates of bubble volume can be made if the water properties and depth at which bubble separation occurs are known.

Our observations demonstrate the mechanism behind previously established correlations between components of underwater soundscapes and the relative abundance of macroalgae on Hawaiian coral reefs [10]. Although ecologically distinct, bubble production is ubiquitous in sea grass beds [16], and can also be observed in marine algal species [17]. Thus, the bubble production mechanism is not specific to S. gracilaria and may be used as a general indicator of photosynthetic activity. Ultrasonic emissions from terrestrial plants are also driven by bubble processes, but the underlying mechanism is water starvation rather than photosynthesis [18].

The contribution of algal sound to soundscapes may be isolated from other ubiquitous sound sources though matched filtering, directional receivers and similar techniques. The observed size distribution and thus corresponding spectral distribution of bubble sound (Fig 3A), and the a-periodic and un-clustered nature of acoustic emission time series (Fig 1), are not shared by any other bubble related process in shallow-water environments, including breaking waves [19]. While active acoustic techniques have been used to estimate marine plant biomass [16, 20], our finding introduces the possibility of noninvasive sensing methods that can quantify relative primary production rates in addition to presence/absence, although further work in disentangling confounding factors is necessary. The relationship between bubble formation processes and gas transfer mechanisms has been well studied [13, 2122] providing a background for the application of bubble science to underwater ecological sensing.

The findings presented here show that algae are capable of producing sound under normal circumstances and that sound may explain a correlative association previously discovered between sound in the 2 to 20 kHz band and algal cover on coral reefs [10]. However, a great deal of further work is required before a passive acoustic tool could be developed for quantifying bubble output from photosynthesis, in-situ primary productivity and algal abundance. A number of caveats in our analysis must be considered and results from soundscape correlation and tank experiments should be considered limited in scope until these caveats are resolved. Underwater soundscapes, especially in the vicinity of a shallow coral reef, are extremely complicated acoustic environments containing sound sources of many types [23]. A number of spatio-temporally variable physical mechanisms that influence the production and propagation of sound from source to receiver also exist.

The most notable interference may be snapping shrimp noise. While the sound generation mechanism of these shrimp involves cavitation [24] and is thus spectrally unique from the relatively narrow-band bubble emissions discussed here, source levels and rate of occurrence can be so high [2527] that they are capable of ‘drowning out’ other sound sources for protracted periods of time. Furthermore, the spectral, temporal frequency, and diurnal behavior of snapping shrimp vary depending on the location and time of year [2829, 2627]. While the frequencies produced by photosynthetic bubbles are within the band of sound produced by snapping shrimp, spectral analysis of sound from shrimp-dominated and coral reef soundscapes reveal spectral differences that suggest other biological contributors can be spectrally evaluated even in the presence of shrimp noise. Bubbles produce transient sounds that contain energy in a relatively narrow frequency band. Due to the higher likelihood of smaller bubbles, it is less likely that bubbles would add significant energy at lower frequencies (i.e., below 5 kHz). Conversely, sound produced by snapping shrimp is broad-band but weighted toward a peak between 4 to 6 kHz [25]. The integration of snapping shrimp sounds over time can produce a low-frequency-weighted spectrum with monotonic decay. Conversely, some reef soundscapes may create higher-frequency-weighted spectra that may include distinct spectral peaks (Fig 4A). Limited comparisons between a snapping-shrimp-dominated soundscape and those correlated with algal dominance in Table 1 show a variation in band levels between 5 and 20 kHz between sites. In this case, the shrimp-dominated soundscape is the Scripps Pier in La Jolla, Ca. Periodic cleaning of encrusting sponges on the pilings reveals large communities of snapping shrimp. Care should be taken when evaluating this comparison as the behavior and sound of snapping shrimp may vary between Hawaii and San Diego. However, the Scripps Pier is relatively unique in that no reefs exist within half a mile of the pier. Recordings were made on a calm day with little swell. The comparison shows that in daytime reef soundscape levels were between 5 to 12 dB higher above approximately 5 kHz. Barring differences in frequency-dependent attenuation or shrimp sound characteristics, the comparison suggests an additional source of higher frequency sound present at the reef sites (Fig 4B). The high source levels emitted by snapping shrimp drive the prevailing thought that they are the overwhelmingly dominant bioacoustic source in coastal underwater ecosystems. However, our suggestion that individually quieter biological sources of sound are detectable in the presence of snapping shrimp is not unprecedented. Chorusing from sea urchin grazing is known to contribute to soundscape spectra off the coast of New Zealand [30] and hermit crabs have been spectrally matched to soundscapes from environments in which both they and snapping shrimp are plentiful [9]. Nevertheless, before accurate quantification of algal photosynthetic activity can be attempted in the field, the behavior and acoustic characteristics of snapping shrimp in the area must be understood to create validated acoustic signal processing algorithms capable of differentiating between shrimp noise and algal sound.

thumbnail
Fig 4. Spectral comparisons between Hawaiian reef soundscapes and the snapping-shrimp dominated soundscape at the Scripps Pier in La Jolla, CA.

(A) Normalized, intensity filtered pressure spectral density levels recorded at midday local time during sunny days on shallow Hawaiian reefs (colored lines) plotted with similarly filtered levels obtained during crepuscular chorusing of snapping shrimp communities on the Scripps pier (black line with grey shading underneath). (B) Logarithmic spectral level differences between crepuscular pier and midday Hawaiian reef soundscapes in Fig 4A. Note that 3 dB represents a doubling or halving of the spectral level difference in each frequency bin. Hawaii locations are identified as follows: FFS–French Frigate Shoals; Kure–Kure Atoll; Big I.–Ke’ei Beach, Big Island Hawaii; Lis.–Lisianski Island; Oahu–Lai’e Beach, Oahu, Hawaii; Maui–La Perouse Bay, Maui, Hawaii; PHR–Pearl and Hermes Reef; Kauai–Tunnels Beach, Kauai, Hawaii.

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

Another confounding factor may be the influence of water movement on bubble formation and retention. Bubbles will not form if water is continually swept from the algae and oxygen saturation state never exceeds 100 percent. Secondly, if bubbles are formed, they may be prematurely removed through wave action, decreasing the size distribution of the bubbles and shifting the frequency distribution higher. The mechanisms governing benthic oxygen saturation state and the effects of water movement on bubble retention are poorly understood. Benthic algae reside at least partially within the fluid boundary layer, especially in rugose environments, and may be somewhat protected from flow. The structure of algal filaments may also assist or retard the removal of bubbles by wave action, meaning that aspects of sound production are likely species-specific. While bubble retainment by algae in current remains to be quantitatively investigated, Fig 5 shows typical turf algae in a shallow reef environment off Hawaii’s Big Island on a sunny day. Fig 5 demonstrates that in very shallow water subject to wave action, bubbles from photosynthesis may be retained if sufficient algal structure exists. Further acoustic, flow, and bubble physics analysis, along with an investigation of biological interaction (both floral and faunal) are required to quantify and understand the impact water movement has on bubble formation and the consequent acoustic emissions by algae in these scenarios.

thumbnail
Fig 5. Bubbles from photosynthesis present on shallow algal turf subject to swell.

Note the movement of suspended particles during the camera exposure period of 1/125 sec. Scale bar indicates 2 cm. Image courtesy Florybeth La Valle, University of Hawaii.

https://doi.org/10.1371/journal.pone.0201766.g005

The probability that too many transient sounds arrive at the receiver during one FFT integration period for spectral analysis on individual sounds to be performed increases with the size of the reef, the nature of propagation pathways from sources to the receiver, and the acoustic productivity of the benthic community. As the range between any given source and the receiver increases, the probability that multipath arrivals are interpreted as separate signals also rises. However, the spectral content of multipath arrivals is less likely to differ significantly from the direct arrival when compared to other sound sources. Future data collection with a spatial filter such as a directional receiver or hydrophone array can limit the area over which passive acoustic surveys are performed, simplifying the challenges associated with multipath and the large number of sources.

Sound is the most efficient radiation underwater and we continue to discover that important biological components of many aquatic ecosystems are both sensitive to [31] and produce acoustic emissions. Contributions to biological soundscapes on coral reefs include sound deliberately produced by fishes and invertebrates for communication [32] and defense [25], or inadvertently through feeding [30] and movement [9]. To date these sounds have been used to study the behavior and distribution of marine mammals and fish [3233], to understand ambient noise levels in the ocean [3435], and to assess the impacts of anthropogenic noise on marine environments [36]. Soundscape measurements are also effective at night when benthic activity is greatest [9] and passive optical techniques are difficult or impossible to implement [37]. In concert with these acoustic observations, a volumetrically integrative, remote and noninvasive method of measuring primary productivity outputs may lead to more accurate, persistent and less expensive surveys of ecological state in shallow water, the majority of which are currently performed in-situ through diver-based optical methods. Quantifying primary productivity through passive acoustic monitoring could become a technique that may not be limited to coral reef environments. Macroalgae are the dominant benthic primary producer in many temperate coastal and freshwater environments that are more difficult to survey optically due to reduced water clarity. Passive acoustic monitoring of ebullition has been utilized in the monitoring of chemical reactions [38]. A semi-real-time and volumetrically integrative approach can be adopted in biofuel-generating algal reactors using similar methods. The sensing technique proposed here is restricted neither to in-situ observations nor to marine environments but may apply basically to all aquatic environments in which macroalgae grow.

Materials and methods

Background data

The summary data shown in Table 1 were collected in the Northwest Hawaiian Islands (NWHI) during a NOAA Remote Areas Monitoring Program Cruise HA-12-04 on the R/V Hi’Ialakai during August 2012 and from sites accessible from shore on the Main Hawaiian Islands during September-October 2012. Permits were obtained for conducting research in the Papahānaumokuākea Marine National Monument (PMNM-2012-029) and the Main Hawaiian Islands (Department of Land and Natural Resources special activity permit 2012–83). Acoustic data were collected using a single hydrophone Loggerhead Instruments DSG Ocean recorder configured with a sampling frequency of 80 kSs-1 and an electronic gain of 20 dB. Simultaneous with acoustic recorder deployments, benthic phototransects were taken in accordance with NOAA Coral Reef Ecosystem Division rapid ecological assessment protocol [39]. Pearson’s correlation coefficient, ρ, between all phototransect-derived environmental variables, and band-limited pressure spectral density estimates, were calculated over data from all field sites for which acoustic and phototransect data were available. Metrics for which ρ values were greater than 0.6 with p < 0.001 (Bonferroni correction) were considered to be sufficiently correlated. Correlations of particular interest were those that appeared to form a distinct pattern across several metrics (i.e., similarly high correlation across adjacent acoustic bands, combined with a consistent temporal pattern). A full description of the data collection and processing methods associated with Table 1 is available in the associated publication [10].

Summary of the experiment

Tank-based experiments were carried out at the Hawaii Institute of Marine Biology at Coconut Island, Kaneohe Bay, Hawaii. Ten kilograms (wet) of the invasive red algae Gracilaria salicornia were collected by the Department of Aquatic Resources (DAR) of the State of Hawaii in Kaneohe Bay and stored in a 1 m diameter seawater holding tank with a steady filtered (100 um) ocean water flow. One kilogram of algae was visually inspected to remove any associated fauna and relocated to a smaller, plastic, opaque-sided rectangular tank with internal dimensions of 550 x 300 x 300 mm. Seawater was used to fill the tank to a depth of 250 mm and a 50 mm layer of algae was positioned on the floor of the tank between 200 and 250 mm depth. The algae remained negatively buoyant throughout the experiment. An aquarium light (Radion XR30w G4 pro) was positioned centrally over the tank and 228 mm (9 in.) above the water surface. All testing was conducted in a dark room with the aquarium light being the only source of photosynthetically available radiation (PAR). The color temperature was set at close to that of natural sunlight on a clear day (5500 K), at which the light operated at 55 percent of maximum output. Software limitations required a warm-up and cool-down period of 30 min before and after full-strength application, during which light levels were linearly ramped up and down between zero and 55 percent. Mean PAR level was obtained from measurements at fifteen equally spaced locations at algae equivalent depth using a calibrated LI-COR Li-193 spherical underwater quantum sensor and a LI-COR Li-250 light meter. Note that the mean PAR level obtained (381 μmol·m-2s-1) is far below what is expected near the surface on a calm sunny day in the tropics [40]. The experiment was started with algae being brought out of a 24-hour period of darkness and exposed to PAR as described above. After 3 h and 45 min, light levels were reduced to zero. After 7 h of darkness, light levels were increased back to the levels described above until the conclusion of the experiment after 10.5 h of light. The algae were then removed from the tank and disposed of in accordance with DAR procedure.

Acoustic recordings

Acoustic data were collected with a High Technology Inc. HTI-92-WB hydrophone, equipped with a low-power, high sensitivity preamplifier. Hydrophone and preamp system sensitivity was -144.8 dB re 1 VμPa-1 over the 2 Hz to 40 kHz band. The hydrophone was positioned to one side of the tank at 125 mm depth, 60 mm from the nearest tank wall. The preamplifier was powered by a variable DC power supply set at 10 V and drawing 2 mA. Data were acquired using a National Instruments® USB-6366 data acquisition module connected to a laptop computer running National Instruments® LabView® 2013 (SP1) and customized script. Time-stamped data were sampled at 100 KSs-1 at 16-bit depth in discrete 10 s intervals, pausing for a short time in between to save the data and verify that the sampling frequency was accurate. The recordings were processed in Mathworks Matlab® software. Before analysis, a high-pass, 5th order Butterworth filter (stop band: 0–0.5 kHz, pass band: 0.8–50 kHz was applied to remove DC offset, 60 Hz electrical noise or low frequency noise from passing vehicles and other sources.

Acoustic data processing

The peaks in amplitude associated with bubble separation ‘pings’ were identified from pre-filtered acoustic data if absolute peak amplitudes were greater than 5 σ from preamp self-noise levels. A 10 ms window around each of these peaks was isolated for further analysis. Information regarding the number of qualifying peaks and the temporal distribution of inter-peak periods within each 10 s file was retained. In order to calculate the sound exposure level (SEL) of bubble sound, the time duration of each air bubble “ping” was determined as the time between the 5th and 95th percentiles of the cumulative energy within the 10 ms window. SEL values were averaged over successive 10-min periods. Spectral analysis was performed on the time series between the 5th and 95th percentiles using a Kaiser-Bessel-windowed (β = 2.5π) 512-point FFT, overlapped 90 percent, meaning a bin width of 195.3 Hz and an overlapped spectral estimate made every 512 μs. The frequency bin containing the highest pressure spectral density level above 1 kHz, temporally matched to the amplitude peak in the corresponding time series, was determined to best represent the dominant frequency of each bubble separation ping. Regression of log-linear relationships and calculation of residuals was conducted using a simple logarithmic (base ten) model for the relationship between SEL and dissolved oxygen, and a natural logarithmic model for the relationship between the number and mean size of bubbles and dissolved oxygen. Residuals were calculated using linearized values.

Determination of bubble size

Estimates of bubble radii (R0) corresponding to each peak frequency were made using the Minnaert equation (1). For R0 to be a good approximation for the true radii of a freely oscillating bubble, the bubble radius must satisfy the following conditions: that the bubbles are “acoustically small”, the geometric mean of bubble radius and acoustic wavelength is large compared to the viscous boundary layer, the bubbles are “thermally large”, and the Laplace pressure is much less than the equilibrium pressure in the liquid (20).

Dissolved oxygen, pH conductivity, and temperature data

Dissolved oxygen (DO), pH, conductivity, and temperature were measured in 60 s intervals over the duration of the experiment. All non-acoustic sensors were part of a Manta+ 20 multi-probe system placed inside the experiment tank opposite the hydrophone. DO was assessed with an optical dissolved oxygen sensor provided by the Hamilton company (HDO) with a resolution of 0.01 mgL-1 or 0.1% saturation. pH was assessed via an electrolyte filled glass sensor next to a reference electrode. For conductivity, a four electrodes sensor was used while the temperature was assessed with a thermistor. Additional specifics are provided by the manufacturer [41]. While essential for the experiment, the oxygen sensor produced two types of acoustic noise that were additionally filtered from the acoustic data before analysis. The first type of noise was a startup transient that occurred upon system power-up every 60 s. A matched filter was used to exclude 0.2 s portions of recordings in which the transient was present. Secondly, a steady-state doubling of high frequency noise levels was emitted by the sensor for approximately 8 s after the startup transient. These recording periods were identified using a Hilbert-transform-based threshold detection algorithm and excluded from analysis.

Optical validation

Bubble radii were directly measured from high resolution photographs taken using a Canon SL1 single lens reflex camera equipped with a Canon 100 mm F2.8 L macro lens (160 mm effective) mounted on a tripod outside of the aquarium and aimed inward through a transparent window. A scale rule bar was inserted into the aquarium and placed across the field of view, orthogonal to the lens axis and within the depth of field. Only bubbles that were in focus within the depth of field of each image were analyzed. Photographs were taken every 150 s for one hour during a period when algae were observed to be actively producing bubbles. Radii measurements were obtained in post-processing using ENVI software, measured by pixel width and converted to μm using the scale bar in individual photographs.

Supporting information

S1 Fig. Chemical and acoustic responses of algae to changes in photosynthetically available radiation (PAR).

All subfigures are time-aligned. (A) Dissolved oxygen time series. (B) Acoustically determined bubble counts per minute. (C) Acoustically determined mean bubble radius per minute.

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

(TIFF)

S2 Fig. Log-linear regression plots of acoustic emissions from photosynthesizing algae.

(A) No. of bubbles per minute against dissolved oxygen levels (R2 = 0.76, exponential coefficients α = -2.54, β = 0.43). (B) Mean bubble radii per minute against dissolved oxygen levels (R2 = 0.47, exponential coefficients α = -3.85, β = 0.19). (C) 10-minute Sound Exposure Level against dissolved oxygen levels (R2 = 0.82, exponential coefficients α = 83.15, β = 1.39). The coefficients may be applied to an exponential regression of linear parameters x and y as follows: (2)

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

(TIFF)

Acknowledgments

We thank G. B. Deane, F. L. Rohwer, B. Bailey and R. F. Keeling for helpful comments and suggestions; K. Lubarsky, B. Neilson, D. Lager, K. Fuller, and K. Tucker for providing algae samples; J. Davidson and W. Au for providing laboratory facilities; S. Calhoun and B. T. Reyes for calibration of laboratory equipment; F. La Valle for the image in Fig 5; J. Tomlinson, P. D. Franck and S. M. Franck for experiment logistical support; G. J. Orris, S. Ackelson, R. D. Norris, M. J. Buckingham and A. J. Miller for funding support. This work was supported by NSF IGERT grant no. 0903551, U.S. Naval Research Laboratory section 219 program funds, National Research Council and American Society for Engineering Education postdoctoral fellowships and NASA grant no. NASA-14HYSP14-0003. Two anonymous reviewers provided comments that greatly improved the quality of the manuscript.

References

  1. 1. Haas AF, Guibert M, Foerschner A, Calhoun S, George E, Hatay M, et al. Can we measure beauty? Computational evaluation of coral reef aesthetics. PeerJ. 2015;3: pe1390. pmid:26587350
  2. 2. Scheffer M, Carpenter S, Foley JA, Folke C, Walker B. Catastrophic shifts in ecosystems. Nature. 2001;413 6856: 591–596.
  3. 3. Jackson JBC. Ecological extinction and evolution in the brave new ocean. Proc Nat Acad Sci. 2008;105: Supplement 11458–11465.
  4. 4. Smith JE, Shaw M, Edwards RA, Obura D, Pantos O, Sala E et al. Indirect effects of algae on coral: algae‐mediated, microbe‐induced coral mortality. Ecol Lett. 2006;9(7): 835–845. pmid:16796574
  5. 5. Pandolfi JM, Jackson JBC, Baron N, Bradbury RH, Guzman HM, Hughes TP et al. Are US coral reefs on the slippery slope to slime? Science. 2005;307(5716): 1725–1726. pmid:15774744
  6. 6. Vargas-Angel B,White D, Storlazzi C, Callender T, Maurin P. Baseline assessments for coral reef community structure and demographics on West Maui National Oceanic and Atmospheric Administration; 2017.
  7. 7. Sueur J, Pavoine S, Hamerlynck O, Duvail S. Rapid acoustic survey for biodiversity appraisal. PLoS ONE. 2008;3(12): e4065. pmid:19115006
  8. 8. Kennedy EV, Holderied MW, Mair JM, Guzman HM, Simpson SD. Spatial patterns in reef-generated noise relate to habitats and communities: evidence from a Panamanian case study. J Exp Mar Biol Ecol. 2010;395(1): 85–92.
  9. 9. Freeman SE, Rohwer FL, D’Spain GL, Friedlander AM, Gregg AK, Sandin SA et al. The origins of ambient biological sound from coral reef ecosystems in the Line Islands archipelago. J Acoust Soc Am. 2014;135: 1775–1788. pmid:25234977
  10. 10. Freeman LA, Freeman SE. Rapidly obtained ecological indicators in coral reef soundscapes. Mar Ecol Prog Ser.2016;501: 69–82.
  11. 11. Radford CA, Tindle CT, Montgomery JC, Jeffs AG. Modeling a reef as an extended sound source increases the predicted range at which reef noise may be heard by fish larvae. Mar Ecol Prog Ser. 2011;428; 167–174.
  12. 12. Akamatsu T, Okumura T, Novarini N Yan HY. Empirical refinements applicable to the recording of fish sounds in small tanks. J Acoust Soc Am. 2002;112(6); 3073–3082. pmid:12509030
  13. 13. Kraines S, Suzuki Y, Yamada K Komiyama H. Separating biological and physical changes in dissolved Oxygen concentration in a coral reef. Limnology and Oceanography. 1996;41(8): 1790–1799.
  14. 14. Minneart M. XVI On musical air-bubbles and the sounds of running water. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 1933;16(104): 235–248.
  15. 15. Leighton TG. The acoustic bubble ISBN: 0-12-44190-8; 1994.
  16. 16. Wilson CJ, Wilson PS Dunton KH. An acoustic investigation of seagrass photosynthesis. Marine Biology. 2012;159(10); 2311–2322.
  17. 17. Medwin H. In situ acoustic measurements of bubble populations in coastal ocean waters. J Geophys Res.1970;75(3): 599–611.
  18. 18. Tyree MT, Sperry JS. Vulnerability of xylem to cavitation and embolism. Ann Rev Plant Biol. 1989;40; 19–38.
  19. 19. Deane GB, Stokes MD. Scale dependence of bubble creation mechanisms in breaking waves. Nature. 2002;418(6900); 839–844. pmid:12192401
  20. 20. Wilson PS, Dunton KH. Laboratory investigation of the acoustic response of seagrass tissue in the frequency band 05–25 kHz. J Acoust Soc Am. 2009;125(4): 1951–1959. pmid:19354371
  21. 21. Ainslie MA, Leighton TG. Review of scattering and extinction cross-sections, damping factors, and resonance frequencies of a spherical gas bubble. J Acoust Soc Am. 2011;130(5): 3184–3208. pmid:22087992
  22. 22. Keeling RF. On the role of large bubbles in air-sea gas exchange and supersaturation in the ocean. J Mar Res. 1993;51(2): 237–271.
  23. 23. Freeman SE, Buckingham MJ, Freeman LA, Lammers MO, D'Spain GL. Cross-correlation, triangulation, and curved-wavefront focusing of coral reef sound using a bi-linear hydrophone array. The Journal of the Acoustical Society of America. 2015 Jan;137(1):30–41. pmid:25618036
  24. 24. Versluis M, Schmitz B, von der Heydt A, Lohse D. How snapping shrimp snap: through cavitating bubbles. Science. 2000 Sep 22;289(5487):2114–7. pmid:11000111
  25. 25. Au WW, Banks K. The acoustics of the snapping shrimp Synalpheus parneomeris in Kaneohe Bay. J Acoust Soc Am. 1998;103: 41–47.
  26. 26. Bohnenstiehl DR, Lillis A, Eggleston DB. The curious acoustic behavior of estuarine snapping shrimp: Temporal patterns of snapping shrimp sound in sub-tidal oyster reef habitat. PLoS ONE 11(1): e0143691. pmid:26761645
  27. 27. Lillis A, Eggleston DB, Bohnenstiehl DR. Estuarine soundscapes: distinct acoustic characteristics of oyster reefs compared to soft-bottom habitats. Marine Ecology Progress Series. 2014 May 28;505:1–7.
  28. 28. Staaterman E, Rice AN, Mann DA, Paris CB. Soundscapes from a Tropical Eastern Pacific reef and a Caribbean Sea reef. Coral Reefs. 2013 Jun 1;32(2):553–7.
  29. 29. Staaterman E, Paris CB, DeFerrari HA, Mann DA, Rice AN, D’Alessandro EK. Celestial patterns in marine soundscapes. Marine Ecology Progress Series. 2014 Aug 4;508:17–32.
  30. 30. Radford C, Jeffs A, Tindle C, Montgomery JC. Resonating sea urchin skeletons create coastal choruses. Mar Ecol Prog Ser. 2008;362: 37–43.
  31. 31. McCauley RD, Day RD, Swadling KM, Fitzgibbon QP, Watson RA, Semens JM. Widely used marine seismic survey air gun operations negatively impact zooplankton. Nature Ecology and Evolution. 2017;1(7); s41559–017.
  32. 32. Luczkovich JJ, Mann DA, Rountree RA. Passive acoustics as a tool in fisheries science. Transactions of the American Fisheries Society. 2008;137: 533–541.
  33. 33. Giorli G, Au WW, Neuheimer A. Differences in foraging activity of deep sea diving odontocetes in the Ligurian Sea as determined by passive acoustic recorders. Deep Sea Research Part I: Oceanographic Research Papers. 2016;107: 1–8.
  34. 34. Chapman NR, Price A. Low frequency deep ocean ambient noise trend in the Northeast Pacific Ocean. J Acoust Soc Am. 2011;129(5): 161–165.
  35. 35. Dahl PH, Miller JH, Cato DH, Andrew RK. Underwater ambient noise. Acoustics Today. 2007;3(1): 23–33.
  36. 36. Dunlop RA, Cato DH, Noad MJ. Your attention please: increasing ambient noise levels elicits a change in communication behaviour in humpback whales (Megaptera novaeangliae). Proceedings of the Royal Society of London B: Biological Sciences. 2010; prspb20092319.
  37. 37. Giorli G, Au WW, Ou H, Jarvis S, Morrissey R, Moretti D. Acoustic detection of biosonar activity of deep diving odontocetes at Josephine Seamount High Seas Marine Protected Area. J Acoust Soc Am. 2015;137(5): 2495–2501. pmid:25994682
  38. 38. Boyd JW, Varley J. The uses of passive measurement of acoustic emissions from chemical engineering processes. Chem Eng Sci. 2001;56(5); 1749–1767.
  39. 39. Preskitt LV, Vroom PS, Smith CM. A rapid ecological assessment (REA) quantitative survey method for benthic algae using photoquadrats with scuba. Pacific Science. 2004;58: 201–209.
  40. 40. Yentch CS, Yentsch CM, Cullen JJ, Lapointe B, Phinney DA, Yentsch SW. Sunlight and water transparency: cornerstones in coral research. J Exp Mar Biol Ecol. 2002;268(2): 171–183.
  41. 41. Sensors for sondes and monitoring. 2017;11:20. Available from: https://wwwwaterprobescom/sensors-for-sondes-and-monitoring.