Background
Sediment dredging
The Dredging project was done in January 1998 to remove about 30 cm of sediment from the lake, covering an area of 3.3 km2 with a total volume of 87×104 m3. The sediment formed in the past 40-100 years was removed according to the measured sedimentation rate of 0.3-0.7 mm/yr. On March 7, 1998, the project on the north lake was finished, and the southeast lake was also by the end of April. Pumping water from the Yangtze River replaced the lake water in late October 1998 and was finished by the end of 1998 (Fig. 2).
The grab dredger removed 40-50cm of sediment on December 27, 2008. The preliminary project was carried out in the north lake and completed in March 2009, with 68×104 m3 removed. Subsequently, it was carried out in the southeast and southwest lakes. The dredging work was carried out in the southwest and southeast lakes, with a total water area of 81×104 m2 in 2014; the final dredging volume reached 43×104 m3. The project was implemented in the southeast lake, with a dredging volume of 3×104m3 in 2017, 9.3×104m3 in 2018, 9×104m3 in 2019, and 14.2×104m3 of sediment in 2020.
Water diversion from the Yangtze River
Since the dredging project was completed in 1998, 8×104 m3 of water has been diverted daily to Xuanwu Lake from the Yangtze River. By October 2002, when the second water diversion project was completed, the daily capacity of water diversion had increased to 18×104m3. The cyanobacteria blooms broke out in Xuanwu Lake in July 2005, and the daily capacity increased to 28×104m3 to control the blooms. Four water outlets introduced the lake water into urban rivers (Fig.1-b), and the annual water replenishment volume was about 50 million m3. The entire lake water can be replaced every 16 days. The changes in water diversion and the dredging project are listed in Fig. 2.
Emergency management with clay flocculation and P. crispus population recovery stage
The first recorded cyanobacteria blooms in Xuanwu Lake occurred in July 2005, with Microcystis dominating the blooms. On September 20, September 30, and October 16, 2005, algal blooms in the north, southeast, and southwest lakes were managed using a clay flocculation technique. The last day of treatment was October 31, 2005. About 300 tons of modified clay (kaolin) were applied for the treatment at a rate of 106 g/m2.
The turions of P. crispus began to sprout naturally in the shallow waters before November 2005, after clay flocculation, and the seedlings of P. crispus were scattered sporadically in shallow water. The average population density was 3-5 plants/m2, and the plant height ranged from 3 to 10 cm. Although P. crispus was still limited to shallow water from December 2005 to February 2006, plant height increased to 30-70 cm, with an average density of 30 plants per m2. As the water temperatures rose rapidly in March, P. crispus flourished and expanded to the whole lake. In April 2006, the P. crispus plants were clipped artificially because of overpropagation. The top of the plants was cut off at about 20-30 cm, and the harvest was completed in mid-April. The life cycle of this plant species in the lake includes a sprouting stage at the end of October, a sporadic growth stage from November 2005 to February 2006, a rapid growth stage from March to April 10th, and the P. crispus population was harvested from the middle to the end of April.
Other measures
The amount of non-point source and point source pollution entering Xuanwu Lake decreased concurrently with the closure of the pollution sources surrounding the lake. Additionally, dredging has been done in a few of the lake's input rivers.
Sampling and analysis methods
We collected water samples jointly with the Nanjing Environmental Monitoring Center from 2003 to 2006, so our sampling sites were consistent, and the analysis results showed that our monitoring results were the same. The monitoring indicators include total nitrogen (TN), ammonium nitrogen (NH3-N), total phosphorus (TP), permanganate index (PI), and five-day biochemical oxygen demand (BOD5). Therefore, the annual water quality index from 1998 to 2020 was based on the data from the Nanjing Environmental Monitoring Center and the specific month data from our monitoring.
Five sampling points are shown in Fig. 1-b. Water temperature, dissolved oxygen (DO) concentration, and pH were measured using a multiparameter meter (YSI 6820, Yellow Spring Instruments, USA) at 9 o'clock, and Secchi disc transparency (SD) was measured using a standard Secchi disk. Water samples were collected at a depth of 5 cm below the surface. The water samples were pre-treated (5% potassium sulfate solution, 121 ºC, 30 min). TN, NH3-N, and TP were determined using an automated wet chemistry analyzer (SAN++, Skalar, Holland). Chlorophyll a content (Chla) was determined by the spectrophotometer method using acetone as extraction, and PI was determined by a titrimetric method by acid digestion with potassium permanganate oxidation (SEPA, 2002). The samplings were conducted at intervals of 7 days with three replicates.
The data on heavy metal concentration (1997–2000) was referenced from Zhu's research in 2004 (Zhu et al. 2004), and the data on heavy metals (2006–2020) was collected from the Nanjing Environmental Quality Report (by the Nanjing Environmental Monitoring Center). The data from the same monitoring points and the same monitoring methods are adopted to ensure the consistency of the data.
The Yangtze River water was precipitated by two waterworks and poured into Xuanwu Lake. The water quality parameters of the Yangtze River were used to compare the correlation with the water quality of Xuanwu Lake because the water quality parameters of the two waterworks were lacking. The water quality indicators of the Yangtze River (Nanjing section) were obtained from the Nanjing Statistical Yearbooks, respectively. The indicators include PI, BOD5, NH3-N, and TP, but TP is from 2003 to 2020; other indicators are from 1999 to 2020.
Geoaccumulation index (Igeo)
The geoaccumulation index (Igeo), designed by Müller (1969), has been widely used to study heavy metal pollution in sediments and soils. The geoaccumulation index (Igeo) is calculated by the following equation:
Igeo= log2(Cn/1.5Bn) (1)
where Cn represents the heavy element's measured content (mg·kg-1), and Bn represents the local geochemical background value of the heavy element (mg·kg-1). The constant, 1.5, was the background matrix correction factor due to the lithogenic effect (Masocha et al. 2022; Mao et al. 2022).
The contamination level is classified into seven classes, where: Class 0: Igeo≤0 uncontaminated; Class 1: 0< Igeo≤1 uncontaminated to moderately contaminated; Class 2: 1 < Igeo≤2 moderately contaminated; Class 3: 2 < Igeo≤3 moderately to heavily contaminated; Class 4: 3 < Igeo≤4 heavily contaminated; Class 5: 4< Igeo≤5 heavily to extremely contaminated; and Class 6: Igeo >5 extremely contaminated (Masocha et al. 2022; Mao et al. 2022).
Statistical analyses
The data were subjected to variance analysis using SPSS 16.0. The significance of differences between treatments was determined at the 0.05 probability level (p<0.05).