Next Article in Journal
Nature-Based Coastal Protection by Large Woody Debris as Compared to Seawalls: A Physical Model Study of Beach Morphology and Wave Reflection
Next Article in Special Issue
Antibiotic Resistance Genes and Potentially Pathogenic Bacteria in the Central Adriatic Sea: Are They Connected to Urban Wastewater Inputs?
Previous Article in Journal
Hybrid Scour Depth Prediction Equations for Reliable Design of Bridge Piers
Previous Article in Special Issue
Modelling the Quality of Bathing Waters in the Adriatic Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial and Temporal Vertical Distribution of Chlorophyll in Relation to Submarine Wastewater Effluent Discharges

1
Center of Excellence for Science and Technology-Integration of Mediterranean Region, University of Split, 21000 Split, Croatia
2
Doctoral Study of Biophysics, Faculty of Science, University of Split, 21000 Split, Croatia
3
Faculty of Civil Engineering, Architecture and Geodesy, University of Split, M. Hrvatske 15, 21000 Split, Croatia
4
Department of Biology, Faculty of Science, University of Split, R. Boškovića 33, 21000 Split, Croatia
*
Author to whom correspondence should be addressed.
Submission received: 15 June 2021 / Revised: 18 July 2021 / Accepted: 21 July 2021 / Published: 23 July 2021

Abstract

:
The vertical distribution of chlorophyll in coastal waters is influenced by a combination of the hydrodynamic environment and different biotic and abiotic processes. The spatial and temporal occurrences of chlorophyll profiles provide a good representation of the changes in the marine environment. The majority of studies in the Adriatic Sea have so far been conducted in areas unaffected by anthropogenic pressure. Our study site is located near two marine outfalls, which are part of the public sewage system. This study aims to characterize the chlorophyll vertical distribution and describe its variability based on the stratification conditions and the presence of a wastewater effluent plume. Based on these conditions, we identified three characteristic scenarios/types of chlorophyll profiles. The first one occurs when the vertical mixing of the water column creates the upwelling of chlorophyll and nutrients to the upper part of the water column. The second and third scenarios occur during stratified conditions and differ by the extent of the effluent plume intrusion. Using modern fluorescence techniques, we identified and described three different vertical chlorophyll profiles, characterizing them according to their physical and biological parameters and processes. For cases with a visible effluent intrusion, we confirmed the importance of the pycnocline formation in keeping the effluent below and maintaining the higher water quality status at the top of the water column.

1. Introduction

Phytoplankton has a major influence on marine ecosystems because of its role in oxygen production and overall significance in supporting and maintaining marine life. Due to the challenges in the direct measuring of phytoplankton, chlorophyll a is mainly used as a proxy for phytoplankton biomass [1,2,3,4,5]. Traditional techniques for chlorophyll analysis are based on bottle sampling, which is limited in terms of obtaining a high spatial and temporal resolution. This was corrected by modern in situ fluorescence techniques and remote sensing observations, which allow for measurements of chlorophyll in a high resolution. The research on chlorophyll temporal and spatial variability in the Mediterranean Sea is mostly restricted to the surface, as most of the collected data come from satellite observations and ex-situ laboratory measurements [5]. The oligotrophic status of the Mediterranean basin [6], with a decreasing trend from west to east [7,8], and seasonal and interannual variability [4,9,10] were observed based on surface chlorophyll analysis.
The vertical distribution of chlorophyll is under the influence of many biotic and abiotic processes, whose interplay affects the depth of the upper mixed layer and controls the vertical chlorophyll distribution [4,11]. So far, information about vertical chlorophyll profiles has been derived from different kinds of measurements [12,13,14,15,16,17,18,19,20,21,22,23].
A similar chlorophyll seasonality was noticed in the Adriatic Sea by remote sensing observations [24,25,26] and laboratory measurements [27,28]. The Adriatic Sea is a semi-enclosed basin of the Mediterranean Sea, where the biogeochemical parameters are influenced mainly by the circulation of the water masses. The chlorophyll trends in the Adriatic Sea have been investigated by remote sensing [29,30,31] and laboratory-based measurements [27,32] to find deep chlorophyll maxima [28,33] and to investigate the anthropogenic impact on the phytoplankton community structure and abundancy [34,35,36].
Coastal waters are highly sensitive ecosystems and may become even more vulnerable when exposed to the negative impact of human activities. Wastewater effluents contain high concentrations of pollutants, which can cause a large amount of deterioration to marine environments [37]. Bacteria are the most abundant and diverse group of organisms living in wastewater [38,39], and bathing water status is mainly determined using fecal indicator bacteria (FIB) as indicators of water quality and related human health issues [40,41]. Besides fecal pollution, eutrophication is another important aspect of coastal water quality. Since primary producers are the first components of an ecosystem to react to an increase in nutrient concentration, phytoplankton biomass is considered a good indicator of the trophic state of coastal water [42]. Changes in phytoplankton biomass abundance occur following nutrient enrichment, which can happen naturally (through riverine influence or the upwelling of nutrients from the bottom) or under anthropogenic influence (urban and industrial wastewater discharges and surface runoff). Anthropogenically introduced nutrients are the major problem of coastal ecosystems [43,44], and their mixing with seawater contributes to the eutrophication of coastal zone waters [45], thus increasing the biomass of macroalgae and seston algae [46].
This paper aims to characterize chlorophyll vertical profiles and describe their spatial and temporal variability based on the stratification status of a water column and the presence of a marine wastewater effluent plume. Due to the nature of fluorescence measurements, which change with the taxonomic species or environmental conditions [47] and can be affected by dissolved materials [48] or by different light irradiation statuses [49,50,51], the present study is focused on the relative distribution of vertical chlorophyll fluorescence in a water column by comparing different segments of the water column for each considered profile. To our knowledge, this is the first classification of chlorophyll vertical profiles in the presence of a wastewater effluent plume. The spatiotemporal variability of the different vertical profiles concerning stratification was also investigated.

2. Materials and Methods

2.1. Study Site

The study area is the eastern central Adriatic Sea, near the city of Split in Croatia (Figure 1). Two main marine wastewater outfalls are located within the study area, about 5 km away from each other, as a part of the Split–Solin public sewage system.
The submarine outfall at Location 1 discharges treated wastewater from the Katalinića brig wastewater treatment plant (WWTP) via a Ø 800 mm pipeline, ending with two diffuser sections 1.3 km away from the shore at a 43 m depth. The Katalinića brig WWTP is designed for the mechanical treatment of both municipal and rainfall-runoff wastewater, with a capacity of 122.000 population equivalents (PE) and an average flow rate equal to 35.000 m3/day approximately.
The marine outfall at Location 2 (Stobreč) discharges the wastewater coming from the Stupe WWTP, where the wastewater is mechanically treated, along with grit and grease removal. The Stupe WWTP facility is designed for municipal wastewater only, with a capacity of 138.000 PE and an average flow rate equal to 30.000 m3/day. The length of the submarine outfall Ø 900 mm pipeline is 2.76 km, with a 200 m long diffuser section ending at a 37 m depth.
Additionally, the monitoring of different physical and chemical parameters is performed several times a year, at stations FP-O14 and FP-O14b, as shown in Figure 1, as a part of the national monitoring program.
To understand the background of the wastewater being discharged into the seawater, we provide data (Table 1) about the seasonal flowrate and available water effluent parameters for two WWTPs.

