Potential Resilience to Ocean Acidiﬁcation of Benthic Foraminifers Living in Posidonia oceanica Meadows: The Case of the Shallow Venting Site of Panarea

: This research shows the results regarding the response to acidic condition of the sediment and Posidonia foraminiferal assemblages collected around the Panarea Island. The Aeolian Archipelago represents a natural laboratory and a much-promising study site for multidisciplinary marine research (carbon capture and storage, geochemistry of hydrothermal ﬂuids and ocean acidiﬁcation vs. benthic and pelagic organisms). The variability and the complexity of the interaction of the ecological factors characterizing extreme environments such as shallow hydrothermal vents did not allow us to carry out a real pattern of biota responses in situ, differently from those observed under controlled laboratory conditions. However, the study provides new insights into foraminiferal response to increasing ocean acidiﬁcation (OA) in terms of biodiversity, faunal density, speciﬁc composition of the assemblages and morphological variations of the shells. The study highlights how the foraminiferal response to different pH conditions can change depending on different environmental conditions and microhabitats (sediments, Posidonia leaves and rhizomes). Indeed, mineral sediments were more impacted by acidiﬁcation, whereas Posidonia microhabitats, thanks to their buffer effect, can offer “refugia” and more mitigated acidic environment. At species level, rosalinids and agglutinated group represent the most abundant taxa showing the most speciﬁc resilience and capability to face acidic conditions.


