Next Article in Journal
An Assessment of Applicability of SNP Chip Developed for Domestic Goats in Genetic Studies of Caucasian Tur (Capra caucasica)
Next Article in Special Issue
Seagrass Meadows Provide a Significant Resource in Support of Avifauna
Previous Article in Journal
Variation in Alpine Plant Diversity and Soil Temperatures in Two Mountain Landscapes of South Patagonia
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 13
Issue 7
10.3390/d13070311
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ocean Acidification and Mollusc Settlement in Posidonia oceanica Meadows: Does the Seagrass Buffer Lower pH Effects at CO2 Vents?

by
Alessandra Barruffo
1,†,
Laura Ciaralli
1,†,
Giandomenico Ardizzone
1,
Maria Cristina Gambi
2,‡ and
Edoardo Casoli
1,*
1
Department of Environmental Biology, University La Sapienza, 00185 Roma, Italy
2
National Institute of Oceanograhy and Applied Geophysics, OGS, 34100 Trieste, Italy
*
Author to whom correspondence should be addressed.
†
These Authors contributed equally to the research.
‡
Prior affiliation: Stazione Zoologica Anton Dohrn, Department Integrative Marine Ecology, Ischia Marine Center, Ischia, 80100 Napoli, Italy.
Diversity 2021, 13(7), 311; https://doi.org/10.3390/d13070311
Submission received: 4 June 2021 / Revised: 3 July 2021 / Accepted: 5 July 2021 / Published: 8 July 2021
(This article belongs to the Special Issue Seagrass Ecosystems, Associated Biodiversity, and Its Management)
Browse Figures
Review Reports Versions Notes

Abstract

:
Ocean acidification has been broadly recognised to have effects on the structure and functioning of marine benthic communities. The selection of tolerant or vulnerable species can also occur during settlement phases, especially for calcifying organisms which are more vulnerable to low pH–high pCO2 conditions. Here, we use three natural CO2 vents (Castello Aragonese north and south sides, and Vullatura, Ischia, Italy) to assess the effect of a decrease of seawater pH on the settlement of Mollusca in Posidonia oceanica meadows, and to test the possible buffering effect provided by the seagrass. Artificial collectors were installed and collected after 33 days, during April–May 2019, in three different microhabitats within the meadow (canopy, bottom/rhizome level, and dead matte without plant cover), following a pH decreasing gradient from an extremely low pH zone (pH < 7.4), to ambient pH conditions (pH = 8.10). A total of 4659 specimens of Mollusca, belonging to 57 different taxa, were collected. The number of taxa was lower in low and extremely low pH conditions. Reduced mollusc assemblages were reported at the acidified stations, where few taxa accounted for a high number of individuals. Multivariate analyses revealed significant differences in mollusc assemblages among pH conditions, microhabitat, and the interaction of these two factors. Acanthocardia echinata, Alvania lineata, Alvania sp. juv, Eatonina fulgida, Hiatella arctica, Mytilys galloprovincialis, Musculus subpictus, Phorcus sp. juv, and Rissoa variabilis were the species mostly found in low and extremely low pH stations, and were all relatively robust to acidified conditions. Samples placed on the dead matte under acidified conditions at the Vullatura vent showed lower diversity and abundances if compared to canopy and bottom/rhizome samples, suggesting a possible buffering role of the Posidonia on mollusc settlement. Our study provides new evidence of shifts in marine benthic communities due to ocean acidification and evidence of how P. oceanica meadows could mitigate its effects on associated biota in light of future climate change.

1. Introduction

In the last two centuries, the atmospheric CO2 concentration levels have steadily increased due to anthropic activities, rising from 280 ppm during the Industrial Revolution to 411.5 ppm today [1]. Oceans act as an essential carbon sink, absorbing part of anthropogenic CO2 emissions; however, this is causing a rise in pCO2 concentrations at the sea surface, leading to a decrease in the ocean pH level and change in carbonate chemistry, a process broadly known as “ocean acidification” (OA) [2]. The Intergovernmental Panel on Climate Change models [3] predict that the pH of surface seawater will drop by 0.4 units by the year 2100, increasing the acidity levels by 170% compared to the pre-industrial ones. OA is an ongoing process expected to intensify in the future, which, together with the other aspects of Global Climate Change, drastically threatens the structure and functioning of marine ecosystems, especially in a semi-closed basin such as the Mediterranean Sea [4].
Changes in pH and carbonate water chemistry have a severe impact at both the individual and ecosystem levels. Several studies and some reviews and meta-analyses highlighted different biological responses to OA across a wide range of taxa with adverse effects on survival, metabolism, calcification, growth, and reproduction ([5,6] and reference therein), especially for calcifiers. The rise of CO2 in ocean waters leads, in fact, to more corrosive conditions, as well as a depletion of calcium carbonate, that affect calcifying organisms hindering them from building and maintain their carbonate shells and skeletons [7,8].
In light of the aforementioned biological responses, the structure and functioning of marine communities will change under acidified conditions: several authors have reported the loss of both species richness and functional diversity in the study of Mediterranean benthic communities associated with CO2 vents systems [9,10,11,12,13], as well as alteration of competitive interactions and an overall trophic simplification [11,14]. These vent systems, characterised by zones of intense CO2 emissions from the seafloor, represent functional natural laboratories to study the effect of medium- and long-term exposure to naturally acidified conditions on the benthic biota, and up to date, several of such natural systems have been described and studied around the world [15,16,17]. From vent studies, changes in benthic community structure due to OA may also depend on the larval pelagic stages or settlement processes that characterise invertebrate recruitment [18,19].
Posidonia oceanica (L.) Delile is the endemic seagrass species in the Mediterranean Sea and the most widespread; it forms extensive meadows in the infralittoral zone, representing one of the most important coastal ecosystems, where it plays key ecological, physical, and economic roles [20,21,22]. Posidonia oceanica meadows host several benthic assemblages within the multitude of microhabitats and the stratified layers (i.e., leaf canopy, rhizomes, and sediment/matte) they provide [23], forming a complex food web [24]. They also represent an important nursery and shelter for many animals and seaweeds, such that Boudouresque and Meinesz [25] named it an “oasis of life”. Posidonia oceanica meadows are threatened all around the Mediterranean Sea by several anthropogenic impacts [24,26] and climate change, since the species is susceptible to the increase of the temperature and has been often impacted by summer heat waves [27,28].
On the other hand, P. oceanica meadows have a pivotal role in mitigating climate change, mainly through carbon sequestration, which is stored in the sediment and matte over millennia (around 27% of the carbon fixed by the plant) [26,29,30,31]. If considered at the scale of the whole Mediterranean, P. oceanica meadows adsorb 11–42% of the CO2 emissions produced around the coasts of this basin [26]. A further aspect of climate mitigation is related to Posidonia and other seagrasses’ photosynthetic performance, which buffers the low pH and reduces the effect of OA on the local benthic communities [32,33,34]. A similar buffering effect has also been demonstrated for some large seaweeds, also in respect to protecting larval stages [35,36].
Mollusca were recognised as one of the most well-known and diverse taxonomic groups in Mediterranean Sea habitats, encompassing around 2000 species [37]: their diversity and abundance are relevant in peculiar and protected Mediterranean habitats [38]. Mollusc settlement, as well as that of other marine invertebrates, is a complex process determined by the interaction of biotic and abiotic factors. Settlement patterns are influenced by the environmental features and habitat type, hydrodynamic process, seasonality, proximity of adult populations, chemical cues, and presence of conspecific/prey or other organisms (e.g., filamentous algae) [39,40]. It is difficult to distinguish among the contribution of each of the factors mentioned above, due to the different spatial and temporal scales at which they act, and because settlement patterns are often inferred from recruitment. The abundance and diversity of molluscs are particularly relevant in P. oceanica beds, where, due to the habitat complexity, they show a wide range of forms, functions, and adaptations. For these reasons, they have been widely investigated along Italian coasts [41,42,43,44,45,46,47,48,49,50]. Molluscs have also been studied comparing two well-distinguished habitats within a P. oceanica meadow: leaf canopy and rhizomes. These two layers present differences in light intensity, hydrodynamics, stability and dynamics, and habitat heterogeneity. These two layers host different mollusc assemblages and zonation patterns within the meadow [51,52].
Experiments and studies carried out at vents of the Castello at Ischia, and in other vent systems, contributed to recognising molluscs as particularly sensitive to OA, especially on the vegetated hard bottoms [11,12,53]. In contrast, they seem less impacted within the Posidonia meadows ([15,51] and reference therein). The decrease in pH affects both larvae and adults’ calcium carbonate secretion to build shells; furthermore, OA threatens metabolic patterns, acid-base status, reproduction, immune response, and the survival of molluscs. Given the diversity and abundance of molluscs, their altered conditions are supposed to have several ecological and economic consequences for the structure and functioning of the entire P. oceanica meadow ecosystem.
This study aims to assess the effect of the pH variation, along a natural gradient of OA in three CO2 vent’s systems, on the settlement of invertebrate fauna associated with P. oceanica, focusing on early-settled stages of Mollusca. This group has been selected due to its high diversity in Posidonia meadows and its potential sensitivity to OA, since most species associated with seagrasses are calcareous. Furthermore, in one of these systems, we tested the hypothesis that there would be significant differences in mollusc’ settlement among stations characterised by living P. oceanica and those with dead matte bottoms, to assess the possible buffering effect provided by the seagrass.

2. Materials and Methods

2.1. Study Area

