Introduction

Benthic foraminifers (Rhizaria) and ostracods (Crustacea) are meiofaunal groups generally provided with calcareous tests and valves, commonly preserved in sea bottom sediments. The composition of their assemblages reflects environmental conditions due to both natural and human causes. Anthropogenic activities produce various effects on shallow marine waters, including organic pollution, changes in sedimentation rates, increase in hydrocarbon and heavy metal concentrations, and eutrophication-induced hypoxia (Gooday et al., 2009; Yasuhara et al., 2012; Wilkinson et al., 2014) that, in turn, lead to an increase in the relative abundance of stress-tolerant foraminiferal species (Hayward et al., 2004; Alve et al., 2009; Frontalini and Coccioni, 2011; Ruiz et al., 2012) and frequently to ostracod diversity decrease (Alve, 1991; Mazzola et al., 1999; Cronin & Vann, 2003; Irizuki et al., 2018). The studies combining the analyses of benthic foraminifers and ostracods in areas where human-induced ecological variations occurred showed the high potential of calcareous meiofaunal assemblages as water quality indicators (Samir, 2000; Triantaphyllou et al., 2003, 2005; Vilela et al., 2003; Bergin et al., 2006; Pascual et al., 2008, Salvi et al., 2015).

In the present study, benthic foraminiferal and ostracod assemblages were studied from eleven samples collected in the infralittoral zone of the Gulf of Pozzuoli, a bay located in the Campania region (Southern Italy) with a narrow continental shelf, a shelf break at about 40 m bsl, a maximum depth of 110 m, and an average depth of ca. 60 m (Fig. 1; Somma et al., 2016). The gulf is mainly exposed to winds and sea waves approaching from the southeast–southwest sector, with a maximum geographic fetch of 665 km for the 205° direction, 0.9–2.2 m average wave height, and a maximum wave height of 4.7 m in winter (De Pippo et al., 2008). The water circulation models (De Maio et al., 1985; Menna et al., 2008; de Ruggiero, 2016) of the gulf generally indicate two main flow patterns: (i) when the open sea currents flow toward the southeast, the inner waters of the bay are cut off in a slow cyclonic gyre; then, the coastal waters slow motion could favor turbidity and a high pollutant concentration; (ii) when the open sea currents flow toward the northwest, some branches enter into the bay; then, a fair renewal of sea waters occurs. Tides are negligible, with a syzygial tide amplitude of 0.35 m (Tammaro et al., 2021); therefore, the gulf is a wave-dominated environment. Salinity, turbidity, and phytoplankton distribution are related to seasonal variation in the sea surface and column temperature, autumn–winter freshwater supply by rainfalls and land runoff, marine currents cell circulation: 37.1–38.6‰ salinity, 27–30° C sea surface temperature, and high phytoplankton biomass (Chl a concentration > 2 µg/L) were registered during the spring–summer season (Bolinesi et al., 2020). The bay was exposed to prolonged anthropogenic disturbance, due to urban and industrial wastes, at least since 1885, when an armaments factory was built by the British company Armstrong Mitchell & Co. The eastern part of the gulf was under the influence of the Bagnoli steel plant from 1910 to 1990, to which the high levels of polycyclic aromatic hydrocarbons (Arienzo et al., 2017; Ferrara et al., 2020), trace metals (Trifuoggi et al., 2017), and rare earth elements (Trifuoggi et al., 2018) in the sediment seem to be linked. Recently, some ecological and paleoecological investigations were performed on Recent (Balassone et al., 2016; Mangoni et al., 2016; Arienzo et al., 2020) and late Quaternary (i.e., from ~ 150 ka to historical times; Aiello et al., 2012; 2018; 2020; 2021; Amato et al., 2019; Petrosino et al., 2021) sediments of the Campania region coastal areas focused on benthic foraminiferal and ostracod assemblages. The present study aims to define the characteristics of the calcareous meiofaunal assemblages in the infralittoral zone (Peres & Picard, 1964; Peres, 1982) of an area showing high geoaccumulation values and to test a possible decrease in benthic faunal abundance and diversity in polluted bottom sediments. Our data, compared with the above-mentioned studies and previous investigations on Campanian infralittoral benthic foraminifers (Moncharmont Zei, 1964; Sgarrella & Barra, 1985; Sgarrella et al., 1985; Sgarrella & Moncharmont, 1993) and ostracods (Müller, 1894; Puri et al., 1964, 1969), may contribute to a more complete understanding of the relationship between meiofaunal assemblages and environmental parameters.

