Next Article in Journal
Hydrochemical Parameters to Assess the Evolutionary Process of CO2-Rich Spring Water: A Suggestion for Evaluating CO2 Leakage Stages in Silicate Rocks
Next Article in Special Issue
Stable Oxygen and Carbon Isotope Composition of Holocene Mytilidae from the Camarones Coast (Chubut, Argentina): Palaeoceanographic Implications
Previous Article in Journal
Materials and Technologies for the Tertiary Treatment of Produced Water Contaminated by Oil Impurities through Nonfibrous Deep-Bed Media: A Review
Previous Article in Special Issue
Spatial and Temporal Distribution of Chemically Characterized Microplastics within the Protected Area of Pelagos Sanctuary (NW Mediterranean Sea): Focus on Natural and Urban Beaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conservation Paleobiology as a Tool to Define Reference Conditions in Naturally Stressed Transitional Settings: Micropaleontological Insights from the Holocene of the Po Coastal Plain (Italy)

1
Department of Biological, Geological and Environmental Sciences, University of Bologna, Via Zamboni 67, 40126 Bologna, Italy
2
Department of Geological Sciences, Jadavpur University, Kolkata 700032, India
*
Author to whom correspondence should be addressed.
Water 2020, 12(12), 3420; https://doi.org/10.3390/w12123420
Submission received: 31 October 2020 / Revised: 28 November 2020 / Accepted: 1 December 2020 / Published: 4 December 2020
(This article belongs to the Special Issue Climate Change and Anthropogenic Impact on Coastal Environments)

Abstract

:
The key role of paralic environments as providers of ecosystem services, associated with their increasingly threatened state, led to the definition of international water management policies aimed to improve ecological quality status (EcoQs). Restoration actions rely on the definition of reference conditions, which is a particularly challenging task in naturally stressed transitional environments. In the present work, we apply the diversity index Exp(H’bc) on benthic foraminifer assemblages from two anthropogenically unimpacted transitional to coastal Holocene sediment successions of the Po coastal plain, in order to assess past EcoQs (PaleoEcoQs). Ostracod ecological groups provided detailed insights on naturally stressful paleoenvironmental conditions. We show that “poor” to “moderate” PaleoEcoQs are recorded by biological indicators at reference conditions under fluctuations of chemical-physical parameters and organic matter enrichment. We emphasize the importance of a site-specific paleobiological approach, as significant differences in diversity occur even on a short spatial scale. This study illustrates that early to mid-Holocene sediment successions resulted to be appropriate for conservation paleobiological purposes, providing a high-resolution paleoecological record under the influence of the Holocene sea-level rise in analogy with the present-day global change.

1. Introduction

At the border between the continental and the marine realm, transitional settings such as lagoons and estuaries are known for their fragile equilibrium [1]. Enrichment in nutrients, heavy metals, and organic contaminants are often determined by the increasing anthropogenic pressure of the last centuries in the form of aquaculture, agricultural, industrial, and commercial activities. The key role of ecosystem services provided by transitional water habitats led to the definition of national and international conservation guidelines back to ca. 50 years ago (e.g., Federal Water Pollution Control Act, USA, 1972; [2]). More recently, European countries have implemented the European Water Framework Directive (WFD − 2000/60/EC) that aims to achieve “good” or “high” ecological quality status (EcoQs) within 2020 [3]. The definition of reference conditions is essential to perform a correct evaluation of the EcoQs, and implement effective actions for the restoration of potentially threatened coastal and transitional water bodies. Specifically, reference conditions are biological quality elements for a type-specific water body at high ecological status, whereas the ecological status includes “biological, physico-chemical and hydromorphological elements” [4]. Reference conditions should reflect the natural variability of the system, but transitional settings are naturally stressed environments subject to strong fluctuations of many chemical-physical parameters, such as water salinity, turbidity, nutrient, and oxygen concentrations. Biotic communities are dominated by taxa adapted to tolerate such intense, rapid excursions of environmental parameters. Nevertheless, communities with similar features in terms of abundance, structure, and composition also occur in anthropogenically impacted areas. The resulting difficulty of discerning natural stressors from anthropogenic impacts in transitional settings determines the so-called “Estuarine Quality Paradox” [5,6].
An accurate definition of reference conditions and appropriate calibration of boundary thresholds between EcoQs could help in the estimation of anthropogenically induced conditions, preventing transitional environments from appearing more impacted than they really are [7,8]. Ideally, reference conditions should be determined at an almost undisturbed site or, in the second instance, on historical data [4]. Conservation paleobiology can provide an appropriate approach to the definition of reference conditions, as it allows to explore the biotic response to both anthropogenic and natural stressors applying “theories and analytical tools of paleontology to biodiversity conservation” [9,10]. Even though most marine conservation biologists considered a centennial or millennial time-scale appropriate for assessing climate change issues only, such paleontological perspective demonstrated to be essential in the definition of (paleo-)environmental ranges of variation in response to a set of environmental stressors [11], also of natural origin [12].
Foraminifers are part of the benthic invertebrate fauna, indicated as a biological quality element for the definition of EcoQs in coastal and transitional waters [4]. Abundant literature demonstrated the efficiency of this group for (bio-)monitoring purposes (e.g., [13,14]). Guidelines for environmental monitoring through benthic foraminifers were produced for the marine realm [15,16,17], and the developed biotic indices have been currently testing [18,19,20]. Nevertheless, only a few works dealt with the same purpose in transitional waters, and essentially focused on modern environments in Italy [21], located on both the Adriatic and the Tyrrhenian coasts. The reference dataset of Italian transitional waters [21] substantially covers environmental conditions for microtidal systems, from marsh to outer lagoon, and related chemical-physical parameters, such as salinity, pH, freshwater contributions, dissolved oxygen, organic carbon, and trace metals. Environmental changes in these areas occur at a local scale, driven by natural (e.g., seasonality and degree of confinement; [22]) and anthropogenic forcings (e.g., eutrophication and heavy metal contaminations; [23,24]). The high environmental coverage of the modern dataset ensures a reliable comparison with the early to mid-Holocene transitional successions of the Po coastal plain, where similar conditions are commonly recorded [25,26,27]. Ostracods also fall within the huge group of the benthic invertebrate fauna, and share with benthic foraminifers many desirable properties, that make these taxa ideal for integrated paleoecological reconstructions, particularly in transitional settings where strongly reduced salinities may prevent the development of foraminifer communities [28,29]. Even though the response of brackish water ostracods to anthropogenic pollution has not been extensively investigated (e.g., [30,31]), they are excellent proxies of freshwater and organic matter inputs, water turbulence, and permanence in a continuum from continental to back-barrier settings [32,33,34].
In this respect, the Holocene back-barrier succession of the Po coastal plain (N Italy, Figure 1) can represent a valuable resource to investigate pre-impacted conditions in transitional environments. The abundance of threatened coastal water bodies in this area [35,36,37,38] makes such a purpose even more urgent to plan efficient restoration actions. Relationships between past and present-day environmental conditions can be efficiently evaluated through a conservation paleobiology perspective, most commonly achieved through invertebrate macrofauna (e.g., mollusks and corals; [12,39]). In the North Adriatic area, such an approach successfully outlined the long-term resilience and persistence of late Pleistocene-Holocene marine molluscan assemblages to major environmental changes under natural forcings in contrast to the abrupt, unprecedented shift in community composition and structure that occurred in the Anthropocene [40,41]. Seasonal hypoxia, sediment pollution, bottom trawling, and warming events were determined to have a great impact on benthic communities in the last centuries [42,43,44,45], but also impacted the sedimentary record affecting the temporal resolution of fossil assemblages [40,46]. Disentangling sedimentological and biological effects on benthic communities appears as a key to quantifying paleoecological changes in a high-resolution chronostratigraphic framework [47]. In this context, ecosystem fluctuations can be effectively assessed from continuous archives such as back-barrier sediment successions [48,49], where the benthic meiofauna is a proxy of PaleoEcoQs of growing importance, thanks to its abundance, biological properties, and preservation potential [50].
In the present work, we evaluate if reference conditions can be assessed within pre-impacted transitional successions applying ecological indices routinely used to detect anthropogenic pollution in modern environments on benthic meiofauna. We reconstruct the paleoecological quality status (PaleoEcoQs) from benthic foraminifers within two Holocene paralic successions (sensu [52]) of the Po coastal plain. The ostracod record is integrated to obtain a more detailed picture of paleoenvironmental conditions under natural stressors from a separate line of evidence.

2. Geological Setting

The Holocene succession of the Po coastal plain is developed atop of late Pleistocene alluvial sediments and consists of a ca. 30 m-thick transgressive-regressive wedge [53,54]. In response to the Holocene transgression, the development of a wave-dominated estuary is recorded by the superposition of back-barrier to offshore deposits (i.e., retrogradational stacking pattern facies), followed by deltaic sediments typifying the progradation of the Po coastal plain after the sea-level stabilization [55,56].
Thick back-barrier successions are mostly developed in proximity to the maximum shoreline ingression, located up to ca. 26 km west from the modern shoreline in the northern part of the study area, whereas an inland migration of a few km only occurred in the southern portion [51] (Figure 1). Lower transgressive intervals are made up of swamp and lagoon fine-grained sediments, enriched in organic material, testifying the first phases of marine ingression that led to the submersion of wide lowland areas where reducing and sediment-starved conditions locally occurred [56]. Swamp-lagoon alternations often recorded within transgressive successions reflect localized progradational episodes connected to autogenic processes under a predominant allogenic-controlled sedimentation [25]. Lagoon sediments with sand intercalations and coarse-grained transgressive barrier sands at seaward locations are indicative of less restricted conditions developed at peak transgression (7.7 to 7.0 cal. ky BP; [57]). Inland, the shift to highstand conditions is represented by thick swamp to alluvial plain sediments indicative of lagoon infilling, whereas beach ridge and delta front sands at seaward locations express coastal progradation [27].
No anthropogenic influence is recorded in the Holocene transgressive and early regressive sediment succession of the Po coastal plain, where transitional to coastal paleoenvironments are considered as pristine and were used to define geochemical background values at a regional scale [58,59,60]. Specifically, Ref. [36] defined natural background values for the southern Po coastal plain including geochemical data from core 240 S5, reported in [25]. Natural enrichment in Cr and Ni determined by the contribution of mafic and ultramafic rocks cropping out extensively in the Po River drainage basin is restricted to the sediments of the northern Po coastal plain, whereas low values are recorded in the southern coastal plain drained by NE Apenninic rivers [61].

