1. Introduction
Coastal waters are fragile ecosystems receiving environmental pressures from human activities. Mussel farming is one of these activities that contributes to the human pressures and potentially changes the natural environment. The
Gulf
of
Thessaloniki
(NW Aegean Sea) is a shallow coastal area influenced by municipal and industrial actions and effluents, shipping, fishing, as well as estuarine inputs from four major rivers. Mussel cultures both by long line and pole systems are located in the area and comprise an important economic activity [1] .
Primary production and assimilation of phytoplankton and detritus as food for mussel growth, through the filtration of water, are the key factors affecting the ecological environment of a mussel farm [2] . The metabolic activity of mussels increases detritus and organic by-products as feces and pseudo feces [3] while nutrients such as ammonium and phosphate are released by the bivalves [4] . The enrichment of the water column in ammonium, phosphate and silicate is also possible due to biodegradation or biodeposition to the bottom [5] . Nutrient cycling is also affected by feces and pseudo feces biodegradation in the water column that deliberates inorganic nutrients [6] . In this way, high density of bivalves in mussel cultures can affect nutrient dynamics in the coastal environment.
The development of indicators of nutrient impacts and eutrophication in coastal waters has received considerable effort [7] [8] [9] . The classification process is based on three major issues: nutrient enrichment, primary production and undesirable effects expressed as dissolved oxygen depletion [10] .
However, mussel farming differs from fish farming concerning the trophic enrichment of coastal waters. Mussels both produce and consume nutrients through filtering of water [4] thus being able of mitigating nutrients and resulting eutrophication in coastal areas. Sediment organic matter enrichment, underneath the mussel culture, though, activates higher nutrient regeneration rates [5] [11] increasing nutrient concentration in the deeper layer of the water column. The identification of potential impacts of dense aggregations of mussels on the trophic status of the coastal waters that host the activity of mussel farming has two aspects: 1) the support of the growth of mussels and 2) the prevention of undesirable effects on the water column and the quality of mussels.
The seasonal nutrient variation in a long-line of a mussel farm in the area of Chalastra (NW Thessaloniki Gulf) was investigated in this work, during the reproduction, growth and harvest cycle of mussels, in twelve samplings from July 2014 to April 2015. Sampling of water took place in four sites of the long-line, with various densities of cultured mussels and in a reference station, outside the culture. Additional samples were collected monthly in the period May-August 2015, in two sites of different mussel density. The results are depicting the variation of nutrients in the culture’s micro-environment in an attempt to identify different impacts of the culture’s density in the trophic status of the coastal waters of
Thessaloniki
’s Gulf. Multivariate statistical analysis and eutrophication indicators [9] [10] were used as tools for the identification of the variation schemes of nutrients, chlorophyll α and dissolved oxygen in the cultures of different mussel density.
2. Materials and Methods
Seawater samples were collected by a van Dorn water sampler, from three depths, 1 m , 5 m and at the bottom layer, close to 10 m , from a long-line in a mussel culture (Figure 1) at sites of different mussel density (Figure 2). The densities were determined by the distance among the suspended mussel socks. Four sampling sites were selected at 30 cm , 50 cm , 70 cm and 90 cm distance (Figure 2). Samples were also collected from a reference site, outside the farm at a distance of 2.5 km (Figure 1). Seawater samples were filtered on-board through membrane 0.45 μm filters [12] and kept in −18˚C until their analysis. Nitrate, nitrite, phosphate and silicate values were measured spectrophotometrically [13] . Dissolved oxygen was measured by the Winkler method on board [14] , while chlorophyll α was determined after the extraction with acetone on a double beam spectrophotometer [13] . Analysis of ammonium was performed by
Figure 1. Map with the location of the experimental long-line of mussels (M) and the reference station (A) at a distance of 2.5 km southeast (adapted from google earth).
Figure 2. Layout of the production line with the different distances between the suspended mussel socks and their vicinity to the bottom.
the phenolhypochlorite method [15] .
Data were treated using a three factor analysis of variance (seasonality, water depth and distance between socks) on the measured parameters (nutrients, chlorophyll α and dissolved oxygen) inside the long line mussel culture. Statistical differences between factor levels were further examined in the graphs of the parameters correlations and sample scores produced by the Principal Component Analysis (PCA) which was applied on the whole data set of the study.
Water Quality Indices
TRIX is an eutrophication index [10] for the evaluation of trophic levels in coastal waters. It is a linear combination of the logarithm of four variables, chlorophyll α, total inorganic nitrogen, inorganic phosphorus and the absolute percentage of the deviation from the dissolved oxygen saturation.
CDIN is the concentration of total dissolved inorganic nitrogen (NO3 + NO2 + NH4) and D%O2 is the absolute percentage of the deviation from the dissolved oxygen saturation. This index was first applied to describe eutrophication status of Italian coastal waters (Adriatic Sea) but has been used in other coastal areas as well [8] .
EI eutrophication index derived from the application of principal component analysis (PCA) on a nutrient and chlorophyll α values data set from Greek oligotrophic waters [9] and is suggested as follows:
3. Results
The variation of phosphates in the period July 2014 to April 2015, in the depths 1, 5 and 10 m , in the reference station and the four different handlings, is presented in Figure 3. The concentration of phosphates ranged from 0.01 to 0.58 μΜ (Table 1). The highest concentrations were measured at the depth of 5 m , where the mussels are suspended. All nutrients showed expectedly great seasonal variability in concentration (Table 2). Phosphates were more sensitive during the sampling period, distance of mussels and deep waters (Table 2).
Nitrate values were similar in the reference station and among the handlings (Figure 4), with a concentration range 0.02 - 8.61 μM (Table 1). Higher concentrations were presented in November and in the end of winter and spring months (Figure 4). Nitrite values ranged from 0.01 to 2.02 μΜ, with a mean value of 0.35 μΜ and were high in October and November as well as in the spring months (Figure 5).
The concentration of ammonium presented statistically important differences between the sampling sites (Table 2). Ammonium values ranged from 0.01 to 42.11 μΜ (Figure 6). The higher concentrations were observed in October and November, in the spring months and in the summer, in the bottom layer.
Figure 3. Variation of phosphates during the sampling period July 2014-April 2015, in the depths 1, 5 and 10 m in the four sites of different distance between the mussel socks and the reference station (distance “ 0” ).
Figure 4. Variation of nitrates during the sampling period July 2014-April 2015, in the depths 1, 5 and 10 m in the four sites of different distance between the mussel socks and the reference station (distance “ 0” ).
Figure 5. Variation of nitrites during the sampling period July 2014-April 2015, in the depths 1, 5 and 10 m in the four sites of different distance between the mussel socks and the reference station (distance “ 0” ).
Figure 6. Variation of ammonium during the sampling period July 2014-April 2015, in the depths 1, 5 and 10 m in the four sites of different distance between the mussel socks and the reference station (distance “ 0” ).
Table 1. Mean and range of values of nutrients and chlorophyll α, from July 2014 to August 2015 in the sampling sites (cultures of different mussel density and reference site).
Table 2. Statistically significant effects (in bold) of seasonal and vertical variation (depth) and sock distance on dissolved oxygen, nutrients and chlorophyll α, at 0.5 probability reference level.
Silicates ranged from 0.28 - 82.76 μΜ (Table 1) and presented low variation among the handlings and the reference station (Figure 7).
Differences in chlorophyll α values were observed in the depths of 5 and 10 m (Figure 8), showing a reduction of about forty per cent during the spring bloom in the handlings compared to the reference site.
Nutrients, dissolved oxygen and chlorophyll α values were used for the calculation of eutrophication indices, TRIX and EI [9] [10] . The environmental quality status was classified in five ranges: high, good, moderate, poor and bad, corresponding to a quality score from 1 to 5, respectively.
Table 3 presents the correlation coefficients between nutrients, dissolved oxygen and chlorophyll a values and PCA axes 1 and 2. The first two PCA axes explain 59.8% of the total variation of the whole data set. Phosphate and N/P show the highest correlation coefficients and therefore are the most important variables for the formation of axes 1 and 2 correspondingly. The variables total N and ammonium, share equal contribution for the formation of both axes. All the variables under study are shown in Figure 9, forming different directions that easily can be grouped into two bundles of variables. Phosphate with chlorophyll α form the first bundle together with nitrate and nitrite, indicating a positive correlation between higher primary productivity and nutrients. The
Figure 7. Variation of silicates during the sampling period July 2014-April 2015, in the depths 1, 5 and 10 m in the four sites of different distance between the mussel socks and the reference station (distance “ 0” ).
Figure 8. Variation of chlorophyll α values during the sampling period July 2014-April 2015, in the depths 1, 5 and 10 m in the four sites of different distance between the mussel socks and the reference station (distance “ 0” ).
Figure 9. PCA biplot of variable correlation with principal axes. Longer arrows indicate higher correlation values, while acute angles indicate positive correlation and obtuse negative ones, between each pair of variables.
Table 3. Correlation coefficients between the physicochemical variables and PCA axes 1 and 2 and percentage variance determination.
second bundle is formed by ammonium and total nitrogen, characterized by high production rates of ammonium. The ratio N:P is negatively correlated with the first bundle of variables, indicative of an irregular seasonal, spatial and vertical variation of the Redfield ratio at the culture’s micro-environment.
In Figure 10, the PCA ordination plot of the variable sample scores, showing the distances between the socks in both sampling periods (July 2014-August 2015) is presented. If we superimpose the Figure 8 and Figure 9, a useful relationship between variables and sample sites reveals. Sites close to particular nutrients show a mutual effect leading to a positive relation: sites with high concentrations of the corresponding nutrients. Marked as “ 10” is the reference station and as “ 100” the dense culture of 30 cm distance. There is a clear distinction of the dense culture (bottom right part) and the less dense culture of 70 cm
Figure 10. Ordination plot of the variable sample scores, showing the distances between the socks in both sampling periods (July 2014-August 2015). The two polygons enclose the samples belonging to the second period of the study.
distance. The less dense culture in the second sampling period comes along with the reference site. In the first sampling period, the densities of 30 and 50 cm distance are not clearly distinguished between each other but still, they have different characteristics from the 70 and 90 cm distances.
Figure 11 shows the position of depth points in the various sampling periods. Depth determines variation of ammonium and total inorganic nitrogen, indicating enrichment in the bottom layer, while chlorophyll α and phosphate concentration correlate well at the upper layer of the water body (1 and 5 m depth).
TRIX index incorporates nutrients as the driving forces of eutrophication, chlorophyll α as the biological response and deviation of dissolved oxygen saturation as the environmental disturbance factor. Eutrophication index (EI), derived from oligotrophic Greek coastal waters [9] is formed by nutrients and chlorophyll α concentration. TRIX values ranged from 0.02 to 4.16, with a mean value of 1.64, while EI index ranged from 0.08 to 12.9, with a mean value of 3.05. The results were classified in five trophic quality ranges, for TRIX, high (<1.6), good (1.6 - 2.8), moderate (2.8 - 4), poor (4 - 5.3) and bad (>5.3) and for EI, high (<0.04), good (0.04 - 0.38), moderate (0.38 - 0.85), poor (0.85 - 1.51) and bad (>1.51), given a quality score from 1 to 5, respectively. For comparison reasons, the trophic range was the same applied for Greek coastal waters [9] [16] .
Figure 12 presents the ordination plot of nutrients, dissolved oxygen, total N, N:P, Chlα, along with TRIX eutrophication index entered as a supplementary
Figure 11. Ordination plot of the variable sample scores, showing the distribution of depth in the mussel culture in both sampling periods (July 2014-August 2015). The two polygons enclose the samples belonging to the second period of the study.
Figure 12. PCA biplot of sampling sites as arrayed by the TRIX classification index (quality scores from 1 - 5). The two polygons enclose the samples belonging to the second period of the study.
Figure 13. PCA biplot of sampling sites as arrayed by the EI classification index (quality values from 1 - 5). The two polygons enclose the samples belonging to the second period of the study.
quality value. The dense culture ( 30 cm distance) is distinguished (at the bottom right) from the 70 cm (upper part), classified in the middle (score 2) eutrophic scale. However, an unclear pattern (upper right part) of the classification of the four densities in the first sampling period is presented.
Figure 13 shows the PCA ordination plot of nutrients, dissolved oxygen, total N, N:P, Chla and EI eutrophication index entered as a supplementary quality value. The application of this index, also confirms that the dense culture ( 30 cm distance) is distinguished (at the bottom right) from the 70 cm (upper part). The classification range that was used, initially selected for oligotrophic waters [9] classifies the mussel’s culture coastal environment in the worse (score 5) eutrophic scale.
4. Discussion
In this work, nutrient values were measured in an experimental long-line mussel culture in the area of Chalastra (NW Thessaloniki Gulf), where intensive mussel culture is an important economic activity. The study was conducted throughout a year (July 2014-August 2015). Two sampling periods can be identified, one assessing the parameters in all handlings (July 2014-April 2015) and the other assessing the parameters in selected handlings (30 and 70 cm between the socks). The study was an attempt to identify the different impacts of the culture density in the trophic conditions of the coastal waters of
Thessaloniki
’s Gulf. The density of the culture was defined by the distance between the mussel socks. In the first study period, four mussel handlings of different sock distance were applied, 30, 50, 70 and 90 cm , while in the second period (May-August 2015), two distances were selected, 30 and 70 cm . Samples from a reference site outside the culture were also collected and analyzed.
Statistically important differences among the different handlings and the reference site were identified in the values of phosphates, nitrites and chlorophyll α as well as in the values of dissolved oxygen. On the contrary, nitrates and ammonium presented unimportant variation in the water column, inside the mussel culture and in the reference site, indicating nutrient mitigation from the mussels [6] . However, the fluctuation of the parameters had a different picture in the deeper layer, close to the bottom, where the prevailing ecological and environmental conditions need further evaluation. Phosphate concentrations increased with depth, suggesting enrichment from the bottom, probably due to resuspension and deliberation of phosphates from the sediments, but remained low in the sampling period. The lowest levels of phosphates were observed in December. Statistically significant differences among the handlings were observed for the phosphate values, indicating correlation between increased metabolic rate and phosphate enrichment [17] [18] . Depth was an important parameter for both phosphate and ammonium accumulation. The hydrodynamic circulation, investigated by Savvidis et al. [19] was characterized by low currents of average speed < 5 cm/s.
The presence of higher concentrations of ammonium and phosphate in the bottom layer of the culture has been also depicted in a monitoring study by Moriki et al. [20] inside the mussel farms at Chalastra area. In the same study, the correlation of elevated nitrite concentrations and low dissolved oxygen values has also been reported [20] , along with low spatial and seasonal variation of nitrate.
The coastal waters of Thessaloniki’s Gulf are affected by the estuaries of three major rivers, Axios, Loudias, Aliakmon as well as municipal, agricultural and industrial effluents from the metropolitan area of Thessaloniki. The nutrient status of the Gulf is constantly under monitoring from the local authorities, in the context of the European Directives, Water FD 2000/60/EC and Marine Strategy FD 2008/56/EC [21] [22] . According to HCMR [23] low dissolved oxygen values along with elevated nutrient levels and low values of current velocities and therefore water circulation have an impact to the eutrophic status of the area, with the contribution of the bivalve cultures in the NW of the Gulf. In this work, it was observed that the mussel culture has a low impact on nitrate and silicate values, while severe oxygen depletion along with elevated nitrite and ammonium values were observed in the bottom layer of the mussel culture. The density of the culture plays an important role in the accumulation of ammonium and nitrites, as observed in the second period of study (May to August 2015).
Multivariate statistical analysis was applied on the data set of both sampling periods, that included seven physicochemical parameters, measured at four sock distances (mussel density) in three vertical points of the water column during the two periods (July 2014-April 2015) and later on (May to August 2015) in two sock distances, corresponding to less dense and dense culture and was a valuable tool in identifying the variation scheme of the measured parameters. The positions of the variables on the PCA axes revealed the correlation of phosphate with chlorophyll α along with nitrate and nitrite, indicating the higher primary productivity and nutrients correlation, while ammonium concentrations played an important role in the variation of total dissolved inorganic nitrogen. Therefore, the driving forces for eutrophication are mainly phosphate, nitrate and nitrite as revealed by the PCA. High ammonium concentrations are accumulated in the mussel culture waters, especially in the deeper layer, close to the bottom. The distinction of different impacts of the mussel density in the eutrophication status was complicated in the first sampling period, when four mussel densities were investigated. However, data analysis, allowed the observation of the seasonal variation of the eutrophication parameters inside the mussel culture, during the mussel life cycle, from reproduction to growth and harvest. Chlorophyll α is an important environmental parameter that represents the phytoplankton biomass both as a eutrophication variable and as a measure of food availability in the mussel’s culture complicated micro-environment. Low chlorophyll α values may have a negative impact in the development of eutrophication and therefore a positive impact in environmental quality, but reflect low food levels to support mussel growth in the culture. PCA revealed higher chlorophyll α values in the reference station, outside the culture, but the identification of the balanced relationship between chlorophyll α and mussel growth needs to be investigated.
Two eutrophication indices, TRIX and EI were used as a supplementary quality value for the classification of the sampling sites at the different mussel densities. Both indices were sensitive in identifying differences seasonally and spatially. However TRIX is a better environmental quality indicator, since clearly showed differences among the two selected mussel densities, 30 and 70 cm distance between the socks, in the period May-August 2015.
Nutrients, chlorophyll α and dissolved oxygen values were used as environmental quality indicators in the coastal waters of
Thessaloniki
’s Gulf, inside and outside an experimental long-line of mussels in suspension. The parameters were undertaken multivariate statistical analysis in order to identify the implication of the culture in the trophic status of the coastal waters. PCA depicted the variation that is influenced by the mussel’s density.
5. Conclusions
The investigation of nutrient variability in a long-line of suspended mussel culture of different densities showed clear differences between the dense ( 30 cm distance of the socks) and the less dense culture ( 70 cm distance of the socks), especially in the bottom layer, while the TRIX and EI eutrophication indicators proved in a confirmatory matter the overall prevailing conditions on a seasonal basis.
The mussel density had implications on water quality, while food availability in terms of chlorophyll α was seriously affected compared to the reference station.
Acknowledgements
This research was partially funded by the operational programme “Fisheries 2007- 2013”.