Mercury cycling in contaminated coastal environments: modeling the benthic-pelagic coupling and microbial resistance in the Venice Lagoon.

Anthropogenic activities have been releasing mercury for centuries, and despite global efforts to control emissions, concentrations in environmental media remain high. Coastal sediments can be a long-term repository for mercury, but also a secondary source, and competing processes in marine ecosystems can lead to the conversion of mercury into the toxic and bioaccumulative species methylmercury, which threatens ecosystem and human health. We investigate the fate and transport of three mercury species in a coastal lagoon affected by historical pollution using a novel high-resolution finite element model that integrates mercury biogeochemistry, sediment dynamics and hydrodynamics. The model resolves mercury dynamics in the seawater and the seabed taking into account partitioning, transport driven by water and sediment, and photochemical and microbial transformations. We simulated three years (early 2000s, 2019, and 2020) to assess the spatio-temporal distribution of mercury species concentrations and performed a sensitivity analysis to account for uncertainties. The modeled mercury species concentrations show high temporal and spatial variability, with water concentrations in some areas of the lagoon exceeding those of the open Mediterranean Sea by two orders of magnitude, consistent with available observations from the early 2000s. The results support conclusions about the importance of different processes in shaping the environmental gradients of mercury species. Due to the past accumulation of mercury in the lagoon sediments, inorganic mercury in the water is closely related to the resuspension of contaminated sediments, which is significantly reduced by the presence of benthic vegetation. The gradients of methylmercury depend on the combination of several factors, of which sediment resuspension and mercury methylation are the most relevant. The results add insights into mercury dynamics at coastal sites characterized by a combination of past pollution (i.e. sediment enrichment) and erosive processes, and suggest possible nature-based mitigation strategies such as the preservation of the integrity of benthic vegetation and morphology.

Why coastal lagoons are so productive? Physical bases of fishing productivity in coastal lagoons.

Coastal lagoons are among the most productive marine ecosystems in the world. Annual primary production varies from 50 to > 500 g C m-2 year-1, being of the same order of magnitude as that of the upwelling areas. Many lagoons lie within the range of eutrophic (300-500 g C m-2 year-1) or hypereutrophic (> 500 g C m-2 year-1) conditions. The high productivity of coastal lagoons makes them subject of exploitation by many marine fishes and invertebrates, that use them as nursery areas and feeding grounds during their early life cycle phases, and most lagoons support important fisheries or maintain aquaculture exploitations. The high levels of their biological production can be explained by some of their common features as shallowness and the strong influence of terrestrial systems. Shallowness favors that the photic zone extends to the lagoon bottom and that wind can promote the resuspension of nutrients and organisms. The interaction with land also introduces significant amounts of nutrients. However, trophic variables can explain < 43 % of the fishing yields, and further than the trophic status of the lagoons, several works showed that the biological productivity of coastal lagoons can be explained by their geomorphological features such as the positive influence of shoreline development and the negative influence of depth. Using the Mar Menor lagoon as a case study, we propose that although nutrient inputs and light can be limiting factors for photosynthetic based productivity, increasing fishing yield up to a certain limit, the productivity of lagoons is mainly promoted by more general forces associated to physical and chemical gradients.

Hot moments and hotspots of cyanobacteria hyperblooms in the Curonian Lagoon (SE Baltic Sea) revealed via remote sensing-based retrospective analysis.

