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

Quantification of chemical and mechanical bioerosion rates of six Caribbean excavating sponge species found on the coral reefs of Curaçao

  • Didier M. de Bakker ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    didier.debakker@wur.nl, didierdebakker@gmail.com

    Affiliations Wageningen Marine Research, Den Helder, Netherlands, Marine Microbiology and Biogeochemistry, NIOZ Royal Netherlands Institute for Sea Research and Utrecht University, Texel, Netherlands

  • Alice E. Webb,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Ocean Systems, NIOZ Royal Netherlands Institute for Sea Research and Utrecht University, Texel, Netherlands

  • Lisanne A. van den Bogaart,

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Marine Biology Research Group, Ghent University, Ghent, Belgium

  • Steven M. A. C. van Heuven,

    Roles Conceptualization, Resources, Writing – review & editing

    Affiliation Ocean Systems, NIOZ Royal Netherlands Institute for Sea Research and Utrecht University, Texel, Netherlands

  • Erik H. Meesters,

    Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing

    Affiliation Wageningen Marine Research, Den Helder, Netherlands

  • Fleur C. van Duyl

    Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Validation, Writing – review & editing

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

Abstract

Excavating sponges are among the most important macro-eroders of carbonate substrates in marine systems. Their capacity to remove substantial amounts of limestone makes these animals significant players that can unbalance the reef carbonate budget of tropical coral reefs. Nevertheless, excavating sponges are currently rarely incorporated in standardized surveys and experimental work is often restricted to a few species. Here were provide chemical and mechanical bioerosion rates for the six excavating sponge species most commonly found on the shallow reef of Curaçao (southern Caribbean): Cliona caribbaea, C. aprica, C. delitrix, C. amplicavata, Siphonodictyon brevitubulatum and Suberea flavolivescens. Chemical, mechanical and total bioerosion rates were estimated based on various experimental approaches applied to sponge infested limestone cores. Conventional standing incubation techniques were shown to strongly influence the chemical dissolution signal. Final rates, based on the change in alkalinity of the incubation water, declined significantly as a function of incubation time. This effect was mitigated by the use of a flow-through incubation system. Additionally, we found that mechanically removed carbonate fragments collected in the flow-through chamber (1 h) as well as a long-term collection method (1 wk) generally yielded comparable estimates for the capacity of these sponges to mechanically remove substratum. Observed interspecific variation could evidently be linked to the adopted boring strategy (i.e. gallery-forming, cavity-forming or network-working) and presence or absence of symbiotic zooxanthellae. Notably, a clear diurnal pattern was found only in species that harbour a dense photosymbiotic community. In these species chemical erosion was substantially higher during the day. Overall, the sum of individually acquired chemical and mechanical erosion using flow-through incubations was comparable to rates obtained gravimetrically. Such consistency is a first in this field of research. These findings support the much needed confirmation that, depending on the scientific demand, the different approaches presented here can be implemented concurrently as standardized methods.

Introduction

The existence of tropical coral reefs relies on an ongoing biogenic precipitation of calcium carbonate (CaCO3) by calcifying organisms (most importantly scleractinian corals) that is exceeding erosional forces [17]. Live corals, however, have been diminishing worldwide at unprecedented rates following phenomena such as reduced water quality, climate change and declined herbivorous pressure (e.g. [8, 914]). Reefs in the Caribbean and Gulf of Mexico, in particular, have been heavily impacted in the past decades. Here coral cover was reduced, on average, to a mere 16% [1517]. Alarmingly, the observed decline can largely be accredited to the mortality of formerly dominant framework-building species (i.e. Acropora spp. and Orbicella spp.) which are especially affected by direct or indirect anthropogenic stressors [1820]. Whilst the main reef-builders appear largely unable to cope with the current level of global and local disturbance [21], the changing reef environment seems to favour, among others, bioeroding organisms) (For an overview see [22]: Table 3, Supplementary data S4).

Within the bioeroding community excavating sponges are frequently the dominant macroboring organisms [2325]. On Caribbean reefs sponges can be responsible for as much as 90% of the total macrobioerosion [26]. Murphy et al. (2016) [27] describe two main strategies by which these sponges erode their substratum: (1) gallery-forming, where the sponge progressively works its way down, and (2) cavity-forming, where the sponge penetrates the limestone and forms chambers inside the substrate with often only the inhalant and exhalant fistules being visible on the surface. Many gallery-forming sponges harbour high densities of dinoflagellate zooxanthellae (i.e. Symbiodinium spp.) and as such their boring strategy allows for optimal light exposure [28]. Excavating sponges erode hard substratum through a combination of chemical dissolution and mechanical removal of CaCO3 fragments (chips, distinctly recognizable by their scalloped surface) [2931]. Previous efforts studying the contribution of the individual components suggest that the mechanical fraction generally accounts for the majority (up to 98%) of total erosion (e.g. [29, 32, 3335]), although an opposite pattern has been described by Zundelevich et al. (2007) [36] for Pione vastifica (β-stage). The product of mechanical erosion by sponges contributes significantly to the total pool of fine silt-sized sediment on coral reefs ([37], and references therein), which is essential to reef cementation and stabilization [38].

There is ample scientific support that excavating sponges respond positively in both abundance and boring activity to factors such as deteriorating water quality (e.g. eutrophication), ocean acidification, elevated sea surface temperature and increased availability of substratum following coral mortality ([22] and references therein). Nevertheless, they are rarely incorporated in standardized reef surveys and experimental work is restricted to a limited number of species. This is in part due the endolithic life-style of the sponge. The few surveys that do include excavating sponges, such as the ReefBudget method [39], have to rely on scarce and often inconsistent data with respect to rates of bioerosion. Currently, techniques to estimate these rates are inconsistent and the procedures involved are tedious and prone to methodological bias.

The available published data on total estimated sponge bioerosion displays tremendous variation with rates ranging from 0.3 kg m-2 y-1 to 29.5 kg m-2 y-1 ([40]: S3 Table]). Density of the invaded substratum has a considerable impact on the attained rates (e.g. [41, 42]), but also factors such as studied species, boring strategy, developmental stage (α, β, γ or δ-stage) and applied methodology can cause significant variability ([40] S3 Table].) In early work erosion rates were generally determined based on weight loss of the host substratum over a designated period [4245], occasionally in combination with chip collection experiments to quantify mechanical erosion. Buoyant weight has been commonly implemented to estimate coral growth rates, but has also proven a good proxy to quantify rates of sponge bioerosion (e.g. [42, 46]). A major advantage of the buoyant weight technique is that the studied material can remain submerged during the weighing process. This approach on itself, however, is not selective to means of erosion (i.e. chemical or mechanical), which is relevant, for instance, when studying the effect of ocean acidification on the chemical fraction of sponge bioerosion. Recent efforts have therefore implemented incubation techniques to study both components separately. By incubating sponge infested substratum it is possible to determine chemical bioerosion based on the increase in alkalinity of the incubation-water over time [36]. This technique has now been widely adopted to quantify effects of environmental change and future climate scenarios on rates of bioerosion (e.g. [34, 35, 47, 48, 49]). A major disadvantage of this method, however, is that sponges are kept in a limited body of water for a relatively long time, which could strongly affect the metabolism of the sponge following food and oxygen depletion and the accumulation of waste products [50]. Moreover, the lack of understanding of the processes underlying mechanical chip removal adds great uncertainty to previously quantified mechanical bioerosion rates. Indeed, periodic expelling of chips and/or contamination of filters by other material may induce error in the quantification of mechanical bioerosion (e.g. [34, 35]). Evidently, the considerable variation in published sponge bioerosion rates stresses the need for consistency in applied procedures.