3. Materials and Methods

Cores 240 S5 and 240 S6 were recovered in the southern Po coastal plain, less than 10 km inland from the modern shoreline and 2.6 km distant from each other (Figure 1). Both cores are 40 m long and were drilled by wire-line perforation, which guaranteed a very high recovery percentage (>90%). Sedimentological characteristics of the cores, including accessory components such as macrofossils and vegetable debris are described in [62]. Further details on core 240 S5, including geochemical and qualitative micropaleontological data are reported in [25]. In the present work, we present the micropaleontological record of core 240 S6 and provide more detailed data from core 240 S5 using quantitative benthic foraminifer and semi-quantitative ostracod fauna.
The lower portion of both cores is made up of continental deposits of late Pleistocene age, passing upward to a ca. 15 m-thick back-barrier succession including a suite of swamp to lagoon facies that evolve into shallow marine and coastal sediments in core 240 S6 only [62]. Subtle changes in the composition of microfossil assemblages within lagoonal intervals are interpreted to reflect short-term paleoenvironmental changes connected to variations in sediment supply at a local scale under a general trend of sea-level rise [25]. Both successions are capped by modern alluvial plain sediments [62].
A total of 129 samples were collected for micropaleontological analysis of core 240 S6, with a thickness of about 1.5 to 2 cm, corresponding to ca. 100 g of sediment. Benthic foraminifer and ostracod fauna was recovered following the standard procedure adopted for sample treatment in core 240 S5 [25] and other reference works dealing with the microfossil content of the Po Plain Holocene successions (e.g., [63]). It involves: (i) oven-drying the samples at 60 °C for at least eight hours; (ii) soaking in water or water + H2O2 (35% vol.) for 24 hours; (iii) wet-sieving at 63 µm; (iv) oven-drying again. Sixty-three samples containing well-preserved and abundant meiofauna were selected for micropaleontological investigations, performed in the sediment fraction >125 µm. The set of 109 samples of core 240 S5 of [25] was reviewed selecting suitable samples for quantitative analyses.
Identification of benthic foraminifers and ostracods followed the original taxonomic descriptions [64,65] and reference works mostly focused on the Mediterranean area as reported in Text S1. When possible, at least 300 well-preserved benthic foraminifer specimens were counted and identified within 41 samples (25 from core 240 S5, 16 from core 240 S6), dividing the residue into small portions with a microsplitter if necessary (abundance data provided in Table S1), otherwise, all the retrieved tests were counted. Relative abundances of ostracods were expressed by four categories (abundant: >30%, common 10% to 30%, scarce 4% to 10%, rare <4%) in a total of 95 samples (54 from core 240 S5, 41 from core 240 S6; Table S2). We considered adult and juvenile ostracods bearing morphological characters that allowed an identification at species level, whereas valves with evident signs of transport or reworking (dark appearance, strong abrasions, fragmented, etc. [66]) were discarded.
To test if natural reference conditions could be assessed by benthic foraminifers, we applied the diversity index proposed by [13] as a proxy for (Paleo-)EcoQs in cores 240 S5 and 240 S6, which represent transitional to marine Holocene paleoenvironments with pristine conditions. The index follows the [67] modification (H’bc) of the Shannon-Weaver index of community entropy (H’; [68]), which fixes the problem of possibly unobserved species in the case of under-sampling. Then, the exponential of the obtained number (exp(H’bc)) allows the shift to a true diversity measure that expresses the number of species with the same abundance required to obtain H’bc [13]. Considering the paleoenvironmental setting, we applied the classes of EcoQs defined for present-day Italian transitional waters [21]. Guidelines for the determination of EcoQs were established by [13], which recommend the use of living assemblages reflecting real-time environmental conditions, nonetheless dead assemblages that represent a longer time frame show a lower, but still significant, degree of correlation with oxygen content. Even though many studies dealing with (Paleo-)EcoQs are performed on the >63 µm size fraction, Ref. [13] suggest that the >125 µm fraction is also adequate due to the negligible information lost, the lower time and taxonomic effort required for the analysis. Fossilizable taxa obtained through sample drying are considered to provide an accurate estimation of (paleo-)environmental conditions to compute (Paleo-)EcoQs, as fragile agglutinated foraminifers preserved by wet-picking do not substantially improve the (paleo-)ecological signal [13].
The Foram−AMBI index was calculated to estimate the stress related to organic carbon concentration on the basis of benthic foraminifer ecological groups with different responses to organic enrichment [15]. Despite the Foram−AMBI being developed for the marine realm, we apply it to the analyzed assemblages to assess the abundance of organic matter sensitive, tolerant, and opportunistic taxa in transitional settings under natural conditions. We used benthic foraminifer groups defined for the Mediterranean area by [17], specifically: GI—sensitive species; GII—indifferent species; GIII—third-order opportunists; GIV—second-order opportunists; GV—first-order opportunists. We computed the exp(H’bc) and Foram−AMBI indices in benthic foraminifer assemblages including >100 tests, which is considered a suitable threshold value for (paleo-)ecological analyses on assemblages that commonly show low abundances in stressed settings [69,70].
Finally, a detrended correspondence analysis (DCA) was computed on benthic foraminifers showing relative frequencies >1% with the aim of exploring the relationship between taxa, assemblages, and the computed ecological indices through paleoenvironmental gradients. DCA is an indirect ordination technique that allows for the assigning of scores to both species and samples in the multivariate space and was successfully applied in back-barrier to shallow marine settings for (paleo-)ecological purposes [71,72,73]. Here, we apply DCA separately for the benthic foraminifer datasets of core 240 S5 and 240 S6 to visualize the correlation between environmental gradients and ecological indices, plotted as additional variables within a more restricted setting (240 S5) and throughout the transition from back-barrier to shallow marine settings (240 S6). A few taxa characteristic of inner to outer shelf environments (Cibicidoides lobatulus, Asterigerinata planorbis, Praeglobobulimina sp.) possibly transported in the restricted back-barrier setting of core 240 S5 that determined strong distortion of the detected paleoecological gradient were removed. The DCA was computed in the R environment (version 3.6.1, R Foundation for Statistical Computing, Vienna, Austria; [74]) using the vegan package [75].
Ostracods provided detailed insights on paleoecological conditions, particularly in back-barrier sediments, integrating the record of benthic foraminifers. Autoecological data on individual taxa were mostly extracted from [76] and additional references when necessary [77,78,79], which allowed to obtain insights on water salinity, oxygen, turbulence, and organic matter inputs (Table S3) also within sediment intervals with no benthic foraminifers. Tolerance to pH oscillations and water permanence was also inferred for assemblages with freshwater to low brackish taxa [76], whereas no data are available for brackish, euryhaline, and marine species.
The chronological framework of the two cores is provided by seven radiocarbon dates performed on samples enriched in organic matter and mollusk shells [25,54]—(Table 1). Calibration of conventional radiocarbon ages was performed with OxCal 4.2 [80], applying the Intcal13 calibration curve or the Marine13, with a ΔR value of 139 ± 28 years [81], and a mixed contribution in back-barrier samples.

4. Results

4.1. Benthic Foraminifers

A total of 15138 benthic foraminifer specimens, belonging to 38 taxa, were recorded in the analyzed successions; in particular, core 240 S5 included 7949 tests represented by 19 taxa from the 17.95 to 11.85 m core depth, and 7189 specimens of 24 taxa were recorded from core 240 S6 in the interval between 23.60 and 10.10 m core depth (Table S1).

4.1.1. Core 240 S5

Exp(H’bc) index showed values <3 in almost all samples, with the exception of the lowermost and uppermost samples (17.95 and 11.85 m core depth) characterized by values slightly >3. The resulting PaleoEcoQs is classified as “poor” from 17.55 to 11.95 m core depth and as “bad” in the remaining samples (Figure 2, Table S1).
The taxonomic composition of fossil assemblages determined moderately high Foram−AMBI values, mostly between 3.5 and 5, excluding the uppermost sample at 11.85 m core depth that shows a value of 3.2. The most abundant ecological group was represented by second-order opportunist taxa, with an average abundance of 88.3% and reaching peak frequencies >95%. Third-order opportunists (tolerant) taxa were relatively scarce (7.6% on average), whereas indifferent and sensitive species were rare (0.6% and 2.3%, respectively) and mostly recorded in the samples between 12.6 and 11.85 m core depth (Figure 2).
The low diversity of samples from the lagoon interval of core 240 S5 is mostly determined by the strong dominance of Ammonia tepida (88.3% on average), associated with lower abundances of Haynesina germanica (7.5%). The former, classified as a second-order opportunist, determined the high Foram−AMBI values throughout the succession in association with tolerant species which were mostly represented by H. germanica. Other taxa, such as the sensitive Ammonia parkinsoniana and the indifferent Haynesina depressula, mostly occur with relatively low abundances in the upper portion of the lagoon sediments (Figure 2).
DCA arranged the benthic foraminifer variables in a wedge-shaped distribution (Figure 3), with DCA axis one accounting for 58% of the total variance. Within this relatively short paleoenvironmental gradient, taxa are distributed along axis one in accordance with their degree of tolerance to organic enrichment. Samples with a higher abundance of sensitive species are placed on the right side of the plot. Foram−AMBI and exp(H’bc) vectors showed a good correlation with axis one.

