Abstract

Endogenous contamination from the newly submerged sediment may have an impact on the water quality of the Danjiangkou Reservoir, the water source of the middle route of the South-to-North Water Diversion Project. In this study, in situ sediment samples were collected at different locations (Guojiashan (DJK1), Kongqueyuan (DJK2), Shijiagou (DJK3), Shiqiao Wharf (DJK4), and Songgang Wharf (DJK5) from the north to the south) of the Danjiangkou Reservoir and the upstream tributary using a columnar sediment sampler and statically cultured in laboratory, and the distribution and release of endogenous N and P from the sediment were investigated based on a double-membrane diffusion model. The results showed that the P contents in the main reservoir and the upstream tributary followed the order of . IP accounted for the largest proportion (47.10-59.70%) of total phosphorus (TP), while LP accounted for only a small proportion (0.10-0.90%) of TP. There is a significant difference in the spatial distribution of P at different sampling points, especially for the OP content (226.90 mg·kg-1) with a coefficient of variation of 26.90%. The contents of different forms of P and NH4+-N decreased from the upstream tributary to the main reservoir, which was mainly attributed to the land use type. In the vertical distribution, the contents of different forms of P in DJK4 and DJK5 decreased with the increase of sediment depth. The contents of NH4+-N and PO43--P in the sediment interstitial water of DJK4 and DJK5 were higher than those in the overlying water. In addition, the contents of NH4+-N and PO43--P were higher in DJK5 than that in DJK4. The static culture experiments showed that N and P were mainly released from the sediment-overlying water interface to the overlying water. The release rate of NH4+-N and PO43--P ranged from 13.08 mg·(m2·d)-1 to 21.39 mg·(m2·d)-1 and from 3.06 mg·(m2·d)-1 to 6.02 mg·(m2·d)-1, and the release amount calculated based on the double-membrane diffusion model was  t·a-1 and  t·a-1, respectively. Thus, endogenous contamination from the newly submerged land is an important factor affecting the water quality of the Danjiangkou Reservoir.

1. Introduction

Reservoir sediment is considered to be an important sink of various nutrients in the overlying water [1] and thus even a small change of reservoir sediment can have a substantial impact on the quality of the overlying water [2]. The nutrients accumulated in the sediment such as nitrogen (N) and phosphorus (P) can be released from the sediment back to the overlying water through the interstitial water as a result of convection, diffusion, and resuspension, causing secondary pollution of the overlying water [3]. For instance, the release rate of NH4+-N from the sediment is still as high as  t·a-1 in the Dianchi Lake of Yunnan Province, China, after dredging and sewage interception [4], and resuspension of nutrient-rich sediment causes eutrophication for several decades even when external nutrient loads are curtailed [5]. In this case, endogenous release of nutrients from the sediment is the main factor affecting the water quality [6]. Therefore, a better understanding of the distribution of nutrients in the sediment and their release fluxes is essential to take steps to reduce endogenous nutrient load and improve reservoir water quality [7].

The Danjiangkou Reservoir is the main water source of the middle route of the South-to-North Water Transfer Project of China, and thus the water quality has been maintained at a high level for years. However, evidence has emerged that there is an increasing risk of endogenous contamination that can seriously affect the water quality of the reservoir [8]. The overflow weir crest of the Danjiangkou Reservoir was heightened in 2012, resulting in an increase in the storage capacity from 17.45 billion m3 to 29.05 billion m3 [9] and consequently inundation of 17333 ha of farmland [10]. The newly formed sediments contain a higher amount of organic matter compared with the original sediments, which implies a higher risk of endogenous release and contamination. N is more likely to be released from the newly submerged farmland and orchard of the Danjiangkou Reservoir [11], and P loss is more prominent in the new riparian zone compared with that in the existing riparian zone [12]. Thus, endogenous release of nutrients such as N and P from the sediment in the newly submerged area may pose a threat to the water quality of the reservoir. A number of studies have been conducted to investigate the factors affecting the release of nutrients from the sediment, such as existing states and seasonal variations of nutrients, dissolved oxygen, and temperature [13], but little is known about the release of nutrients from the sediment in the newly submerged area.

