Subscriber access provided by UNIV DI BOLOGNA
Article
Time-Domain NMR study of Mediterranean scleractinian corals reveals skeletal-porosity sensitivity to environmental changes Paola Fantazzini, Stefano Mengoli, Stefania Evangelisti, Luca Pasquini, Manuel Mariani, Leonardo Brizi, Stefano Goffredo, Erik Caroselli, Fiorella Prada, Giuseppe Falini, Oren Levy, and Zvy Dubinsky Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 21 Oct 2013 Downloaded from http://pubs.acs.org on October 29, 2013
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 27
Environmental Science & Technology
TOC 47x26mm (300 x 300 DPI)
ACS Paragon Plus Environment
Environmental Science & Technology
1 2 3 4
Page 2 of 27
Time-Domain NMR study of Mediterranean scleractinian corals reveals skeletal-porosity sensitivity to environmental changes Paola Fantazzinia,b*, Stefano Mengolic, Stefania Evangelistia, Luca Pasquinia, Manuel Mariania,b,
5
Leonardo Brizia,b, Stefano Goffredod, Erik Carosellid, Fiorella Pradad, Giuseppe Falinie, Oren Levyf,
6
Zvy Dubinskyf
7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
a
Department of Physics and Astronomy, University of Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy, b Centro Enrico Fermi, Roma, Italy, c Management Department, University of Bologna, Bologna, Italy, d Marine Science Group, Department of Biological, Geological and Environmental Sciences, University of Bologna, Bologna, Italy, e Department of Chemistry, University of Bologna, Bologna, Italy, f The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan, Israel. *
Corresponding author: Paola Fantazzini; tel: +390512095119; fax: +390512095047;
[email protected]
22
1 ACS Paragon Plus Environment
Page 3 of 27
Environmental Science & Technology
23 24 25 26 27
Abstract
28
Mediterranean corals are a natural model for studying global warming, as the Mediterranean
29
basin is expected to be one of the most affected regions and the increase of temperature is
30
one of the greatest threats for coral survival. We have analyzed for the first time with Time
31
Domain Nuclear Magnetic Resonance (TD-NMR) porosity and pore-space structure,
32
important aspects of coral skeletons, of two scleractinian corals, Balanophyllia europaea
33
(zooxanthellate) and Leptopsammia pruvoti (non-zooxanthellate), taken from three different
34
sites on the western Italian coast along a temperature gradient. Comparisons have been made
35
with Mercury Intrusion Porosimetry and SEM images. TD-NMR parameters are sensitive to
36
changes of the pore-structure of the two coral species. A parameter, related to the porosity,
37
is larger for Leptopsammia pruvoti than for Balanophyllia europaea, confirming previous
38
non-NMR results. Another parameter representing the fraction of the pore-volume with
39
pore-sizes less than 10-20 µm is inversely related, with high statistical significance, to the
40
mass of the specimen, and, for Balanophyllia europaea, to the temperature of the growing
41
site. This effect in the zooxanthellate species, that could reduce its resistance to mechanical
42
stresses, may depend on an inhibition of the photosynthetic process at elevated temperatures
43
and could have particular consequences in determining the effects of global warming on
44
these species.
45 46 47 48 2 ACS Paragon Plus Environment
Environmental Science & Technology
49
Page 4 of 27
INTRODUCTION
50
Corals and global warming. Global climate change is the defining environmental issue of
51
our times and is expected to profoundly affect all levels of ecological hierarchies and a broad array
52
of terrestrial and marine ecosystems.1-5 Marine communities are expected to be affected more than
53
terrestrial ones by the effects of climate change,6 especially in temperate areas.7 Thus, the
54
Mediterranean basin8 represents a natural focus of interest for researchers and at the same time a
55
natural laboratory to model and predict climate change and its ecological effects. In particular, the
56
increase of temperature is one of the greatest threats for corals, which can be considered as a probe
57
of global warming effects, as it triggers bleaching events and widespread mortality.9-11 Several
58
recent mass mortality events of Mediterranean corals have been reported as related to high
59
temperatures.12-16
60
This study focuses on two scleractinian species of the Mediterranean Sea, already studied as a
61
model for climate change: Balanophyllia europaea (Risso, 1826) (B. europaea) and Leptopsammia
62
pruvoti Lacaze-Duthiers (1897) (L. pruvoti)17 (Fig. 1). B. europaea is a solitary, zooxanthellate (i.e.
63
symbiotic with unicellular algae named zooxanthellae) coral, endemic to the Mediterranean Sea. Its
64
distribution is limited to depths of 0-50 m because of its symbiosis with zooxanthellae, which
65
require light.17-19 L. pruvoti is a non-zooxanthellate and solitary scleractinian coral, distributed in the
66
Mediterranean basin and along the European Atlantic coast from Portugal to Southern England and
67
Ireland. Its distribution is limited to semi-enclosed rocky habitats, under overhangs, in caverns, and
68
small crevices at 0–70 m depth.17,18 Corals were collected in three different Italian sites (see Fig. 2),
69
where the porosity of the two species has been studied previously,17along a latitude and sea surface
70
temperature (SST) gradient. Temperature, the variation of which is mainly influenced by latitude,20
71
is linked to coral biometry, physiology, and demography.21-23
72
Coral porosity and pore-size distribution. The porosity (pore-volume to sample-volume
73
ratio) and the pore-space structure of mineralized tissues are of crucial importance in determining 3 ACS Paragon Plus Environment
Page 5 of 27
Environmental Science & Technology
74
overall properties and biological functions, such as the coral skeleton resistance to natural and
75
anthropogenic breakage. They are important parameters for studying scleractinian coral growth and
76
the effects of abiotic and anthropogenic influences on coral reefs.24 Of additional interest is a good
77
knowledge of the role of porosity in diagenesis.
78
The measured values of these parameters depend on the measurement methods24,25 and can
79
present spatial variations.26 This could be the reason why their dependence on environmental
80
conditions remains largely unstudied, notwithstanding their importance and their variation with
81
factors such as exposure, temperature, latitude depth and species.
82
Recent investigations17,19,27 have shown that along the Italian coast porosity of B. europaea is
83
positively correlated with SST. There is concern for the future of this species17,27,28 in relation to the
84
current predictions of global warming by the Intergovernmental Panel on Climate Change.7 Instead,
85
for L. pruvoti, both SST and solar radiation do not seem to influence significantly the porosity and
86
the space colonization potential.17,19,29 An important aspect is the possible hierarchical structure of
87
the porosity. Recently30 the flaw-tolerance in nacre has been ascribed to the nanoparticle-
88
architecture of the aragonite platelet, which makes a crack propagate in an intergranular manner.
89
The structure of porous media at different length scales is of great importance also in the use of
90
corals as potential bone graft substitute material.31
91
Porosity and pore-size distributions can be investigated by many methods.32-36 The results
92
strongly depend on the physical principles adopted and on the assumptions on pore shape and
93
connectivity (see Supporting Information). Porosity and pore size distributions by Mercury
94
Intrusion Porosimetry (MIP) on eight different coral species36 showed large differences, with
95
diameters ranging from 0.2 to 100 µm. However, it has been emphasized that particle compression
96
and rupture can result from the high Hg pressure used. Time-Domain Nuclear Magnetic Resonance
97
(TD-NMR)37 has the advantage of being non-destructive and non-invasive.
98
particular Magnetic Resonance Relaxometry of 1H nuclei of water saturating the pore-space is an
99
efficient tool to investigate pore-space structure. Known from the ’50s,38,39 and validated in the 4 ACS Paragon Plus Environment
TD-NMR and in
Environmental Science & Technology
Page 6 of 27
100
course of time by comparison with MIP and, for specific surface, with the Brunauer, Emmett, Teller
101
(BET) method, it is now widely applied.40-52 It is particularly useful for porous media with wide
102
pore-sizes distributions, as corals have. In this paper the distribution of the local transverse
103
relaxation time (T2) of 1H of water saturating the pore-space of the cleaned coral skeletons,
104
corresponding to distributions of “NMR pore-sizes”, are used (more details on NMR and surface
105
effects in Supporting Information). To the best of our knowledge, this is the first time that this
106
technique is applied for this kind of investigation of corals.
107 108
MATERIALS AND METHODS
109
The corals. Specimens of B. europaea and L. pruvoti (54 specimens in total) were randomly
110
collected from three sites: Calafuria (CL), Palinuro (PL) and Pantelleria Isle (PN) (see Fig. 2). Coral
111
tissue was totally removed and corals cleaned as described in Ref. 17. The skeletons were weighed
112
to determine the mass (m). The total volume (VT) was determined19 including the volume of the oral
113
cavity. Then the specimens were saturated with water for NMR measurements. Further details are
114
reported in Supporting Information.
115
Total NMR Signal, microporosity and cut-off definitions. The total NMR signal (SNMR),
116
represented by the area below each T2 distribution, is proportional to the volume of water saturating
117
the pore-space volume (VP). This signal divided by the total sample volume gives a value
118
proportional to the total porosity of the specimen (see Supporting Information). The fraction of
119
water with relaxation times over a given interval of the distribution corresponds to the pore volume
120
fraction over a corresponding pore-size range. “NMR microporosity", “microporosity" for short,
121
will indicate the fraction of VP where the smaller pores are weakly coupled by water diffusion to the
122
large ones at the local relaxation time scale. This can be accomplished if the slope of the
123
distribution shows a strong increase at a certain T2 value, to be chosen as the point of separation
124
between “smaller” and “larger” pores. This relaxation time will be called the “cut-off”. The 5 ACS Paragon Plus Environment
Page 7 of 27
Environmental Science & Technology
125
“microporosity” is then defined as the fraction of 1H signal with T2 smaller than the cut-off, divided
126
by the total 1H signal. Operatively, it is the ratio of the area under the distribution for T2 smaller
127
than the cut-off, to the total area under the distribution.
128
Statistical Analysis. Statistical analysis was performed using Statistical Package STATA 9.0
129
(StataCorp LP). To test the significance of the differences among species and growth site
130
parametric and non-parametric tests were performed. Multivariate analyses were made using both
131
Ordinary Least Squares (OLS) robust to outliers and non-parametric bootstrapping regression
132
procedure following Efron,53 applied to check the robustness of the results, that could be affected by
133
small sample bias. The models are described by the function:
134
yi = a + b1 ⋅ mi + b2 ⋅ SSTi + ε i (1),
135
the index i refers to the n-observations, yi is the value of the dependent variable and εi the
136
corresponding error. The parameters a and bj (j=1,2) are the best fit parameters, to be determined by
137
OLS referring to the independent variables m and SST. More details in Supporting Information.
138 139
RESULTS AND DISCUSSION
140
Figure 3 shows the T2 distributions of all B. europaea and L. pruvoti specimens. All the
141
distributions show a main peak at long relaxation times and a long tail, of smaller amplitude, about
142
three order of magnitudes wide. The major shape difference between the two species is the length of
143
the tail. For B. europaea (a, c, e) the tails go down to T2 = (0.1 - 0.2) ms, values shorter than for L.
144
pruvoti (b, d, f). The principal differences among each group are given by the total areas and are due
145
to the wide range of masses and volumes of the specimens. In principle, 1H signal can be produced
146
from other sources than water, namely the intraskeletal organic matrix, consisting of proteins,
147
polysaccharides and lipids. In order to check this possibility, T2 distributions for a dry coral and the
148
same after complete water saturation were obtained and are shown in Fig.4. The discussion reported
149
in Supplementary Information (referring to Ref. 54-61) leads us to conclude that in the distribution 6 ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 27
150
for the fully water saturated sample there is no contribution to the signal from macromolecular
151
nuclei and only a maximum on the order of 2% could be attributed to lipids.
152
In order to discuss the distributions in Fig. 3 in terms of pore-sizes, one should get an
153
approximate value for the radius of the pores inside the coral. The inset in Fig. 4 reports the MIP
154
distribution for the same specimen. The major fraction of the pore-volume is given by pores whose
155
entrance radii are on the order of tens of µm, while a minor fraction corresponds to a three order of
156
magnitude long tail of very small pores, down to tens of nm. NMR and MIP are exceptionally
157
similar and consistent. The two classes of pores are easily distinguished also in the distributions in
158
Fig. 3, so that the two parameters “microporosity” and “cut-off” were determined for each
159
distribution. The sharp boundary between the two classes, with a cut-off in the T2 range 200-400
160
ms, suggests that the two classes of pores are not well connected by water diffusion during a local
161
relaxation time. Also, the long tail indicates that these pores are poorly connected both to the other
162
small pores and to the large ones in the major class. On the basis of the comparison with MIP,
163
microporosity should correspond to pore sizes in the range from ∼10 nm to ∼10-20 µm. The
164
existence of a wide class of pores with sizes less than 10-20 µm is well described in the SEM
165
images of both species reported in Fig. 5.
166
Table 1 shows means, standard errors and statistical significances of the differences between
167
the two species by both parametric and non-parametric tests for all the variables considered:
168
microporosity, cut-off, mass, total volume, SNMR, and SNMR/VT. The two species behave differently
169
with high statistical significance (p< 0.01 for all variables, including NMR parameters). In
170
particular mass and total volume for B. europaea are much larger than for L. pruvoti, and,
171
viceversa, total porosity estimated by the ratio SNMR/VT is larger for L. pruvoti than for B. europaea.
172
This result is consistent with the higher porosity obtained for L. pruvoti by previous non-NMR
173
analysis,17 where the difference was considered as a likely consequence of the different habitat of
174
the two species.
7 ACS Paragon Plus Environment
Page 9 of 27
Environmental Science & Technology
175
Table 2 shows the differences of the means of the same species among growing sites. For B.
176
europaea the differences have high statistical significance for almost all variables. Microporosity
177
has the highest significance (microporosity, p < 0.01; cut-off, p < 0.05; mass and total volume, p <
178
0.1; SNMR, p < 0.05 ). For L. pruvoti, only NMR parameters cut-off and microporosity show
179
significant differences among sites (p< 0.05-0.1).
180
Table S1 of Supporting Information shows the results of the correlation among variables
181
performed separately for B. europaea (Panel A) and L. pruvoti (Panel B). For both species, the pairs
182
of variables mass and total volume, mass and SNMR, total volume and SNMR are significantly
183
correlated (p < 0.01). The correlations between microporosity and cut-off, and microporosity and
184
mass are statistically significant for both species (p < 0.01 for B. europaea, p < 0.05 for L. pruvoti).
185
Figure 6 reveals a counter-intuitive behavior: the longer the cut-offs, the smaller the
186
microporosities, while, at the first glance, one would expect the contrary. The pore-space
187
architecture differs between samples with higher or lower cut-offs (between smaller microporosity
188
or higher microporosity). The scatterplots in Figure S2 of Supporting Information and in Fig. 7,
189
showing cut-off versus mass and microporosity versus mass, respectively, suggest that the observed
190
correlation between cut-off and mass are governed by the mass: the smaller the mass, the higher the
191
microporosity and the shorter the cut-off. As the corals increase their mass, both the cut-off (which
192
separates the two main pore classes) and the ratio between the fraction of the two pore classes
193
change, with larger pores becoming more abundant. This effect is shown also by L. pruvoti, but it is
194
not as marked as in B. europaea due to the smaller range of masses of the corals. This is consistent
195
with the gradual “filling up” of the smaller pores with the growth of the coral. A secondary infilling
196
of skeletal pores in the older portion of the skeleton is a consistent characteristic of the skeletal
197
density of branches of tall branching corals, in which growing tips are very porous, while basal
198
regions are extremely dense.62
199
In order to study how microporosity, considered as a dependent variable, is affected by mass
200
and SST, multiple regression analysis was performed (Eq. 1) both for all the specimens and for the 8 ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 27
201
two species separately. Panel A of Table 3 summarizes the results. First of all, it is important to
202
observe that from a statistical point of view a potential correlation among the variables m and SST
203
does not invalidate the results, as all the values of the variance inflation factor (VIF) are 0.1). Overall, it emerges that the mass significantly affects microporosity at least at the 5
217
percent level (p< 0.05 for L. pruvoti and p< 0.01 for B. europaea). Results for microporosity vs
218
SST show that a significant relationship exists for B. europaea (p< 0.05) but not for L. pruvoti (p>
219
0.1).
220
In a previous study17 it has been shown that porosity depends on temperature for B. europaea,
221
but not for L. pruvoti. It has been hypothesized that the increase in porosity with temperature in the
222
zooxanthellate species could depend on an inhibition of the photosynthetic process at elevated
223
temperatures,23,63 causing an attenuation of calcification64 with possible negative consequences also
224
on space colonization and population density.19,27 The NMR results point in the same direction and
225
seem to indicate also that this effect could be accompanied by a reduction of microporosity,
226
meaning an increase of the fraction of the largest pores in the pore-space. 9 ACS Paragon Plus Environment
Page 11 of 27
Environmental Science & Technology
227
TD-NMR is a quick, non-invasive, non-destructive method that does not use ionizing radiation,
228
that can be applied to gain insight regarding the pore-space architecture of scleractinian corals,
229
showing differences between species and growing sites, and sensitivity to environmental changes.
230
Of course, this method, as well as MIP, BET and the hydrostatic balance method, gives information
231
on the connected porosity only and, as such, can be applied in systems with low fractions of isolated
232
pores. This method can give information not attainable in other ways, like changes of the internal
233
architecture of corals described by microporosity and cut-off with increasing mass and growing
234
temperature. Even if this method cannot spatially locate the heterogeneity of the pore space, the
235
existence of a clear cut-off in almost all the distributions (a very high slope at a certain point of the
236
distribution) means that the smallest pores are not well connected by diffusion at the NMR time
237
scale (corresponding to the local value of T2) to the largest ones. Moreover the NMR defined
238
parameter “microporosity” can quantify the ratio between the volume of the smallest pores (sizes
239
less than 10-20 µm) and the total pore-volume.
240
The increased fraction of larger pores in the zooxanthellate corals with increasing SST, that could
241
reduce their resistance to mechanical stresses, could have particular consequences in determining
242
effects of global warming on these species.65 The described method will be applied in future work
243
to the effects of ocean acidification on the skeletal properties of corals66 and other calcifying
244
organisms.
245 246
ACKNOWLEDGEMENTS
247
This research has received funding from the European Research Council under the European
248
Union's Seventh Framework Programme (FP7/2007-2013) / ERC grant agreement n° 249930 –
249
CoralWarm: Corals and global warming: the Mediterranean versus the Red Sea. We also thank F.
250
Gizzi, C. Marchini and S. Prantoni for help in collecting the samples, the diving centers Centro
251
Immersioni Pantelleria and Il Pesciolino for logistic assistance in the field, the Scientific Diving 10 ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 27
252
School for collaboration in the underwater activities, Fausto Peddis for MIP measurements, Gianni
253
Neto for pictures of living specimens and Robert James Sidford Brown for carefully reading the
254
manuscript.
255 256
Supporting Information Available
257
All references numbers in Supporting Information refer to the Reference list in the main text.
258
This information is available free of charge via the Internet at http://pubs.acs.org/.
259 260 261
11 ACS Paragon Plus Environment
Page 13 of 27
Environmental Science & Technology
262
B. europaea
Parameter
microporosity cut-off m VT SNMR SNMR/ VT
268 269 270 271 272 273 274 275
263
L. pruvoti
n
mean
se
n
mean
se
27 27 27 27 27 27
31.3 337 0.88 0.92 4060 4800
1.6 15 0.15 0.17 820 275
26 26 26 27 27 27
38.6 249 0.25 0.20 1233 6993
2.1 20 0.04 0.03 165 349
t
Z
2.79*** 3.59*** 3.92*** 4.13*** 3.38*** 4.94***
264 6.66*** 12.92*** 265 11.19*** 17.46*** 266 11.21*** 17.60*** 267
Table 1. Descriptive and test statistics split by species *** p