Fig. 1
figure 1

Location of the sampling stations (black solid circle) in the Gulf of Pozzuoli. Legend: 1, pyroclastics of the Phlegrean Fields (Late Pleistocene–Holocene); 2, deposits of transitional environments (Quaternary); 3, isobath(-m); 4, edge of the continental shelf break, from a depth of about 25 down to 40 m. Depth is in meters b.s.l. (after Somma et al., 2016 and the morphobathymetric and sedimentological surveys carried out for this research). The geographic coordinate system is WGS84

Material and methods

Eleven samples of very fine to coarse sands, and very fine gravels, were collected by a Van Veen grab above the shelf break of the Gulf of Pozzuoli (~ 40 m bsl), in a water depth range between 7.5 m and 38 m, within the infralittoral zone. The grab collected superficial sediments, including the first ca. 5 cm of the seabed, related to texture and consistency of silt or sand deposits. Sharp & Nardi (1987) calculated, in this area, a sedimentation rate of about 4 mm/year and consequently the sampled sediments deposited, at most, in the last 15 years. The surface of the sampler is about 150 square cm (10 × 15 cm), while the volume of the sampled sea bottom surface sediment generally is about 750–1000 cubic cm. Samples were taken along one campaign in the spring of 2017 along transects and aboard a motor vessel. Bathymetry, grain size, number, and location of samples are reported in Fig. 1 and Table 1. For meiofaunal analyses, all the samples were oven-dried, and 100 g of dry sediment was taken. They were washed through 230-mesh (63 μm) and 120-mesh (125 μm) sieves, and the residues were oven-dried and examined with a reflected light microscope. A microsplitter was used to obtain subsamples when necessary. About 300 benthic foraminiferal tests and 300 ostracod valves were picked from the coarsest fraction (> 125 μm), classified, and counted. Abundance and diversity indices were calculated using the number of foraminiferal specimens, the ostracod Minimum Number of Individuals (MNI), and the Total Number of Valves (TNV). MNI is the greater number between right and left adult valves plus the number of adult carapace; when only juvenile shells (j) were recorded, the MNI equals one. TNV includes all the adult and young instar valves. Assemblage composition as well as diversity indices was considered for environmental discussion. The following indices were calculated: S (taxa richness), I (individuals per 100 g of sediment), D (dominance), H' (Shannon’s diversity index, using natural logarithm), and J (equitability). The species were identified according to classic and modern literature both for benthic foraminifers and ostracods (Aiello & Barra, 2010; Aiello et al., 2018, and references therein). The studied specimens are housed in the Aiello Barra Micropaleontological Collection (A.B.M.C.), Dipartimento di Scienze della Terra, dell’Ambiente e delle Risorse, Università degli Studi di Napoli Federico II. Statistical analyses were performed using abundance values of foraminiferal (I = number of individuals per 100 g) and ostracod (both MNI = minimum number of individuals per 100 g, and TNV = total number of valves per 100 g) assemblages. Q-mode cluster analysis (paired group as an algorithm, Rho as similarity measure) was performed to determine groups of samples with similar meiofaunal composition, using abundance values of all foraminiferal and ostracod (both MNI and TNV) species. Pearson’s correlation coefficient was used to test for correlation between assemblage features, depth, major and trace elements, and polycyclic aromatic hydrocarbons of eight fine-grained samples; benthic foraminiferal and ostracod species with relative abundance greater than 5% in at least two samples were considered. The abiotic variables were subject to z-standardization. Analyses were carried out on the same set of samples used by Arienzo et al. (2017), Trifuoggi et al. (2017, 2018), and Ferrara et al. (2020) in their investigations on the distribution of polycyclic aromatic hydrocarbons (PAHs), trace metals (HMs), and rare earth elements (REEs). All the analytical determinations were performed in triplicate for each sample taken at each site. The quality of the analytical results is assured by participation in ring tests for the determination of HMs, PAHs, and REEs from sediments and similar matrices. Mean recoveries ranged from a minimum of 85% to a maximum of 97%. Grain size analyses were performed following the standard methodology of Folk and Ward (1957). Full methodological details on sampling techniques, geochemical and grain size analyses were reported in Arienzo et al. (2017), Trifuoggi et al. (2017, 2018), and Ferrara et al. (2020).

Table 1 Coordinates of sampling stations, grain size, and water depth

Computation of diversity indices and statistical analysis were performed with STATISTICA 5 (StatSoft Inc., Tulsa, OK, USA).

Results