A number of models have been proposed to describe the release of nutrients from the lake sediment, such as the permeation model, the double-membrane diffusion model, and the surface renewal model. Among them, the double-membrane diffusion model is particularly suitable to characterize the flux of endogenous nutrients from the kinetic perspective, as it uses mathematical conversion to obtain the double film thickness that is difficult to measure by other models and generalizes the actual complex release process. However, there are three theoretical premises: both water and sediment phases are single homogeneous phases, transmission takes place by means of molecular diffusion, and the total amount of sediment nutrients remains unchanged. It avoids the influence of complicated environmental factors and allows for the estimation of the release of sediment nutrients in a large reservoir area. The value is calculated from the data obtained from the simulated release experiment, which is the net release coefficient derived from the actual sediment adsorption, precipitation, and microbial action [14].

In this study, in situ sediment samples were collected at different locations of the Danjiangkou Reservoir and statically cultured in laboratory, and the distribution and release of endogenous N and P from sediments were investigated based on the double-membrane diffusion model. The results of this study may provide important insights into how to improve the water quality of the reservoir.

2. Materials and Methods

2.1. Study Area and Sampling Sites

The Danjiangkou Reservoir is located in the northwest of Hubei province and the southwest of Henan province and consists of Hanjiang and Danjiang reservoirs. The annual average inflow is about 39.48 billion m3 from the Hanjiang River and its tributary Danjiang River, and the terrain is high in the northwest but low in the southeast. The reservoir was put into operation in 1973, and the overflow weir crest was heightened in 2012. Now, it is the largest artificial freshwater lake in Asia, with a water area of about 102275 ha and a storage capacity of about 29.05 billion m3 [15].

In this study, five representative sampling sites were established at Guojiashan (DJK1; 33°0239N and 111°2846E), Kongqueyuan (DJK2; 33°051N and 111°2711E), Shijiagou (DJK3; 32°5442N and 111°3047E), Shiqiao Wharf (DJK4; 32°4923N and 111°3428E), and Songgang Wharf (DJK5; 32°4602N and 111°3850E), as shown in Figure 1. These sampling sites were about 0.5 km away from the reservoir shore with a water depth of about 2 m. The overlying water and sediment samples were collected in March and August, 2018. The longitude and latitude of the sampling sites were determined by GPS (see Table 1). AD, DJK 4, and DJK 5 were continuous farmland since 2012 where a large amount (188-356 kg·ha-1) of chemical fertilizers (mainly potash and phosphate fertilizers) were used with a maximum of 646 kg·ha-1 [16]. Thus, the double-membrane diffusion model was used to investigate the effect of the release of endogenous N and P on the overlying water at these two sampling sites based on the data collected in March, 2018.

2.2. Sample Collection and Analysis

Columnar sediment samples were collected using a columnar sampler with an organic glass tube (), and three parallel samples were collected at each sampling site for static culture experiments. The overlying water samples were also collected at each sampling site without disturbing the sediment. A portion of water samples was encapsulated in 5 L plastic bottles to supplement the water in static culture experiments, while the rest was encapsulated in 500 ml polyethylene bottles for the measurement of N and P concentrations in the overlying water, and several drops of concentrated sulfuric acid were added to reach a pH value of 2 and shaken uniformly in order to inhibit the microbial activity. Statically cultured samples were cut according to the depth (0-1 cm, 1-2 cm, 2-3 cm, 3-4 cm, 4-5 cm, 5-6 cm, 6-7 cm, 7-8 cm, 8-9 cm, and 9-10 cm) and centrifuged at 5000 r/min for 20 min. Sediment interstitial water was obtained and filtered using a 0.45 μm fiber filtration membrane for the measurement of NH4+-N and PO43--P concentrations in the sediment interstitial water.

2.3. Static Culture Experiments