2.2. Sampling Methodology

This study was performed during the implementation of the Interreg Italy–Croatia project AdSWiM (Managed use of treated urban wastewater for the quality of the Adriatic Sea). The measurements were performed with the CTD fast profiling probe (The Sea & Sun Technology) and C3 fluorometer during sampling sessions conducted with a boat (Nautica 600) every month from February to September 2020 (with the exception of March and August due to the COVID restrictions).
CTD fast profiling probe, with an additional sensor for dissolved oxygen (model: 48M), was used for temperature, salinity, density anomaly, and dissolved oxygen measurements. The temperature is measured by a classical resistance thermometer using a platinum resistor, while the salinity is calculated from the seawater conductivity using a 7-poll cell sensor. The density anomalies were calculated automatically by the CTD/DO probe internal software using the international thermodynamic equation of seawater-2010 [52]. The dissolved oxygen sensor works on the principle of a Clark electrode. The depth is calculated from the pressure, measured by a strain gauge, and is converted to depth as 1dbar = 1 m. The turbidity, chlorophyll, and colored dissolved organic matter (CDOM) were measured using a C3 Submersible fluorometer (Turner Design). The chlorophyll measurements were performed with a blue mercury lamp, with a peak emission at 460 nm, and the fluorescence was collected at 680 nm. The CDOM was measured using a UV LED (CWL: 365 nm), with a peak emission at 350 nm and fluorescence collection at 450 nm. The turbidity was measured with the IR lamp, with a peak emission at 850 nm, and the collection of scattered light at 90°. Since the relationship between phytoplankton biomass and chlorophyll fluorescence is highly variable and depends on the physiological state and community composition [53,54,55], in this study, we focused on the general distribution of chlorophyll profiles in the water column concerning the stratification and effluent discharge.
The fecal pollution was estimated by the level of FIB. The FIB concentration was assessed by the membrane filtration method and the enumeration of fecal coliforms (E. coli) [56,57] and intestinal enterococci (ENT) [58], as indicated in the Croatian legislation (Regulation on the Quality of Marine Bathing Waters; OG 73/08). The number of viable heterotrophic bacteria was determined using the spread plate method. An automated Protos3 colony counting system (Synbiosis, UK) was used to determine the bacteriological concentration. Sampling was conducted from the boat using a Niskin sampler. The samples were collected in sterile bottles, protected from light in cool boxes, and transported to the microbiology laboratory of the Faculty of Science in Split to be processed within 4 h.
The chlorophyll concentration data at locations FP-O14 and FP-O14b were collected during the regular monitoring performed by the water management agency (Croatia Waters). The chlorophyll a concentrations were determined by the laboratory fluorometric method, with filtration through Whatman GF/C glass filters, at 0, 5, 10, 20, and 50 m depths at the FP-O14 location and at 0, 5, 10, and 43 m depths at FP-O14b.

2.3. Data Analysis

The data collection was performed by the vertical profiling of the water column down to 40 m and 35 m depths at location 1 and location 2, respectively. A moving average filter was applied to raw data sampled at a maximum spatial resolution of 3 cm to achieve spatial uniformity. The data were smoothed by a moving average filter, with a 20 cm window.
The vertical distributions of the different layers were examined by calculating the vertical gradients of the density anomalies and derived buoyancy frequency. The buoyancy frequency or Brunt-Väisälä frequency (N) represents the strength of the density stratification and, consequently, the stability of the water column. It is defined as:
N = g ρ 0 d ρ dz
where g is the gravitational acceleration, ρ 0 is the averaged density of the whole water profile, and d ρ /dz is the water density gradient at each sampling step. To identify events when stratification occurred, data points are selected where N2 exceeds 0.001 s−2 as possible pycnocline candidates [59]. The data points are then aggregated if the distance between them was less than 0.8 m, forming a possible pycnocline layer. Layers thinner than 0.8 m are neglected to exclude the impact of anomalous extremes. Aggregated layers are classified as a pycnocline and their profiles as stratified, which potentially aids in the creation of chlorophyll layers. To support the analysis of forming chlorophyll layers under the influence of marine wastewater discharge, the vertical profiles of dissolved oxygen, turbidity, CDOM, temperature, and salinity are also examined.

3. Results

The main results from six sampling cruises are presented in Figure 2 and Figure 3 for location 1 and location 2, respectively. Each Figure contains six vertical profiles, corresponding to the sampling month with the set of parameters measured. Analysis of the buoyancy frequency showed that stratification occurred in May, June, July, and September for both locations. In February and April, the winter and early spring, lower temperatures induced vertical mixing in the water column. As a result of the strong sea surface warming, the shallowest stratification occurred in July, while the deepest occurred in September. At both locations 1 and 2, the temperature and salinity followed their seasonal pattern, with low gradients in winter and steeper gradients during stratified periods.

3.1. Location 1

3.1.1. February 2020

During the February sampling cruise, the winter mixing of the water column, induced by wind and low temperatures, allowed the upwelling of the bottom seawater rich with nutrients, causing the rise of chlorophyll concentrations from the bottom to the upper layers of the water column (Figure 2, February). The February vertical profiles shown in Figure 2 also depict a slight increase in the turbidity and CDOM around a 25 m depth, probably due to the phytoplankton activity, which can contribute to the organic matter by excreting photosynthetic products or by lysing [60]. The homogeneous vertical dissolved oxygen profile indicates a high oxygen enrichment throughout the water column, mostly as a consequence of the higher oxygen production and water column aeration during the winter weather conditions.

3.1.2. April 2020

The April vertical profiles shown in Figure 2 show a higher concentration of chlorophyll at the bottom, along with higher values of CDOM and turbidity. The observed dip in the dissolved oxygen levels near the bottom indicates the presence of a wastewater plume. To support this hypothesis, the observed high level of fecal pollution is presented in Figure 4. The high levels of E. coli and ENT support the dissolved oxygen dip, as bacteria consume oxygen from seawater during the process of decomposition of organic matter from the wastewater plume. The homogeneous vertical profile of dissolved oxygen in the upper part of the water column indicates that the oxygen production exceeds the consumption in that part of the water column.

3.1.3. May 2020

The stratification of the water column started in May. The chlorophyll vertical profile, shown in Figure 2, has no gradients in the upper layer, with a steady increase below the forming stratification point. The increase in the chlorophyll concentration is followed by an increase in CDOM and turbidity. Since the dissolved oxygen profile was homogeneous, we presume that the phytoplankton production of oxygen dominated over the consumption, without any observed effluent effect.

3.1.4. June 2020

During the June sampling cruise (Figure 2), we observe a broad chlorophyll peak with three plateaus, starting below a depth of 15 m. The chlorophyll peaks coincide with the CDOM and turbidity peaks, indicating an effluent plume layer. Consequently, the phytoplankton biomass increases, further contributing to the organic matter concentration and turbidity rise. This increase is matched by a decrease in salinity and dissolved oxygen. The oxygen decrease is a consequence of the bacterial oxidation of organic matter, which is flushed into the seawater by the marine outfall discharge. In the June sampling cruise, the stratification is much steeper, compared to May, and it occurs at a depth of around 18 m.