All the samples yielded both benthic foraminiferal and ostracod shells (no barren samples) (Tables 28). A total of 4262 foraminiferal individuals and 3607 ostracod valves were collected. The good state of preservation, the distribution data, and the presence of all developmental stages (in ostracods: different young instars and adults) suggested that the calcareous meiofaunal assemblages could be considered entirely autochthonous. The benthic foraminiferal assemblages included 142 species assigned to 74 genera; 127 ostracod species in 49 genera were recorded (Appendix 1; Figs. 23). Five benthic foraminiferal species and eight ostracod species were tentatively identified or left in open nomenclature and nine with affinitive status due to the absence of adult specimens, or because of poorly preserved shells. The good state of preservation, the distribution data, and the presence of all developmental stages (in ostracods: different young instars and adults) suggested that the calcareous meiofaunal assemblages could be considered entirely autochthonous.

Table 2 Benthic foraminiferal absolute abundance (I = individuals per 100 g of sediment)
Table 3 Benthic foraminiferal relative abundance (RA, %) samples
Table 4 Ostracod absolute abundance [I(MNI) = minimal number of individuals per 100 g of sediment]; j indicates juvenile specimens
Table 5 Ostracod relative abundance (RA, %) (MNI = Minimal Number of Individuals)
Table 6 Ostracod absolute abundance [I(TNV) = Total Number of Valves per 100 g of sediment]
Table 7 Ostracod relative abundance (RA, %) (TNV = Total Number of Valves)
Table 8 Benthic foraminiferal and ostracod assemblage indices
Fig. 2
figure 2

1 Eggerelloides scaber (Williamson, 1858), lateral view, sample TP2—2, ABMC 2019/042 2 Quinqueloculina lata Terquem, 1876, four chamber side, sample TP1—1, ABMC 2019/054 3 Quinqueloculina seminulum (Linnaeus, 1758), four chamber side, sample TP2—2, ABMC 2019/045 4 Quinqueloculina stelligera Schlumberger, 1893, peripheral view, sample TP1—2, ABMC 2019/044 5 Siphonaperta aspera (d’Orbigny, 1826), side view, sample TP2—1, ABMC 2019/040 6 Cycloforina contorta (d´Orbigny, 1846), four chamber side, sample TP1—2, ABMC 2019/058 7 Cycloforina contorta (d´Orbigny, 1846), peripheral view, sample TP1—2, ABMC 2019/059 8 Triloculina trigonula (Lamarck, 1804), peripheral view, sample TP1—2, ABMC 2019/046 9 Miliolinella semicostata (Wiesner, 1923), side view, sample TP1—1, ABMC 2019/048 10 Bulimina elongata d’Orbigny, 1846, lateral view, sample TP3—2, ABMC 2019/047 11 Rosalina macropora (Hofker, 1951), spiral side, sample TP1—2, ABMC 2019/055 12 Rosalina macropora (Hofker, 1951), umbilical side, sample TP1—2, ABMC 2019/056 13 Asterigerinata mamilla (Williamsom, 1858), spiral side, sample TP1—2, ABMC 2019/051 14 Ammonia aberdoveyensis Haynes, 1973, spiral side, sample TP1—1, ABMC 2019/039 15 Ammonia beccarii (Linnaeus, 1758), spiral side, sample TP2—3, ABMC 2019/060 16 Ammonia beccarii (Linnaeus, 1758), umbilical side, sample TP2—3, ABMC 2019/061 17 Buccella granulata (Di Napoli Alliata, 1952), spiral side, sample TP2—3, ABMC 2019/043 18 Tretomphalus concinnus (Brady, 1884), spiral side, sample TP1—1, ABMC 2019/052 19 Cibicides lobatulus (Walker & Jacob, 1798), umbilical side, sample TP1—2, ABMC 2019/053 20 Haynesina depressula (Walker & Jacob, 1798), side view, sample TP1—1, ABMC 2019/050 21 Elphidium crispum (Linnaeus, 1758), side view, sample TP2—3, ABMC 2019/041 22 Elphidium maioricense Colom, 1942, side view, sample TP3—2, ABMC 2019/049 23 Planorbulina mediterranensis d’Orbigny, 1826, unattached side, sample TP3—2, ABMC 2019/057 Scale bar 100 µm

Fig. 3
figure 3