The overlying water was removed from columnar sediment samples by siphonage, and then the original overlying water was filtered using the 0.45 μm fiber filtration membrane and carefully introduced into the columnar sediment samples along the wall by siphonage. In order to better simulate the field environment of the Danjiangkou Reservoir (the annual average temperature is about 16°C), all columnar sediment samples were placed vertically in a culture vessel with a water depth of 0.40 m in dark at 16°C for 10 days. Water samples were collected at the middle of the sample at the beginning of the experiment and then every 8 h during the culture. Prior to each sampling, the overlying water was stirred slowly using a slender glass rod in order to ensure uniform sampling and to avoid sediment resuspension. After sampling, an equal amount of original overlying water was supplemented. All samples were kept in a constant temperature incubator at 4°C.

2.4. Analysis Methods

The contents of NH4+-N in the overlying water, sediment, and sediment interstitial water were measured using Nessler’s reagent colorimetric method (GB7479-87, National Standards of People’s Republic of China), KCl extraction-distillation method [17], and Nessler’s reagent colorimetric method [17], respectively, while those of PO43--P were measured using molybdenum-antimony anticolorimetry (GB11893-89, National Standards of People’s Republic of China) and the molybdenum blue coloration method [17], respectively. The states of P were determined as described in using the Standards, Measurements, and Testing (SMT) method (Standards Measurements and Testing Program of the European Commission) [18] (see Figure 2). Labile phosphorus (LP) was extracted using 1 mol/L NH4Cl [19]: 1 g of sediment was placed in a centrifuge tube, and 25 ml of NH4Cl (1 mol/L) was added, oscillated for 0.5 h, and centrifuged. The LP content in the supernatant was measured using molybdenum blue colorimetry. All measurements were repeated three times, and the average was reported (the error was within <5%). The significance of spatial difference was analyzed by one-way ANOVA using SPSS 19.0, and figures were plotted using Origin 9.0 and ArcGIS 10.2. Error bars represented the standard deviations (SD) throughout the paper.

2.5. Calculation of Experimental Results
2.5.1. Double-Membrane Diffusion Model

A number of models have been proposed to describe the release of contaminants from the lake sediment, such as the permeation model, the double-membrane diffusion model, and the surface renewal model. In this study, the double-membrane diffusion model is used, in which the molecular diffusion coefficient and the ratio of the thickness of double membranes () are used instead of the thickness of double membranes that is difficult to determine [20]. Three assumptions are made in this model: (1) both the overlying water and sediment are single homogeneous phases, and the overlying water is well mixed, (2) the mass transfer between the sediment and the overlying water takes place only through molecular diffusion, and (3) the nutrient content is relatively high in the sediment compared with that in the overlying water, and thus nutrient loss due to sediment release is negligible throughout the study. The release processes of nutrients from the sediment into the overlying water consist of (1) migration of nutrients from the interior to the surface of the sediment, (2) interfacial transfer of nutrients to the water side, (3) diffusion through the liquid membrane, and (4) uniform mixing in the overlying water. Given the above three assumptions, steps 1, 2, and 4 are assumed to occur instantaneously, and thus step 2 is the controlling step of the release processes [21].

According to the double-membrane diffusion model, the transfer flux should be proportional to the molecular diffusion coefficient and the concentration gradient in the double membrane, the latter of which is the concentration difference between the two sides of the membrane divided by the thickness of the double membrane : where C is the contaminant concentration on the water side (mg·L-1), and is the equilibrium concentration of contaminants in the water phase (mg·L-1).

Let the transfer velocity constant be . According to the double-membrane diffusion model and the dynamics theory,

Thus, the dynamics of contaminant concentration in the water phase can be described as follows: where is the release area (m2), and is the volume of the overlying water (m3).

Eq. (3) can be transformed into where is the depth of the overlying water in the experimental column [22].

The plot of against time is a straight line through the origin of coordinates, and the slope is the value [21]. The values calculated from experimental results are shown in Table 2.

According to the double-membrane diffusion model, the release flux of contaminants from the sediment is

The release rate of contaminants is where is the sediment area (m2).

Thus, the annual diffusion amount is where is the volume of the overlying water (m3), , is the depth of the overlying water (m), and is the diffusion period [21].

2.5.2. Release Rate of Nutrients in Statically Cultured Sediments

