Integrated Monitoring with Moss-Bag and Mussel Transplants in Reservoirs

: For the ﬁrst time, transplants with moss-bags and mussels together were applied to study the water quality in standing water bodies. The tested species: Fontinalis antipyretica Hedw. and Sinanodonta woodiana (Lea, 1834) were collected from unpolluted sites and analyzed to obtain background levels. Then, the moss and mussels were left in cages for a period of 30 days in three reservoirs where both are not present naturally. Two of the reservoirs su ﬀ er from old industrial contamination and one is a ﬀ ected by untreated wastes. Twenty-four compounds were studied, among them trace elements Al, As, Cd, Co, Cr, Cu, Fe, Hg, Mn, Ni, Pb, Zn and organic priority substances: six polybrominated diphenyl ethers (PBDEs) congeners and short-chain chlorinated para ﬃ ns (SCCPs). The trace element accumulation was signiﬁcant after the exposition period in all studied stations. PBDEs and SCCPs were also accumulated up to two times more in the moss tissues. PBDEs in the mussels exceeded the environmental quality standard (EQS). The applied combined transplants, and especially the moss-bags, revealed severe contamination with heavy metals not detected by the water samples. The moss and the mussel followed a di ﬀ erent model of trace element and PBDEs accumulation. The SCCPs levels were alarmingly high in all plant samples. The study conﬁrmed PBDEs and SCCPs as bioaccumulative compounds and suggested that an EQS for SCCPs in biota needs to be established.


Introduction
Heavy metals are among the most common contaminants in the aquatic environment [1] and due to their excessive accumulation, biomagnification, and toxicity, high levels of heavy metals in surface waters have evoked significant concern from governments and the public [2]. Polybrominated diphenyl ethers (PBDEs) due to their persistence, bioaccumulation, long-range transport, toxicity, and adverse effects on wildlife and humans have been restricted or banned in most countries [3,4]. On the other hand, PBDEs are mainly used in plastic, textile, and electronic applications as fire retardants [5,6]. Short-chain chlorinated paraffins (SCCPs) are used for various industrial applications, such as flame retardants, plasticizers, metal-working fluids, lubricant additives, paints, sealants, and leather fat liquors. Due to their environmental persistence, potentials for bioaccumulation, and long-range atmospheric transport, in 2017, the Stockholm Convention listed SCCPs as a new group of persistent organic pollutants (POPs) [7,8]. Moreover, compared with medium-chain CPs (MCCPs, C14-17) and long-chain CPs (LCCPs, C18-30), short-chain CPs (SCCPs, C10-13) have higher toxicity effects on aquatic organisms, slower biodegradation rates, and, therefore, the most significant potential for long-range transport, persistence, and bioaccumulation [9,10]. There is very limited data on the environmental exposure of PBDEs and SCCPs in aquatic environments, especially on the potential accumulation in aquatic mosses and mussels. A review of PBDEs in the global environment recommended that more research needs to be performed on PBDEs in a variety of matrices and locations in order to assess the current status of these compounds [11]. SCCPs were found to be bioaccumulated by invertebrates and fish from Lake Ontario and Lake Michigan [12], and their possible biomagnification was also indicated [13].
Trace element and persistent organic pollution in aquatic ecosystems is traditionally monitored using chemical analyses of water and sediments [14,15]. Moreover, heavy metals and organic contaminants can bioaccumulate in aquatic organisms depending on their concentration in the waterfront, interaction with suspended organic matter, and bioavailability in sediments [16]. Therefore, biomonitoring can be used to determine the toxicant fraction in the environmental matrices (e.g., water and sediment) that is of ecotoxicological significance [17]. According to Azizi et al. [18], monitoring programs use a high number of bioindicators known as "sentinel organisms" to detect the temporal and spatial variation of contaminants and to contribute to the knowledge of trends in aquatic contamination.
The moss-bag technique was widely used to monitor air pollution. However, this simple and economical biomonitoring approach is rarely used to evaluate the aquatic environment quality, especially in lentic ecosystems. In comparison with the chemical monitoring of the pollutant concentrations in the water, the bryomonitoring technique allows the detection of intermittent or sporadic contamination, precise location of contamination sources, and simultaneous monitoring of a large number of contaminants (inorganic and organic) via analysis of a single sample [19]. Furthermore, active bryomonitoring can be used in environments where native mosses are absent, and it improves interpretation of the temporal variation in contaminants by taking into consideration the initial concentrations and time of exposure. The first transplants were applied fifty years ago [20]. Transplants of Fontinalis antipyretica Hedw. have been used to assess river pollution, mainly with heavy metals [21,22]. Fontinalis antipyretica is widespread throughout the northern hemisphere, and it is the most easily recognized aquatic species, as well as it is the most commonly used moss to monitor water quality [23].
Due to their sedentary filter-feeding, geographical distribution, size, and high contaminant bioconcentration facts, mussels are also widely used for assessing the ecological risk of different contaminants [24]. In fact, the Mussel Watch [25,26] is the oldest biomonitoring program in progress worldwide, and it has been successfully developed in many countries to determine the occurrence and the toxic effects of heavy metals in aquatic ecosystems and also human health [27][28][29][30][31]. The Chinese pond mussel, Sinanodonta woodiana (Bivalvia: Unionida: Unionidae) is a widely distributed benthic, freshwater species, with a native range that expands into East Asia from the Russian Far East (Amur Basin) to Taiwan, Indo-China, and Malaysia [32,33]. Currently, the species has widely spread across many regions, e.g., Europe, the USA, Hispaniola, Costa Rica, Indonesia, and the Philippines [34,35]. Moreover, according to Kolarević et al. [36], it fulfills the main criteria required for a bioindicator organism. Earlier studies have revealed the ability of Sinanodonta woodiana (S. woodiana) to accumulate trace elements [37] and pesticides [38], as well as its potential to detect genotoxicity [39], lysosomal membrane destabilization, and respiration rate [40,41]. Lastly, S. woodiana is a traditionally edible species in its native range, including China [42], and the monitoring of different toxicants is essential in terms of human health.
Although active biomonitoring has proven its effectiveness, until now, moss-bags were not applied in standing water bodies, except for a research in a Cu-contaminated lake in Italy [43]. Furthermore, to our knowledge, no studies on simultaneously transplanted mosses and mussels have been performed except a single study on trace metal contamination in rivers in France [44].
Therefore, in the present paper, we aimed to see whether contamination with trace elements and organic contaminants in reservoirs can be detected with the moss-bags technique and caged mussels. For this reason, alive individuals of Fontinalis antipyretica (F. antipyretica) and Sinanodonta woodiana from unpolluted sites were exposed during the summer of 2019 at three reservoirs influenced by industry and untreated wastes in order to study the contamination patterns in both, water and biota.

