The qualitative and quantitative analysis of the coupled C, N, P and Si retention in complex of water reservoirs

The Solina–Myczkowce complex of reservoirs (SMCR) accounts about 15 % of the water storage in Poland. On the base of historical (2004–2006 years) data, the mass balance of nitrogen, phosphorus, total organic carbon and dissolved silicon were calculated. Large, natural affluents were the main source of the biogenic compounds in the studied ecosystem, delivering 90 % of TOC, 87 % of TN and 81 % of TP and DSi load. Moreover, results show that SMCR is an important sink for all the analysed biogenic elements. About 15–30 % of external loads were retained in the reservoir mainly in upper Solina. Due to the intensive processes of primary production, inorganic forms of nitrogen and phosphorus were mainly retained. Internal production of organic matter lead to an amount of the organic matter deposited in the sediments greater than was anticipated on the basis of the mass balance calculations. A constant load of dissolved silicon originating only from natural sources did not contribute to supplement deficits of Si present in the body of water in the reservoirs, promoting disturbances in N:C:P:Si ratios and another growth condition for other types of algae.

reservoirs and knowledge of these element cycles is essential for our understanding of the ecosystem biogeochemistry (Bouwman et al. 2013). Additionally, this information is important for interpretation of greenhouse gasses (e.g. N 2 O, CO 2 and CH 4 ) emission from reservoirs causes and intensity (Gruca-Rokosz and Tomaszek 2015).
Retention of biogenic elements in water ecosystems is frequently described as a function of the morphometric and hydrologic parameters of reservoirs. The most frequently considered factor in model studies is a hydraulic retention time (HRT); also depth and area of the reservoir and loads (Behrend and Opitz 2000;Seitzinger et al. 2002;Tomaszek and Koszelnik 2003;Hejzlar et al. 2009) are utilised. In general, HRT changes proportionally, as longer stoppage of waters in slow flow zones promotes two basic mechanisms of retention, i.e. uptake by organisms and sedimentation. The behaviour of nitrogen is slightly different, as its retention may be high also in shallow and flowing basins due to conditions promoting denitrification, which in the literature is referred to as a third nitrogen retention mechanism, understood as a difference between inflow and outflow of load (Seitzinger et al. 2002;Tomaszek and Koszelnik 2003).
The purpose of the study is interpretation of how a mountain complex of reservoirs can modify natural biogeochemical fluxes of four major biogenic elements in the river waterway. The Solina-Myczkowce complex of mountain reservoirs (SMCR) is a perfect training field for this interpretation. During the 2004-2006 period reservoir chemistry, hydrology and catchment area management were studied intensely and independently by different institutions.

Study site
The SMCR is located in the upstream part of the River San (SE Poland) within the River Vistula system (Fig. 1). The upper, Solina reservoir is the biggest man-made lake in the Vistula basin and accounts for about 15 % of the total water storage capacity in Poland (volume: 502 mln m 3 , mean depth: 22 m, hydraulic retention time: 215 days, mean discharge: 35 m 3 s −1 ). The lower Myczkowce Reservoir, as a compensatory water body (volume: 10 mln m 3 , mean depth: 5 m, hydraulic retention time: 6 days), is supplied by the hypolimnetic waters of the upper one (90 % of total supply), and by minor tributaries. The upper reservoir has three major inflows (accounting for 90 % of total water supply) and three minor ones. Reversible pumping takes place sporadically. The outflow from the complex involves bottom water from the lower reservoir flowing through a hydroelectric-power plant.
The greater part (c. 75 %) of the 1250 km 2 catchment area is covered by forest, followed by meadows and pastures. Arable land accounts for only a small fraction of the area. The drainage basin has a low population density of about 6 inhabitants per km 2 , while about half of the households concerned are not connected to either sewerage or septic systems. A relatively steep (6 %) slope favours the leaching of soil and ground cover, especially during periods of intensive atmospheric precipitation and snow melting.