The release rate of contaminants can be calculated according to the mass conservation principle: where is the release rate of contaminants from sediment samples, mg/(m2·d); is the volume of the overlying water (L);, , and are the contents of contaminants collected at the beginning of the experiment, time and , respectively (mg/L); is the contaminant content in the added water (mg·L-1), is the volume of samples (L); is the contact area between the sediment and the overlying water (m2); and is the release time () [22].

3. Results and Discussion

3.1. Distribution of Different Forms of P and N in the Sediment

There is a significant difference in total phosphorus (TP) contents (402.40-560.30 mg·kg-1) among the five sampling sites (DJK 1-5) with a coefficient of variation of 14.4% (), which is consistent with the previous finding in the Danjiangkou Reservoir (530 mg·kg-1) [15] but lower than that in the Three Gorges Reservoir (1189 mg·kg-1) [23]. In sampling sites from DJK1 to DJK5, a great amount of P in chemical fertilizers and crop residues can be fixed in the soil by adsorption, sedimentation, and microbial immobilization [24]. However, the TP contents in DJK 1-5 seem to be less affected by agricultural practices compared with that in the Three Gorges Reservoir. The highest TP content is found in DJK 1, which is 1.34 times that in DJK 4. In general, the TP content decreases from the upstream tributary to the main reservoir, which may be associated with the original land use type [25] and the amount of chemical fertilizers used (see Table 1).

The contents of different forms of P follow the order of . Thus, IP accounts for the largest proportion (47.10-59.70%) of TP in the five sampling sites, while the LP content is only 0.98-3.50 mg·kg-1, accounting for a small proportion (0.10-0.90%) of TP (see Figure 3). However, it is noted that LP is a weakly bound phosphate on the sediment surface that is readily exploitable by aquatic organisms, and P is most easily released from the sediment in the form of LP [26]. The vertical distribution of various forms of P in DJK 4 shows a decreasing trend with sediment depth and the contents of Fe-P, which can also be easily utilized by aquatic organisms, range from 1.70 mg·kg-1 to 4.70 mg·kg-1, while that in DJK 5 shows a W-shaped decreasing trend and the Ca-P contents range from 9.30 to 28.50 mg·kg-1. The vertical distribution of P may be associated with the environment, sedimentation rate, and age. The greater the sediment depth is, the lower the contents of various forms of P will be [27]. Fertilizers containing P were not allowed in early years, thus resulting in less accumulation of nutrients in the sediment; whereas an increasing amount of P fertilizers are used in recent years, thus resulting in an increase in P contents in the sediment.

Inorganic N exists in the sediment and overlying water mainly in the form of NH4+-N [28]. Table 3 shows that the NH4+-N content ranges from 173.00 mg·kg-1 to 208.10 mg·kg-1 in the five sampling sites, which is slightly higher than that obtained in a previous study (142.30 mg·kg-1) [29]. In line with the study of the Yuqiao Reservoir [30], the NH4+-N content also decreases from the upstream tributary to the main reservoir. Compared with the reservoir, the tributary river is characterized by smaller bed slopes and flow velocities, making it easier to accumulate organic matter [31]. The spatial distribution of P is also affected by agricultural practices at the five sampling sites due to extensive use of chemical fertilizers. In the vertical distribution, the NH4+-N content increases as the sediment depth increases (see Figure 3), because NH4+-N can be easily diffused from the sediment surface to the overlying water through complex physical, chemical, and biological exchanges [32]. However, the NH4+-N content is increased at 4-12 cm. Organic N is mineralized at an increasing rate in the deep sediment under the effect of ammonifiers, and more importantly, less oxygen is available as the sediment depth increases, resulting in inhibition of nitrification and thus transformation of less NH4+-N into other forms of N. As a consequence, an increase in the NH4+-N content is observed at 4-12 cm [33].