Field Design
Individuals of the aquatic moss F. antipyretica were collected from a small unpolluted stream (42.119625 • N, 24.556055 • E WGS 84) in Bulgaria. The material was rinsed on-site with the stream water before being transferred to the laboratory. Basal material, material in poor condition, and epiphytes and other foreign objects were removed. The material was divided into three parts for the three studied reservoirs. Flat bags (30 × 20 cm) were prepared with plastic mesh (aperture 1 mm), which was previously rinsed in distilled water and air-dried.
Mussels from one size class (154 ± 5.5 g; 11 ± 3.5 cm) were collected manually from a pool of the Institute of Fisheries and Aquaculture (Bulgaria), where abiotic factors are monitored regularly, and there are no sources of anthropogenic contamination. The mussels were transported to the laboratory on the same day in clean plastic containers filled with pool water. They were left in an aquarium filled with dechlorinated water and pool water (50:50) and fitted with air pumps until the field experiment started. Biometric analyses were performed before and after the exposition according to Gasmi et al. [45]. The total weight (Ww, g) of each mussel was measured after cleaning the outer shells from particles with an analytical balance (Kern, Germany) in g. The shells were measured with a caliper in mm.
The moss-bags (~100 g of F. antipyretica wet weight) together with the mussels (n = 10) were placed in cages (30 × 15 × 10 cm; galvanized steel) at a depth of 2 m in each of the reservoirs. The cages were removed from the reservoirs 30 days after the start of the experiment. Water samples from a depth of 2 m were taken simultaneously for inorganic and organic chemical analyses. The mussels' total soft tissues were dissected out for further chemical analyses according to the adapted for mussels EMERGE protocol (European Mountain lake Ecosystems: Regionalisation, diaGnostic and socioeconomic Evaluation) [46]. The exposure sites were located at three large, deep reservoirs (maximum depth between 48 and 98 m) with different uses (energy and water supply, irrigation, fish farming) and affected by a different type of pollution ( Figure 1). In such conditions, aquatic bryophytes are absent. Kardzhali (41. 25.519459 • E WGS 84) reservoirs suffer from old industrial contamination, especially ore extraction and mining activities [47], and have sediment loads of priority substances. Zhrebchevo Reservoir (42.585571 • N, 25.885592 • E WGS 84) is located in a region with light industries but is under the influence of illegal dumping sites (solid waste) and untreated wastewaters. Another reason for selection of the three heavily modified water bodies is that biomonitoring with S. woodiana in human-made lakes minimizes the potential threats caused by accidentally introduced species. As was mentioned above, all three studied reservoirs belong to the type large, deep reservoirs with concrete walls surrounding the prevailing part of their area. The only possible stations for positioning of the cages were shallow areas, which are used for fish farming.