A temporally and spatially detailed historical (1985-2018) analysis of cyanobacteria blooms was performed in the Curonian Lagoon (Lithuania, Russia), the largest coastal lagoon in the Baltic Sea. Satellite data allowed the mapping of cyanobacteria surface accumulations, so-called "scums", and of chlorophyll-a concentration. The 34-year time series shows a tendency towards later occurrence (October-November) of the cyanobacteria scum presence, whereas the period of its onset (June-July) remains relatively constant. The periods when scums are present, "hot moments", have been consistently increasing in duration since 2008. The differences in the starting, ending and annual duration of cyanobacteria blooms have been significantly altered by hydro-meteorological conditions (river discharge, water temperature, and wind conditions) and their year-round patterns. The most important environmental factors that determined the temporal changes of the scum presence and area were the standing stock of cyanobacteria and the ambient wind conditions. The "hotspots", the areas where the blooms most likely occur, were distributed in the south-southwestern and central parts of the lagoon. The least affected areas were the northern part, which is connected to the coastal waters of the Baltic Sea, and the Nemunas River delta region. The longstanding, well-established spatial patterns of cyanobacteria blooms were linked to hydrodynamic features, namely water renewal time and current patterns, and to potential nutrient sources that included muddy sediments and the locations of colonies of piscivorous birds. Our findings confirmed that the annual and seasonal variations of cyanobacteria blooms and their regulation are a complex issue due to interactions between multiple factors over spatially and temporally broad scales. Despite great progress in the prevention and control of eutrophication and cyanobacteria blooms, the lagoon is still considered to be in a poor ecological status. This work provides a new and missing understanding on the spatial and temporal extent of cyanobacteria blooms and the factors that govern them. Such an understanding can help in planning management strategies, forecasting the magnitude and severity of blooms under changing nutrient loads and potential climate scenarios.

Modelling as decision support for the localisation of submarine urban wastewater outfall: Venice lagoon (Italy) as a case study.

Microbiological impact is critical in coastal areas where tourism is particularly important for both the local and regional economy. Submarine outfalls are commonly used to enhance the dispersion of treated sewage thus avoiding pollution along the coast. The Venice lagoon (North Italy) has a very sensitive ecosystem, due to the morphological and natural characteristics of the basin and the co-existence of human activities. To preserve the lagoon, the discharge from the treatment plant for urban wastewater collected from the Venezia-Mestre agglomeration, neighbouring areas and local industries (total of 400,000 population equivalent-PE) has been moved from the lagoon to the open Adriatic Sea since November 2013 by means of an approximately 20-km pipeline. Microbiological pollution inside the lagoon can affect shellfish breeding areas instead, along the coast it affects the quality of bathing waters. In this study, and for the first time, a 3D hydrodynamic SHYFEM model (shallow water finite element model) with high spatial resolution coupled with a microbiological module has been applied to the lagoon and to the Adriatic Sea, to evaluate the effectiveness of the location of the submarine outfall. Microbiological data have been produced by the control Authority according to official analytic methods and by the plant operator. The module of survival of free Escherichia coli follows a variable rate in dependence of UV radiation, temperature and salinity in the water. Two scenarios were modelled: final discharge into the lagoon before November 2013 and after into the open sea. In the latter case, two situations have been considered, one with "Bora" and the other with "Scirocco" winds. Our results indicate that the model correctly simulates microbiological decay and dispersion. The transferral of the final discharge point far from the shoreline improves pollution dispersion, thus preserving the lagoon without evidence of impacts on the bathing waters in all meteorological conditions.

Natural resources and climate change: A study of the potential impact on Manila clam in the Venice lagoon.

A crucial aspect in climate change is to understand how an ecosystem will adapt under different environmental conditions and how it will influence the ecological resources and the connected human activities. In this study, a numerical model reproduces the growth dynamics, dispersion and settlement of clam's larvae in the Venice lagoon. On the basis of the last IPCC scenarios for the years 2050 and 2100, the model simulates the changes in larval settlement, showing how the geographical distribution and, consequently, the nursery area changes over time. Our results indicate that climate change will modify, not only the timing of the settlements (from spring-summer to winter autumn) and the spatial distribution of nursery areas (from central to southern lagoon), but also the absolute quantity of settled larvae in the lagoon. This can strongly affect aquaculture in terms of availability of seed and farming practice. Given that these changes are due to the variations in temperature and circulation, similar processes are likely to happen in other transitional environments all over the world affecting the global aquaculture resources. In this regard, the tool we developed could support local policymakers in the knowledge-based planning and sustainable management of clam aquaculture in vulnerable environments.

GIS-based multi-criteria site selection for zebra mussel cultivation: Addressing end-of-pipe remediation of a eutrophic coastal lagoon ecosystem.