Among the vent’s system, the Castello Aragonese off the island of Ischia (Tyrrhenian Sea, Italy) represents model conditions to investigate how benthic invertebrates respond to OA in complex coastal systems, such as seagrass meadows. Another shallow system named “Vullatura”, located near the Castello and recently investigated [13,54], is characterised by extended shallow P. oceanica meadows. These sites have been widely used, since they mimic future environmental conditions with a decrease in pH, even though they cannot be considered as perfect proxies of OA at a larger scale [15]. The study was carried out in three vent sites: two located in proximity to the Castello Aragonese (40°43′53.53″ N, 13°57′50.98″ E, north-eastern coast of Ischia, Gulf of Naples, Italy), on the north and south sides of the Castello (N and S), and a third site, named “Vullatura” (Vu) [55], located approx. 500 m north to the Castello (40°44′9.53″ N; 13°57′28.19″ E) (Figure 1). The gas emissions at the three vent systems are almost entirely composed of CO2 (90.1–95.3%), and no differences in temperature and salinity occur if compared to surrounding waters, as well as no toxic sulphur compounds being associated with the emissions [9,15,55]. At the Castello, the vents cover about 5000 m2, with differences in bubble emissions/m2 among N and S sides: 0.7 × 106 L of gas emitted/day in the former and 1.4 × 106 L of gas emitted/day in the latter [9]. At the Vullatura, bubble emissions cover approximately 1200 m2 with high densities of bubbles/m2, as suggested by the site name, which means “boiling water” in the local vernacular [56]. These CO2 flows modify pH and alkalinity of seawater at each of the sites mentioned above, creating a pH gradient developing parallel to the coastline.
Dense P. oceanica beds develop between 1 and 6 m depth at Castello and Vullatura, with the highest shoot densities (up to 1000 shoots/m2) reported under the most acidified conditions at the Castello [5,9,57,58].
In all of the vent sites here considered, P. oceanica meadows settle on the sandy bottom. On the north and south sides of Castello, Posidonia forms a matte approx. 60 cm high, while at the Vullatura there are high pinnacles and ridges of matte up to 3 m high [13,56]. The shoot density at all stations have been measured in various previous studies, and show the highest values in the acidified stations both on the north and south sides of the Castello, and at the Vullatura (>1000 shoots m2 at S3, 866 ± 148 shoots/m2 at Vu3, and 858 ± 85 shoots m2 at N3), while they are lower at the control stations (438 ± 88 shoots m2 at N1; 494 ± 125 shoots m2 at SC) [13,58].
At each of the sites, three stations characterised by different pH conditions were selected: control stations with ambient pH (8.2–8.0 pH units; N1, SC, and Vu1), low pH stations (7.8–7.4 pH units; N2, S2, and Vu2), and extremely low pH stations (less than 7.4 pH units; N3, S3, and Vu3) (Figure 1). The pH variations among and within the studied stations were assessed in previous studies, and here are reported graphically modified data available in the literature [13,53]. The pH data were measured in situ through the hourly acquisition of Sea-FET pH-meters, which was performed over more than 23 days at each of the investigated stations [59]. The SeaFET pH-meters were arranged close to the P. oceanica meadows: due to the size of the instrument and its mooring, and the canopy features as well, it was not possible to place them directly into the Posidonia canopy; therefore, the pH measurements refer to the meadow’s surrounding waters.

2.2. Sampling and Data Collection

Artificial collectors are represented by 8 cm diameter rounded non-toxic nylon scouring pads; they were placed underwater (by scientific SCUBA divers) at each of the stations and collected after 33 days during the spring season of 2019 (from 28 April to 30 May). The same type of artificial collectors have been largely used in previous studies carried out at the Castello Aragonese vent systems [18,19,60], because they represent an effective tool for the collection of juvenile and adult invertebrate specimens. One month-long sampling period was chosen to allow sufficient time for invertebrate settlement. The April–May period was selected because it represents an active reproduction and settlement phase for several invertebrate species. At each vent and station, 4 scouring pad replicates were fixed to 1 m long bars and placed above the seagrass canopy (C), and an additional 4 replicates were fixed to 20 cm long bars and placed close to the bottom (B), at rhizome levels. Thus, a total of 8 artificial collectors were installed at each station. The replicate samples were placed at least 3 m distant from each other to avoid the pseudo-replicate related issue.
At the Vullatura, an additional 4 replicates of scouring pads for each station were placed on the dead matte substratum (M), without Posidonia cover, to assess P. oceanica buffering effects on mollusc settlement.
The artificial collectors were collected on the 30th of May 2019 and fixed in an alcohol 70% solution. The scouring pads were opened in the laboratory, and the material was washed in running water and filtered with a 120 μm mesh sieve size to retain all of the organisms and detritus collected. Organisms were sorted into large homogeneous taxonomic groups (i.e., Polychaeta, Gastropoda, Bivalvia, Amphipoda Tanaidacea, Isopoda, Cumacea, Decapoda, Echinodermata) and counted under a stereomicroscope. As for Mollusca, all the specimens were identified to the lowest possible taxonomical level, and their names crosschecked with online resources (WORMS; www.marinespecies.org, accessed on 30 April 2021).
During both the installation and collection of the scouring pads, 15 orthotropic P. oceanica shoots were randomly collected at each station and then transported inside small plastic bags filled with seawater in a cooler box to the laboratory. For each shoot, the lengths of the longest leaf were measured. Furthermore, to characterise the epiphytic assemblages, we selected the external side of the two outermost adult leaves of each shoot (generally the oldest and therefore the most colonised by epiphytes). The percentage coverage of the main taxonomic groups was visually assessed: encrusting brown algae, encrusting red algae (mainly Corallinales), erect algae (both brown and green algae), bryozoans, hydrozoans, spirorbid worms, and foraminifers were taken into account. The percentage coverage of the main taxonomic groups was visually estimated through the semi-quantitative Braun–Blanquet categories: 0 = absent; 1 = occasionally present, <5%; 2 = 5–25%; 3 = 26–50%; 4 = 51–75%; 5 = >75%.

2.3. Statistical Analyses

Student’s t-tests were carried out to test significant differences of P. oceanica mean leaf length according to pH conditions and sampling time (April–May). The number of taxa and number of specimens counted in each artificial collector were used as primary descriptors of mollusc assemblages, and were analysed using Generalized Linear Mixed Effect Models (GLMMs) to assess their variation along the OA gradient. We selected a Poisson distribution with a logarithmic link (typical for count data), and the site was set as a random effect. Then, the number of taxa and the number of specimens were tested in relation to two fixed nested factors: pH conditions (3 levels: control pH, low pH, and extremely low pH) and Microhabitat (2 levels: Canopy and Bottom). Conversely, Generalized Linear Models (GLMs) were used to assess changes in the number of taxa and the number of individuals at the Vullatura in relation to the presence or the absence of the seagrass: pH conditions (3 levels) and Microhabitat [3 levels: C, B, and Matte(M)] were included in the model formulas as fixed and nested independent variables.
Multivariate permutational analyses of variance (PERMANOVA) based on Bray–Curtis similarity were used to test differences in the mollusc assemblages among sites, pH conditions, and microhabitats. Furthermore, PERMANOVA was applied to test differences in the mollusc assemblages among pH conditions and the three microhabitats sampled at the Vullatura. Furthermore, similarity percentage–species contribution analyses (SIMPER) were carried out to identify the taxa that mainly contributed to the dissimilarity among stations previously highlighted by PERMANOVA. In order to visualise dissimilarities among mollusc assemblages sampled at different pH conditions and microhabitats, Principal Coordinates Analyses (PCoA) based on Bray–Curtis similarity were applied.
For GLMMs, GLMs, and PERMANOVA, a significance level of 0.05 (p-value < 0.05) was chosen. All of the statistical analyses were performed in the R platform (version 3.5.2).

3. Results

3.1. pH Measures

The pH characterisation of the sites and stations, derived from previous measurements using in situ continuous SeaFET pH meters [13,53,59], showed a high daily variability at both low pH and extremely low pH stations if compared to the ambient pH condition of control stations (Supplementary Materials, Figure S2). In-situ pH measures highlighted severe acidification at the extremely low pH stations (N3: 7.2 ± 0.3; S3: 6.6 ± 0.4; Vu3: 7.2 ± 0.5), if compared to low pH stations (N2: 7.5 ± 0.2; S2: 7.7 ± 0.3; Vu2: 7.8 ± 0.3). Minimum pH values of 6.0 and 6.1 were recorded at stations S3 and Vu3, respectively.

3.2. Posidonia Oceanica Leaf Phenology and Epiphytic Assemblages

The mean length of Posidonia leaves was significantly higher in the control stations (t-test, p-value < 0.001), in respect to the acidified station, and no significant differences (t-test, p-value > 0.05) were reported between low and extremely low pH conditions (Supplementary Materials, Table S1 and Figure S3). The mean leaf length significantly differed among sampling months, with values in May higher than in April.
Bryozoans and hydrozoans mainly characterised the epiphytic assemblages in April 2019; erect fleshy algae were reported in the acidified stations (both low and extremely low pH), whereas encrusting red algae (Corallinales) and calcareous invertebrates (spirorbids and foraminifers) were found mainly in control stations (Supplementary Materials, Figure S4). In May 2019 as well, the percentage coverage of calcareous organisms was low in the acidified stations (Figure S4).

3.3. Effects of OA on Mollusc’ Settlement