Analytical Procedures
All samples were transported in coolers to the laboratories on the day of the sampling. The plant material was rinsed on-site with the reservoir water before being transferred to the laboratory, and material in poor condition and foreign objects were removed. Two to three-centimeter tips of F. antipyretica and whole-tissue samples from the mussels were analyzed. The moss samples were mineralized in a microwave digestion system (Ethos 1, Milestone, Denmark) following the procedure described in a previous work [48]. Briefly, the homogenized air-dried material (~1 g) was treated with 8 mL of HNO3 (65%, Merck, Darmstadt, Germany supra pure) overnight, then 2 mL of H2O2 (30%, Merck, Darmstadt, Germany, pa) were added, containers were sealed and irradiated. Fresh frozen mussel samples were homogenized by grinding, and ~2 g of the material was directly subjected to MB dissolution with the same mixture of HNO3 and H2O2 (8:2). After digestion, the samples were transferred in volumetric flasks (50 mL) and diluted with double distilled water (Fison Scientific Equipment, Fi-stream). For determination of trace elements, the sample solutions were measured directly, while for the elements with higher concentrations, such as Al, Ca, Fe, K, Mg, Mn, Na, and P, the sample solutions were additionally diluted by a factor of 5. The water samples were stabilized in 0.5% HNO3. Blank samples were prepared for each batch of samples. The elements Na (589.592 nm); K (766.490 nm); Mg (285.213 nm); Ca (422.673 nm), and Fe (238.204 nm) were determined by ICP-AES using iCAP 6300 Duo, S (Thermo Scientific, Waltham, MA, USA) using radial observation mode.
The elements and corresponding monitored isotopes 31 P; 52, 53 Cr; 59 Co; 60, 62 Ni; 63, 65 Cu; 64,66 Zn; 75 As; 111,113 Cd; 206, 208 Pb; and 201, 202 Hg were determined by ICP-MS (Agilent 7700, Agilent Technologies, Tokyo, Japan). Two isotopes per element (where possible) were monitored for dynamic evaluation of potential spectral matrix interference, and He (4.8 mL/min) was used as collision gas in order to prevent polyatomic ion formation. Additionally, two tungsten isotopes (182, 184 W) were permanently monitored in order to ensure dynamic assessment of the risk of overlapping the mercury signals by tungsten oxide radicals. In order to assure quantification of very