3.1.5. July 2020

The July vertical profiles shown in Figure 2 show that the chlorophyll vertical profile is not comprised of several plateaus, like the one in June, but has a steep slope on the deeper side and two shoulders on the shallow side. Furthermore, it is present much deeper than the stratification layer, which occurs at a depth of 7 m. The increase in chlorophyll is matched by an increase in CDOM and turbidity values and a decrease in dissolved oxygen and salinity. The observed peaks of different parameters during the July cruise indicate the presence of an effluent plume below the pycnocline. The dissolved oxygen increases down to the 24 m depth and then begins to decline near the chlorophyll peak.

3.1.6. September 2020

The September vertical profiles shown in Figure 2 show the pycnocline presence down to the 22–23 m depth, which was the deepest pycnocline observed during our cruises. Below the pycnocline, at a 25 m depth, an increase in chlorophyll, CDOM, and turbidity is observed, together with a small decrease in dissolved oxygen at the same depth. Near the bottom of the water column, the oxidation of organic matter reduced the dissolved oxygen concentration, indicating that more oxygen was being used than produced. In this case, the deep pycnocline probably submerged the effluent plume, keeping it near the bottom of the water column.

3.2. Location 2

3.2.1. February 2020

During the February sampling cruise, a similar vertical chlorophyll profile occurred at location 2 as that at location 1. The February chlorophyll profile shown in Figure 3 shows a rise in the chlorophyll concentrations from the bottom to the top of the water column. The homogeneous vertical turbidity, CDOM, salinity, and dissolved oxygen profiles indicate the water column mixing. Furthermore, the elevated levels of FIB (Figure 5) at the surface of the water column indicate the effluent rising to the surface of the water column due to the mixing process under the unstratified water column conditions.

3.2.2. April 2020

The April vertical profiles shown in Figure 3 indicate that the mixing of the water column creates the upwelling of the bottom water, causing a similar vertical chlorophyll profile to the one observed in February. The observed peaks in CDOM and turbidity, within the mixing process, are due to the available nutrients for primary production, resulting in the phytoplankton activity. The homogeneous dissolved oxygen vertical profile indicates a high oxygen enrichment throughout the water column, reflecting surplus phytoplankton production and an aerated water column due to the lack of stratification.

3.2.3. May 2020

Stratification of the water column at location 2 was first observed in May 2020. The observed chlorophyll, CDOM, and turbidity profiles in May 2020 (Figure 3) reveals no significant concentration and/or gradients. Such situations are mostly a consequence of the complex physical, biological, and chemical interactions, together with hydrodynamic processes, particularly due to the fact that the wastewater effluent is only periodically discharged into the sea. However, the May chlorophyll profile, shown in Figure 3, shows an increase in values near the bottom of the water column. The increase in the chlorophyll concentration is matched by a small increase in CDOM and turbidity near the bottom, without a decline in the dissolved oxygen vertical profile. Since the dissolved oxygen profile was homogeneous, we presume that the phytoplankton production of oxygen exceeded the consumption by bacteria.

3.2.4. June 2020

The June vertical profiles presented in Figure 3 show the pycnocline at a 14 m depth, after which an increase in chlorophyll, CDOM, and turbidity is observed. Higher chlorophyll concentrations below a 15 m depth, in conjunction with the CDOM and turbidity increase, indicate the presence of an effluent plume layer accumulated below the pycnocline. The presence of organic matter, together with a higher phytoplankton biomass, leads to an increase in turbidity. Furthermore, the input of organic material contributes to the reduction of available oxygen through bacteria oxidation processes, which is supported by the shape of the vertical dissolved oxygen profile. The reduced salinity at 15 to 20 m depths supports the conclusion that these changes are a direct consequence of the wastewater plume. In comparison to the other months, in June 2020 the highest chlorophyll, CDOM, and turbidity values are measured at the bottom of the water column. The concentration of 220 CFU/100 mL of E. coli indicates the presence of an effluent, in agreement with other observed parameters.

3.2.5. July 2020

During the July sampling cruise (Figure 3, July), pycnocline occurred at a very shallow depth of 4 m. Increased values of chlorophyll, CDOM, and turbidity are observed from 15 m to 35 m depths. The observed peaks of chlorophyll, turbidity, and CDOM, with a homogeneously dissolved oxygen profile, down to the 25 m depth, indicate the presence of phytoplankton biomass, producing dissolved oxygen, which prevailed over the consumption. A sharp CDOM peak at the 15 m depth was due to the phytoplankton activity, as no dissolved oxygen decrease was observed at that depth. The decrease in the dissolved oxygen below the 25 m depth indicates that the plume effluent had a more significant effect in this layer. The oxygen consumption during the organic matter decomposition at this depth prevailed over the production, causing an increase of CDOM and turbidity in the bottom layer. The observed concentration of 103 CFU/100 mL for E. coli indicates an intrusion of fecal material by the marine outfall. Thus, the July profile showed the characteristics of two competing processes. One is a higher oxygen production from the phytoplankton activity and the other is a higher oxygen consumption from the oxidation of organic material entering the marine system by the wastewater discharge.

3.2.6. September 2020

September vertical profiles in Figure 3 show a homogeneous chlorophyll profile within the upper layer up to the pycnocline at a 15 m depth. Below the stratification, the chlorophyll concentration is increasing towards the bottom of the water column. The organic matter peaks were observed below the pycnocline, while the turbidity increases only near the bottom. As the vertical dissolved oxygen profile shows little variability, we presume that the organic matter reflects phytoplankton activity, indicating that the production dominated over the consumption.

4. Discussion