4.1.2. Core 240 S6

The first abundant benthic foraminifer assemblages are recorded in two single samples at 23.60 and 21.45 m core depth within the lagoon-swamp deposits, which showed exp(H’bc) values of 3.1 and 3.2, respectively. Above an interval barren in benthic foraminifers attributed to a swamp paleoenvironment from 21.20 to 18.80 m core depth, lagoon sediments were mostly characterized by values of the diversity index exp(H’bc) < 3, with a slightly higher value in one sample at 16.6 m core depth. An increase of the exp(H’bc), with values between 5.8 and 8.4, is recorded within shoreface deposits above 13.10 m core depth. PaleoEcoQs in the lowermost lagoon samples is determined as “bad”, whereas it resulted to be “poor” in the lagoon sediments between 18.80–13.35 m core depth, with the exception of one sample classified as “bad” at 16.60 m core depth. PaleoEcoQs within the uppermost coastal sediments resulted to be oscillating between “bad” and “moderate” (Figure 4, Table S1).
Foram−AMBI was relatively high in the lower assemblages, which showed values of 3.4 and 3.2 at 23.60 and 21.45 m core depth, respectively, and ranged from 3.3 to 3.8 up to 13.70 m core depth. Foram−AMBI dropped at 1.1 at 13.00 m core depth, and relatively low values between 1.3 and 1.5 persisted in the uppermost assemblages (Figure 4).
The low diversity and the moderately high values of Foram−AMBI in the lagoon assemblages are mostly determined by the dominance of A. tepida (with relative abundances even >70%), classified as second-order opportunists, which represented the most abundant taxon throughout the analyzed succession (45.7% on average). In addition, moderately high frequencies of tolerant taxa, such as H. germanica (with peak frequencies of ca. 45%) and Aubignyna perlucida (up to 5.7%; Table S1) contributed to an increase in the Foram−AMBI. Upcore, the higher diversity and the drop of Foram−AMBI within coastal sediments is mostly produced by the high abundance of indifferent species, as Cribroelphidium lidoense and Ammonia beccarii (up to 29% and 51%, respectively). These taxa are associated with a variety of sensitive species with moderate abundances, mostly represented by Ammonia parkinsoniana (up to >35%) and lower frequencies of Adelosina cliarensis, Buccella granulata, Elphidium crispum, Elphidium macellum, and Siphonaperta aspera (relative frequencies of individual taxa <5%; Table S1).
Taxa and samples in the DCA are mostly distributed along DC1 that explained 76% of the total data variance (Figure 5). Along axis one, most marine taxa sensitive or indifferent to organic enrichment are placed on the right of the plot. Whereas, the left side gathered the second-order opportunist A. tepida and other species, such as H. germanica and Aubignyna perlucida, tolerant to fluctuations in salinity and organic matter and thriving in restricted, fluvial-influenced conditions [82,83,84]. This is confirmed by the depositional environment attributed to the samples, as negative scores characterized lagoon sediments, whereas shoreface samples showed positive scores. The orientation of environmental vectors suggests a highly positive and negative correlation of exp(H’bc) and Foram−AMBI, respectively, with ordination axis one.

4.2. Ostracods

Twenty-five ostracod taxa were identified from the two study successions. Specifically, 11 taxa were found between 18.95 and 4.30 m core depth in core 240 S5, whereas core 240 S6 included 21 taxa from 24.50 to 11.60 m core depth (Table S2).

4.2.1. Core 240 S5

Ostracod taxa tolerant to low oxygen concentrations are abundant throughout the core, mostly represented by Pseudocandona albicans, Candona neglecta, and Ilyocypris spp. within swamp sediments, whereas Cyprideis torosa is dominant within lagoon sediments (Figure 6). Scattered valves of Cytherois fischeri recorded from 16.00 to 15.50 m core depth determined the rare occurrence of taxa characteristic of highly oxygenated conditions. According to salinity, freshwater to low brackish taxa are most abundant within swamp sediments, whereas the euryhaline C. torosa is dominant within lagoon deposits, that also included variable amounts of taxa tolerant to hypohaline conditions (mostly P. albicans and Ilyocypris spp.) in the lower and upper portions. Rare brackish-marine and marine species are recorded in scattered lagoon samples between 16.00 and 14.20 m core depth, represented by C. fischeri and Palmoconcha turbida. Concerning water turbulence, most taxa within swamp sediments are polyrheophilic and mesorehoephilic, and rheoeuryplastic within lagoon deposits, respectively, driven by P. albicans and C. torosa. Within swamp sediments, most taxa are titanoeuryplastic and able to live in both permanent and temporary waters.

4.2.2. Core 240 S6

The high abundance of ostracods tolerant to reduced oxygen concentrations throughout the succession of core 240 S6 is mostly determined by the high abundance of P. albicans and C. torosa in swamp and lagoon sediments, respectively. The presence of taxa associated with high oxygen concentrations is mostly recorded in lagoon deposits from 17.70 to 13.75 m core depth due to the common to scarce occurrence of Leptocythere ex gr. castanea. Freshwater to low brackish species are mostly recorded in the lower part of the analyzed succession, from 24.50 to 18.80 m core depth, where the fauna is predominantly composed of P. albicans associated with lower occurrences of Candona candida, and Ilyocypris spp. This interval also included rare brackish and brackish-marine taxa, and very rare occurrences of the marine Pontocythere turbida, possibly transported from nearby marine settings. The euryhaline C. torosa is abundantly found within scattered intervals from 24.30 to 21.20 m core depth. Lagoon sediments from 18.80 to 13.10 m core depth include associations of euryhaline, brackish, and brackish-marine taxa (e.g., C. torosa, Loxoconcha elliptica, Leptocythere bacescoi), whereas exclusively marine species occur in shoreface sands. Freshwater to brackish ostracod communities are mostly composed of polyrheophilic to mesorheophilic and rheoeuryplastic taxa, whereas rheophobic to oligorheophilic species (Leptocythere ex gr. castanea, L. elliptica) occurred above 17.70 m core depth. Up to 18.80 m core depth, assemblages are composed of titanoeuryplastic taxa, characteristic of both permanent and temporary waters (Figure 7).

5. Discussion

5.1. PaleoEcoQs through Short Paleoenvironmental Gradients: Restricted Back-Barrier Settings (Core 240 S5)

A relatively short paleoenvironmental gradient is depicted by the DCA samples of core 240 S5, deposited in a back-barrier, restricted lagoon setting (Figure 2 and Figure 3). The composition of benthic foraminifer assemblages is rather stable with dominant A. tepida and subordinate H. germanica, taxa characteristic of partially barred environments enriched in organic matter and with reduced salinity [85,86]. Specifically, the former is able to rapidly consume high rates of both fresh and refractory organic matter, whereas the latter tolerates moderate organic enrichments, showing lower consumption rates of organic material [87]. The upward increasing concentrations of taxa commonly reported from lagoons but with a relatively low tolerance to organic enrichment, such as H. depressula and A. parkinsoniana [86], suggests lower organic matter levels. Despite the relatively short paleoenvironmental gradient, the Foram−AMBI index confirms the organic matter trend derived from benthic foraminifers underlain by DCA axis one (Figure 3), even under restricted conditions.
Dominant ostracod taxa (e.g., P. albicans, C. torosa) confirm abundant organic matter and reduced oxygenation, in addition to oscillating, reduced salinities and possible water turbulence during the development of lagoon conditions (Figure 6). Indeed, the dominant C. torosa is one of the most opportunistic ostracod species, able to thrive on extremely variable conditions, especially under high organic matter concentrations [88,89]. The scattered freshwater to low brackish ostracods within lagoon deposits, such as P. albicans and C. fragilis, also tolerate fluctuating pH and water table [76] (Figure 6) and likely suggest the development of unstable paleoenvironmental conditions. High amounts of organic materials are also suggested by geochemical indicators [25] based on the ratio between loss of ignition (LOI) and CaO of carbonates (Figure 8). LOI quantifies the total volatile content in a sediment sample, including humidity, H2O in the lattice of clay minerals, CO2 in the carbonate minerals, and organic matter. Even though the individual contribution of each element cannot be discerned, the integration with selected major elements still provides a relevant assessment of partitioning among them. The majority of lagoon and swamp samples are located below the LOI-carbonate mixing line towards the LOI apex, suggesting organic enrichment and low carbonate concentrations, and possibly lowered amounts of oxygen and oscillating pH in relationship with organic matter consumption [90].
Unstable paleoenvironmental conditions, confirmed by ostracod assemblages, determined the low diversity of the benthic foraminifer fauna within lagoon sediments of core 240 S5 (Figure 2). Despite the low geochemical background levels indicative of unimpacted, natural conditions, very low values of exp(H’bc) resulted in the definition of “poor” or “bad” PaleoEcoQs, on the basis of threshold values defined by [21] for Italian transitional waters. This is in accordance with the Estuarine Quality Paradox, which prevents discrimination between natural and anthropogenic stresses [5,6] relying on biotic indicators, including the macrofauna [2,91].

5.2. PaleoEcoQs through Long Paleoenvironmental Gradients: From Transitional to Coastal Settings (Core 240 S6)