Analytical Procedures
All samples were transported in coolers to the laboratories on the day of the sampling. The plant material was rinsed on-site with the reservoir water before being transferred to the laboratory, and material in poor condition and foreign objects were removed. Two to three-centimeter tips of F. antipyretica and whole-tissue samples from the mussels were analyzed. The moss samples were mineralized in a microwave digestion system (Ethos 1, Milestone, Denmark) following the procedure described in a previous work [48]. Briefly, the homogenized air-dried material (~1 g) was treated with 8 mL of HNO 3 (65%, Merck, Darmstadt, Germany supra pure) overnight, then 2 mL of H 2 O 2 (30%, Merck, Darmstadt, Germany, pa) were added, containers were sealed and irradiated. Fresh frozen mussel samples were homogenized by grinding, and~2 g of the material was directly subjected to MB dissolution with the same mixture of HNO 3 and H 2 O 2 (8:2). After digestion, the samples were transferred in volumetric flasks (50 mL) and diluted with double distilled water (Fison Scientific Equipment, Fi-stream). For determination of trace elements, the sample solutions were measured directly, while for the elements with higher concentrations, such as Al, Ca, Fe, K, Mg, Mn, Na, and P, the sample solutions were additionally diluted by a factor of 5. The water samples were stabilized in 0.5% HNO 3 . Blank samples were prepared for each batch of samples. The elements Na (589.592 nm); K (766.490 nm); Mg (285.213 nm); Ca (422.673 nm), and Fe (238.204 nm) were determined by ICP-AES using iCAP 6300 Duo, S (Thermo Scientific, Waltham, MA, USA) using radial observation mode.
The elements and corresponding monitored isotopes 31 P; 52, 53 Cr; 59 Co; 60, 62 Ni; 63, 65 Cu; 64, 66 Zn; 75 As; 111,113 Cd; 206, 208 Pb; and 201, 202 Hg were determined by ICP-MS (Agilent 7700, Agilent Technologies, Tokyo, Japan). Two isotopes per element (where possible) were monitored for dynamic evaluation of potential spectral matrix interference, and He (4.8 mL/min) was used as collision gas in order to prevent polyatomic ion formation. Additionally, two tungsten isotopes (182, 184 W) were permanently monitored in order to ensure dynamic assessment of the risk of overlapping the mercury signals by tungsten oxide radicals. In order to assure quantification of very low concentration of Hg in the tested samples, the observation time for both isotopes was increased up to 0.5 s.
The elements Al and Mn were measured by both methods (ICP-AES and ICP-MS) at the corresponding emission lines (Mn 259.373 nm and Al 396.152 nm) and isotopes (27 Al and 55 Mn).
The non-spectral matrix effect and sensitivity drift inherent for the ICP-MS method were corrected by internal standardization. For this purpose, 103 Rh was selected as the most appropriate internal standard and was introduced to the calibration solutions as well as to the measured samples.
For calibration, the following traceable to NIST standard solutions were used after appropriate dilution: ICP multi-element standard solution IV, Certipur (Merck, Darmstadt, Germany); single-element standard solutions for K, Na, P, As, Hg, and Rh (CPAChem, Stara Zagora, Bulgaria-France).
For validation of the developed analytical methods, a certified reference material NCS DC 73,348-"Bush Branches and Leaves" was applied, and the obtained recoveries for all tested elements varied between 93.8% for Cr and 109% for Al.
The limits of detection (LOD) for Al and Hg were <0.005, for As <1, for Cd and Fe <0.1 and for Co < 0.01 mg L −1 in water. The concentrations in water and biota were presented in %, µg L −1 , and mg L −1 for water samples and in mg kg −1 for biota, respectively. Six PBDE congeners (BDEs 28, 47, 99, 100, 153, and 154) and SCCPs were analyzed in the extracted matrices (water and biota) by gas chromatography-mass spectrometry (GC-MS) with electron capture negative ionization (Thermo TSQ  8000  Before the analysis of PBDEs and SCCPs in water, 0.1 µg mL −1 deuterated standard was added to the samples, and solid-phase extraction was applied. The extracts were concentrated with a rotary vacuum evaporator to 1 mL. The injection volume was 1 µL, and the injection was performed in the splitless mode. For PBDEs, helium was employed as carrier gas (0.9 mL min −1 ), and Ar was used as the collision gas. The GC oven temperature was programmed at 80 • C (1 min), from 20 to 140 • C (0 min), and from 8 to 310 • C (5 min); while for the determination of SCCPs, methane was used as carrier gas (1.3 mL min −1 ) at the following temperature program: 100 • C (0.5 min), from 10 to 150 • C (0.5 min), and from 30 to 280 • C (7 min).
The sample preparation of biota samples (moss and mussels) involved a grinding step followed by extraction with ethyl acetate. The extract was concentrated and purified by adsorption chromatography, resulting in PBDEs being separated from the chlorinated compounds. The PBDE-containing fraction was concentrated and subjected to GC-MS analysis at the same conditions as mentioned above.
The analyzed sum of SCCPs includes the variety of SCCPs with their differing chlorine content and C-number distribution patterns following ISO 12010:2019.
The LOD for the PBDE congeners were 0.004 in water and 0.003 µg L −1 in biota; and for SCCPs 0.15 in water and 1.5 µg L −1 in biota. The concentrations in water and biota were presented in µg L −1 and µg kg -1 , respectively.

Contamination Factors and Metal Pollution Index
Metal enrichment in moss tissues was established based on contamination factor (CF), calculated as the ratio between the metal content at a given reservoir and the unpolluted stream. The following scale based on the geometric progression adopted by Mouvet [49] was applied: The threshold discriminating between suspected and certain pollution (CF = 6) is high enough to minimize the effects of possible variability sources (e.g., interspecific differences in cation exchange capacity, intraspecific differences due to the vertical position in the water column or to the horizontal position in the streambed, seasonal changes, water chemistry, sample treatment, analytical errors).
An evaluation of the overall metal contamination in moss samples was carried out by calculating a metal pollution index (MPI) derived as a linear sum of the contamination factors weighted to take into account the differences in toxicity of the various metals [50], i.e., where CFi is the contamination factor for metal i; wi is the weight for metal i; and wt = Σiwi.
The bioaccumulation factor (BAF) for trace elements and organic contaminants for the mussels was calculated after Nikanorov [53] according to the formula shown below: BAF = (toxicant concentration in mussel total soft tissue)/(toxicant concentrations in water).
In general, the BAF calculations show how many times more the concentration of the toxicant in biota is compared to that in the water. According to its values, the biota could be categorized as macro (BAF > 2), micro (BAF = 1 to 2), and deconcentrators (BAF < 1).

Statistical Analysis
The similarity of the accumulated pollutants in the aquatic moss and mussels was assessed using cluster analysis (unweighted pair-group average, Bray-Curtis similarity index). All statistical tests were performed using the computer software "PAST" [54]. A level of probability p ≤ 0.05 was accepted as statistically significant. The relationship between the analyzed toxicants in water, moss, and mussel samples was studied with principal component analysis (PCA). The data were transformed (x' = log (x + 1)), automatically centered, and standardized by the Canoco v.5 program [55].

Water
The water in all the reservoirs was alkaline, with pH values between 6.84 and 8.29. The concentrations in water samples were below the LOD for Hg (<0.05 µg L −1 ) and Fe (<0.1 mg L −1 ) at all sampling locations ( Table 1) All studied PBDE congeners and SCCPs were below the LOD in water samples from the background stations, except for BDE-99 and 100 (Table 1). The maximums were analyzed in waters from Kardzhali and Studen Kladenets reservoirs. BDE-47 and -153 were the dominant congeners determined in all water samples. Total PBDEs in water samples did not exceed the EQS, while SCCPs were above the EQS in Zhrebchevo Reservoir (about three times greater). The later confirmed the impact of untreated wastes discharged directly into the reservoir and along the banks.

Moss-Bags
In this study, Cd, Cr, Cu, Ni, and Pb were accumulated by the moss between 17 and 130 times more than their respective background levels (Table 1). Significant contaminants analyzed in moss tissues were Al, As, Cd, Cr, Cu, Ni, and Pb ( Table 2). The calculated CFs indicated certain pollution for As, Al, Cr; strong for Cu and Ni; and certain to extreme for Cd and Pb. As a consequence, the MPIs derived revealed heavy contamination in both Kardzhali and Studen Kladenets reservoirs, which could be related to the former activity of metallurgical plant and industrial impact. The Cu and Pb concentrations measured in the moss-bags in the present study were similar to those reported for river transplants in areas contaminated due to past mining industries [56]. The comparison with transplanted F. antipyretica for 28 days in copper-contaminated Orta Lake in Italy [43] showed similar accumulation for Cr and Pb; however, we determined higher values in the present study for Cr and lower ones for Cu, respectively. The analyzed Pb concentrations in the transplants corresponded to the ambient water levels and confirmed that the accumulation of the elements with no metabolic function, increased with increasing the metal supply in the environment [57].
The transplanted moss CFs for Cd, Cu, and Pb were several times higher than CFs from a study of trace metals from past mining activities with F. antipyretica transplants in rivers in France [58]. Furthermore, the highest CFs in the present study were more than 30 times higher for Cd, 20 times for Cu, and 13 times for Pb, which confirms the high toxic bioaccumulation potential of the studied heavy metals.
Three main groups of elements were revealed by cluster analysis: the first separated Hg with the highest accumulation in the Zhrebchevo Reservoir ( Figure 2). The second cluster included the anthropogenic elements Cu, Cr, Mn, Ni, Pb, and Zn, whose highest concentrations were measured mainly in the Studen Kladenets Reservoir. Manganese was separated from the rest of the elements early (at about 15% similarity). The third cluster grouped elements that could be interpreted as natural elements (Al and Fe), as well as As, Cd, and Co, the highest concentrations of which were analyzed in moss-bags predominantly from the Kardzhali Reservoir. The moss-bags displayed highest overall concentrations of PBDEs (from 0.09 to 0.118 µg kg −1 ). Moreover, the SCCPs levels were alarmingly high in all plant samples (Table 1).

Mussels
The minimum and maximum concentrations for the studied compounds in mussels were presented in Table 1. Overall, the levels were similar between the studied stations. However, the highest values were determined for Al, As, Co, Fe, Mn in Zhrebchevo; Cr, Cu, Ni in Kardzhali; Cd, Pb, Zn in Studen Kladenets Reservoir. On the other hand, the lowest concentrations were measured in the mussels from the background site, except for Zn.
The BAF calculations showed variations between the three reservoirs ( Table 3). The highest values were found for Cu, Mn, and Pb for the mussels from Zhrebchevo; for Al, Cd, Co, Cr, P, and Zn for the mussels from Studen Kladenets, and As, M, Na, and Ni for the mussels from Kardzhali Reservoir. In addition, the BAF for Al, Cr, Cu, Ni, and Pb was extremely high for the mussels from all three studied sites, reaching thousand times more its values in the mussels compared to its corresponding levels in the water. Thus, according to the calculated BAF values, we can categorize the mussels as macroconcentrators (>2). A result, which is worth mentioning is that the BAF for As and Zn differed between Kardzhali, Studen Kladenets, and Zhrebchevo reservoirs (Table 3), which is most likely due to the different type and degree of contamination. The moss-bags displayed highest overall concentrations of PBDEs (from 0.09 to 0.118 µg kg −1 ). Moreover, the SCCPs levels were alarmingly high in all plant samples (Table 1).

Mussels
The minimum and maximum concentrations for the studied compounds in mussels were presented in Table 1. Overall, the levels were similar between the studied stations. However, the highest values were determined for Al, As, Co, Fe, Mn in Zhrebchevo; Cr, Cu, Ni in Kardzhali; Cd, Pb, Zn in Studen Kladenets Reservoir. On the other hand, the lowest concentrations were measured in the mussels from the background site, except for Zn.
The BAF calculations showed variations between the three reservoirs ( Table 3). The highest values were found for Cu, Mn, and Pb for the mussels from Zhrebchevo; for Al, Cd, Co, Cr, P, and Zn for the mussels from Studen Kladenets, and As, M, Na, and Ni for the mussels from Kardzhali Reservoir. In addition, the BAF for Al, Cr, Cu, Ni, and Pb was extremely high for the mussels from all three studied sites, reaching thousand times more its values in the mussels compared to its corresponding levels in the water. Thus, according to the calculated BAF values, we can categorize the mussels as macroconcentrators (>2). A result, which is worth mentioning is that the BAF for As and Zn differed between Kardzhali, Studen Kladenets, and Zhrebchevo reservoirs (Table 3), which is most likely due to the different type and degree of contamination. Table 3. Bioaccumulation factors (BAF) for 11 elements and organic contaminants in mussels from the studied stations. The analyzed elements (Ca, Fe, Hg, K), which were below the limit of detection (LOD) in all studied stations were excluded. Bold-highest BAF value for each element/substance. Our results for bioaccumulation were in line with other studies on trace elements in the mussels from anthropogenic sites. We consider that such studies are essential not only for ecological monitoring but also for health risk assessment of heavy metal exposure via consumption of mussels. For example, our results for Cu and Fe concentrations in the total soft tissue of mussels corresponded to those of Yap et al. [59] (2016); however, their results for Cd, Pb, Ni, and Zn were higher. Furthermore, Azizi et al. [60] analyzed the metal concentrations (Cd, Cr, Cu, Fe, Ni, Zn, Co, and Pb) in the whole soft tissues of mussels from Cala Iris offshore (Northern Morocco) in different seasons. Unlike us, the authors did not detect Co and Pb; the Cu concentrations were lower; the Cd concentrations were similar, but the Cr, Fe, Ni, and Zn concentrations were times higher.

Kardzhali Studen Kladenets Zhrebchevo
The applied cluster analysis divided the trace elements in S. woodiana into three groups. The first cluster, separated early at about 5% similarity, further divides into two sub-groups-one, containing only Fe, which is a naturally occurring trace element, but along with the elements in the second sub-cluster (Zn and Al, Cu), they also could be of anthropogenic origin, mainly due to ore extraction and mining activities. The second cluster, at only about 10% similarity contains one element-Hg, whose toxic properties and behavior are different form all other studied contaminants, and this may explain its separation in a substantive cluster. The third cluster was also further divided into two sub-groups-the first one contained As and Pb, and the second Cd, Ni, Mn, Co, and Cr, which are typical contaminants, probably grouped together, based on synergetic mechanisms in their interactions in the environment (Figure 3). In general, the sum of PBDEs in the mussels from Kardzhali and Studen Kladenets showed a similar trend as in the water-the levels were close, 0.041 and 0.045 µg kg −1 , respectively (Table 1). Furthermore, they were also higher than the measured ∑PBDE concentration in the mussels from the Zhrebchevo Reservoir (0.012 µg kg −1 ). The predominant PBDE congener was BDE-47, accounting for 50% of the total PBDEs in mussels from all three reservoirs. This finding is in line with Debruyn et al. [61], who studied marine mussels (Modiolus modiolus Linnaeus, 1758) from 14 stations near a municipal outfall. A notable result, which should be pointed out, is that the SCCPs (6.12 µg kg −1 ) in the mussels from the Studen Kladenets Reservoir were 27 times higher than those measured in the Zhrebchevo Reservoir and 10 times higher than those in the Kardzhali Reservoir, respectively, which again confirmed the impact of human activity in the studied stations.
The BAF calculations for PBDEs were similar for all three reservoirs (Table 3). They were found to be between 0.46 and 0.59, which, according to Nikanorov [53], determines the mussels as deconcentrators (BAF < 1). However, the BAF calculations for SCCPs varied between reservoirs with highest values of 7.12 for the Studen Kladenets Reservoir. This finding describes S. woodiana as a macroconcentrator (BAF > 2) for this organic contaminant, even though for Kardzhali and Zhrebchevo Reservoir it ranges between 0.05 and 0.47. The calculated BAFs indicated that PBDEs have a pattern of bioaccumulative behavior in the mussels, but SCCPs may be more bioaccumulative in the mussels than PBDEs.
Our results for S. woodiana showed lower values for PBDEs compared to those for Mytilus edulis (Linnaeus, 1758) along the Dutch coast [62], where BAF were measured to be reaching up to 10 for some PBDEs congeners. Our BAF results (except for Studen Kladenets Reservoir) were also lower for SCCPs compared to the calculated BAF for fish and invertebrates from China by Sun et al. [63], which varied from 2.46 to 3.49. They were also lower than the BAF result for SCCPs for bivalves (1.08-1.61) in the study of Ma et al. [64]. We consider that in our case, this result also depends on the tested species.

Integrated Transplant Monitoring
Both bryophytes and mussels are in a close relationship with their immediate environment. The selected bryomonitor, F. antipyretica, has large biomass and high cation exchange capacity. Mussels In general, the sum of PBDEs in the mussels from Kardzhali and Studen Kladenets showed a similar trend as in the water-the levels were close, 0.041 and 0.045 µg kg −1 , respectively (Table 1). Furthermore, they were also higher than the measured PBDE concentration in the mussels from the Zhrebchevo Reservoir (0.012 µg kg −1 ). The predominant PBDE congener was BDE-47, accounting for 50% of the total PBDEs in mussels from all three reservoirs. This finding is in line with Debruyn et al. [61], who studied marine mussels (Modiolus modiolus Linnaeus, 1758) from 14 stations near a municipal outfall. A notable result, which should be pointed out, is that the SCCPs (6.12 µg kg −1 ) in the mussels from the Studen Kladenets Reservoir were 27 times higher than those measured in the Zhrebchevo Reservoir and 10 times higher than those in the Kardzhali Reservoir, respectively, which again confirmed the impact of human activity in the studied stations.
The BAF calculations for PBDEs were similar for all three reservoirs (Table 3). They were found to be between 0.46 and 0.59, which, according to Nikanorov [53], determines the mussels as deconcentrators (BAF < 1). However, the BAF calculations for SCCPs varied between reservoirs with highest values of 7.12 for the Studen Kladenets Reservoir. This finding describes S. woodiana as a macroconcentrator (BAF > 2) for this organic contaminant, even though for Kardzhali and Zhrebchevo Reservoir it ranges between 0.05 and 0.47. The calculated BAFs indicated that PBDEs have a pattern of bioaccumulative behavior in the mussels, but SCCPs may be more bioaccumulative in the mussels than PBDEs.
Our results for S. woodiana showed lower values for PBDEs compared to those for Mytilus edulis (Linnaeus, 1758) along the Dutch coast [62], where BAF were measured to be reaching up to 10 for some PBDEs congeners. Our BAF results (except for Studen Kladenets Reservoir) were also lower for SCCPs compared to the calculated BAF for fish and invertebrates from China by Sun et al. [63], which varied from 2.46 to 3.49. They were also lower than the BAF result for SCCPs for bivalves (1.08-1.61) in the study of Ma et al. [64]. We consider that in our case, this result also depends on the tested species.

Integrated Transplant Monitoring
Both bryophytes and mussels are in a close relationship with their immediate environment. The selected bryomonitor, F. antipyretica, has large biomass and high cation exchange capacity. Mussels as filter feeders were known to reflect levels of organic and inorganic pollutants occurring in the water, food, or sediment and thus also raise public health concerns.
Mosses contain elevated amount of N and P and other essential nutrients, and their uptake is influenced by the trace elements. A loss of K in the transplanted moss (between 2 and >4 times, Table 1), which is a very soluble nutrient, could be a result of heavy metal accumulation as suggested by Vázquez et al. [65]. The trace elements probably influenced also the content of Na and P, most displaced in the Kardzhali Reservoir. The concentrations of mainly intracellular macroelements (e.g., K and P) decreased in moss tissue. In contrast, both moss and mussel tissues were enriched with Ca and Mg at the end of the exposition period in all studied stations, which confirmed the elements as mainly exchangeable ones [66]. Mussels retained similar levels of K and Na during the experiment.
Principal component analysis (PCA) logically divided the water samples from biota along the first axis. Eigenvalues of the first two axes were 0.824 and 0.136, together explaining 95.9% of the samples-pollutants relation. The analysis indicated that in general, mosses and mussels followed a different model of trace element accumulation (Figure 4). F. antipyretica accumulated higher amounts of most of the elements (on the bottom right: Al, As, Cd, Co, Cr, Ni, Pb, Zn), while S. woodiana was associated with elevated concentrations of Mn. as filter feeders were known to reflect levels of organic and inorganic pollutants occurring in the water, food, or sediment and thus also raise public health concerns. Mosses contain elevated amount of N and P and other essential nutrients, and their uptake is influenced by the trace elements. A loss of K in the transplanted moss (between 2 and >4 times, Table  1), which is a very soluble nutrient, could be a result of heavy metal accumulation as suggested by Vázquez et al. [65]. The trace elements probably influenced also the content of Na and P, most displaced in the Kardzhali Reservoir. The concentrations of mainly intracellular macroelements (e.g., K and P) decreased in moss tissue. In contrast, both moss and mussel tissues were enriched with Ca and Mg at the end of the exposition period in all studied stations, which confirmed the elements as mainly exchangeable ones [66]. Mussels retained similar levels of K and Na during the experiment.
Principal component analysis (PCA) logically divided the water samples from biota along the first axis. Eigenvalues of the first two axes were 0.824 and 0.136, together explaining 95.9% of the samples-pollutants relation. The analysis indicated that in general, mosses and mussels followed a different model of trace element accumulation (Figure 4). F. antipyretica accumulated higher amounts of most of the elements (on the bottom right: Al, As, Cd, Co, Cr, Ni, Pb, Zn), while S. woodiana was associated with elevated concentrations of Mn. PBDEs in the biota exceeded the EQS (∑PBDEs in moss from 0.09 to 0.12 µg kg −1 ; mussel 0.01-0.05 µg kg −1 ). There was a predominance of BDE-99, 100, and 153 in the moss tissues, whereas BDE- PBDEs in the biota exceeded the EQS ( PBDEs in moss from 0.09 to 0.12 µg kg −1 ; mussel 0.01-0.05 µg kg −1 ). There was a predominance of BDE-99, 100, and 153 in the moss tissues, whereas BDE-47 dominated in the mussels. BDE-47 dominated both in terrestrial mosses and herbivores in Norway [67], and its dominance in mussels probably reflects different bioaccumulation properties or metabolic transformation processes of the BDE-47 in aquatic environment. BDE-47 was again the most frequently found in Swedish terrestrial mosses with concentrations between <0.06 and 0.46 ng/g DW, while BDE 28, 47, 99, 100, 153, 154 were generally low or below the quantification limit for the analysis [68].
The levels of PBDEs in the moss showed similar pollution in the three reservoirs, while the mussels assessed Kardzhali and Studen Kladenets as more contaminated. There were no statistically significant linear correlations between the PBDEs in the water and in the biota.
As for the SCCPs, the moss samples accumulated high levels, which correlated to those in the water. Short-chain chlorinated paraffins in mussels have a highest calculated BAF of 7.12, thus, this compound could be considered to be bioaccumulative.

Conclusions
The active monitoring applied here, which combines moss and mussel data, seems adequate to study the persistent impact of industrial and untreated waste pollution on aquatic environment. The transplants with Fontinalis antipyretica and Sinanodonta woodiana in this study revealed severe contamination with heavy metals in two of the reservoirs related to former industrial impact, which on the other hand was not detected by water samples. Another general finding was that PBDEs and SCCPs were found in all biotic samples and were confirmed as bioaccumulative compounds.
Although inorganic and organic pollutants accumulated by the moss were higher than those in the mussels, both monitors are complementary. Mosses and mussels followed a different model of trace element accumulation. In addition, the domination of BDE-99, 100, and 153 in the moss and BDE-47 in the mussels reflects different bioaccumulation properties. Thus, integrated monitoring with moss-bags and transplanted mussels should be applied in reservoirs, especially those used for aquaculture, in order to control the water status, to observe pollution trends, and to prevent a potential hazard to human health.
Still, many questions persist: how many monitoring sites are necessary and where to place them. Furthermore, in this pilot study we assume that the uptake of the contaminants is rapid when concentrations are high, but the influence of the negligible water flow velocity in reservoirs has to be studied further. We believe that in the future, such monitoring will be increasingly used to find the answer to the above questions and allow assessment of different pollutants in standing water bodies.