Farming of shellfish and seaweeds is a tested tool for mitigating eutrophication consequences in coastal environments, however as many other marine economic activities it should be a subject of marine spatial planning for designating suitable sites. The present study proposes site selection framework for provisional zebra mussel farming in a eutrophic lagoon ecosystem, aimed primarily at remediation purposes. GIS-based multi-criteria approach was applied, combining data from empirical maps, numerical models and remote sensing to estimate suitability parameters. Site selection and prioritisation of suitable areas considered 15 environmental and socio-economic criteria, which contributed to 4 optimisation models (settlement, growth and survival of mussels, environmental and socio-economic) and 3 predefined scenarios representing provisional goals of mussel cultivation: spat production, biomass production and bioremediation. The relative importance of each criterion was assessed utilizing the Analytical Hierarchy Process. Site suitability index was calculated and the final result of the site selection analysis was summarized for 3 scenarios and overall suitability map. Four suitability classes (unsuitable, least, moderately and most suitable) were applied, and 3 most suitable zones for provisional zebra mussel cultivation with 12 candidate sites were selected accordingly. The integrated approach presented in this study can be adjusted for designating zebra mussel farming sites in other estuarine lagoon ecosystems, or cultivation of other mussel species for bioremediation purposes. The analytical framework and the workflow designed in this study are also adoptable for addressing other aquaculture-related spatial planning issues.

High resolution multibeam and hydrodynamic datasets of tidal channels and inlets of the Venice Lagoon.

Tidal channels are crucial for the functioning of wetlands, though their morphological properties, which are relevant for seafloor habitats and flow, have been understudied so far. Here, we release a dataset composed of Digital Terrain Models (DTMs) extracted from a total of 2,500 linear kilometres of high-resolution multibeam echosounder (MBES) data collected in 2013 covering the entire network of tidal channels and inlets of the Venice Lagoon, Italy. The dataset comprises also the backscatter (BS) data, which reflect the acoustic properties of the seafloor, and the tidal current fields simulated by means of a high-resolution three-dimensional unstructured hydrodynamic model. The DTMs and the current fields help define how morphological and benthic properties of tidal channels are affected by the action of currents. These data are of potential broad interest not only to geomorphologists, oceanographers and ecologists studying the morphology, hydrodynamics, sediment transport and benthic habitats of tidal environments, but also to coastal engineers and stakeholders for cost-effective monitoring and sustainable management of this peculiar shallow coastal system.

The Trapping Index: How to integrate the Eulerian and the Lagrangian approach for the computation of the transport time scales of semi-enclosed basins.

In this work, we investigated if the Eulerian and the Lagrangian approaches for the computation of the Transport Time Scales (TTS) of semi-enclosed water bodies can be used univocally to define the spatial variability of basin flushing features. The Eulerian and Lagrangian TTS were computed for both simplified test cases and a realistic domain: the Venice Lagoon. The results confirmed the two approaches cannot be adopted univocally and that the spatial variability of the water renewal capacity can be investigated only through the computation of both the TTS. A specific analysis, based on the computation of a so-called Trapping Index, was then suggested to integrate the information provided by the two different approaches. The obtained results proved the Trapping Index to be useful to avoid any misleading interpretation due to the evaluation of the basin renewal features just from an Eulerian only or from a Lagrangian only perspective.

Assessing confinement in coastal lagoons.

Measures of transport scale in aquatic systems can contribute to the formulation of definitions of indicators of the system's ecological properties. This paper addresses confinement, a specific transport scale proposed by biological scientists as a parameter that can capture and synthesize the principal properties that determine the spatial structure of biological communities in transitional environments. Currently, there is no direct experimental measure of confinement. In this study, a methodology based on the accumulation rate within a lagoon of a passive tracer of marine origin is proposed, the influences of different factors in the calculation of confinement are analyzed, and general recommendations are derived. In particular, we analyze the spatial and the temporal variability of confinement and its sensitivity to the seasonal variability of climatic forcing, the inputs from rivers and the parameterization of the tidal exchanges. The Lagoon of Venice is used as a case study.

Investigation of residence time and groundwater flux in Venice Lagoon: comparing radium isotope and hydrodynamic models.