The majority of the vertical chlorophyll profile studies conducted in coastal region and estuaries are based on laboratory measurements or satellite techniques, lacking direct in situ measurements, especially the influence of marine outfalls. Therefore, this study aimed to identify the chlorophyll vertical profiles under the influence of the water column stratification status and the presence of an effluent plume. Based on the measurement results, we classify the chlorophyll profiles into three types (Table 2). The profiles observed in February and April 2020, with the observed rise of the chlorophyll concentrations from the bottom to the top of the water column, are classified as the Upper Chlorophyll Profile (UCP). The other two profile types are divided between the sampling sessions when stratification of the water column was observed, and they differ by the presence of an effluent plume. In May at both locations and in September at location 2, the vertical chlorophyll profile was characterized by a higher concentration near the bottom, without a decrease in dissolved oxygen. These profiles are classified as the Bottom Chlorophyll Profile (BCP). The third profile type is characterized by peaks of chlorophyll coinciding with CDOM and turbidity peaks and a matching decrease in salinity and dissolved oxygen, which is a consequence of an effluent plume intrusion. These profiles are classified as the Effluent Chlorophyll Profile (ECP). A summary of this classification is given in Table 2.
In unstratified conditions, observed in February and April 2020, at both locations, the upwelling of the bottom, nutrient-rich waters induced a rise in the chlorophyll concentrations from the bottom to the top of the water column, which we classified as UCP. This rise is supported by the CDOM and turbidity vertical profiles, which show only small variations along the vertical column. The bottom water upwelling also resulted in a higher FIB level at the surface, indicating a rise in fecal pollution due to marine outfalls. In unstratified conditions, the dissolved oxygen vertical profiles were homogeneous. This is due to vertical mixing, which allows air–sea gas exchange with the bottom waters and the influence of external forces, such as wind [61]. An exception occurred in April 2020 at location 1 due to an effluent plume intrusion, which decreased the dissolved oxygen values. These conditions induce the characteristics of the chlorophyll profile of both the ECP and UCP (Figure 2). The oxygenation of the water column down to a 30 m depth was due to the primary production of phytoplankton, but it was also due to the unstratified water column, which caused vertical mixing and aeration of the entire water column. At the bottom of the water column, the dissolved oxygen concentration and chlorophyll are inversely correlated, suggesting that the consumption of oxygen dominates over its production. Similar observations have been made in many estuaries with nutrient inputs from rivers [62,63,64,65,66,67]. Marine outfalls discharge nutrients into the sea, resulting in similar biological processes occurring as those when rivers discharge nutrients into the estuary.
From May to September, stratification occurred at both locations, and the observed profiles are classified as BCP (May 2020, September 2020) and ECP (June 2020, July 2020, September 2020) (Table 2). BCPs are characterized by a higher chlorophyll concentration near the bottom of the water column. This increase is likely due to the phytoplankton activity and resuspension of the bottom sediment, where dead plankton and other sources of organic matter fertilizer, unconsumed food, and fecal matter settle. In May 2020, despite the water column stratification, the vertical distribution of dissolved oxygen was homogeneous, similar to a phenomenon observed in some estuaries, where, under stratified conditions and higher temperatures, hypoxia occurred [68,69,70]. In our case, hypoxia was never observed, even in summer, but a higher dissolved oxygen variability was present.
To confirm the occurrence of UCP and BCP types in the absence of an effluent plume, we used data from measurements conducted by the Croatian water management agency in the period from 2012–2019, sampled at two locations (FP-O14 and FP-O14b). Both stations are used to represent coastal seawater unaffected by marine outfalls. At both of these locations, during months when no stratification is present, the vertical chlorophyll profile shows upwelling, indicating the existence of UCP (Figure 6 and Figure 7). Under stratified conditions, narrower ranges of the chlorophyll concentrations are observed in the upper part, with the highest values measured at the bottom, which are classified as BCP (Figure 6 and Figure 7). These data show the existence of UCP and BCP profiles in coastal, unpolluted waters, which follow a seasonal cycle of sea warming and the formation of stratified conditions, and differentiate the naturally occurring chlorophyll profiles from the ones induced by anthropogenic influence.
The vertical profiles of CDOM, turbidity, dissolved oxygen, and salinity further highlight the difference between ECP and BCP. The nutrient enrichment by marine outfalls induces phytoplankton biomass growth and increases the chlorophyll concentration. The highest chlorophyll concentration observed was at location 1 in July 2020. Such a surge in the chlorophyll concentration could be due to the accumulation of living algae, supported by raw sewage, containing higher amounts of nutrients, especially nitrogen and phosphorus. Besides nutrients, the marine outfalls also discharge organic matter into the coastal waters, resulting in increased organic matter concentrations and turbidity. The increased turbidity levels are the result of the input of organic matter from the wastewater discharge, creating a good condition for the presence of fecal bacteria. The effluent intrusion creates conditions where the oxygen consumption, during the bacterial oxidation of organic matter and decomposition of dead algae and plants, exceeds the production of dissolved oxygen by photosynthesis. The observed decrease in the dissolved oxygen profiles furthermore separates the BCP and ECP types. The gradient of dissolved oxygen with depth (higher concentration above the pycnocline) indicates that stratification may play a role in lowering the oxygen concentration by reducing the gas exchange between the air and sea surface layer [61]. While stratification can influence the dissolved oxygen profile, the homogeneously dissolved oxygen profiles observed in May suggest that the highest variability of dissolved oxygen concentration is due to the influence of a wastewater effluent plume and not stratification. The variations in salinity also supported our hypothesis about the nature of ECP. Similar changes in salinity were previously observed for other submarine outfall sites and provide further evidence of a freshwater input through the outfall system [71,72,73,74].
Due to the fact that both WWTPs are discharging a municipal wastewater effluent, the seasonal variations are due to the variations in the residential population of the Split municipal area, which increases during the summer months. Consequently, such changes can influence the wastewater volume. Thus, the summer discharges are approximately 12% higher than those during the winter period, as presented in Table 1. This is also shown in Figure 2 and Figure 3 in terms of an increased concentration and/or gradients of measured parameters. While the discharged wastewater and the seawater start mixing quickly and reaching equilibrium with the surrounding seawater, the salinity profiles may serve as additional information, aiding in distinguishing the chlorophyll profiles. During the observed stratification, the effluent plume was held below the pycnocline, which is generally ecologically favorable, as it maintains a higher quality of the surface waters.
In addition, the circulation in the wider area surrounding the two studied outfalls has a typical barocline current pattern, with relatively larger surface layer velocities due to the regular wind intensities. Furthermore, during the unstratified conditions, the bottom upwelling is intensified, creating the potential for effluent plume surfacing. In such cases, it is expected that surface currents will further contribute to the mixing and spreading of chlorophyl.

5. Conclusions

To study the wastewater impact on the chlorophyll vertical distribution, measurements were performed in the eastern central Adriatic Sea at two locations of marine outfalls, near the city of Split in Croatia. We describe the three vertical chlorophyll profile types, based on the stratification conditions and the effluent plume presence. The water column mixing induced the upwelling of the bottom water, generating UCP and enabling the rise of the effluent towards the surface. In the stratified conditions, categorized as the ECP type, the pycnocline prevents the plume from rising, keeping it in the lower part of the water column, while the BCP type of the profile represents a case when the water column is stratified, but no effluent plume is present. Our categorization of the different vertical chlorophyll profiles is supported by data collected by the Croatian water agency at the locations unaffected by the marine outfalls. All three observed chlorophyll profile types are the consequence of the interplay between the physical, biological, and chemical processes. To distinguish these processes and confidently identify the presence of the effluent plume we use the FIB measurements, along with dissolved oxygen gradients, to examine the conditions when the oxygen consumption exceeds the production, indicating organic matter oxidation. This, together with the salinity, CDOM, and turbidity measurements, allows us to separate the wastewater influence from the natural processes. We confirm the importance of the stratification in submerging the effluent below the pycnocline, which is highlighted in February 2020, when column mixing allowed the FIB to rise to the surface of the water column. To our knowledge, this is a first study involving in situ vertical measurements coupled with the submarine outfall, analyzing its influence on the vertical chlorophyll distribution.

Author Contributions

Conceptualization, M.K.; Methodology, M.K., M.V., T.K., A.M. and R.A.; Investigations: M.K., M.V., T.K. and R.A.; Writing—Original Draft Preparation, M.K.; Writing—Review and Editing, M.V., T.K., A.M. and R.A.; Visualization, M.K. and T.K.; Project Administration, R.A.; Funding Acquisition, R.A. All authors have read and agreed to the published version of the manuscript.

Funding