Introduction
Ocean acidification (OA) is one of the main consequences of climate change linked to the increasing atmospheric CO 2 concentration, affecting current and future marine ecosystems [1][2][3]. OA associated with global warming of 1.5 • C to 2 • C [3], is projected to amplify its adverse effects on the growth, development, calcification and survival of marine organisms [4,5]. Reduced calcification, reduced rates of repair, and weakened calcified structures are among the most significant effects to the future levels of OA expected by 2100. However, the responses are species-specific and some organisms have not shown high sensitivity to changes in CO 2 , pH and carbonate concentrations [6][7][8][9][10]. Among meiofaunal protists, foraminifers, single-celled protists, are by far the most useful meiofaunal group utilized by geoscientists in environmental and climatological reconstructions of recent and past geological record. Based upon their widespread distribution, short life and reproductive cycle, high biodiversity, and specific ecological requirements, foraminifers may quickly respond to environmental changes and thus can be used as bio-indicators in a wide range of marine environments, e.g., [11][12][13][14][15]. The calcification process in benthic foraminifers is sensitive both to ocean acidification and local chemical and physical features of seawater. Therefore, they can record specific environmental conditions and they can provide a possible response to the ongoing acidification trend giving a picture of the present day and future changes [16,17]. Shallow-water (less than 200 m water depth) hydrothermal fields provide the opportunity to study fluid impacts and the effects of elevated pCO 2 on the benthic foraminiferal communities [1, [18][19][20]. Among them, the seafloor around Panarea Island (Aeolian Archipelago) represents a natural laboratory and a much-promising study site for multidisciplinary marine research (carbon capture and storage, geochemistry of hydrothermal fluids and ocean acidification vs benthic and pelagic organisms). The seabottoms of Panarea Island offer a wide environmental variability linked to different degrees of hydrothermal activity varying from extreme environments to normal marine conditions. In addition, luxuriant Posidonia oceanica meadows host many ecological niches and different microhabitats useful to investigate the response of epiphytic community in which foraminifers represent one of the most conspicuous components [21][22][23][24]. Previous studies on benthic foraminiferal assemblages in Mediterranean areas characterized by underwater hydrothermal vents and CO 2 fluid emissions, have indicated biodiversity loss, decrease in faunal density and variations in the test-type abundances (agglutinated, porcelaneous, hyaline) [25][26][27]. Moreover, experimental data showed how P. oceanica meadows metabolism can modulate the responses of epiphytic micro/meio-invertebrates to future OA, conferring more uniform pH regimes [28][29][30].
The main part of the experimental studies was carried out under controlled conditions in laboratory where it is possible to observe the response of microfauna (as a whole in foraminiferal assemblages or at single-species level) to the variations of each parameter (temperature, pH) [28,31]. The relationships between the environmental parameters (e.g., pH, temperature, salinity, CO 2 ) and foraminiferal distribution in such environments in situ are yet poorly studied. Actually, the results obtained in the laboratory are very difficult to apply in nature because of the complexity of interaction of abiotic and biotic factors added to the temporal and spatial variability of these same factors.
On this perspective, the aim of this research is to present the results of the study of foraminiferal assemblages recorded both in sediments and Posidonia leaves off Panarea Island. Our purpose is to give new insights into the overall pattern of foraminiferal community response to the extreme environment of shallow hydrothermal vents in terms of biodiversity, specific composition of the assemblages, and morphological variations of the foraminiferal shells. These results could represent a background for biomonitoring the area in the future and to acquire data to elaborate short-time scenarios.
Panarea Island is part of the Aeolian Archipelago, a 200 km-long, ring-shaped partly submerged volcanic ridge, located in the SE sector of the Tyrrhenian Sea. The latter is a back-arc basin developed since the Neogene, as a consequence of the NW-directed subduction of the Ionian oceanic lithosphere below the Calabrian Arc [32][33][34][35][36]. Therefore, the archipelago is a subduction related volcanic arc [37,38], composed by stratovolcanoes formed by both explosive and effusive activity [39][40][41]. The volcanic products composition display calcalkaline, shoshonitic and potassic affinity, with an age ranging from 1.3 Ma to present [40][41][42]. Panarea Island and the surrounding islets represent the emerged portions of a wide, mostly submerged, volcanic edifice, more than 2000 m high, mainly dismantled by erosion and by neo-and volcano-tectonics activity [43][44][45]. The volcanic complex shows a conic shape, with a flat summit formed by a submerged shelf located at depths from 80 to 130 m b.s.l. To the east of Panarea, the islets of Basiluzzo, Panarelli, Dattilo, Bottaro, Lisca Bianca, Lisca Nera and Formiche represent an archipelago partly emerging from the flat eroded plateau, at a depth of less than 30 m [43]. The eruptive history of the volcanic system is formed by successive eruptive cycles of volcanic activity, separated by quiescent stages [46]. Panarea island is mostly composed of andesitic to dacitic lava domes, dated between 149 ± 5 and 127 ± 1.5 ka [46][47][48], interbedded with pyroclastic deposits. The islets are made of dacite and andesite lavas, the last of which were emplaced with the Basiluzzo endogenous dome (54 ± 8 ka) [47,49]. Currently, Panarea volcanic complex is characterized by local subsidence, observed in the Basiluzzo area [50] and by degassing from several fumarolic areas, both inland and offshore. Gas venting is frequent along a NE-SW trending extensional fault system, located NE of Panarea Island. Here a series of fault scarps and fresh volcanic rocks and mineralization are exposed, where gas venting has been observed [51,52]. These hydrothermal processes produced sediment-hosted sulphide deposits, Fe-rich crusts, bacterial mats and chimneys [51,[53][54][55]. Hydrothermalism is also attested by diffuse degassing at the seafloor, and gas emissions and boiling of the seawater has been described since historical times at Panarea Island [51,56]. The active vents are characterized by discharge of CO 2 -dominated gases and thermal fluids with temperatures up to 140 • C [57,58]. An intense hydrothermal activity was reached in November 2002: a gas eruption occurred 125 km east of Panarea island, between 2 and 30 m b.s.l., localized on top of a shallow rise, surrounded by the islets of Panarelli, Lisca Bianca, Bottaro, Lisca Nera and Dattilo [56]. The most intense gas flux was measured to the west of the islet of Bottaro, where one vent was erupting at 108-109 lt/day [59,60]. According to [57], the total gas output under normal conditions is about 9 × 106 lt/day over the whole area. Direct measurements at the site of eruption reported values up to 110 • C, with pH between 5.0 and 5.5 [61,62]. The active degassing phase lasted several months, up to mid-2003, and had great influence on the biota and water properties in the proximity of vents, such as the death of benthic communities, high bacteria abundance and diversity, water salinity increasing, pH lowering, strong Cl-excess with respect to seawater [63][64][65]. Progressively, the hydrothermal activity reduced its intensity, and gas vents shifted their position with time [56]. Nowadays, most of the venting activity is concentrated within the submarine crater formed close to Bottaro islet during the paroxysmal phase of November 2002 [56].

Sampling Strategy
In the present paper the authors show the results of a sampling survey carried out in the frame of the Scientific Diving Summer School hosted in the ECCSEL-NatLab Italy laboratory of Panarea Island on September 2019, an ECCSEL-ERIC infrastructure of the National Institute of Oceanography and applied Geophysics. The samples were collected, by scuba divers of the diving center Amphibia, from a depth interval ranging from 4 to 24 m off the NE coastal sector of the Panarea Island. Four areas affected by submarine hydrothermal manifestations with CO 2 dominated emissions (>90%) located at shallow depth around Panarea island were considered: Ditella areas (Hot and Cold), Bottaro crater, P21 and Black Point (Figure 1). These areas are well-known and documented and they are characterized by the presence of wide seagrasses patchy meadows, e.g., [57,58,[65][66][67][68][69]. Study on the Panarea bottoms showed that P. oceanica shoot and leaf density seems to be unaffected by OA [70].
Two sites (Raya and B1) displaying normal marine pH and temperature values (pH = 8.15, 23-26 • C) were considered as control sites. Temperature and pH values reported in Table 1 are unpublished data (personal communication). For Ditella stations T and pH values reply the values reported in literature [68].
Ditella area is located NE of Panarea Island at 10-12 m depth. This area is characterized by spots with very different temperatures at the distance of few meters from each other, which probably influence the distribution of the patchy meadow of P. oceanica seagrass. The hot spots are characterized by CO 2 emissions with a temperature up to 60 • and pH values ranging from 7 to 5.6 (parameters were measured at the emission points). The cold spots are characterized by reduced venting activity and display temperature of about 26 • C and pH of about 7.9. In this site Posidonia coverage appears continuous. Data collected during 2016 reported higher shoot density near the cold areas (365.6 ± 89.7 shoots/m 2 ) than near the hot plots (273.2 ± 99.2 shoots/m 2 ) [71].  Table 1 are unpublished data (personal communication). For Ditella stations T and pH values reply the values reported in literature [68].
Ditella area is located NE of Panarea Island at 10-12 m depth. This area is characterized by spots with very different temperatures at the distance of few meters from each other, which probably influence the distribution of the patchy meadow of P. oceanica seagrass. The hot spots are characterized by CO2 emissions with a temperature up to 60° and pH values ranging from 7 to 5.6 (parameters were measured at the emission points). The cold spots are characterized by reduced venting activity and display temperature of about 26 °C and pH of about 7.9. In this site Posidonia coverage appears continuous. Data collected during 2016 reported higher shoot density near the cold areas (365.6 ± 89.7 shoots/m 2 ) than near the hot plots (273.2 ± 99.2 shoots/m 2 ) [71].
Posidonia leaves and soft sediment samples were collected from each station to highlight differences between sediment and leaves foraminiferal assemblages.
The Bottaro crater is an area characterized by a depression (about 14 m wide, over 20 m long and about 11 m deep) located 3.3 km to the east of Panarea close to the Bottaro islet. This crater is the result of the massive gas eruption occurred in November 2002 [56,58]. In the Bottaro crater zone only Posidonia samples were considered (13 m depth) coming from two of the three sites (B1-B2). Literature data reported shoot density ranging from 674.7 ± 26.6 to 609.3 ± 91.8 shoots/m 2 [70]. Considering a profile of the crater, starting from the central site inside the crater, to the margin of the depression (B3) any Posidonia coverage was recorded. Inside the crater pH values range from 6.0 to 7.8 with a temperature of about 22 °C. B2 is 7 m from the edge of the crater and represents the first site in which Posidonia were recorded along the crater profile; temperature is 22.6 °C, while pH is about 7.9. Finally, B1 is about 30 m from the edge of the crater and is characterized by temperature of 22.7 °C and pH of 8.15. It was considered representative of normal marine conditions.  Posidonia leaves and soft sediment samples were collected from each station to highlight differences between sediment and leaves foraminiferal assemblages.
The Bottaro crater is an area characterized by a depression (about 14 m wide, over 20 m long and about 11 m deep) located 3.3 km to the east of Panarea close to the Bottaro islet. This crater is the result of the massive gas eruption occurred in November 2002 [56,58]. In the Bottaro crater zone only Posidonia samples were considered (13 m depth) coming from two of the three sites (B1-B2). Literature data reported shoot density ranging from 674.7 ± 26.6 to 609.3 ± 91.8 shoots/m 2 [70]. Considering a profile of the crater, starting from the central site inside the crater, to the margin of the depression (B3) any Posidonia coverage was recorded. Inside the crater pH values range from 6.0 to 7.8 with a temperature of about 22 • C. B2 is 7 m from the edge of the crater and represents the first site in which Posidonia were recorded along the crater profile; temperature is 22. by a rocky bottom with strong localized gas emissions adjacent to a tectonic fault [63,72]. The Posidonia cover is very patchy. These two sites represent the stations with the highest venting activity of all study sites with temperatures ranging from 55.7 • C (P21) to 133 • C (Black Point); pH value is 5.21 for P21 and ranges from 2.78 to 4.9 in the Black Point site. At these stations no sediment samples were collected due to the rocky feature of the bottom.
Finally, sediment samples coming from both the hot and cold sites and from an area far from the sites, have been mineralogically characterized by diffractometric analysis.

Benthic Microfauna
Two typologies of samples were considered: sediment and Posidonia (rhizomes and leaves). For each type of samples different treatments were adopted following, as far as possible, the standard procedures [21,24,73].
Sediment samples-Sediment was collected only in Ditella and Raya stations with three replicates for each site. Sediment samples were sampled by mean of the syringe technique (diameter 3 cm, 50 mL related to 10 cm top layer equivalent to 100 cc of sediment). The whole syringe content was analyzed because of the coarse grain size of the sediment that made difficult to slice 0.5 or 1 cm subsamples thick. To distinguish the living fauna, all sediment samples were stained and preserved in a solution of 90% ethanol with 2 g/L of Rose Bengal [73][74][75]. After 15 days, the samples were wet sieved through a 63 µm sieve and then dried at 40 • C. For each sample, stained foraminifers were counted, hand-picked, and identified using a binocular microscope. The Rose Bengal staining method has been widely used in ecological studies to distinguish live from dead foraminifers because it is inexpensive and easy to use [22,76,77]. However, under specific conditions (i.e., anoxic environments) the accuracy of this method may be affected by the presence of undecayed protoplasm, which can persist for weeks or months after death [76][77][78][79][80]. While the staining criteria are confidently applied to the superficial samples, ambiguities may arise in the case of deeper intervals [81], commonly consisting of a slight overestimation of the living assemblages [82]. However, the Rose Bengal method when used with care, permits the obtainment of results as reliable as those obtained by other techniques [75,80,83]. The quantitative analysis of benthic foraminifers was based on the count of all specimens present in the whole sample. In addition, the frequency of groups distinguished by their basic wall structure (agglutinated, porcelaneous and hyaline) was calculated in order to display their distribution in different pH conditions.
The scarce and irregular quantitative data didn't allow us to carry out statistical analyses. Density in the sediment samples was calculated as number of individuals in 100 cc of sediment, while diversity was calculated as total number of taxa in the sediments (S) and also α-Fisher [84] and Shannon (H) indices [85]. All these parameters were calculated by mean of the PAST software (ver.4.03) [86].
Posidonia samples-Posidonia rhizomes were cut off the substrate approximately one centimeter above the sediment surface. Leaves and rhizomes were immediately stowed away under water in plastic bags, and later carefully washed with seawater over 63 µm sieve. Epiphytes were washed into larger bowls, washed with fresh water and dried. Plant fragments were examined under the microscope to remove living specimens that may have remained glued to the plant surface. All epiphytic foraminifers were picked from each sample and identified at species level. Posidonia samples comprise material from approximately 27 rhizomes and 172 leaves ( Table 2). Similar length leaves were considered, and the ratio between number of specimens recorded and number of leaves analyzed (F/P) was calculated to have an indicative estimation of epifaunal density. A scheme for all the samples is reported in Table 2. The classification at the genus level was made according to the most used taxonomical study on foraminiferal genera [87], while species were determined according to some important studies on the Mediterranean area [25,88,89], and to the World Modern Foraminifera Database. Neoconorbina spp. includes mainly N. posidonicola and little specimens (possible juvenile stage) very difficult to distinguish at specific level. Frequent agglutinated subcircular encrusted tests attached to phytal surface (rhizomes or leaves) were grouped under Daitrona sp. [88].
The foraminifers were hand-picked and photographed using FEIQUANTA400 scanning electron microscope at the SEM Laboratory of Earth Sciences Department, Sapienza University of Rome.

General Characters of Sediment and Posidonia Samples
All sediment samples consist mainly of coarse sandy volcanic grains and are characterized by the dominance of plagioclase (andesine) and amphibole (pargasite), which are the most common phenocryst phases in the andesitic and dacitic lava widely outcropping in Panarea and in the islets [40] (Figures 2 and 3).    Additional phases recognized in the sediments are dolomite, which is present in nor amounts in almost all samples and an argillitic component (illite) that occurs mai in the hot sample HFT (Figure 3).
The biogenic fraction of the sediment floor is characterized by scarce or absent ve tal debris, mollusks (bivalves and gastropods), bryozoans, serpulids, diatoms, echino and ostracods in addition to foraminifers. Brachiopods, Joania cordata and Argyrotheca neata, two common species in the Mediterranean Sea [91] are recorded mainly in the r zomes at the cold and Bottaro stations; in the control samples, rare individuals are fou both in the sediment and attached on the Posidonia rhizomes.
Posidonia leaves and rhizomes are characterized by the presence of encrusting org isms such as serpulids, coralline algae, bryozoans, sponges and diatoms. They are v scarce or totally absent at the P21 station leaves, except for diatoms that are always ab dant and highly diversified associated with siliceous sponge spicules.

Foraminiferal Content: General Features
Sediment samples-Dead assemblage is totally or almost absent in most of the sa ples. A total of 40 foraminiferal specimens were found in the sediment samples wit reduction in the number of species and individuals along a pH gradient from 8 to 6 (Ta 3a). The most acidic conditions (pH < 6) were characterized by rare and exclusively agg tinated taxa.
Posidonia samples-A total of 1125 individuals were counted on the total Posido leaves and rhizomes. The number of individuals tends to increase along the pH gradi Additional phases recognized in the sediments are dolomite, which is present in minor amounts in almost all samples and an argillitic component (illite) that occurs mainly in the hot sample HFT (Figure 3).
The biogenic fraction of the sediment floor is characterized by scarce or absent vegetal debris, mollusks (bivalves and gastropods), bryozoans, serpulids, diatoms, echinoids and ostracods in addition to foraminifers. Brachiopods, Joania cordata and Argyrotheca cuneata, two common species in the Mediterranean Sea [91] are recorded mainly in the rhizomes at the cold and Bottaro stations; in the control samples, rare individuals are found both in the sediment and attached on the Posidonia rhizomes.
Posidonia leaves and rhizomes are characterized by the presence of encrusting organisms such as serpulids, coralline algae, bryozoans, sponges and diatoms. They are very scarce or totally absent at the P21 station leaves, except for diatoms that are always abundant and highly diversified associated with siliceous sponge spicules.

Foraminiferal Content: General Features
Sediment samples-Dead assemblage is totally or almost absent in most of the samples. A total of 40 foraminiferal specimens were found in the sediment samples with a reduction in the number of species and individuals along a pH gradient from 8 to 6 (Table 3a). The most acidic conditions (pH < 6) were characterized by rare and exclusively agglutinated taxa.
Posidonia samples-A total of 1125 individuals were counted on the total Posidonia leaves and rhizomes. The number of individuals tends to increase along the pH gradient (Table 3b,c) except for the hot site, where an increase of density is recorded both in the leaves and in the rhizomal microhabitats. Generally, the highest density is concentrated on leaves particularly on the adult portion of them. Only in the most acidic areas (P21, BP and B2 sites) rhizomal density (mainly agglutinated taxa, Table 3b) is similar or higher than leaves one. The biodiversity is low in all sites. The greatest number of taxa (S) is recorded at the control site B1 (Table 3c). Diversity indices display α-Fisher index < 6 except for the Cold site where it reaches 6.45 value, and Shannon index (H) always <2. Agglutinated tests are clearly dominant in more acidic sites, while porcelaneous and hyaline shells reach their highest frequencies in normal marine conditions (control sites: Raya and B1). However, it is noteworthy that the hyaline taxa abundance and diversity increase on both Posidonia leaves and rhizomes at the Hot spot sites (Table 3a).   The most abundant species is Rosalina bradyi (25.4% of epiphytic total) followed by Neoconorbina spp. (18.9%) and Planorbulina mediterranensis (9.8%). All these hyaline species accounted the highest frequencies at Hot spots stations while decrease until to be totally absent at P21 and Black Point sites. The list of total species is reported in the Table S1. In all sites with pH values < 8, the SEM analyses highlight strong signs of decalcifications and alterations (variations in the pore distribution, cracks and fractures) resulting an increase of shell fragility for the most part of the tests (Figures 4 and 5 The most abundant species is Rosalina bradyi (25.4% of epiphytic total) followed by Neoconorbina spp. (18.9%) and Planorbulina mediterranensis (9.8%). All these hyaline species accounted the highest frequencies at Hot spots stations while decrease until to be totally absent at P21 and Black Point sites. The list of total species is reported in the Table S1. In all sites with pH values < 8, the SEM analyses highlight strong signs of decalcifications and alterations (variations in the pore distribution, cracks and fractures) resulting an increase of shell fragility for the most part of the tests (Figures 4 and 5).
The most abundant species is Rosalina bradyi (25.4% of epiphytic total) followed by Neoconorbina spp. (18.9%) and Planorbulina mediterranensis (9.8%). All these hyaline species accounted the highest frequencies at Hot spots stations while decrease until to be totally absent at P21 and Black Point sites. The list of total species is reported in the Table S1. In all sites with pH values < 8, the SEM analyses highlight strong signs of decalcifications and alterations (variations in the pore distribution, cracks and fractures) resulting an increase of shell fragility for the most part of the tests (Figures 4 and 5).   (1,2) Cold site at 7 < pH < 8, similarly to Rosalina specimens, the test shows morphological variations mainly in the pore size and distribution on the surface; (3)(4)(5) Hot site at acidic conditions 6 < pH < 7.9 with high temperature, signs of decalcification (3,5), increase of pore size and morphological abnormalities are visible (4).

Ditella Station Hot Samples
Sediment samples-The foraminiferal assemblage is constituted exclusively of rare specimens of Lepidodeuterammina ochracea. Only diatoms and partly dissolved carbonate tests (serpulids, ostracods and foraminifers) are found in bottom samples.
Posidonia leaves-A total of 292 living foraminifers were found on 11 leaves (2 leaves were totally barren in foraminifers) with a F/P ratio of 24.43. The most part of specimens was attached in the mature part of the leaves, while only 9 specimens (Lobatula lobatula, Asterigerinata mamilla, P. mediterranensis, R. bradyi) were found also on the juvenile portion of them. The foraminiferal community is represented mainly by Neoconorbina spp. (40%) and R. bradyi (33%) which are clearly the most abundant taxa, P. mediterranensis follows with values of 15%, while the remaining species are <4% ( Figure 6).
With reference to the shell composition, the main part of foraminifers belongs to hyaline group (99%), while agglutinated and porcelaneous (miliolids) groups are about 1% (Figure 7).
Posidonia rhizomes-A total of 50 specimens were recorded on the rhizomes. Rosalina brady (38%) and Neoconorbina spp. (20%) dominate the assemblage followed by miliolids with 18%. The remaining taxa display values < 10% except for Spirillina vivipara which is recorded with frequencies of 14%. The agglutinated taxa are totally absent.
With reference to the shell composition, the main part of foraminifers belongs to hyaline group (99%), while agglutinated and porcelaneous (miliolids) groups are about 1% (Figure 7).
Posidonia rhizomes-A total of 50 specimens were recorded on the rhizomes. Rosalina brady (38%) and Neoconorbina spp. (20%) dominate the assemblage followed by miliolids with 18%. The remaining taxa display values <10% except for Spirillina vivipara which is recorded with frequencies of 14%. The agglutinated taxa are totally absent.    Figure 6).
Posidonia rhizomes-B1 site shows 96 specimens, among them the most abundant are R. bradyi (38%) followed by miliolids (21%) that slightly increase, while Neoconorbina spp. (8%) decrease with respect to the leaf's assemblages. In B2 rhizomes foraminiferal content displays 45 specimens among them the most abundant species is R. bradyi (69%). Porcelaneous are present with percentages of 17% while other species are very scarce with values <3% (Figures 6 and 7). In both sites agglutinated taxa are present with low percentages <3%.

P21 and Black Point Site
Posidonia leaves-At the P21 site a total of 39 leaves was analyzed. The F/P ratio records its lowest value (0.18). The most part of them is barren in foraminiferal content; only 5 specimens of Daitrona sp., an agglutinated taxon, and one specimen of R. bradyi were recorded. At the Black Point site Posidonia leaves developed no closer than 10 m from emission points and exhibit very scarce foraminiferal content. The F/P is 2.5 and on 12 leaves analyzed, 30 specimens were found. They are represented mainly by agglutinated taxa such as A. globigeriniformis (33%), Eggerelloides scaber (20%), A. planorbis (7%), Trochammina inflata (7%) and Daitrona sp. (13%).

Summary of Different pH Conditions
In Panarea sites, our observations allowed us to detect three main different pH conditions with respect to normal marine conditions with pH = 8.1 (control Raya and B1 sites): low acidic conditions (pH = 7.9), acidic conditions (Hot sites; pH=7.0-5.6) and strong acidic conditions (Black Point and P21 sites; pH = 3.8-5.2) The normal marine conditions display assemblages typical of shallow water oligotrophic environment. On rhizomes, the assemblages are mainly characterized by relatively high frequencies of miliolids and R. brady, whereas on leaves L. lobatula and P. mediterranensis dominate. No deformities or abnormal test were observed. In the low acidic conditions, rosalinids increase significantly with R. bradyi mainly concentrated on the rhizomes, whereas Neoconorbina spp. are more frequent on sediments and leaves. In acidic conditions Rosalinids continue to dominate the epifaunal assemblages, while under strong acidic conditions, moving from the spot emission at not less than 10 m (5.2 < pH < 7), the first Posidonia shoots host only a typical agglutinated assemblage both on leaves and rhizomes. In addition, the rare miliolids recorded in the rhizomes, are covered by an agglutinated protective layer. Frequent deformities or abnormal test were evidenced starting from pH values < 7.

Response of Foraminiferal Assemblages to CO 2 Fluid Emission
The abiotic parameter pH is well known to influence the distribution, biodiversity and abundance of benthic foraminifers and, in case its value is far from the optimum, it acts as environmental stressor [22,26,31]. The microfaunal analysis conducted in sites under different pH and T conditions around the Panarea Island, allow us to make some considerations. The benthic foraminiferal response to increasing water acidification consists of reduction of biomineralization potentiality, loss of biodiversity and morphological modifications. In addition, environmental contexts may affect differently these responses to future OA [28,92,93]. It can be very variable on the basis of complex environmental interactions. Experimental studies showed as the exposure times and the temporal variations of the acidic conditions are among the most significant factors [28]. In the Panarea sites, the available geophysical data show a continuous hydrothermal acidic activity affecting the benthic ecosystem [55,69] with a pH decrease ranging from 0.2 to 4.3 units from natural marine conditions (pH = 8).
Differently from controlled laboratory conditions, the great variability of physical and chemical parameters of the Panarea sites makes very difficult to define a pattern of biota responses in situ. This is probably the reason of the observed different response of benthic fauna (e.g., Cold and B2 sites) under similar pH values. From this perspective, the obtained data should be considered as the results of lasting foraminiferal adaptation in response to the overall environmental conditions rather than the immediate response to the values measured at the sampling time.
Our observations show that the response of the Panarea foraminiferal assemblages displays different patterns in sediment and Posidonia (rhizomes and leaves) samples. Regarding the sediment assemblages, the scarcity of the microfaunal content (dead and living) and the lack of data from all analyzed sites, does not allow a reliable assessment of the foraminiferal response to the different acidic conditions. Based on our data and likewise to Pontine vent activity (central Tyrrhenian Sea [27]), strong and persistent CO 2 emissions probably cause the complete dissolution of post-mortem calcareous tests at the sea bottom and the inhibition of the living assemblage. On the other hand, the calcareous component in the sediment is very low except for the presence of rare dolomite (low-temperature dolomite) that may be precipitate from the hydrothermal fluids which permeate the host rock ( Figure 3) [94,95].
On the contrary, Posidonia epiphytic assemblages show a more complete picture. The distinction between leaves and/or rhizomes microhabitats allow us to carry out considerations not only on the microhabitat preference of foraminifers at the species level but also to highlight their functional aspect. The highest density is concentrated on leaves particularly on the adult portion of them, confirming the typical pattern of epiphytic colonization under normal marine condition [21]. Leaves represent "high" substrates (e.g., [96][97][98]) on which suspension feeders enjoy of a better exploitation of nutrients in the water mass around the plant, a greater degree of oxygenation and probably a better advantage of the Posidonia buffer effect [21,99,100]. Literature data demonstrated that effect of the phytal metabolism related to Posidonia activity may favor the development of foraminifers under acidic conditions [101][102][103][104][105]. Both controlled experiments conducted in laboratories and in situ studies, showed how little and fluctuating pH decrease (<0.4 unit) can be modulated by Posidonia oceanica conferring to the environment a major pH stability [28].
Under significant pH decrease (P21, BP and B2 sites), rhizomal foraminiferal density (mainly agglutinated taxa) become higher than foraminiferal density on leaves. These observations compel us to consider rhizomes as "refugia" with respect to more acidic surrounding environments. Rhizomes can trap sediment and with-it nutrients, and they can host fungi and bacteria providing food sources for many organisms including foraminifers [88,[105][106][107].
No taxon at the species level is exclusive of the two Posidonia microhabitats (leaves and rhizomes). Rosalinids (R. bradyi and Neoconorbina spp.) are abundant in both microhabitats although Neoconorbina spp. show a preference for leaves. They are temporarily attached foraminifers with a great ability for detachment in response to changing environmental conditions moving from epifaunal (substrate and seagrass) to sediment microhabitats (or conversely) to find more advantageous life conditions [24,88,[108][109][110][111][112]. Probably the motility and relatively short life span of rosalinids make these taxa more easily adapted to stress conditions, showing a typical opportunistic behavior (B morphotype after [21,24]).
On leaves rosalinids are frequently associated with cibicidids (L. lobatula, C. refulgens) which have similar ecological features. Although they are characterized by an attachment surface, they are predominantly motile species. These taxa are always associated with P. mediterranensis that, on the contrary, is considered a permanently attached species with long life spans and needs of optimal environmental conditions (A morphotype after [21,24,[113][114][115]). This is testified by the high frequencies of this species in the control site (Raya site) and its clear decrease under acidic condition (except for warm hot site) indicating its preference of normal marine conditions and low tolerance towards stressed environment. The presence of P. mediterranensis in Panarea acidic sites may be explained by the mitigating effect of Posidonia meadow. The photosynthetic activity of the plant ensures the survival conditions for P. mediterranensis although evident test abnormalities and dissolution signs testify non-optimal conditions. In the rhizomes, rosalinids are associated with miliolids and small agglutinated species that are permanently motile and have a very short life cycle (<3 months, D morphotype after [21,24], well adapted to stressed environment. Moreover, it is noted that many agglutinated species (mainly A. globigeriniformis, L. ochracea, A. planorbis) are characterized by a very small (about few microns), agglutinated grains notwithstanding the coarse grain size of the sediments. This highlights their selective capability to exploit the scarce fine fraction deriving from the alteration of the volcanoclastic material. Under acidic/neutral conditions these alteration processes may result in the formation of clay minerals such as kaolinite and illite/smectite, respectively ( Figure 3) [94,95]. Detailed considerations regarding different pH conditions follow.

Normal Marine Conditions (pH > 8)
The normal marine conditions sites (B1 and Raya) are characterized by epiphytic and epifaunal assemblages typical of shallow water oligotrophic environments [22,116]. In the sediment, although the sites are not affected by venting activity, the recorded faunal density is low, suggesting a possible indirect influence of CO 2 emissions (testified by frequent agglutinated taxa, Figure 6), or other adverse factors (e.g., scarce food availability on and into the sediments, the nature of the bottom or water current), that may negatively influence the foraminiferal development in these areas. This is evidenced also by the low diversity indices of the epifaunal assemblages (leaves and rhizomes) that, although more diversified than sediment assemblage, is more similar to marginal environments rather than typical shelf ecosystems [22] Our observations carried out from low acidic conditions confirm literature data with slight changes in the foraminiferal assemblages [26]. The two sites associated with low acidic conditions, Cold and Bottaro B2, show mainly differences in terms of quantitative data. However, at both sites, the transition from normal marine (pH = 8.1) to low acidic conditions (decreasing of 0.2 unit) is marked by a general decrease in Posidonia foraminiferal content and diversity (mainly highlighted on leaves), although the carbonate tests are also well represented. A first sign of stressed conditions can be inferred by the significant reduction of the species with limited adaptative capability such as P. mediterranensis on leaves. Similarly, on the rhizomes, the porcelaneous group decreases significantly. Literature data suggest that their high Mg content make them more susceptible to dissolution in acidic conditions and to undergoing tests deformities in response to seawater pollution [26,[117][118][119]. On the contrary, the sediment community (Cold site) displays an increase in terms of both density and diversity (H > 2 and Fisher indices > 9) proportionally for each test typology group with respect to the scarce foraminiferal content at the control site. As above mentioned, based on the available data, it is difficult to give an exhaustive explanation of that. However, a possible cause may be found in the more protected conditions of the sites that favor trapping of organic matter, a food source for foraminifers. The transition to acidic conditions is marked by a clear drop of foraminiferal content in the sediments with rare agglutinated taxa (L. ochracea). This finding confirms literature data indicating critical threshold pH values around 7.8 and 7.6 that prevent the building of foraminiferal carbonate test [26,31].
Regarding epiphytic foraminiferal content, contrary to what is expected, it accounts for the highest density values of all studied sites (high values of F/P ratio) and highly diversified assemblages (similar to control site and higher than low acidic conditions) probably explained by the buffer effect of Posidonia. This is testified by the increase of P. mediterranensis on leaves, and the re-occurrence of miliolids in the rhizomes although the most part of the tests recorded show clear signs of carbonate dissolution in terms of fragility, reductions of number of pores and increasing size of the single pores (L. lobatula, Rosalina bradyi, and P. mediterranensis, Figures 4 and 5).
In addition, the presence of warm water due to the proximity of Hot sites in this zone seems to enhance the foraminiferal development amplifying the beneficial effect of vegetable activity. Temperature is an important factor controlling not only geographical distributions but also metabolism, growth rates and reproductive processes of foraminifers [120], however little is known about the relationship between a single species and this parameter. Literature data and laboratory experiments demonstrated that 45 • C is the lethal value threshold for most part of macroforaminifers [22,121]. In the analyzed samples, among cibicidids, L. lobatula shows a higher tolerance to increasing temperatures than C. refulgens that is more frequent at cold sites.
Rosalinids continue to dominate the assemblage showing their opportunistic behaviour. Rosalina species are considered to have a high tolerance to elevated temperatures resulting in greater adaptability to the future warming [122]. Its strong adaptative behavior is also demonstrated under a wide range of stressed conditions such as high euthrophic environments [123], pH fluctuations [28] or other extreme environments [124]. Little is known about the ecological characteristics of Neoconorbina species. Based on our observations their behavior is very similar to that of Rosalina spp. According to literature data these conditions (low pH and very high temperatures up to 133 • C) do not allow any carbonate biomineralization process and prevent the development of carbonate foraminiferal tests [26,31,123] even in the presence of Posidonia seagrass.
At the species level, the association of A. globigeriniformis and E. scaber seems to be more resistant in strong acidic conditions confirming their opportunistic behavior in extreme environments such as submarine caves [123], other venting areas [27] or oil spill regions [124]. In addition, specimens attributable to Daitrona cf. Daitrona sp., a permanently attached taxon [21], seem to show tolerance to acidic conditions, because they agglutinate sediment particles to a proteinaceous or non-calcite matrix [125], a potential advantage under acidic conditions (Figure 8). nently attached taxon [21], seem to show tolerance to acidic conditions, because they agglutinate sediment particles to a proteinaceous or non-calcite matrix [125], a potential advantage under acidic conditions ( Figure 8).

Conclusions
The microfaunal analysis conducted at sites affected by CO2 venting activity around the Panarea Island, allow us to make some considerations on benthic foraminiferal response to increasing ocean acidification (OA). Under different pH conditions (from strong to low acidic conditions) the foraminiferal assemblage's response can be very different depending on several environmental conditions and microhabitats. A decrease of pH from the normal marine values ranging between 2.9 and 1.1 units (strong acidic conditions) determines a strong decline of faunal density and biodiversity with the presence of only rare agglutinated taxa.
Under acidic conditions, although the sediment niches proved prohibitive for foraminiferal life, Posidonia leaves and rhizomes can offer a better condition for the epiphytic foraminiferal fauna that can take advantage from the buffer effect of Posidonia seagrass.

Conclusions
The microfaunal analysis conducted at sites affected by CO 2 venting activity around the Panarea Island, allow us to make some considerations on benthic foraminiferal response to increasing ocean acidification (OA). Under different pH conditions (from strong to low acidic conditions) the foraminiferal assemblage's response can be very different depending on several environmental conditions and microhabitats. A decrease of pH from the normal marine values ranging between 2.9 and 1.1 units (strong acidic conditions) determines a strong decline of faunal density and biodiversity with the presence of only rare agglutinated taxa.
Under acidic conditions, although the sediment niches proved prohibitive for foraminiferal life, Posidonia leaves and rhizomes can offer a better condition for the epiphytic foraminiferal fauna that can take advantage from the buffer effect of Posidonia seagrass. Warm water conditions seem to amplify the beneficial effect of the Posidonia. However evident morphological abnormalities with reduction of biomineralization potentiality affect the foraminifers' tests.
The threshold of pH = 7.9 determines the restoration of the foraminiferal community in the sediments and a more abundant and diversified epiphytic assemblages. Among the foraminiferal assemblages rosalinids and agglutinated taxa represent the most abundant taxa and the more resilient species capable of facing future OA scenarios.
Similar response could be expected by all epibiont calcifying community (mollusks, brachiopods but also weakly calcified such as crustaceans) that could be less disadvantaged from the effect of OA thanks to the photosynthesis processes and the consequent buffering role of seagrass. This study confirms the role of Posidonia oceanica as refuge and in the mitigation contrasting the future climate change.