A well-developed paleoenvironmental gradient is evidenced by the DCA in core 240 S6, spanning from the lagoon to shoreface conditions (Figure 5). Species are distributed along axis one in accordance with their response to organic enrichment, which in turn is linked to the distribution of taxa in the back-barrier setting. The opportunists A. tepida, A. perlucida, and H. germanica are abundant in Mediterranean transitional waters, from inner lagoon to shallow prodelta settings in proximity to river outlets [22,82,86,92]. Cribroelphidium oceanense is also reported from inner lagoon environments enriched in organic matter, and H. depressula seems to tolerate variable inputs of organic material in back-barrier settings, where it is commonly recorded with moderate frequencies [86,93]. The remaining taxa, mostly restricted to shoreface sediments, show a trend of increasing sensitivity to organic matter towards positive values of DCA axis one, with A. parkinsoniana located in an intermediate position due to its common occurrence in both outer lagoon and shallow marine settings [86,94,95]. In addition, vector orientation of the Foram−AMBI index (Figure 5) confirms a strong influence of organic matter on the benthic foraminifer distribution from back-barrier to coastal settings, which is an important factor controlling the degree of confinement in paralic environments, along with salinity, water circulation, turbidity, etc. [83].
Ostracod ecological groups confirm stressful conditions with relatively low oxygen concentration, fluctuating salinity, and water turbulence throughout back-barrier sediments (Figure 7). Above mid-Holocene lagoon-swamp alternations recording strong paleoenvironmental fluctuations (in terms of oxygen, salinity, pH, and water table; Figure 7), thick lagoon deposits with low to moderate proportions of Leptocythere ex gr. castanea possibly suggests a tendency towards more stable, higher oxygenated conditions [78]. Nevertheless, the abundance of brackish taxa tolerant to organic enrichments such as C. torosa, L. elliptica, L. stellifera, and L. bacescoi [76,96,97] are indicative of high organic matter concentrations at the bottom and, concurrently, of increasing marine influence.
The taxonomic composition of benthic assemblages suggests a lower enrichment in organic matter and, in turn, higher oxygen conditions in core 240 S6, where exp(H’bc) values are higher than in core 240 S5 (Figure 2 and Figure 4). Similarly, the general upward increase of the diversity benthic foraminifer index exp(H’bc) reflects a decrease in organic matter and the development of less restricted conditions in response to the Holocene transgression. Both back-barrier and coastal deposits are affected by an estimated PaleoEcoQs lower than it really is, which suggests a re-evaluation of thresholds values also in low confined (paleo-)environments.

5.3. Temporal Resolution of the Fossil Record

Living assemblages, used to define present-day EcoQs, inevitably differ from dead assemblages in taxonomic composition and structure as multiple generations are accumulated over time and space within the fossil record: this phenomenon, defined as time-averaging, has a key importance in paleoecology and, in turn, in conservation paleobiology [98,99]. In addition, taphonomic action occurring within the surface mixed-sediment layers further dampen the biological signal recorded by benthic communities [100], increasing the temporal scale of single assemblages to the decadal, centennial, or millennial time-scales [101]. Time-averaging tends to increase the diversity of assemblages, with a higher proportion of rare stenotopic taxa at the expense of abundant species [98]. Bioturbation is reported within back-barrier sediments of the Po coastal plain [51] but the presented study successions still exhibit extremely low diversity values (Figure 2 and Figure 4), despite the potential time-averaging effects that could have mixed communities with different composition. Paralic environments were affected by rapid sediment accumulation rates, due to the high accommodation space available after the Holocene transgression (<7.0 ky BP) [57] and the high subsidence rates reached in specific depocenters, such as our study area, connected to the structural setting of the Po basin [51]. These characteristics make the transitional succession of the Po coastal plain appropriate for paleobiological investigations, especially through benthic foraminifers which are affected by a lower time-averaging than the bigger mollusk shells [102].

5.4. How to Define Reference Conditions?

Directions of the European Commission [4] affirmed “an urgent need” of data that would allow for the defining of reference conditions in naturally variable settings, such as transitional and coastal environments. Among benthic taxa, foraminifers resulted to be suitable biological quality elements in modern [21] and fossil [69] settings. A conservation paleobiology perspective extended through the Holocene age would cover many aspects of the natural variability of coastal systems in time and space and would also allow to investigate extended (paleo-)environmental gradients that are difficultly encountered in short-time intervals [103], ensuring an accurate assessment of PaleoEcoQs.
Geochemical data from our study material are considered as reference background values for the southern Po coastal plain [36] representing natural conditions at “high” PaleoEcoQs, nevertheless, the low-diversity associated fauna indicates low PaleoEcoQs. In the presence of anthropogenically unimpacted environmental conditions, a possible solution to overcome the Estuarine Quality Paradox issue is represented by the definition of threshold values for (Paleo-)EcoQs at a local scale, as also evidenced in pre-industrial to present-day sediment successions of the Norwegian coasts [50,69]. The principle of a location-specific approach is widely adopted in paleoecology, as dissimilarities between modern and fossil assemblages highly impact on final results (e.g., [104]). Such a location-specific approach traditionally implies an adequate sampling of modern data to develop quantitative paleoenvironmental reconstructions. In the field of conservation paleobiology and, in turn, for biomonitoring purposes, the problem is reversed: a high-resolution chronostratigraphic sampling at a local scale would allow to cover a wide range of past analogs (i.e., reference conditions) to properly assess present-day EcoQs. In our case, exp(H’bc) thresholds defined for Italian transitional waters are too high to correctly assess reference conditions in the analyzed transitional successions, despite some modern sites being located only ca. 100 km away from our study cores (e.g., Venice Lagoon; [21]). In addition, the high lateral and cross variability of transitional settings determined significant differences in exp(H’bc) values, low for the lagoon assemblages of core 240 S6, and even lower in the more restricted, organic enriched setting of core 240 S5 (Figure 2 and Figure 4) in spite of the few km among them. Taking into account such aspects, biotic indices should be lowered in accordance with faunal diversity at a very local scale.
Assemblages at high ecological status should, theoretically, decrease their diversity in the case of anthropogenic impact. Nevertheless, the assessment of EcoQs could be challenging when low diversity assemblages are already recorded under pristine conditions, as for core 240 S5. In such a case, would the first stages of anthropogenic impacts be detected? Low diversity assemblages impacted by additional stressors either locally disappear when environmental conditions overcome organisms’ adaptive means or do not show changes in community structure due to the eurytopic nature of the dominant opportunistic taxa [69,105]. In such cases, diversity indices do not work properly and inferences about ecological status can be performed through geochemical indicators [69]. Indeed, further studies linking the early to mid-Holocene unimpacted assemblages with benthic communities subject to the increasing anthropogenic pressures of the last centuries are required for the Mediterranean area, in order to test benthic foraminifers as biological quality elements even in transitional environments.

6. Conclusions

The anthropogenically undisturbed Holocene succession of the Po coastal plain allows for the assessment of ecological conditions in transitional settings by means of benthic foraminifers, which represent biological quality elements of growing importance [13,15]. Additional paleoenvironmental insights were efficiently provided by benthic ostracods ecological groups, demonstrating excellent integration with the benthic foraminifer record.
At the locations of core 240 S5 and 240 S6, we showed that reference conditions estimated from the benthic foraminifer exp(H’bc) are much lower than the “high” PaleoEcoQs expected in the pristine conditions of the Holocene sediments. In particular, the diversity decreased from coastal towards restricted back-barrier settings in accordance with the Estuarine Quality Paradox.
Within Holocene deposits naturally enriched in organic matter, where also oxygen deficiency potentially occurred, threshold values of (Paleo-)EcoQs should be established at a local scale. Values determined for Italian transitional waters underestimated the PaleoEcoQs at our study sites, which also showed significant differences among them despite the short geographical distance, due to the influence of sea-level rise.
We demonstrated that the conservation paleobiology approach offers a great opportunity to explore high ranges of natural variability throughout both short and long (paleo-)environmental gradients, in accordance with the suggestions of the European WFD [4]. The throwback to early or mid-Holocene provides “long-term” changes not necessarily recorded within recent sediment successions traditionally adopted for paleoecological purposes (i.e., back to the pre-industrial era), but possibly developed in present-day environments subject to the increasing pressure of global change.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/12/12/3420/s1, Text S1: List of references adopted for taxonomic attribution of benthic foraminifers and ostracods. Table S1: Benthic foraminifer abundance, Shannon H, Chao, and Shen (2003) H’bc, exp(H’bc) values with related PaleoEcoQs, Foram−AMBI values with related classification and percentages of individual ecological groups in cores 240 S5 and 240 S6. Table S2: Ostracod semi-quantitative abundance in cores 240 S5 and 240 S6. Table S3: Autoecological characteristics of ostracod taxa from cores 240 S5 and 240 S6, based on [76,78,79,106,107,108,109] and references therein.

Author Contributions