The contents of NH4+-N and PO43--P in the sediment interstitial water of DJK 4 and DJK 5 increase with increasing depth (see Figure 4), and they are significantly higher in DJK4 than in DJK5 (). The average NH4+-N content in the interstitial water at a depth of 10 cm in samples from DJK4 and DJK5 is 16.00 mg·L-1 and 3.50 mg·L-1 with a coefficient of variation of 42.50% and 107.70%, respectively. Thus, the NH4+-N content varies more dramatically in DJK5 [34]. In DJK5, the NH4+-N content increases as the depth increases until a maximum of 12.75 mg·L-1is reached at 3 cm, after which it decreases with further increase of the depth. This agrees well with the study of Liu et al. about the Beili Lake of the Hangzhou West Lake [35] and the study of Wang et al. about the Nansi Lake [36]. The significant difference in the NH4+-N content at 4-10 cm between DJK 4 and DJK 5 can be attributed to the presence of more iron and manganese oxides in DJK5, which act as electron acceptors of NH4+ for conversion of NH4+ into NH3 gas under anaerobic conditions [37]. It is also noted that the higher NH4+-N content in the interstitial water of DJK 4 is consistent with the NH4+-N content in the sediment.

The average PO43--P content in the interstitial water at a depth of 10 cm in samples from DJK4 and DJK5 is 5.26 mg·L-1 and 4.30 mg·L-1 with a coefficient of variation of 20.50% and 36.20%, respectively, which is related to the high P content in DJK 4. In DJK 5, the PO43--P content first increases and then decreases with increasing depth, and the maximum of 8.16 mg·L-1 is reached at 2 cm; while in DJK 4, the PO43--P content is maximized at 4 cm (7.45 mg·L-1) and 10 cm (5.47 mg·L-1). This is because the PO43--P in the surface interstitial water can be easily diffused into the overlying water, resulting in a lower PO43--P content compared with that in the deep interstitial water [38]. In both DJK 4 and DJK 5, the PO43--P content is decreased at 4-7 cm, at which there is a transition from aerobic to anaerobic conditions that can inhibit the effect of aerobic microbes on the active organic detritus and increase the dissolution of iron and manganese oxides under anaerobic conditions, and consequently result in the release of adsorbed P.

3.2. Release of NH4+-N from Statically Cultured Sediments

NH4+-N is released from the sediment-overlying water interface into the overlying water in March in DJK1-5 (see Figure 5), which is in good agreement with that of the Danjiangkou Reservoir in August [11]. The release rate ranges from 13.08 mg·(m2·d)-1 to 21.39 mg·(m2·d)-1, and the maximum is observed in DJK3 (21.39 mg·(m2·d)-1), which is 1.32 times of the average. It is higher than that in the Meiliang Bay of the Tai Lake at 15°C (12.60 mg·(m2·d)-1) [39], but lower than that in the Fubao Bay of the Dianci Lake at 14-16°C (93.03 mg·(m2·d)-1) [40], because organic matter is released in a larger amount from the sediment to the overlying water in DJK3. In the spatial distribution, the coefficient of variation of the NH4+-N content is 50.60% in DJK 4 and 74.30% in DJK 5, both of which are higher than that in DJK1-DJK3 (18.70%, 20.90%, and 20.10%, respectively) which were predominantly used for agricultural purposes.

The flow velocity is generally low at the bottom of the reservoir, making it difficult for nutrients to migrate to the overlying water through the sediment-interstitial water-overlying water interface. Under oxygen-rich conditions, the organic N in the sediment can be converted into inorganic ions such as NH4+ and NO3- by nitrification, which in turn can be diffused to the interstitial water and then to the overlying water and, as a result, have an effect on the NH4+-N content in the overlying water [37]. It is also noted that more NH4+ would be released under anaerobic conditions compared with that under aerobic conditions [41]. In addition, the decomposition of organic N results in the formation of NH3 under hypoxia conditions but NH4+-N under acidic conditions. The NH4+ in the overlying water can also be diffused back to the sediment. The activities of microorganisms in the sediment and subsequent ammonification are reduced in an anoxic environment, and the reducing environment can also promote the denitrification, and finally, N is released into the atmosphere.

3.3. Release of PO43--P from Statically Cultured Sediments

