Research on Water Environment Regulation of Artificial Playground Lake Interconnected Yangtze River

With the rapid development of China, water pollution is still a serious problem despite implementation of control measures. Reasonable water environment management measures are very important for improving water quality and controlling eutrophication. In this study, the coupled models of hydrodynamics, water quality, and eutrophication were used to predict artificial Playground Lake water quality in the Zhenjiang, China. Recommended “unilateral” and “bilateral” river numerical models were constructed to simulate the water quality in the Playground Lake without or with water diversion by pump, sluice and push pump. Under “unilateral” and “bilateral” river layouts, total nitrogen and total phosphorus meet the landscape water requirement through water diversion. Tourist season in spring and summer and its suitable temperature result in heavier eutrophication, while winter is lighter. Under pumping condition, water quality and eutrophication of “unilateral” river is better than “bilateral” rivers. Under sluice diversion, the central landscape lake of “unilateral river” is not smooth, and water quality and eutrophication is inferior to the “bilateral”. When the water level exceeds the flood control level (4.1 m), priority 1 is launched to discharge water from the Playground Lake. During operation of playground, when water level is less than the minimum level (3.3 m), priority 2 is turned on for pumping diversion or sluice diversion to Playground Lake. After opening the Yangtze river diversion channel sluice, priority 3 is launched for sluice diversion to the Playground Lake. When the temperature is less than 15 °C, from 15 °C to 25 °C and higher than 25 °C, the water quality can be maintained for 15 days, 10 days and 7 days, respectively. Corresponding to the conditions of different priority levels, reasonable choices of scheduling measures under different conditions to improve the water quality and control eutrophication of the Playground Lake. This article is relevant for the environmental management of the artificial Playground Lake, and similar lakes elsewhere.

The key pollution areas that may exist in the plant area are the original glyphosate production workshop, the phosphonic acid production workshop, the 10% glyphosate aqueous solution production workshop and the sewage treatment station. Landfill covers an area of about 200,000 m 2 , mainly used for domestic waste and a small amount of construction waste, with emphasis on groundwater pollution monitoring. This place contains toxic and hazardous substances, mainly volatile organic compounds, semi-volatile organic compounds, heavy metals, and total petroleum hydrocarbons. The shipyard covers an area of about 27,000 m 2 , and the plant was mainly for machining, and the production of small vessels. The key toxic and hazardous substances of the shipyard are mainly volatile organic compounds, semi-volatile organic compounds, heavy metals and total petroleum hydrocarbons. Based on the results of the preliminary sampling, a total of 139 soil test sites were found in the chemical plant area, and pollutants  The chemical plant, demolished from 2008 to 2009, covers an area of about 170,000 m 2 . From 2011 to 2015, the construction of an asphalt pavement was carried out, and soil exchange and agricultural planting were carried out in some areas. The key pollution areas that may exist in the plant area are the original glyphosate production workshop, the phosphonic acid production workshop, the 10% glyphosate aqueous solution production workshop and the sewage treatment station. Landfill covers an area of about 200,000 m 2 , mainly used for domestic waste and a small amount of construction waste, with emphasis on groundwater pollution monitoring. This place contains toxic and hazardous substances, mainly volatile organic compounds, semi-volatile organic compounds, heavy metals, and total petroleum hydrocarbons. The shipyard covers an area of about 27,000 m 2 , and the plant was mainly for machining, and the production of small vessels. The key toxic and hazardous substances of the shipyard are mainly volatile organic compounds, semi-volatile organic compounds, heavy metals and total petroleum hydrocarbons. Based on the results of the preliminary sampling, a total of 139 soil test sites were found in the chemical plant area, and pollutants to 2015, the construction of an asphalt pavement was carried out, and soil exchange and agricultural planting were carried out in some areas. The key pollution areas that may exist in the plant area are the original glyphosate production workshop, the phosphonic acid production workshop, the 10% glyphosate aqueous solution production workshop and the sewage treatment station. Landfill covers an area of about 200,000 m 2 , mainly used for domestic waste and a small amount of construction waste, with emphasis on groundwater pollution monitoring. This place contains toxic and hazardous substances, mainly volatile organic compounds, semi-volatile organic compounds, heavy metals, and total petroleum hydrocarbons. The shipyard covers an area of about 27,000 m 2 , and the plant was mainly for machining, and the production of small vessels. The key toxic and hazardous substances of the shipyard are mainly volatile organic compounds, semi-volatile organic compounds, heavy metals and total petroleum hydrocarbons. Based on the results of the preliminary sampling, a total of 139 soil test sites were found in the chemical plant area, and pollutants were detected at all sites. A total of 95 soil contaminants were detected in soil testing, with a detection rate of 72%. There are 38 pollutants that exceeded the screening value. A total of 62 kinds of groundwater pollutants were detected, and 25 kinds of pollutants exceeded the screening value. After the implementation of the playground project, all chemical plants will be cleaned up. Therefore, this study does not consider the impact of the chemical plant on the project. According to the preliminary investigation and analysis of the site environment, the detection indexes of domestic waste landfill soil are generally good. Most of the indicators did not exceed the selection of environmental risk evaluation screening standards (the garbage layer was not analysed), groundwater environmental quality standards greatly exceed the standard of Class III [47,48], subject to a greater degree of pollution. After demolition was completed in 2009, the shipyard was covered by imported earth filling and construction waste and has remained idle till now. The wharf is far away from the wading area of the project and is free from water pollution according to this study. According to the plan, the total excavation earthwork within the scope of soil treatment is 300,000 m 3 , the earthwork to be treated is 220,000 m 3 , and the remaining 80,000 m 3 shall be earthwork without treatment.