A total of 19,613 invertebrate specimens were collected and sorted at the three investigated sites and under different pH conditions. Invertebrate assemblages were quite similar among sites in the proportions of the different taxonomic groups. The structure of the assemblages found in the artificial collectors (Table 1) revealed a dominance of Polychaeta and Amphipoda: the former showed lower abundance at the acidified stations (from 37% to 33%) if compared to the ambient pH stations. On the other hand, the Amphipoda were more abundant at low and extremely low pH stations (from 19% to 31% and 25%). Mollusca (Gastropoda, Bivalvia, and Polyplacophora pooled together) accounted for 30% of the number of individuals at the three investigated sites. Each of the other taxonomic groups (i.e., Tanaidacea, Decapoda, Cumacea, Isopoda, and Echinodermata) never exceeded 10% of the total number of individuals.
The mollusc specimens sorted and identified were 4659, belonging to 57 different taxa (15 Bivalvia, 41 Gastropoda, and 1 Polyplacophora; Table 2). Many of the collected specimens were juveniles, demonstrating the effectiveness of scouring pads in collecting these stages and assessing invertebrate settlement. Mytilus galloprovincialis juv. was the main representative species, accounting for 63% of Bivalvia and 33% among all Mollusca (1540 individuals); in extremely low pH stations, this percentage reached 45% of the total mollusc individuals. Other common taxa found in the samples were Alvania lineata (360 individuals: 7.7%) Phorcus sp. juv. (probably Phorcus turbinatus) (310, 6.6%), Rissoa variabilis (276, 5.9%), Venerupis sp. juv. (232, 5.0%), Hiatella arctica (191, 4.0%), Eatonina fulgida (147 specimens, 3.2%), and Ecrobia ventrosa (132 specimens, 2.9%) (Supplementary Materials, Figure S5).
The number of taxa showed higher values from extremely low to ambient pH at the investigated sites, while the number of individuals showed comparable values, except at the S extremely low pH station where a peak was observed (Figure 2). The mean values of the two descriptors were lower at N and Vu stations if compared with S ones, where variability among replicates was also wider. The maximum number of individuals (159.0 ± 161.0) was reached at the extremely low pH S station, whereas the minimum value (4.7 ± 3.0) was reported for low pH station at the Vu site. The mean number of taxa did not show any relevant difference between canopy (C) and bottom (B) layers, especially between low and extremely low pH conditions. The intercepts of the GLMMs carried out on both the number of taxa and individuals were deemed significant and were defined by ambient pH and microhabitat B (Table 3). The number of taxa was significantly lower under extremely low and low pH conditions; all of the factors considered, and their relative interactions, yielded significant p-values for the number of individuals. Mollusc abundance was significantly lower under low and extremely low pH, especially in microhabitat C. Conversely, the interaction between low and extremely low pH conditions and microhabitat C showed a significant increase in the number of individuals if compared to the combination defined by the intercept.
The PCoA ordination plots of the mollusc assemblages revealed distinct clusters according to pH conditions and microhabitats (Figure 3): the three sites (N, S, and Vu) were shown in separate plots to avoid site-related effects. The separation between acidified and ambient stations was appreciable in Vu and N plots, whereas the S plot showed replicate dispersion, indicating higher variability among assemblages. Furthermore, mollusc assemblages differed between C and B microhabitats, which revealed dissimilarity between the layers within the meadows. Significant differences were also highlighted through PERMANOVA (Table 4) when analysing mollusc assemblage composition for all the factors considered and their relative interactions.
The similarity percentages–species contribution analysis (SIMPER) indicated high dissimilarity among pH conditions at the three sites, ranging from 51.4% to 88.4%. The taxa most responsible for the 60% dissimilarity among stations (SIMPER analyses) are shown in Table S2 (Supplementary Materials). Alvania lineata, Alvania sp. juv, E. fulgida, H. arctica, M. galloprovincialis, Phorcus sp. juv, and R. variabilis mostly contributed for among pH conditions dissimilarity and were mainly found in samples exposed to extremely low pH conditions.

3.4. Does Posidonia Oceanica Buffer the Effect of OA on Mollusc Settlement?

The number of taxa and the number of individuals did not show a visible differentiation pattern among the three investigated microhabitats (C, B, and M) at Vu site (Figure 4). Both the descriptors, however, showed higher mean values in samples collected at control pH stations. The maximum value of the number of taxa was reported (11.5 ± 3.1) in C, whereas the maximum number of individuals was collected in M (40.2 ± 18.6), both under control pH conditions. The intercepts of both the GLMs were significant and defined by control pH conditions and B microhabitat (Table 5). The number of taxa significantly differed at low pH exclusively; the other factors and their combination analysed in the model did not yield significant p-values. The GLM carried out on the number of individuals found on artificial collectors showed significant changes according to pH conditions, microhabitats, and their combination. In particular, the highest significant reduction in mollusc abundance was reported for the combination defined by extremely low pH condition and microhabitat M.
The PCoA revealed diversification of the mollusc assemblages according to microhabitats (Figure 5). For example, the observations from the canopy (C) produced a separate cluster if compared to the samples of B and M. Within each of the considered microhabitats, the four replicates exposed to a given pH condition were closer to each other, although mollusc assemblages collected on dead matte (M) were characterised by higher dispersion in the Euclidean plane, reflecting considerable variability.
The observed differences in the composition of mollusc assemblages and relative abundance of the taxa were significant according to the pH condition, microhabitat, and pH × microhabitat interaction (PERMANOVA analysis; Table 6). Alvania lineata, R. variabilis, E. ventrosa, and H. arctica were the species most responsible for dissimilarity among pH conditions within each microhabitat, and those characterised by higher abundances under acidified conditions (SIMPER, Supplementary Materials, Table S3).

4. Discussion