Conceptualization, G.B. and S.C.V.; sampling, sample treatment and methodology, G.B., S.C.V., V.R., and A.G.; multivariate elaborations, G.B.; interpretation, G.B., S.C.V., and V.R.; writing—Original draft preparation, G.B.; writing—Review and editing, S.C.V., V.R., and A.G.; figure preparation, G.B.; supervision, S.C.V. and V.R.; project administration, S.C.V.; funding acquisition, S.C.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Regione Emilia-Romagna—RER—Geological, Seismic, and Soil Survey is acknowledged for providing core material. Anupam Ghosh acknowledges Erasmus Mundus Lot 13 Academic Exchange Programme for this collaborative research work. Authors thank the two anonymous reviewers for their comments on the early version of the manuscript, with a special acknowledgment to the first reviewer for the constructive suggestions that greatly enhanced the quality of the work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lloret, J.; Marín, A.; Marín-Guirao, L. Is coastal lagoon eutrophication likely to be aggravated by global climate change? Estuar. Coast. Shelf Sci. 2008, 78, 403–412. [Google Scholar] [CrossRef]
  2. Munari, C.; Mistri, M. The performance of benthic indicators of ecological change in Adriatic coastal lagoons: Throwing the baby with the water? Mar. Pollut. Bull. 2008, 56, 95–105. [Google Scholar] [CrossRef] [PubMed]
  3. European Parliament. Directive 2008/56/EC of the European Parliament and of the Council of 17 June 2008 establishing a framework for community action in the field of marine environmental policy (Marine Strategy Framework Directive). Off. J. Eur. Union 2008, L164, 1–51. [Google Scholar]
  4. European Commission. Common Implementation Strategy for the Water Framework Directive (2000/60/EC). Transitional and Coastal Waters: Typology, Reference Conditions and Classification Systems; Guidance Document. No 5; Office for Official Publications of the European Communities: Luxembourg, 2003. [Google Scholar]
  5. Dauvin, J.-C. Paradox of estuarine quality: Benthic indicators and indices, consensus or debate for the future. Mar. Pollut. Bull. 2007, 55, 271–281. [Google Scholar] [CrossRef]
  6. Elliott, M.; Quintino, V. The Estuarine Quality Paradox, Environmental Homeostasis and the difficulty of detecting anthropogenic stress in naturally stressed areas. Mar. Pollut. Bull. 2007, 54, 640–645. [Google Scholar] [CrossRef]
  7. Ruellet, T.; Dauvin, J.-C. Benthic indicators: Analysis of the threshold values of ecological quality classifications for transitional waters. Mar. Pollut. Bull. 2007, 54, 1707–1714. [Google Scholar] [CrossRef]
  8. Fitch, J.E.; Cooper, K.M.; Crowe, T.P.; Hall-Spencer, J.M.; Phillips, G. Response of multi-metric indices to anthropogenic pressures in distinct marine habitats: The need for recalibration to allow wider applicability. Mar. Pollut. Bull. 2014, 87, 220–229. [Google Scholar] [CrossRef]
  9. Dietl, G.P.; Flessa, K.W. Conservation paleobiology: Putting the dead to work. Trends Ecol. Evol. 2011, 26, 30–37. [Google Scholar] [CrossRef]
  10. Dietl, G.P.; Kidwell, S.M.; Brenner, M.; Burney, D.A.; Flessa, K.W.; Jackson, S.T.; Koch, P.L. Conservation paleobiology: Leveraging knowledge of the past to inform conservation and restoration. Annu. Rev. Earth Planet. Sci 2015, 43, 79–103. [Google Scholar] [CrossRef]
  11. Dietl, G.P.; Smith, J.A.; Durham, S.R. Discounting the past: The undervaluing of paleontological data in conservation science. Front. Ecol. Evol. 2019, 7, 108. [Google Scholar] [CrossRef] [Green Version]
  12. O’Dea, A.; Lepore, M.; Altieri, A.H.; Chan, M.; Morales-Saldaña, J.M.; Muñoz, N.-H.; Pandolfi, J.M.; Toscano, M.A.; Zhao, J.; Dillon, E.M. Defining variation in pre-human ecosystems can guide conservation: An example from a Caribbean coral reef. Sci. Rep. 2020, 10, 2922. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Bouchet, V.M.P.; Alve, E.; Rygg, B.; Telford, R.J. Benthic foraminifera provide a promising tool for ecological quality assessment of marine waters. Ecol. Indic. 2012, 23, 66–75. [Google Scholar] [CrossRef]
  14. Dimiza, M.D.; Triantaphyllou, M.V.; Koukousioura, O.; Hallock, P.; Simboura, N.; Karageorgis, A.P.; Papathanasiou, E. The Foram Stress Index: A new tool for environmental assessment of soft-bottom environments using benthic foraminifera. A case study from the Saronikos Gulf, Greece, Eastern Mediterranean. Ecol. Indic. 2016, 60, 611–621. [Google Scholar] [CrossRef]
  15. Alve, E.; Korsun, S.; Schönfeld, J.; Dijkstra, N.; Golikova, E.; Hess, S.; Husum, K.; Panieri, G. Foram−AMBI: A sensitivity index based on benthic foraminiferal faunas from North-East Atlantic and Arctic fjords, continental shelves and slopes. Mar. Micropaleontol. 2016, 122, 1–12. [Google Scholar] [CrossRef]
  16. Alve, E.; Hess, S.; Bouchet, V.M.P.; Dolven, J.K.; Rygg, B. Intercalibration of benthic foraminiferal and macrofaunal biotic indices: An example from the Norwegian Skagerrak coast (NE North Sea). Ecol. Indic. 2019, 96, 107–115. [Google Scholar] [CrossRef]
  17. Jorissen, F.; Nardelli, M.P.; Almogi-Labin, A.; Barras, C.; Bergamin, L.; Bicchi, E.; El Kateb, A.; Ferraro, L.; McGann, M.; Morigi, C.; et al. Developing Foram−AMBI for biomonitoring in the Mediterranean: Species assignments to ecological categories. Mar. Micropaleontol. 2018, 140, 33–45. [Google Scholar] [CrossRef] [Green Version]
  18. Melis, R.; Celio, M.; Bouchet, V.M.P.; Varagona, G.; Bazzaro, M.; Crosera, M.; Pugliese, N. Seasonal response of benthic foraminifera to anthropogenic pressure in two stations of the Gulf of Trieste (northern Adriatic Sea, Italy): The marine protected area of Miramare versus the Servola water sewage outfall. Medit. Mar. Sci. 2019, 20, 120–141. [Google Scholar] [CrossRef]
  19. El Kateb, A.; Stalder, C.; Martínez-Colón, M.; Mateu-Vicens, G.; Francescangeli, F.; Coletti, G.; Stainbank, S.; Spezzaferri, S. Foraminiferal-based biotic indices to assess the ecological quality status of the Gulf of Gabes (Tunisia): Present limitations and future perspectives. Ecol. Indic. 2020, 111, 105962. [Google Scholar] [CrossRef]
  20. Klootwijk, A.T.; Alve, E.; Hess, S.; Renaud, P.E.; Sørlie, C.; Dolven, J.K. Monitoring environmental impacts of fish farms: Comparing reference conditions of sediment geochemistry and benthic foraminifera with the present. Ecol. Indic. 2021, 120, 106818. [Google Scholar] [CrossRef]
  21. Bouchet, V.M.P.; Goberville, E.; Frontalini, F. Benthic foraminifera to assess Ecological Quality Statuses in Italian transitional waters. Ecol. Indic. 2018, 84, 130–139. [Google Scholar] [CrossRef]
  22. Frontalini, F.; Margaritelli, G.; Francescangeli, F.; Rettori, R.; Armynot du Châtelet, E.; Coccioni, R. Benthic foraminiferal assemblages and biotopes in a coastal lake: The case study of Lake Varano (Southern Italy). Acta Protozool. 2013, 52, 147–160. [Google Scholar]
  23. Coccioni, R.; Frontalini, F.; Marsili, A.; Mana, D. Benthic foraminifera and trace element distribution: A case-study from the heavily polluted lagoon of Venice (Italy). Mar. Pollut. Bull. 2009, 59, 257–267. [Google Scholar] [CrossRef] [PubMed]
  24. Frontalini, F.; Coccioni, R.; Bucci, C. Benthic foraminiferal assemblages and trace element contents from the lagoons of Orbetello and Lesina. Environ. Monit. Assess. 2010, 170, 245–260. [Google Scholar] [CrossRef] [PubMed]
  25. Dinelli, E.; Ghosh, A.; Rossi, V.; Vaiani, S.C. Multiproxy reconstruction of late Pleistocene-Holocene environmental changes in coastal successions: Microfossil and geochemical evidences from the Po Plain (Northern Italy). Stratigraphy 2012, 9, 153–167. [Google Scholar]
  26. Vaiani, S.C.; Pennisi, M. Tracing freshwater provenance in palaeo-lagoons by boron isotopes and relationship with benthic foraminiferal assemblages. A comparison from late Quaternary subsurface successions in Northern and Central Italy. Boll. Soc. Paleontol. Ital. 2014, 53, 55–66. [Google Scholar] [CrossRef]
  27. Amorosi, A.; Bruno, L.; Campo, B.; Morelli, A.; Rossi, V.; Scarponi, D.; Hong, W.; Bohacs, K.M.; Drexler, T.M. Global sea-level control on local parasequence architecture from the Holocene record of the Po Plain, Italy. Mar. Pet. Geol. 2017, 87, 99–111. [Google Scholar] [CrossRef]
  28. Pint, A.; Engel, M.; Melzer, S.; Frenzel, P.; Plessen, B.; Brückner, H. How to discriminate athalassic and marginal marine microfaunas: Foraminifera and other fossils from an early Holocene continental lake in northern Saudi Arabia. J. Foraminifer. Res. 2017, 47, 175–187. [Google Scholar] [CrossRef]
  29. Barbieri, G.; Vaiani, S.C. Benthic foraminifera or Ostracoda? Comparing the accuracy of palaeoenvironmental indicators from a Pleistocene lagoon of the Romagna coastal plain (Italy). J. Micropalaeontol. 2018, 37, 203–230. [Google Scholar] [CrossRef]
  30. Pascual, A.; Rodriguez-Lazaro, J.; Weber, O.; Jouanneau, J.M. Late Holocene pollution in the Gernika estuary (southern Bay of Biscay) evidenced by the study of Foraminifera and Ostracoda. Hydrobiologia 2002, 475, 477–491. [Google Scholar]
  31. Ruiz, F.; Abad, M.; Bodergat, A.M.; Carbonel, P.; Rodríguez-Lázaro, J.; Yasuhara, M. Marine and brackish-water ostracods as sentinels of anthropogenic impacts. Earth Sci. Rev. 2005, 72, 89–111. [Google Scholar] [CrossRef]
  32. Pint, A.; Frenzel, P.; Horne, D.J.; Franke, J.; Daniel, T.; Burghardt, A.; Funai, B.; Lippold, K.; Daut, G.; Wennrich, V. Ostracoda from inland waterbodies with saline influence in Central Germany: Implications for palaeoenvironmental reconstruction. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2015, 419, 37–46. [Google Scholar] [CrossRef]
  33. Rossi, V.; Amorosi, A.; Sarti, G.; Mariotti, S. Late Quaternary multiple incised valley systems: An unusually well-preserved stratigraphic record of two interglacial valley-fill successions from the Arno Plain (northern Tuscany, Italy). Sedimentology 2017, 64, 1901–1928. [Google Scholar] [CrossRef]
  34. D’Orefice, M.; Bellotti, P.; Bertini, A.; Calderoni, G.; Censi Neri, P.; Di Bella, L.; Fiorenza, D.; Foresi, L.M.; Louvari, M.A.; Rainone, L.; et al. Holocene evolution of the Burano Paleo-Lagoon (Southern Tuscany, Italy). Water 2020, 12, 1007. [Google Scholar] [CrossRef] [Green Version]
  35. Donnini, F.; Dinelli, E.; Sangiorgi, F.; Fabbri, E. A biological and geochemical integrated approach to assess the environmental quality of a coastal lagoon (Ravenna, Italy). Environ. Int. 2007, 33, 919–928. [Google Scholar] [CrossRef]
  36. Migani, F.; Borghesi, F.; Dinelli, E. Geochemical characterization of surface sediments from the northern Adriatic wetlands around the Po river delta. Part I: Bulk composition and relation to local background. J. Geochem. Explor. 2015, 156, 72–88. [Google Scholar] [CrossRef]
  37. Zonta, R.; Cassin, D.; Pini, R.; Dominik, J. Assessment of heavy metal and As contamination in the surface sediments of Po delta lagoons (Italy). Estuar. Coast. Shelf Sci. 2019, 225, 106235. [Google Scholar] [CrossRef]
  38. Pitacco, V.; Mistri, M.; Ferrari, C.R.; Sfriso, A.; Sfriso, A.A.; Munari, C. Multiannual trend of micro-pollutants in sediments and benthic community response in a Mediterranean Lagoon (Sacca di Goro, Italy). Water 2020, 12, 1074. [Google Scholar] [CrossRef] [Green Version]
  39. Kusnerik, K.M.; Lockwood, R.; Grant, A.N. Using the fossil record to establish a baseline and recommendations for oyster mitigation in the Mid-Atlantic U.S. In Marine Conservation Paleobiology; Tyler, C.L., Schneider, C.L., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 75–103. ISBN 978-3-319-73795-9. [Google Scholar]
  40. Tomašových, A.; Albano, P.G.; Fuksi, T.; Gallmetzer, I.; Haselmair, A.; Kowalewski, M.; Nawrot, R.; Nerlović, V.; Scarponi, D.; Zuschin, M. Ecological regime shift preserved in the Anthropocene stratigraphic record. Proc. R. Soc. B Biol. Sci. 2020, 287, 20200695. [Google Scholar] [CrossRef]
  41. Kowalewski, M.; Wittmer, J.M.; Dexter, T.A.; Amorosi, A.; Scarponi, D. Differential responses of marine communities to natural and anthropogenic changes. Proc. R. Soc. B 2015, 282, 20142990. [Google Scholar] [CrossRef]
  42. Barmawidjaja, D.M.; van der Zwaan, G.J.; Jorissen, F.J.; Puskaric, S. 150 years of eutrophication in the northern Adriatic Sea: Evidence from a benthic foraminiferal record. Mar. Geol. 1995, 122, 367–384. [Google Scholar] [CrossRef]
  43. Pranovi, F.; Giovanardi, O.; Franceschini, G. Recolonization dynamics in areas disturbed by bottom fishing gears. In Recruitment, Colonization and Physical-Chemical Forcing in Marine Biological Systems; Baden, S., Pihl, L., Rosenberg, R., Strömberg, J.O., Svane, I., Tiselius, P., Eds.; Springer: Heidelberg, Germany, 1998; Volume 375/376, pp. 125–135. [Google Scholar]
  44. Gallmetzer, I.; Haselmair, A.; Tomašových, A.; Stachowitsch, M.; Zuschin, M. Responses of molluscan communities to centuries of human impact in the northern Adriatic Sea. PLoS ONE 2017, 12, e0180820. [Google Scholar] [CrossRef] [PubMed]
  45. Pitacco, V.; Mistri, M.; Munari, C. Long-term variability of macrobenthic community in a shallow coastal lagoon (Valli di Comacchio, northern Adriatic): Is community resistant to climate changes? Mar. Environ. Res. 2018, 137, 73–87. [Google Scholar] [CrossRef] [PubMed]
  46. Tomašových, A.; Gallmetzer, I.; Haselmair, A.; Kaufman, D.S.; Kralj, M.; Cassin, D.; Zonta, R.; Zuschin, M. Tracing the effects of eutrophication on molluscan communities in sediment cores: Outbreaks of an opportunistic species coincide with reduced bioturbation and high frequency of hypoxia in the Adriatic Sea. Paleobiology 2018, 44, 575–602. [Google Scholar] [CrossRef] [Green Version]
  47. Albano, P.G.; Gallmetzer, I.; Haselmair, A.; Tomašových, A.; Stachowitsch, M.; Zuschin, M. Historical ecology of a biological invasion: The interplay of eutrophication and pollution determines time lags in establishment and detection. Biol. Invasions 2018, 20, 1417–1430. [Google Scholar] [CrossRef] [Green Version]
  48. Seddon, A.W.R.; Froyd, C.A.; Witkowski, A.; Willis, K.J. A quantitative framework for analysis of regime shifts in a Galápagos coastal lagoon. Ecology 2014, 95, 3046–3055. [Google Scholar] [CrossRef] [Green Version]
  49. Stegner, M.A.; Ratajczak, Z.; Carpenter, S.R.; Williams, J.W. Inferring critical transitions in paleoecological time series with irregular sampling and variable time-averaging. Quat. Sci. Rev. 2019, 207, 49–63. [Google Scholar] [CrossRef]
  50. Dolven, J.K.; Alve, E.; Rygg, B.; Magnusson, J. Defining past ecological status and in situ reference conditions using benthic foraminifera: A case study from the Oslofjord, Norway. Ecol. Indic. 2013, 29, 219–233. [Google Scholar] [CrossRef]
  51. Campo, B.; Bruno, L.; Amorosi, A. Basin-scale stratigraphic correlation of late Pleistocene-Holocene (MIS 5e-MIS 1) strata across the rapidly subsiding Po Basin (northern Italy). Quat. Sci. Rev. 2020, 237, 106300. [Google Scholar] [CrossRef]
  52. Jones, G.D. Paralic. In Beaches and Coastal Geology; Schwartz, M., Ed.; Encyclopedia of Earth Science Series I; Springer: Boston, MA, USA, 1984; p. 631. ISBN 978-0-387-30843-2. [Google Scholar]
  53. Amorosi, A.; Colalongo, M.L.; Pasini, G.; Preti, D. Sedimentary response to Late Quaternary sea-level changes in the Romagna coastal plain (northern Italy). Sedimentology 1999, 46, 99–121. [Google Scholar] [CrossRef]
  54. Campo, B.; Amorosi, A.; Vaiani, S.C. Sequence stratigraphy and late Quaternary paleoenvironmental evolution of the Northern Adriatic coastal plain (Italy). Palaeogeogr. Palaeoclimatol. Palaeoecol. 2017, 466, 265–278. [Google Scholar] [CrossRef]
  55. Amorosi, A.; Colalongo, M.L.; Fiorini, F.; Fusco, F.; Pasini, G.; Vaiani, S.C.; Sarti, G. Palaeogeographic and palaeoclimatic evolution of the Po Plain from 150-ky core records. Glob. Planet. Chang. 2004, 40, 55–78. [Google Scholar] [CrossRef]
  56. Bruno, L.; Amorosi, A.; Severi, P.; Costagli, B. Late Quaternary aggradation rates and stratigraphic architecture of the southern Po Plain, Italy. Basin Res. 2017, 29, 234–248. [Google Scholar] [CrossRef]
  57. Amorosi, A.; Barbieri, G.; Bruno, L.; Campo, B.; Drexler, T.M.; Hong, W.; Rossi, V.; Sammartino, I.; Scarponi, D.; Vaiani, S.C.; et al. Three-fold nature of coastal progradation during the Holocene eustatic highstand, Po Plain, Italy—Close correspondence of stratal character with distribution patterns. Sedimentology 2019, 66, 3029–3052. [Google Scholar] [CrossRef]
  58. Amorosi, A.; Sammartino, I. Influence of sediment provenance on background values of potentially toxic metals from near-surface sediments of Po coastal plain (Italy). Int. J. Earth Sci. Geol. Rundsch. 2007, 96, 389–396. [Google Scholar] [CrossRef]
  59. Amorosi, A.; Centineo, M.C.; Dinelli, E.; Lucchini, F.; Tateo, F. Geochemical and mineralogical variations as indicators of provenance changes in Late Quaternary deposits of SE Po Plain. Sediment. Geol. 2002, 151, 273–292. [Google Scholar] [CrossRef]
  60. Borghesi, F.; Migani, F.; Dinelli, E. Geochemical characterization of surface sediments from the northern Adriatic wetlands around the Po River delta. Part II: Aqua regia results. J. Geochem. Explor. 2016, 169, 13–29. [Google Scholar] [CrossRef]
  61. Amorosi, A. Chromium and nickel as indicators of source-to-sink sediment transfer in a Holocene alluvial and coastal system (Po Plain, Italy). Sediment. Geol. 2012, 280, 260–269. [Google Scholar] [CrossRef]
  62. Cibin, U.; Severi, P.; Correggiari, A.; Roveri, M. Note Illustrative della Carta Geologica d’Italia alla Scala 1:50.000—Foglio 240–241 Forlì–Cervia; ISPRA—Servizio Geologico d’Italia: Roma, Italy, 2005. [Google Scholar]
  63. Rossi, V.; Vaiani, S.C. Benthic foraminiferal evidence of sediment supply changes and fluvial drainage reorganization in Holocene deposits of the Po Delta, Italy. Mar. Micropaleontol. 2008, 69, 106–118. [Google Scholar] [CrossRef]
  64. Ellis, B.F.; Messina, A.R. Catalogue of Foraminifera; Micropaleontology Press: New York, NY, USA, 1940. [Google Scholar]
  65. Ellis, B.F.; Messina, A.R. Catalogue of Ostracoda; The American Museum of Natural History Special Publications: New York, NY, USA, 1952. [Google Scholar]
  66. Boomer, I.; Horne, D.J.; Slipper, I. The use of ostracods in palaeoenvironmental studies or what can you do with an Ostracod shell? Paleontol. Soc. Pap. 2003, 9, 153–179. [Google Scholar] [CrossRef]
  67. Chao, A.; Shen, T.-J. Nonparametric estimation of Shannon’s index of diversity when there are unseen species in sample. Environ. Ecol. Stat. 2003, 10, 429–443. [Google Scholar] [CrossRef]
  68. Shannon, C.E.; Weaver, W.W. The Mathematical Theory of Communications; University of Illinois Press: Urbana, IL, USA, 1963. [Google Scholar]
  69. Hess, S.; Alve, E.; Andersen, T.J.; Joranger, T. Defining ecological reference conditions in naturally stressed environments—How difficult is it? Mar. Environ. Res. 2020, 156, 104885. [Google Scholar] [CrossRef] [PubMed]
  70. Bergamin, L.; Romano, E. Suitable sediment fraction for paleoenvironmental reconstruction and assessment of contaminated coastal areas based on benthic foraminifera: A case study from Augusta Harbour (Eastern Sicily, Italy). Ecol. Indic. 2016, 71, 66–78. [Google Scholar] [CrossRef]
  71. Gliozzi, E.; Grossi, F. Late Messinian lago-mare ostracod palaeoecology: A correspondence analysis approach. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2008, 264, 288–295. [Google Scholar] [CrossRef]
  72. Patzkowsky, M.E.; Holland, S.M. Stratigraphic Paleobiology: Understanding the Distribution of Fossil Taxa in Time and Space; University of Chicago Press: Chicago, IL, USA, 2012. [Google Scholar]
  73. Francescangeli, F.; Portela, M.; Armynot du Chatelet, E.; Billon, G.; Andersen, T.J.; Bouchet, V.M.P.; Trentesaux, A. Infilling of the Canche Estuary (eastern English Channel, France): Insight from benthic foraminifera and historical pictures. Mar. Micropaleontol. 2018, 142, 1–12. [Google Scholar] [CrossRef]
  74. The R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  75. Oksanen, J.; Blanchet, G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Package Version 2.5-6 2019. Available online: https://rdrr.io/cran/vegan/ (accessed on 10 May 2020).
  76. Meisch, C. Freshwater Ostracoda of Western and Central Europe (Süsswasserfauna von Mitteleuropa); Spektrum Akademischer Verlag: Heidelberg, Germany, 2000; Volume 3. [Google Scholar]
  77. Mazzini, I.; Anadon, P.; Barbieri, M.; Castorina, F.; Ferreli, L.; Gliozzi, E.; Mola, M.; Vittori, E. Late Quaternary sea-level changes along the Tyrrhenian coast near Orbetello (Tuscany, central Italy): Palaeoenvironmental reconstruction using ostracods. Mar. Micropaleontol. 1999, 37, 289–311. [Google Scholar] [CrossRef]
  78. Frenzel, P.; Keyser, D.; Viehberg, F.A. An illustrated key and (palaeo)ecological primer for Postglacial to Recent Ostracoda (Crustacea) of the Baltic Sea. Boreas 2010, 39, 567–575. [Google Scholar] [CrossRef]
  79. Salel, T.; Bruneton, H.; Lefèvre, D. Ostracods and environmental variability in lagoons and deltas along the north-western Mediterranean coast (Gulf of Lions, France and Ebro delta, Spain). Rev. Micropaléontol. 2016, 59, 425–444. [Google Scholar] [CrossRef]
  80. Bronk Ramsey, C. Bayesian Analysis of Radiocarbon Dates. Radiocarbon 2009, 51, 337–360. [Google Scholar] [CrossRef] [Green Version]
  81. Langone, L.; Asioli, A.; Correggiari, A.; Trincardi, F. Age-depth modelling through the late Quaternary deposits of the central Adriatic basin. Mem. Ist. Ital. Idrobiol. 1996, 55, 177–196. [Google Scholar]
  82. Coccioni, R. Benthic foraminifera as bioindicators of heavy metal pollution—A case study from the Goro Lagoon (Italy). In Environmental Micropaleontology: The Application of Microfossils to Environmental Geology; Martin, R.E., Ed.; Kluwer Academic/Plenum Publishers: New York, NY, USA, 2000; pp. 71–103. [Google Scholar]
  83. Debenay, J.-P.; Guillou, J.-J. Ecological transitions indicated by foraminiferal assemblages in paralic environments. Estuaries 2002, 25, 1107–1120. [Google Scholar] [CrossRef]
  84. Mojtahid, M.; Geslin, E.; Coynel, A.; Gorse, L.; Vella, C.; Davranche, A.; Zozzolo, L.; Blanchet, L.; Bénéteau, E.; Maillet, G. Spatial distribution of living (Rose Bengal stained) benthic foraminifera in the Loire estuary (western France). J. Sea Res. 2016, 118, 1–16. [Google Scholar] [CrossRef]
  85. Debenay, J.-P. Distribution trends of foraminiferal assemblages in paralic environments: A base for using foraminifera as bioindicators. In Environmental Micropaleontology: The Application of Microfossils to Environmental Geology; Martin, R.E., Ed.; Kluwer Academic/Plenum Publishers: New York, NY, USA, 2000; pp. 39–67. [Google Scholar]
  86. Melis, R.; Covelli, S. Distribution and morphological abnormalities of recent foraminifera in the Marano and Grado Lagoon (North Adriatic Sea, Italy). Medit. Mar. Sci. 2013, 14, 432–450. [Google Scholar] [CrossRef] [Green Version]
  87. Wukovits, J.; Bukenberger, P.; Enge, A.J.; Gerg, M.; Wanek, W.; Watzka, M.; Heinz, P. Food supply and size class depending variations in phytodetritus intake in the benthic foraminifer Ammonia tepida. Biol. Open 2018, 7, bio030056. [Google Scholar] [CrossRef] [Green Version]
  88. Vesper, B. Zur Morphologie und Ökologie von Cyprideis Torosa (Jones, 1850) (Crustacea, Ostracoda, Cytheridae) unter Besonderer Berücksichtigung Seiner Biometrie; Mitteilungen Hamburger Zoologisches Museum und Institut: Hamburg, Germany, 1972; pp. 21–77. [Google Scholar]
  89. Athersuch, J.; Horne, D.J.; Whittaker, J.E. Marine and Brackish Water Ostracods; Kermack, D.M., Barnes, R.S.K., Eds.; Synopses of the British Fauna (New Series); Brill E. J.: Leiden, The Netherlands, 1989. [Google Scholar]
  90. Arndt, S.; Jørgensen, B.B.; LaRowe, D.E.; Middelburg, J.J.; Pancost, R.D.; Regnier, P. Quantifying the degradation of organic matter in marine sediments: A review and synthesis. Earth Sci. Rev. 2013, 123, 53–86. [Google Scholar] [CrossRef]
  91. Pitacco, V.; Reizopoulou, S.; Sfriso, A.; Sfriso, A.; Mistri, M.; Munari, C. The difficulty of disentangling natural from anthropogenic forcing factors makes the evaluation of ecological quality problematic: A case study from Adriatic lagoons. Mar. Environ. Res. 2019, 150, 104756. [Google Scholar] [CrossRef]
  92. Pérez-Asensio, J.N.; Rodríguez-Ramírez, A. Benthic foraminiferal salinity index in marginal-marine environments: A case study from the Holocene Guadalquivir estuary, SW Spain. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2020, 560, 110021. [Google Scholar] [CrossRef]
  93. Carboni, M.G.; Succi, M.C.; Bergamin, L.; Bella, L.D.; Frezza, V.; Landini, B. Benthic foraminifera from two coastal lakes of southern Latium (Italy). Preliminary evaluation of environmental quality. Mar. Pollut. Bull. 2009, 59, 268–280. [Google Scholar] [CrossRef]
  94. Jorissen, F.J. Benthic Foraminifera from the Adriatic Sea; Principles of Phenotypic Variation; Utrecht Micropaleontological Bulletins: Utrecht, The Netherlands, 1988; Volume 37. [Google Scholar]
  95. Di Bella, L.; Bellotti, P.; Frezza, V.; Bergamin, L.; Carboni, M.G. Benthic foraminiferal assemblages of the imperial harbor of Claudius (Rome): Further paleoenvironmental and geoarcheological evidences. Holocene 2011, 21, 1245–1259. [Google Scholar] [CrossRef]
  96. Barbieri, G.; Rossi, V.; Vaiani, S.C.; Horton, B.P. Benthic Ostracoda and foraminifera from the North Adriatic Sea (Italy, Mediterranean Sea): A proxy for the depositional characterisation of river-influenced shelves. Mar. Micropaleontol. 2019, 153, 101772. [Google Scholar] [CrossRef]
  97. Montenegro, M.E.; Pugliese, N. Autecological remarks on the ostracod distribution in the Marano and Grado Lagoons (Northern Adriatic Sea, Italy). Boll. Soc. Paleontol. Ital. 1996, 3, 123–132. [Google Scholar]
  98. Tomašových, A.; Kidwell, S.M. Predicting the effects of increasing temporal scale on species composition, diversity, and rank-abundance distributions. Paleobiology 2010, 36, 672–695. [Google Scholar] [CrossRef]
  99. Johnson, R.G. Models and methods for analysis of the mode of formation of fossil assemblages. Bull. Geol. Soc. Am. 1960, 71, 1075–1086. [Google Scholar] [CrossRef]
  100. Davies, D.J.; Powell, E.N.; Stanton, R.J. Relative rates of shell dissolution and net sediment accumulation—A commentary: Can shell beds form by the gradual accumulation of biogenic debris on the sea floor? Lethaia 1989, 22, 207–212. [Google Scholar] [CrossRef]
  101. Martin, R.E. Taphonomy and temporal resolution of foraminiferal assemblages. In Modern Foraminifera; Sen Gupta, B.K., Ed.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1999; pp. 281–298. [Google Scholar]
  102. Martin, R.E.; Harris, M.S.; Liddell, W.D. Taphonomy and time-averaging of foraminiferal assemblages in Holocene tidal flat sediments, Bahia la Choya, Sonora, Mexico (northern Gulf of California). Mar. Micropaleontol. 1995, 26, 187–206. [Google Scholar] [CrossRef]
  103. Kowalewski, M. The youngest fossil record and conservation paleobiology: Holocene shells as eco-environmental recorders. In Conservation Paleobiology: Using the Past to Manage for the Future, Paleontological Society Short Course, Portland, OR, USA, 17 October 2009; Dietl, G.P., Flessa, K.W., Eds.; The Paleontological Society: New York, NY, USA, 2009; Volume 15, pp. 1–25. [Google Scholar]
  104. Horton, B.P.; Edwards, R.J. Quantifying Holocene Sea-Level Change Using Intertidal Foraminifera: Lessons from the British Isles; Special Publication No. 40.; Cushman Foundation for Foraminiferal Research: Lawrence, KS, USA, 2006. [Google Scholar]
  105. Levinton, J.S. The paleoecological significance of oppurtunistic species. Lethaia 1970, 3, 69–78. [Google Scholar] [CrossRef]
  106. Uffenorde, H. Ökologie und jahreszeitliche Verteilung rezenter benthonischer Ostracoden des Limski kanal bei Rovinj (nördliche Adria). Göttinger Arb. Geol. Paläont. 1972, 13, 142. [Google Scholar]
  107. Mazzini, I.; Rossi, V.; Da Prato, S.; Ruscito, V. Ostracods in archaeological sites along the Mediterranean coastlines: Three case studies from the Italian peninsula. In The Archaeological and Forensic Applications of Microfossils: A Deeper Understanding of Human History; Williams, M., Hill, T., Boomer, I., Wilkinson, I.P., Eds.; The Micropalaeontological Society Special Publications; Geological Society: London, UK, 2017; pp. 121–142. [Google Scholar]
  108. Gobert, S. Freshwater Ostracods as Palaeoevironmental Proxies in the Moervaart Depression, a Palaeolake in Sandy Flanders (NW Belgium). Master’s Thesis, University of Ghent, Ghent, Belgium, 2012. [Google Scholar]
  109. Ruiz, F.; González-Regalado, M.L.; Baceta, J.I.; Menegazzo-Vitturi, L.; Pistolato, M.; Rampazzo, G.; Molinaroli, E. Los ostrácodos actualesde la laguna de Venecia (NE de Italia). Geobios 2000, 33, 447–454. [Google Scholar] [CrossRef]
Figure 1. Location of cores 240 S5 and 240 S6. Maximum marine ingression from [51]. Satellite image: © 2020 Google.
Figure 1. Location of cores 240 S5 and 240 S6. Maximum marine ingression from [51]. Satellite image: © 2020 Google.
Water 12 03420 g001
Figure 2. Stratigraphy, vertical distribution of benthic foraminifer taxa, exp(H’bc), and related PaleoEcoQs, Foram−AMBI ecological groups, and Foram−AMBI values of core 240 S5. Ages are reported as calibrated years BP. Samples barren in benthic foraminifers are evidenced. Arrows are colored according to sample classification to improve the graphical visibility of selected samples.
Figure 2. Stratigraphy, vertical distribution of benthic foraminifer taxa, exp(H’bc), and related PaleoEcoQs, Foram−AMBI ecological groups, and Foram−AMBI values of core 240 S5. Ages are reported as calibrated years BP. Samples barren in benthic foraminifers are evidenced. Arrows are colored according to sample classification to improve the graphical visibility of selected samples.
Water 12 03420 g002
Figure 3. Detrended correspondence analysis (DCA) applied on the benthic foraminifer dataset of core 240 S5, showing samples and taxa colored in accordance with the Foram−AMBI ecological group. Total variance explained by each axis is reported. Exp(H’bc) and Foram−AMBI indices are plotted as additional variables.
Figure 3. Detrended correspondence analysis (DCA) applied on the benthic foraminifer dataset of core 240 S5, showing samples and taxa colored in accordance with the Foram−AMBI ecological group. Total variance explained by each axis is reported. Exp(H’bc) and Foram−AMBI indices are plotted as additional variables.
Water 12 03420 g003
Figure 4. Stratigraphy, vertical distribution of benthic foraminifer taxa, exp(H’bc), and related PaleoEcoQs, Foram−AMBI ecological groups, and Foram−AMBI values of core 240 S6. Ages are reported as calibrated years BP. Samples barren in benthic foraminifers, including rare or transported tests are evidenced. Arrows are colored according to sample classification to improve the graphical visibility of selected samples.
Figure 4. Stratigraphy, vertical distribution of benthic foraminifer taxa, exp(H’bc), and related PaleoEcoQs, Foram−AMBI ecological groups, and Foram−AMBI values of core 240 S6. Ages are reported as calibrated years BP. Samples barren in benthic foraminifers, including rare or transported tests are evidenced. Arrows are colored according to sample classification to improve the graphical visibility of selected samples.
Water 12 03420 g004
Figure 5. Detrended correspondence analysis (DCA) applied on the benthic foraminifer dataset of core 240 S6, showing samples and taxa. Samples are colored according to the facies association; taxa are colored in accordance with the Foram−AMBI ecological group. Total variance explained by each axis is reported. Exp(H’bc) and Foram−AMBI indices are plotted as additional variables.
Figure 5. Detrended correspondence analysis (DCA) applied on the benthic foraminifer dataset of core 240 S6, showing samples and taxa. Samples are colored according to the facies association; taxa are colored in accordance with the Foram−AMBI ecological group. Total variance explained by each axis is reported. Exp(H’bc) and Foram−AMBI indices are plotted as additional variables.
Water 12 03420 g005
Figure 6. Vertical distribution patterns of benthic ostracods in core 240 S5. Relative abundances of taxa and ecological groups, reported as abundance categories are shown. Samples with no or transported ostracod valves are evidenced.
Figure 6. Vertical distribution patterns of benthic ostracods in core 240 S5. Relative abundances of taxa and ecological groups, reported as abundance categories are shown. Samples with no or transported ostracod valves are evidenced.
Water 12 03420 g006
Figure 7. Vertical distribution patterns of benthic ostracods in core 240 S5. Relative abundances of taxa and ecological groups, reported as abundances categories, are shown. Samples with no or transported ostracod valves are evidenced.
Figure 7. Vertical distribution patterns of benthic ostracods in core 240 S5. Relative abundances of taxa and ecological groups, reported as abundances categories, are shown. Samples with no or transported ostracod valves are evidenced.
Water 12 03420 g007
Figure 8. Ternary plot of CaO, Al2O3, and LOI in accordance with [36] used to assess organic matter enrichment from core 240 S5 [25]. Points located below the average pelite (lower star)–carbonate (upper star) composition line suggest high organic matter content.
Figure 8. Ternary plot of CaO, Al2O3, and LOI in accordance with [36] used to assess organic matter enrichment from core 240 S5 [25]. Points located below the average pelite (lower star)–carbonate (upper star) composition line suggest high organic matter content.
Water 12 03420 g008
Table 1. List of radiocarbon ages from cores 240 S5 and 240 S6.
Table 1. List of radiocarbon ages from cores 240 S5 and 240 S6.
CoreCore DepthMaterialConventional AgeCalibrated 2σ AgeCalibrated Age Mean (yr BP)σDating LaboratoryReference
240 S59.5Organic clay5290 ± 1006300–57706080110ENEA (Bologna, Italy)[62]
240 S512.9Organic clay6770 ± 1057600–7230741090
240 S518.5Organic clay10650 ± 13012770–1214012540150
240 S523.4Organic clay12500 ± 22015470–1396014690390
240 S531.6Organic clay19050 ± 23023520–2245022980280
240 S621.55Mollusk shell7384 ± 458037–7856795845CEDAD (Lecce, Italy)[54]
240 S623.6Mollusk shell7358 ± 458013–7843793245
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barbieri, G.; Rossi, V.; Ghosh, A.; Vaiani, S.C. Conservation Paleobiology as a Tool to Define Reference Conditions in Naturally Stressed Transitional Settings: Micropaleontological Insights from the Holocene of the Po Coastal Plain (Italy). Water 2020, 12, 3420. https://doi.org/10.3390/w12123420

AMA Style

Barbieri G, Rossi V, Ghosh A, Vaiani SC. Conservation Paleobiology as a Tool to Define Reference Conditions in Naturally Stressed Transitional Settings: Micropaleontological Insights from the Holocene of the Po Coastal Plain (Italy). Water. 2020; 12(12):3420. https://doi.org/10.3390/w12123420

Chicago/Turabian Style

Barbieri, Giulia, Veronica Rossi, Anupam Ghosh, and Stefano Claudio Vaiani. 2020. "Conservation Paleobiology as a Tool to Define Reference Conditions in Naturally Stressed Transitional Settings: Micropaleontological Insights from the Holocene of the Po Coastal Plain (Italy)" Water 12, no. 12: 3420. https://doi.org/10.3390/w12123420

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