Hydrodynamic Two-Dimensional Model
The Cartesian coordinate system of two-dimensional hydrodynamic governing equations is the continuity equations and momentum equations for the integral of the three-dimensional Renault Navier-Stokes equations [2,3] of the incompressible fluid along the direction of water depth, which can be expressed as follows: + ∂ ∂x hT xy + ∂ ∂y hT yy + hv s Q (3) where: t represents time; x, y represent Cartesian coordinates; h represents total water depth; η represents water level; ρ represents water density; u and v represent average water depth; f = 2Ωsinϕ denotes the Coriolis factor (Ω is the angular velocity of the Earth's rotation, ϕ is the geographical latitude); S xx , S xy and S yy are the radiation stress tensors; P a is the atmospheric pressure; Q is the point source emissions; g is the gravitational acceleration: where: ρ 0 represents the relative density of water; (u s , v s ) represents the rate at which the outside world is released into the waterbody. Transverse stress T ij includes viscous resistance, turbulent frictional resistance, and differential advection frictional resistance, which can be calculated using the eddy viscosity equation of the mean vertical velocity:

Water Quality and Eutrophication Two-Dimensional Model
(1) Basic equations of water quality model The water quality equation is based on the mass balance equation. The three-dimensional water quality transport equation contains a lot of uncertain parameters. Under the existing conditions, the verification of the model is difficult. Considering the factors such as data and model calculation workload, the average vertical two-dimensional water quality model is adopted [49,50]. Two-dimensional water quality transport equation is: where: C i is the pollutant concentration; u, v are the flow velocity components in the x and y directions, respectively; E x and E y are the diffusion coefficients in the x and y directions, respectively; K i is the pollutant degradation coefficient; S i is the pollutant sediment release item.
In order to introduce a quantitative relationship between sediment resuspension flux and hydrodynamic conditions in the model and reflect the change of resuspension flux of each pollutant in the sediment with the flow velocity. The sediment resuspension flux is calculated using the relationship obtained from sediment resuspension experiments when establishing the mathematical model [49,50], which mainly reflects the handling of the source sink term Si, as follows: where: α i is the sediment resuspension flux (g/(m 2 ·d)), α i = ζ i ·β i exp(ξ i ·P); H represents water depth (m); β i is the proportion of sediment pollutants in SS (%); P represents co-velocity (cm/s), P = √ u 2 + v 2 ; ζ i , ξ i are the sediment resuspension parameters.
(2) Basic equations of Ecolab eutrophication model The content of Chl-a in lakes is the major parameter for evaluating the water trophic status. There are many factors affecting Chl-a content in lakes. It is generally acknowledged that light, temperature, precipitation, nutrients and pH affect it. In this paper, the impact of nutritive salts (total nitrogen, total phosphorus) concentrations on algae growth were investigated. Chl-a concentration was regarded as the evaluation index. According to the principle of mass conservation, the basic equation of eutrophication variables is [36][37][38][39][40][41]: where: In the equation, C chl-a represents the concentration of Chl-a; U and V respectively represent the flow velocity components in the x and y axis directions, which can be calculated from the water volume model; E x and E y represent the lateral and longitudinal diffusion coefficients of algae, respectively; S chl-a represents the conversion of Chl-a; G PI represents the algal growth; D PI represents the algal death; V s represents the algal sedimentation; D represents the water depth; K 1 refers to the correlation coefficient between Chl-a content and photosynthesis of phytoplankton; Phtsyn refers to the photosynthesis value of plants in unit water volume; µ refers to the mortality rate under optimal nutrition conditions; F (N, P) indicates the nutrient limit function: where PN min and PN max are respectively the minimum and maximum nitrogen content of algae (gN/gC). PP min and PP max are the minimum and maximum values of phosphorus content in algae, respectively (gP/gC). KC is the phytoplankton half-saturation content of phosphorus (gP/gC).

Pollution Source Prediction after the Completion of the Playground
After the chemical plant clean-up completion, landfills are also all clean. The ground is also cleared of pollutants in the Shuijie commercial area, Sea World and Polar Dry Skiing areas. Only travelers will be exposed to some source of pollution in all aspects of tourism. The exposed population in this project mainly comes from two sources: one is the tourist population from other places; the other is the residents nearby. Based on the analysis of market penetration of the same type of tourism products, it is concluded that the capture rate of the project to the residential market is about 15% and that of the tourist market is about 6%. The number of tourists in this project is about 1.7 million/a. According to the prediction formula of tourist quantity, the traffic flow of local residents is about 1.03 million/a. The number of foreign visitors is about 670,000/a and the peak daily passenger flow is 6000/a.
(1) Index construction Emission of pollutants from tourism activities (S): affected by such factors as the per capita water consumption of tourism, sewage discharge coefficient, sewage treatment rate, pollutant concentration and other factors. Emissions of tourism activities can be calculated using the following formula: where S is the discharge of pollutants (such as tourism accommodation activities, tourism catering activities, excursions, etc.); W is the per capita water consumption; r is the sewage discharge coefficient; and u is the sewage treatment rate; C D represents the concentration of pollutants after sewage treatment; C w represents the concentration of pollutants before sewage treatment. Tourist Activities Disturbance Index, TADI: The indexes of total nitrogen (TN), total phosphorus (TP), five days biochemical oxygen demand (BOD 5 ) and chemical oxygen demand (COD) were selected to construct a comprehensive evaluation index of pollutant emissions of tourists. Due to the different evaluation criteria of various types of pollutant emissions, in order to synthesize various types of pollutant emissions evaluation results, each type of pollutant emissions must be non-dimensional treatment. First, the indicators are standardized mean difference, and then use the weighted arithmetic average method to construct tourist comprehensive evaluation index of pollutant emissions. Among them, the range of standardization formula is as follows: where M i is the standard value of pollutants; S i is the pollutant index value; S min is the minimum pollutant index value; S max is the maximum pollutant index value. The TADI is calculated as follows: where M is the standard value of pollutants; i is the type of tourism activity subscript, i = 1, 2, 3 . . . n, that there are n kinds of tourism activities; j is the pollutant type subscript j = 1, 2, 3, 4. TADI, is in the range of [0, 1]. Based on the range of pollutant emission intensity and application, TADI can be divided into five basic evaluation grades of 1, 0.8, 0.6, 0.4, 0.2. Grade 1 value below 0.2, low-interference type; grade 2 value between 0.2 and 0.4, lower interference type; grade 3 value between 0.4 to 0.6, moderate interference type; grade 4 value between 0.6 and 0.8, higher interference type; grade 5 value above 0.8, highly interference type. The calculated concentration of pollutants in the main part of the playground is shown in Table 1.

Model Setup
In the model calculation, the Magical Marine World in Zhenjiang is divided into a three-quadrangle mixed grid with a grid spacing of about 8-10 m. Suppose the initial time the lake is stationary, there is no disturbance, the time step ∆t = 60 s. The original planning "bilateral" river, the proposed "bilateral" river grid model and topographic elevation shown in Figures 3 and 4.
where is the standard value of pollutants; is the type of tourism activity subscript, = 1, 2, 3 … , that there are kinds of tourism activities; is the pollutant type subscript = 1, 2, 3, 4. TADI, is in the range of [0, 1]. Based on the range of pollutant emission intensity and application, TADI can be divided into five basic evaluation grades of 1, 0.8, 0.6, 0.4, 0.2. Grade 1 value below 0.2, low-interference type; grade 2 value between 0.2 and 0.4, lower interference type; grade 3 value between 0.4 to 0.6, moderate interference type; grade 4 value between 0.6 and 0.8, higher interference type; grade 5 value above 0.8, highly interference type. The calculated concentration of pollutants in the main part of the playground is shown in Table 1.

Model Setup
In the model calculation, the Magical Marine World in Zhenjiang is divided into a threequadrangle mixed grid with a grid spacing of about 8-10 m. Suppose the initial time the lake is stationary, there is no disturbance, the time step t = 60 s. The original planning "bilateral" river, the proposed "bilateral" river grid model and topographic elevation shown in Figures 3 and 4.   The wading area elevation as a whole changed within 0-3 m, of which the deepest lake center, the central elevation of 0 m, Moya area at about 1.5 m, Shuijie area at about 2 m. According to the actual topography of the bwading area, the geology and geographical location of the model, the Manning coefficient is 38 m 1/3 /s and the eddy parameter is 0.28. See the Appendix A for initialization parameters.

Parameters Selection
Due to the current planning state and no excavation operations on the playground, the model calculation results were checked according to Wang [35][36][37][38], in situ data and routine monitoring data of Neijiang Lake and diversion channel, in order to meet more the strong relation between the real status of the park and the model to estimate the effect of changes. The model water quality and eutrophication parameters are listed in Table 2. During the high tide of spring-tide of the Yangtze River, the water level of the Yangtze River is higher than that of the Neijiang Lake, with a 0.3 m sluice. Water will be diverted into the Playground Lake by gravity to improve water quality and eutrophication by the use of water head between the #1 inlet, #2 inlet and outlet ( Figure 5).

Parameters Selection
Due to the current planning state and no excavation operations on the playground, the model calculation results were checked according to Wang [35][36][37][38], in situ data and routine monitoring data of Neijiang Lake and diversion channel, in order to meet more the strong relation between the real status of the park and the model to estimate the effect of changes. The model water quality and eutrophication parameters are listed in Table 2. During the high tide of spring-tide of the Yangtze River, the water level of the Yangtze River is higher than that of the Neijiang Lake, with a 0.3 m sluice. Water will be diverted into the Playground Lake by gravity to improve water quality and eutrophication by the use of water head between the #1 inlet, #2 inlet and outlet ( Figure 5). Model boundary conditions: The initial water level was set at 2.67 m. The temperature was 28 °C. The initial flow rate was set to zero [42][43][44]. The initial water head Δh (the water head between the Yangtze River and the Playground Lake) was 0.2-0.3 m. The flow rate of water diversion through sluice at #1 inlet and #2 inlet were 8.27 and 1.82 m 3 /s, respectively.
Weather conditions: The actual measurement data show that the average temperature in 2016 is 16.8 °C; the coldest month is January, the average monthly temperature is 2.97 °C; the hottest month is August, the average monthly temperature is 29.39 °C. The average annual precipitation in recent two years is 1082.7 mm, which is unevenly distributed. The precipitation is mostly concentrated in the spring, summer and autumn seasons [45,46]. Especially, the precipitation is the highest in the summer, exceeding 45% of the total annual precipitation. The average annual wind speed is 3.4 m/s. Initial nutrient salts and pollution sources load: The amount of emissions by tourists and sewage produced by tourism facilities were calculated as pollution sources load, according to the measured data and index construction method [48].
(2) Water diversion through Pump Through the #1 pumping station, the diversion water will enter the playground to improve its water quality and eutrophication status. The position of "three sluices and one culvert" is shown in Figure 6. The specification parameters are given in Table 3. Table 3. Specification parameters of "three sluices and one culvert". The Playground Lake has a storage capacity of about 350,000 m 3 . The design flow of a single pump for the #1 pump is 1.85 m 3 /s. water can be changed once within 27 h by a double pump diversion, according to the initial calculation. The pump diversion will not affect recreational activities such as watercraft due to its low flow, so that continuous pump diversion can be used. In order to ensure that the water is completely replaced, the proposed pumping duration is 30 h. If we turn off the #1 sluice, turn on the #1 pump, #2 sluice, #3 sluice and #4 culvert to divert water, most of the water will flow out from #2 sluice, and the water in the Moya area cannot be effectively changed. Therefore, it is necessary to properly allocate "three sluices and one culvert" to make the water Model boundary conditions: The initial water level was set at 2.67 m. The temperature was 28 • C. The initial flow rate was set to zero [42][43][44]. The initial water head ∆h (the water head between the Yangtze River and the Playground Lake) was 0.2-0.3 m. The flow rate of water diversion through sluice at #1 inlet and #2 inlet were 8.27 and 1.82 m 3 /s, respectively.

Number Size
Weather conditions: The actual measurement data show that the average temperature in 2016 is 16.8 • C; the coldest month is January, the average monthly temperature is 2.97 • C; the hottest month is August, the average monthly temperature is 29.39 • C. The average annual precipitation in recent two years is 1082.7 mm, which is unevenly distributed. The precipitation is mostly concentrated in the spring, summer and autumn seasons [45,46]. Especially, the precipitation is the highest in the summer, exceeding 45% of the total annual precipitation. The average annual wind speed is 3.4 m/s. Initial nutrient salts and pollution sources load: The amount of emissions by tourists and sewage produced by tourism facilities were calculated as pollution sources load, according to the measured data and index construction method [48].
(2) Water diversion through Pump Through the #1 pumping station, the diversion water will enter the playground to improve its water quality and eutrophication status. The position of "three sluices and one culvert" is shown in Figure 6. The specification parameters are given in Table 3. Table 3. Specification parameters of "three sluices and one culvert". The Playground Lake has a storage capacity of about 350,000 m 3 . The design flow of a single pump for the #1 pump is 1.85 m 3 /s. water can be changed once within 27 h by a double pump diversion, according to the initial calculation. The pump diversion will not affect recreational activities such as watercraft due to its low flow, so that continuous pump diversion can be used. In order to ensure that the water is completely replaced, the proposed pumping duration is 30 h. If we turn off the #1 sluice, turn on the #1 pump, #2 sluice, #3 sluice and #4 culvert to divert water, most of the water will flow out from #2 sluice, and the water in the Moya area cannot be effectively changed. Therefore, it is necessary to properly allocate "three sluices and one culvert" to make the water completely replaced. As the ratio of the storage capacity of the Shuijie area and central landscape lake region to the Moya area is about 3:1, the pumping method design is from 0:00 to 22:00 on the first day (lasting 22 h) to change water in the Shuijie area and central landscape lake region, and from 22:00 on the first day to 6:00 on the second day (lasting 8 h) to change water in the Moya area (Table 4). completely replaced. As the ratio of the storage capacity of the Shuijie area and central landscape lake region to the Moya area is about 3:1, the pumping method design is from 0:00 to 22:00 on the first day (lasting 22 h) to change water in the Shuijie area and central landscape lake region, and from 22:00 on the first day to 6:00 on the second day (lasting 8 h) to change water in the Moya area (Table 4).   Table 5.   Table 5.  The model calculation schemes are shown in Table 6. According to the measured meteorological data, the southeaster wind with the highest wind frequency was selected in summer and the average temperature was 28 °C. The flow field, water quality and eutrophication of the wading area under no water diversion, pump diversion and sluice diversion were simulated, respectively. The northeast wind with the highest wind frequency was selected in winter and the average temperature was 7 °C. The flow field, water quality and eutrophication of the wading area under no water diversion and pump diversion were simulated respectively.    The model calculation schemes are shown in Table 6. According to the measured meteorological data, the southeaster wind with the highest wind frequency was selected in summer and the average temperature was 28 • C. The flow field, water quality and eutrophication of the wading area under no water diversion, pump diversion and sluice diversion were simulated, respectively. The northeast wind with the highest wind frequency was selected in winter and the average temperature was 7 • C. The flow field, water quality and eutrophication of the wading area under no water diversion and pump diversion were simulated respectively.  The model calculation schemes are shown in Table 6. According to the measured meteorological data, the southeaster wind with the highest wind frequency was selected in summer and the average temperature was 28 °C. The flow field, water quality and eutrophication of the wading area under no water diversion, pump diversion and sluice diversion were simulated, respectively. The northeast wind with the highest wind frequency was selected in winter and the average temperature was 7 °C. The flow field, water quality and eutrophication of the wading area under no water diversion and pump diversion were simulated respectively.   It can be seen from the figures that the overall flow rate of the entire wading area is small under the condition of no water diversion, with an average velocity of 0.003 m/s. Only the center of the lake is shallow at the edge of the lake and the flow velocity is relatively fast. Pump and sluice diversion of the overall flow field have significantly improved. However, due to the meandering distribution of the waterway in the Shuijie 3 zone and the narrow width of the estuary (12~20 m), the flow field in the Shuijie 3 area is not smooth. Further, a large area of stagnant area will be formed in the Shuijie 3 part connected to the central landscape lake. Coupled with the large flow of people here this caused a greater point source pollution.

Number Size
(2) The new plan of the river layout program In order to improve the flow field in wading area, three kinds of river construction schemes are proposed considering the original river layout and project construction progress: 1 Original Planning river: the original planning river layout shown in Figure 9a; 2 "Roughed" river: roughly handle the Shuijie 2 and the ring lake (planted with plants, retaining walls and other "soft isolation" devices). To improve the flow ratio of Shuijie 2 and 3, so as to allow more water to enter the Shuijie 3. The layout of the river is shown in Figure 9b. 3 "Bilateral" rivers: smoothing the river channel in the Shuijie 3 to improve the flow ratio of the Shuijie 2 and 3 so as to allow the water to pass through the two rivers evenly. The layout of the river is shown in Figure 9c; 4 "Unilateral" river: remove the wading area north of Shuijie 3, so that all water flows into the Shuijie 2, the river layout shown in Figure 9d. It can be seen from the figures that the overall flow rate of the entire wading area is small under the condition of no water diversion, with an average velocity of 0.003 m/s. Only the center of the lake is shallow at the edge of the lake and the flow velocity is relatively fast. Pump and sluice diversion of the overall flow field have significantly improved. However, due to the meandering distribution of the waterway in the Shuijie 3 zone and the narrow width of the estuary (12~20 m), the flow field in the Shuijie 3 area is not smooth. Further, a large area of stagnant area will be formed in the Shuijie 3 part connected to the central landscape lake. Coupled with the large flow of people here this caused a greater point source pollution.
(2) The new plan of the river layout program In order to improve the flow field in wading area, three kinds of river construction schemes are proposed considering the original river layout and project construction progress: ① Original Planning river: the original planning river layout shown in Figure 9a; ② "Roughed" river: roughly handle the Shuijie 2 and the ring lake (planted with plants, retaining walls and other "soft isolation" devices). To improve the flow ratio of Shuijie 2 and 3, so as to allow more water to enter the Shuijie 3. The layout of the river is shown in Figure 9b. ③ "Bilateral" rivers: smoothing the river channel in the Shuijie 3 to improve the flow ratio of the Shuijie 2 and 3 so as to allow the water to pass through the two rivers evenly. The layout of the river is shown in Figure 9c; ④ "Unilateral" river: remove the wading area north of Shuijie 3, so that all water flows into the Shuijie 2, the river layout shown in Figure 9d. "Roughed" river project complex, being costly, it is not recommended so the "bilateral" river, and "unilateral" river layout plans are recommended.

River Seepage Prevention Area Determination
After the project is implemented, the landfill where the river is located is excavated. At the chemical plant location the old soil is replaced with new, non-polluting soil backfill. The remaining landfill scope and bottom elevation at the end of the project are shown in Figure 10a. Playground design ground elevations after infrastructure foundation and river excavation are shown in Figure 10b.
The range and elevation of the remaining landfill overlaying the recommended "bilateral" rivers are obtained by linear interpolation, as shown in Figure 10c. It can be seen from the figure that there is no rubbish in the area where the bottom of the river is located, and only some of the river banks are close to the remaining rubbish. Therefore, there is a risk of pollutants infiltrating into the earth, which may cause water pollution in the wading area. "Roughed" river project complex, being costly, it is not recommended so the "bilateral" river, and "unilateral" river layout plans are recommended.

River Seepage Prevention Area Determination
After the project is implemented, the landfill where the river is located is excavated. At the chemical plant location the old soil is replaced with new, non-polluting soil backfill. The remaining landfill scope and bottom elevation at the end of the project are shown in Figure 10a. Playground design ground elevations after infrastructure foundation and river excavation are shown in Figure 10b.
The range and elevation of the remaining landfill overlaying the recommended "bilateral" rivers are obtained by linear interpolation, as shown in Figure 10c. It can be seen from the figure that there is no rubbish in the area where the bottom of the river is located, and only some of the river banks are close to the remaining rubbish. Therefore, there is a risk of pollutants infiltrating into the earth, which may cause water pollution in the wading area. After the implementation of the project, the chemical plant will be completely cleaned up and After the implementation of the project, the chemical plant will be completely cleaned up and the landfill site under the river will be completely cleaned up. The Shuijie business district, Sea World and Polar Dry Skiing ground contaminants have also been cleared. In order to prevent landfill leachate approaching the slope from entering the river water body, seepage prevention is required when the distance between the remaining garbage and the riverway is small (less than 20 m), and the seepage prevention zone is finally determined, as shown in Figure 11a. River slope protection type is shown in Figure 11b.

Results and Discussion
According to the definition of eutrophication, it is well known that abundant nitrogen, phosphorus and other nutrients, slow flow velocity and suitable climatic temperature are the three essential conditions for the eutrophication. Although the climatic conditions cannot be controlled by man and cannot be described quantitatively, the nutrient content and hydrodynamic conditions in the water body can be controlled. The following two factors are used as variables to predict the eutrophication of water body. Nitrogen and phosphorus nutrients are the most important material basis for algae growth. Chlorophyll (Chl-a) can reflect the status of algae breeding. Therefore, this study proposed the combination of three single indicators of total nitrogen (TN), total phosphorus (TP) and chlorophyll (Chl-a) as the evaluation factors of urban water body landscape water quality.

Analysis of Flow Field without Water Diversion
Using the constructed two-dimensional model of the hydrodynamic, the comparison of flow fields in "unilateral" and "bilateral" rivers under no water diversion conditions is shown in Figure  12. It can be seen in the absence of water diversion, the flow velocity in the two rivers is very slow, with little difference.

Results and Discussion
According to the definition of eutrophication, it is well known that abundant nitrogen, phosphorus and other nutrients, slow flow velocity and suitable climatic temperature are the three essential conditions for the eutrophication. Although the climatic conditions cannot be controlled by man and cannot be described quantitatively, the nutrient content and hydrodynamic conditions in the water body can be controlled. The following two factors are used as variables to predict the eutrophication of water body. Nitrogen and phosphorus nutrients are the most important material basis for algae growth. Chlorophyll (Chl-a) can reflect the status of algae breeding. Therefore, this study proposed the combination of three single indicators of total nitrogen (TN), total phosphorus (TP) and chlorophyll (Chl-a) as the evaluation factors of urban water body landscape water quality.

Analysis of Flow Field without Water Diversion
Using the constructed two-dimensional model of the hydrodynamic, the comparison of flow fields in "unilateral" and "bilateral" rivers under no water diversion conditions is shown in Figure 12. It can be seen in the absence of water diversion, the flow velocity in the two rivers is very slow, with little difference.

Analysis of Flow Field without Water Diversion
Using the constructed two-dimensional model of the hydrodynamic, the comparison of flow fields in "unilateral" and "bilateral" rivers under no water diversion conditions is shown in Figure  12. It can be seen in the absence of water diversion, the flow velocity in the two rivers is very slow, with little difference.

Analysis of Flow Field with Water Diversion through Pump
Using the constructed two-dimensional hydrodynamic model, the comparison of the flow field between the "unilateral" and "bilateral" channels after pumping for 30 h is shown in Figure 13.

Analysis of Flow Field with Water Diversion through Pump
Using the constructed two-dimensional hydrodynamic model, the comparison of the flow field between the "unilateral" and "bilateral" channels after pumping for 30 h is shown in Figure 13. As can be seen from the figure, under the condition of pumping, due to the reduction of "unilateral" river bank capacity, the pump effect is better than that of "bilateral" river channel, especially in the marked red box mark.

Analysis of Flow Field with Water Diversion through Sluice
Using the constructed two-dimensional hydrodynamic model, we can calculate the flow field of "unilateral" and "bilateral" rivers under the sluice diversion conditions, as shown in Figure 14. As can be seen from the figure, due to the removal of the river in the Shuijie 3 area, the water flow at the location indicated by the red box in the "unilateral" channel is not smooth, and forms a large stagnant area, which is not conducive to the smooth flow of water in the wading area. Due to the sudden widening of the water surface in the lake area from the water block to the central area as shown by the yellow box in the "bilateral" rivers, the stagnation of the reciprocating flow occurs, but it does not affect the overall flow field. This shows that the distribution of flow field in "bilateral" rivers is better than "unilateral" river, so we recommend "bilateral" rivers. As can be seen from the figure, under the condition of pumping, due to the reduction of "unilateral" river bank capacity, the pump effect is better than that of "bilateral" river channel, especially in the marked red box mark.

Analysis of Flow Field with Water Diversion through Sluice
Using the constructed two-dimensional hydrodynamic model, we can calculate the flow field of "unilateral" and "bilateral" rivers under the sluice diversion conditions, as shown in Figure 14. As can be seen from the figure, due to the removal of the river in the Shuijie 3 area, the water flow at the location indicated by the red box in the "unilateral" channel is not smooth, and forms a large stagnant area, which is not conducive to the smooth flow of water in the wading area. Due to the sudden widening of the water surface in the lake area from the water block to the central area as shown by the yellow box in the "bilateral" rivers, the stagnation of the reciprocating flow occurs, but it does not affect the overall flow field. This shows that the distribution of flow field in "bilateral" rivers is better than "unilateral" river, so we recommend "bilateral" rivers.

Analysis of Water Quality and Eutrophication without Water Diversion in Summer
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication comprehensive scores without water diversion in summer (Figure 15a). The mean value of TP, TN, Chl-a and eutrophication score of the wading area obtained after the "unilateral" river reached a completely steady state after 10 days under a southeastern wind were 0.107, 1.385 and 0.04 mg/L and 57.1, and almost all areas show mild eutrophication; Spatial distribution of "bilateral" rivers TP, TN, Chl-a and eutrophication comprehensive scores without water diversion in summer (Figure 15b). After southeastern wind exposure for 10 days, the TP, TN, Chl-a and eutrophication comprehensive scores obtained after the "bilateral" rivers reached a completely steady state have mean values of 0.106, 1.383 and 0.039 mg/L and 56.7. In addition to the exchange of water inlet and outlet area #2, almost all areas show mild eutrophication.

Analysis of Water Quality and Eutrophication without Water Diversion in Summer
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication comprehensive scores without water diversion in summer (Figure 15a). The mean value of TP, TN, Chl-a and eutrophication score of the wading area obtained after the "unilateral" river reached a completely steady state after 10 days under a southeastern wind were 0.107, 1.385 and 0.04 mg/L and 57.1, and almost all areas show mild eutrophication; Spatial distribution of "bilateral" rivers TP, TN, Chl-a and eutrophication comprehensive scores without water diversion in summer (Figure 15b). After southeastern wind exposure for 10 days, the TP, TN, Chl-a and eutrophication comprehensive scores obtained after the "bilateral" rivers reached a completely steady state have mean values of 0.106, 1.383 and 0.039 mg/L and 56.7. In addition to the exchange of water inlet and outlet area #2, almost all areas show mild eutrophication.

Analysis of Water Quality and Eutrophication without Water Diversion in Summer
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication comprehensive scores without water diversion in summer (Figure 15a). The mean value of TP, TN, Chl-a and eutrophication score of the wading area obtained after the "unilateral" river reached a completely steady state after 10 days under a southeastern wind were 0.107, 1.385 and 0.04 mg/L and 57.1, and almost all areas show mild eutrophication; Spatial distribution of "bilateral" rivers TP, TN, Chl-a and eutrophication comprehensive scores without water diversion in summer (Figure 15b). After southeastern wind exposure for 10 days, the TP, TN, Chl-a and eutrophication comprehensive scores obtained after the "bilateral" rivers reached a completely steady state have mean values of 0.106, 1.383 and 0.039 mg/L and 56.7. In addition to the exchange of water inlet and outlet area #2, almost all areas show mild eutrophication.

Analysis of Water Quality and Eutrophication with Pumping Diversion in Summer
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores after pumping 22 h in summer, after which the "unilateral" river model reached a completely steady state ( Figure  16a). The mean values of TP, TN, Chl-a and eutrophication in the wading area were 0.087, 1.296 and 0.028 mg/L and 52.9 respectively; Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer, whereby the "unilateral" river model reached a completely steady state (Figure 16b

Analysis of Water Quality and Eutrophication with Pumping Diversion in Summer
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores after pumping 22 h in summer, after which the "unilateral" river model reached a completely steady state (Figure 16a). The mean values of TP, TN, Chl-a and eutrophication in the wading area were 0.087, 1.296 and 0.028 mg/L and 52.9 respectively; Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer, whereby the "unilateral" river model reached a completely steady state (Figure 16b

Analysis of Water Quality and Eutrophication with Pumping Diversion in Summer
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores after pumping 22 h in summer, after which the "unilateral" river model reached a completely steady state ( Figure  16a). The mean values of TP, TN, Chl-a and eutrophication in the wading area were 0.087, 1.296 and 0.028 mg/L and 52.9 respectively; Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer, whereby the "unilateral" river model reached a completely steady state (Figure 16b). The mean values of TP, TN, Chl-a and eutrophication in wading area were 0.05, 1.203 and 0.017 mg/L and 48.5; Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication score after pumping for 22 h in summer (Figure 16c). After the 22 h pumping, the "bilateral" river model reaches a completely steady state. The mean values of TP, TN, Chl-a and eutrophication score in wading area were 0.086, 1.294 and 0.027 mg/L and 52.7, respectively; Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer (Figure 16d) Chl-a(mg/L) Score (d) Figure 16. (a) Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores after pumping 22 h in summer; (b) Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer, whereby the "unilateral" river model reached a completely steady state; (c) Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication score after pumping for 22 h in summer; (d) Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer.

Analysis of Water Quality and Eutrophication with Sluice Diversion
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores by sluice diversion in summer (Figure 17a). After the "unilateral" river model has reached a completely steady Figure 16. (a) Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores after pumping 22 h in summer; (b) Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer, whereby the "unilateral" river model reached a completely steady state; (c) Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication score after pumping for 22 h in summer; (d) Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication score after pumping for 30 h in summer.

Analysis of Water Quality and Eutrophication with Sluice Diversion
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores by sluice diversion in summer (Figure 17a). After the "unilateral" river model has reached a completely steady state under the sluice diversion conditions, the mean values of TP, TN, Chl-a and eutrophication score in wading area were 0.04, 1.202 and 0.016 mg/L and 46.3, respectively. Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication scores by sluice diversion in summer (Figure 17b). After the "bilateral" rivers model has reached a completely steady state under the sluice diversion conditions, the mean values of TP, TN, Chl-a and eutrophication score in wading area were 0.04, 1.191 and 0.012 mg/L and 45.3, respectively.

Analysis of Water Quality and Eutrophication without Water Diversion in Winter
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores without water diversion in winter (Figure 18a). The mean values of TP, TN, Chl-a and eutrophication score of the wading area after the "unilateral" channel model reached the completely steady state were 0.103, 1.36, and 0.026 mg/L and 51.6. Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication comprehensive scores in winter without water diversion (Figure 18b). The mean values of TP, TN, Chl-a and eutrophication score of the wading area after the "bilateral" channel model reached the completely steady state were 0.101, 1.355 and 0.02 mg/L and 50.6, respectively.

Analysis of Water Quality and Eutrophication without Water Diversion in Winter
Spatial distribution of "unilateral" river TP, TN, Chl-a and eutrophication scores without water diversion in winter (Figure 18a). The mean values of TP, TN, Chl-a and eutrophication score of the wading area after the "unilateral" channel model reached the completely steady state were 0.103, 1.36, and 0.026 mg/L and 51.6. Spatial distribution of "bilateral" river TP, TN, Chl-a and eutrophication comprehensive scores in winter without water diversion (Figure 18b). The mean values of TP, TN, Chl-a and eutrophication score of the wading area after the "bilateral" channel model reached the completely steady state were 0.101, 1.355 and 0.02 mg/L and 50.6, respectively.  Based on the TN, TP and Chl-a predicted concentrations calculated by the model, the assessment is made according to the current domestic lake eutrophication score and classification criteria as stipulated in the "Assessment Methods of Surface Water Environmental Quality". The monthly eutrophication scores are given in Table 7. The table shows that throughout the year there was mild eutrophication in June to October, of which the most serious occurred in August.  Based on the TN, TP and Chl-a predicted concentrations calculated by the model, the assessment is made according to the current domestic lake eutrophication score and classification criteria as stipulated in the "Assessment Methods of Surface Water Environmental Quality". The monthly eutrophication scores are given in Table 7. The table shows that throughout the year there was mild eutrophication in June to October, of which the most serious occurred in August.  Based on the TN, TP and Chl-a predicted concentrations calculated by the model, the assessment is made according to the current domestic lake eutrophication score and classification criteria as stipulated in the "Assessment Methods of Surface Water Environmental Quality". The monthly eutrophication scores are given in Table 7. The table shows that throughout the year there was mild eutrophication in June to October, of which the most serious occurred in August.
From the water quality and eutrophication model space and time prediction shows: (1) In the layout of the "unilateral river" and "bilateral rivers" schemes, TP and TN meet the requirements of landscape and recreation with the water diversion, and the overall water quality meets the requirements. Spring and summer correspond to the tourist season and the light temperature is appropriate, the degree of eutrophication heavier, while in winter it is lighter; (2) In the absence of water diversion in summer, due to the slow flow rate in the wading area, the entire wading area under the "unilateral river" and "bilateral rivers" layout is in a mild eutrophication status, in particular in the southeast corner of the Shuijie 1-2 corner, due to dense passenger flow, resulting in heavier human point source pollution. As the winter temperature is low and the light intensity is low and short, is not conducive to algae growth, chlorophyll a concentration is low, and there is less eutrophication. (3) In the condition of pumping, the water quality and eutrophication in the location where #1 push pump is located is better than that of "bilateral" rivers, because "unilateral" river lacks the Shuijie Section 3 and more water flows through Shuijie Sections 1 and 2. (4) In the case of sluice diversion, the water quality and eutrophication in the west side of center landscape lake area the "unilateral" river was inferior to the "bilateral" river because of poor water flow. (5) In overview of the eutrophication occurred region, for the use of three kinds of water diversion: sluice diversion, pumping and push pump, sluice diversion and pumping improved the entire wading area as a whole. For the commercial area of Shuijie where cyanobacteria are prone to occur, as well as the central landscape lake and the southeast corner of the wading area, push pump or algae salvaging are provided for local water quality improvement.

Comprehensive Benefit Analysis
According to the comprehensive benefit analysis of "unilateral" and "bilateral" rivers from four aspects of comprehensive project cost, tourism benefit, environmental benefit and river layout, we can see: (1) From the analysis of environmental benefits: there is not much difference between the two river layouts. (2) From the economic and social benefits analysis: from the short-term benefits point of view, we recommended "unilateral" river. From a long-term benefit point of view, "bilateral" rivers will have more domestic and international market share of the source, so that both social benefits and economic benefits. See Table 8 for details.

Overall Scheduling Scheme
The Neijiang lake normal water level is 3.9 m, the minimum control level when changing water is 3.6 m, the highest control level is 4.1 m, the highest flood control water level is 5.9 m. The internal water system and Neijiang lake water levels are basically the same, and the designed playground water level is between 3.60~4.10 m. In the rainy season, once the water level in the playground area exceeds 4.0 m, the pumping station will start to drain water. When the water level in the river reaches 3.8 m, the pumping station will be stopped. Based on the comparison and selection of programs, the final overall scheduling solution is determined. Flood control and drainage scheduling are implemented to ensure that floods do not occur in the wading area. When the approach channel is opened, the sluice will start water diversion and the pump station will start diversion when it cannot be opened. According to the temperature of the corresponding scheduling days, we control the pump diversion days. When the temperature is higher than 15 • C and no water diversion using push pump to speed up the retention, the flow area slows. The overall schedule is shown in Table 9.
Priority 1: prioritize flood control and drainage schedule to ensure that no disaster occurs in the wading area. The specific drainage program is as follows: (1) storms when the water level over 4.0 m, close three sluices and one culvert, open 2 sets of pumps for drainage. (2) when the water level is lower than 3.9 m, stop one pump (when the water level exceeds 3.9, all will resume). (3) when the water level is lower than 3.8 m stop all pumps.
Priority 2: when the playground water level is low(generally less than 3.3 m), affecting the operation of the water project, shut down the three sluices and one culvert and open the 2 pumps of #1 to divert water to stop at 4.1 m, or by sluice diversion hold the wading area operating water level.
Priority 3: when the tidal level of the Yangtze River is higher than the water level of Neijiang Lake (3.9~4.1 m), the approach channel will open its sluice and divert water to open the three sluices and one culvert to change the playground water.
Priority 4: when there is no sluice diversion use pumping: according to the temperature of the corresponding scheduling days, open the #1 pump diversion. When the temperature is higher than 15 • C and no water diversion when the push pump is used to speed up the retention, slow the flow area. In rainy season, when the rainfall level is higher than 4.0 m.
During the project operation period, the river water level is lower than 3.3.
Open the sluice in the diversion channel Note: Sluice diversion is generally carried out in June, July and August (temperature above 23 • C), twice per month.

Conclusions
The purpose of this article was to regulate the water environment of the artificial Playground Lake-interconnected Yangtze River system. Playground location choice, construction, mathematical modelling, water quality forecasting and other aspects were analyzed, and after a comprehensive benefit analysis a scheduling scheme was proposed. The eutrophication model presented in the paper was used to simulate under different water diversion scenarios. According to the simulation and numerical analysis, the main conclusions are as follows: (1) The soil and groundwater pollutants were analyzed for the site selection and the effects of pollutants in the chemical plant, domestic waste landfill and shipyard were known. According to the excavation scope of the playground, the impermeable range of the artificial lake is obtained, and different types of slope protection (scenery stone revetment, soil slope protection, wood pile slope protection) are confirmed in combination with the landscape. The pollution source during the operation of playground is predicted by the Tourist Activities Disturbance Index method. (2) In order to ensure that the water quality and eutrophication of the original planning river meet the water requirements for landscape entertainment, three water diversion schemes (pumping, sluice, push pump) are proposed to control the eutrophication of the water in the wading area. In order to improve the flow field in the wading area, three new river programs ("roughed" river, "bilateral" river, "unilateral" river) are proposed. Due to the complexity of the "roughed" river engineering, two river layout schemes, "bilateral" and "unilateral", are recommended. (3) In the layout of the "unilateral river" and "bilateral rivers" schemes, TP and TN meet the requirements of landscape and recreation by the water diversion, and the overall water quality meets the requirements. Spring and summer correspond to the tourist season and the light temperature is appropriate, the degree of eutrophication heavier, while in winter it is lighter.
In the absence of water diversion during summer, the entire wading area under the "unilateral" and "bilateral" rivers is under mild eutrophication. In winter, due to the low temperature, light intensity and shorter duration, conditions are not conducive to the growth of algae, the concentration of chlorophyll a is low, and the degree of eutrophication is lighter. In the condition of pumping, the water quality and eutrophication in the location where #1 push pump is located is better than that of "bilateral" rivers. In the case of sluice diversion, the water quality and eutrophication in the west side of center landscape lake area the "unilateral" river was inferior to the "bilateral" river because of poor water flow. In view of the eutrophication area, three types of water transfer methods were used: sluice diversion, pump diversion, and push pump cycle. sluice and pump diversion leading to an overall improvement over the entire wading area, while push pump provided partial improvements. (4) Comparing the comprehensive schemes, we determine the final overall scheduling schemes: Priority 1: prioritize flood control and drainage schedule to ensure that no disaster occurs in the wading area. The specific drainage program is as follows: 1 storms when the water level over 4.0 m, close three sluices and one culvert, open 2 sets of pumps for drainage. 2 when the water level is lower than 3.9 m, stop one pump (when the water level exceeds 3.9, it will resume pumping). 3 when the water level is lower than 3.8 m stop all pumps. Priority 2: when playground water level is low, affecting the operation of the water project (generally less than 3.3 m), shut down the three sluices and one culvert and open the 2 pumps of #1 to divert water to stop at 4.1 m, or by sluice diversion holding the wading area operating water level. Priority 3: when the tidal level of the Yangtze River is higher than the water level of Neijiang Acknowledgments: In this section you can acknowledge any support given which is not covered by the author contribution or funding sections. This may include administrative and technical support, or donations in kind (e.g., materials used for experiments).

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A