Ocean acidification modifies the structure and composition of marine benthic communities; these aspects and processes can be addressed in studying naturally acidified systems, such as CO2 vents. Although vents are not perfect predictors of changes that will occur at larger scale in future oceans, they provide information of the possible effects on local marine communities [61], otherwise impossible to achieve through mesocosms or Free Ocean Carbon Dioxide Enrichment (FOCE) systems [62]. Vents acidify surrounding waters and, as open systems, allow habitat connectivity through natural larval exchange, allowing researchers to study ecosystem processes, such as benthic community dynamics and changes under OA conditions [9].
In the three systems studied here (Castello north and south, and Vullatura), the pH exhibited wide daily fluctuations, especially at low and extremely low pH stations. This variability reflects local conditions at a small–medium spatial scale and in time ([15] and references therein), influenced by hydrodynamics, vents emissions, and biological activity (i.e., photosynthesis and respiration rates).
Our results show that settlement of mollusc taxa in P. oceanica meadows subjected to OA is highly influenced by the decrease of seawater pH at the three vents systems studied, although some species are less affected by these environmental conditions. Furthermore, this study highlights the role of P. oceanica meadows as a nursery area, buffering the effects of OA on the associated biota.
Posidonia oceanica leaf length showed significant variations among pH conditions, with lower values reported for low and extremely low pH stations. This pattern can be linked to the greater vulnerability of the seagrass under acidified conditions, due to the biological erosion performed by Sarpa salpa (Linnaeus, 1758) through grazing [13,58]. In fact, rcent studies reported a higher abundance of the herbivorous fish S. salpa in the vents area if compared to adjacent ambient pH stations [63,64]. Posidonia oceanica leaves seem more palatable under acidified conditions, due to their nutritional quality [65,66,67] and the reduction or disappearance of calcareous epiphytes [13]. The alterations in the epiphytic community structure were evident at the acidified stations at the three studied sites. Calcareous algae and invertebrate taxa (i.e., encrusting red algae, spirorbids, and foraminifers) were the dominant groups in control stations, whereas they strongly reduced their abundance at acidified stations. At low and extremely low pH stations, the epiphytic communities shifted towards encrusting and erect non-calcareous algae and non-calcareous invertebrates. The results here provided are consistent with similar studies carried out in the same area, highlighting the reduction of the diversity of the P. oceanica epiphytic communities exposed to OA [13,58].
The scouring pads, used as artificial collectors, were effective tools in collecting both juvenile and adult invertebrate stages, and therefore assessing settlement in P. oceanica meadows. The artificial collectors provided the different composition of the assemblages according to the pH gradient, and the structure of the assemblages. As for the main taxonomic group collected and their relative abundance, this was consistent with results reported by Ricevuto et al. [60], a study that was performed on the Castello south and north side vents.
The molluscs sampled through artificial collectors showed a relatively high number of species, especially considering the relatively shallow depth examined and compared with other Posidonia beds studied around Ischia [42,50]. Molluscs sampled in several studies, and different habitats around the Castello vents (including control areas outside the CO2 bubbling), accounted for 87 species [15,68], while molluscs collected in a previous study at the Castello Posidonia meadows accounted for 51 species [54,69]. Therefore, the diversity observed in our artificial collectors is comparable to the natural habitat surrounding the vents. Most of the species found are typical of P. oceanica meadows, especially among the gastropods [70], while many have already been recorded in the Castello area, associated with both Posidonia and macroalgae bottoms [15,54,71]. Most of the individuals collected were composed mainly of juvenile stages. Posidonia oceanica meadows are considered nursery areas for many invertebrates and fishes [72]; however, this role has not been well evaluated and supported in specific studies, especially for mobile macro-invertebrates. Posidonia vents systems have been recently indicated as nursery grounds by Mirasole et al. [64] for the fish fauna of the Castello vents, while juveniles of invertebrates have been largely sampled by similar artificial collectors placed next to hard vegetated bottoms off the Castello vents [19,60].
Our results provide evidence that, although molluscs actively recruit in all of the Posidonia stations here studied, they exhibited a significant reduction in species settlement under OA conditions, showing sensitivity to pH reduction. The lower number of taxa and number of individuals were observed according to pH gradient, although at the Vullatura and North Castello site, extremely low pH conditions showed slightly higher values of both parameters if compared to low pH stations. As for the number of individuals, this could be due to the reduced competition among species under more severe OA conditions. In addition, the high number of individuals recorded on the South Castello side could be related to a storm (SE wind direction) that occurred in May 2019, and which may have contributed to the resuspension of organisms from the bottom. However, the south side of the Castello, which is more sheltered to the dominant winds on the area (N–W), showed the highest Posidonia shoot density at the vents area [54], and it was characterised by a higher number of individuals collected in previous studies in the hard bottoms [53], Posidonia meadow [54] and artificial collectors [19,60]. In general, such between-site differences in species diversity and abundance, although showing a similar distributional pattern across pH conditions, suggest that various local factors may have as much of an effect on mollusc settlement. However, mollusc settlement changed according to decreasing pH, although with different magnitudes according to site. The general pattern indicates that the number of individuals differed by 0% to 25%, and the number of taxa differed by 10% to 40% from ambient to extremely low pH. In particular, the mollusc diversity reduction according to OA conditions was evaluated in previous studies, considering all habitats and data available up to 2018 [15,73]. Data are consistent with those reported in this study: on the north side, mollusc’s diversity changed by 39% between ambient and low pH, and by 48% between ambient and extremely low pH conditions, while on the south side the number of mollusc taxa was reduced by 47% between ambient and low pH, and by 53% between ambient and extremely low pH conditions.
This general pattern reflects what was previously observed in the same area in studies where molluscs and other calcifying invertebrates were considered, and it is related to the sensitivity of calcareous organisms to OA [9,18,53,60]. Furthermore, our results are consistent with similar studies carried out in Kueishan Island (Taiwan) vent systems [74,75], where water physio-chemical features, modified by vent emissions (although different from those characterising our Ischia vents), shape the composition and structure of mollusc assemblages.
However, some species showed variable trends according to pH reduction, with the gastropods A. lineata, E. fulgida, Phorcus sp. juv, R. variabilis, and the bivalves A. echinata, H. arctica, M. subpictus, and M. galloprovincialis, as the most frequent species occurring at low and extremely low pH. Mytilus galloprovincialis juvenile stages (2–3 mm length) showed a discrete tolerance to acidified conditions and accounted for 45% of individuals sampled in extremely low pH stations. High M. galloprovincialis settlement under acidified conditions was described previously at the Castello north and south sides using artificial collectors placed next to hard vegetated bottoms [60]. Although the juveniles of M. galloprovincialis were dominant in samples at the acidified stations, adults of the species are absent in the extremely low pH zones and very scarce at low pH, and occur mainly on the control areas far from the CO2 emissions [71]. The absence of adults of several taxa in the adjacent areas suggests that selection produced by OA happens subsequently, during the adult phase rather than in the settlement and post-settlement stage. As evidence of a post-settlement selection, the comparison of the adult mollusc assemblages found in the Castello Aragonese vent areas are provided as Supplementary Materials (Table S4; [15,54,69]); in our acidified stations, 30 taxa (52.6%), over the total of 57 taxa identified in this study, had been reported in adult stages in Posidonia meadows or in the vegetated hard bottoms inside or near the vent’s area. Also, natural/control stations showed a high rate of post-settlement selection (54.9% of the total identified taxa). However, several natural factors, besides OA, cause the mortality of early juvenile invertebrates, such as biological and physical disturbance, delay in metamorphosis, predation, and competition [76]. Likely, OA increases the early life stages and adult mortality rate, acting as a steady physiological stress inducing the loss of functional and taxonomic diversity at the vent’s sites [10].
The occurrence of A. echinata, H. arctica, and M. subpictus individuals at acidified stations seems to confirm the robustness of heavily calcified bivalves to a decrease of pH [54]. No visible sign of shell erosion, or differences in shell thickness among control and acidified conditions, were reported for the mollusc specimens found in the artificial collectors exposed to OA conditions. Likely, the short period of exposure to low/extremely low pH (33 days), as well as the natural fluctuation of pH values at the studied sites, can explain the observed pattern. Shell erosion has been observed in different adult mollusc species, such as Patella spp., Phorcus turbinatus (Born, 1778) and Columbella rustica (Linnaeus, 1758), in both vegetated hard bottoms and Posidonia meadows of the Castello vents [9,69].
Among gastropods, many herbivorous taxa can benefit from CO2 increases at vent sites through the provision of higher habitat complexity and food, as an effect of the increase in primary production [77]. In particular, R. variabilis deserves particular mention, since this species showed sensitivity to pH reduction in laboratory trials, with alteration in its ability to perceive chemical cues from prey [78]; however, R. variabilis did not show changes in settlement along the pH gradient.
No relevant differences in species composition, richness, and abundance were reported between canopy (C) and bottom (B) samples, except for stations at low and extremely low pH on the south side. At this site, the bottom samples hosted a higher number of taxa if compared to the canopy ones, while abundance showed a peak in the canopy at extremely low pH. This pattern may be related to high shoot density values reported in the south side vents (>1000 shoots/m2). In this zone, a similar pattern has been observed on the invertebrate community associated with Posidonia by Garrard et al. [54], where the abundance recorded in the acidified zone (low pH) was doubled if compared to the control zone.
The mollusc species showing a certain degree of robustness to OA may take advantage of the absence of the other species selected against water acidification, and, therefore, reduced competition. A similar mechanism of community compensation has been observed for the benthic community associated with the vegetated hard bottom of both sides of the Castello vents [53].
The samplings performed, especially at the Vullatura vent and control stations, considering three different microhabitats, offered an interesting perspective about the buffering effect provided by the P. oceanica cover. Mollusc assemblages showed qualitative and quantitative differences among microhabitats, due to the change of multiple environmental parameters, such as hydrodynamics and light, and the plant cover, which occurs from the canopy to the bottom of the meadow and in the dead matte. In our samples, and in all three of the sites considered, canopy and bottom samples always are separated (see PCoA models), likely suggesting a certain degree of buffering of the canopy. At the Vullatura, samples placed on the dead matte were more similar to those collected at the bottom level (B, among rhizomes). They showed lower diversity and abundances in the acidified stations if compared to the bottom and canopy samples. At the same time, these two descriptors were similar to those reported for the canopy level and even higher than the bottom samples in the Vullatura control station. The mollusc assemblage settled on the dead matte also showed less consistency in the PCoA plane according to pH gradient. This pattern suggests that when living Posidonia shoots are present, they can buffer the effect of pH reduction on associated fauna, likely through photosynthesis processes, oxygen production, and carbon dioxide removal [62,79,80]. Buffering effects on the local pH have been highlighted recently in Posidonia and other seagrasses by an increasing number of studies [32,33,34]; however, the direct or indirect effects on the local fauna have been poorly considered. Also, in the study carried out by Garrard et al. [54], the highest abundances and richness recorded in Posidonia under acidified conditions have been attributed more to the indirect effect of the local increase of the habitat complexity (i.e., shoot density) than to the direct effect of buffering low pH. Especially for small invertebrates with cryptic life habits, as most of the invertebrates associated with Posidonia are, dense seagrass and macroalgal cover may represent a “chemical” refuge from local pH fluctuations and other stressing conditions. Therefore, even calcifying organisms, as molluscs, could be less disadvantaged from the effect of the alteration of carbonate chemistry within the canopy of this seagrass.
Although preliminary to properly considering this aspect, the present study suggests that P. oceanica plays a pivotal role in the survival of several associated species, both as a nursery area and also as a refuge, via the mitigation of future climate change through contrasting OA, by oxygen production, CO2 removal, and by carbon storage. Therefore, both high habitat complexity and climate mitigation effects contribute to making P. oceanica meadows an increasingly essential marine ecosystem for the future of the Mediterranean Sea, so that protection and conservation measures, such as the EU Biodiversity Strategy for 2030, should be encouraged to maintain this habitat to be as healthy, diverse, and productive as possible.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/d13070311/s1, Figure S1: picture of P. oceanica subject to CO2 bubbling at the S3 station (Castello south side vents), Figure S2: pH trend at the studied stations, Table S1: leaves length values and t-test pairwise comparisons, Figure S3: leaf length of P. oceanica measured during the two sampling periods, Figure S4: composition of epiphytic assemblages on P. oceanica leaves, Figure S5: common taxa found in the artificial collectors, Table S2: results of the SIMPER analyses performed among pH conditions at the investigated vent sites, Table S3: results of the SIMPER analyses performed among the three investigated microhabitats at the Vullatura vent site, Table S4: comparison among post-settlement molluscs described in the present study and the checklist of the mollusc taxa recorded in the CO2 vent sites of the Castello Aragonese.

Author Contributions

Conceptualisation, M.C.G. and E.C.; methodology, M.C.G. and E.C.; formal analysis, A.B. and L.C.; data curation, A.B., L.C. and E.C.; writing—original draft preparation, A.B., L.C., M.C.G. and E.C.; writing—review and editing, G.A., M.C.G. and E.C.; supervision, G.A., M.C.G. and E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study is available on request from the corresponding author.

Acknowledgments