In the present study, the chemical and mechanical bioerosion rate of six boring sponge species that dominate bioeroding community on the southern Caribbean reefs of Curaçao was determined. Chemical rates were estimated using flow-through incubation techniques, which were compared to previously used standing incubations. For mechanical rates a distinction was made between actual chips and other CaCO3 material removed by the eroding activities of the sponge over a designated time span. These approaches were then compared to conventional buoyant weighing procedures. We aim to optimize techniques to quantify bioerosion and to determine variation in chemical and mechanical bioerosion of sponge species with different functional traits (e.g. presence/absence of photosymbionts, boring strategies and developmental stages) and in light of a changing reef environment.

Materials and methods

Ethics statement

All research was conducted under research permit (#2012/48584) issued by the Curaçaoan Ministry of Health, Environment and Nature (GMN) to the CARMABI Foundation.

Study area and sample collection

Sponges were collected on reefs at the leeward side of the island of Curaçao (southern Caribbean) in February 2017. Reefs here are characterized by a reef flat and a drop-off that gradually slopes down from approximately 10 m depth [51]. Samples were collected at Snake Bay (12°8’N, 68°59’W), Piscadera Bay (12°7’N, 68°58’ W) and Directors Bay (12°3’N, 68°51’W) between 5 and 15 m depth. Sponge infested cores of the most dominant species were removed from old dead coral substratum using a pneumatic drill with hole-saw (inner diameter 45 mm). Collected cores were cleaned and non-sponge organisms removed. The collected material was kept in large flow-through aquaria in the research facilities of the Carmabi Research Station. Cores were left to acclimatize and recover for at least seven days. Full tissue regeneration was observed within this period (see also [52]).

Species description

A total of six excavating sponge species were collected including two gallery-forming species (Fig 1): Cliona caribbaea (Carter 1882) and C. aprica [28], three cavity-forming species: C. delitrix [28], C. amplicavata [53] and Siphonodictyon brevitubulatum [28] and Suberea flavolivescens, a species that formed a dense network through the carbonate substratum (from here onwards referred to as network-forming). This is, to our knowledge, the first scientific description of C. amplicavata and S. flavolivescens on the reefs of Curaçao. Both gallery-forming species belong to the Cliona viridis species complex [54, 55] and harbour high densities of dinoflagellate zooxanthellae [28]. C. delitrix has also been reported to hold zooxanthellae, but in considerably lower densities [56]. The other three species do not appear to harbour symbiotic zooxanthellae [28, 53]. α-Stage specimens were collected for C. aprica, C. delitrix, C. amplicavata, and S. brevitubulatum and β-stage specimens for C. caribbaea and S. flavolivescens. The α-stage of C. delitrix (often confused with C. laticavicola) [57] was chosen since the β-stage was only recently observed for the first time on Curaçaoan reefs [58]. To allow for inter-specific comparison, all specimens, with the exception of C. amplicavata, were collected from infested dead Orbicella spp. coral, which is the most common substratum on many Caribbean reefs [59, 60] and of average skeletal density (~1.8 g cm-3) [61]. C. amplicavata was only found in Acropora spp. rubble. Same sized uncolonized cores were collected from old Orbicella spp. and Acropora spp. skeleton to serve as control substratum throughout all experiments allowing adjustment for possible dissolution through micro-bioerosion, abrasion caused by handling and accretion due to the potential presence of calcifying organisms (e.g. [35]). In order to prevent decomposition of S. brevitubulatum, the majority of the soft internal tissue was removed and cores were left to recover for several additional days before being used experimentally. Species identification was confirmed by spicule morphology analysis.

thumbnail
Fig 1. In situ close-up of the six studied excavating sponge species.

A. Cliona amplicavata (α-stage) B. Cliona aprica (α-stage) C. Cliona caribbaea (β-stage) D. Cliona delitrix (α-stage) E. Siphonodictyon brevitubulatum (α-stage), please note that the yellow fistules of S. brevitubulatum are surrounded by a different sponge species (pink coloured) F. Suberea flavolivescens (α-stage). Black arrows in A, E and F point out oscula and ostia of the studied sponges.

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

Experimental set-ups

Sponge bioerosion was measured in different experimental settings in large flow-through (~5 L min-1) holding tanks. We used 500 mL incubation chambers consisting of a transparent body with detachable top and bottom. The top contained an electrically driven magnetic stirrer to ensure thorough homogenisation of the incubation water. The incubation chambers were used in a standing mode (ST) without water flowing through and in a continuous flow-through mode (FT). In addition and for comparison total bioerosion was also determined gravimetrically by applying the buoyant weight (BW) technique [62, 63] to hanging cores. Debris from these cores was collected to determine the mechanical bioeroding capacity of the sponges.

With the conventional standing incubation the effect of incubation time on rates of chemical erosion was tested for two sponge species. C. caribbaea (symbiont-rich) and C. delitrix (symbiont-poor) infested cores were incubated for 1, 2, 4 and 6 h in fully enclosed incubation chambers. Two cores were placed in each chamber which was sealed air-tight to avoid gaseous exchange with the external environment. Chambers were largely submerged to prevent temperature anomalies. Temperature (°C) and oxygen (%) were monitored continuously over the course of the incubation using a PreSens O2 sensor (Fibox 4, PSt3). Incubations at night were carried out for a maximum of 3 hours to avoid oxygen levels declining too much (> 80% loss of O2 saturated water). Samples for alkalinity (~250 mL) and nutrient (PO43-, NH4+, NO3-, NO2-) analyses were collected at the start and the end of each run.

For the flow-through set-up the incubation chamber was modified to allow a constant through flow of water (Fig 2). Each chamber contained two sponge-infested or control cores. By means of a peristaltic pump, with adjustable flow-speed, fresh sea water was transported from a reservoir into the chamber where the magnetic stirrer ensured homogenized mixing. The excess water overflowed through flexible Teflon® tubing into a collection vessel. A filter (0.7 μm pore size) was placed on the out-flow tube to prevent loss of particulate material from the chamber. A flow-speed of 1.5 L h-1 allowed for continuous refreshment of the incubation water approximately every 20 minutes while still obtaining a reliable signal for water chemistry analyses. This flow-through incubation method ensured a constant supply of food and prevented excessive build-up of waste products, therewith simulating more adequately in situ reef conditions. Water samples for alkalinity (AT, μmol kg-1) and nutrient analyses were collected at the inlet and outlet of the incubation chamber. The system was run for 30 minutes before sample collection started, to allow for an initial full refreshment of the incubation water. This was done to prevent catching a signal that might be affected by the initial incubation water the cores were placed in. Additionally, a possible effect of sponge handling and subsequent acclimatization was herewith largely excluded. Final dissolution rates were calculated from the observed difference in AT which was corrected for the retention time of the water in the incubation chamber. Cores were carefully removed from the FT chamber immediately after water sample collection. Subsequently, all particulate material produced during the entire duration of the incubation (~1 h) was collected as a measure for mechanically removed substratum (mechanical erosion). Chemical and mechanical erosion of sponge infested and control cores was determined at day and night in the FT incubation because erosion rates of sponges, particularly symbiont-bearing species, were found to be affected by light availability (diurnal variation) [34, 64, 65]. Flow through incubations were run for all six sponge species. For two of them, C. caribbaea and C. delitrix, results (chemical dissolution rates) of the FT were compared to results of the ST incubation. In order to test the consistency in time of the FT methodology, additional water samples were collected for C. caribbaea 1, 2, 3 and 4 h after the start of the incubation while the FT system kept running.

thumbnail
Fig 2. Flow-through set-up.

A. Overview of the flow-through incubation set-up including the peristaltic pump with adjustable flow-speed continuously pumping fresh sea water through the incubation chambers. The excess water subsequently overflows into collection vessels for nutrient and total alkalinity analyses. Chambers are partly submerged in a holding tank to prevent temperature fluctuation. B. Close-up of a flow-through incubation chamber (500 mL) holding two sponge infested limestone cores.

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

Mechanical bioerosion by the six species was also determined over a longer period (7 consecutive days) as part of the buoyant weight approach. Three infested cores per species and blank cores were positioned upside down directly above glass funnels (Ø 10 cm) in a large holding tank. Particulate material caught in each funnel (placed in a collection vial) as well as in the FT incubation chamber was collected on pre-combusted (550°C, 3 h), pre-weighed, 0.7 μm Whatman® glass microfiber filters. Filters were subsequently re-combusted (550°C, 3 h) to remove any organic material and re-weighed resulting in an estimate of the particulate material produced over the course of each experimental approach. As an additional analysis, material collected on the filters of the one week experiment was homogenized and subsampled. Scanning Electron Microscopy (SEM) of the subsample was used to estimate the actual contribution of sponge chips, non-chip CaCO3 fragments and other material (e.g. spicules, foreign materials). FT filters contained too little material for adequate SEM analysis. Results of both experiments were corrected for the weight of fragments collected from uncolonized control cores to ensure that only the fraction removed by the sponge would be included in estimates of total mechanical erosion.

To determine total bioerosion based on the loss in buoyant mass, the hanging cores were first weighed at the beginning of the 7-day collection period following the BW technique as described by Fang et al. (2013b) [35]. Cores were re-weighed after hanging for three full weeks. All measurements were conducted by means of an electronic weighing apparatus (0.1 mg accuracy). Organic sponge tissue was assumed to equal the density of the ambient sea water so to eliminate the effect of sponge growth. Buoyant weight of the cores was calibrated for fluctuations in seawater density. The change in weight of each core over the three week period was corrected by the change in buoyant mass of the uncolonized control cores and subsequently converted into rates of total sponge bioerosion.

Chemical erosion

The mass of CaCO3 dissolved by each sponge (ΔM) was calculated from the change in total alkalinity over the course of the incubation following Eq (1) [32, 36]. Fluctuations in PO43-, NH4+, NO3- and NO2- can affect ΔAT and are therefore included in the equation as well [47, 66]. (1) where ΔAT is the change in total alkalinity. Multiplication factor ‘0.5’ reflects the 2:1 relationship between the increase in AT during CaCO3 dissolution. VSW is the volume (L) of seawater in the incubation chamber and ρSW is seawater density (1.023 kg L-1). The multiplier ‘100’ incorporates the molecular mass of CaCO3.

AT was measured using an Automated Spectrophotometric Alkalinity System (ASAS) [67] following the method outlined by Breland and Byrne (1993) [68] and Yao and Byrne (1998) [69]. This optical titration procedure has a remarkable high precision (± 0.7 μmol kg-1) making it possible to detect minor alkalinity fluctuations. AT was measured immediately after sample collection. To compensate for a drift in AT over the course of the experimental period, certified reference material (CRM; supplied by Dr. A. Dickson, Scripps Institution of Oceanography) was analyzed every 20 samples. Filtered (acrodisc: 0.2 μm pore size) nutrient samples were stored at -20°C until transport to the Royal Netherlands Institute for Sea Research (NIOZ) where they were analysed on a QuAAtro continuous flow analyzer (SEAL Analytical, GmbH, Norderstedt, Germany) following GO-SHIP protocol [70]. Final chemical erosion rates were expressed in mg CaCO3 cm-2 d-1 or convertible units and nutrient dynamics in μmol L-1 h-1.

Determining sponge surface area and biomass

At the end of each experiment all cores were photographed to estimate the surface area covered by external sponge tissue. For α-stage (only single inhalant and exhalant papillae are visible on the surface) sponges, the area in between the outermost papillae was taken as a measure for infested surface area. For β-stage specimens (papillae are fused into an encrusting sheet) the total sponge covered surface on both the top and side of each core was used. Additionally, the cores used in the BW experiment were dissolved in acid (2 M HCl) after the final weighing to isolate the sponge tissue for wet and dry-weight (24 h at 60°C) measurements and to verify the absence of other macroborers (e.g. worms and molluscs).

Statistics

All statistical testing was conducted in the R programming environment R 3.3.2 [71] using the packages “stats” and “gam”. Assumptions of normality and homogeny were checked visually and all collected data were fourth-root transformed to stabilize variances. Generalized additive modelling was implemented to predict the fluctuations in rates of chemical erosion as a function of time in ST incubations. Two-way ANOVA (species and diurnal cycle) testing complemented by post-hoc Tukey-HSD paired comparisons were used to examine variation in patterns of excavation among and within the six studied sponge species.

Results

Chemical erosion

Comparison of chemical dissolution rates acquired through ST and FT incubation showed that duration of the ST incubation had a distinct effect on sponge chemical bioerosion. Dissolution rates by C. caribbaea (F(3,8) = 33.9, p < 0.001) and C. delitrix (F(3,8) = 27.3, p < 0.001) declined significantly as a function of the time cores were retained in the confined water body of the incubation chamber (Fig 3). After 2 hours (T2), rates had already decreased by 31% (C. caribbaea) and 70% (C. delitrix) compared to the modelled rate at T0, which is the hypothetical point at which the sponge should not yet experience any negative effect of being incubated. Comparison of dissolution rates between the ST and the FT during day-time (Fig 3) shows that the chemical dissolution extrapolated to T0 of the ST incubation, C. caribbaea: 0.014 mg cm-2 h-1 (95% CI: 0.011–0.017) and C. delitrix: 0.021 mg cm-2 h-1 (95% CI: 0.012–0.032), agrees with dissolution values obtained in FT incubation (Table 1, daily rates). A negative effect of incubation time was not observed for the two species in the flow-through methodology, at least not up to 4 h after the start of the incubation (Fig 3). This supports the long-term consistency of the FT methodology introduced here. Consequently it was decided to focus on measurements with the FT set-up.

thumbnail
Fig 3. Modelled effect of incubation time on sponge chemical bioerosion rates.

Black regression lines represent the modelled (GAM) mean dissolution in mg cm-2 h-1 by Cliona caribbaea (left panel) and Cliona delitrix (right panel) in standing incubations. 95% confidence limits are indicated by the grey band. Solid dots show the calculated mean dissolution rates at T1, T2, T4 and T6 (ST). Open dots represent rates derived from running flow-through incubations at T1, T2, T3 and T4, dotted lines here indicate 95% confidence limits. All rates are obtained using alkalinity anomaly techniques.

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

thumbnail
Table 1. Summary of all calculated sponge bioerosion rates.

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

Chemical dissolution rates acquired through FT incubation revealed significant (F(38) = 4.11, p < 0.01) variation among sponge species (Fig 4 and Table 1). Interspecific variation was predominantly associated with morphology resulting in significantly slower erosion rates for both gallery-forming species: C. caribbaea and C. aprica. These high density zooxanthellate sponges also showed clear diurnal disparity with significantly higher (p < 0.01 for both species) dissolution rates during the day (Fig 4). In C. amplicavata, C. delitrix, S. brevitubulatum and S. flavolivescens night erosion appeared to be higher than day erosion, although this was not significant for any of these species. Notably, the inflow water was slightly more basic during day-time (pH = 8.08 ± 0.001) compared to night-time (pH = 8.02 ± 0.007).

thumbnail
Fig 4. Chemical bioerosion.

Hourly rates of chemical dissolution (mg CaCO3 cm-2 d-1) are given for the six studied sponge species (top). Both day (light grey) and night (dark grey) erosion rates are included. Rates are derived from alkalinity change in flow-through incubation. Error bars indicate the 95% confidence limits.

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

Nutrient dynamics

Considerable nutrient accumulation (especially NH4+ and NO3-) occurred at day-time in ST (1 h) incubations holding C. caribbaea (Table 2). Accumulation of nutrients after 1 hour was largely prevented in the FT system (Table 2). Observed fluctuations in phosphate concentrations were largely negligible, there was essentially no indication for net uptake. In FT incubations a minor elevation in levels of ammonium and nitrate was observed during most incubations at night. During the day, these nutrients appear to be consumed in the majority of the incubations, including those containing non-infested cores. Nitrogen can be taken up by the sponge itself or by its zooxanthellae (i.e. C. caribbea, C. aprica and to a lesser extend C. delitrix). The observed decline in control incubations (containing uncolonized cores) indicates that there is likely also an effect on nutrient levels caused by micro-organisms in the incubation water. When considering the nutrient trajectories in the control water many of the observed dynamics in chambers containing sponges are largely dampened (Table 2).

thumbnail
Table 2. Summary of nutrient dynamics in the incubation chambers.

https://doi.org/10.1371/journal.pone.0197824.t002

Mechanical erosion

The mechanical fraction of sponge bioerosion was based on the weight of the particulate material collected in the incubation chamber after 1 h FT incubation and in the funnels for the long-term (7 day) collection experiment, yielding two estimates per species (Fig 5 and Table 1). Although both estimates are in a similar range, rates derived from FT incubation are consistently higher. SEM analyses of 7-day filters did revealed that, depending on the species, 17% to 45% of the collected material (long-term) could be attributed to actual sponge chips (Fig 6). The remaining material consisted mostly of non-chip CaCO3 fragments with some occurrences of foreign planktonic material or spicules. The fraction of non-chip material that was dislodged by sponge could not be traced back directly to boring activity and thus represents an estimate (i.e. total non-chip CaCO3 fragments minus the fragment production by non-infested cores). CaCO3 fragment removal differed significantly among species (F(30) = 4.43, p = 0.003), with most material removed by C. delitrix and C. amplicavata. Similar to the chemical fraction, there appears to be a strong diurnal pattern in mechanical substrate removal in both symbiont-rich species as well as the symbiont-poor C. delitrix (Table 1). ANOVA testing of the FT results showed that mechanical erosion at night was significantly lower only for C. caribbaea (F(1,5), = 59.6, p = 0.016). It should be noted, however, that only one sample was available for C. aprica at day which strongly hampered proper statistical testing. Mechanical erosion by the other four species did not differ significantly at day or night (p > 0.05).

thumbnail
Fig 5. Total sponge bioerosion.

Final rates of total erosion by the six sponge species (top) acquired through flow-through incubation (FT) and buoyant weight loss over a three week period (BW). The different bioerosive components (FT) are based on a 12 h day/night diurnal cycle. Mechanical bioerosion rates based on long-term collection are provided separately (LT). Here, a distinction is made between actual sponge chips and other CaCO3 fragments that have been dislodged through sponge bioeroding activity. Rates are given in mg CaCO3 cm-2 d-1 and kg CaCO3 m-2 y-1. Error bars indicate the 95% confidence intervals for each individual element. Species and boring strategies are provided at the top of each panel, CF: cavity-forming, GF: gallery-forming; NF: network-forming.

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

thumbnail
Fig 6. Scanning Electron Microscopy (SEM) image of material produced by C. caribbaea.

The provided image shows only a fraction of the total material produced over the 7-day collection experiment. Black arrows: some characteristic sponge chips, white arrows: other CaCO3 fragments. Scale is provided in lower-right corner.

https://doi.org/10.1371/journal.pone.0197824.g006

Total bioerosion

Total bioerosion (i.e. chemical and mechanical) quantified from FT incubations was consistently higher than BW acquired rates (Table 1 and Fig 5). For C. amplicavata, S. brevitubulatum and S. flavolivescens both rates are still within a similar range. However, the discrepancy between methods was particularly pronounced in both symbiont-rich species (C. caribbaea and C. aprica) and S. flavolivescens. Nevertheless, the observed patterns of inter-specific variation were highly comparable among both approaches, with highest erosion rates found for the three cavity-forming species. Lowest rates are consistently observed for network-forming S. flavolivescens. There appears to be considerable variation in BW rates among cores of the same species (suggested by the large confidence intervals), in particular among the cavity-formers (Fig 5). BW rates converted to mg CaCO3 g DWsponge-1 d-1 based on the dry weight of the sponge tissue in each core yielded considerable similarity in patterns of inter-specific variation to rates based on surface area (Table 1 and S1 Fig). The cavity-forming C. amplicavata, however, appear to be an exception to this apparent surface area–biomass relation. FT estimates indicate that for five out of the six species the mechanical fraction accounted for the vast majority of total bioerosion (86%–96%). Interestingly, for the network-forming S. flavolivescens chemical dissolution (46%) and mechanical fragment removal (54%) contributed almost equally.

Discussion

We quantified the bioerosive capacity of six excavating sponge species commonly found on the reefs of Curaçao. This is the first time individual chemical and mechanical rates and total bioerosion rates are provided simultaneously for the majority of these species. Overall, Cliona delitrix and Cliona amplicavata (both cavity-forming species) were found to have the highest boring rates. C. delitrix is indeed known to be among the most destructive bioeroding sponges on many Caribbean reefs [72]. It should be pointed out that C. amplicavata was found solely in Acropora spp. rubble which is of higher skeletal density than Orbicella spp. [61]. Substratum density can strongly affect bioerosion, where denser substratum generally leads to higher erosion rates (e.g. [42, 43]). Erosion by C. caribbaea and C. aprica, both adopting a gallery-forming strategy, was found to be considerably lower compared to cavity-formers. Nonetheless, these species, together with C. delitrix, are the dominant excavating sponges on the shallow reefs of Curaçao and many sites in wider Caribbean region [39, 57, 58] and thus contribute significantly to total macrobioerosion. S. flavolivescens, previously undescribed to occur on Curaçaoan reefs, appears to be the weakest eroder of all sponge species studied here. This network-forming species adapts a different boring strategy (i.e. network-forming) from those outlined by Murphy et al. (2016) [27].

Methodological considerations

Comparison of the chemical bioerosion in the flow-through (FT) and standing (ST) approach reveals that an incubation time of 1 h or longer in the ST incubation chambers was too long to estimate the purported “real” bioerosion rate. Bioerosion steeply dropped in time, suggesting that food depletion and or waste accumulation inhibited chemical bioerosion. Neumann (1966) [45] suggests that sponges can bore more actively in a high water flow environment, possibly due to a larger supply of fresh sea water. Being known as highly efficient filter-feeding organisms (e.g. [73]) it is not surprising that these sponge have the capacity to rapidly consume the available food sources, both particulate and dissolved [74] in any limited body of water. Consequential stress could significantly affect the metabolism of the sponge, possibly redirecting energy away from secondary processes such as boring and focussing more on vital functions like maintaining a high pumping velocity. Extrapolation of rates determined in the ST to T0 revealed chemical dissolution rates comparable to rates found in the FT incubators. This implies that whenever the ST set-up is preferable, the optimal incubation times should be determined in advance in the devices used and individually for each sponge species. To date, standing incubation techniques generally resulted in lower rates for chemical erosion, but have also predominantly been used to test the effect of environmental change on the boring efficiency through sea water manipulation ([40] and references therein). Although in most of these studies the authors looked at relative variation in rates as a response to treatment, the acquired rates should be interpreted with some degree of caution. The processes acting inside these enclosed chambers are largely unknown and unpredictable. More desirable, therefore, would be the use of FT incubation, since here the effect of incubation time on measured variables appears to be negligible (Fig 3 and Table 2).

CaCO3 fragment production (including sponge chips) in the FT incubation chambers was generally higher than the fragment production by hanging cores in a large aquarium. Webb et al. (2017) [34] propose that elevated levels of Ca2+ in incubation water (ST) following chemical dissolution of CaCO3 may cause contraction of the conductive pathway of the sponge. Such a contraction could result in an initial pulse ejection of the chips present within the sponge tissue, resulting in excessively high production over the relatively short incubation. However, in our FT incubations accumulation of Ca2+ will have been limited because the incubation water was constantly refreshed. Since the discrepancy between methods was most pronounced in both symbiont-rich species (C. caribbaea and C. aprica) and S. flavolivescens, it is more likely that the reversed orientation of the hanging cores may have reduced the boring capacity of these sponges following limited light availability. In addition, this methodology (hanging cores) might be more prone to loss of a fraction of the produced silt-sized fragments due to water movement in the holding tank [35]. Overall, there is still hardly any information available on the mechanisms underlying mechanical erosion or the expelling of dislodged material by excavating sponges. Both approaches (short and long term collection), however, can provide better insight in the mechanisms underlying mechanical erosion and the expelling of dislodged material by excavating sponges.

For most species, final FT rates are comparable to rates acquired gravimetrically based on the loss in buoyant weight (BW) of the invaded cores (Fig 5). Strongest variation between methodologies was found in C. caribbaea and C. aprica and S. flavolivescens. Similar to the long-term fragment collection, the relatively low BW estimates are likely the result of the reversed orientation of the cores (see previous paragraph). This might also explain why for some hanging species mechanical erosion exceeded total bioerosion. CaCO3 fragment production was only measured in the first week. Possibly the boring capacity of the hanging sponges had further decreased in the two following weeks. We therefore argue that BW estimates in this study may reflect a slight underestimation of sponge bioeroding capacity and FT obtained rates should be considered more reliable. Up to now, however, conflicting studies generally found up to five times higher rates obtained through buoyant weighing compared to standing incubation derived estimates [34, 35]. The observed discrepancy has caused scepticism regarding the use of BW technique. Our results, however, support the use of the buoyant weight methodology, as described in previous studies, as a reliable approach to estimate rates of bioerosion and suggest that the previous implementation of standing incubation may often have resulted in underestimated rates. Nevertheless, to examine the effect that changing environmental conditions might have on bioerosion it is essential to consider both fractions independently. Particularly in view of the effect factors such as elevated levels of CO2 in ocean water have on the chemical dissolution capacity of the various sponge species. Preferably this is done in an incubation set-up where chemical erosion can be determined individually as well. Consequently, the methods presented here, with some minor improvements where necessary (read BW, long-term fragment collection), can and should be used concurrently depending on the underlying scientific aim.

Bioerosion by different species

FT incubation revealed that the chemical bioerosive capacity of C. caribbaea and C. aprica was notably higher during the day (Fig 4). This is likely related to the presence of photosynthetic symbionts in these viridis-complex species. The photosynthetic capacity of intracellular symbiotic zooxanthellae was pinpointed as a possible primary source of carbon for these sponge species [65, 75, 76]. A consequential diurnal variation in rates of bioerosion following light availability has previously been described for C. caribbaea in St. Eustatius (Lesser Antilles) [34], Cliona varians in the Florida Keys [65] and Cliona orientalis in Australia [64, 77], all of which are species in the C. viridis species complex. Symbiodinium spp. in C. delitrix, present in much lower densities [56], did not seem to cause significant variation in chemical rates between day and night. Accordingly, this species appears to rely predominantly on heterotrophic feeding to fuel bioerosion. Similarly, no significant diurnal pattern was found for any of the other species, although chemical dissolution by these species does appear to be somewhat lower during the day. Possibly, these species benefited from the small decrease in sea water pH at night. Indeed, photosynthesis during the day-time and respiration at night by primary producers is known to cause diurnal variation in pH [78]. The latter suggests that these sponges might respond more firmly than symbiont-rich species to a decreasing sea water pH as a result of ocean acidification. This corresponds to the findings of Webb et al. (2017) [34] that at night dissolution by C. caribbaea only started increasing at a strongly elevated pCO2 of the incubation water. The absence of a clear diurnal pattern in azooxanthellate species is in compliance with findings of Schönberg (2008) [79] who studied rates of Cliona celata in the North Sea, but also with Zundelevich et al. (2007) [36] who studied the zooxanthellate Pione vastifica in the Red Sea.

Rates obtained through FT incubation show that mechanical removal of invaded substratum contributed most significantly (86–95%) to total estimated bioerosion rates. Only for S. flavolivescens (network-forming) chemical and mechanical erosion was approximately equal. A detailed visual inspection of S. flavolivescens specimens, however, revealed a large quantity of CaCO3 particles to be present in its dense tissue network. It is possible that this species retains most of the mechanically removed fragments, causing an underestimation of mechanical erosion when using methods as presented here. With the exception of S. flavolivescens the ratio between dissolution and mechanical substrate removal resembles that described in the majority of the published literature, where it is suggested that mechanical dissolution contributes between 83% to 98% to total bioerosion by various studied excavating sponge species (e.g. [29, 32, 3335]). Based on our findings there appears to be no obvious variation in this ratio among species adapting different boring strategies. However, we did observe some diurnal variation in the mech:chem ratio, with chemical erosion having a slightly more pronounced contribution to total erosion at night. This was especially clear for C. delitrix (chemical contribution: day: 8%, night 26%) and might be a result of decreased pH at night-time. This suggests that the mech:chem ration might shift with the predicted ocean acidification scenarios. The portion of mechanically removed fragments that could be attributed to the characteristic sponge chips was limited compared to the total collected CaCO3 for all six species. This was particularly evident in the three cavity-forming species (C. amplicavata C. delitrix and S. brevitubulatum) and S. flavolivescens, where chips accounted for a mere 17%–24%. CaCO3 material produced by C. aprica (74%) and C. caribbaea (64%) contained considerably more sponge chips, suggesting chip production may be linked to the adopted excavating strategy. Although chips are with great certainty removed by the sponge itself, many of the remaining CaCO3 fragments are most likely also the dislodged as a result of sponge boring activity since substantially less of this material was present on filters of uncolonized cores.

The sum of chemical and mechanical rates (both FT) realized final estimates of total annual bioerosion for all six studies species: 10.95 kg m-2 y-1 (C. amplicavata); 4.96 kg m-2 y-1 (C. aprica); 8.25 kg m-2 y-1 (C. caribbaea); 12.52 kg m-2 y-1 (C. delitrix); 7.23 kg m-2 y-1 (S. brevitubulatum) and 2.01 kg m-2 y-1 (S. flavolivescens). With the exception of C. amplicavata all rates presented here is based on bioerosion of Orbicella spp. coral substratum which is of ‘average’ density and the most ubiquitous substratum on many reefs throughout the Caribbean. The considerable variation in published bioerosion rates ([40]: S3 Table) can largely be attributed to factors such as the density of the substratum (e.g. [41, 42]), but also studied species, developmental stage and applied methodology. Acker and Risk (1985) [33] provided total erosion rates for C. caribbaea in Orbicella spp. (8.0 kg m-2 y-1) that are comparable to rates presented here. More recently, a considerably lower rate (1.7 kg m-2 y-1) for total erosion by C. caribbaea was described for C. caribbaea, but in less dense Diploria spp. substratum and based on 6 h standing incubation [34]. Similar rates (2.2–2.5 kg m-2 y-1) were described for its ‘Pacific equivalent’ C. orientalis in massive Porites spp. coral [35, 52]. Both Fang et al. (2013b) [35] and Webb et al. (2017) [34], however, found 3 to 5 times higher rates using the buoyant weight technique which results in rates that are similar to our findings. Estimates described for C. aprica range higher (7.0 kg m-2 y-1) than our findings but it should be noted that this is based on erosion of extremely dense Conch shells [44]. Substantially higher bioerosion rates were found for other C. viridis-complex species in substratum that is substantially denser than that of Orbicella spp. skeleton: C. orientalis (up to 17.6 kg m-2 y-1) [42] and C. varians (up to 22.8 kg m-2 y-1) [65]. Of the sponges studied here, C. delitrix is generally considered among the most aggressive eroding species [72, 80] and our data appears to support this claim. The observed eroding capacity of C. delitrix in the α-stage (this study) infers the bioeroding potential of C. delitrix since the β-stage is deemed considerably more destructive [81]. Next to C. delitrix, substantially high rates were only found for C. amplicavata, but in Acropora spp. rubble which has a considerably denser skeletal structure than Orbicella spp. Murphy et al. (2016) [27] are the only study to have assessed the bioerosive capacity of S. brevitubulatum, adapting a completely different approach including substrate density and lateral expansion rate. Following their analogy a hypothetical sponge would bore with a rate of 4.4 kg m-2 y-1 in carbonate substrate with a density of 1.7 g cm-3. This is in range with our estimates in substratum of comparable density. It is noteworthy that the estimated rates for S. brevitubulatum as presented in this work may not fully reflect the boring capacity of this species considering that most of the soft internal tissue needed to be removed to prevent tissue decomposition. To our knowledge, the bioeroding capacity of C. amplicavata has not previously been determined. Also, no representative data is available for S. flavolivescens or any species that adapts a similar excavating strategy.

Overall, patterns of inter-specific variation were strikingly similar when related to either surface area or biomass of the sponge. This would suggest a consistent link between surface area and internal biomass of the sponge with regards to bioerosive activity, which is in contrast to the few existing preceding studies addressing this issue [40, 44, 82]. Nonetheless, it remains difficult to properly relate surface area to internal biomass of the sponge and we only quantified sponge biomass in three cores per species. Elaborate research on this specific relationship is crucial to properly estimate the contribution of excavating sponges to the reef carbonate budget.

We provide a comprehensive comparison of methods to quantify chemical and mechanical erosion and tentatively endorse methodologies that could be implemented widely in an attempt to support standardized data collection. The ‘novel’ FT technique described here is assumed to more adequately mimic natural in situ conditions, thus providing a more reliable representation of the bioerosive capacity of the excavating sponge community. We therefore encourage the use of flow-through incubation techniques next to buoyant weighing techniques, in particular when the aim is to determine the chemical fraction of bioerosion separately. If for some reason standing incubation techniques are preferred, the length of the incubation should be minimized or at least two time points should be included so that the signal at T0 can be predicted. It should be noted that determining mechanical erosion based on collected particulate material will always be prone to underestimation because loss of particles in the smallest size class can not be fully excluded. This is the first study to estimate rates for such an extensive range of species encompassing the relevant majority of the sponge bioeroding community of Curaçao and the wider Caribbean. This includes multiple α-stage species which, to date, have largely been ignored because their distinct papillate-morphology complicates linking assessed rates to the infested surface area [40]. Although we recognize that these data cover only part of the possible variation in rates of bioerosion caused by inter-specific and intra-specific variation, developmental stage, successional stage, density of the invaded substrate, reef-zone [43] and demography it can be considered an important addition to the currently still limited knowledge on the bioerosive impact of the excavating sponge community on coral reefs. The data presented here, combined with previously published and future work can also be included in standardized surveys, such as the ReefBudget methodology [39], resulting in a more complete and correct image of the carbonate budget of Caribbean coral reefs.

Supporting information

S1 Fig. Daily removed CaCO3 (mg) per gram of dry weight sponge tissue.

Provided rates are based on the loss in buoyant weight of sponge infested cores (n = 3 per species) over a three week period. Sponge tissue was dried for 24 h at 60°C. Rates are given for all six studied species: Cliona amplicavata (cavity-forming), Cliona aprica (gallery-forming), Cliona caribbaea (gallery-forming), Cliona delitrix (cavity-forming), Siphonodictyon brevitubulatum (cavity-forming) and Suberea flavolivescens (network-forming). All rates include a correction for the outcomes of uncolonized control core experiments. Error bars represent the 95% confidence limits.

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

(EPS)

S1 Table. Overview of the data supporting the results presented in this article.

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

(XLSX)

Acknowledgments

The authors would like to extend their sincere gratitude to Dr. Sven Zea for helping to identify the collected sponge species. We also wish to thank Karel Bakker and Sharyn Ossebaar for analysing nutrient samples, the Carmabi foundation for facilitating the fieldwork on Curaçao and the EMBC Master Programme. External funding for A.E. Webb and S.M.A.C. van Heuven was provided by the Netherlands Organization for Scientific Research (NWO grants 858.14.021 and 858.14.022). All research was conducted under research permit (#2012/48584) issued by the Curaçaoan Ministry of Health, Environment and Nature (GMN) to the CARMABI Foundation.

References

  1. 1. Glynn PW, Manzello DP. Bioerosion and coral reef growth: a dynamic balance. Coral reefs in the Anthropocene: Springer; 2015. p. 67–97.
  2. 2. Kleypas JA, Buddemeier RW, Gattuso J-P. The future of coral reefs in an age of global change. International Journal of Earth Sciences. 2001;90(2):426–37.
  3. 3. Goreau TF. The physiology of skeleton formation in corals. I. A method for measuring the rate of calcium deposition by corals under different conditions. The Biological Bulletin. 1959;116(1):59–75.
  4. 4. Hutchings P. Biological destruction of coral reefs. Coral reefs. 1986;4(4):239–52.
  5. 5. Andersson AJ, Gledhill D. Ocean acidification and coral reefs: effects on breakdown, dissolution, and net ecosystem calcification. Annual Review of Marine Science. 2013;5:321–48. pmid:22881351
  6. 6. Erez J, Reynaud S, Silverman J, Schneider K, Allemand D. Coral calcification under ocean acidification and global change. Coral reefs: an ecosystem in transition: Springer; 2011. p. 151–76.
  7. 7. Glynn PW, Morales GEL. Coral reefs of Huatulco, West Mexico: reef development in upwelling Gulf of Tehuantepec. Revista de Biología Tropical. 1997;45(3):1033–47.
  8. 8. Bellwood DR, Hughes TP, Folke C, Nyström M. Confronting the coral reef crisis. Nature. 2004;429(6994):827–33. pmid:15215854
  9. 9. Hughes TP, Baird AH, Bellwood DR, Card M, Connolly SR, Folke C, et al. Climate change, human impacts, and the resilience of coral reefs. science. 2003;301(5635):929–33. pmid:12920289
  10. 10. Bak RPM, Nieuwland G, Meesters EHWG. Coral reef crisis in deep and shallow reefs: 30 years of constancy and change in reefs of Curacao and Bonaire. Coral reefs. 2005;24(3):475–9.
  11. 11. De Bakker DM, Van Duyl FC, Bak RP, Nugues MM, Nieuwland G, Meesters EH. 40 Years of benthic community change on the Caribbean reefs of Curaçao and Bonaire: the rise of slimy cyanobacterial mats. Coral Reefs. 2017;36(2):355–67.
  12. 12. Hoegh-Guldberg O, Mumby PJ, Hooten AJ, Steneck RS, Greenfield P, Gomez E, et al. Coral reefs under rapid climate change and ocean acidification. science. 2007;318(5857):1737–42. pmid:18079392
  13. 13. Carpenter KE, Abrar M, Aeby G, Aronson RB, Banks S, Bruckner A, et al. One-third of reef-building corals face elevated extinction risk from climate change and local impacts. Science. 2008;321(5888):560–3. pmid:18653892
  14. 14. Jackson JBC, Kirby MX, Berger WH, Bjorndal KA, Botsford LW, Bourque BJ, et al. Historical overfishing and the recent collapse of coastal ecosystems. science. 2001;293(5530):629–37. pmid:11474098
  15. 15. Jackson JBC, Donovan M, Cramer K, Lam V. Status and trends of Caribbean coral reefs: 1970–2012. Global Coral Reef Monitoring Network, 2014.
  16. 16. Wilkinson C. Status of coral reefs of the world. 2008.
  17. 17. Gardner TA, Côté IM, Gill JA, Grant A, Watkinson AR. Long-term region-wide declines in Caribbean corals. Science. 2003;301(5635):958–60. pmid:12869698
  18. 18. Alvarez-Filip L, Dulvy NK, Côté IM, Watkinson AR, Gill JA. Coral identity underpins architectural complexity on Caribbean reefs. Ecol Appl. 2011;21(6):2223–31. pmid:21939056
  19. 19. Perry CT, Steneck RS, Murphy GN, Kench PS, Edinger EN, Smithers SG, et al. Regional‐scale dominance of non‐framework building corals on Caribbean reefs affects carbonate production and future reef growth. Glob Chang Biol. 2015;21(3):1153–64. pmid:25537577
  20. 20. Aronson RB, Precht WF. White-band disease and the changing face of Caribbean coral reefs. The ecology and etiology of newly emerging marine diseases: Springer; 2001. p. 25–38.
  21. 21. De Bakker DM, Meesters EH, Bak RPM, Nieuwland G, Van Duyl FC. Long-term Shifts in Coral Communities On Shallow to Deep Reef Slopes of Curaçao and Bonaire: Are There Any Winners? Frontiers in Marine Science. 2016;3(247).
  22. 22. Schönberg CH, Fang JK, Carreiro-Silva M, Tribollet A, Wisshak M. Bioerosion: the other ocean acidification problem. ICES Journal of Marine Science. 2017;74(4):895–925.
  23. 23. Sheppard CR, Spalding M, Bradshaw C, Wilson S. Erosion vs. recovery of coral reefs after 1998 El Niño: Chagos reefs, Indian Ocean. AMBIO: A Journal of the Human Environment. 2002;31(1):40–8.
  24. 24. Hudson J, editor Long-term bioerosion rates on a Florida reef: a new method. Proceedings of the 3rd International Coral Reef Symposium; 1977.
  25. 25. Rützler K. Impact of crustose clionid sponges on Caribbean reef corals. Acta Geol Hisp. 2002;37(1):61–72.
  26. 26. Macgeachy JK, Stearn CW. Boring by macro‐organisms in the coral Montastrea annularis on Barbados Reefs. International Review of Hydrobiology. 1976;61(6):715–45.
  27. 27. Murphy G, Perry C, Chin P, McCoy C. New approaches to quantifying bioerosion by endolithic sponge populations: applications to the coral reefs of Grand Cayman. Coral Reefs. 2016;35(3):1109–21.
  28. 28. Pang RK. The ecology of some jamaican excavating sponges. Bull Mar Sci. 1973;23.
  29. 29. Rützler K, Rieger G. Sponge burrowing: fine structure of Cliona lampa penetrating calcareous substrata. Marine Biology. 1973;21(2):144–62.
  30. 30. Pomponi SA. Cytological mechanisms of calcium carbonate excavation by boring sponges. International review of Cytology. 1980;65:301–19.
  31. 31. Cobb WR. Penetration of calcium carbonate substrates by the boring sponge, Cliona. American Zoologist. 1969;9(3):783–90.
  32. 32. Nava H, Carballo JL. Chemical and mechanical bioerosion of boring sponges from Mexican Pacific coral reefs. Journal of experimental Biology. 2008;211(17):2827–31.
  33. 33. Acker KL, Risk MJ. Substrate destruction and sediment production by the boring sponge Cliona caribbaea on Grand Cayman Island. Journal of Sedimentary Research. 1985;55(5).
  34. 34. Webb AE, van Heuven SMAC, de Bakker DM, van Duyl FC, Reichart G-J, de Nooijer LJ. Combined Effects of Experimental Acidification and Eutrophication on Reef Sponge Bioerosion Rates. Frontiers in Marine Science. 2017;4(311).
  35. 35. Fang JK, Mello‐Athayde MA, Schönberg CH, Kline DI, Hoegh‐Guldberg O, Dove S. Sponge biomass and bioerosion rates increase under ocean warming and acidification. Global change biology. 2013;19(12):3581–91. pmid:23893528
  36. 36. Zundelevich A, Lazar B, Ilan M. Chemical versus mechanical bioerosion of coral reefs by boring sponges-lessons from Pione cf. vastifica. Journal of experimental biology. 2007;210(1):91–6.
  37. 37. Carballo JL, Ovalle‐Beltrán H, Yáñez B, Bautista‐Guerrero E, Nava‐Bravo H. Assessment of the distribution of sponge chips in the sediment of East Pacific Ocean reefs. Marine Ecology. 2017;38(1).
  38. 38. Ogden JC. Carbonate-sediment production by parrot fish and sea urchins on Caribbean reefs. Studies in geology. 1977;4:281–8.
  39. 39. Perry C, Edinger E, Kench P, Murphy G, Smithers S, Steneck R, et al. Estimating rates of biologically driven coral reef framework production and erosion: a new census-based carbonate budget methodology and applications to the reefs of Bonaire. Coral Reefs. 2012;31(3):853–68.
  40. 40. Schönberg CH, Fang JK-H, Carballo JL. Bioeroding sponges and the future of coral reefs. Climate Change, Ocean Acidification and Sponges: Springer; 2017. p. 179–372.
  41. 41. Calcinai B, Azzini F, Bavestrello G, Gaggero L, Cerrano C. Excavating rates and boring pattern of Cliona albimarginata (Porifera: Clionaidae) in different substrata. Porifera Research: Biodiversity, Innovation & Sustainability. 2007:203–10.
  42. 42. Schönberg CH. Substrate effects on the bioeroding demosponge Cliona orientalis. 1. Bioerosion rates. Marine Ecology. 2002;23(4):313–26.
  43. 43. Bak R. The growth of coral colonies and the importance of crustose coralline algae and burrowing sponges in relation with carbonate accumulation. Netherlands Journal of Sea Research. 1976;10(3):285IN1293–292IN4337.
  44. 44. Rützler K. The role of burrowing sponges in bioerosion. Oecologia. 1975;19(3):203–16. pmid:28309235
  45. 45. Neumann AC. Observations on coastal erosion in Bermuda and measurements of the boring rate of the sponge, Cliona lampa. Limnology and Oceanography. 1966;11(1):92–108.
  46. 46. Fang JK, Schoenberg CH, Kline DI, Hoegh‐Guldberg O, Dove S. Methods to quantify components of the excavating sponge Cliona orientalis Thiele, 1900. Marine Ecology. 2013;34(2):193–206.
  47. 47. Wisshak M, Schönberg C, Form AU, Freiwald A. Effects of ocean acidification and global warming on reef bioerosion—lessons from a clionaid sponge. Aquatic Biology. 2013;19(2):111–27.
  48. 48. Achlatis M, van der Zande RM, Schönberg CH, Fang JK, Hoegh-Guldberg O, Dove S. Sponge bioerosion on changing reefs: ocean warming poses physiological constraints to the success of a photosymbiotic excavating sponge. Scientific Reports. 2017;7.
  49. 49. Wisshak M, Schönberg C, Form A, Freiwald A. Sponge bioerosion accelerated by ocean acidification across species and latitudes? Helgoland Marine Research. 2014;68(2):253.
  50. 50. Rix L, Goeij JM, Oevelen D, Struck U, Al‐Horani FA, Wild C, et al. Differential recycling of coral and algal dissolved organic matter via the sponge loop. Functional Ecology. 2017;31(3):778–89.
  51. 51. Van Duyl FC. Atlas of the Living Reefs of Cura «; ao and Bonaire (Netherlands Antilles). Foundation for Scientific Research in Surinam and the Netherlands Antilles. 1985;117.
  52. 52. Wisshak M, Schönberg CH, Form A, Freiwald A. Ocean acidification accelerates reef bioerosion. PloS one. 2012;7(9):e45124. pmid:23028797
  53. 53. Rützler K. The burrowing sponges of Bermuda: Smithsonian Institution Press; 1974.
  54. 54. Schönberg C, editor Sponges of the ‘Cliona viridis complex’–a key for species identification. Proceedings of 9th international coral reef symposium; 2002.
  55. 55. Zea S, Weil E. Taxonomy of the Caribbean excavating sponge species complex Cliona caribaea-C. aprica-C. langae (Porifera, Hadromerida, Clionaidae). Caribbean Journal of Science. 2003;39(3):348–70.
  56. 56. Granados C, Camargo C, Zea S, Sánchez J. Phylogenetic relationships among zooxanthellae (Symbiodinium) associated to excavating sponges (Cliona spp.) reveal an unexpected lineage in the Caribbean. Molecular phylogenetics and evolution. 2008;49(2):554–60. pmid:18725307
  57. 57. Chaves‐Fonnegra A, Feldheim KA, Secord J, Lopez JV. Population structure and dispersal of the coral‐excavating sponge Cliona delitrix. Molecular ecology. 2015;24(7):1447–66. pmid:25732628
  58. 58. Mueller B. First documentation of encrusting specimen of Cliona delitrix on Curaçao: a cause for concern? Marine Biodiversity Records. 2017;10(1):10.
  59. 59. Goreau TF. The ecology of Jamaican coral reefs I. Species composition and zonation. Ecology. 1959;40(1):67–90.
  60. 60. Bak RP. Coral reefs and their zonation in Netherlands Antilles: modern and ancient reefs. 1977.
  61. 61. Hughes TP. Skeletal density and growth form of corals. Marine Ecology Progress Series. 1987:259–66.
  62. 62. Davies PS. Short-term growth measurements of corals using an accurate buoyant weighing technique. Marine biology. 1989;101(3):389–95.
  63. 63. Bak R. Coral weight increment in situ. A new method to determine coral growth. Marine Biology. 1973;20(1):45–9.
  64. 64. Fang JK, Schönberg CH, Hoegh-Guldberg O, Dove S. Day–night ecophysiology of the photosymbiotic bioeroding sponge Cliona orientalis Thiele, 1900. Marine Biology. 2016;163(5):1–12.
  65. 65. Hill M. Symbiotic zooxanthellae enhance boring and growth rates of the tropical sponge Anthosigmella varians forma varians. Marine Biology. 1996;125(4):649–54.
  66. 66. Jacques T, Pilson M. Experimental ecology of the temperate scleractinian coral Astrangia danae I. Partition of respiration, photosynthesis and calcification between host and symbionts. Marine Biology. 1980;60(2–3):167–78.
  67. 67. Liu X, Byrne RH, Lindemuth M, Easley R, Mathis JT. An automated procedure for laboratory and shipboard spectrophotometric measurements of seawater alkalinity: Continuously monitored single-step acid additions. Marine Chemistry. 2015;174:141–6.
  68. 68. Breland JA, Byrne RH. Spectrophotometric procedures for determination of sea water alkalinity using bromocresol green. Deep Sea Research Part I: Oceanographic Research Papers. 1993;40(3):629–41.
  69. 69. Yao W, Byrne RH. Simplified seawater alkalinity analysis: use of linear array spectrometers. Deep Sea Research Part I: Oceanographic Research Papers. 1998;45(8):1383–92.
  70. 70. Hydes D, Aoyama M, Aminot A, Bakker K, Becker S, Coverly S, et al. Determination of dissolved nutrients (N, P, Si) in seawater with high precision and inter-comparability using gas-segmented continuous flow analysers. 2010.
  71. 71. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2013. ISBN 3-900051-07-0; 2013.
  72. 72. Chaves-Fonnegra A, Zea S. Observations on reef coral undermining by the Caribbean excavating sponge Cliona delitrix (Demospongiae, Hadromerida) Porifera research: biodiversity, innovation and sustainability. Série Livros. 2007;28.
  73. 73. Vacelet J, Boury-Esnault N. Carnivorous sponges. Nature. 1995;373(6512):333–5.
  74. 74. Mueller B, de Goeij JM, Vermeij MJ, Mulders Y, van der Ent E, Ribes M, et al. Natural diet of coral-excavating sponges consists mainly of dissolved organic carbon (DOC). PloS one. 2014;9(2):e90152. pmid:24587253
  75. 75. Fang JK, Schönberg CH, Mello‐Athayde MA, Hoegh‐Guldberg O, Dove S. Effects of ocean warming and acidification on the energy budget of an excavating sponge. Global change biology. 2014;20(4):1043–54. pmid:23966358
  76. 76. Weisz JB, Massaro AJ, Ramsby BD, Hill MS. Zooxanthellar symbionts shape host sponge trophic status through translocation of carbon. The Biological Bulletin. 2010;219(3):189–97. pmid:21183440
  77. 77. Schönberg C, editor Growth and erosion of the zooxanthellate Australian bioeroding sponge Cliona orientalis are enhanced in light. Proceedings of the 10th International Coral Reef Symposium; 2006.
  78. 78. Smith JE, Price NN, Nelson CE, Haas AF. Coupled changes in oxygen concentration and pH caused by metabolism of benthic coral reef organisms. Marine Biology. 2013;160(9):2437–47.
  79. 79. Schönberg CH. A history of sponge erosion: from past myths and hypotheses to recent approaches. Current developments in bioerosion: Springer; 2008. p. 165–202.
  80. 80. Rose CS, Risk MJ. Increase in Cliona delitrix infestation of Montastrea cavernosa heads on an organically polluted portion of the Grand Cayman fringing reef. Marine Ecology. 1985;6(4):345–63.
  81. 81. Chaves-Fonnegra A, Zea S. Coral colonization by the encrusting excavating Caribbean sponge Cliona delitrix. Mar Ecol. 2011;32.
  82. 82. Holmes G, Ortiz J-C, Schönberg CHL. Bioerosion rates of the sponge Cliona orientalis Thiele, 1900: spatial variation over short distances. Facies. 2009;55(2):203–11.