https://orcid.org/0000-0002-1868-9630
Figures
Abstract
The North Atlantic Oscillation (NAO) has been hypothesized to drive interannual variability in Bermudan coral extension rates and reef-scale calcification through the provisioning of nutritional pulses associated with negative NAO winters. However, the direct influence of the NAO on Bermudan coral calcification rates remains to be determined and may vary between species and reef sites owing to implicit differences in coral life history strategies and environmental gradients across the Bermuda reef platform. In this study, we investigated the connection between negative NAO winters and Bermudan Diploria labyrinthiformis, Pseudodiploria strigosa, and Orbicella franksi coral calcification rates across rim reef, lagoon, and nearshore reef sites. Linear mixed effects modeling detected an inverse correlation between D. labyrinthiformis calcification rates and the winter NAO index, with higher rates associated with increasingly negative NAO winters. Conversely, there were no detectable correlations between P. strigosa or O. franksi calcification rates and the winter NAO index suggesting that coral calcification responses associated with negative NAO winters could be species-specific. The correlation between coral calcification rates and winter NAO index was significantly more negative at the outer rim of the reef (Hog Reef) compared to a nearshore reef site (Whalebone Bay), possibly indicating differential influence of the NAO as a function of the distance from the reef edge. Furthermore, a negative calcification anomaly was observed in 100% of D. labyrinthiformis cores in association with the 1988 coral bleaching event with a subsequent positive calcification anomaly in 1989 indicating a post-bleaching recovery in calcification rates. These results highlight the importance of assessing variable interannual coral calcification responses between species and across inshore-offshore gradients to interannual atmospheric modes such as the NAO, thermal stress events, and potential interactions between ocean warming and availability of coral nutrition to improve projections for future coral calcification rates under climate change.
Citation: Courtney TA, Kindeberg T, Andersson AJ (2020) Coral calcification responses to the North Atlantic Oscillation and coral bleaching in Bermuda. PLoS ONE 15(11): e0241854. https://doi.org/10.1371/journal.pone.0241854
Editor: Erik Caroselli, University of Bologna, ITALY
Received: May 20, 2020; Accepted: October 21, 2020; Published: November 11, 2020
Copyright: © 2020 Courtney et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The online code and data repository is publicly accessible from GitHub: https://github.com/traviscourtney/BermudaCoring.
Funding: This work was funded by the University of California, San Diego Academic Senate via a Marine Sciences Research Grant (AJA) and the National Science Foundation via OCE–1829778 (AJA) and DGE–1650112 (TAC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
For millions of years, calcification by scleractinian corals has accounted for the majority of the growth and maintenance of tropical to subtropical shallow-water coral reef structures [1–3], which currently sustain approximately 10% of humanity with critical ecosystem services including shoreline protection, fisheries provisioning, and cultural significance [4–6]. The current range of optimal environmental conditions have limited coral dominated reef systems to approximately just 0.2% of the global oceans [7, 8]. Within this limited geographic range, considerable variation in coral calcification rates have been observed between species owing to varying ecological traits and within species owing to gradients in the environmental drivers of calcification through space and time [9]. However, ocean warming, acidification, eutrophication, and deoxygenation currently threaten the capacity for corals and coral reefs to maintain their calcium carbonate structures [9, 10]. Understanding the drivers of variability in calcification rates over broad spatial and temporal scales is therefore critical to inform effective management and conservation strategies [10]. Calcification rate time series data can be quantified from the product of linear extension and skeletal density of annual skeletal growth bands to evaluate the potential environmental drivers of coral calcification for many massive coral taxa [9, 11–14]. The high latitude coral reefs of Bermuda were chosen as the focus for this study because they contain coral taxa with observable growth bands and a naturally variable environment that may have conferred an increased resiliency of Bermudan coral communities to environmental change relative to their more southerly conspecifics [15–17].
Relatively few species of scleractinian corals occupy the high latitude coral reefs of Bermuda, but there are nonetheless considerable differences in calcification rates by up to an order of magnitude between species [18]. Within species, extension rates for Diploria labyrinthiformis, Porites astreoides, and Pseudodiploria strigosa have been observed to vary across the Bermuda platform with a tendency towards more rapid extension rates in nearshore environments relative to rim reef environments, which may be due to increased temperature extrema, elevated nutrition, or reduced wave activity closer to shore [19–21]. Further variability in D. labyrinthiformis and P. astreoides extension rates have been observed across depth gradients in Bermuda with lower extension rates observed for corals residing in deeper waters relative to shallower locations likely owing to reduced light intensity at deeper depths [20]. At the regional spatial scale, Bermudan coral growth rates are generally lower than their Caribbean conspecifics (e.g., Diploria labyrinthiformis, Pseudodiploria strigosa, Orbicella franksi, Porites astreoides, and Scolymia cubensis), which has been attributed to the on average cooler waters and reduced light intensity at higher latitudes [18–23].
Coral calcification rates in Bermuda have also been found to vary considerably over time. For example, P. astreoides and P. strigosa calcification rates have been shown to vary on diel time scales in association with the beneficial effects of light intensity on calcification rates [24]. Diploria labyrinthiformis and P. astreoides calcification rates were observed to vary on seasonal time scales primarily due to the positive effects of warmer sea surface temperature (SST) on calcification rates [16]. Diploria labyrinthiformis, P. strigosa, and Montastraea cavernosa extension rates and skeletal densities have been further observed to vary on interannual time scales in apparent synchrony with the North Atlantic Oscillation (NAO) [19, 25–28]. The NAO refers to the variability in the atmospheric pressure gradient observed at sea level between the subtropical high pressure system near the Azores and the arctic low pressure system near Iceland that drives climatic and oceanographic variability across the North Atlantic Ocean [29–32]. Notably, winters with a negative NAO index (i.e., hereafter referred to as negative NAO winters), are generally characterized by intensified water column mixing that entrain cold, nutrient-rich waters that drive enhanced phytoplankton blooms and community composition shifts in the North Atlantic subtropical gyre near Bermuda [29–34]. Elevated extension rates and reduced skeletal densities have been ascribed to increased coral nutrition (e.g., particulate organic matter [35]) coinciding with cooler SST during negative NAO winters [19, 25–28] with coral calcification rates assumed to remain approximately constant owing to the offsetting increases in extension rates and decreases in skeletal density [27]. However, summertime coral reef net ecosystem calcification on the Bermuda platform was found to increase following negative NAO winters [36] and food availability was positively correlated with increasing Bermudan coral calcification rates on seasonal time scales [16] suggesting negative NAO winters may also benefit Bermudan coral calcification rates.
The primary objective of this study was to test the hypothesis that negative NAO winters increase annual coral calcification rates across the Bermuda platform, which may be driven by a postulated increase in coral nutrition stimulating coral calcification rates on the reef [37]. We further hypothesize that negative NAO winters differ between species owing to differences in life history strategies, which can modulate calcification responses to environmental conditions [38]. Finally, we hypothesize that coral calcification rates and their responses to negative NAO winters vary across the platform, which may be due to the differential influence of the NAO interacting with environmental gradients across the Bermuda platform. To test these hypotheses, annual calcification rates from Diploria labyrinthiformis, Orbicella franksi, and Pseudodiploria strigosa coral cores collected across the Bermuda platform were analyzed with respect to the winter NAO phase. However, it is important to emphasize that the data presented here only assess the links between Bermudan coral calcification rates and the winter NAO index, and not the actual interannual nutritional availability on the reef, which remains a working hypothesis of the potential connection between the winter NAO index and coral calcification rates.
Methods
Species and site descriptions
Diploria labyrinthiformis, O. franksi, and P. strigosa were selected for coring in this study because they form massive skeletal morphologies with documented density bands [21, 39–41], use both tentacle and mucous-based feeding strategies [42], and represent some of the dominant reef-building corals in Bermuda [18]. Additionally, D. labyrinthiformis and P. strigosa are both considered as having stress-tolerant life histories whereas O. franksi is considered as having a generalist life history strategy, which suggests that their responses to environmental conditions may differ [38]. All coring sites were within the boundaries of the rim reefs surrounding Bermuda along a gradient of more oceanic to terrestrial proximity (Fig 1). The North Lagoon coring locations include an outer rim reef site (Hog Reef), two centrally located lagoon patch reef locations (Halfway Flat and Three Hill Shoals), and a nearshore reef site (Whalebone Bay, Fig 1). The south shore Gurnet Rock coring site is proximal to the open ocean and is also located at the tidal outflows of higher turbidity Castle Harbour inshore waters [43, 44] suggesting it may experience a mixture of nearshore and oceanic characteristics (Fig 1).
Colored shapes correspond to locations where three cores were collected from the respective species at each of the following reef sites: HR = Hog Reef; TH = Three Hill Shoals; WB = Whalebone Bay; GR = Gurnet Rock; and HF = Halfway Flat. Bathymetry data are from [45].
Limited observations support a gradient of decreasing seasonal temperature extrema, decreasing coral nutrition, and increasing wave energy from the nearshore to lagoon to rim reef locations [16, 19–21, 37, 43, 46–49]. For example, in situ seawater temperature ranges in 2006 decreased slightly from nearshore to rim reef locations with observations of 15.6–30.1°C for Whalebone Bay, 17.3–29.6°C for Crescent Reef (located near Halfway Flat), 17.1–29.4°C for Three Hill Shoals, and 17.9–29.0°C for Hog Reef [46]. Recorded salinity values have been reported to be relatively constant at 36.5 ppt across the majority of the Bermuda platform [50]. The Bermuda Inshore Waters Investigation provides some evidence for decreasing coral nutrition from nearshore to lagoon waters with higher nitrogen content of particulate matter (Great Sound: 4.0 μmol l–1 vs. North Lagoon: 2.8 μmol l–1), higher phosphorous content of particulate matter (Great Sound: 0.11 μmol l–1 vs. North Lagoon: 0.06 μmol l–1), and higher chlorophyll α concentrations (Great Sound: 0.78±0.37 μg l–1 vs. North Lagoon: 0.43±0.19 μg l–1) observed at the more nearshore Great Sound sampling location relative to the North Lagoon sampling location [49]. Higher turbidity was also observed for the more nearshore reef sites compared to the rim reef sites, suggesting a greater abundance of suspended particulate matter available for coral heterotrophy at nearshore relative to rim reef sites [37]. Mean flow speeds (ū) and turbulence intensities (q) on a calm day were lowest nearshore at Tynes Bay (ū = 1.64–2.03 cm s–1, q = 0.85–1.09 cm s–1), higher in the central lagoon at Crescent Reef (ū = 2.22–2.52 cm s–1, q = 1.03–1.11 cm s–1), and highest in the rim reef at Hog Reef (ū = 5.58–8.42 cm s–1, q = 2.49–4.49 cm s–1), which in part supports a gradient of increasing wave energy from nearshore to lagoon to rim reef locations [37]. We caution that further measurements across space and time are necessary to more strongly elucidate the persistence and variability in environmental parameters across the Bermuda platform while the data presented here nonetheless provide some evidence to support the presence of environmental gradients from the nearshore to rim reef environments.
Sample collection
All coral cores were collected in August 2016 using a pneumatic drill with 5-cm diameter wet concrete diamond core drill bits powered by a boat-based air compressor or underwater diving cylinders. Following extraction, cement plugs were placed in each colony and sealed with underwater epoxy to promote coral recovery. Each coral colony was selected for coring based on similar colony heights of approximately 30 cm to limit potential ontogenetic effects on growth rates [39–42], within a mean (±SD) depth of 4.5±1.8 m to limit depth effects [20], and with seemingly healthy appearances to reduce the effects of other potential confounding variables on growth rates [43]. Each coral core was collected along the vertical growth axis up to 30 cm in length from three replicate colonies of D. labyrinthiformis (n = 15), P. strigosa (n = 12), and O. franksi (n = 15) from four to five distinct reef locations for a total of 42 cores from 42 distinct coral colonies (Fig 1). Regrettably, fewer cores were collected from P. strigosa because no colonies meeting our search criteria were observed at the central platform reef site (Fig 1). The Government of Bermuda Department of Environment and Natural Resources is acknowledged for permitting collection and approving exports of the coral cores under license number 2015–05–26–46.
Quantification of linear extension, skeletal density, and calcification
Coral cores and density reference materials from the laboratory of Dr. Karl Castillo at the University of North Carolina at Chapel Hill were imaged via computed tomography (CT) with a GE Lightspeed VCT at San Diego Imagery. Individual cores and density standards were visualized as the mean voxel intensity of a ~3.00 mm thick digital slice using standardized window widths and window levels in the software Horos (version 2.4.1) and exported as an image for measuring growth rates in Coral XDS [51]. To quantify density from CT image intensity following established methods [52–55], linear models between known density and CT intensity of the density reference materials for three independent scans were derived using the function lm in the statistical package R version 4.0.2 [56]. Transects were analyzed along the major growth axis of individual polyps to quantify mean extension and intensity of subsequent high and low intensity bands using the half range mode and converted to skeletal density using the density calibration curves for each scan. Annual extension was determined as the sum of extension for adjacent high and low intensity bands while annual skeletal density was determined as the weighted average of skeletal densities for adjacent high and low skeletal density bands to account for any potential differences in extension rates between bands. Annual calcification rates (g cm–2 yr–1) were determined from the product of annual extension (cm yr–1) and skeletal density (g cm–3) [11, 13, 14, 57]. The 2016 growth bands were incompletely formed due to the timing of collection and were removed from the data prior to subsequent analyses of annual calcification rates.
Time series construction and pointer year analyses
To visualize temporal trends, mean chronologies (±95% confidence intervals) were developed for the annual linear extension, skeletal density, and calcification data from each coral core (S2 Fig) and the monthly winter NAO (DJFM = December, January, February, and March) index from the NOAA National Weather Service (i.e., 1950–2015, URL: https://www.cpc.ncep.noaa.gov/products/precip/CWlink/pna/nao.shtml) using the function MeanCI in the R package DescTools [58]. A pointer year analysis was used to test whether individual negative NAO winters corresponded to detectable deviations in annual coral linear extension, skeletal density, and calcification. In this study, pointer years were defined as individual years that exhibited ±10% positive or negative growth rate anomalies relative to the previous year in at least 75% of the cores for each year [59, 60] using the R package dplR [61]. Pointer years with less than one third of the full series of cores from each species were excluded to avoid potential false detection of growth rate anomalies based on a limited number of cores. The pointer year analyses were conducted separately for all cores collected from each species (e.g., D. labyrinthiformis [n = 15], P. strigosa [n = 12], and O. franksi [n = 15]) to test for species-level effects of negative NAO winters on annual linear extension, skeletal density, and calcification.
Linear mixed effects modeling
Linear mixed effects modeling was used to detect the fixed effects of winter NAO index, species, site, and the interactions between the winter NAO index and species or site on annual calcification rates for the duration of the NAO index (i.e., 1950–2015) using the R package nlme [62]. This analytical approach was used to account for nesting of calcification rates within each core and data imbalances generated by variable lengths of core time series data with the capacity to detect small trends in coral calcification rates following established methods [53, 55, 63]. Briefly, each coral core was treated as an individual sampling unit with annual calcification rates nested within each core. Random effects were incorporated to allow the estimated slopes and intercepts to vary by core. A stepwise modeling approach building from simple null models to more complex models utilizing random and fixed effects was followed with minimized AIC values used to select the optimal model. Significance tests of the main effects and interactions for the optimal model were evaluated with ANOVA (Type III) in the anova.lme function in R package nlme [62]. Simple slopes for site and species that incorporate the main effects and interactions and Tukey Honest Significant Difference (HSD) pairwise comparisons of species and site were determined using the function emtrends in R package emmeans [64].
Results
Density calibrations
The measured densities of coral core standards were strongly correlated with CT scan intensity for three independent scans (Scan 1: Slope = 0.019±0.002, adjusted R2 = 0.92, p<0.001; Scan 2: Slope = 0.018±0.001, adjusted R2 = 0.97, p<0.001; Scan 3: Slope = 0.0205±0.0009, adjusted R2 = 0.98, p<0.001). The 95% confidence intervals of the predicted mean skeletal densities for all cores based on these regressions generates a mean calculated skeletal density uncertainty of ±0.06 g cm–3 (S1 Fig).
Mean chronologies and pointer years
The winter NAO index exhibited interannual variability with a tendency towards more negative phases from approximately 1950–1980 and more positive phases thereafter with the exception of a notably negative NAO winter in 2010 (Fig 2a). There was a high degree of variability in calcification rates between cores and years that tended to increase towards the base of the time series (Fig 2b–2d) as the number of cores with growth rate data for the respective year declined (Fig 2e). No apparent long-term calcification trends were observed except for D. labyrinthiformis, which appeared to be slightly elevated from approximately 1950–1985 (Fig 2b–2d). Coral core chronologies with respect to the mean annual extension, skeletal density, and calcification data for each species reflected similar trends with calcification appearing to more closely follow patterns in linear extension than skeletal density (S2 Fig). There was a recorded negative pointer year in 1988 (100% of cores exhibited negative growth relative to the previous year) and a positive pointer year in 1989 (83% of cores exhibited positive growth relative to the previous year) in the annual calcification time series data for D. labyrinthiformis (Fig 2). No other pointer years were detected in the linear extension, skeletal density, or calcification rate data (S2 Fig).
(a) Black line indicates mean of NAO index for DJFM with shaded ±95% confidence intervals in light blue and dashed-red line indicating zero. (b,c,d) Gray lines indicate calcification rates (g cm–2 yr–1) for individual coral cores and the black lines with light blue shading represent mean calcification rates ±95% confidence intervals for each species. The vertical red line indicates the negative pointer year and the vertical light blue line indicates the positive pointer year. (e) Lines indicate the number of cores with calcification rate data for the respective year from each species.
Linear mixed effects modeling
The random slopes and intercepts model with winter NAO index, species, site, interaction of winter NAO index with species, and interaction of winter NAO index with site as fixed effects was selected as the best model to fit the annual calcification rates (Table 1). Type III ANOVA revealed that species (p = 0.0002) and the interactions of winter NAO index with species (p = 0.0130) and site (p = 0.0319) were significant predictors of calcification rates while the main effects of winter NAO index (p = 0.0990) and site (p = 0.0960) were not (Table 2). The simple slope (±SE) for calcification per unit winter NAO index was negative for D. labyrinthiformis (–0.013±0.004 g cm–2 per NAO) and significantly lower than the slopes for O. franksi (0.003±0.005 g cm–2 yr–1 per NAO) and P. strigosa (0.003±0.005 g cm–2 yr–1 per NAO), which exhibited no detectable trends or differences between each other (Fig 3a, Table 3). Similarly, there were no detectable trends in simple slopes in calcification with respect to the NAO across all sites (i.e., 95% confidence intervals overlap zero), but simple slopes at Hog Reef were significantly more negative than for Whalebone Bay (Fig 3b, Table 3).
Squares represent mean slope value with shaded bars representing ± 95% confidence intervals and the red dashed line indicates zero slope. The letters overlying each confidence interval represents statistical significance from post-hoc Tukey Honest Significant Difference (HSD) tests between species (a) and sites (b). In this notation, slopes identified by the same letters are not significantly different from each other and slopes with different letters are significantly different from each other.
Discussion
The negative correlation between D. labyrinthiformis calcification rates and winter NAO index partially supports the hypothesis that stronger negative winter NAO phases influence annual Bermudan coral (Fig 3a) calcification rates [36]. This extends the established links between coral extension rates and skeletal densities [19, 25, 26, 28] to include enhanced coral calcification rates associated with negative NAO winters, which are characterized by cooler SST, elevated nutrients, and enhanced phytoplankton blooms in offshore waters [29–34]. However, there were no detectable correlations between O. franksi or P. strigosa calcification rates and the winter NAO index (Fig 3a), which may suggest that calcification responses to NAO linked changes vary by species. While we cannot unequivocally rule out the possibility that the increased coverage of calcification time series data in earlier years preferentially enabled detection of the winter NAO index impacts on D. labyrinthiformis calcification rates, all three species have calcification rates overlapping with nearly the entire 1950–2015 winter NAO index suggesting we could still expect to detect any real correlations with the NAO (Fig 2). Negative NAO winters may indeed be more beneficial to D. labyrinthiformis calcification rates than O. franksi and P. strigosa (Fig 3a), but the mechanisms for this potential difference such as potential provisioning of nutritional pulses to the Bermudan corals during negative NAO winters remains to be elucidated. For example, heterotrophic status may be more similar for D. labyrinthiformis and P. strigosa on the basis of similar suspension feeding rates [65] and tissue thickness (S1 Table) that were greater than for O. franksi. However, O. franksi and P. strigosa exhibited statistically similar nitrogen uptake of particulate matter in a controlled feeding experiment [37] and had similar calcification responses to the winter NAO index in this study (Fig 3a). Future studies are needed to more firmly establish the mechanisms for the species-specific calcification responses to the winter NAO observed herein (Fig 3a). Critically, such studies would further elucidate the species-specific effects of heterotrophy on Bermudan coral calcification rates and the inshore-offshore coupling in phyto- and zooplankton biomass with respect to the NAO.
Spatial gradients in nutrition with respect to coral growth rates have long been hypothesized for the Bermuda platform [19, 21, 37]; however, environmental data to support the higher concentration of available coral nutrition on nearshore reef sites is sparse (e.g., see Species and site descriptions section) and warrants further investigation [16, 19–21, 37, 43, 46–49]. While no individual sites had detectable correlations between calcification and the winter NAO index, the correlation for the most outer rim reef coring site, Hog Reef, was significantly more negative than for the most nearshore coring site, Whalebone Bay (Fig 3, Table 3). There is therefore some tentative evidence that the negative NAO index may be more beneficial to coral calcification rates at the outer rim reefs than nearshore reefs. This may be due to the stronger influence of postulated nutritional pulses in negative NAO winters along the rim reef compared to nearshore because of lower background concentrations of nutrition at the former [16, 43, 46–48]. However, further research is needed to establish the mechanisms for any potential spatial differences in the response of coral clacification rates to the winter NAO index across the Bermuda platform. For example, site was not a significant predictor of calcification rates in this study (p = 0.0960, Fig 3b, Table 2), but previous studies with greater sample sizes (e.g., n = 7–18; [21]; n = 1–34 [19]) were able to detect elevated extension rates at sites more proximal to shore relative to sites closer to the outer rim reef [19, 21]. It is important to note that while calcification rates may respond more strongly to environmental conditions than linear extension [66], uncertainties in the measurement of calcification rates may be greater than for linear extension [67]. Therefore, an increased sampling size beyond the three corals per species per site in this study may be needed to further elucidate site-level differences in Bermudan coral calcification rates and their responses to the winter NAO index.
In contrast to the overall correlation between D. labyrinthiformis calcification rates and winter NAO, the identified negative pointer year in 1988 (winter NAO index = 0.48) and positive pointer year in 1989 (winter NAO index = 1.41) (Fig 2b), do not support enhanced coral growth rates associated with negative NAO winters. Instead, the reduced calcification rates for 100% of D. labyrinthiformis cores in 1988 coincided with the year of the first reported coral bleaching event in Bermuda [68] and suggest that the 1988 bleaching event impacted D. labyrinthiformis calcification rates more negatively than O. franksi and P. strigosa. The positive pointer year in 1989 and absence of any noticeable post-bleaching depression of annual calcification rates [69, 70] in the time series data (Fig 2b) suggests a rapid recovery in D. labyrinthiformis calcification rates. Furthermore, high density skeletal growth anomalies have been associated with coral bleaching events [70–79], but no positive pointer years were detectable in the skeletal density data of this study (S2 Fig). While subsequent minor coral bleaching events have occurred throughout the Bermuda platform in the years following the 1988 bleaching event [15], there were no detectable impacts of these additional events on calcification rates in the pointer year analyses of this study (Fig 2b–2d). The single negative pointer year observed in 1988 suggests that the coincident coral bleaching event may have been the most intense bleaching event on record. Further research is needed to determine whether this is due to potential differences in the magnitude of subsequent marine heat waves in Bermuda or potential long-term acclimatization of Bermudan corals following the 1988 bleaching event.
The high degree of stability of Bermudan coral communities over recent decades is evidenced by the lack of any long-term declines in coral cover or changes in the coral community composition [15–17] that have been observed for many of the lower-latitude Western Atlantic coral reefs [17, 38, 80–83]. However, it remains unclear whether this stability is owing to Bermudan corals thus far escaping or resisting thermal stress events [16] as unprecedented marine heatwaves have driven extensive coral mortality in other high-latitude coral reef systems with previously stable coral cover [84]. Interestingly, previous work in Bermuda has indicated that the thermal optima for D. labyrinthiformis [16] and O. franksi [85] are not currently exceeded by typical summer SST, but that calcification rates under future warming and acidification may become increasingly dependent on the continued availability of adequate nutrition [16, 85]. Warming winters have been hypothesized to benefit high latitude corals through warming calcification rate-limiting winter SST [9, 16]. However, increases in SST (i.e., +0.021°C yr–1) have coincided with declines in offshore primary productivity (i.e., –5.6 mg C m–2 d–1 per year) from 1990–2016 in the offshore waters of Bermuda suggesting these sources of nutrition may decrease with continued warming [86]. The mechanistic links between greenhouse gas emissions and the NAO remain equivocal [29, 31, 87, 88] making it difficult to project how the NAO may interact with climate change to influence offshore primary productivity and Bermudan coral calcification rates. We posit that if warming-induced declines in offshore nutrition persist and they are indeed an important source of Bermudan coral nutrition, then this will negatively affect future Bermudan and other high-latitude coral calcification rates, but the relative balance between the positive effects of winter warming and negative effects of reduced nutrition on calcification rates remains unclear (Fig 4). Meanwhile, the effects of warming on calcification rates in lower latitude corals is hypothesized to be more negative owing to summer SST exceeding coral thermal optima [9] with coinciding declines in nutrition likely to reduce the capacity of corals to resist and recover from the projected increases in coral bleaching events [89–91] (but see also the potential role of nutrients in coral bleaching events [92]). Regardless, reductions in the rates of greenhouse gas emissions under the Paris Agreement remain critical to reduce the probability of frequent and severe coral bleaching events in Bermuda [16] and worldwide [93].
The figure is modeled after [36] and the brain coral illustration (not drawn to scale) is courtesy of the integration and application network, university of Maryland center for environmental science (ian.umces.edu/symbols/).
Conclusions
Collectively the findings of this study highlight the complex interplay between coral species and site on calcification responses associated with the NAO and a thermal stress event. Calcification rates of D. labyrinthiformis responded most strongly to the winter NAO index and the marine heatwave associated with the 1988 coral bleaching event, while O. franksi and P. strigosa exhibited no detectable interannual calcification responses to the NAO or the 1988 bleaching event. Coral calcification responses were more sensitive to negative NAO winters at Hog Reef compared to Whalebone Bay, lending provisional support of potential platform-scale trends in calcification responses to the NAO. In addition to improved Bermudan coral futures under lower CO2 emission pathways [16, 85], further research should therefore be conducted to elucidate the mechanisms for the variable calcification responses to the NAO, marine heatwaves, and the potential role of coral nutrition between species and sites in Bermuda and throughout corals’ geographic ranges.
Supporting information
S1 Fig. Coral core standards density curves are plotted for each scan.
Each circle represents the mean CT-scan luminance and measured density for each coral standard. Solid black line represents linear regression of measured density vs. CT-scan luminance with shaded gray regions indicating ±95% confidence intervals of regression data.
https://doi.org/10.1371/journal.pone.0241854.s001
(TIFF)
S2 Fig. Coral core linear extension, skeletal density, and calcification rate time series data are plotted for each species.
Gray lines indicate growth measurements for individual coral cores, black lines represent mean coral growth measurements for each species, and shaded blue regions indicate ±95% confidence intervals of the species-level mean time series data. The vertical red line indicates a negative pointer year and vertical light blue line indicates a positive pointer year.
https://doi.org/10.1371/journal.pone.0241854.s002
(TIFF)
S1 Table. Mean (± standard deviation) core physiological measurements of the entire cored material and the number of cores collected are summarized for all cores from each species and separately for each species from each site.
Tissue thickness was measured on the day of collection for each core using Vernier calipers following established methods [94]. Mean coral extension for each core was determined as the mean of the top five years of the coral growth analysis (i.e., 2010–2015) to compare the most recent overlapping growth intervals. Core skeletal densities were determined as the mean skeletal density for the entire cored material to account for variations in skeletal densities between skeletal components. Mean calcification rates for each core were determined as the product of the most recent 5-year mean extension rate and mean core skeletal density [9, 11–14]. Mean annual linear extension, skeletal density, and calcification rates were computed for the entire cored skeletal material to provide context for future physiological studies and carbonate production estimates.
https://doi.org/10.1371/journal.pone.0241854.s003
(DOCX)
Acknowledgments
We thank Dr. Jessica Carilli for guidance on coral drilling methods, Richard Walsh for assistance with compressed air fittings, Dr. Justin Baumann for guidance on analyzing coral skeletal densities and discussions regarding coral heterotrophy, Dr. Karl Castillo for loaning coral density standards, Donna Hutchinson at San Diego Imaging for scanning the cores, and Lark Starkey for copy editing the manuscript. We thank the staff at the Bermuda Institute of Ocean Sciences for logistical support and assistance with collection permits. The Government of Bermuda Department of Environment and Natural Resources is acknowledged for permitting collection and approving exports of the coral cores under license number 2015–05–26–46. We further thank the reviewers for improving this manuscript with their constructive comments and insights.
References
- 1. Stanley GD. Early history of scleractinian corals and its geological consequences. Geology. 1981;9: 507–511.
- 2. Goreau T. Calcium carbonate deposition by coralline algae and corals in relation to their roles as reef-builders. Annals of the New York Academy of Sciences. 1963;109: 127–167. pmid:13949254
- 3. Hart DE, Kench PS. Carbonate production of an emergent reef platform, Warraber Island, Torres Strait, Australia. Coral Reefs. 2007;26: 53–68.
- 4. Moberg F, Folke C. Ecological goods and services of coral reef ecosystems. Ecological Economics. 1999;29: 215–233.
- 5. Donner SD, Potere D. The Inequity of the Global Threat to Coral Reefs. BioScience. 2007;57: 214–215.
- 6. Woodhead AJ, Hicks CC, Norström AV, Williams GJ, Graham NAJ. Coral reef ecosystem services in the Anthropocene. Functional Ecology. 2019;33: 1023–1034.
- 7. Reaka-Kudla ML. The Global Biodiversity of Coral Reefs: A Comparison with Rain Forests. Biodiversity II: Understanding and protecting our biological resources. 1997. pp. 83–108.
- 8. Kleypas JA, McManus JW, Meñez LAB. Environmental limits to coral reef development: where do we draw the line? American Zoologist. 1999;39: 146–159.
- 9. Pratchett MS, Anderson KD, Hoogenboom MO, Widman E, Baird AH, Pandolfi JM, et al. Spatial, temporal and taxonomic variation in coral growth—implications for the structure and function of coral reef ecosystems. Oceanography and Marine Biology: An Annual Review. 2015;53: 215–296.
- 10. Andersson AJ, Kline DI, Edmunds PJ, Archer SD, Bednaršek N, Carpenter RC, et al. Understanding Ocean Acidification Impacts on Organismal to Ecological Scales. Oceanography. 2015;28: 16–27.
- 11. Knutson DW, Buddemeier RW, Smith SV. Coral Chronometers: Seasonal Growth Bands in Reef Corals. Science. 1972;177: 270–272. pmid:17815626
- 12. Buddemeier RW, Maragos JE, Knutson DW. Radiographic studies of reef coral exoskeletons: Rates and patterns of coral growth. Journal of Experimental Marine Biology and Ecology. 1974;14: 179–199.
- 13. Dodge RE, Brass GE. Skeletal extension, density, and calcification of the reef coral, Montastrea annularis: St. Croix, U.S. Virgin Islands. Bulletin of Marine Science. 1984;34: 288–307.
- 14. Barnes DJ, Lough JM. On the nature and causes of density banding in massive coral skeletons. Journal of Experimental Marine Biology and Ecology. 1993;167: 91–108.
- 15.
Smith S, Sarkis S, Murdoch T, Weil E, Croquer A, Bates N, et al. Threats to Coral Reefs of Bermuda. In: Sheppard C, editor. Coral Reefs of the United Kingdom Overseas Territories. Dordrecht: Springer Science+Business Media; 2013. pp. 173–188. https://doi.org/10.1007/978-94-007-5965-7
- 16. Courtney TA, Lebrato M, Bates NR, Collins A, de Putron SJ, Garley R, et al. Environmental controls on modern scleractinian coral and reef-scale calcification. Science Advances. 2017;3: e1701356. pmid:29134196
- 17.
Jackson JBC, Donovan MK, Cramer KL, Lam V, Lam W. Status and Trends of Caribbean Coral Reefs: 1970–2012. Gland: Global Coral Reef Monitoring Network, IUCN; 2014 p. 306.
- 18. Courtney TA, Andersson AJ, Bates NR, Collins A, Cyronak T, de Putron SJ, et al. Comparing Chemistry and Census-Based Estimates of Net Ecosystem Calcification on a Rim Reef in Bermuda. Frontiers in Marine Science. 2016;3: 181.
- 19. Dodge R. The Natural Growth Records of Reef-Building Corals. Doctoral dissertation Nova Southeastern University Retrieved from NSUWorks. 1978; 237.
- 20. Logan A, Tomascik T. Extension growth rates in two coral species from high-latitude reefs of Bermuda. Coral Reefs. 1991;10: 155–160.
- 21. Logan A, Yang L, Tomascik T. Linear skeletal extension rates in two species of Diploria from high-latitude reefs in Bermuda. Coral Reefs. 1994;13: 225–230.
- 22.
Iams WJ. New methods for studying the growth rates of reef-building organisms. Bermuda Biological Station Special Publication 2. 1969; 65–75.
- 23. Tomascik T, Logan A. A comparison of peripheral growth rates in the recent solitary coral Scolymia cubensis (Milne-Edwards and Haime) from Barbados and Bermuda. Bulletin of Marine Science. 1990;46: 8.
- 24. Venti A, Andersson A, Langdon C. Multiple driving factors explain spatial and temporal variability in coral calcification rates on the Bermuda platform. Coral Reefs. 2014;33: 979–997.
- 25. Dodge RE, Vaišnys JR. Hermatypic coral growth banding as environmental recorder. Nature. 1975;258: 706–708.
- 26.
Pätzold J, Wefer G. Bermuda (32°N, 65°W): SST in the North Atlantic Gyre. In: Dunbar RB, Cole JE (Organizers) Coral records of ocean-atmosphere variability. Report from the workshop on coral paleoclimate reconstruction. La Parguera, Puerto Rico. 1992; 24–26.
- 27. Draschba S, Pätzold J, Wefer G. North Atlantic climate variability since AD 1350 recorded in δ18O and skeletal density of Bermuda corals. International Journal of Earth Sciences. 2000;88: 733–741.
- 28. Cohen AL, Smith SR, McCartney MS, Van Etten J. How brain corals record climate: An integration of skeletal structure, growth and chemistry of Diploria labyrinthiformis from Bermuda. Marine Ecology Progress Series. 2004;271: 147–158.
- 29. Visbeck MH, Hurrell JW, Polvani L, Cullen HM. The North Atlantic Oscillation: Past, present, and future. Proceedings of the National Academy of Sciences. 2001;98: 12876–12877. pmid:11687629
- 30. Gruber N, Keeling CD, Bates NR. Interannual Variability in the North Atlantic Ocean Carbon Sink. Science. 2002;298: 2374–2378. pmid:12493911
- 31. Hurrell JW, Deser C. North Atlantic climate variability: The role of the North Atlantic Oscillation. Journal of Marine Systems. 2009;78: 28–41.
- 32. Delworth TL, Zeng F, Vecchi GA, Yang X, Zhang L, Zhang R. The North Atlantic Oscillation as a driver of rapid climate change in the Northern Hemisphere. Nature Geoscience 2016;9: 509–512.
- 33. Menzel DW, Ryther JH. Annual variations in primary production of the Sargasso sea off Bermuda. Deep Sea Research (1953). 1961;7: 282–288.
- 34. Casey JR, Aucan JP, Goldberg SR, Lomas MW. Changes in partitioning of carbon amongst photosynthetic pico- and nano-plankton groups in the Sargasso Sea in response to changes in the North Atlantic Oscillation. Deep Sea Research Part II: Topical Studies in Oceanography. 2013;93: 58–70.
- 35. Houlbrèque F, Ferrier-Pagès C. Heterotrophy in tropical scleractinian corals. Biological Reviews. 2009;84: 1–17. pmid:19046402
- 36. Yeakel KL, Andersson AJ, Bates NR, Noyes TJ, Collins A, Garley R. Shifts in coral reef biogeochemistry and resulting acidification linked to offshore productivity. Proceedings of the National Academy of Sciences. 2015;112: 14512–14517. pmid:26553977
- 37. Mills MM, Lipschultz F, Sebens KP. Particulate matter ingestion and associated nitrogen uptake by four species of scleractinian corals. Coral Reefs. 2004;23: 311–323.
- 38. Darling ES, Alvarez-Filip L, Oliver TA, Mcclanahan TR, Cote IM. Evaluating life-history strategies of reef corals from species traits. Ecology Letters. 2012;15: 1378–1386. pmid:22938190
- 39. Saenger C, Cohen AL, Oppo DW, Hubbard D. Interpreting sea surface temperature from strontium/calcium ratios in Montastrea corals: Link with growth rate and implications for proxy reconstructions. Paleoceanography. 2008;23.
- 40. Dodge RE, Thomson J. The natural radiochemical and growth records in contemporary hermatypic corals from the Atlantic and Caribbean. Earth and Planetary Science Letters. 1974;23: 313–322.
- 41. Guzmán HM, Jackson JBC, Weil E. Short-term ecological consequences of a major oil spill on Panamanian subtidal reef corals. Coral Reefs. 1991;10: 1–12.
- 42. Lewis JB, Price WS. Feeding mechanisms and feeding strategies of Atlantic reef corals. Journal of Zoology. 1975;176: 527–544.
- 43.
Morris B, Barnes J, Brown F, Markham J. The Bermuda marine environment: a report of the Bermuda inshore waters investigation 1976–1977. Bermuda Biological Station Special Publication 15. 1977.
- 44. Flood VS, Pitt JM, Smith SR. Historical and ecological analysis of coral communities in Castle Harbour (Bermuda) after more than a century of environmental perturbation. Marine Pollution Bulletin. 2005;51: 545–557. pmid:16242159
- 45.
Sutherland MG, McLean SJ, Love MR, Carignan KS, Eakins BW. Digital Elevation Models of Bermuda: Data Sources, Processing and Analysis. Boulder, CO: NOAA National Geophysical Data Center, US Dept of Commerce. 2014;7.
- 46.
Jones RJ. Bermuda Institute of Ocean Sciences (BIOS). Marine Environmental Program (MEP) Annual Report (2006–2007), submitted to the Bermuda Government Department of Environmental Protection, Ministry of the Environment. 2006 p. 79.
- 47. Wang XT, Sigman DM, Cohen AL, Sinclair DJ, Sherrell RM, Weigand MA, et al. Isotopic composition of skeleton-bound organic nitrogen in reef-building symbiotic corals: A new method and proxy evaluation at Bermuda. Geochimica et Cosmochimica Acta. 2015;148: 179–190.
- 48. Sims ZC, Cohen AL, Luu VH, Wang XT, Sigman DM. Uptake of groundwater nitrogen by a near-shore coral reef community on Bermuda. Coral Reefs. 2020;39: 215–228.
- 49. Jickells TD, Knap AH, Smith SR. Trace metal and nutrient fluxes through the Bermuda inshore waters. Rapports et procès-verbaux des réunions—Conseil international pour l’exploration de la mer. 1986;186: 251–262.
- 50.
Smith S. Bermuda. CARICOMP—Carribean Coral Reef, Seagrass and Mangrove Sites (edited by B Kjerfve). UNESCO; 1998. pp. 247–257.
- 51.
Helmle K, Kohler K, Dodge R. Relative optical densitometry and the Coral X-radiograph densitometry system: Coral XDS. Int Soc Reef Studies 2002 European Meeting, Cambridge. 2002; www.nova.edu/ocean/coralxds/index.html.
- 52. DeCarlo TM, Cohen AL, Barkley HC, Cobban Q, Young C, Shamberger KE, et al. Coral macrobioerosion is accelerated by ocean acidification and nutrients. Geology. 2015;43: 7–10.
- 53. Rippe JP, Baumann JH, Leener DND, Aichelman HE, Friedlander EB, Davies SW, et al. Corals sustain growth but not skeletal density across the Florida Keys Reef Tract despite ongoing warming. Global Change Biology. 2018;24: 5205–5217. pmid:30102827
- 54. Groves SH, Holstein DM, Enochs IC, Kolodzeij G, Manzello DP, Brandt ME, et al. Growth rates of Porites astreoides and Orbicella franksi in mesophotic habitats surrounding St. Thomas, US Virgin Islands. Coral Reefs. 2018;37: 345–354.
- 55. Baumann JH, Ries JB, Rippe JP, Courtney TA, Aichelman HE, Westfield I, et al. Nearshore coral growth declining on the Mesoamerican Barrier Reef System. Global Change Biology. 2019;25: 3932–3945. pmid:31456305
- 56.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2020; www.R-project.org/.
- 57. Buddemeier R, Kinzie R. Coral growth. Oceanography and Marine Biology: An Annual Review. 1976;14: 183–225.
- 58.
Signorell A et. mult. al. DescTools: Tools for descriptive statistics. R package version 0.99.34. https://cran.r-project.org/package=DescTools. 2020.
- 59. Becker M, Nieminen TM, Gérémia F. Short-term variations and long-term changes in oak productivity in northeastern France. The role of climate and atmospheric CO2. Annals of Forest Science. 1994;51: 477–492.
- 60. Mérian P, Lebourgeois F. Size-mediated climate—growth relationships in temperate forests: A multi-species analysis. Forest Ecology and Management. 2011;261: 1382–1391.
- 61.
Bunn A, Korpela M, Biondi F, Campelo F, Mérian P, Qeadon F, et al. dplR: Dendrochronology Program Library in R. R package version 1.7.0. https://CRAN.R-project.org/package=dplR. 2019.
- 62.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1–139. https://CRAN.R-project.org/package=nlme. 2019.
- 63. Castillo KD, Ries JB, Weiss JM. Declining Coral Skeletal Extension for Forereef Colonies of Siderastrea siderea on the Mesoamerican Barrier Reef System, Southern Belize. PLOS ONE. 2011;6: e14615. pmid:21359203
- 64.
Lenth R. emmeans: Estimated Marginal Means, aka Least-Squares Means. R package version 1.4.5. https://CRAN.R-project.org/package=emmeans. 2020.
- 65. Lewis JB. Experimental tests of suspension feeding in Atlantic reef corals. Marine Biology. 1976;36: 147–150.
- 66. Carricart-Ganivet JP. Coral skeletal extension rate: An environmental signal or a subject to inaccuracies? Journal of Experimental Marine Biology and Ecology. 2011;405: 73–79.
- 67. Rico-Esenaro SD, Sanchez-Cabeza J-A, Carricart-Ganivet JP, Montagna P, Ruiz-Fernández AC. Uncertainty and variability of extension rate, density and calcification rate of a hermatypic coral (Orbicella faveolata). Science of The Total Environment. 2019;650: 1576–1581. pmid:30308843
- 68. Cook CB, Logan A, Ward J, Luckhurst B, Berg CJ. Elevated temperatures and bleaching on a high latitude coral reef: the 1988 Bermuda event. Coral Reefs. 1990;9: 45–49.
- 69. Carilli JE, Norris RD, Black BA, Walsh SM, McField M. Local Stressors Reduce Coral Resilience to Bleaching. PLOS ONE. 2009;4: e6324. pmid:19623250
- 70. Cantin NE, Lough JM. Surviving Coral Bleaching Events: Porites Growth Anomalies on the Great Barrier Reef. PLOS ONE. 2014;9: e88720. pmid:24586377
- 71. Carilli JE, Norris RD, Black B, Walsh SM, Mcfield M. Century-scale records of coral growth rates indicate that local stressors reduce coral thermal tolerance threshold. Global Change Biology. 2010;16: 1247–1257.
- 72. Mollica NR, Cohen AL, Alpert AE, Barkley HC, Brainard RE, Carilli JE, et al. Skeletal records of bleaching reveal different thermal thresholds of Pacific coral reef assemblages. Coral Reefs. 2019;38: 743–757.
- 73. Barkley HC, Cohen AL. Skeletal records of community-level bleaching in Porites corals from Palau. Coral Reefs. 2016;35: 1407–1417.
- 74. Smithers SG, Woodroffe CD. Coral microatolls and 20th century sea level in the eastern Indian Ocean. Earth and Planetary Science Letters. 2001;191: 173–184.
- 75. Hendy EJ, Lough JM, Gagan MK. Historical mortality in massive Porites from the central Great Barrier Reef, Australia: evidence for past environmental stress? Coral Reefs. 2003;22: 207–215.
- 76. Cantin NE, Cohen AL, Karnauskas KB, Tarrant AM, McCorkle DC. Ocean Warming Slows Coral Growth in the Central Red Sea. Science. 2010;329: 322–325. pmid:20647466
- 77. Mallela J, Hetzinger S, Halfar J. Thermal stress markers in Colpophyllia natans provide an archive of site-specific bleaching events. Coral Reefs. 2016;35: 181–186.
- 78. Barkley HC, Cohen AL, Mollica NR, Brainard RE, Rivera HE, DeCarlo TM, et al. Repeat bleaching of a central Pacific coral reef over the past six decades (1960–2016). Communications Biology. 2018;1: 1–10.
- 79. Carilli J, Donner SD, Hartmann AC. Historical Temperature Variability Affects Coral Response to Heat Stress. PLOS ONE. 2012;7: e34418. pmid:22479626
- 80. Gardner TA, Cote IM, Gill JA, Grant A, Watkinson AR, Côté IM, et al. Long-Term Region-Wide Declines in Caribbean Corals. Science. 2003;301: 958–960. pmid:12869698
- 81. Edmunds PJ. Decadal-scale changes in the community structure of coral reefs of St. John, US Virgin Islands. Marine Ecology Progress Series. 2013;489: 107–123.
- 82. Green DH, Edmunds PJ, Carpenter RC. Increasing relative abundance of Porites astreoides on Caribbean reefs mediated by an overall decline in coral cover. Marine Ecology Progress Series. 2008;359: 1–10.
- 83. Pandolfi J. Coral community dynamics at multiple scales. Coral Reefs. 2002;21: 13–23.
- 84. Moore JAY, Bellchambers LM, Depczynski MR, Evans RD, Evans SN, Field SN, et al. Unprecedented Mass Bleaching and Loss of Coral across 12°of Latitude in Western Australia in 2010–11. PLoS ONE. 2012;7. pmid:23284773
- 85. Silbiger NJ, Goodbody-Gringley G, Bruno JF, Putnam HM. Comparative thermal performance of the reef-building coral Orbicella franksi at its latitudinal range limits. Marine Biology. 2019;166: 126.
- 86. D’Alelio D, Rampone S, Cusano LM, Morfino V, Russo L, Sanseverino N, et al. Machine learning identifies a strong association between warming and reduced primary productivity in an oligotrophic ocean gyre. Scientific Reports. 2020;10: 3287. pmid:32098970
- 87. Pinto JG, Raible CC. Past and recent changes in the North Atlantic oscillation. WIREs Climate Change. 2012;3: 79–90.
- 88.
Gillett NP, Graf HF, Osborn TJ. Climate Change and the North Atlantic Oscillation. The North Atlantic Oscillation: Climatic Significance and Environmental Impact. American Geophysical Union (AGU); 2013. pp. 193–209.
- 89. Grottoli AG, Rodrigues LJ, Palardy JE. Heterotrophic plasticity and resilience in bleached corals. Nature. 2006;440: 1186–1189. pmid:16641995
- 90. Hughes AD, Grottoli AG. Heterotrophic Compensation: A Possible Mechanism for Resilience of Coral Reefs to Global Warming or a Sign of Prolonged Stress? PLoS ONE. 2013;8: e81172. pmid:24278392
- 91. Conti-Jerpe IE, Thompson PD, Wong CWM, Oliveira NL, Duprey NN, Moynihan MA, et al. Trophic strategy and bleaching resistance in reef-building corals. Science Advances. 2020;6: eaaz5443. pmid:32300659
- 92. DeCarlo TM, Gajdzik L, Ellis J, Coker DJ, Roberts MB, Hammerman NM, et al. Nutrient-supplying ocean currents modulate coral bleaching susceptibility. Science Advances. 2020;6: eabc5493. pmid:32937375
- 93. Van Hooidonk R, Maynard J, Tamelander J, Gove J, Ahmadia G, Raymundo L, et al. Local-scale projections of coral reef futures and implications of the Paris Agreement. Scientific Reports. 2016;6: 39666. pmid:28000782
- 94. Barnes DJ, Lough JM. Systematic variations in the depth of skeleton occupied by coral tissue in massive colonies of Porites from the Great barrier reef. Journal of Experimental Marine Biology and Ecology. 1992;159: 113–128.