Thanks are due to Claudia Meneghesso (supported by a fellowship of the Stazione Zoologica Anton Dohrn of Napoli) and Davide Forestiere for sorting part of the benthic material considered in this paper, and Martina Gaglioti for help at sea in April 2019. We are extremely grateful to Luigia Donnarumma for help in the identification of some species. We wish to thank Juhyung Lee for letting us show the trend of the pH measurements recorded in some of the studied stations with the SeaFET pH-meters. Finally, we are grateful to Pietro Sorvino, and the ANS diving staff of Ischia, for the thorough support during the field activities at sea.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shukla, P.R.; Skea, J.; Calvo Buendia, E.; Masson-Delmotte, V.; Pörtner, H.O.; Roberts, D.C.; Zhai, P.; Slade, R.; Connors, S.; Van Diemen, R. IPCC, 2019: Climate Change and Land: An IPCC Special Report on Climate Change, Desertification, Land Degradation, Sustainable Land Management, Food Security, and Greenhouse Gas Fluxes in Terrestrial Ecosystems; Intergovernmental Panel on Climate Change (IPCC): Geneva, Switzerland, 2019. [Google Scholar]
  2. Caldeira, K.; Wickett, M.E. Anthropogenic carbon and ocean pH. Nature 2003, 425, 365. [Google Scholar] [CrossRef]
  3. Solomon, S.; Qin, D.; Manning, M.; Chen, Z.; Marquis, M.; Averyt, K.B.; Tignor, M.; Miller, H.L. Climate Change 2007: The Physical Science Basis; Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  4. Micheli, F.; Halpern, B.S.; Walbridge, S.; Ciriaco, S.; Ferretti, F.; Fraschetti, S.; Lewison, R.; Nykjaer, L.; Rosenberg, A.A. Cumulative human impacts on Mediterranean and Black Sea marine ecosystems: Assessing current pressures and opportunities. PLoS ONE 2013, 8, e79889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Fabry, V.J.; Seibel, B.A.; Feely, R.A.; Orr, J.C. Impacts of ocean acidification on marine fauna and ecosystem processes. ICES J. Mar. Sci. 2008, 65, 414–432. [Google Scholar] [CrossRef]
  6. Kroeker, K.J.; Kordas, R.L.; Crim, R.N.; Singh, G.G. Meta-analysis reveals negative yet variable effects of ocean acidification on marine organisms. Ecol. Lett. 2010, 13, 1419–1434. [Google Scholar] [CrossRef]
  7. Raven, J.; Caldeira, K.; Elderfield, H.; Hoegh-Guldberg, O.; Liss, P.; Riebesell, U.; Shepherd, J.; Turley, C.; Watson, A. Ocean Acidification Due to Increasing Atmospheric Carbon Dioxide; The Royal Society: London, UK, 2005; ISBN 0854036172. [Google Scholar]
  8. Doney, S.C.; Fabry, V.J.; Feely, R.A.; Kleypas, J.A. Ocean acidification: The other CO2 problem. Ann. Rev. Mar. Sci. 2009, 1, 169–192. [Google Scholar] [CrossRef] [Green Version]
  9. Hall-Spencer, J.M.; Rodolfo-Metalpa, R.; Martin, S.; Ransome, E.; Fine, M.; Turner, S.M.; Rowley, S.J.; Tedesco, D.; Buia, M.C. Volcanic carbon dioxide vents show ecosystem effects of ocean acidification. Nature 2008, 454, 96–99. [Google Scholar] [CrossRef] [Green Version]
  10. Teixidó, N.; Gambi, M.C.; Parravacini, V.; Kroeker, K.; Micheli, F.; Villéger, S.; Ballesteros, E. Functional biodiversity loss along natural CO2 gradients. Nat. Commun. 2018, 9, 5149. [Google Scholar] [CrossRef] [PubMed]
  11. Vizzini, S.; Martínez-Crego, B.; Andolina, C.; Massa-Gallucci, A.; Connell, S.D.; Gambi, M.C. Ocean acidification as a driver of community simplification via the collapse of higher-order and rise of lower-order consumers. Sci. Rep. 2017, 7, 4018. [Google Scholar] [CrossRef] [Green Version]
  12. Auriemma, R.; De Vittor, C.; Esposito, V.; Gaglioti, M.; Gambi, M.C. Motile Fauna associated to Cystoseira brachycarpa along a gradient of Ocean Acidification at a vent system off Panarea (Aeolian Islands, Italy). Biol. Mar. Mediterr. 2019, 26, 216–219. [Google Scholar]
  13. Mecca, S.; Casoli, E.; Ardizzone, G.D.; Gambi, M.C. Effects of ocean acidification on phenology and epiphytes of the seagrass Posidonia oceanica at two CO2 vent systems of Ischia (Italy). Mediterr. Mar. Sci. 2020, 21, 70–83. [Google Scholar] [CrossRef]
  14. Kroeker, K.J.; Micheli, F.; Gambi, M.C. Ocean acidification causes ecosystem shifts via altered competitive interactions. Nat. Clim. Chang. 2013, 3, 156–159. [Google Scholar] [CrossRef]
  15. Foo, S.A.; Byrne, M.; Ricevuto, E.; Gambi, M.C. The carbon dioxide vents of Ischia, Italy, a natural system to assess impacts of ocean acidification on marine ecosystems: An overview of research and comparisons with other vent systems. Oceanogr. Mar. Biol. 2018, 237–310. [Google Scholar]
  16. Gonzalez-Delgado, S.; Hernández, J.C. The importance of natural acidified systems in the study of ocean acidification: What have we learned? Adv. Mar. Biol. 2018, 80, 57–99. [Google Scholar] [PubMed]
  17. Rastrick, S.S.P.; Graham, H.; Azetsu-Scott, K.; Calosi, P.; Chierici, M.; Fransson, A.; Hop, H.; Hall-Spencer, J.; Milazzo, M.; Thor, P. Using natural analogues to investigate the effects of climate change and ocean acidification on Northern ecosystems. ICES J. Mar. Sci. 2018, 75, 2299–2311. [Google Scholar] [CrossRef]
  18. Cigliano, M.; Gambi, M.C.; Rodolfo-Metalpa, R.; Patti, F.P.; Hall-Spencer, J.M. Effects of ocean acidification on invertebrate settlement at volcanic CO2 vents. Mar. Biol. 2010, 157, 2489–2502. [Google Scholar] [CrossRef] [Green Version]
  19. Ricevuto, E.; Kroeker, K.J.; Ferrigno, F.; Micheli, F.; Gambi, M.C. Spatio-temporal variability of polychaete colonization at volcanic CO2 vents indicates high tolerance to ocean acidification. Mar. Biol. 2014, 161, 2909–2919. [Google Scholar] [CrossRef]
  20. Boudouresque, C.F.; Bernard, G.; Bonhomme, P.; Charbonnel, E.; Diviacco, G.; Meinesz, A.; Pergent, G.; Pergent-Martini, C.; Ruitton, S.; Tunesi, L. Protection and Conservation of Posidonia oceanica Meadows; RAMOGE, RAC/SPA and GIS Posidonie Publications: Marseille, France, 2012; ISBN 2905540311. [Google Scholar]
  21. Salomidi, M.; Katsanevakis, S.; Borja, Á.; Braeckman, U.; Damalas, D.; Galparsoro, I.; Mifsud, R.; Mirto, S.; Pascual, M.; Pipitone, C.; et al. Assessment of goods and services, vulnerability, and conservation status of European seabed biotopes: A stepping stone towards ecosystem-based marine spatial management. Mediterr. Mar. Sci. 2012, 13, 49–88. [Google Scholar] [CrossRef] [Green Version]
  22. Vassallo, P.; Paoli, C.; Rovere, A.; Montefalcone, M.; Morri, C.; Bianchi, C.N. The value of the seagrass Posidonia oceanica: A natural capital assessment. Mar. Pollut. Bull. 2013, 75, 157–167. [Google Scholar] [CrossRef]
  23. Mazzella, L.; Buia, M.C.; Gambi, M.C.; Lorenti, M.; Russo, G.F.; Scipione, M.B.; Zupo, V. Plant-animal trophic relationships in the Posidonia oceanica ecosystem of the Mediterranean Sea: A review. Plant-Anim. Interact. Mar. Benthos 1992, 46, 165–187. [Google Scholar]
  24. Giakoumi, S.; Halpern, B.S.; Michel, L.N.; Gobert, S.; Sini, M.; Boudouresque, C.F.; Gambi, M.C.; Katsanevakis, S.; Lejeune, P.; Montefalcone, M.; et al. Towards a framework for assessment and management of cumulative human impacts on marine food webs. Conserv. Biol. 2015, 29, 1228–1234. [Google Scholar] [CrossRef]
  25. Boudouresque, C.F.; Mainesz, A. Découverte de l’herbier de Posidonie. In Cahier du Parc National du Port-Cros; Le-Livre: Sablons, France, 1982; Volume 4, pp. 1–81. [Google Scholar]
  26. Pergent, G.; Bazairi, H.; Bianchi, C.N.; Boudouresque, C.F.; Buia, M.C.; Calvo, S.; Clabaut, P.; Harmelin-Vivien, M.; Mateo, M.A.; Montefalcone, M. Climate change and Mediterranean seagrass meadows: A synopsis for environmental managers. Mediterr. Mar. Sci. 2014, 15, 462–473. [Google Scholar] [CrossRef] [Green Version]
  27. Mayot, N.; Boudouresque, C.F.; Leriche, A. Unexpected response of the seagrass Posidonia oceanica to a warm-water episode in the North Western Mediterranean Sea. C. R. Biol. 2005, 328, 291–296. [Google Scholar] [CrossRef] [PubMed]
  28. Marbà, N.; Duarte, C.M. Mediterranean warming triggers seagrass (Posidonia oceanica) shoot mortality. Glob. Chang. Biol. 2010, 16, 2366–2375. [Google Scholar] [CrossRef]
  29. Mateo, M.A.; Romero, J.; Pérez, M.; Littler, M.M.; Littler, D.S. Dynamics of millenary organic deposits resulting from the growth of the Mediterranean seagrass Posidonia oceanica. Estuar. Coast. Shelf Sci. 1997, 44, 103–110. [Google Scholar] [CrossRef]
  30. Mateo, M.A.; Serrano, O. The carbon sink associated to Posidonia oceanica. In Mediterranean Seagrass Meadows: Resilience and Contribution to Climate Change Mitigation; IUCN: Gland, Switzerland; Málaga, Spain, 2012. [Google Scholar]
  31. Pergent-Martini, C.; Pergent, G.; Monnier, B.; Boudouresque, C.F.; Mori, C.; Valette-Sansevin, A. Contribution of Posidonia oceanica meadows in the context of climate change mitigation in the Mediterranean Sea. Mar. Environ. Res. 2021, 165, 105236. [Google Scholar] [CrossRef]
  32. Hendriks, I.E.; Olsen, Y.S.; Ramajo, L.; Basso, L.; Steckbauer, A.; Moore, T.S.; Howard, J.; Duarte, C.M. Photosynthetic activity buffers ocean acidification in seagrass meadows. Biogeosciences 2014, 11, 333–346. [Google Scholar] [CrossRef] [Green Version]
  33. Bergstrom, E.; Silva, J.; Martins, C.; Horta, P. Seagrass can mitigate negative ocean acidification effects on calcifying algae. Sci. Rep. 2019, 9, 1932. [Google Scholar] [CrossRef]
  34. Ricart, A.M.; Ward, M.; Hill, T.M.; Sanford, E.; Kroeker, K.J.; Takeshita, Y.; Merolla, S.; Shukla, P.; Ninokawa, A.T.; Elsmore, K. Coast-wide evidence of low pH amelioration by seagrass ecosystems. Glob. Chang. Biol. 2021, 27, 2580–2591. [Google Scholar] [CrossRef]
  35. Wahl, M.; Schneider Covachã, S.; Saderne, V.; Hiebenthal, C.; Müller, J.D.; Pansch, C.; Sawall, Y. Macroalgae may mitigate ocean acidification effects on mussel calcification by increasing pH and its fluctuations. Limnol. Oceanogr. 2018, 63, 3–21. [Google Scholar] [CrossRef] [Green Version]
  36. Greiner, C.M.; Klinger, T.; Ruesink, J.L.; Barber, J.S.; Horwith, M. Habitat effects of macrophytes and shell on carbonate chemistry and juvenile clam recruitment, survival, and growth. J. Exp. Mar. Bio. Ecol. 2018, 509, 8–15. [Google Scholar] [CrossRef]
  37. Coll, M.; Piroddi, C.; Steenbeek, J.; Kaschner, K.; Ben Rais Lasram, F.; Aguzzi, J.; Ballesteros, E.; Bianchi, C.N.; Corbera, J.; Dailianis, T.; et al. The Biodiversity of the Mediterranean Sea: Estimates, Patterns, and Threats. PLoS ONE 2010, 5, e11842. [Google Scholar] [CrossRef] [Green Version]
  38. Casoli, E.; Bonifazi, A.; Ardizzone, G.; Gravina, M.F.; Russo, G.F.; Sandulli, R.; Donnarumma, L. Comparative analysis of mollusc assemblages from different hard bottom habitats in the central Tyrrhenian sea. Diversity 2019, 11, 74. [Google Scholar] [CrossRef] [Green Version]
  39. Rodriguez, S.R.; Ojeda, F.P.; Inestrosa, N.C. Settlement of benthic marine invertebrates. Mar. Ecol. Prog. Ser. 1993, 97, 193–207. [Google Scholar] [CrossRef]
  40. Ardizzone, G.D.; Belluscio, A.; Gravina, M.F.; Somaschini, A. Colonization and disappearance of Mytilus galloprovincialis Lam. on an artificial habitat in the Mediterranean Sea. Estuar. Coast. Shelf Sci. 1996, 43, 665–676. [Google Scholar] [CrossRef]
  41. Idato, E. Zonazione verticale della fauna vagile di strato foliare in una prateria di Posidonia oceanica (L.) Delile: 1. Molluschi. Boll. Malacol. 1983, 19, 109–120. [Google Scholar]
  42. Russo, G.F. Mollusk syntaxon of foliar stratum along a depth gradient in a Posidonia oceanica (L.) Delile meadow: Diel variability. In Proceedings of the International Workshop on Posidonia oceanica Beds, Ischia, Italie, 7–11 October 1985; GIS Posidonie Publishing: Marseille, France, 1985. [Google Scholar]
  43. Russo, G.F.; Fresi, E.; Buia, M.C.; Vinci, D. Malacofauna delle praterie a Posidonia oceanica della zona di Capo Passero (Sicilia Sud-Orientale): Analisi comparativa con i popolamenti dell’isola d’Ischia. Oebalia 1985, 22, 339–348. [Google Scholar]
  44. Gambi, M.C.; Lorenti, M.; Russo, G.F.; Scipione, M.B.; Zupo, V. Depth and seasonal distribution of some groups of the vagile fauna of the Posidonia oceanica leaf stratum: Structural and trophic analyses. Mar. Ecol. 1992, 13, 17–39. [Google Scholar] [CrossRef]
  45. Terlizzi, A.; Russo, G.F. Variabilità strutturale del popolamento a molluschi di alcune praterie di Posidonia oceanica del Mediterraneo. Biol. Mar. Mediterr. 1995, 2, 423–426. [Google Scholar]
  46. Terlizzi, A.; Russo, G.F. Analisi della dinamica annuale del taxocene a molluschi di una prateria superficiale di Posidonia oceanica: Confronto tra due diverse metodiche di campionamento. Biol. Mar. Mediterr. 1996, 3, 489–492. [Google Scholar]
  47. Bonfitto, A.; Fellegara, I.; Gillone, G. Sampling techniques and structure of the malacofauna associated to the rhizome zone in Posidonia oceanica (L.) Delile. Boll. Malacol. 1997, 33, 83–88. [Google Scholar]
  48. Martini, N.; Gillone, G.; Lombardi, C.; Sabelli, B. Mollusc community of a Posidonia oceanica (L.) Delile bed: Annual variability. Boll. Malacol. 2000, 36, 191–194. [Google Scholar]
  49. Beqiraj, S.; Kashta, L.; Kuci, M.; Kasemi, D.; Mato, X.; Gace, A. Benthic macrofauna of Posidonia oceanica meadows in the Albanian coast. Nat. Montenegrina 2008, 7, 55–69. [Google Scholar]
  50. Russo, G.F.; Terlizzi, A. Structural patterns in the mollusc assemblages of Posidonia oceanica beds: Methodological, edaphic or biogeographical product? Boll. Malacol. 1997, 33, 89–94. [Google Scholar]
  51. Albano, P.G.; Sabelli, B. The molluscan assemblages inhabiting the leaves and rhizomes of a deep water Posidonia oceanica settlement in the central Tyrrhenian Sea. Sci. Mar. 2012, 0, 721–732. [Google Scholar] [CrossRef] [Green Version]
  52. Urra, J.; Ramírez, Á.M.; Marina, P.; Salas, C.; Gofas, S.; Rueda, J.L. Highly diverse molluscan assemblages of Posidonia oceanica meadows in northwestern Alboran Sea (W Mediterranean): Seasonal dynamics and environmental drivers. Estuar. Coast. Shelf Sci. 2013, 117, 136–147. [Google Scholar] [CrossRef]
  53. Kroeker, K.J.; Micheli, F.; Gambi, M.C.; Martz, T.R. Divergent ecosystem responses within a benthic marine community to ocean acidification. Proc. Natl. Acad. Sci. USA 2011, 108, 14515–14520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Gambi, M.C.; Gaglioti, M.; Teixidó, N. The CO2 vent’s systems off the island of Ischia (Tyrrhenian Sea). Mem. Descr. Cart. Geol. Italia 2020, 205, 55–64. (In Italian) [Google Scholar]
  55. Garrard, S.L.; Gambi, M.C.; Scipione, M.B.; Patti, F.P.; Lorenti, M.; Zupo, V.; Paterson, D.M.; Buia, M.C. Indirect effects may buffer negative responses of seagrass invertebrate communities to ocean acidification. J. Exp. Mar. Bio. Ecol. 2014, 461, 31–38. [Google Scholar] [CrossRef]
  56. Gambi, M.C. Emissioni sommerse di CO2 lungo le coste dell’isola d’Ischia. Rilievi su altre aree come possibili laboratori naturali per lo studio dell’acidificazione e cambiamento climatico a mare. Not. SIBM 2014, 66, 67–79. [Google Scholar]
  57. Buia, M.C.; Gambi, M.C.; Lorenti, M.; Dappiano, M.; Zupo, V. Aggiornamento sulla distribuzione e sullo stato ambientale dei sistemi a fanerogame marine (Posidonia oceanica e Cymodocea nodosa) delle isole Flegree. Acc. Sc. Lett. Arti Napoli Mem. Soc. Sc. Fis. Mat. 2003, 5, 163–186. [Google Scholar]
  58. Donnarumma, L.; Lombardi, C.; Cocito, S.; Gambi, M.C. Settlement pattern of Posidonia oceanica epibionts along a gradient of ocean acidification: An approach with mimics. Mediterr. Mar. Sci. 2014, 15, 498–509. [Google Scholar] [CrossRef] [Green Version]
  59. Lee, J. Global Change, Response Variability, and Their Implications for Coastal Ecosystem Resilience and Functioning; Stanford University: Stanford, CA, USA, 2020. [Google Scholar]
  60. Ricevuto, E.; Lorenti, M.; Patti, F.P.; Scipione, M.B.; Gambi, M.C. Temporal trends of benthic invertebrate settlement along a gradient of ocean acidification at natural CO2 vents (Tyrrhenian Sea). Biol. Mar. Mediterr. 2012, 19, 49–52. [Google Scholar]
  61. Riebesell, U. Acid test for marine biodiversity. Nature 2008, 454, 46–47. [Google Scholar] [CrossRef]
  62. Cox, T.E.; Gazeau, F.; Alliouane, S.; Hendriks, I.E.; Mahacek, P.; Le Fur, A.; Gattuso, J.P. Effects of in situ CO2 enrichment on structural characteristics, photosynthesis, and growth of the Mediterranean seagrass Posidonia oceanica. Biogeosciences 2016, 13, 2179–2194. [Google Scholar] [CrossRef] [Green Version]
  63. Mirasole, A.; Signa, G.; Gianguzza, P.; Bonaviri, C.; Mazzola, A.; Vizzini, S. Fish assemblages cope with ocean acidification in a shallow volcanic CO2 vent benefiting from an adjacent recovery area. Mar. Environ. Res. 2020, 157, 104851. [Google Scholar] [CrossRef] [PubMed]
  64. Mirasole, A.; Badalamenti, F.; Di Franco, A.; Gambi, M.C.; Teixidó, N. Boosted fish abundance associated with Posidonia oceanica meadows in temperate shallow CO2 vents. Sci. Total Environ. 2021, 771, 145438. [Google Scholar] [CrossRef]
  65. Arnold, T.; Mealey, C.; Leahey, H.; Miller, A.W.; Hall-Spencer, J.M.; Milazzo, M.; Maers, K. Ocean acidification and the loss of phenolic substances in marine plants. PLoS ONE 2012, 7, e35107. [Google Scholar] [CrossRef] [Green Version]
  66. Ricevuto, E.; Vizzini, S.; Gambi, M.C. Ocean acidification effects on stable isotope signatures and trophic interactions of polychaete consumers and organic matter sources at a CO2 shallow vent system. J. Exp. Mar. Bio. Ecol. 2015, 468, 105–117. [Google Scholar]
  67. Scartazza, A.; Moscatello, S.; Gavrichkova, O.; Buia, M.C.; Lauteri, M.; Battistelli, A.; Lorenti, M.; Garrard, S.L.; Calfapietra, C.; Brugnoli, E. Carbon and nitrogen allocation strategy in Posidonia oceanica is altered by seawater acidification. Sci. Total Environ. 2017, 607, 954–964. [Google Scholar] [CrossRef]
  68. Gaglioti, M.; Mirasole, A.; Teixidó, N.; Gambi, M.C. Finestre sul futuro climatico dei nostri mari: Sintesi sulla biodiversità marina del sistema di emissione di CO2 del Castello d’Ischia. Ambiente Rischio Comun. 2019, 16, 42–47. [Google Scholar]
  69. Garrard, S.L. The Effect of Ocean Acidification on Plant-Animal Interactions in a Posidonia Oceanica Meadow; Open University: Milton Keynes, UK, 2013. [Google Scholar]
  70. Donnarumma, L.; Sandulli, R.; Appolloni, L.; Russo, G.F. Assessing molluscs functional diversity within different coastal habitats of Mediterranean marine protected areas. Ecol. Quest. 2018, 29, 35–51. [Google Scholar] [CrossRef] [Green Version]
  71. Chiarore, A.; Bertocci, I.; Fioretti, S.; Meccariello, A.; Saccone, G.; Crocetta, F.; Patti, F.P. Syntopic Cystoseira taxa support different molluscan assemblages in the Gulf of Naples (southern Tyrrhenian Sea). Mar. Freshw. Res. 2019, 70, 1561–1575. [Google Scholar] [CrossRef]
  72. Guidetti, P. Differences among fish assemblages associated with nearshore Posidonia oceanica seagrass beds, rocky–algal reefs and unvegetated sand habitats in the Adriatic Sea. Estuar. Coast. Shelf Sci. 2000, 50, 515–529. [Google Scholar] [CrossRef]
  73. Gaglioti, M.; Ricevuto, E.; Gambi, M.C. Pattern and map of biodiversity related to ocean acidification in CO2 vents of Ischia. Biol. Mar. Mediterr. 2018, 25, 235–236. [Google Scholar]
  74. Chan, B.K.K.; Wang, T.W.; Chen, P.C.; Lin, C.W.; Chan, T.Y.; Tsang, L.M. Community structure of macrobiota and environmental parameters in shallow water hydrothermal vents off Kueishan Island, Taiwan. PLoS ONE 2016, 11, e148675. [Google Scholar] [CrossRef] [Green Version]
  75. Chen, C.; Chan, T.Y.; Chan, B.K.K. Molluscan diversity in shallow water hydrothermal vents off Kueishan Island, Taiwan. Mar. Biodivers. 2018, 48, 709–714. [Google Scholar] [CrossRef]
  76. Hunt, H.L.; Scheibling, R.E. Role of early post-settlement mortality in recruitment of benthic marine invertebrates. Mar. Ecol. Prog. Ser. 1997, 155, 269–301. [Google Scholar] [CrossRef] [Green Version]
  77. Connell, S.D.; Doubleday, Z.A.; Hamlyn, S.B.; Foster, N.R.; Harley, C.D.G.; Helmuth, B.; Kelaher, B.P.; Nagelkerken, I.; Sarà, G.; Russell, B.D. How ocean acidification can benefit calcifiers. Curr. Biol. 2017, 27, R95–R96. [Google Scholar] [CrossRef] [Green Version]
  78. Mutalipassi, M.; Fink, P.; Maibam, C.; Porzio, L.; Buia, M.C.; Gambi, M.C.; Patti, F.P.; Scipione, M.B.; Lorenti, M.; Zupo, V. Ocean acidification alters the responses of invertebrates to wound-activated infochemicals produced by epiphytes of the seagrass Posidonia oceanica. J. Exp. Mar. Bio. Ecol. 2020, 530–531, 151435. [Google Scholar] [CrossRef]
  79. Invers, O.; Zimmerman, R.C.; Alberte, R.S.; Pérez, M.; Romero, J. Inorganic carbon sources for seagrass photosynthesis: An experimental evaluation of bicarbonate use in species inhabiting temperate waters. J. Exp. Mar. Bio. Ecol. 2001, 265, 203–217. [Google Scholar] [CrossRef]
  80. Koch, M.; Bowes, G.; Ross, C.; Zhang, X.H. Climate change and ocean acidification effects on seagrasses and marine macroalgae. Glob. Chang. Biol. 2013, 19, 103–132. [Google Scholar] [CrossRef] [PubMed]