1 Aurila convexa (Baird, 1850), left valve, sample TP2—3, ABMC 2019/065 2 Callistocythere lobiancoi (Müller, 1894), left valve, sample TP1—2, ABMC 2019/078 3Callistocythere flavidofusca (Ruggieri, 1950), right valve, sample TP2—3, ABMC 2019/067 4 Urocythereis margaritifera (Müller, 1894), right valve, sample TP2—1, ABMC 2019/068 5 Cistacythereis turbida (Müller, 1894), left valve, sample TP2—3, ABMC 2019/074 6 Pontocythere turbida (Müller, 1894), right valve, sample TP2—2, ABMC 2019/072 7 Procytherideis retifera Ruggieri, 1978, right valve, sample TP1—1, ABMC 2019/081 8 Procytherideis retifera Ruggieri, 1978, carapace in dorsal view, sample TP1—1, ABMC 2019/082 9 Sagmatocythere napoliana, left valve, sample TP3—2, ABMC 2019/075 10 Paracytheridea triquetra (Reuss, 1850), left valve, sample TP1—1, ABMC 2019/066 11 Paracytheridea paulii Dubowsky, 1939, right valve, sample TP1—1, ABMC 2019/062 12 Semicytherura robusta Bonaduce, Ciampo & Masoli,1976, left valve, sample TP2—3, ABMC 2019/077 13 Semicytherura incongruens (Müller, 1894), right valve, sample TP2—2, ABMC 2019/070 14 Semicytherura rarecostata Bonaduce, Ciampo & Masoli, 1976, right valve, sample TP1—2, ABMC 2019/076 15 Semicytherura ruggierii (Pucci, 1955), right valve, sample TP1—1, ABMC 2019/080 16 Loxoconcha affinis (Brady, 1866), left valve, sample TP1—1, ABMC 2019/063 17 Loxoconcha rhomboidea (Fischer, 1855), left valve, sample TP1—1, ABMC 2019/069 18 Loxoconcha ovulata (Costa, 1853), left valve, sample TP2—2, ABMC 2019/071 19 Cytherois uffenordei Ruggieri, 1974, left valve, sample TP1—2, ABMC 2019/064 20 Xestoleberis communis Müller, 1894, left valve, sample TP1—2, ABMC 2019/073 21 Xestoleberis dispar Müller, 1894, right valve, sample TP2—2, ABMC 2019/079 Scale bar 100 µm

Benthic foraminifers.

Six benthic foraminiferal species were present in all the samples, that is, Ammonia aberdoveyensis, Buccella granulata, Cibicides lobatulus, Elphidium crispum, Quinqueloculina seminulum and Triloculina schreiberiana. Assemblages were characterized by the genera Quinqueloculina (19 species) and Elphidium (10 species). Cibicides lobatulus was the most abundant species, with a Medium Relative Abundance (MRA) of 12.8%, followed by Tretomphalus concinnus (MRA = 5.54%), Siphonaperta aspera (MRA = 4.30%), Elphidium crispum (MRA = 3.97%), Asterigerinata mamilla (MRA = 3.75%) and Q. seminulum (MRA = 3.53%).

The number of species (S) was between 17 and 82, with discrimination between the three coarse-grained samples (TP2-1 = coarse sand; TP3-1 = medium sand; TP5-1 = very fine gravel) and the remaining eight samples, all made of fine or very fine sands. The former displayed a S range from 17 to 48, with the mean value of 34.33; the eight fine-grained samples had S between 52 and 82 (mean value = 61.25).

The number of specimens (I) showed a wide range, from 69 to 105,984. The three coarse-grained samples displayed a mean value of 484, whereas a I mean value of 51,540 was recorded for the remaining samples.