The four naturally-occurring isotopes of radium were coupled with a previously evaluated hydrodynamic model to determine the apparent age of surface waters and to quantify submarine groundwater discharge (SGD) into the Venice Lagoon, Italy. Mean apparent age of water in the Venice Lagoon was calculated using the ratio of 224Ra to 228Ra determined from 30 monitoring stations and a mean pore water end member. Average apparent age was calculated to be 6.0 d using Ra ratios. This calculated age was very similar to average residence time calculated for the same period using a hydrodynamic model (5.8 d). A mass balance of Ra was accomplished by quantifying each of the sources and sinks of Ra in the lagoon, with the unknown variable being attributed to SGD. Total SGD were calculated to be 4.1 +/- 1.5, 3.8 +/- 0.7, 3.0 +/- 1.3, and 3.5 +/- 1.0 x 10(10) L d(-1) for (223,224,226, 228)Ra, respectively, which are an order of magnitude larger than total mean fluvial discharge into the Venice Lagoon (3.1 x 10(9) L d(-1)). The SGD as a source of nutrients in the Venice Lagoon is also discussed and, though significant to the nutrient budget, is likely to be less important as the dominant control on SGD is recirculated seawater rather than freshwater.

Submarine wastewater discharges: dispersion modelling in the Northern Adriatic Sea.

BACKGROUND, AIM AND SCOPE: Opposite interests must coexist in coastal areas: the presence of significant cities and urban centres, of touristic and recreational areas, and of extensive shellfish farming. To avoid local pollution caused by treated wastewaters along the Northern Adriatic coast (Friuli Venezia-Giulia and Veneto regions), marine outfall systems have been constructed. In this study, the application of a numerical dispersion model is used to support the traditional monitoring methods in order to link information concerning the hydrodynamic circulation and the microbiological features, to evaluate possible health risks associated with recreational and coastal shellfish farming activities. The study is a preliminary analysis of the environmental impact of wastewater treatment plants (WWTPs) with submarine discharge outfalls. It also could be useful for the water profile definition according to the Directive 2006/7/EC on the quality of bathing water and for the integrated areal analysis (Ostoich et al. 2006), to define the area of influence of each submarine discharge point. MATERIALS AND METHODS: Historical data on discharges of the considered WWTPs were recovered and evaluated. Data on discharges' control for Veneto region (WWTPs of Lido and Cavallino) were produced by the WWTPs' manager Veritas Laboratory service, while data for the WWTPs of Friuli Venezia-Giulia region were produced by the regional environmental protection agency in the institutional control activity following official methods. The hydrodynamic model used in this work is the three-dimensional version of the finite element model SHYFEM, developed at ISMAR-CNR (Marine Science Institute of the Italian National Research Council) in Venice (Umgiesser et al. J Mar Syst 51:123-145, 2008). RESULTS AND DISCUSSION: Numerical simulations have been carried out with the 3D version of the finite element model SHYFEM for 3 months during autumn 2007 to evaluate the bacterial pollution dispersion along the coasts of Veneto and Friuli Venezia-Giulia regions, prescribing meteo-marine forcings and concentration values at the points corresponding to the positions of the submarine outfalls. Model results show that during autumn 2007 the discharges of the submarine outfalls of the Venice province seem to have no impact on the surface water quality, while there are some visible effects in the Gulf of Trieste. This reflects the behaviour of the experimental data collected by ARPAV and ARPA FVG and monitoring campaigns both on water and shellfish quality. Further results have been elaborated to identify the area of influence of each discharge point; scenarios were developed with imposed concentrations. The results seem to highlight that the two discharges of the Veneto region are not noticeable, while the discharges of the Gulf of Trieste (in particular the Servola and Barcola ones) are perceptible. CONCLUSIONS: This study represents a new step towards the study of the microbiological pollution dispersion and impact due to the discharges of the submarine outfalls of the Veneto and Friuli Venezia-Giulia regions (nine considered discharge points). With the 3D version of the finite element model SHYFEM, the information obtained from the hydrodynamic circulation has been linked to the classical methods of analysis, to assess possible risks connected to the microbiological parameter Escherichia coli. RECOMMENDATIONS AND PERSPECTIVES: In future studies the time scale for microbiological parameters' decay could be linked to various environmental parameters such as light climate, temperature, and salinity. Interesting information would come from the study of new scenarios with different configurations of the discharge of the pipelines and/or the treatment plants and in particular from the improvements of the 3D version of the SHYFEM model, to take the stratification process into account which occurs during spring-summer, since the Northern Adriatic Sea is a very complex ecosystem, both as physical and ecological processes.