Impacts of multiple environmental stressors on coral reef erosion and secondary accretion Nyssa J Silbiger, Òscar Guadayol, Florence I.M. Thomas, Megan J Donahue
Ocean acidification threatens to shift coral reefs from net accreting to net eroding. While corals build reefs through accretion of calcium carbonate (CaCO3) skeletons, net reef
PrePrints
growth depends on bioerosion by grazers and borers and on secondary calcification by crustose coralline algae and other calcifying invertebrates. Primary calcification, secondary calcification, and erosion processes respond differently to climate change stressors; therefore, the combined accretion-erosion response is uncertain. Using a new microcomputed tomography (µCT) method, we measured the simultaneous response of secondary accretion and bioerosion along an environmental gradient: bioerosion rates ranged from 0.02 to 0.91 kg m-2 yr-1 and secondary accretion rates ranged from 0.01 to 0.4 mm yr-1 across a 32m transect. Co-located measures of secondary accretion and bioerosion had different environmental drivers: bioerosion rates were highly sensitive to ocean acidity while secondary accretion rates were most sensitive to physical drivers. These results suggest that bioerosion plays a significant role in the shift from net accretion to net erosion on coral reefs.
PeerJ PrePrints | https://dx.doi.org/10.7287/peerj.preprints.1190v1 | CC-BY 4.0 Open Access | rec: 19 Jun 2015, publ: 19 Jun 2015
1
Title: Impacts of multiple environmental stressors on coral reef erosion and secondary accretion
2
3
Authors: Nyssa J Silbiger1 *, Òscar Guadayol2 , Florence I.M. Thomas1 and Megan J. Donahue1
4
5
6
7
Author affiliations: 1
Hawai‘i Institute of Marine Biology, University of Hawai‘i at M¯anoa, K¯ane‘ohe, Hawai‘i, USA
2
School of Life Sciences, Joseph Banks Laboratories, University of Lincoln, Green Lane, Lincoln, UK
8
PrePrints
9
Corresponding author:
10
Nyssa J Silbiger
11
Hawai‘i Institute of Marine Biology
12
University of Hawai‘i at M¯anoa
13
PO Box 1346
14
K¯ane‘ohe, Hawai‘i, USA 96744
15
Ph: (808) 236-7424
16
email:
[email protected]
17
18
KEYWORDS: Ocean acidification, accretion, bioerosion, natural gradients, coral reefs
19
20
ABSTRACT
21
Ocean acidification threatens to shift coral reefs from net accreting to net eroding. While corals build reefs
22
through accretion of calcium carbonate (CaCO3 ) skeletons, net reef growth depends on bioerosion by graz-
23
ers and borers and on secondary calcification by crustose coralline algae and other calcifying invertebrates.
24
Primary calcification, secondary calcification, and erosion processes respond differently to climate change
25
stressors; therefore, the combined accretion-erosion response is uncertain. Using a new micro-computed
26
tomography (µCT) method, we measured the simultaneous response of secondary accretion and bioerosion
27
along an environmental gradient: bioerosion rates ranged from 0.02 to 0.91 kg m−2 yr−1 and secondary
28
accretion rates ranged from 0.01 to 0.4 mm yr−1 across a 32m transect. Co-located measures of secondary
29
accretion and bioerosion had different environmental drivers: bioerosion rates were highly sensitive to ocean
30
acidity while secondary accretion rates were most sensitive to physical drivers. These results suggest that
31
bioerosion plays a significant role in the shift from net accretion to net erosion on coral reefs.
PeerJ PrePrints | https://dx.doi.org/10.7287/peerj.preprints.1190v1 | CC-BY 4.0 Open Access | rec: 19 Jun 2015, publ: 19 Jun 2015
PrePrints
32
INTRODUCTION
33
Human-induced changes in ocean chemistry (Hoegh-Guldberg et al., 2007; Andersson and Gledhill, 2013; DeCarlo
34
et al., 2014; Silbiger et al., 2014; Silbiger and Donahue, 2015; Wisshak et al., 2012; Tribollet et al., 2009; Reyes-
35
Nivia et al., 2013; Fang et al., 2013), temperature (Hoegh-Guldberg et al., 2007; Silbiger and Donahue, 2015;
36
Davidson et al., 2013; Fang et al., 2013), and water quality (Fabricius, 2005; DeCarlo et al., 2014; Risk et al.,
37
1995b; Holmes et al., 2000; Tribollet and Golubic, 2005; Edinger et al., 2000; Le Grand and Fabricius, 2011) are
38
threatening coral reefs (Pandolfi et al., 2011; Hoegh-Guldberg et al., 2007; Fabricius, 2005). Predictions of reef
39
response to changing ocean conditions are often based on the response of reef building corals alone (Silverman
40
et al., 2009; Pandolfi et al., 2011); however, coral reef bioerosion from borers and grazers and secondary accretion
41
from crustose coralline algae and other encrusting invertebrates are also critical processes for reef sustainability
42
(Przeslawski et al., 2008). There are several gaps in our knowledge about the response of coral reefs to future ocean
43
conditions and whether reef accretion will continue to exceed reef erosion. 1) Will primary accretion, secondary
44
accretion, and reef erosion respond similarly to environmental drivers and will their responses combine to accel-
45
erate reef loss? Studies that examine accretion or erosion processes individually have found different responses
46
to environmental stress. In a field experiment in Indonesia, coral calcification and bioerosion had different func-
47
tional relationships with land-based pollution (Edinger et al., 2000). Laboratory experiments focusing on climate
48
stressors (i.e. temperature and ocean acidity) have found that bioerosion is linearly related to ocean acidity and
49
temperature (Silbiger and Donahue, 2015; Wisshak et al., 2012; Tribollet et al., 2009; Reyes-Nivia et al., 2013;
50
Fang et al., 2013), but that calcification exhibits both linear (Pandolfi et al., 2011; Comeau et al., 2013) and non-
51
linear (Castillo et al., 2014; Silbiger and Donahue, 2015; Comeau et al., 2013; Pandolfi et al., 2011) responses.
52
These differential responses of primary and secondary accretion and bioerosion challenge our ability to predict
53
the net response of coral reefs to environmental change. 2) How will multiple environmental stressors impact
54
individual reef processes? Many environmental parameters interact to drive patterns in accretion and erosion, in-
55
cluding ocean acidity (Hoegh-Guldberg et al., 2007; Andersson and Gledhill, 2013; DeCarlo et al., 2014; Silbiger
56
et al., 2014; Silbiger and Donahue, 2015; Wisshak et al., 2012; Tribollet et al., 2009; Reyes-Nivia et al., 2013;
57
Fang et al., 2013), temperature (Hoegh-Guldberg et al., 2007; Silbiger and Donahue, 2015; Davidson et al., 2013;
58
Fang et al., 2013), nutrients (Fabricius, 2005; DeCarlo et al., 2014; Risk et al., 1995b; Holmes et al., 2000; Tri-
59
bollet and Golubic, 2005), and gradients of human influence (e.g., chlorophyll, turbidity, sedimentation) (Edinger
60
et al., 2000; Le Grand and Fabricius, 2011). This myriad of drivers complicates the predictions of reef response
61
to climate change. 3) How does local, natural variability contextualize global changes in the accretion-erosion
62
balance? Shallow coral reef ecosystems persist in highly variable physicochemical environments (Guadayol et al.,
63
2014; Hofmann et al., 2011; Rivest and Gouhier, 2015) driven by tidal flushing, photosynthesis and respiration,
64
ground-water inputs, and other benthic biological processes (Shamberger et al., 2014; Yates et al., 2007; Drupp
65
et al., 2011; Massaro et al., 2012; Duarte et al., 2013; Smith et al., 2013). Little is known about how climate change
PeerJ PrePrints | https://dx.doi.org/10.7287/peerj.preprints.1190v1 | CC-BY 4.0 Open Access | rec: 19 Jun 2015, publ: 19 Jun 2015
PrePrints
66
interacts with this naturally variable environment to drive patterns accretion and erosion. Here, we address each
67
of these knowledge gaps by leveraging a natural environmental gradient in K¯ane‘ohe Bay to rank the impact of
68
multiple environmental stressors on simultaneous measures of secondary accretion and reef bioerosion.
69
The physicochemical variability of reefs in K¯ane‘ohe Bay, Hawai‘i has been extensively investigated (Guadayol
70
et al., 2014; Drupp et al., 2011; Lowe et al., 2009). Like many reefs around the globe (Hofmann et al., 2011;
71
Rivest and Gouhier, 2015; Shamberger et al., 2014; Manzello et al., 2008; Crook et al., 2013), K¯ane‘ohe Bay has
72
persistent areas of natural acidification that reach the low open-ocean pH values expected by the end of the 21st
73
century (Silbiger et al., 2014; Guadayol et al., 2014). Using a space-for-time framework, we can leverage this
74
natural variability to understand how reefs will respond to ocean acidification in the context of a naturally variable
75
environment. Our prior work in K¯ane‘ohe Bay demonstrates that net reef erosion (calculated as the percent change
76
in volume of experimental blocks) is driven by natural variation in pH and that reefs could shift from net accretion
77
to net erosion with increasing ocean acidity (Silbiger et al., 2014). Here, we separate and simultaneously measure
78
secondary accretion and bioerosion along the same environmental gradient (Figure 1) and then compare the relative
79
importance of a suite of environmental drivers to each of these processes.
80
In prior studies using experimental substrates, pre- and post-deployment buoyant weights have been used to
81
calculate rates of change in density, mass, or volume (reviewed in Table 1), but this method cannot distinguish
82
between secondary accretion and bioerosion processes. In prior studies using natural substrates, imaging method-
83
ologies in 2-dimensions and, more recently, 3-dimensions (CT and µCT) can be used to separate accretion and
84
bioerosion on slabs or cores of reef (Table 1), but rates are difficult to estimate because the time the substrate
85
became available to bioeroders and secondary calcifiers is unknown. Here, we describe a new analysis using µCT
86
to separate secondary accretion and bioerosion from the same experimental substrate exposed to the same envi-
87
ronmental variation over the same time-scale (Figure 2). Our µCT method also allows for a 3D visualization of
88
the experimental blocks that highlights specific areas of secondary accretion and bioerosion (See Movies S1 and
89
S2 in Supplement). Using µCT, we calculate secondary accretion and bioerosion rates from experimental blocks
90
deployed along a 32 m reef transect (Figure S1). Patterns in carbonate chemistry, nutrients, chlorophyll a, tempera-
91
ture, and depth were characterized along this same transect (Figure 1; described in Silbiger et al. (2014)), and used
92
to compare five specific hypotheses about drivers of the accretion-erosion balance: carbonate chemistry, resource
93
availability, temperature, depth, and hydrodynamics (distance from shore and depth). We use a model selection
94
approach to test which of these drivers has the strongest relationship to secondary accretion and bioerosion.
95
Materials and Methods
96
(a) Experimental Design: Our study site is located in K¯ane‘ohe Bay, O‘ahu on the windward (eastern) side of
97
Moku o Lo‘e. We used the same environmental gradient as Silbiger et al. (2014), which we briefly describe here
PeerJ PrePrints | https://dx.doi.org/10.7287/peerj.preprints.1190v1 | CC-BY 4.0 Open Access | rec: 19 Jun 2015, publ: 19 Jun 2015
PrePrints
98
(see supplement for full description of methods). Twenty-one experimental blocks were deployed along a 32 m
99
transect, stratified between reef flat and reef slope (Figure S1). Experimental blocks were cut from dead pieces of
100
massive Porites sp. skeleton into 5cm x 5cm x 2cm blocks, soaked in freshwater, and then autoclaved to remove
101
any living organisms. Collections were made under special activity permit # SAP2011-1 to the Hawai‘i Institute
102
of Marine Biology at the University of Hawai‘i at M¯anoa. The average skeletal density of each experimental
103
block was 1.57 ±0.07(sd) g cm−3 . Blocks were deployed from March 31, 2011 to April 10, 2012. We collected
104
both discrete water samples (pH, TA, nitrate, nitrite, ammonium, phosphate, and chlorophyll a) and data from
105
continuous sensors (temperature and depth) along the transect. Water samples were collected directly above each
106
block four times within 24 hours in September, December, and April in order to capture both diel and seasonal
107
variability in the environment. Continuous sensors were stationed over each block for a minimum of two weeks to
108
calculate high frequency (0.1 min−1 ) variation in temperature and depth. These short time series were normalized
109
to a continuous time series from a permanent station positioned adjacent to the transect, allowing comparison of
110
the micro-environments at each block (Guadayol et al., 2014).
111
(b) µCT: Secondary accretion and bioerosion rates were calculated using µCT (Figure 2). µCT is an X-ray
112
technology that non-destructively images the external and internal structures of solid objects, resulting in a three-
113
dimensional array of object densities. We used an eXplore CT120 µCT (GE Healthcare Xradia, Inc) at the Cornell
114
University Imaging Multiscale CT Facility (Figure 2.1-2) to scan blocks before and after deployment (voltage =
115
100kV, current = 50mA). A three-dimensional array of isotropic voxels at 50 µm3 resolution was generated using
116
the GE Console Software and were averaged to 100 µm3 for data analysis. We used a threshold of 200 Hounsfield
117
Units to separate coral from air (Silbiger et al., 2014) (Figure 2.3). The pre and post-deployment scans were aligned
118
using an intensity-based registration technique (Figure 2.4), converted to binary (Figure 2.5), and subtracted from
119
one another resulting in a matrix of 0’s, 1’s, and -1’s (Figure 2.5). All positive values were new pixels added to the
120
post-deployment scan which indicate secondary accretion, negative values were pixels that were lost and indicate
121
bioerosion, and zeros meant there was no change at that pixel between the two scans. All values were summed and
122
multiplied by the resolution of the scan to obtain the volume lost or gained per block (Figure 2.6).
123
(c) Rates: Bioerosion and secondary accretion rates were calculated using the following equations: Bioerosion
124
Rate (kg m−2 yr−1 ) = (V oli × ρi )/(SAi × T ime) and Secondary Accretion Rate (mm yr−1 ) = V oli /(SAi ×
125
T ime) , where i represents an individual block, V ol is the volume lost (bioerosion) or gained (secondary accretion)
126
in m3 , SA is the surface area of the pre-deployment blocks (m2 ), ρ is the skeletal density of the pre-deployment
127
block (kg m−3 ), and T ime is the deployment time (years). Secondary accretion rates were converted from m to
128
mm per year to compare with literature values.
129
(d) Model Selection: Our goal was to compare pH with other known drivers and correlates of the accretion-
130
erosion balance. Many of the environmental variables were collinear along the transect; thus, we removed collinear-
131
ity by using the residuals of a regression of each environmental variable against log(depth) and distance from shore
PeerJ PrePrints | https://dx.doi.org/10.7287/peerj.preprints.1190v1 | CC-BY 4.0 Open Access | rec: 19 Jun 2015, publ: 19 Jun 2015
PrePrints
132
(Silbiger et al., 2014). Correlation coefficients for raw environmental data and the residual environmental data are
133
available in Silbiger et al. (2014).
134
We used a model selection framework to compare models for five specific hypotheses about the accretion-
135
erosion balance and test which of these drivers had the strongest relationship to secondary accretion and bioerosion
136
(Table 2). In a model selection framework, Akaike Information Criterion (AIC) values are used to rank candidate
137
models, accounting for both fit and complexity. Carefully constructed model selection avoids problems associated
138
with multiple hypothesis testing that are common in stepwise regression, such as arbitrary α levels and uninter-
139
pretable functional relationships (Johnson and Omland, 2004; Anderson, 2008). Here, we used the corrected AIC
140
(AICc), which is recommended for sample sizes