The dominance (D) was between 0.04 and 0.20, with high values in the samples TP2-1 (D = 0.20) and TP5-1 (D = 0.16) and low (D(0.06) the remaining ones.

TP2-1 and TP5-1 assemblages showed a Shannon diversity index (H ') less than 3; in the other assemblages H' >3. Mean H' is 3.28. Mean H' is 3.28. A similar trend was observed for equitability (J), low in the samples TP2-1 and TP5-1 (J<0.8) and high (J>0.8) in the other samples.

Ostracods.

The most diversified genera were Semicytherura (18 species) and Xestoleberis (12 species). Characteristic species were Urocythereis margaritifera [MRA(MNI) = 8.10%; MRA(TNV) = 9.94%], Pontocythere turbida [MRA(MNI) = 6.00%; MRA(TNV) = 5.64%], Semicytherura rarecostata [MRA(MNI) = 5.56%; MRA(TNV) = 3.87%], Loxoconcha rhomboidea [MRA(MNI) = 5.08%; MRA(TNV) = 4.83%] and Loxoconcha ovulata [MRA(MNI) = 3.68%; MRA(TNV) = 5.55%]. Loxoconcha affinis [MRA(MNI) = 3.45%; MRA(TNV) = 3.64%], Xestoleberis dispar [MRA(MNI) = 2.93%; MRA(TNV) = 5.26%], Xestoleberis communis [MRA(MNI) = 2.25%; MRA(TNV) = 5.22%] and Aurila convexa [MRA(MNI) = 1.93%; MRA(TNV) = 4.95%] were considered accessory species.

The three samples with coarser granulometry (TP2-1, TP3-1, TP5-1) yielded relatively poor ostracod assemblages, showing low diversity and high dominance. Conversely, in the fine-grained samples, diversity and abundance were high and the dominance low.

Taxa richness (S) ranged from 4 to 12 in the coarse-grained samples and from 31 to 57 in the remaining ones. In the former samples, abundance (I) was between 9 and 27 (MNI) and between 18 and 40 (TNV); in the latter samples, the mean value of I was 2065.5 (MNI) and 7781 (TNV). Shannon index H' followed a similar trend: in TP2-1, TP3-1, TP5-1 mean H' (MNI) was 1.68 and medium H' (TNV) was 1.56; in the assemblages occurring in the fine-grained sediments, medium H' (MNI) was 3.35 and medium H' (TNV) was 3.10.

Dominance (D) values were high in the assemblages of samples TP2-1, TP3-1, TP5-1 [D (MNI) range = 0.09–0.48; D(TNV) range = 0.14–0.70]; in the fine-grained samples the average D was 0.05 (MNI) and 0.07 (TNV).

Equitability (J) ranged from 0.72 to 0.99 (MNI) and from 0.46 to 0.94 (TNV). The mean J values were 0.86 (MNI) and 0.79 (TNV).

Statistics.

The cluster analysis (Fig. 4) revealed two clusters of samples, obtained at a similarity cut-off level of 0.45. Cluster B consists of the coarse-grained samples TP2-1, TP3-1, and TP5-1, with low diversity–low abundance assemblages; cluster A includes the remaining eight fine-grained (fine and very fine sands) samples, characterized by high-diversity—high abundance assemblages. The three sediment samples grouped in Cluster B, consisting of medium sand, coarse sand, and very fine gravel, showed low geochemical accumulation. Since both low meiofaunal abundance/diversity values and low pollutant concentrations are highly associated with grain size (v. Discussion section), a correlation analysis including all the samples would provide results strongly influenced by granulometry. Consequently, we opted for performing the Pearson’s correlation coefficient analysis on the eight fine-grained samples included in Cluster A. Results of Pearson’s correlation coefficient analysis using meiofaunal assemblages, depth, major and trace elements, total organic carbon, and polycyclic aromatic hydrocarbons (Table 9) are reported in Table 10. The foraminiferal species Cibicides lobatulus and Elphidium crispum are common in all the samples. The assemblages included in Cluster A are characterized by the foraminifers Tretomphalus concinnus, Asterigerinata mamilla, Triloculina trigonula, and Elphidium punctatum and by the ostracods Semicytherura rarecostata, Loxoconcha ovulata, L. rhomboidea, Aurila convexa, and Xestoleberis communis. In Cluster B, the foraminiferal species Siphonaperta aspera, Quinqueloculina seminulum, Elphidium macellum, and the ostracods Urocythereis margaritifera and Pontocythere turbida characterized the assemblages.

Fig. 4
figure 4

Dendrogram based on cluster analysis of benthic foraminiferal and ostracod (MNI and TNV) abundance data (I)

Our results show that the anthropogenic impact in the infralittoral zone of the Gulf of Pozzuoli, recorded in the geochemical accumulation (v. Table 9), was not reflected by diversity indices of the calcareous meiofaunal assemblages. The present findings were compared with the results of a previous study by Moncharmont (1964), based on a sampling carried out in 1961 by Harbans S. Puri and the Stazione Zoologica Anton Dohrn, where the characteristic species of the infralittoral zone were Ammonia beccarii, A. mamilla, C. lobatulus and T. concinnus (splitted by Moncharmont in T. concinnus and Rosalina globularis); conversely, Q. seminulum, S. aspera, Bulimina elongata were very rare. Our data suggested that Q. seminulum and B. elongata, stress-tolerant species (Aiello et al., 2018; Debenay et al., 2009), and S. aspera may have increased their abundance in the last decades.

Table 9 Concentrations (%) of TOC = total organic carbon; concentrations (mg/kg) of heavy metals, REE + Y = rare earth elements + Y, LREE (light rare earth elements), HREE (heavy rare earth elements), total PAHs (polycyclic aromatic hydrocarbons) and total PD PAHs (PD: Priority Dangerous PAHs: naphthalene (NAP), anthracene (ANT), benzo(b)fluoranthene (BbF), benzo(k)fluoranthene (BkF), benzo(ghi)perylene (BgP), benzo(a)pyrene (BaP); (Arienzo et al., 2017; Trifuoggi et al., 2017; 2018; Ferrara et al., 2020); concentrations of national regulatory guidelines (Ministero dell'Ambiente e della Tutela del Territorio, 2003; Ministero dell'Ambiente e della Tutela del Territorio e del Mare, 2009); background concentrations for the study area (Damiani et al., 1987)

The 1961 sampling campaign also provided bottom sediments of the adjacent Gulf of Naples. Their ostracofaunas were analyzed by Puri et al. (1964). In the infralittoral zone, the ostracod assemblages were dominated by Aurila, Urocythereis, Carinocythereis, Costa and Pontocythere (Cushmanidea in Puri et al., 1964) species, Tenedocythere prava (= Quadracythere (?) prava plus «Cythereis» polygonata, the latter name including the juveniles, in Puri et al., 1964) and S. incongruens. In the present infralittoral ostracofauna of the Gulf of Pozzuoli U. margaritifera and P. turbida are very common, but overall the assemblages are not dominated by trachyleberid species, being characterized by the genera Loxoconcha and Xestoleberis and by the species S. rarecostata.

It can be hypothesized that industrial and urban contaminants caused a meiofaunal change, more apparent in ostracod than in benthic foraminifers. Conversely, the high-diversity and low-dominance values showed that in oligotrophic and well-oxygenated waters (e.g., Hyams-Kaphzan et al., 2009, in a coastal area under the influence of sewage sludge with organic matter input), the diversity indices were not negatively influenced by the high concentration of pollutants in bottom sediments.

The sensitivity of ostracod and benthic foraminifers was also displayed by statistical analysis. Noteworthy correlations between ostracod dominance and sediment contaminants (polycyclic aromatic hydrocarbons, rare earth elements) encourage the use of calcareous meiofaunal assemblages as environmental bioindicators. The analysis confirmed the results of previous investigations, such as the tolerance of Q. lata to high heavy metal concentrations (Romano et al., 2009) testified by high correlations (≥ 0.76) with Ni, Pb, and Zn. Quinqueloculina lata was considered a pollution-tolerant species, by Elshanawany et al. (2011; 2018) and by Romano et al. (2013) in environment under strong anthropogenic pressure. The study of Mangoni et al. (2016) showed a correlation of this species with high concentrations of inorganic nutrients.

The strong anticorrelations (≤ -0.80) between both TNV and MNI abundances of the typical shallow marine (Aiello et al., 2016, 2018) species P. paulii and U. margaritifera, and water depth suggested the preferences of these taxa for the uppermost part of the upper infralittoral zone. From a paleoecological point of view, the refining of the distribution range of paleobathymetric indicators is a primary objective. The paleodepth estimates of sedimentary successions, or levels, using benthic foraminifers and ostracods, may contribute to the reconstruction of the dynamics of volcanic areas (Aiello et al., 2007, 2012; Marturano et al., 2009; 2011a; 2011b; 2013; 2018; Di Vito et al., 2016; Isaia et al., 2019). The present results showed that high abundance values of the genera Paracytheridea and Urocythereis, and specifically of P. paulii and U. margaritifera, are characteristic of upper shoreface environments in waters shallower than 23 m.

Table 10 Bivariate correlation with the Pearson’s correlation coefficient. TOC = total organic carbon; TTPAH = total polycyclic aromatic hydrocarbons; TTPD = total priority dangerous polycyclic aromatic hydrocarbons; REE + Y = rare earth elements + Y; LREE = light rare earth elements; HREE = heavy rare earth elements

Pearson’s correlation coefficient analysis of benthic foraminifers revealed that T. concinnus and taxa richness (S) correlated with As; foraminiferal dominance with Fe and Zn; Q. lata correlated with Ni, Pb, Zn, TOC and anticorrelated with water depth.

The Pearson’s correlation coefficient analysis of ostracods showed that Cr and Cu displayed high correlations with Carinocythereis whitei, Semicytherura incongruens (both MNI and TNV), L. rhomboidea (MNI), L. ovulata, and X. dispar and an inverse correlation with equitability J (TNV). Water depth correlated with C. turbida and taxa richness S (MNI and TNV), L. ovulata and Shannon Index H' (MNI), and anticorrelated with Paracytheridea paulii and U. margaritifera (MNI and TNV) and with Procytherideis retifera (MNI).

Total PAHs and total PD PAHs (Priority Dangerous PAHs) correlated with P. retifera (MNI), A. convexa and D and showed an inverse correlation with J (TNV). Rare earth elements correlated with A. convexa and D (TNV).

Ni correlated with P. paulii, P. retifera, and U. margaritifera (MNI); As showed correlation with L. affinis (MNI and TNV), Cytherois uffenordei, Paracytheridea triquetra (MNI); Hg anticorrelated with H' (TNV); TOC correlated with D(TNV).

Discussion

The study of shallow water samples of the Gulf of Pozzuoli allowed us to investigate the distribution of meiofaunal calcareous assemblages in an area subject to long-term industrial and urban pollution. All the sediments were collected within the infralittoral zone where previous researches displayed high geoaccumulation levels, especially in the eastern part of the bay (Arienzo et al., 2017; Trifuoggi et al., 2017, 2018). Anthropogenic and natural environmental pressure may influence both the meiofaunal features (abundance, diversity indices, dominance) and the taxonomic composition of foraminiferal and ostracod assemblages. In the sediments of the Gulf of Pozzuoli, their diversity was strongly related to the grain size of the bottom sediments, being low in the coarse-grained samples and high in fine and very fine sands. The relationship between granulometry and meiofaunal diversity in shallow marine waters was investigated by several researchers, who identified indicative trends. Pokorný (1978) stated, as a general rule, that coarse sediments (e.g., oolites and clean sands), can hold only a limited number of ostracod species, whereas on pelitic bottoms and mud-mixed sands the diversity of the assemblages increases. A number of studies corroborated this relationship, showing high ostracod diversity in fine-grained sands and low diversity in coarser sediments (Hazel, 1975; Aiello et al., 2006), higher diversity on silts or mud-mixed sands, and lower in clean sands (Puri 1966, 1971) and negative correlation with granulometry in Quaternary sandy successions (Aiello et al., 2020). On the other hand, on muddy bottoms, ostracod diversity frequently diminishes (Benson & Maddocks, 1964; Hong, 2016); consequently, it was suggested that the relationship between "ostracod diversity and particle size fractionation is not unimodal but rather hump shape" (Hong, 2016).

A similar trend was observed concerning benthic foraminifers. Investigations on assemblages collected on muddy bottom recorded low diversity, increasing in sandy sediments (Diz et al., 2004; Debenay et al., 2005); in large-grain sands and coarser bottom sediments a decrease was observed (Samir & El-Din, 2001; Temelkov, 2008; Delaine et al., 2015). It has to be noted that some studies reported apparently contradictory results, showing, for example, high foraminiferal diversity in sheltered areas with fine-grained sediments and high ostracod diversity in more exposed coarser-grained sediments (Morley & Hayward, 2014).

The characteristic of the assemblages occurring in the study samples supported the link between granulometry and calcareous meiofaunal diversity. Coarse-grained sediments were present in the samples TP3-1 (medium sand), TP2-1 (coarse sand), and TP5-1 (very fine gravel), where the ostracod dominance is high, whereas abundance and diversity are low; foraminiferal taxa richness and abundance are low. In the infralittoral zone of the bay, muddy sediments were virtually absent, and Pokorný’s statement was confirmed by higher diversity and lower dominance displayed by ostracod, and, to a lesser extent, by benthic foraminiferal assemblages in fine and very fine sandy samples.

Chemical analyses revealed lower concentrations of pollutants in coarse-grained samples. Arsenic, copper, mercury, lead, zinc, total polycyclic aromatic hydrocarbons (pahs; in particular anthracene, benzo(a)anthracene, benzo(b)fluoranthene, and chrysene) and priority dangerous pahs were significantly higher in fine and very fine sands. The recognized inverse relationship between grain size and anthropogenic chemical contaminants (Horowitz, 1985, 1991; Herut & Sandler, 2006) and the above-mentioned link between meiofaunal remains and granulometry, could erroneously suggest the preference of ostracod and benthic foraminiferal assemblages for polluted bottom sediments. Consequently, to achieve reliable results, we considered separately the fine sandy samples and the coarser-grained sediments.

Calcareous meiofaunal assemblages exhibit different responses to anthropogenic inputs, including the decrease in diversity and increase in the abundance of tolerant taxa or morphotypes typical of polluted waters (Frontalini & Coccioni, 2008, 2011; Yasuhara et al., 2012; Wilkinson et al., 2014). As regards the presence of pollutants in the bottom sediments, researchers assumed that the geochemical nature of the substrate exerts only a modest influence on benthic foraminiferal and ostracod assemblages (Albani et al., 1998; Eagar, 1999; Irizuki et al., 2015); conversely, some investigations suggested negative effects (Schornikov, 2000; Mostafawi, 2001; Martins et al., 2015) or displayed controversial results (Coccioni, 2000; Debenay & Fernandez, 2009; Choi & An, 2012).

Our point of view is that neither the contaminants accumulated in the bottom sediments nor the dead assemblages are representative of the ecological conditions at the time of sampling. Instead, they are indicative of the environmental history of the bay during the last years or decades. On the base of the present data, it is hypothesized that the continued anthropogenic disturbance, primarily due to industrial wastes, had not resulted in a decrease in ostracod and benthic foraminiferal diversity and had its effect by changing the taxonomic composition of the assemblages. In Table 11, the mean values of diversity indices of the infralittoral assemblages of Monte di Procida (Mangoni et al., 2016) and Falerno-Domitio (Balassone et al., 2016) were reported. Both ostracod and benthic foraminiferal assemblages of the Gulf of Pozzuoli showed higher diversity and lower dominance in comparison with the nearby areas. The expected response to contaminant input was the decrease in meiofaunal diversity, nonetheless, some investigations on benthic foraminifers reported inverse trends (Debenay & Fernandez, 2009 and Li et al., 2015, in metal-contaminated waters) and highly diversified assemblages in polluted waters (Romano et al., 2009 and Choi & An, 2012, in metal-contaminated waters; Armynot du Châtelet et al., 2011, in metal-contaminated, rich organic carbon waters). Barras et al. (2014) suggested that in oligotrophic environments the diversity indices are not appropriate to describe water quality. Moreover, the common presence of miliolids implies coastal waters supersaturated in calcium carbonate (Aiello et al., 2018, and references therein) where the presence of industrial wastes did not result in a pH lowering.

Table 11 Comparative table of diversity/dominance indices mean values of benthic foraminiferal and ostracod infralittoral assemblages of northwestern Campanian coastal areas (only fine and very fine sands); MDP = Monte di Procida (Mangoni et al., 2016), FD = Litorale Falerno Domitio (Balassone et al., 2016), GPI = infralittoral assemblages of the Gulf of Pozzuoli (this paper)

Reports of high-diversity ostracod assemblages in stressed environments are very rare (Amato et al., 2019; Aiello et al., 2020) which hampers to perform a proper comparison with the present study.

Conclusions

A quantitative study of benthic foraminiferal and ostracod assemblages, integrated with chemical and sedimentological parameters, was carried out along the coast of the Gulf of Pozzuoli. A total of 11 samples were collected within the infralittoral zone, to investigate the relationship between meiofaunal calcareous remains and pollution indicators preserved in the shallow marine bottom sediments of the bay. Assemblages were characterized by species typical of infralittoral Mediterranean environments such as the benthic foraminifers A. aberdoveyensis, B. granulata, and C. lobatulus, and the ostracods U. margaritifera and P. turbida. The genera Elphidium and Quinqueloculina (Foraminifera), and Semicytherura and Xestoleberis (Ostracoda) were highly diversified.

Despite the continued anthropogenic disturbance, testified by high geoaccumulation levels, the shallow bottom sediments of the Gulf of Pozzuoli yielded high-diversity, low-dominance assemblages. It was here hypothesized that the oligotrophic, well-oxygenated, and supersaturated in CaCO3 shallow Tyrrhenian waters may promote "complex trophic relationship" and "full exploitation of ecological niches" as stated by Holbourn et al. (2013) in their investigations on foraminiferal deep-water assemblages. The diversity was linked to the grain size of the bottom sediments, being higher on fine and very fine sands, and lower on coarser sediments.

On the other hand, the comparison with the assemblages collected in 1961 in the Gulf of Pozzuoli and the Gulf of Naples showed an increase in the abundance of taxa that are suggested to tolerate to geochemical pollution (Q. seminulum, B. elongata, Xestoleberis, Loxoconcha, S. rarecostata). Statistical analysis confirmed the tolerance of Q. lata to a high level of heavy metals contaminants and the preference of U. margaritifera and P. paulii for the shallow part of the infralittoral zone.