Diversity 13 00311 g001 550
Figure 1. Location of the study area and disposition of the sampling stations along a pH gradient. Control/ambient conditions (Vu1, N1, and SC) are highlighted in green, low pH conditions (Vu2, N2, and S2) are indicated with orange, and extremely low pH conditions (Vu3, N3, and S3) are highlighted in red. Study sites were distinguished as follows: Vu-Vullatura; N-North of Castello Aragonese; S-South of Castello Aragonese. Basemap source: ESRI.
Figure 1. Location of the study area and disposition of the sampling stations along a pH gradient. Control/ambient conditions (Vu1, N1, and SC) are highlighted in green, low pH conditions (Vu2, N2, and S2) are indicated with orange, and extremely low pH conditions (Vu3, N3, and S3) are highlighted in red. Study sites were distinguished as follows: Vu-Vullatura; N-North of Castello Aragonese; S-South of Castello Aragonese. Basemap source: ESRI.
Diversity 13 00311 g001
Diversity 13 00311 g002 550
Figure 2. Mean number of individuals (cyan) and mean number of taxa (dark blue) with SD bars according to site (N, S, and Vu), pH stations, and microhabitats within the meadow (C = canopy; B = bottom). Number of individuals’ axis break for the South (S) site is indicated in red.
Figure 2. Mean number of individuals (cyan) and mean number of taxa (dark blue) with SD bars according to site (N, S, and Vu), pH stations, and microhabitats within the meadow (C = canopy; B = bottom). Number of individuals’ axis break for the South (S) site is indicated in red.
Diversity 13 00311 g002
Diversity 13 00311 g003 550
Figure 3. Principal Coordinates Analysis (PCoA) of the mollusc assemblages sampled in N, S, and Vu sites. Colours indicate pH conditions: green = ambient pH; orange = low pH; red = extremely low pH. Microhabitats are indicated through circles (Canopy) and triangles (Bottom).
Figure 3. Principal Coordinates Analysis (PCoA) of the mollusc assemblages sampled in N, S, and Vu sites. Colours indicate pH conditions: green = ambient pH; orange = low pH; red = extremely low pH. Microhabitats are indicated through circles (Canopy) and triangles (Bottom).
Diversity 13 00311 g003
Diversity 13 00311 g004 550
Figure 4. Mean value of number of individuals (cyan) and mean number of taxa (dark blue) with SD bars in the three pH conditions, and microhabitats within the meadow (C = canopy; B = bottom, M = dead matte) at the Vu site.
Figure 4. Mean value of number of individuals (cyan) and mean number of taxa (dark blue) with SD bars in the three pH conditions, and microhabitats within the meadow (C = canopy; B = bottom, M = dead matte) at the Vu site.
Diversity 13 00311 g004
Diversity 13 00311 g005 550
Figure 5. Principal Coordinates Analysis (PCoA) ordination plot of the mollusc assemblages sampled at the Vu site. Colours indicate pH conditions: green = control pH; orange = low pH; red = extremely low pH. Microhabitats are indicated through circles (Canopy), triangles (Bottom), and squares (Matte).
Figure 5. Principal Coordinates Analysis (PCoA) ordination plot of the mollusc assemblages sampled at the Vu site. Colours indicate pH conditions: green = control pH; orange = low pH; red = extremely low pH. Microhabitats are indicated through circles (Canopy), triangles (Bottom), and squares (Matte).
Diversity 13 00311 g005
Table
Table 1. Composition of the benthic organisms collected in the samples in different sites and pH conditions and expressed as percentages (%) of each taxonomic group on the total number of individuals sampled.
Table 1. Composition of the benthic organisms collected in the samples in different sites and pH conditions and expressed as percentages (%) of each taxonomic group on the total number of individuals sampled.
SitespH Conditions
ClassVuNSControlLow pHExtremely Low pH
Polychaeta363633373233
Amphipoda242226193125
Tanaidacea974556
Gastropoda172210171212
Bivalvia7820151318
Decapoda000000
Cumacea111211
Isopoda125443
Echinodermata521122
Table
Table 2. Number of individuals divided per species in each site (Vu, N, S), pH gradient (control, low, extremely low) and microhabitat (C = “canopy”, B = “bottom”, M = “dead matte”).
Table 2. Number of individuals divided per species in each site (Vu, N, S), pH gradient (control, low, extremely low) and microhabitat (C = “canopy”, B = “bottom”, M = “dead matte”).
Control/Ambient pH Low pH Extremely Low pH
VuCVuBVuMNCNBSCSBVuCVuBVuMNCNBSCSBVuCVuBVuMNCNBSCSB
Bivalvia
Abra alba (W. Wood, 1802) 1 47
Acanthocardia echinata (Linnaeus, 1758)1125341814 3214210 1228
Anomia ephippium (Linnaeus, 1758)9 21 36 3 14
Arca noae (Linnaeus, 1758)91 1 5 1 1 2
Gregariella semigranata (Reeve, 1858)11 2
Haminoea sp. 1
Hiatella arctica (Linnaeus, 1767)1521016 2129 42819991 7 282
Lima sp. 104684691 13 535 3 1
Mimachlamys varia (Linneo, 1758)4 21752 1 81 1 3
Musculus subpictus (Cantraine, 1835)1232121621412 8143 311324
Mytilus galloprovincialis juv. (Lamarck, 1819)234111846140187 12322223121 1310340277
Parvicardium sp. 1 3 23 5
Pecten sp. 5 2 1
Venerupis sp.3242143050 5 2534121 3614
Veneridae gen. sp.7 4 21213 11715
Gastropoda
Alvania cimex (Linneo, 1758) 76 2 6
Alvania discors (T.Brown, 1818) 3 2 5
Alvania lineata (Risso, 1826)84724217120261815 62 7625 18226
Alvania pagodula (Bucquoy, Dautzenberg & Dollfus, 1884)220
Alvania piersmai (Moolenbeek e Hoenselaar, 1989) 3 1 2 2 8
Alvania scabra (Filippi, 1844) 5 1 9 3 4 8 1 1
Alvania sp. juv. 12 13 1 916 487 4
Alvania subareolata(Monterosato, 1869) 13 1
Ammonicera fischeriana (Monterosato, 1869)2115123 4 81 51 1
Bittium latreillii (Payraudeau, 1826) 4 2 1 1 1
Bolma rugosa (Linneo, 1767) 15
Cerithium vulgatum (Bruguière, 1792) 521 112 212
Crisilla beniamina (Monterosato, 1884)8 3 27
Eatonina fulgida (J. Adams, 1797) 27 1131 3164147 411922924
Ecrobia ventrosa (Montagu, 1803)3 541313 20 11015 11 71
Gibberula philippii (Monterosato, 1878) 1 7 2 2
Gibbula sp. 2411 1 2120 3 3428
Hydrobia acuta neglecta (Muus, 1963) 2 2 1
Jujubinus exasperatus (Pennant, 1777) 1 1
Jujubinus gravinae (Dautzenberg, 1881) 421923111 4 22 1 6 2
Jujubinus striatus (Linnaeus, 1758)1 34 65 26 44 1212 1
Melanella petitiana (Brusina, 1869) 1 118 21 1 1153
Mitrella sp. 1 2
Nodulus contortus (Jeffreys, 1856)1145219239 2 73 1
Nudibranchia gen. sp. 1 12 1
Phorcus sp. juv. 16945101 2644442 3813620
Pisinna glabrata (Megerle von Mühlfeld, 1824)2 11
Pusillina marginata (Michaud, 1830)314952094 115214 13 1085
Retusa truncatula (Bruguière, 1792)2 1
Rissoa italiensis (Verduin, 1985) 12 1 3 1
Rissoa lilacina (Récluz, 1843) 10 4 1 1 1 1
Rissoa splendida (Eichwald, 1830) 1625 112 7
Rissoa variabilis (Megerle von Mühlfeld, 1824)4 342401161121 3612 23133018168
Rissoa ventricosa (Desmarest, 1814) 1
Rissoa violacea (Desmarest, 1814) 1
Rissoella inflata (Locard, 1892) 633 6 45
Scissurella costata (d’Orbigny, 1824) 2 1
Tricolia entomochelia (Gofas, 1994) 1 2
Tricolia landinii (Bogi & Campani, 2007) 1 3 3 1
Tricolia pullus (Linneo, 1758) 7 5 8 15 5 2
Vitreolina philippii (Rayneval & Ponzi, 1854) 1 1 10131 2114
Polyplacophora
Leptochiton scabridus (Jeffreys, 1880) 1
Table
Table 3. Outputs of the GLMMs carried out on the number of taxa and number of individuals. Significant values are highlighted in bold.
Table 3. Outputs of the GLMMs carried out on the number of taxa and number of individuals. Significant values are highlighted in bold.
Source of VariationValueSD Errort-Valuep-Value
Number of taxa
Intercept2.6290.19413.567<0.001
Low pH−0.3550.118−3.0090.002
Extremely low pH−0.4140.120−3.449<0.001
Microhabitat (C)−0.0170.107−0.1620.871
Low pH × Microhabitat (C)0.0960.1650.5820.561
Extremely low pH × Microhabitat (C)−0.0260.174−0.1530.878
Number of individuals
Intercept4.0560.4199.667<0.001
Low pH−0.3820.052−7.307<0.001
Extremely low pH−0.1360.048−2.8000.005
Microhabitat (C)−0.1540.049−3.1530.001
Low pH × Microhabitat (C)0.2480.0743.338<0.001
Extremely low pH × Microhabitat (C)0.1710.0722.3790.017
Table
Table 4. Results of permutational multivariate analysis of variance (PERMANOVA) on mollusc assemblages found on P. oceanica meadows. Significant values are highlighted in bold.
Table 4. Results of permutational multivariate analysis of variance (PERMANOVA) on mollusc assemblages found on P. oceanica meadows. Significant values are highlighted in bold.
Source of VariationDfSSMSFp-Value
Site25.3760.23314.8400.001
pH21.0650.0462.9320.001
Microhabitat12.3290.10112.8330.001
Site × pH41.8140.0782.4980.001
Site × Microhabitat20.7760.0332.1360.003
pH × Microhabitat20.7160.0311.9720.010
Site × pH × Microhabitat41.3530.0581.8630.002
Residuals539.6230.417
Table
Table 5. Output of the GLMs fitted to the number of taxa and number of individuals in the three microhabitats at the Vu site. Significant values are highlighted in bold.
Table 5. Output of the GLMs fitted to the number of taxa and number of individuals in the three microhabitats at the Vu site. Significant values are highlighted in bold.
Source of VariationEstimateSD Errort-Valuep-Value
Number of species
Intercept2.0470.17911.401<0.001
Low pH−0.7950.322−2.4690.013
Extremely low pH−0.1380.263−0.5250.599
Microhabitat (C)0.3950.2321.6980.089
Microhabitat (M)0.2790.2381.1750.24
Low pH × Microhabitat (C)0.2620.4030.6500.515
Extremely low pH × Microhabitat (C)−0.4320.359−1.2020.229
Low pH × Microhabitat (M)−0.2110.441−0.4770.633
Extremely low pH × Microhabitat (M)−0.4390.370−1.1880.235
Null deviance: 62.86 on 35 df
Residual deviance: 30.58 on 27 df
Number of individuals
Intercept3.1670.10230.874<0.001
Low pH−1.6090.251−6.404<0.001
Extremely low pH0.3800.1332.8590.004
Microhabitat (C)0.3730.1332.8010.005
Microhabitat (M)0.5270.1294.078<0.001
Low pH × Microhabitat (C)0.7930.2942.6920.007
Extremely low pH × Microhabitat (C)−1.1180.200−5.583<0.001
Low pH × Microhabitat (M)−0.2140.328−0.6510.515
Extremely low pH × Microhabitat (M)−1.3350.200−6.669<0.001
Null deviance: 513.63 on 35 df
Residual deviance: 251.88 on 27 df
Table
Table 6. Results of PERMANOVA on mollusc assemblages found on P. oceanica meadow at the Vu site. Significant values are highlighted in bold.
Table 6. Results of PERMANOVA on mollusc assemblages found on P. oceanica meadow at the Vu site. Significant values are highlighted in bold.
Source of VariationDfSSMSFp-Value
pH21.6010.1222.9430.001
Microhabitat21.9780.1513.6380.001
pH × Microhabitat42.1470.1641.9740.010
Residual277.3410.561
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barruffo, A.; Ciaralli, L.; Ardizzone, G.; Gambi, M.C.; Casoli, E. Ocean Acidification and Mollusc Settlement in Posidonia oceanica Meadows: Does the Seagrass Buffer Lower pH Effects at CO2 Vents? Diversity 2021, 13, 311. https://doi.org/10.3390/d13070311

AMA Style

Barruffo A, Ciaralli L, Ardizzone G, Gambi MC, Casoli E. Ocean Acidification and Mollusc Settlement in Posidonia oceanica Meadows: Does the Seagrass Buffer Lower pH Effects at CO2 Vents? Diversity. 2021; 13(7):311. https://doi.org/10.3390/d13070311

Chicago/Turabian Style

Barruffo, Alessandra, Laura Ciaralli, Giandomenico Ardizzone, Maria Cristina Gambi, and Edoardo Casoli. 2021. "Ocean Acidification and Mollusc Settlement in Posidonia oceanica Meadows: Does the Seagrass Buffer Lower pH Effects at CO2 Vents?" Diversity 13, no. 7: 311. https://doi.org/10.3390/d13070311

APA Style

Barruffo, A., Ciaralli, L., Ardizzone, G., Gambi, M. C., & Casoli, E. (2021). Ocean Acidification and Mollusc Settlement in Posidonia oceanica Meadows: Does the Seagrass Buffer Lower pH Effects at CO2 Vents? Diversity, 13(7), 311. https://doi.org/10.3390/d13070311

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop