Next Article in Journal
Spatiotemporal Dynamics of Nitrogen Transport in the Qiandao Lake Basin, a Large Hilly Monsoon Basin of Southeastern China
Next Article in Special Issue
Improved Principal Component-Fuzzy Comprehensive Assessment Coupling Model for Urban River Water Quality: A Case Study in Chongqing, China
Previous Article in Journal
An Overview of Managed Aquifer Recharge in Brazil
Previous Article in Special Issue
Hydrochemical Characteristics, Controlling Factors, and Solute Sources of Streamflow and Groundwater in the Hei River Catchment, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiannual Trend of Micro-Pollutants in Sediments and Benthic Community Response in a Mediterranean Lagoon (Sacca di Goro, Italy)

1
National Institute of Biology, Marine Biology Station, Fornače 41, 6330 Piran, Slovenia
2
Department of Chemical and Pharmaceutical Sciences, University of Ferrara, Via Fossato di Mortara 17, 44121 Ferrara, Italy
3
ARPAE, Struttura Oceanografica Daphne, Viale Vespucci 2, 47042 Cesenatico (FC), Italy
4
Department of Environmental Sciences, Informatics and Statistics, University Ca’ Foscari, Via Torino 155, 30127 Mestre, Venice, Italy
*
Author to whom correspondence should be addressed.
Water 2020, 12(4), 1074; https://doi.org/10.3390/w12041074
Submission received: 24 February 2020 / Revised: 30 March 2020 / Accepted: 7 April 2020 / Published: 9 April 2020

Abstract

:
Long-term variations of ecological status in a Mediterranean coastal lagoon (Sacca di Goro, Northern Adriatic) were investigated, combining data on the concentration of surface sediment contaminants and on the structure of the macrobenthic community. The aim was to assess any amount of chemical contamination and check the response of the macrobenthic community to sediment contamination. Over the studied period, the sediments of the lagoon showed contamination by trace metals and organochlorine pesticides, with most of them exceeding the thresholds indicated by the Italian legislation in many samples. Contamination by polychlorinated dibenzodioxins and dibenzofurans (PCDD/Fs) and polycyclic aromatic hydrocarbons (PAHs) instead never exceeded the threshold. The ecological status based on the macrobenthic community, evaluated through biotic indices (AMBI and M-AMBI), fell below the Good/Moderate threshold in most samples. The results indicate a possible influence of toxic compounds in sediment on benthic organisms, but most of the variability shown by the macrobenthic community is probably due to other factors. The difficulty in establishing a cause/effect relationship was due to the co-occurrence and variability of various stressors (both natural and anthropogenic) and their interactions. The methods currently used for monitoring transitional waters thus seem insufficient to disentangle the effect of pollutants and other environmental variables on the benthos. Integrated approaches (e.g., bioaccumulation and toxicity tests) are thus needed for a more precise identification of the risk posed by a high concentration of pollutants in such environments.

1. Introduction

Contaminated sediments are a priority issue in aquatic systems impacted by anthropogenic use, industry, and urban development. Sediments act as a "sink" for contaminants, making the analytical determination of their concentrations easier with respect to the water column and providing time-integrated information about the ecosystem health. The evaluation of temporal changes in sediment contamination is crucial to answer questions about the state of the system, whether it is improving in time or to determine if remediation activities have benefited the system [1]. Unlike in the marine environment, where long-term studies are more frequent [2], coastal lagoon studies on the effect of contaminants on the biota are often a snapshot in time, with very few samples in a short study period. In such studies, the influence of temporal variation cannot be evaluated. The knowledge on temporal variation in contaminant availability is desirable, as recommended by the Water Framework Directive [3], implemented in Italy through the National Act 260/10, which focuses mostly on trend analysis of time series data obtained in sediment monitoring programs.
Coastal lagoons are subjected to multiple stressors, from the high rate of natural dynamic changes (fluctuating temperature, salinity, and instability of sediment surface) to degradation due to human population growth and economic development [4]. Together with climatic changes, the major threats are related to the increase in pollutants produced by anthropogenic activities that influence both the geomorphology and geochemical backgrounds [5]. Coastal lagoons receive large amounts of pollutants derived from urban, agricultural, and industrial effluents, and can accumulate them, temporarily or permanently [6]. Interpreting the effects of anthropogenic contamination in coastal lagoons is often complex and confounding, since dynamic physical, chemical, and geologic conditions may interfere with the assessment of anthropogenic impacts on lagoon biotic integrity [7]. Macrobenthic communities play an important functional role in lagoons and other aquatic ecosystems. They alter geochemical conditions at the sediment–water interface, promote decomposition and nutrient cycling, and transfer energy to other food web components. Community composition and structure are influenced by an array of anthropogenic and non-anthropogenic stressors, and direct causal relationships are not evident [7]. The evaluation of the relative influence and interaction of different factors, including the effects of sediment-associated contaminants, poses a major challenge. The analysis of long-term data series could be a powerful tool to analyze those dynamics, and to assess the ecological risk posed by sediment pollutants in transitional waters.
The Sacca di Goro, a large Mediterranean lagoon, provides a good case study for examining spatial/temporal trends of contaminated sediments. It is a natural heritage (wetland of international importance according to the Ramsar Convention, UNESCO heritage site, and the Site of Interest of the European Union) and is part of one of the largest Mediterranean deltaic systems (Po River Delta). This gives the study a supra-local perspective: the need to better understand field-observed contaminant effect on the biota in coastal lagoons is relevant for many places on the globe. Chemical contamination from industry, urban, and agricultural runoff of the Po Valley are known to affect coastal lagoon ecosystems at the Po Delta [8,9]. Many of those micro-pollutants are known to have adverse biological effects on benthic organisms [10,11]. In the present work, the ecological status (ES) of the lagoon was assessed through the combined analysis of chemical elements (the presence of anthropogenic pollutants in sediments) and biological elements (the structure and diversity of macrobenthic community as a proxy for healthy ecosystem functioning). To assess the level of chemical pollution, we targeted major anthropogenic micro-pollutants required by Italian legislation (National Act 260/10), such as: organochlorine pesticides (OCPs), trace elements (TEs), polychlorinated dibenzodioxins (PCDDs), dibenzofurans (PCDFs), and polycyclic aromatic hydrocarbons (PAHs). For the assessment of the integrity of the macrobenthic community, structural (richness and diversity) and biotic indices (AMBI and M-AMBI) were used (as required by Act 260/10). The aim of the present work was threefold: (i) to evaluate the level and temporal trend of sediment contamination of the lagoon and to compare contaminants’ concentration with national and international thresholds; (ii) to test if the macrobenthic community showed a response to multiannual sediment contamination; and (iii) to combine the temporal trend of the ecological quality based on chemical and biological indices of the Sacca di Goro.

2. Materials and Methods

2.1. Study Area

The Sacca di Goro is a shallow coastal lagoon located in the southern part of the Po River Delta, with a surface area of 2600 ha and an average depth of about 1.2–1.5 m. Freshwater inputs come mainly from the Po di Volano river (about 3.5 × 108 m3/y) and from canals with similar flows (2.0–5.5 × 107 m3/y), named Giralda, Romanina, and Canal Bianco, that flow directly into the western part of the lagoon, and from the Po di Goro deltaic branch, which is artificially regulated through a dam in its eastern part. Marine inflows, varying according to the tidal dynamics, come from two mouths connecting the lagoon to the Northern Adriatic Sea. Consequently, the area is characterized by large daily fluctuations, linked to the height of the tides, and seasonal fluctuations in the physicochemical parameters of the lagoon’s water (temperatures: 2–33°C, salinity: 6–30 psu, and pH: 7–8.8; [12]. Most of the lagoon’s floor consists of silty-clay sediments, carried to the sea by rivers, but there are also sandy areas, particularly near the lagoon mouths and behind the spit and the barrier island [12]. Sampling was overall performed from 2004 to 2010 at 4 sampling sites (Figure 1), representative of the different hydrological characteristics of the lagoon. Site GOR1 (44°49’39.21″ N, 12°16’54.68″ E) was located in the western part. It was influenced mainly by freshwater discharged from Po di Volano and Giralda and was therefore characterized by variable salinity. Sites GOR2 (44°49’07.19″ N, 12°19’11.93″ E) and GOR4 (44°49’47.11″ N, 12°18’18.30″ E) were located in the central area and were influenced more by tidal exchange, and showed higher salinity. Site GOR3 (44°48’43.67″ N, 12°20’20.09″ E) was located in the eastern part, at the edge of the Valle di Gorino, the shallower and most confined part of the lagoon.

2.2. Sample Collection and Chemical Analysis Methods

One sediment sample (top 5 cm) was collected seasonally at each site by means of a box corer, distributed into separate pre-labeled, acid-washed, and solvent-rinsed glass jars, stored on board a ship in coolers filled with ice and then stored at −20 °C immediately upon return to the laboratory. All samples were analyzed for grain size and organic matter. Particle size distribution was determined by ASTM D 422–63 method. Organic matter was estimated by loss on ignition (LOI: 105 °C for 24 h followed by ashing at 550 °C for 4 h). Trace elements (Cr, Ni, Pb, Hg, Cd, and As) were analyzed by inductively coupled plasma mass spectrometry (ICP-MS) according to the EPA method 6020B [13]. Organochlorine pesticides (OCPs) were extracted and cleaned using a QuEChERS method [14], followed by gas chromatography mass spectrometry (GC-MS) analysis. This QuEChERS extraction method involves mixing the sample with acetonitrile and permits the salt out liquid–liquid partitioning step using anhydrous MgSO4 and sodium acetate. After shaking and centrifugation, cleanup is done by dispersive solid phase extraction (d-SPE). Sulfur cleanup followed the EPA method 3660B (USEPA, 1996). Analyses of PCDD/Fs were performed following the EPA method 1613B and quantification through high-resolution gas chromatography high-resolution mass spectrometer (HRGC-HRMS). Briefly, sediment samples were extracted in a Soxhlet apparatus using 300 mL of 10% acetone in toluene for 20 h. The extract was divided into portions prior to the purification for PCDD/Fs, dioxin-like PCBs, and PBDEs. The extracts for PCDD/F, PCBs, and PBDEs analysis were cleaned by passing through a multi-layer silica gel column. Analytical accuracy was achieved by the use of blanks and certified reference material for sediment included in all series of analysis. The analyses of polycyclic aromatic hydrocarbons (PAHs) was performed following EPA methods 8270D [15], and quantification was performed through high-resolution gas chromatography high-resolution mass spectrometer (HRGC-HRMS).
Three replicates were collected seasonally for macrofaunal community analysis using a 4.0 L Van Veen grab at three sampling sites (GOR1, GOR3, and GOR4; the site GOR2 was not sampled for macrofauna). These samples were sieved at 0.5 mm and preserved in 8% formalin. Benthic organisms were sorted and identified to the species level.

2.3. Sediment and Benthic Data Analysis

In order to obtain a global picture and a temporal trend of the general status of the lagoon, sampling sites were used as replicates. Given the hydro-morphological differences within the lagoon, Trellis graphs, using Lattice package for R [16], were used as explorative analysis to check whether the trend of each variable was consistent among sampling sites. If the prerequisite of no differences among sampling sites was not met, the temporal trend was analysed separately for each site. To illustrate the main characteristics of sediments in Sacca di Goro, descriptive statistics were computed. In order to define the chemical status of the lagoon, chemical data were confronted to the existing national thresholds, when available, according to the Italian Law (National Act 260/10). This law fixes the threshold concentration of harmful elements and compounds to discriminate between Good and Moderate status, as requested by the Water Framework Directive (WFD). For trace metals, local background values should be taken into consideration, together with national thresholds, so values calculated from Holocene lagoon deposits sampled from several deep cores drilled in the Po river coastal plain [17,18] were considered as well, as suggested by the WFD. Moreover, in order to evaluate the potential risk for benthic organisms, chemical concentrations were compared with thresholds defined by Sediment quality guidelines [19], which define two thresholds for a chemical: Effects Range-Low (ERL) and Effects Range-Median (ERM). These values defined three ranges, rarely (below the ERL), occasionally (above the ERL, but below ERM), and frequently (above the ERM) associated with adverse effects on acquatic organisms. To test for differences in sediment characteristics and contaminant concentrations among years, a non-parametric Friedman test [20] was used. When significant differences were encountered, a post-hoc test according to Nemenyi [21] was also carried out. If a linear increase or decrease was observed, its significance was tested with Spearman rank correlation coefficient (rs).
For benthic communities, richness (S), Shannon index (H), AMBI index [22], M-AMBI index [23], and their Ecological Quality Ratios (EQR) were calculated for each sample using the free software (http://www.azti.es, v.5.0). Then, their annual averages were calculated. AMBI and its multivariate version M-AMBI are probably the most widely used biotic indices all over the world [24]. The percentage of invertebrates belonging to the different ecological groups (EGs) according to AMBI library was also calculated for each sample and then annual averages were calculated. The reference values (Hlog2 = 3.4, S = 28 and AMBI = 2.14) and thresholds for ecological status (ES) classification were those reported by Italian legislation (National Act 260/10), for the typology “M-AT-2, oligo/meso/poly microtidal lagoon”.
Differences in macrofaunal richness, diversity, AMBI and M-AMBI values, and EGs among the different years, were investigated through permutational multivariate analysis of variance (PERMANOVA) [25], using a single factor design with 7 levels (corresponding to years of sampling) and the “unrestricted permutation of row data”, as recommended in [25].
In order to explore the relations between contaminants and benthic indices, the sites were considered separately. The Spearman rank correlation coefficient (rs) was used in order to check for correlations: (i) among different pollutants, (ii) among sediment texture and pollutants, and (iii) among biotic index and abiotic parameters (sediment texture and pollutants). Distance-based linear models (DISTLM) were used to test and quantify the variation in macrobenthic community explained by abiotic variables. DISTLM does a partitioning of variation in a data cloud described by a resemblance matrix, according to a multiple regression model [26]. First, the marginal test shows the amount of variation explained by each variable when taken alone, ignoring all other variables; then, with the sequential test, individual variables are selected with a forward procedure (R2 selection criteria), and the proportion of explained variation attributed to each variable that is added to the model is a function of other variables already present in the model [25]. Chemical and sediment variables were analysed separately and selected in order to meet the assumption required by DISTLM analysis: the number of samples higher than the number of variables [25]. Abiotic variables were chosen in order to avoid multi-collinearity, keeping variables that showed a significant relation with biotic variables. Matrices with (i) sediment and (ii) chemical variables for year for site were used as predictor variables. As response variables, the following matrices were used: (i) ‘taxa/abundance’ matrix; (ii) ‘richness’ matrix; (iii) ‘diversity’ matrix; (iv) ‘AMBI’ matrix; (v) ‘M-AMBI’ matrix; and (vi) ‘EGs’ matrix. Bray Curtis similarity was used for ‘taxa’ and ‘EGs’ matrices, while Euclidean distance was used for ‘richness’, ‘diversity’, ‘AMBI’, and ‘M-AMBI’ matrices. All those calculations were performed with PRIMER v6 + PERMANOVA software package [25,27]. For a more detailed analysis of the relationships between chemical pollutants and macrobenthic community, factors showing a significant correlation according to the marginal test (DISTLM analysis) where analysed with the Spearman rank correlation coefficient (rs) and plotted to check differences among sites. A p-value < 0.05 was chosen as a significant threshold. Those analyses were run using R version 2.4.0 [28].

3. Results

3.1. Sediment Characteristics

Sediments of Sacca di Goro lagoon at all four sampling sites were typically sandy mud (from 50% to 75% of silt). Not significant differences were observed among years in the percentage of OM, sand, silt, and clay in the sediment (Friedman test, p < 0.05) (Figure 2). Organic matter (OM) ranged between 37.6 ± 5.4 mg kg−1 in 2007 and 50.4 ± 8.6 mg kg−1 in 2009 (Figure 2).

3.2. Trace Elements

Sacca di Goro lagoon showed signs of contamination by TEs from 2004 to 2010 (Table 1), particularly by Pb and Ni (Table 1), with average concentrations exceeding both their EQS and their background value in most of the studied years (Table 1). Concentrations of Cr and Cd exceeded their EQS almost every year as well, but Cr concentration was always below its background level (Table 1). Average concentrations of Hg reached its EQS only once in 2008, but exceeded its background values in most samples, while mean annual concentrations of As were always below the EQS values (Table 1). Mean annual concentrations of Ni exceeded both the ERL and ERM threshold every year (Table 1). As and Cr exceeded their ERL but not their ERM thresholds (Table 1). The concentrations of Pb and Cd were always below their ERL threshold, and Hg exceeded its ERL only in 2006, 2008, and 2009. No linear trend was observed during the studied period (rs, p > 0.05 for all TEs), but Ni, Cr, and Hg showed significant differences among the studied years (Friedman test, p < 0.05). The lowest values of Ni and Cr were observed in 2004 and 2005 and the highest values in 2007 (Figure 3). The highest value of Hg was observed in 2008, whereas in 2004 it was under the limit of detection. The concentration of As, Pb, and Cd did not differ among years (Friedman test, p > 0.05).

3.3. Organochlorine Pesticides (OCPs)

A moderate contamination by OCPs was observed, mainly related to dichlorodiphenyltrichloroethane (DDT) metabolites. Total concentrations of OCPs ranged from 0.09 ± 0.18 µg/kg in 2004 to 4.14 ± 2.44 µg/kg in 2007, showing two main peaks in 2007 and 2010 (ΣOCPs, Table 2). DDT metabolites (DDE and DDD), were recorded almost every year (Table 2). DDE was the metabolite present with the highest concentrations almost every year and exceeded its national threshold in 2007 (1.91 ± 1.04 µg/kg) and 2010 (2.47 ± 1.34 µg/kg). The only exception was in 2008 when DDD concentration was higher (1.35 ± 2.20 µg/kg) and exceeded its threshold value. The concentration of DDE in 2010 and DDD in 2008 exceeded even the 20% of tolerance permitted by Italian legislation. Some samples exceeded the ERL value for DDE, particularly in 2010, but not the ERM value (Table 2). DDE also showed a slight increase (rs = 0.37, p < 0.05) during the studied period. Conversely, concentrations of DDT were under the limit of detection every year, with the exception of 2007, when it was recorded with values not exceeding the national EQS. Other OCPs, namely Hexachlorobenzene, Simazine, Terbuthylazine, and Alaclor, were under the level of detection every year except in 2005 and 2007, with concentrations not exceeding their national EQS. Therefore, DDT metabolites were the only OCPs preventing the achievement of Good chemical status in 2007, 2008, and 2010.

3.4. Polychlorinated Dibenzodioxins (PCDDs) and Dibenzofurans (PCDFs)

PCDDs and PCDFs were measured from 2007 to 2010 (Table 3). Total PCDD/Fs ranged from 128.15 ng/kg in 2009 to 178.46 ng/kg in 2007. Concentrations of PCDD/Fs did not vary significantly during the period (Friedman testKW p > 0.05). Octachlorodibenzo-p-dioxin (O8CDD) was predominant every year, with concentrations between 88.4 ng/kg in 2009 and 123.9 ng/kg in 2007. The mean annual toxic equivalents (TEQs) for PCDD/Fs varied from 0.32 ng/kg in 2009 to 0.37 ng/kg in 2007, never exceeding the threshold (2 ng/kg) imposed by Italian legislation (Table 3).

3.5. Polycyclic Aromatic Hydrocarbons (PAHs)

Contamination by PAHs was observed every year (Figure 3). PAHs’ concentrations varied during the period (Friedman test p < 0.05). The lowest values of the sum of PAH (Σ16PAHs) were observed in 2004 (4.05 ± 5.7 µg/kg) and 2005 (4.05 ± 52.6 µg/kg). PAH concentrations markedly increased in 2006 (349.92 ± 144.4 µg/kg), reaching a peak (370.8 ± 151.1 µg/kg) in 2007 (Figure 3). The sum of PAHs never exceeded the threshold imposed by Italian legislation (Figure 3) for a Good chemical status, and was markedly below the threshold established by the Sediment Quality Guidelines (ERL = 4022 ng/g; ERM = 44,792 ng/g). Four-ring PAHs were the most abundant every year (38.9%–54.3% of total PAHs), followed by five-ring PAHs (9.8%–29.7%) and three-rings PAHs (15.7%–23.6%). Four-ring PAHs were mainly represented by pyrene (Py) and benzofluoranthene (BFl), five-rings were mainly benzopyrene (BaP) and indenopyrene (InPy), and three-ring PAHs were mainly fluoranthene (Fl) and phenanthrene (Phe). Two-ring PAHs (0.2%–9.9%) and six-ring PAHs (2.5%–13.5%) accounted for small percentages of the total PAHs. Concentrations of high molecular weight (HMW) PAHs (4- to 6-ring), clearly predominated along the studied period (from 66.5% to 83% of total PAHs) compared to those of low molecular weight (LMW) PAHs (2- to 3-ring). The LMW/HMW ratio ranged from 0.2 to 0.5, the Phe/Ant ratio was low (always <30), and the Flu/Py ratio was high (always >0.8, except in 2005 and 2009).

3.6. Ecological Status Based on Macrobenthic Community

The macrobenthic community was dominated by tolerant species (EGIII) along the entire studied period (Figure 4), with percentages ranging from 61.6% in 2008 to 73.4% in 2006. Among them, the most abundant species was the polychaete Streblospio shrubsolii (32.6%–75.2%). The percentages of opportunistic species were more variable but generally quite high. Second order opportunistic species (EGIV) ranged from 3.9% in 2006 to 29.3% in 2008, and were mainly represented by the polychaete Polydora ciliata (2.4%–32.4%). First order opportunistic species (EGV) ranged from 0.8% in 2006 to 23.4% in 2004 and were almost entirely represented by the polychaete Capitella capitata (0.04%–6.6%) and by oligochaetes (0%–17.2%). Contrastingly, the percentages of sensitive species (EGI) were generally low, ranging from 1.1% in 2008 to 10.9% in 2010, and they were mainly represented by amphipods, such as Ampelisca sarsi (0%–1.5%), Gammarus aequicauda (0%–4%), and Melita palmata (0.1%–2.3%). The percentages of indifferent species (EGII) were generally even lower (0.5%–4.1%), with the only exception of the high percentage observed in 2006 (18.3%), mainly due to the polychaete Salvatoria clavata (13.6%).
According to the AMBI index, ecological status was Good/High only in 2006, 2009 and 2010 (AMBI < 3.3), while M-AMBI index exceeded the threshold for Good Ecological status (M-AMBI > 0.71) only in 2006 and 2009. Taxa/abundance matrix varied significantly among years (PERMANOVA, p > 0.05). The percentage of each ecological group (EG), richness, diversity, AMBI and M-AMBI average values did not showed significant differences among years (PERMANOVA, p > 0.05). The percentage of sensitive species (EGI), S, H and M-AMBI showed an increase of within-year variation along the studied period (increasing SD, Figure 4), because the temporal trend of biotic indices varied among sites, as well. Oscillations between Good and Moderate ecological statuses were observed at each site, for both the AMBI and M-AMBI index. The temporal trend deviated more or less markedly from linearity in all cases and differed among sites. The M-AMBI index showed a general improving trend of ecological condition with increasing M-AMBI values at sites GOR3 and GOR4 but a general worsening trend, with decreasing M-AMBI values, at site GOR1. The AMBI index showed a similar trend: a general increasing in AMBI values (corresponding to worsening ecological conditions) at GOR1, and a decreasing at GOR3 (corresponding to improving ecological conditions), whereas, at site GOR4, values did not vary so much around the Good/Moderate threshold. S showed an increasing trend at site GOR4, but a decreasing one at sites GOR1 and GOR3. H values were generally stable at site GOR3, while a marked non-linear trend was observed at sites GOR1, with the highest values in 2007, and at GOR4, with the highest value in 2009.

3.7. Combined Trend of Abiotic and Biotic Variables

Collinearity was observed between some pollutants: Ni was positively correlated with Cr and As and negatively correlated with Hg (Table 4). A positive correlation was observed also between Cr and Pb, between PCDFs and PCDDs (Table 4). Few pollutants showed a correlation with sediment texture (Table 4). Organic matter (OM) showed a significant negative relationship with concentrations of Ni and Cr, and silt was negatively correlated with DDE (Table 4).
Collinearity was observed also among biotic indices, with a positive correlation between M-AMBI and H, M-AMBI and S, and a negative correlation between M-AMBI and AMBI, and between AMBI and H (Table 4). Among the biotic indices, only AMBI values were significantly related with sediment texture, and only with clay (Table 4). Among the pollutants, only As and DDE were negatively related with species richness (S) (Table 4). DISTLM analysis showed no correlation of macrobenthic community (‘taxa/abundance’ matrix) and biotic indices with sediment variables, i.e., sand, silt, clay and OM (DISTLM, p > 0.05, marginal test). Ten variables (PCDD/Fs, PAH tot, DDE, DDD, Ni, Hg, Pb, Cd, Cr, As) were chosen among chemical pollutants, to meet the assumptions required by DISTLIM analysis. Those chemical variables did not show significant relationships with macrobenthic community in terms of the ‘taxa/abundance’ matrix (DISTLM, p > 0.05, marginal test). Considering biotic indices, species richness (S) was the only one showing a correlation with pollutants, namely As and DDE (DISTLM, p < 0.05, marginal test). The highest percentage of variability within the S index was explained by As alone (52% of variability, sequential test). Site GOR1 was the most contaminated, with concentrations of both As and DDE exceeding the threshold for good chemical status in most samples, and with the lowest values of S (Figure 5A). Conversely, in most samples at sites GOR3 and GOR4, concentrations of As and DDE did not exceed their threshold and S values were higher (Figure 5B).

4. Discussion

4.1. Trace Elements Contamination

Evidences of contamination by TEs were observed in sediments of Sacca di Goro from 2004 to 2010. Concentrations of most TEs exceeded their national EQS in many samples, with the exception of Hg and As, showing annual means always below their EQS values. According to the WFD and Italian Legislation, for a correct estimation of the chemical status, the local background levels of trace metals of the area should be considered as well. Cr and Ni, even exceeding their national EQS values in most samples, were below background levels [17,18] in most samples. Similar situations with high values of Cr and Ni, were reported also from other lagoons of the Po River Delta: Pialassa Baiona (Ni: 35–57 mg/kg; Cr: 53–142 mg/kg; [18] and Comacchio lagoon (Ni: 59–75.4 mg/kg; Cr: 72.5–97.5 mg/kg [28]. Sediments of the Po River Delta are known to be enriched in Cr and Ni, due to geological reasons [9,17], and such concentrations are attributed to sediment provenance from ultramorfic source rocks of the Po catchment basin [29]. Therefore, such an enrichment of Cr and Ni in Sacca di Goro lagoon was likely of natural origin, as already pointed out by previous works [5]. Conversely, the high levels of Pb recorded in Sacca di Goro sediments exceeded both its national EQS and its background level, suggesting an anthropogenic enrichment. Concentrations of Pb were higher than those reported for other areas of the Po River Delta, such as Comacchio lagoon (21–28 mg/kg, [30], Pialassa Baiona lagoon (17 ± 9.1 mg/kg, [18], and Pialassa dei Piomboni lagoon (mean 18 mg/kg [17]). High levels of Pb in Sacca di Goro, were reported also by [5]. Concentrations of Hg also suggest an anthropogenic enrichment, even if this element did not exceeded threshold for Good chemical status. In Sacca di Goro lagoon, the mean annual concentrations of Pb, As, and Cr exceeded the ERL but not the ERM threshold defined by the Sediment quality guidelines, indicating a "Possible-effects" within which effects on benthic organisms would occasionally occur [19]. Cd and Hg concentrations were mostly below their ERL threshold, so an effect on benthic organisms was unlikely, whereas Ni exceeded its ERM threshold every year, indicating potentially adverse biological effects.

4.2. Organochlorine Pesticides Contamination

We found the presence of OCP in the Sacca di Goro: DDE was the dominant metabolite, persistent with increasing concentration along the entire studied period. Its mean concentration was comparable with values reported for other Mediterranean lagoons, such as Comacchio lagoon (0.68 ± 0.42 SD, [30], and Berre lagoon (DDE: 1.7 μg/kg; [31], but it exceeded its EQS in some samples. The watershed of Sacca di Goro is predominantly agricultural, with only a few main crop typologies: corn and wheat, rice, sugar beets, soybean, and vegetables [8]. The main freshwater input and therefore the major input of pollutants of agricultural origin in the lagoon comes from Po di Volano canal. The use of dichlorodiphenyltrichloroethane (DDT) in agriculture was banned by the European Union in the 1970s but is still very common in many ecosystems because of its high persistence [32]. Both DDE and DDD existed as by-products in commercial DDT formulations, and both may be formed by environmental degradation of DDT [33]. The dominance of DDE, with low concentrations of DDT, in Sacca di Goro indicates that the contamination was likely not recent. One among the probable causes for the persistence of those contaminants is the storage in the watershed soils and/or aquatic sediments and delayed delivery [8]. DDT compounds are of particular concerned, since their toxicity for benthic organisms has been proved by several studies (e.g., for amphipods) [34,35]. According to the Sediment quality guidelines [31], some samples in Sacca di Goro exceeded the ERL value for DDE, even if not the ERM value; therefore, an effect on benthic organisms was probable.

4.3. Polychlorinated Dibenzodioxins and Dibenzofurans Contamination

Concentrations of PCDDs and PCDFs in Sacca di Goro were comparable with those recorded for other Adriatic lagoons, such as Pialassa Baiona (PCDDs: 96–19,641 ng/kg and PCDFs: 16–1313 ng/kg) [18], Venice lagoon (PCDDs: 16–13,642 ng/kg and PCDFs: 49–126,561 ng/kg) [36], Comacchio lagoon (PCDD/Fs: 15.31 ng/kg–52.76 ng/kg) [30], and other locations, such as the Spanish northern Atlantic coast (PCDD/Fs: 0.15–3.99 ng/kg) [37], and the intertidal zone of the North Sea (PCDD/Fs: 0.124–3.156 ng/kg) [38]. Differences in concentration could be explained with the different distances from industrial and/or thermal power plants. In fact, PCDD/Fs do not occur naturally, but they are formed as by-products during the production of other chemicals or during combustion and incineration processes [39]. They are highly toxic, and very persistent, so they are now globally distributed. These compounds have low solubility and are highly hydrophobic; thus, they tend to accumulate in organisms and to adsorb to particles [39]. In Sacca di Goro lagoon, the mean annual TEQs for PCDD/Fs was always well below the threshold imposed by Italian legislation and therefore no adverse effects on benthic organisms were expected.

4.4. Polycyclic Aromatic Hydrocarbons Contamination

Sediments of Sacca di Goro were contaminated by PAHs along the studied period, but total PAHs concentrations (range Σ16PAHs: 4.05–370.8 µg/kg) never exceeded the threshold for a Good chemical status. Total concentration of PAHs was higher than those reported from other Adriatic lagoons, such as Lesina (8.51–70.41 ng/g) and Varano (6.61–55.06 ng/g) lagoons [40] but lower than Comacchio lagoon (9.1–2100.2 µg/kg) [30], Pialassa Baiona lagoon (2500–120000 ng/g) [41], Grado and Marano lagoon (<limit of detection-1056 mg/kg) [42], and other Mediterranean lagoons, such as Faro and Ganzirri lake (74–5755 ng/g) [43] and Stagnone lagoon (72–18381 ng/g) [44]. PAHs are an important class of persistent organic pollutants, introduced into the environment mainly via anthropogenic input [45]. Their potential sources could be estimated analysing the relative concentrations of PAHs by aromatic groups (two-ring, three-ring, four-ring, five-ring, and six-ring PAHs). The dominance of High Molecular Weight PAHs (=low LMW/HMW ratio), together with the low Phe/Ant ratio and high Flu/Py ratio observed suggested that the main sources of PAHs in Sacca di Goro were combustion processes. The predominant pyrolitic origin of PAHs in sediments has also been commonly observed in other Mediterranean lagoons [6,42]. PAHs could result from a large number of possible sources, such as industrial wastewater, sewage, road runoff/street dust, and petroleum-related activities [6], but a more accurate identification of the origins of PAH is beyond the scope of this work. Some PAHs and their metabolites are recognized to be carcinogenic; therefore, different tools have been proposed in order to predict their potential toxicity. Total PAHs’ concentrations in Sacca di Goro were markedly below thresholds established by the Sediment Quality Guidelines [19] and therefore no adverse effects on benthic organisms were expected.

4.5. Ecological Status Based on Macrobenthic Community

The macrobenthic community in Sacca di Goro was dominated by tolerant species, and values of M-AMBI and AMBI indices resulted below the threshold for Good ecological status in most samples. Low diversity and high abundances of tolerant species is a common condition of transitional waters, which are known to be harsh environments because of shallow waters, limited water renewal and marked seasonal variations of oxygen, temperature, and salinity. Those naturally stressed conditions coincide with anthropogenic disturbance according to the models based on Pearson–Rosenberg paradigm [46], potentially reducing the robustness of the AMBI index in transitional ecosystems [47,48]. Annual averages of EGs, richness, and diversity did not show a general clear pattern along the studied period, because of the marked differences among sampling sites. Those differences lead to the increase in within-year variation along the studied period (SD). The general trend of most biotic indices differed among sampling sites, suggesting a general tendency towards an improvement of ecological conditions at sites GOR3 and GOR4 (increasing M-AMBI), and a tendency towards a deterioration at site GOR1 (decreasing M-AMBI). At site GOR4, the improvement of ecological conditions was mainly related to an increase in richness and diversity, whereas, at site GOR3, it was mainly related to a decrease of AMBI index and therefore to a proportional decrease in tolerant and opportunistic species. At GOR1, the worsening of ecological condition was related to the decrease in richness and abundance and an increase in AMBI values. Those differences represent a potential impediment for a correct evaluation of the ES, since the positive trend of two sites could mask the negative trend in one site. An incorrect estimation of the temporal trend of ecological status could limit the possibility of identifying drivers of changes and consequently plan appropriate restoration measures.

4.6. Coupling Abiotic and Biotic Variables

The fluctuations observed in the present work in terms of both pollutant concentrations and biotic indices reflected the marked fluctuations of environmental parameters characterizing those ecosystems. In this framework, the relationship between biotic and physico-chemical variables is always complex, and disentangling the effect of natural and anthropogenic disturbances is a major challenge [49,50].
It is well known that salinity and confinement (defined as the degree of connection to the sea) [51] are the main structuring factors of the biotic lagoon communities; however, this study suggests that other factors, such as some contaminants, may also play a role. Macrobenthic community responded, in terms of species richness, to high levels of two micro-pollutants: As and DDE. The highest concentration of As and DDE was observed at the site strongly influenced by the watershed through the Po di Volano and minor canals. This influence of watershed and irrigation canals was reported also by [8]: peaks of pesticides were observed in areas of the lagoon coinciding with freshwater plumes, whereas water exchange with the adjacent sea favour herbicide dilution and export, keeping their concentrations at lower levels in other parts of the lagoon. The problem is the difficulty of separating the effect of contaminants from other factors, especially salinity, which is typically lower close to river mouths. Salinity itself can directly influence species’ richness: a decrease in richness with decreasing salinity has been reported as a common pattern from many transitional areas, even if the relationship is not necessarily strong and linear [51,52]. Moreover, salinity is also recognised as an important driving factor in metal toxicity, since free metal ion are more concentrated, with related bioavailability and toxicity, in conditions of low salinity [53]. Therefore, it is not possible to disentangle the effects of these two pollutants and salinity on the macrobenthic community, but it is likely that salinity and pollutants played a synergic role. Pollutants themselves can act synergistically when they co-occur. Organisms in the aquatic environment are commonly exposed to chemical mixtures, and the toxicity of a mixture is usually higher than each of the individual substances [54]. Toxicological investigations on the polychaete Laeonereis acuta showed that prolonged exposure to high levels of As resulted in a conspicuous reduction in the activity of an enzyme involved in the detoxification of several compounds, including organophosphorus pesticides [55]. The effects of contaminated freshwater released in wetlands was observed also in other lagoons, such as the Yellow River Delta [56], where heavy metal contamination of sediments and macrobenthic organisms, together with a dominance of pollution-resistant species were recorded in the affected area.
Interestingly, there was a lack of response of benthic organisms to other metals (namely Ni, Cr, Pb), present in even higher concentrations in relation to their EQS and their ERL/ERM thresholds [19]. The complexity of the response of the macrobenthic community to metal pollution reflects the fact that the measurement of contaminant levels alone does not directly indicate the bioavailability of pollutants to organisms [57,58]. The way in which high levels of pollutants in sediments affect the macrobenthic community, could vary in relation to abiotic processes, causing metal mobilizations, and to the physiological or genetic response of benthic organisms. Recent field and laboratory studies have documented the fact that environmental history can change the sensitivity of benthic organisms to stressful conditions [54,59]. The chronic exposure of animals to trace metals can considerably affect their physiological and biochemical responses, leading to increased tolerance to chronic contamination [60]. This adaptation ability was observed in different acquatic species, such as the bivalves R. philippinarum [61], and Macoma balthica [62], and the isopod Platynympha longicaudata [63].
Our results show that the methods currently used for monitoring transitional waters are insufficient to disentangle the effect of pollutants and other environmental variables, such as salinity or confinement, on benthic organisms and community structures. Integrated interdisciplinary approaches with new additional analyses, for instance, focusing on bioaccumulation and toxicity tests, are needed for a more precise identification of the risk posed by high concentrations of pollutants in transitional waters.

Author Contributions

Conceptualization, M.M.; methodology, A.S., C.R.F.; formal analysis, C.R.F., C.M.; investigation, C.M., A.A.S.; writing—original draft preparation, V.P.; writing—review and editing, M.M., V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Munari, C.; Mistri, M. Spatio-temporal pattern of community development in dredged material used for habitat enhancement: A study case in a brackish lagoon. Mar. Pollut. Bull. 2014, 89, 340–347. [Google Scholar] [CrossRef] [PubMed]
  2. Friedman, A.; Reverdin, G.; Khodri, M.; Gastineau, G. A new record of Atlantic sea surface salinity from 1896 to 2013 reveals the signatures of climate variability and long-term trends. Geophys. Res. Lett. 2017, 44, 1866–1876. [Google Scholar] [CrossRef] [Green Version]
  3. European Community. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 establishing a framework for community action in the field of water policy. Off. J. Eur. Communities 2000, 43, 1–72. [Google Scholar]
  4. Newton, A.; Icely, J.; Cristina, S.; Brito, A.; Cardoso, A.C.; Colijn, F. An overview of ecological status, vulnerability and future perspectives of European large shallow, semi-enclosed coastal systems, lagoons and transitional waters. Estuar. Coast. Shelf Sci. 2014, 140, 95–122. [Google Scholar] [CrossRef]
  5. Natali, C.; Fogli, R.; Bianchini, G.; Tassinari, R.; Tessari, U. Heavy Metals Backgrounds in Sediments From the Sacca di Goro (NE, Italy). EQA-Int. J. Environ. Qual. 2016, 20, 15–26. [Google Scholar]
  6. Barhoumi, B.; LeMenach, K.; Devier, M.-H.; Ameur, W.B.; Etcheber, H.; Budzinski, H. Polycyclic aromatic hydrocarbons (PAHs) in surface sediments from the Bizerte Lagoon, Tunisia: Levels, sources, and toxicological significance. Environ. Monit. Assess. 2014, 186, 2653–2669. [Google Scholar] [CrossRef] [PubMed]
  7. Borja, A.; Barbone, E.; Basset, A.; Borgersen, G.; Brkljacic, M.; Elliott, M.; Garmendia, J.M.; Marques, J.C.; Mazik, K.; Muxika, I.; et al. Response of single benthic metrics and multi-metric methods to anthropogenic pressure gradients, in five distinct European coastal and transitional ecosystems. Mar. Pollut. Bull. 2011, 62, 499–513. [Google Scholar] [CrossRef] [PubMed]
  8. Carafa, R.; Marinov, D.; Dueri, S.; Wollgast, J.; Ligthart, J.; Canuti, E.; Viaroli, P.; Zaldívar, J.M. A 3D hydrodynamic fate and transport model for herbicides in Sacca di Goro coastal lagoon (Northern Adriatic). Mar. Pollut. Bull. 2006, 52, 1231–1248. [Google Scholar] [CrossRef]
  9. 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]
  10. Rainbow, P.S. Trace metal bioaccumulation: Models, metabolic availability and toxicity. Environ Intl 2007, 33, 576–582. [Google Scholar] [CrossRef]
  11. Overmans, S.; Nordborg, M.; Rua, R.D.; Brinkman, D.L.; Negri, A.P.; Agustí, S. Phototoxic effects of PAH and UVA exposure on molecular responses and developmental success in coral larvae. Aquat. Toxicol. 2018, 198, 165–174. [Google Scholar] [CrossRef] [PubMed]
  12. Corbau, C.; Munari, C.; Mistri, M.; Lovo, S.; Simeoni, U. Application of the Principles of ICZM for Restoring the Goro Lagoon. Coast. Manage. 2016, 44, 350–365. [Google Scholar] [CrossRef]
  13. USEPA Method 6020B: Inductively Coupled Plasma-Mass Spectrometry; USEPA, United States Environmental Protection Agency: Washington, DC, USA, 2014.
  14. Anastassiades, M.; Lehotay, S.; Stajnbaher, D.; Schenck, F. Fast and easy multiresidue method employing acetonitrile extraction/partitioning and “dispersive solid-phase extraction” for the determination of pesticide residues in produce. J. AOAC Int. 2003, 86, 412–431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. USEPA. Method 8270D: Semivolatile organic compounds by Gas Chromatography/Mass Spectrometry (GC/MS); USEPA (United States Environmental Protection Agency): Washington, DC, USA, 2007. [Google Scholar]
  16. Sarkar, D. Lattice: Multivariate Data Visualization with R; Springer: New York, NY, USA, 2008. [Google Scholar]
  17. 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]
  18. Guerra, R.; Pasteris, A.; Lee, S.; Park, N.; Ok, G. Spatial patterns of metals, PCDDs/Fs, PCBs, PBDEs and chemical status of sediments from a coastal lagoon (Pialassa Baiona, NW Adriatic, Italy). Mar. Pollut. Bull. 2014, 89, 407–416. [Google Scholar] [CrossRef] [PubMed]
  19. 18 Long, E.R.; Macdonald, D.D.; Smith, S.L.; Calder, F.D. Incidence of adverse biological effects within ranges of chemical concentrations in marine and estuarine sediments. Environ. Manage. 1995, 19, 81–97. [Google Scholar] [CrossRef]
  20. Myles, H.; Wolfe, D.A. Nonparametric Statistical Methods; John Wiley & Sons: New York, NY, USA, 1973; pp. 139–146. [Google Scholar]
  21. Nemenyi, P. Distribution-free Multiple Comparisons. Biometrics 1963, 18, 263. [Google Scholar]
  22. Borja, A.; Franco, J.; Pérez, V. A marine biotic index to establish the ecological quality of soft-bottom benthos within European estuarine and coastal environments. Mar. Pollut. Bull. 2000, 40, 1100–1114. [Google Scholar] [CrossRef]
  23. Muxika, I.; Borja, A.; Bald, J. Using historical data, expert judgement and multivariate analysis in assessing reference conditions and benthic ecological status, according to the European Water Framework Directive. Mar. Pollut. Bull. 2007, 55, 16–29. [Google Scholar] [CrossRef]
  24. Borja, A.; Marín, S.; Muxika, I.; Pino, L.; Rodríguez, J. Is there a possibility of ranking benthic quality assessment indices to select the most responsive to different human pressures? Mar. Pollut. Bull. 2015, 97, 85–94. [Google Scholar] [CrossRef]
  25. Anderson, M.J.; Gorley, R.N.; Clarke, K.R. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods; PRIMER-E: Plymouth, UK, 2008. [Google Scholar]
  26. Legendre, P.; Anderson, M.J. Distance-based redundancy analysis: Testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 1999, 69, 1–24. [Google Scholar] [CrossRef]
  27. Clarke, K.R.; Gorley, R.N. PRIMER v6: User Manual/Tutorial; PRIMER-E: Plymouth, UK, 2006. [Google Scholar]
  28. RDevelopmentCoreTeam “R: A Language and Environment for Statistical Computing”; R Foundation for Statistical Computing: Vienna, Austria, 2008; Available online: https://www.R-project.org/ (accessed on 12 March 2020).
  29. 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. 2007, 96, 389–396. [Google Scholar] [CrossRef]
  30. Pitacco, V.; Mistri, M.; Ferrari, C.R.; Munari, C. Heavy metals, OCPs, PAHs, and PCDD/Fs contamination in surface sediments of a coastal lagoon (Valli di Comacchio, NW Adriatic, Italy): Long term trend (2002−2013) and effect on benthic community. Mar. Pollut. Bull. 2018, 135, 1221–1229. [Google Scholar] [CrossRef]
  31. Arienzo, M.; Masuccio, A.; Ferrara, L. Evaluation of sediment contamination by heavy metals, organochlorinated pesticides, and polycyclic aromatic hydrocarbons in the Berre coastal lagoon (southeast France). Arch. Environ. Contam. Toxicol. 2013, 65, 396–406. [Google Scholar] [CrossRef] [PubMed]
  32. Casatta, N.; Mascolo, G.; Roscioli, C.; Viganò, L. Tracing endocrine disrupting chemicals in a coastal lagoon (Sacca di Goro, Italy): Sediment contamination and bioaccumulation in Manila clams. Sci. Total Environ. 2015, 511, 214–222. [Google Scholar] [CrossRef] [PubMed]
  33. Quensen, J.F.; Mueller, S.A.; Jain, M.K.; Tiedje, J.M. Reductive dechlorination of DDE to DDMU in marine sediment microcosms. Science 1998, 280, 722–724. [Google Scholar] [CrossRef] [Green Version]
  34. Hoke, R.; Cotter, A.; Goldenstein, T.; Kosian, P.; Ankley, G.; Phipps, G. Evaluation of equilibrium partitioning theory for predicting acute toxicity of field-collected sediments contaminated with DDT, DDE and DDD to the amphipod Hyalella azteca. Environ. Toxicol. Chem. 1994, 13, 157–166. [Google Scholar] [CrossRef]
  35. Swartz, R.C.; Cole, F.A.; Lamberson, J.O.; Ferraro, S.P.; Schults, D.W.; Deben, W.A.; Lee Ii, H.; Ozretich, R.J. Sediment toxicity, contamination and amphipod abundance at a DDT-and dieldrin-contaminated site in San Francisco Bay. Environ. Toxicol. Chem. 1994, 13, 949–962. [Google Scholar] [CrossRef]
  36. Bellucci, L.; Frignani, M.; Raccanelli, S.; Carraro, C. Polychlorinated dibenzo-p-dioxins and dibenzofurans in surficial sediments of the Venice Lagoon (Italy). Mar. Pollut. Bull. 2000, 40, 65–76. [Google Scholar] [CrossRef]
  37. Gómez-Lavín, S.; Gorri, D.; Irabien, Á. (Assessment of PCDD/Fs and PCBs in sediments from the Spanish northern Atlantic coast. Water Air Soil Pollut. 2011, 221, 287–299. [Google Scholar] [CrossRef]
  38. Danis, B.; Debacker, V.; Miranda, C.T.; Dubois, P. Levels and effects of PCDD/Fs and co-PCBs in sediments, mussels, and sea stars of the intertidal zone in the southern North Sea and the English Channel. Ecotoxicol. Environ. Saf. 2006, 65, 188–200. [Google Scholar] [CrossRef] [PubMed]
  39. Dueri, S.; Marinov, D.; Fiandrino, A.; Tronczyński, J.; Zaldívar, J.M. Implementation of a 3D coupled hydrodynamic and contaminant fate model for PCDD/Fs in Thau Lagoon (France): The importance of atmospheric sources of contamination. Int. J. Env. Res. Public Health 2010, 7, 1467–1485. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Specchiulli, A.; Renzi, M.; Scirocco, T.; Cilenti, L.; Florio, M.; Breber, P. Comparative study based on sediment characteristics and macrobenthic communities in two Italian lagoons. Environ. Monit. Assess. 2010, 160, 237–256. [Google Scholar] [CrossRef] [PubMed]
  41. Fabbri, D.; Baravelli, V.; Giannotti, K.; Donnini, F.; Fabbri, E. Bioaccumulation of cyclopenta [cd] pyrene and benzo[ghi]fluoranthene by mussels transplanted in a coastal lagoon. Chemosphere 2006, 64, 1083–1092. [Google Scholar] [CrossRef]
  42. Acquavita, A.; Falomo, J.; Predonzani, S.; Tamberlich, F.; Bettoso, N.; Mattassi, G. The PAH level, distribution and composition in surface sediments from a Mediterranean Lagoon: The Marano and Grado Lagoon (Northern Adriatic Sea, Italy). Mar. Pollut. Bull. 2014, 81, 234–241. [Google Scholar] [CrossRef]
  43. Giacalone, A.; Gianguzza, A.; Mannino, M.R.; Orecchio, S.; Piazzese, D. Polycyclic aromatic hydrocarbons in sediments of marine coastal lagoons in Messina, Italy: Extraction and GC/MS analysis, distribution and sources. Polycyc Aromat. Compd. 2004, 24, 135–149. [Google Scholar] [CrossRef]
  44. Culotta, L.; De Stefano, C.; Gianguzza, A.; Mannino, M.R.; Orecchio, S. The PAH composition of surface sediments from Stagnone coastal lagoon, Marsala (Italy). Mar. Chem. 2006, 99, 117–127. [Google Scholar] [CrossRef]
  45. Countway, R.E.; Dickhut, R.M.; Canuel, E.A. Polycyclic aromatic hydrocarbon (PAH) distributions and associations with organic matter in surface waters of the York River, VA Estuary. Org. Geochem. 2003, 34, 209–224. [Google Scholar] [CrossRef]
  46. Pearson, T.; Rosenberg, R. Macrobenthic succession in relation to organic enrichment and pollution of the marine environment. Oceanogr. Mar. Biol. Ann. Rev. 1978, 16, 229–311. [Google Scholar]
  47. Borja, A.; Muxika, I. Guidelines for the use of AMBI (AZTI’s Marine Biotic Index) in the assessment of the benthic ecological quality. Mar. Pollut. Bull. 2005, 50, 787–789. [Google Scholar] [CrossRef]
  48. 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] [PubMed]
  49. Prato, S.; La Valle, P.; De Luca, E.; Lattanzi, L.; Migliore, G.; Morgana, J.G.; Munari, C.; Nicoletti, L.; Izzo, G.; Mistri, M. The "one-out, all-out" principle entails the risk of imposing unnecessary restoration costs: A study case in two Mediterranean coastal lakes. Mar. Pollut. Bull. 2014, 80, 30–40. [Google Scholar] [CrossRef] [PubMed]
  50. Pitacco, V.; Reizopoulou, S.; Sfriso, A.; Sfriso, A.A.; Mistri, M.; Munari, C. The difficulty of disentangling natural from anthropic forcing factors makes the evaluation of ecological quality problematic: A case study from Adriatic lagoons. Mar. Environ. Res. 2019, 150, 104756. [Google Scholar] [CrossRef] [PubMed]
  51. Canu, D.M.; Solidoro, C.; Umgiesser, G.; Cucco, A.; Ferrarin, C. Assessing confinement in coastal lagoons. Mar. Pollut. Bull. 2012, 64, 2391–2398. [Google Scholar] [CrossRef] [PubMed]
  52. Jorcin, A. Temporal and spatial variability in the macrozoobenthic community along a salinity gradient in the Castillos Lagoon (Uruguay). Archiv. Hydrobiol. 1999, 146, 369–384. [Google Scholar] [CrossRef]
  53. Telesh, I.V.; Khlebovich, V.V. Principal processes within the estuarine salinity gradient: A review. Mar. Pollut. Bull. 2010, 61, 149–155. [Google Scholar] [CrossRef]
  54. Blackmore, G.; Wang, W.X. Uptake and efflux of Cd and Zn by the green mussel Perna viridis after metal preexposure. Environ. Sci. Technol. 2002, 36, 989–995. [Google Scholar] [CrossRef]
  55. Relyea, R.A. A cocktail of contaminants: How mixtures of pesticides at low concentrations affect aquatic communities. Oecologia 2008, 159, 363–376. [Google Scholar] [CrossRef]
  56. Li, M.; Yang, W.; Sun, T.; Jin, Y. Potential ecological risk of heavy metal contamination in sediments and macrobenthos in coastal wetlands induced by freshwater releases: A case study in the Yellow River Delta, China. Mar. Pollut. Bull. 2016, 103, 227–239. [Google Scholar] [CrossRef] [Green Version]
  57. Ventura-Lima, J.; Sandrini, J.Z.; Cravo, M.F.; Piedras, F.R.; Moraes, T.B.; Fattorini, D. Toxicological responses in Laeonereis acuta (Annelida, Polychaeta) after arsenic exposure. Environ. Int. 2007, 33, 559–564. [Google Scholar] [CrossRef]
  58. Depledge, M.H.; Galloway, T.S. Healthy animals, healthy ecosystems. Front Ecol Environ 2005, 3, 251–258. [Google Scholar] [CrossRef]
  59. Yang, W.; Li, X.; Pei, J.; Sun, T.; Shao, D.; Bai, J. Bioavailability of trace metals in sediments of a recovering freshwater coastal wetland in China’s Yellow River Delta, and risk assessment for the macrobenthic community. Chemosphere 2017, 189, 661–671. [Google Scholar] [CrossRef] [PubMed]
  60. Shi, D.; Blackmore, G.; Wang, W.X. Effects of aqueous and dietary preexposure and resulting body burden on silver biokinetics in the green mussel Perna viridis. Environ. Sci. Technol. 2003, 37, 936–943. [Google Scholar] [CrossRef] [PubMed]
  61. Paul-Pont, I.; De Montaudouin, X.; Gonzalez, P.; Soudant, P.; Baudrimont, M. How life history contributes to stress response in the Manila clam Ruditapes philippinarum. Environ. Sci. Pollut. Res. 2010, 17, 987–998. [Google Scholar] [CrossRef]
  62. Boisson, F.; Hartl, M.G.; Fowler, S.W.; Amiard-Triquet, C. Influence of chronic exposure to silver and mercury in the field on the bioaccumulation potential of the bivalve Macoma balthica. Mar. Environ. Res. 1998, 45, 325–340. [Google Scholar] [CrossRef]
  63. Ross, K.; Cooper, N.; Bidwell, J.R.; Elder, J. Genetic diversity and metal tolerance of two marine species: A comparison between populations from contaminated and reference sites. Mar. Pollut Bull. 2002, 44, 671–679. [Google Scholar] [CrossRef]
Figure 1. Map of the studied area with sampling sites.
Figure 1. Map of the studied area with sampling sites.
Water 12 01074 g001
Figure 2. Sediment parameters of samples collected from 2006 to 2010 (n = 4). The average annual values (±SD) of sediment classification (sand > 50µm, silt from 50 to 2 µm, clay < 2µm). The values of organic matter are expressed in mg/kg, values of sand, silt, and clay are expressed in %.
Figure 2. Sediment parameters of samples collected from 2006 to 2010 (n = 4). The average annual values (±SD) of sediment classification (sand > 50µm, silt from 50 to 2 µm, clay < 2µm). The values of organic matter are expressed in mg/kg, values of sand, silt, and clay are expressed in %.
Water 12 01074 g002
Figure 3. Annual average of polycyclic aromatic hydrocarbons (PAHs) (Σ16PAHs in µg/kg) for 2004 to 2010 (n = 4) and group profiles of LPAHs (two or three-fused rings) and of HPAHs (four to six-fused rings) in the sediments from Sacca di Goro. The EQS limit (Good/Moderate) for the sum of PAHs according to the WFD is indicated.
Figure 3. Annual average of polycyclic aromatic hydrocarbons (PAHs) (Σ16PAHs in µg/kg) for 2004 to 2010 (n = 4) and group profiles of LPAHs (two or three-fused rings) and of HPAHs (four to six-fused rings) in the sediments from Sacca di Goro. The EQS limit (Good/Moderate) for the sum of PAHs according to the WFD is indicated.
Water 12 01074 g003
Figure 4. Annual mean (±SD) of percentage of Ecological Groups (sensitive species = EGI, indifferent = EGII, tolerant = EGIII, second order opportunistic = EGIV, first order opportunistic = EGV), structural indices (richness = S, Shannon index = H) and ecological indices (AMBI and M-AMBI) for macrobenthic community in Sacca di Goro from 2004 to 2010 (n = 3). Grey line = Good/Moderate threshold.
Figure 4. Annual mean (±SD) of percentage of Ecological Groups (sensitive species = EGI, indifferent = EGII, tolerant = EGIII, second order opportunistic = EGIV, first order opportunistic = EGV), structural indices (richness = S, Shannon index = H) and ecological indices (AMBI and M-AMBI) for macrobenthic community in Sacca di Goro from 2004 to 2010 (n = 3). Grey line = Good/Moderate threshold.
Water 12 01074 g004
Figure 5. Values of species richness (S) against the concentration of As (A) and DDE (B), with threshold values for S (high/good, horizontal dashed line), As, and DDE (good/moderate, vertical dotted line), defined by National Law (n = 3).
Figure 5. Values of species richness (S) against the concentration of As (A) and DDE (B), with threshold values for S (high/good, horizontal dashed line), As, and DDE (good/moderate, vertical dotted line), defined by National Law (n = 3).
Water 12 01074 g005
Table 1. Yearly average (±SD), maximum and minimum yearly concentrations of trace elements (TEs) from 2004 to 2010 in sediments of the Sacca di Goro (n = 4), threshold limits, and the number of years from 2004 to 2010 in which TEs’ concentrations exceeded those limits expressed in percentage. EQS = limits for ecological quality status (Good/Moderate) according to the Water Framework Directive (WFD); Ref = background levels of Ni, Pb, Cr and Hg for Po River Delta [18]; ERM = Effects Range-Low values [19]; ERL = Effects Range-Low values [19].
Table 1. Yearly average (±SD), maximum and minimum yearly concentrations of trace elements (TEs) from 2004 to 2010 in sediments of the Sacca di Goro (n = 4), threshold limits, and the number of years from 2004 to 2010 in which TEs’ concentrations exceeded those limits expressed in percentage. EQS = limits for ecological quality status (Good/Moderate) according to the Water Framework Directive (WFD); Ref = background levels of Ni, Pb, Cr and Hg for Po River Delta [18]; ERM = Effects Range-Low values [19]; ERL = Effects Range-Low values [19].
Ni (mg/kg)Cr (mg/kg)Hg (mg/kg)Pb (mg/kg)Cd (mg/kg)As (mg/kg)
200495.5 ± 7.880.6 ± 8.90.0 ± 0.029.6 ± 7.20.4 ± 0.110.0 ± 3.3
200583.9 ± 2.477.4 ± 9.50.1 ± 0.228.8 ± 10.10.3 ± 0.29.3 ± 1.2
2006110.3 ± 14.2149.9 ± 17.70.2 ± 0.136.0 ± 14.60.4 ± 0.29.6 ± 3.5
2007113.9 ± 7.2157.5 ± 8.30.1 ± 0.143.5 ± 13.90.5 ± 0.29.7 ± 3.6
200893.4 ± 7.5143.6 ± 13.10.3 ± 0.133.6 ± 8.40.4 ± 0.16.7 ± 1.5
2009100.1 ± 16.0129.6 ± 28.40.2 ± 0.132.0 ± 10.20.4 ± 0.18.3 ± 2.8
2010106.7 ± 10.3137.2 ± 14.60.1 ± 0.128.6 ± 13.90.4 ± 0.19.2 ± 3.5
Max113.9 ± 7.2157.5 ± 8.30.0 ± 0.043.5 ± 13.90.5 ± 0.210.0 ± 3.3
Min83.9 ± 2.477.4 ± 9.60.3 ± 0.128.6 ± 13.90.3 ± 0.26.7 ± 1.5
EQS (mg/kg)30500.3300.312
>EQS100%100%14%57%100%0%
Ref (mg/kg)971580.1224--
>Ref71%0%43%100%--
ERL (mg/kg)20.9810.1546.71.28.2
>ERL10.710.43000.86
ERM (mg/kg)51.63700.712189.670
>ERM100%0%0%0%0%0%
Table 2. Mean (±SD) annual organochlorine pesticides (OCPs) concentrations (µg/kg) in sediments of Sacca di Goro (n = 4). The EQS according to the WFD and D.M. 260/2010.
Table 2. Mean (±SD) annual organochlorine pesticides (OCPs) concentrations (µg/kg) in sediments of Sacca di Goro (n = 4). The EQS according to the WFD and D.M. 260/2010.
HexachlorobenzeneDDEDDDDDTSimazineTerbuthylazineAlaclorΣOCPs
2004<0.10.06 ± 0.130.03 ± 0.05<0.1<0.1<0.1<0.10.09 ± 0.18
20050.15 ± 0.301.04 ± 1.290.40 ± 0.80<0.1<0.1<0.1<0.11.59 ± 1.88
2006<0.11.41 ± 1.21<0.1<0.1<0.1<0.1<0.11.41 ± 1.21
20070.09 ± 0.071.91 ± 1.040.28 ± 0.490.08 ± 0.230.15 ± 0.421.25 ± 1.670.38 ± 0.604.14 ± 2.44
2008<0.11.09 ± 0.441.35 ± 2.20<0.1<0.1<0.1<0.12.44 ± 2.19
2009<0.11.51 ± 0.750.35 ± 0.60<0.1<0.1<0.1<0.11.86 ± 1.29
2010<0.12.47 ± 1.340.69 ± 0.88<0.1<0.1<0.1<0.13.16 ± 2.07
EQS0.41.80.81
ERL 2.21.58
ERM 2746.1
Table 3. Mean (±SD) annual concentrations (ng/kg) PCCDs and dibenzofurans (PCDFs) (n = 4), and TEQs (toxic equivalents) calculated for total PCCD/Fs, based on Toxic Equivalency Factors (TEFs) according to the WHO (World Health Organization).
Table 3. Mean (±SD) annual concentrations (ng/kg) PCCDs and dibenzofurans (PCDFs) (n = 4), and TEQs (toxic equivalents) calculated for total PCCD/Fs, based on Toxic Equivalency Factors (TEFs) according to the WHO (World Health Organization).
H7CDDO8CDDH7CDFO8CDFΣPCDD/FsΣTEQ
2007Mean ± SD19.58 ± 22.60123.91 ± 91.5113.30 ± 7.5321.67 ± 16.71178.460.47
TEQ0.200.120.130.02
2008Mean ± SD18.83 ± 16.8989.14 ± 82.676.66 ± 4.7822.14 ± 22.92136.760.37
TEQ0.190.090.070.02
2009Mean ± SD12.74 ± 9.4188.40 ± 58.458.59 ± 6.4318.41 ± 13.68128.150.32
TEQ0.130.090.090.02
2010Mean ± SD17.23 ± 11.1094.18 ± 72.888.58 ± 4.7921.15 ± 10.62141.130.37
TEQ0.170.090.090.02
Table 4. Collinearity (rs) among abiotic and biotic variables (n = 3). Only variables showing significant correlations were displayed. NS = p > 0.05.
Table 4. Collinearity (rs) among abiotic and biotic variables (n = 3). Only variables showing significant correlations were displayed. NS = p > 0.05.
NiCrAsPbPCDFsOMDDEClayM-AMBI AMBIH
Ni
Cr0.56
As0.69NS
Hg−0.54NSNS
PbNS0.62NS
PCDFsNSNSNS0.7
PCDDsNSNSNS0.90.89
OM−0.61−0.58NSNSNS
DDE 0.69NS0.89NSNS−0.54
SandNSNSNSNSNS0.55NS
Silt0.55NSNSNSNSNSNS
AMBINSNSNSNSNSNSNS−0.69−0.59
HNSNSNSNSNSNSNSNS0.89−0.64
SNSNS−0.67NSNSNS−0.62NS0.71NSNS

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/w12041074

AMA Style

Pitacco V, Mistri M, Ferrari CR, Sfriso A, Sfriso AA, Munari C. Multiannual Trend of Micro-Pollutants in Sediments and Benthic Community Response in a Mediterranean Lagoon (Sacca di Goro, Italy). Water. 2020; 12(4):1074. https://doi.org/10.3390/w12041074

Chicago/Turabian Style

Pitacco, Valentina, Michele Mistri, Carla Rita Ferrari, Adriano Sfriso, Andrea Augusto Sfriso, and Cristina Munari. 2020. "Multiannual Trend of Micro-Pollutants in Sediments and Benthic Community Response in a Mediterranean Lagoon (Sacca di Goro, Italy)" Water 12, no. 4: 1074. https://doi.org/10.3390/w12041074

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