Sampling strategy and methods
In order to assess the mass balance of N, P, DSi and TOC, water was sampled from the mouth parts of rivers and streams feeding the Solina reservoir, as well as from the outflow. Also water samples obtained from four stations on the Solina reservoir and two on the Myczkowce reservoir were subjected to analysis. Water was sampled 33 times every 1-6 weeks. Sampling dates were adjusted on the basis of the meteorological and hydrological factors. About 0.5 dm 3 of glass-fiber-membrane-filtered water was made subject to spectrophotometric determinations for concentrations of nitratenitrogen (N-NO 3 − , salicylate method, coefficient of variation for the procedure-CVP: ±1.5 %), nitrite-nitrogen (N-NO 2 − , Griess reaction, CVP: ±1.7 %), ammonium-nitrogen (N-NH 4 + , Berthelot's reaction, CVP: ±1.4 %), phosphate-phosphorus (P-PO 4 3− , Molybdate method, CVP: ±1.8 %) dissolved silicon (the Molybdate method, CVP: ±1.6 %) and chlorophyll a (only in stagnant water, in the methanol) using a PhotoLab S12 spectrophotometer (WTW GmBH). Moreover, Kjeldahl nitrogen (N Kjeld , distillation), Total Organic Carbon (TOC, Shimadzu TOC analyzer, CVP: ±0.6 %) as well as total phosphorus (TP, oxidation to phosphate), were analyzed in non-filtered samples. Total nitrogen (TN) was calculated as the sum of nitrate-and nitrite-nitrogen and N Kjeld . Balances of respective elements, both for the entire complex and the two reservoirs, were calculated as follows: where (mass time −1 ): L R -load inflowing from the drainage basin in affluent waters; L DD -load from the direct drainage basin; L At -load supplied from atmosphere along (1) L R + L DD + L At = L Out + R(E) Fig. 1 The bathymetric map of the Solina-Myczkowce complex of reservoirs. Location of sampling stations from stagnant water (M-Myczkowce, S-Solina) as well as studied affluents was shown with precipitation; L Out -load removed along with run-off; R(E)-retention (elimination) of an element in the ecosystem.
Load inflowing from the drainage basin in affluent waters (L R ) constituted a sum of loads inflowing with all the rivers and streams feeding the balanced ecosystems. Loads of elements for the particular sections were calculated as a product of concentration and water flow rate. Daily flow rate (Q) values of the analysed sections, necessary for the calculation of loads from the three main affluents (about 90 % of total supply) and run-offs from the reservoirs were obtained from Solina-Myczkowce Power Plant S.A (continuous stage measurements. Q value for smaller watercourses was calculated on the day of water sampling, using installed staff gauge readings. Concentrations between the sampling days were calculated using the statistical approach, according to Mukhopadhyay and Smith (2000). Uncertainty of calculated loads was approximated on the base of Harmel et al. (2006) as cumulative uncertainties of potential sources of error. We assumed errors during sampling (uncertainty of 15 %), chemical procedures and analyses (above mentioned CVP, ca. 2 %) and flows measurement (continuous stage measurements, uncertainty of 2 %). Therefore, C, N, P and Si loads are calculation with errors of ±19 %. To estimate the cumulative probable uncertainty of calculated retentions the root mean square error propagation method were used (Harmel et al. 2006). However, the analysis and discussion of data are based on the most probabilistic, average values due to clarity of the paper and per analogiam to many other authors (e.g. Mengis et al. 1997;Garnier et al. 1999;Torres et al. 2007;Hejzlar et al. 2009).
L DD was calculated with as a sum of (1) load inflowing from point sources within the direct drainage basin (including the WWTP); (2) load inflowing from nonpoint sources within the direct drainage basin; (3) load introduced by bathers. Respective summands were calculated from the available models (Giercuszkiewicz-Bajtlik 1990; Jørgensen 2011) including data on the direct drainage basin area development, touristic burden and amount of wastewater discharged from the WWTP obtained from the Solina Municipality with its seat in Polańczyk. L At was obtained from parallel study (Urbanik 2007, MSc thesis, unpublished data) utilising precipitation rate measurements carried out by the hydrological survey station in Lesko. As the reservoirs are not located in the area affected by ground waters, this source of feed was disregarded. In addition, because of the mesotrophic nature of waters, the effect of atmospheric nitrogen binding by cyanobacteria was recognised as insignificant for the balance of this element (Ferber et al. 2004;Koszelnik et al. 2007).
The sedimentation rate of TN, TOC was calculated on the basis of the phosphorus mass balance (which, in contrast to nitrogen and carbon, is not present in a gaseous state in the biogeochemical cycle), using the following formulae (Dudel and Kohl 1992): where: N sed (C sed )-nitrogen (carbon) sedimentation rate (t year −1 ); P ret -phosphorus retention in the reservoir between the sampling time points (t year −1 ); N:P-ratio of total N:P concentration in benthic deposits (from Koszelnik 2009a, b); C:N-ratio of organic C:total N concentration in benthic deposits (from Koszelnik 2009a, b); ∆P-change of (2) N s ed = P r et−∆P(N : P) (3) C s ed = N s ed(C : N ) mean phosphorus content in the water body between the sampling time days (from Koszelnik 2009a, b; t year −1 ); calculated from the following formula: where P n+1 and P n correspond to TP concentration at nth and n + 1th day of sampling (g P m −3 ), V R is a reservoir volume (m 3 ), ∆t is a time period between nth and n + 1th day of sampling (day), and 10 6 and 365 are conversion factors to obtain the t year −1 unit;

Results
The results of the mass balance for the dammed reservoir complex Solina-Myczkowce are listed in Table 1. The complied data show that the inflow of TN to the reservoir from all the sources in 2004 and 2005 amounted to approx. 1700 t. In the following, last year of study, the value was ca. 50 % greater, and amounted to 2489 t. Soluble forms, in particular nitrate(V) nitrogen (ca. 60 %) were the main contributors to the load. For the other contributors loads calculated for each year of balance were less variable. Annual inflow of TP was within a range of 76-101 t; RSi 2056-2244 t; and TOC 2517-2837 t. Predominantly the loads fed the Solina reservoir. Natural affluents of the Myczkowce reservoir played a minor part in biogenic compound feed (3-7 %). Large, natural affluents ( Fig. 1) were the main source of the biogenic compounds in the studied ecosystem, delivering 90 % of TOC, 87 % of TN and 81 % of TP and DSi load (Fig. 2). The share of inflows from the direct drainage basin is more evident in the case of DSi balance (5 %), whereas the load originating from atmosphere is immaterial in the annual balance.
Except for TOC (R 2 = 0.85-0.90; p < 0.001), no significant relations between hydraulic flows and concentration of total forms of the analysed elements in water were discovered. No seasonal changes in N and P concentration indicate that both N and P originate from point and non-point sources. Occurrence of these changes in the case of regionaloriginating Si can be related to its assimilation in river waters upstream of the reservoirs (Humborg et al. 2000).
Calculated values of loads, retention/elimination rates of N, P, Si and C are listed in Table 1. An apparent discrepancy between results from both the reservoirs and the total balance, results from inclusion of reverse pumping of water from the Myczkowce to Solina reservoir in partial balances, which influences particular balances, but has no effect on the balance for the reservoir complex. Significant masses of the balanced elements were found retained within the analysed ecosystem. This indicates that the majority of biogenic compounds feeding the reservoirs is incorporated into the trophic chain and/or accumulated in the benthic deposits (Behrend and Opitz 2000;Koszelnik et al. 2007) by various retention mechanisms. The overall reservoir balance reveals that only for nitrogen was the outflowing load higher than the inflowing load in 2006. In general, retention of biogenic compounds in the entire cascade depended on the amount of elements retained in the upper reservoir. The lower reservoir, due to a short HRT and unfavourable thermal conditions normally has neither retained nor eliminated significant amounts of biogenic compounds. Only in 2005 the silicon retention was predominantly present in the Myczkowce reservoir, with the similar relation true for N-NO 3 − and TOC for the entire study period.  The determined TP loads and percentages of element retention in the reservoir complex during the balanced years were slightly less diversified than nitrogen. 19-33 % of the inflowing TP was retained in the reservoirs. Despite the fact that phosphates constitute slightly more than 50 % of the load inflowing to the reservoirs, the calculated balance reveals, that retention of the easily assimilative phosphorus form was prevailing.
Determined retention of the dissolved silicon (DSi ret% ) in the reservoir complex varied from 5 % (2006) to 24 % (2005) of the annual external load. It was observed that in the first and the last year of the study the majority of the supplied load was accumulated in the Solina reservoir. On the contrary, in 2005 from 492 t of the retained silicon, as much as 352 t was retained in the Myczkowce reservoir. Retention of TOC was observed in both the reservoirs and ranged from 11 to 22 % of the supplied load. In a way similar to the totals of remaining elements, a major amount of carbon was retained in the studied reservoir complex in 2005.
Upon analysis of the morphometric factors on retention, a significant influence of HRT on retention of N, P and Si in waters of the entire reservoir complex was observed (Table 2). Flow rate of inflowing waters was correlated only with N ret% in the Myczkowce reservoir, while the load of elements supplied to the reservoir influenced the N retention in the Myczkowce reservoir and DSi ret% in the Solina reservoir, as well as in the entire complex. Much better correlations were observed upon analysis of the influence of hydraulic retention on the elements retention (W ret -water inflow reduced by water outflow). Increase in W ret led to significant rise of N ret% in both the reservoirs and the complex, DSi ret% and OWO ret% in the Solina reservoir and the complex and P ret% , but only in the balance of the entire reservoir complex (Table 1). Significant correlations between retention values of N, P and Si for the Solina reservoir and the entire complex  (Table 2). Such relations were not seen for the Myczkowce reservoir, there was no influence of TOC ret% on N ret% , and P ret% was correlated only with DSi ret% .
Average annual sedimentation rate of TN, TOC in the Solina reservoir, calculated from the mass balance and levels of individual elements in sediments is provided in Table 3. Sedimentation in the Myczkowce reservoir was negligible due to the very short HRT in this reservoir (approximately 2 days). Production of sediments in this reservoir is mainly related to macrophyte production and supply of external matter. The TOC sedimentation calculated for the entire Solina reservoir amounts to 1300 t per year, and TN 110 t per year. Thickness analysis of the matter deposited from the beginning of reservoir's existence shows that approximately twice as much of the sediments were formed in backwaters (S1, S2 stations) of the reservoir when compared to lacustrine deep areas (S3, S4 stations). This information was taken into consideration upon calculation-it was anticipated that 2/3 of the total sedimentation takes place in the said area.

Discussion
The hydrologic balance (Urbanik 2007-unpublished data) shows that within the 3 years of study the hydraulic inflow to the reservoir was compensated by the outflow, so the calculated element retention is not the result of hydrologic factors, but biogeochemical  factors. Mass balance of N, C, P and Si, conducted separately for both the reservoirs and the entire complex enabled identification of the classical biogeochemical cycle of conversion of the analysed elements. The results of the mass balance show that the SMCR retains a significant amount of biogenic elements (Table 1). The major part of elements was retained mostly in the Solina reservoir. The biogenic compounds were retained sporadically in the Myczkowce reservoir due to the hydrologic factors, i.e. feeding of a hypolimnion with N, C, P and Si-reach waters. The relationships presented in Table 2 show that retention of biogenic elements within the reservoir complex (mostly the Solina reservoir) results not only from the hydrologic or functional factors of the power station (storage of waters during the spring period and discharging during low water periods in summer and autumn), but also from inclusion of easily assimilative forms to the trophic chain and various chemical transformations. The fact that correlation between TN, TP and DSi, which are retained mostly in the easily assimilative forms (Table 1), is significant, it confirms that mechanisms of assimilation by water organisms are crucial for retention of elements in the studied reservoirs. The distinctive feature of the mass balance was DSi depletion from the water body in the reservoirs, mostly in the Solina reservoir (Fig. 2). Approximately 20 % of the inflowing load of DSi was retained in the studied reservoirs. Humborg et al. (2006) reports that the DSi load flowing off from 1 km 2 of the River Vistula basin to the Baltic Sea amounts to 0.8 t per year. Average annual load of DSi feeding the SMCR amounts to 1947 t, which is equivalent to 1.5 t of flow off from 1 km 2 of the basin. Hence, anticipating that silicon originates only from sources connected with soil erosion (Garnier et al. 1999;Humborg et al. 2002) and that the DSi load level produced in other areas of the Vistula basin is similar, ca. 50 % of dissolved silicon is retained within the area of the Vistula basin, unfavourably reducing loads feeding the Baltic sea. This decrease leads to deterioration of seawater quality due to the deficiency of DSi when compared to other biogenic compounds from anthropogenic sources. The said phenomenon leads to imbalance between diatoms and other algae (Humborg et al. 2000). A similar decrease in DSi load in dammed reservoirs was described by Garnier et al. (1999), while Humborg et al. (2006) conclude that the cascade design of dammed reservoirs on rivers increases HRT and favours retention of dissolved silicon, depleting it from downstream waters. 20 % of DSi retention is a distinctive feature of oligotrophic waters (Garnier et al. 1999), but the data available in the literature (Garnier et al. 1999;Humborg et al. 2002Humborg et al. , 2006 also show that similar levels of DSi retention were seen in both oligo-and eutrophic ecosystems. In the studied case, depletion of dissolved silicon in the surface lake body, related in turn to DSi ret affects the water quality, stimulating growth of algae. The correlations shown on Fig. 3 may confirm that increase of chlorophyll level can be related with emergence of non-diatomic (green) algae. Diatoms are seen in lake and reservoir waters mostly in spring, but also even in late winter (Humborg et al. 2000;Lehmann et al. 2004). By analogy to the condition of Lake Lugano (Lehmann et al. 2004) it can be concluded that DSi level exceeding 0.7 g m −3 in the epilimnion of the Solina reservoir contributes to chl a level related with the presence of diatoms and green algae. Below this level, in summer, a rapid increase in chl a, reaching even as much as 12 mg m −3 , is observed which, in turn, may lead to occurrence of thermophilic cyanobacteria with concomitant disappearance of diatoms. Despite low water temperature in the Myczkowce reservoir, an elevated level of chl a was seen, but the index of >2.5 mg m −3 was noted only in 2005, while in 2006 it was low. In general, phytoplankton production in this reservoir is minor. However, due to the poor silicon feed, silicon shortages can occur, which lead to minor tides of algae, mostly in the warmer water area of the dam, analogous to those present in the Solina reservoir (Koszelnik 2013). A decrease in water DSi below 1 g m −3 was seen in 2005, but not in 2006. In addition, when silicon was almost completely depleted from the water body, a decrease in Si:N and Si:P ratios (see Fig. 4) was seen, and DSi became the limiting element. With silicon shortage present, phosphorus and nitrogen are the main substrates utilised in production of organic matter in the reservoirs. In turn, the value of N:P molar ratio (Fig. 4) significantly exceeding 16:1 proves the stoichiometric excess of nitrogen versus phosphorus.

Abbreviations as in the text
Approximately 28 % of phosphorus supplied to the studied reservoirs is retained. The process is occurring mostly in the Solina reservoir. This result depends on various factors allowing storage of phosphorus from the water body in the benthic deposits, related in general with an affinity to specific metals, presence of aerobic conditions or with pH (Golterman 1998). Phosphorus retention is a result of sedimentation of solid particles introduced to the reservoir with affluents and assimilative forms incorporated to the biomass of phytoplankton and transferred to sediments (Hejzlar et al. 2009;Dunalska et al. 2013). A part of such retained phosphorus can be released again to the water body due to resuspension or decomposition of bonds with iron or other metals in anaerobic conditions. The benthic deposits of the Solina reservoir are rich in metals with affinity to phosphorus, mostly in iron. Retention volume of these deposits is very high (Bartoszek et al. 2009), and favourable aerobic conditions make release of phosphorus from the deposit practically absent in both the reservoirs. On that basis it can be concluded that the principal phosphorus retention mechanism in the reservoir is a direct sedimentation of phosphorus contained in suspension and intermediate sedimentation of mineral forms, after their assimilation into the trophic chain. Storage properties of the benthic deposits are large enough so that the calculated P ret value could be greater, but the morphometry of the Solina reservoir (high depth, low area of an active bed) determines the identified level. Consequently, upon balancing of average phosphorus mass retained in the reservoir it was anticipated that the phosphate phosphorus retention will be equal to the amount of this element utilised by phytoplankton, and the difference between retention of P-PO 4 3− and TP will be equal to the amount of element originating from external sources and accumulated in the deposit: A complexity of biogeochemical nitrogen changes in the water environment affects the retention level of this element. In contrast to Si and P, nitrogen has a larger gaseous phase, and denitrification leading to change in state of aggregation affects the mass balance of this element in water ecosystems. In the studied period of 3 years, decrement of load inflowing to the reservoirs amounted only to 9 % per year. This value was affected by a significant element elimination in 2006. In 2005 N ret% was equal to ca. 22 % of the introduced load. Studies on nitrogen retention, carried out between 1999 and 2003, reveal that the retention of this element was varying significantly, amounting to from 12 to 36 % of the annual load (Tomaszek and Koszelnik 2003). Previous studies show that the rate of denitrification in the Solina reservoir is stable and amounts to ca. 5 g m −2 year −1 (Koszelnik et al. 2007), which corresponds to ca. 20 % of retention and 5 % of nitrogen load. Empirical models of denitrification in conditions occurring in the Solina reservoir, contingent on presence of nitrates and temperature (Gruca-Rokosz 2005, PhD thesis, unpublished data), were utilised to analyse various nitrogen retention mechanism. Estimated denitrification rate was 4 g N m −2 year −1 , and on that basis it was calculated that, on average, 70 t of N is denitrified annually. Nitrogen sedimentation calculated from the balance mass equals to 110 t year −1 , hence the balance of retained nitrogen is as follows: Share of denitrification process in nitrogen retention amounted to 40 % and was twice as high as that calculated for previous years (Koszelnik et al. 2007). However, load reduction is significantly lesser (similar to N ret% ) and was only 3.6 %.
Influence of denitrification on the nitrogen mass balance depends on various factors. Seitzinger et al. (2002) describes that for North American estuaries 50 % of nitrogen is supplied from denitrification. Estuaries are bodies of water similar in many cases to dammed reservoirs; mostly due to the ratio between areas of basin and water table, retention time or biogenic compounds load. The said value can be real, but in many cases-including estuaries-significantly lower values are seen, i.e. 5-30 % (Dudel and Kohl 1992;Koszelnik et al. 2007;Povilaitis et al. 2012), but also values as high as 70 % are reported (Mengis et al. 1997). A comprehensive analysis of the available data presented in previous papers (Koszelnik et al. 2007) enables one to conclude that in water regions with high hydraulic dynamics contribution of nitrification to the mass balance is minor when compared to natural lakes, where this process can be significant.
Nitrogen sedimentation, in a way similar to phosphorus, is a result of its consumption. Thus, calculation of N sed and N cons contributions can be difficult. Jickells et al. (2000) states that in inland waters, retention mainly consists of storage of organic nitrogen produced within the ecosystems in the benthic deposits. In the studied case retention of assimilative forms, mostly nitrates, is equal to ca. 100 t year −1 . Nevertheless, it should not be anticipated that the overall mass of retained NO 3 − will be assimilated. Some part of it will be denitrified, as the main substrate for the said process, occurring in the anoxic layer of the benthic deposits, are nitrates(V) diffusing from water . A surplus nitrogen seen in the above balance (−6 t) can result from utilisation of nitrogen stored in the deposits in the nitrification process, which, after various transformations, can be denitrified. The mass balance of total organic carbon calculated for the reservoir complex shows that annually approximately 442 t of TOC is retained. The calculated TOC sedimentation is three times higher and amounts to 1300 t year −1 . Hence, it should be recognised that a significant part of sedimentation matter is produced within the ecosystems: Lentic waters are characterised by a high carbon retention capability, as the major part of TOC supplied to lakes and reservoirs is respirated and included in the trophic chain (Garnier et al. 1999). The above is true for both deep and shallow reservoirs. Anderson and Sobek (2006) provide an example of a shallow lake in Sweden, in which the annual carbon load amounts to approximately 3 t. Calculated phytoplankton production for the lake is as high as 53 t C per year, and the macrophyte production-16 t C per year, while carbon sedimentation is three times greater than the carbon inflow to the ecosystem. Unlike the other elements, TOC retention in the Myczkowce reservoir was fairly high (7-8 % of the load), which can be explained by carbon uptake by macrophytes after its respiration. Rate of decomposition for these forms is low, and annually they release only 40 % of the retained organic carbon (Gessner 2001).

Conclusions
Both the SMCR reservoirs are loaded with nitrogen and phosphorus in amounts significantly exceeding theoretical values considered to be allowable. Although the easily assimilative inorganic forms are dominating in the supplied mass of elements, concentration of both of the biogenic compounds and the amount of chlorophyll a fall within the level specific for mesotrophy. However, the inflow of biogenic elements is so significant that no distinct, seasonal variations of nitrogen and phosphorus concentrations were seen in the reservoirs. The major part of the biogenic compound load supplied to the studied reservoirs is retained therein via utilisation as substrates in the primary production process.
During the summer, dissolved silicon deficits were observed in waters in both the reservoirs. This phenomenon was present due to silicon consumption within the water body and reduced inflow from the basin. In this case silicon became a limiting element for the production of diatomic organic matter, especially in the warmer Solina reservoir. This effect was accompanied by an increase in chlorophyll a concentration, sporadically reaching the value specific for eutrophy, which can be related to production of other species of (non-diatomic) algae.
Retention values of particular elements, calculated from the mass balance, prove the intensity of element uptake process carried out by organisms. Approximately 20 % of inflowing total forms of N, P and Si are accumulated, mostly in the Solina reservoir. Dissolved inorganic forms are retained in more than 50 % of cases. Sedimentation of autochthonous biogenic forms is a result of inclusion of supplied elements to the trophic chain. The denitrification rate of nitrogen amounts to 20 % of retention and only to 5 % of the supplied load. TOC retention at the level of 30 % proves the allochthonous matter is accumulated in the deposit.
Sedimentation of allochthonous organic matter calculated from the mass balance is three times lower than the value estimated on the basis of the overall TOC sedimentation in the Solina reservoir. The remaining amount results from an intra-ecosystemic production stimulated by an external inflow of biogenic compounds and elements circulation within the reservoirs.