The contributions of M.K., M.V. and R.A. were supported by the project, STIM-REI, Contract Number: KK.01.1.1.01.0003, funded by the European Union through the European Regional Development Fund—the Operational Programme Competitiveness and Cohesion 2014–2020 (KK.01.1.1.01). T.K. and R.A. acknowledge the support from the project, CAAT (Coastal Auto-purification Assessment Technology), also funded by European Union from European Structural and Investment Funds 2014–2020, Contract Number: KK.01.1.1.04.0064. The contribution of R.A. was supported by the project, COMON (COastal zone MONitoring using multi-scaling methods), funded by European Union from European Structural and Investment Funds 2014–2020, Contract Number: KK.01.1.07.0033. This research is partially supported through project KK.01.1.1.02.0027, a project co-financed by the Croatian Government and the European Union through the European Regional Development Fund—the Competitiveness and Cohesion Operational Programme.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available from the corresponding author, upon a request.

Acknowledgments

The cruise monitoring was completed through the implementation of the Interreg CBC Italy–Croatia project, “Managed use of treated urban wastewater for the quality of the Adriatic Sea” (AdSWiM), by the Faculty of Civil Engineering, Architecture and Geodesy, University of Split. Microbiological data from April 2020 are obtained by the Institute of Public Health Zadar. Seawater data gathered at stations FP-O14 and FP-O14b was provided by the Croatian water management administration (Croatia Waters). We thank Nika Ugrin for her valuable help during the measurement campaigns.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lorenzen, C.J. A method for the continuous measurement of in vivo chlorophyll concentration. Deep Sea Res. Oceanogr. Abstr. 1966, 13, 223–227. [Google Scholar] [CrossRef]
  2. Cullen, J.J. The Deep Chlorophyll Maximum: Comparing Vertical Profiles of Chlorophyll a. Can. J. Fish. Aquat. Sci. 1982, 39, 791–803. [Google Scholar] [CrossRef]
  3. Palter, J.; Coto, S.L.; Ballestero, D. The distribution of nutrients, dissolved oxygen and chlorophyll a in the upper Gulf of Nicoya, Costa Rica, a tropical estuary. Rev. Biol. Trop. 2007, 55, 427–436. [Google Scholar] [CrossRef] [Green Version]
  4. Lavigne, H.; D’Ortenzio, F.; Migon, C.; Claustre, H.; Testor, P.; D’Alcalà, M.R.; Lavezza, R.; Houpert, L.; Prieur, L. Enhancing the comprehension of mixed layer depth control on the Mediterranean phytoplankton phenology. J. Geophys. Res. Ocean. 2013, 118, 3416–3430. [Google Scholar] [CrossRef] [Green Version]
  5. Lavigne, H.; D’Ortenzio, F.; Ribera D’Alcalà, M.; Claustre, H.; Sauzède, R.; Gacic, M. On the vertical distribution of the chlorophyll a concentration in the Mediterranean Sea: A basin-scale and seasonal approach. Biogeosciences 2015, 12, 5021–5039. [Google Scholar] [CrossRef] [Green Version]
  6. Dugdale, R.C.; Wilkerson, F.P. Nutrient sources and primary production in the Eastern Mediterranean. Oceanologica Acta. In Proceedings of the Oceanographie Pelagique Mediterraneene, Villefranche-sur-Mer, France, 16–20 September 1988; pp. 179–184. [Google Scholar]
  7. Bosc, E.; Bricaud, A.; Antoine, D. Seasonal and interannual variability in algal biomass and primary production in the Mediterranean Sea, as derived from 4 years of SeaWiFS observations. Glob. Biogeochem. Cycles 2004, 18. [Google Scholar] [CrossRef] [Green Version]
  8. Barale, V.; Jaquet, J.M.; Ndiaye, M. Algal blooming patterns and anomalies in the Mediterranean Sea as derived from the SeaWiFS data set (1998–2003). Remote Sens. Environ. 2008, 112, 3300–3313. [Google Scholar] [CrossRef]
  9. D’Ortenzio, F.; Ribera d’Alcalà, M. On the trophic regimes of the Mediterranean Sea: A satellite analysis. Biogeosciences 2009, 6, 139–148. [Google Scholar] [CrossRef] [Green Version]
  10. Volpe, G.; Nardelli, B.B.; Cipollini, P.; Santoleri, R.; Robinson, I.S. Seasonal to interannual phytoplankton response to physical processes in the Mediterranean Sea from satellite observations. Remote Sens. Environ. 2012, 117, 223–235. [Google Scholar] [CrossRef]
  11. Falkowski, P.; Woodhead, A. Primary Productivity and Biogeochemical Cycles in the Sea; Springer: Heidelberg, Germany, 1992; Volume 43, ISBN 0306441926. [Google Scholar]
  12. Macías, D.; Stips, A.; Garcia-Gorriz, E. The relevance of deep chlorophyll maximum in the open Mediterranean Sea evaluated through 3D hydrodynamic-biogeochemical coupled simulations. Ecol. Modell. 2014, 281, 26–37. [Google Scholar] [CrossRef]
  13. Crise, A.; Allen, J.I.; Baretta, J.; Crispi, G.; Mosetti, R.; Solidoro, C. The Mediterranean pelagic ecosystem response to physical forcing. Prog. Oceanogr. 1999, 44, 219–243. [Google Scholar] [CrossRef]
  14. Moutin, T.; Raimbault, P. Primary production, carbon export and nutrients availability in western and eastern Mediterranean Sea in early summer 1996 (MINOS cruise). J. Mar. Syst. 2002, 33–34, 273–288. [Google Scholar] [CrossRef]
  15. Crombet, Y.; Leblanc, K.; Quéuiner, B.; Moutin, T.; Rimmelin, P.; Ras, J.; Claustre, H.; Leblond, N.; Oriol, L.; Pujo-Pay, M. Deep silicon maxima in the stratified oligotrophic Mediterranean Sea. Biogeosciences 2011, 8, 459–475. [Google Scholar] [CrossRef] [Green Version]
  16. Dolan, J.R.; Claustre, H.; Carlotti, F.; Plounevez, S.; Moutin, T. Microzooplankton diversity: Relationships of tintinnid ciliates with resources, competitors and predators from the Atlantic Coast of Morocco to the Eastern Mediterranean. Deep. Res. Part I Oceanogr. Res. Pap. 2002, 49, 1217–1232. [Google Scholar] [CrossRef]
  17. Christaki, U.; Giannakourou, A.; Van Wambeke, F.; Grégori, G. Nanoflagellate predation on auto- and heterotrophic picoplankton in the oligotrophic Mediterranean Sea. J. Plankton Res. 2001, 23, 1297–1310. [Google Scholar] [CrossRef] [Green Version]
  18. Estrada, M.; Marrase, C.; Latasa, M.; Berdalet, E.; Delgado, M.; Riera, T. Variability of deep chlorophyll maximum characteristics in the northwestern Mediterranean. Mar. Ecol. Prog. Ser. 1993, 92, 289–300. [Google Scholar] [CrossRef]
  19. Casotti, R.; Landolfi, A.; Brunet, C.; D’Ortenzio, F.; Mangoni, O.; Ribera d’Alcalà, M.; Denis, M. Composition and dynamics of the phytoplankton of the Ionian Sea (eastern Mediterranean). J. Geophys. Res. Ocean. 2003, 108. [Google Scholar] [CrossRef]
  20. Marty, J.C.; Chiavérini, J.; Pizay, M.D.; Avril, B. Seasonal and interannual dynamics of nutrients and phytoplankton pigments in the western Mediterranean Sea at the DYFAMED time-series station (1991–1999). Deep. Res. Part II Top. Stud. Oceanogr. 2002, 49, 1965–1985. [Google Scholar] [CrossRef]
  21. Psarra, S.; Tselepides, A.; Ignatiades, L. Primary productivity in the oligotrophic Cretan Sea (NE Mediterranean): Seasonal and interannual variability. Prog. Oceanogr. 2000, 46, 187–204. [Google Scholar] [CrossRef]
  22. Krom, M.D.; Brenner, S.; Kress, N.; Neori, A.; Gordon, L.I. Nutrient dynamics and new production in a warm-core eddy from the Eastern Mediterranean Sea. Deep Sea Res. Part A Oceanogr. Res. Pap. 1992, 39, 467–480. [Google Scholar] [CrossRef]
  23. Siokou-Frangou, I.; Christaki, U.; Mazzocchi, M.G.; Montresor, M.; Ribera D’Alcala, M.; Vaque, D.; Zingone, A. Plankton in the open mediterranean Sea: A review. Biogeosciences 2010, 7, 1543–1586. [Google Scholar] [CrossRef] [Green Version]
  24. Barale, V.; Rizzoli, P.M.; Hendershott, M.C. Remotely sensing the surface dynamics of the Adriatic Sea. Deep Sea Res. Part A Oceanogr. Res. Pap. 1984, 31, 1433–1459. [Google Scholar] [CrossRef]
  25. Barale, V.; McClain, C.R.; Malanotte-Rizzoli, P. Space and time variability of the surface color field in the northern Adriatic Sea. J. Geophys. Res. 1986, 91, 12957. [Google Scholar] [CrossRef]
  26. Böhm, E.; Banzon, V.; D’Acunzo, E.; D’Ortenzio, F.; Santoleri, R. Adriatic Sea surface temperature and ocean colour variability during the MFSPP. Ann. Geophys. 2003, 21, 137–149. [Google Scholar] [CrossRef] [Green Version]
  27. Krivokapić, S.; Stanković, Ž.; Vuksanović, N. Seasonal variations of phytoplankton biomass and environmental conditions in the inner Boka Kotorska Bay (Eastern Adriatic sea). Acta Bot. Croat. 2009, 68, 45–55. [Google Scholar]
  28. Ninčević, Ž.; Marasović, I.; Kušpilić, G. Deep chlorophyll-a maximum at one station in the middle Adriatic Sea. J. Mar. Biol. Assoc. UK 2002, 82, 9–19. [Google Scholar] [CrossRef]
  29. Mauri, E.; Poulain, P.M.; Južnič-Zonta, Ž. MODIS chlorophyll variability in the northern Adriatic Sea and relationship with forcing parameters. J. Geophys. Res. Ocean. 2007, 112. [Google Scholar] [CrossRef]
  30. Kišević, M.; Morović, M.; Smailbegović, A.; Andričević, R. Above and in-water measurements of reflectance and chlorophyll algorithms in the Kaštela Bay in the Adriatic Sea. In Proceedings of the 2nd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing, WHISPERS 2010—Workshop Program, Reykjavik, Iceland, 14–16 June 2010; pp. 1–4. [Google Scholar] [CrossRef]
  31. Mozetič, P.; Solidoro, C.; Cossarini, G.; Socal, G.; Precali, R.; Francé, J.; Bianchi, F.; De Vittor, C.; Smodlaka, N.; Fonda Umani, S. Recent trends towards oligotrophication of the northern adriatic: Evidence from chlorophyll a time series. Estuaries Coasts 2010, 33, 362–375. [Google Scholar] [CrossRef] [Green Version]
  32. Cozzi, S.; Cabrini, M.; Kralj, M.; De Vittor, C.; Celio, M.; Giani, M. Climatic and anthropogenic impacts on environmental conditions and phytoplankton community in the gulf of trieste (Northern adriatic sea). Water 2020, 12, 2652. [Google Scholar] [CrossRef]
  33. Viličić, D.; Orlić, M.; Jasprica, N. The deep chlorophyll maximum in the coastal north eastern Adriatic Sea, July 2007. Acta Bot. Croat. 2008, 67, 33–43. [Google Scholar]
  34. Mozetič, P.; Malačič, V.; Turk, V. A case study of sewage discharge in the shallow coastal area of the Northern Adriatic Sea (Gulf of Trieste). Mar. Ecol. 2008, 29, 483–494. [Google Scholar] [CrossRef]
  35. Ninčević-Gladan, Ž.; Bužančić, M.; Kušpilić, G.; Grbec, B.; Matijević, S.; Skejić, S.; Marasović, I.; Morović, M. The response of phytoplankton community to anthropogenic pressure gradient in the coastal waters of the eastern Adriatic Sea. Ecol. Indic. 2015, 56, 106–115. [Google Scholar] [CrossRef] [Green Version]
  36. Bužančić, M.; Ninčević Gladan, Ž.; Marasović, I.; Kušpilić, G.; Grbec, B. Eutrophication influence on phytoplankton community composition in three bays on the eastern Adriatic coast. Oceanologia 2016, 58, 302–316. [Google Scholar] [CrossRef] [Green Version]
  37. Teodoro, A.C.; Duleba, W.; Gubitoso, S.; Prada, S.M.; Lamparelli, C.C.; Bevilacqua, J.E. Analysis of foraminifera assemblages and sediment geochemical properties to characterise the environment near Araçá and Saco da Capela domestic sewage submarine outfalls of São Sebastião Channel, São Paulo State, Brazil. Mar. Pollut. Bull. 2010, 60, 536–553. [Google Scholar] [CrossRef] [PubMed]
  38. Despland, L.M.; Vancov, T.; Aragno, M.; Clark, M.W. Diversity of microbial communities in an attached-growth system using BauxsolTM pellets for wastewater treatment. Sci. Total Environ. 2012, 433, 383–389. [Google Scholar] [CrossRef] [PubMed]
  39. Yu, J.; Seon, J.; Park, Y.; Cho, S.; Lee, T. Electricity generation and microbial community in a submerged-exchangeable microbial fuel cell system for low-strength domestic wastewater treatment. Bioresour. Technol. 2012, 117, 172–179. [Google Scholar] [CrossRef]
  40. Decision, C. Directive 2006/7/EC of the European Parliament and of the Council of 15 February 2006 concerning the management of bathing water quality and repealing Directive 76/160/EEC. Off. J. Eur. Union 2006, L064, 37–51. [Google Scholar]
  41. Di Dato, M.; Galešić, M.; Šimundić, P.; Andričević, R. A novel screening tool for the health risk in recreational waters near estuary: The Carrying Capacity indicator. Sci. Total Environ. 2019, 694. [Google Scholar] [CrossRef] [PubMed]
  42. Brettum, P.; Andersen, T. The use of phytoplankton as indicators of water quality. In Norwegian Institute for Water Research SNO Report; NIVA: Oslo, Norway, 2004; p. 33. [Google Scholar]
  43. Howarth, R.W.; Sharpley, A.; Walker, D. Sources of nutrient pollution to coastal waters in the United States: Implications for achieving coastal water quality goals. Estuaries 2002, 25, 656–676. [Google Scholar] [CrossRef]
  44. Halpern, B.S.; Walbridge, S.; Selkoe, K.A.; Kappel, C.V.; Micheli, F.; D’Agrosa, C.; Bruno, J.F.; Casey, K.S.; Ebert, C.; Fox, H.E.; et al. A global map of human impact on marine ecosystems. Science 2008, 319, 948–952. [Google Scholar] [CrossRef] [Green Version]
  45. Windom, H.L. Contamination of the marine environment from land-based sources. Mar. Pollut. Bull. 1992, 25, 32–36. [Google Scholar] [CrossRef]
  46. Lapointe, B.E.; Thacker, K.; Hanson, C.; Getten, L. Sewage pollution in Negril, Jamaica: Effects on nutrition and ecology of coral reef macroalgae. Chin. J. Oceanol. Limnol. 2011, 29, 775–789. [Google Scholar] [CrossRef]
  47. Kiefer, D.A. Fluorescence properties of natural phytoplankton populations. Mar. Biol. 1973, 22, 263–269. [Google Scholar] [CrossRef]
  48. Röttgers, R.; Koch, B.P. Spectroscopic detection of a ubiquitous dissolved pigment degradation product in subsurface waters of the global ocean. Biogeosciences 2012, 9, 2585–2596. [Google Scholar] [CrossRef] [Green Version]
  49. Sosik, H.M.; Chisholm, S.W.; Olson, R.J. Chlorophyll fluorescence from single cells: Interpretation of flow cytometric signals. Limnol. Oceanogr. 1989, 34, 1749–1761. [Google Scholar] [CrossRef]
  50. Kolber, Z.; Falkowski, P.G. Use of active fluorescence to estimate phytoplankton photosynthesis in situ. Limnol. Oceanogr. 1993, 38, 1646–1665. [Google Scholar] [CrossRef]
  51. Müller, P.; Li, X.P.; Niyogi, K.K. Update on Photosynthesis Non-Photochemical Quenching. A Response to Excess Light Energy. Plant Physiol. 2001, 125, 1558–1566. [Google Scholar] [CrossRef] [Green Version]
  52. Intergovernmental Oceanographic Commission; Scientific Committee on Oceanic Research; International Association for the Physical Sciences of the Ocean. The international thermodynamic equation of seawater—2010: Calculation and use of thermodynamic properties. In Intergovernmental Oceanographic Commission, Manuals and Guides No. 56; UNESCO: Paris, France, 2010; p. 196. [Google Scholar]
  53. Petersen, W.; Schroeder, F.; Bockelmann, F.-D. FerryBox—Application of continuous water quality observations along transects in the North Sea. Ocean Dyn. 2011, 61, 1541–1554. [Google Scholar] [CrossRef]
  54. Cunningham, A. Variability of In-Vivo Chlorophyll Fluorescence and Its Implication for Instrument Development in Bio-Optical Oceanography. Sci. Mar. 1996, 60, 309–315. [Google Scholar]
  55. Falkowski, P.; Kiefer, D.A. Chlorophyll a fluorescence in phytoplankton: Relationship to photosynthesis and biomass. J. Plankton Res. 1985, 7, 715–731. [Google Scholar] [CrossRef]
  56. Croatian Standards Institute (HRN EN ISO 9308-1:2014/A1:2017). Water Quality—Enumeration of Escherichia Coli and Coliform Bacteria—Part 1: Membrane Filtration Method for Waters with Low Bacterial Background Flora; Croatian Standards Institute: Zagreb, Croatia, 2017. [Google Scholar]
  57. Jozić, S.; Vukić Lušić, D.; Ordulj, M.; Frlan, E.; Cenov, A.; Diković, S.; Kauzlarić, V.; Fiorido Đurković, L.; Stilinović Totić, J.; Ivšinović, D.; et al. Performance characteristics of the temperature-modified ISO 9308-1 method for the enumeration of Escherichia coli in marine and inland bathing waters. Mar. Pollut. Bull. 2018, 135, 150–158. [Google Scholar] [CrossRef] [PubMed]
  58. Croatian Standards Institute (HR EN ISO 7899-2:2000). Water Quality—Detection and Enumeration of Intestinal Enterococci—Part 2: Membrane Filtration Method; Croatian Standards Institute: Zagreb, Croatia, 2000. [Google Scholar]
  59. Zhao, C.; Maerz, J.; Hofmeister, R.; Röttgers, R.; Wirtz, K.; Riethmüller, R.; Schrum, C. Characterizing the vertical distribution of chlorophyll a in the German Bight. Cont. Shelf Res. 2019, 175, 127–146. [Google Scholar] [CrossRef]
  60. Goosen, N.K.; Kromkamp, J.; Peene, J.; Van Rijswijk, P.; Van Breugel, P. Bacterial and phytoplankton production in the maximum turbidity zone of three European estuaries: The Elbe, Westerschelde and Gironde. J. Mar. Syst. 1999, 22, 151–171. [Google Scholar] [CrossRef]
  61. Ni, X.; Huang, D.; Zeng, D.; Zhang, T.; Li, H.; Chen, J. The impact of wind mixing on the variation of bottom dissolved oxygen off the Changjiang Estuary during summer. J. Mar. Syst. 2016, 154, 122–130. [Google Scholar] [CrossRef] [Green Version]
  62. Epifanio, C.E.; Maurer, D.; Dittel, A.I. Seasonal changes in nutrients and dissolved oxygen in the Gulf of Nicoya, a tropical estuary on the Pacific coast of Central America. Hydrobiologia 1983, 101, 231–238. [Google Scholar] [CrossRef]
  63. Frankignoulle, M.; Abril, G.; Borges, A.; Bourge, I.; Canon, C.; Delille, B.; Libert, E.; Théate, J.M. Carbon dioxide emission from European estuaries. Science 1998, 282, 434–436. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Kress, N.; Coto, S.L.; Brenes, C.L.; Brenner, S.; Arroyo, G. Horizontal transport and seasonal distribution of nutrients, dissolved oxygen and chlorophyll-a in the Gulf of Nicoya, Costa Rica: A tropical estuary. Cont. Shelf Res. 2002, 22, 51–66. [Google Scholar] [CrossRef]
  65. Gazeau, F.; Smith, S.V.; Gentili, B.; Frankignoulle, M.; Gattuso, J.P. The European coastal zone: Characterization and first assessment of ecosystem metabolism. Estuar. Coast. Shelf Sci. 2004, 60, 673–694. [Google Scholar] [CrossRef]
  66. Gazeau, F.; Gattuso, J.P.; Middelburg, J.J.; Brion, N.; Schiettecatte, L.S.; Frankignoulle, M.; Borges, A.V. Planktonic and whole system metabolism in a nutrient-rich estuary (the Scheldt estuary). Estuaries 2005, 28, 868–883. [Google Scholar] [CrossRef] [Green Version]
  67. Yuan, X.; Yin, K.; Harrison, P.J.; Cai, W.J.; He, L.; Xu, J. Bacterial production and respiration in subtropical Hong Kong waters: Influence of the Pearl River discharge and sewage effluent. Aquat. Microb. Ecol. 2010, 58, 167–179. [Google Scholar] [CrossRef] [Green Version]
  68. Stanley, D.W.; Nixon, S.W. Stratification and bottom-water hypoxia in the Pamlico River estuary. Estuaries 1992, 15, 270–281. [Google Scholar] [CrossRef]
  69. Breitburg, D.L.; Adamack, A.; Rose, K.A.; Kolesar, S.E.; Decker, M.B.; Purcell, J.E.; Keister, J.E.; Cowan, J.H. The pattern and influence of low dissolved oxygen in the Patuxent River, a seasonally hypoxic estuary. Estuaries 2003, 26, 280–297. [Google Scholar] [CrossRef]
  70. Bergondo, D.L.; Kester, D.R.; Stoffel, H.E.; Woods, W.L. Time-series observations during the low sub-surface oxygen events in Narragansett Bay during summer 2001. Mar. Chem. 2005, 97, 90–103. [Google Scholar] [CrossRef]
  71. Roth, F.; Lessa, G.C.; Wild, C.; Kikuchi, R.K.P.; Naumann, M.S. Impacts of a high-discharge submarine sewage outfall on water quality in the coastal zone of Salvador (Bahia, Brazil). Mar. Pollut. Bull. 2016, 106, 43–48. [Google Scholar] [CrossRef] [Green Version]
  72. Ramos, P.A.; Neves, M.V.; Pereira, F.L. Mapping and initial dilution estimation of an ocean outfall plume using an autonomous underwater vehicle. Cont. Shelf Res. 2007, 27, 583–593. [Google Scholar] [CrossRef]
  73. Petrenko, A.A. Effects of a sewage plume on the biology, optical characteristics, and particle size distributions of coastal waters. J. Geophys. Res. C Ocean. 1997, 102, 25061–25071. [Google Scholar] [CrossRef] [Green Version]
  74. Washburn, B.L.; Jones, B.H.; Bratkovich, A.; Dickey, D.; Chen, M. Mixing, Dispersion, and resuspension in vicinity of ocean wastewater plume. J. Hydraul. Eng. 1992, 118, 38–58. [Google Scholar] [CrossRef]