The PO43--P concentration in the overlying water of five sampling points (DJK1-5) tends to increase over time in March, and the coefficient of variation in the reservoir (35.40% and 64.70% in DJK 4 and DJK 5, respectively) is significantly higher than that in the upstream tributary (21.40%, 15.80%, and 26.70% in DJK1-DJK3, respectively) (see Figure 6). The same phenomenon is also observed for the NH4+-N content in the overlying water. The release rate of PO43--P ranges from 3.06 mg·(m2·d)-1 to 6.02 mg·(m2·d)-1 with an average of 4.65 mg·(m2·d)-1, which is lower than that in the Meiliang Bay of the Tai Lake in July (76.55 mg·(m2·d)-1) [42], but slightly higher than that in the Wuli Lake of the Tai Lake (2.05 mg·(m2·d)-1) [43]. The highest release rate is found in DJK1, which is 1.29 times the average. The release rate of PO43--P decreases from the upstream tributary to the main reservoir, which is consistent with the contamination extent of the sediment. However, the release potential of P from the sediment is not entirely determined by the TP amount in the sediment. The release rate of P is high in water with severe organic contamination, which is related to the presence of a large amount of easily degradable matter [44]. At the beginning of the experiment when there is rich oxygen, the active organic debris in the sediment can release a large amount of PO43- under the effect of microorganisms. PO43- can be diffused to the interstitial water and then to the overlying water driven by the concentration gradient, forming the first release peak of PO43-. On the other hand, an anaerobic environment would be formed under the effect of aerobic microbes, and iron and manganese oxides can be easily reduced to increase the solubility, resulting in the release of adsorbed P and thus the formation of the second release peak of PO43-.

3.4. Annual Release Amount of N and P from the Sediment and Risk Evaluation

The increase of dam height in 2012 results in a total reservoir area of 1023 km2, a storage capacity of 29.05 billion m3 and an average water depth of about 30 m. After interception of contaminants, the annual release amount of PO43--P and NH4+-N from the newly submerged sediment into the overlying water is calculated to be 0 t and  t, respectively, using formula (7) based on the migration coefficient () obtained in experiments and the double-membrane diffusion model at an average temperature of 16°C for the Danjiangkou Reservoir. This differs substantially from the release amount of TP based on the SWAT model ( t) [45], probably because some management practices such as moving the topsoil to cultivate fertility can reduce the endogenous contaminants in the sediment and consequently the release amount of TP. The amounts of NH4+-N and PO43--P would be increased from 0.44 mg·L-1 and 0.02 mg·L-1 to 0.57 mg·L-1 and 0.02 mg·L-1 according to the average annual inflow (39.48 billion m3) and to 0.62 mg·L-1 and 0.03 mg·L-1 according to the total reservoir storage (29.05 billion m3), respectively. In theory, the water quality is between mesotrophy and eutrophication (0.03-0.05 mg·L-1 for TP and 0.31-0.65 mg·L-1 for TN) [46]. The reservoir storage is increased in the wet season from May to October, and the increase of temperature leads to an increase in dissolved carbon and microbial activities in the soil. As a result, a large amount of N and P accumulated in the soil would migrate to the overlying water, which is now the main factor affecting the water quality of the Danjiangkou Reservoir. In addition, the newly submerged area is often in alternative submergence-dry states, and land is cultivated for agricultural production under dry conditions [16], thus resulting in an increase in the diffusion of P and N from the sediment into the overlying water.

4. Conclusions

(1)The TP contents in the 0-10 cm sediment samples collected in the Danjiangkou Reservoir and its upstream tributary range from 402.40 mg·kg-1 to 560.30 mg·kg-1, and the contents of different forms of P follow the order of . The spatial difference in P contents among the five sampling points may be related to the original land use type(2)The average PO43--P and NH4+-N contents in the interstitial water of DJK 4 and DJK 5 are higher than those in the overlying water, and thus PO43--P and NH4+-N are likely to be released from the sediment to the overlying water(3)Both experimental and calculation results show that N and P are easily released from the sediment to the overlying water, and thus endogenous contamination from the newly submerged sediment is an important factor affecting the water quality of the Danjiangkou Reservoir

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This study was supported jointly by the National Key Research and Development Project of China (No. 2019YFC0409202), the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (No. 51721006), High Level Talent Support Program of North China University of Water Resources and Electric Power, and Special Support for Innovative Scientific and Technological Team of Water Ecological Security in the water source area of the middle route of South-to-North Water Diversion in Henan Province.