Figure 1. Locations of marine wastewater outfalls and national monitoring program stations.
Figure 1. Locations of marine wastewater outfalls and national monitoring program stations.
Water 13 02016 g001
Figure 2. Vertical profiles sampled during 2020 at location 1. The measured parameters are density anomaly (red line), buoyancy frequency (blue line), chlorophyll (light blue line), dissolved oxygen (magenta line), turbidity (green line), CDOM (black line), temperature (orange line), and salinity (brown line).
Figure 2. Vertical profiles sampled during 2020 at location 1. The measured parameters are density anomaly (red line), buoyancy frequency (blue line), chlorophyll (light blue line), dissolved oxygen (magenta line), turbidity (green line), CDOM (black line), temperature (orange line), and salinity (brown line).
Water 13 02016 g002
Figure 3. Vertical profiles sampled during 2020 at location 2. The measured parameters are density anomaly (red line), buoyancy frequency (blue line), chlorophyll (light blue line), dissolved oxygen (magenta line), turbidity (green line), CDOM (black line), temperature (orange line), and salinity (brown line).
Figure 3. Vertical profiles sampled during 2020 at location 2. The measured parameters are density anomaly (red line), buoyancy frequency (blue line), chlorophyll (light blue line), dissolved oxygen (magenta line), turbidity (green line), CDOM (black line), temperature (orange line), and salinity (brown line).
Water 13 02016 g003
Figure 4. Microbial water quality ((A) Escherichia coli, (B) Intestinal enterococci) at location 1. L1-B refers to location 1 bottom (blue), while L1-S refers to location 1 surface location (orange).
Figure 4. Microbial water quality ((A) Escherichia coli, (B) Intestinal enterococci) at location 1. L1-B refers to location 1 bottom (blue), while L1-S refers to location 1 surface location (orange).
Water 13 02016 g004
Figure 5. Microbial water quality ((A) Escherichia coli, (B) Intestinal enterococci) at location 2. L2-B refers to location 2 bottom (blue), while L2-S refers to location 2 surface location (orange).
Figure 5. Microbial water quality ((A) Escherichia coli, (B) Intestinal enterococci) at location 2. L2-B refers to location 2 bottom (blue), while L2-S refers to location 2 surface location (orange).
Water 13 02016 g005
Figure 6. Geometric mean vertical chlorophyll (A) and density anomaly (B) concentrations observed during the 2012–2019 period at the measuring station, FP-O14.
Figure 6. Geometric mean vertical chlorophyll (A) and density anomaly (B) concentrations observed during the 2012–2019 period at the measuring station, FP-O14.
Water 13 02016 g006
Figure 7. Geometric mean vertical chlorophyll (A) and density anomaly (B) concentrations observed during the 2012–2019 period at the measuring station, FP-O14b.
Figure 7. Geometric mean vertical chlorophyll (A) and density anomaly (B) concentrations observed during the 2012–2019 period at the measuring station, FP-O14b.
Water 13 02016 g007
Table 1. Available data for two WWTP locations, considering the average discharge (m3/day) and different wastewater effluent parameter concentrations (mg/L), presented as the mean ± standard deviation.
Table 1. Available data for two WWTP locations, considering the average discharge (m3/day) and different wastewater effluent parameter concentrations (mg/L), presented as the mean ± standard deviation.
Location/ParameterQSUMMER (m3/Day)QWINTER (m3/Day)COD (mg/L)BOD (mg/L)TSS (mg/L)TN (mg/L)TP (mg/L)
Stupe30.00026.000372 ± 235199 ± 11389 ± 5952 ± 156 ± 1.5
Katalinića brig38.00034.000246 ± 107135 ± 5869 ± 3832 ± 54 ± 2
Table 2. Vertical profile classification summary.
Table 2. Vertical profile classification summary.
Location 1Location 2
MonthStratificationEffluentTypeStratificationEffluentType
FebruaryXXUCPXXUCP
AprilXUCP-ECPXXUCP
MayXBCPXBCP
JuneECPECP
JulyECPECP
SeptemberECPXBCP
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kvesić, M.; Vojković, M.; Kekez, T.; Maravić, A.; Andričević, R. Spatial and Temporal Vertical Distribution of Chlorophyll in Relation to Submarine Wastewater Effluent Discharges. Water 2021, 13, 2016. https://0-doi-org.brum.beds.ac.uk/10.3390/w13152016

AMA Style

Kvesić M, Vojković M, Kekez T, Maravić A, Andričević R. Spatial and Temporal Vertical Distribution of Chlorophyll in Relation to Submarine Wastewater Effluent Discharges. Water. 2021; 13(15):2016. https://0-doi-org.brum.beds.ac.uk/10.3390/w13152016

Chicago/Turabian Style

Kvesić, Marija, Marin Vojković, Toni Kekez, Ana Maravić, and Roko Andričević. 2021. "Spatial and Temporal Vertical Distribution of Chlorophyll in Relation to Submarine Wastewater Effluent Discharges" Water 13, no. 15: 2016. https://0-doi-org.brum.beds.ac.uk/10.3390/w13152016

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