Optimization Model of the Ecological Water Replenishment Scheme for Boluo Lake National Nature Reserve Based on Interval Two-Stage Stochastic Programming

: In the present paper, an optimization model of the ecological water replenishment quantity of Boluo Lake wetland was constructed for ecological water replenishment in the Boluo Lake National Nature Reserve based on the interval two-stage stochastic programming (ITSP) method, which scientiﬁcally allocates the ecological water resources and enhances the utilization of ﬂood resources in order to meet the minimum ecological water replenishment quantity requirement to meet the local economic and social water demand, to restore the wetland function of the reserve, and to improve the ecological environment. In addition, it considers the ecological service value of the reserve in order to achieve a shared increase in the ecological and economic beneﬁts. The optimization model of the ecological water replenishment of Boluo Lake wetland considered the minimum ecological water replenishment as the objective function, while the water diversion capacity, water supply capacity, water diversion sequence, functional area, and ecological service value of the Lake bubble were utilized as constraints in order to restore the ecosystem function of the Boluo Lake National Nature Reserve. The results from the model revealed that the amount of ecological water replenishment was signiﬁcantly reduced after optimization simulation, with the maximum reduction range in the ecological water replenishment being − 100.00% and − 74.58%. In addition, the total amount of ﬂood diversion was signiﬁcantly increased, and the ﬂood resources could be fully utilized. Moreover, the recovery effect on the lake and pond functional areas was signiﬁcant and compared much better to that prior to the ecological water replenishment, which was up to 2300.00% and 1987.59%. The ecological service value also increased signiﬁcantly, and the rate of this increase was as high as 23.90% and 21.58%. In the present study, an optimization model was constructed for the ecological water replenishment of the Boluo Lake wetland supplement project based on the interval two-stage stochastic programming method, which would achieve the entire scope of ecological and economic beneﬁts of the ecological water replenishment project by realizing the ecological system reconstruction and providing a feasible and reliable plan to the decision-makers.


Introduction
Wetlands are rich in biodiversity and are considered to be one of the most productive and economically valuable ecosystems in the world [1][2][3] because they provide the ecosystem products and services necessary for human activities, such as fish breeding, plant-carbon sequestration, biodiversity protection, cultural and entertainment services, etc. [4][5][6][7][8]. Despite their importance, over half of the wetland ecosystems in North America, Europe, Australia, and China have disappeared in the early 20th century due to the production activities conducted by humans [9,10]. Currently, there is a severe decline in the wetland ecosystems. The Boluo Lake National Nature Reserve, which is the only large natural wetland in the central part of the Jilin Province, belongs to the "Inland Wetland and Water Ecosystem" type of natural ecosystem. In recent years, this wetland ecosystem has been critically damaged due to the arid climate and over-exploitation by humans.
The amount of available ecological water is a key factor determining the function of the wetland ecosystem [11]. Therefore, the ecological supplement of the wetland is one of the most direct and effective methods for the restoration of a wetland area and improves the wetland ecosystem [12,13]. The main objectives of ecological recharge are to restore the hydrological and soil conditions of degraded wetlands, restore vegetation and faunal habitats, increase species composition, restore biodiversity and ultimately achieve the sustainable development of wetland ecosystems. The main objectives of ecological recharge are to restore hydrological and soil conditions in degraded wetlands, restore vegetation and fauna, increase species composition and restore biodiversity, and ultimately to achieve sustainable development of wetland ecosystems, which are beneficial to the local economy. Ecological water replenishment is applied widely in wetland restoration with good outcomes. Maleki et al. [14] employed a multicriteria-spatial decision support system (MC-SDSS) to design an ecological water replenishment program for the restoration of the Harmon wetland and to reduce the negative impact of drought on this wetland area. Xu et al. [15] constructed a random water allocation optimization model based on the watershed scale to conduct the ecological supplement of wetlands in the Xiaoqing River Basin and to restore the wetland activities in this region. Yang et al. [16] adopted an interval fuzzy double-sided chance-constrained stochastic programming method for designing the wetland ecological water replenishment project in the western region of the Jilin Province, which in addition to improving the ecological water supply, improved the local economic benefits as well. While the ecological water replenishment method is significantly effective in restoring the wetland ecosystem, the excessive allocation of ecological water replenishment may lead to water shortages in the reservoirs and rivers. Therefore, a scientific allocation of the ecological water replenishment is crucial in ecological water replenishment projects.
The interval two-stage stochastic programming method (ITSP) is a novel mathematical programming method that combines the interval programming method [17] with the twostage stochastic programming method (TSP) [18,19], thereby reflecting the uncertainty of parameters in both interval and probability distribution forms. The ITSP method is applied widely to resolve the uncertainty problem and has achieved good results in various research fields. Maqsood et al. [20] employed the ITSP method in the uncertain conditions under the solid waste management system and formulated the waste flow model with the smallest system cost and the largest system feasibility. Meng et al. [21] employed the ITSP method to allocate water resources from the perspective of the pollution capacity of the Yinma River basin, and the wastewater carrying improvement project was adopted, which promoted the improvement in the quality of the water environment to a certain extent. Qiu et al. [22] constructed an optimal allocation model of flood resources based on the ITSP method and incorporated the corresponding flood resource effectiveness under the different discharge scenarios of the Nenjiang River into the total water resources, which improved the flood utilization rate.
In summary, the present study involved the construction of an optimization model based on the ITSP method for the ecological water replenishment of the Baluo Lake wetland supplement project. By utilizing the minimum ecological water replenishment amount of the Boluo Lake National Nature Reserve as the objective function, the optimal allocation of the ecological water replenishment quantity of each lake bubble was performed to ensure the complete utilization of the flood resources, to increase the amount of flood resources, and to effectively restore the lake bubble and the wetland area in reserve. Meanwhile, the ecological service value of the protected area was considered in order to realize the improvement in both ecological and economic benefits.

Study Area
The Boluo Lake National Nature Reserve is located in the northwest region of Nong'an County, Changchun City, Jilin Province. Figure 1 is the map of the Boluo Lake National Nature Reserve. It spans an area of 24,915 hectares and is the only large natural wetland in the central region of the Jilin Province [23]. It is located at the western edge of the depression in the southeast region of the Songliao Plain, with the terrain higher in the southeast and lower in the northwest. Boluo Lake is an inland closed-flow alkaline freshwater lake [24,25]. The Boluo Lake wetland has remained dry for a long time, the rivers in the reserve have been cut-off, and the groundwater level has dropped, which has resulted in a decrease in the lake water, the phenomenon of drying up, a serious reduction in the wetland area, and aggravation of soil erosion. Owing to the unique terrain of this region, the circulation in the water resources is relatively difficult, while soil desertification and land salinization are aggravated, which has caused a serious impact on the local ecosystem security and gradual degradation of the wetland ecosystem function.

Construction of the Ecological Water Replenishment Optimization Model for the Boluo Lake Wetland Based on the ITSP Method
The objective of the present study was to realize the restoration of wetland function and the reconstruction of the wetland ecosystem in the Boluo Lake wetland water replenishment project using minimum ecological water supply through the scientific allocation of the ecological water resources. Therefore, in the present study, the minimum ecological water replenishment was considered to be the objective, while the water diversion capacity, water replenishment quantity of the lake bubbles, and water supply capacity were utilized as constraints. Finally, an ecological water replenishment optimization model, based on the ITSP method, was constructed for the Boluo Lake wetland water replenishment project. The model may be expressed as follows: where f represents the total quantity of ecological water replenishment in the Boluo Lake National Nature Reserve; i = 1 to 5 represent the Toudaogang Reservoir, Boluo Lake, Mobopao, Yuanbaowapao, and Aobaotupao, respectively, in the Nature Reserve; j = 1 to 4 represent the four different ecosystem types, namely, fish ponds, crab ponds, reed wetlands, and marsh wetlands, respectively; h = 1 to 3 represent the flood amount levels in the low flow years, normal flow years, and high flow years, respectively; p h denotes the scenario probability; EAT ± i represents the area of water replenishment recommended by the project (hm 2 ); QT ± i represents the total water replenishment amount for the lake and pond (t/hm 2 ); DEA ± ijh represents the area adjustment amount for the lake and pond i in the ecosystem j under scenario h (hm 2 ); and QLT ± ih represents the total water transport loss for the lake and pond i under scenario h.
(1) Constraints for water diversion and supplementation amount [23]: (2) Constraints for the water supply capacity [26]: where QAI ± ih represents the local water amount in the lake and pond i under scenario h (m 3 ); QI ± i represents the local water in the lake and pond (m 3 ); QAN ± ih represents the normal supplement of the lake and pond i under scenario h (m 3 ); and QN ± i represents the normal supplement of the lake and pond (m 3 ).
(3) Constraints for the functional area [23]: where PLA ± i represents the upper limit of the area of the lake and pond i (hm 2 ). (4) Constraints for the diversion and supplement of the water amount [24]: where QTF ± ih represents the amount of flood resources in the lake and pond i under scenario h (m 3 ). (5) Constraints for water diversion and the supplementation order [24]: where QLI ± ih represents the loss of local water transport for the lake and pond i under scenario h (m 3 ); QLN ± ih represents the loss of water transport for the normal supplement of the lake and pond i under scenario h (m 3 ); QLF ± ih represents the loss of water transport for the flood diversion and supplement of the lake and pond i under scenario h (m 3 ); Y i3 denotes the presence or absence of the reed wetland in the lake and pond i; Y i4 denotes the presence or absence of marsh wetland in the lake and pond i; Y = 0 denotes that the type of wetland is present; and Y = 1 denotes that there is no such type of wetland.
(6) Constraints for the ecological value: where k is the ecosystem service function of the reserve; TEBE ± kh represents the total ecological value of the indicator system under scenario h (10 4 CNY); EBW ± jk represents the ecological value compensation of the ecosystem j (10 4 CNY); and EBW ± jk represents the total ecological value of the index system of the project scheme (10 4 CNY). (7) The non-negative constraint: In the optimization scheme, the water replenishment quantity of each lake and pond is non-negative.

Model Parameters
The present work considered the Toudaogang Reservoir, Boluo Lake, Mobopao, Yuanbaowapao, and Aobaotupao in the Boluo Lake Wetland National Nature Reserve as the study areas. Each of these five study areas were divided into four types, namely fish pond, crab pond, reed wetland, and marsh wetland. Three scenarios, namely, low flow years, normal flow years, and high flow years, were set for the analysis. The scenario probability for each of these scenarios was 0.55, 0.3, and 0.15, respectively. The original data used in the present paper are from the interconnected river systems in the western region of the Jilin Province that had to be implemented.
Lingo18.0 software was employed to solve the optimization model based on the ITSP method for the ecological water replenishment of the Boluo Lake wetland supplement project. The model solution was obtained in two steps. First, the upper bound sub-model was established and solved. The initial area of the different land types in each study area in the first stage and the lower limit of the reduced area of the different land types in each study area in the second stage were calculated. Next, the initial area calculated in the first step was utilized as the constraint of the lower boundary sub-model to further solve the lower boundary sub-model. The objective of the present study was to alleviate the water shortage in the lake wetlands in the Boluo Lake Wetland National Nature Reserve through the complete utilization of the water resources via ecological replenishment to achieve the restoration of the lake wetland area in the reserve, as well as the reconstruction of the ecosystem, while ensuring the ecological service value of the reserve. The ecological replenishment quantity for each lake in the Boluo Lake Wetland National Nature Reserve is presented in Table 1. Additionally, Figure 2 shows a comparison of the recommended ecological water replenishment project and the ITSP model optimization scheme for each lake in the Boluo Lake Wetland National Nature Reserve.  After the optimization, the total ecological water replenishment of the Toudaogang Reservoir was altered to −13.67% and 6.37%, with a decrease in the lower limit of the total ecological water replenishment in the area and an increase in the upper limit compared to the recommended scheme for the project. Meanwhile, the total ecological water replenishment for the Boluo Lake was altered to −86.81% and 22.43%, with a significant decrease in the lower limit of total ecological water replenishment and an increase in the upper limit of the area compared to the recommended scheme for the project. Moreover, the total ecological water replenishment for Mobopao was altered to −94.59% and −9.12%, with both upper and lower limits of the total ecological water replenishment in the area reduced, as compared to the recommended scheme for the project and a significant reduction in the lower limit. Furthermore, the total ecological water replenishment for the Yuanbaowapao was altered to −100.00% and −74.58%, with both upper and lower limits of the total ecological water replenishment in the area significantly reduced compared to the recommended scheme for the project. In addition, the total ecological water replenishment for the Aobaotupao was altered to −80.65% and 48.42%, with a significant decrease in the lower limit of the total ecological water replenishment in the area and a significant increase in the upper limit compared to the recommended scheme for the project. A comprehensive analysis of the changes in the ecological water replenishment in the study area revealed that the changes in the ecological water replenishment varied from one lake bubble to another, although the lower limit of the ecological water replenishment was significantly low in all five bubbles, while the upper and lower limits of the ecological water replenishment for Mobopao and Yuanbaowapao were relatively lower. Overall, the total ecological water replenishment in the Boluo Lake Wetland National Nature Reserve was reduced, and the simulation results for this optimization scheme are presented in the form of intervals in order to enable the decision-makers to allocate and adjust the ecological water replenishment of each lake bubble according to the actual situation while undertaking decisions.

Variation in the Flood Diversion Volume of the Water Replenishment Project for the Boluo Lake Wetland Based on the ITSP Method
In order to ensure safety, flood waters could be diverted and stored in the form of usable inland water to be used as a resource as required [27,28]. In the ecological water replenishment of wetlands, increasing the amount of flood resources utilized is conducive to reducing the amount of water replenishment required in the other methods of water replenishment, thereby achieving the optimal utilization of the water resources of ecological water replenishment [29]. The values of the flood diversion volume from each lake bubble under different scenarios are presented in Table 2. In the recommended scheme for the project, the Toudaogang Reservoir region did not utilize the flood resources, while in the ITSP optimization scheme, the Boluo Lake region did not utilize the flood resources. This was a consequence of optimizing the flood diversion volume for each lake bubble in the Boluo Lake Wetland National Nature Reserve, considering the water diversion sequence and water supply capacity as constraints. In Mobopao, the flood diversion volumes in the three scenarios were in the range of (0.00, 431.10), (0.00, 478.80), and (0.00, 12.60), respectively, while in the Yuanbaowapao, the flood diversion volume ranges in the three scenarios were (0.00, 289.80), (0.00, 289.80), and (0.00, 56.70), respectively. It could be observed that the flood diversion volumes in each of the lake bubbles were varied in the different scenarios, which was mainly because of the high flow years, the water resources being plentiful, and the ecological water replenishment being accordingly adjusted downwards. Since floods are usually diverted by storing the floodwater in reservoirs during the flood season, they are converted into usable water resources [30]. Therefore, it is possible to control the allocation of the flood resources in ecological water replenishment and rationalize it according to the different scenarios. While the flood diversion volume ranges vary from region to region, the upper limit of the total flood diversions remained substantially higher compared to that in the recommended scheme for the project, due to which the results of the optimized scheme were presented in the form of intervals, providing a wider scope for decision-making.

Analysis of the Ecological Water Replenishment Variation in the Boluo Lake Wetland Water Replenishment Project Based on the ITSP Method
Ecological water replenishment is one of the most direct and effective ways to reconstruct the wetland ecosystem and to restore its functions [31,32], thereby resolving the issues of wetland area reduction and effectively restoring the wetland areas. The initial areas of the different functional areas within each lake bubble in the Boluo Lake Wetland National Nature Reserve, the areas of the recommended scheme for the project, and the areas of each ecosystem type after optimization using the ITSP method are presented in Table 3, and the changes in the areas of each ecosystem type after optimization using the ITSP method compared to those in the areas of the original scheme are presented in Figure 3. As presented in Table 3, in the Toudaogang Reservoir, the largest area after optimization was the fish pond restoration area, which reached a value of (0.72, 1.04) × 10 4 hm 2 , accounting for an increase of (2300.00%, 1987.59%), as compared to the initial area. While the lower limits of the reed wetland area and the marsh wetland area were lower than those of the corresponding initial areas, the upper limits were increased by 700.00% and 716.00%, respectively. In the Boluo Lake, the largest after the optimization was the area of fish ponds, with a value of (2.94, 8.00) × 10 4 hm 2 , which was an increase of (326.76%, 515.35%) compared to the initial area. In the Boluo Lake region, except for the crab ponds, the other three functional regions greatly exhibited increased upper and lower limits of their respective areas. In Mobopao, except for the crab pond area that remained unchanged, the other three functional regions exhibited decreased lower limits of their respective areas compared to the initial area, while the upper limits of their areas increased significantly by 449.72%, 700.00%, and 700.00%, respectively. In Yuanbaowapao, no change occurred in the lower limits of the fish pond area and crab pond area after the optimization. The lower limits of two areas, the reed wetland, and the marsh wetland, were reduced compared to the respective initial areas, and the upper limits of three functional regions, except for the fish pond, were increased substantially by 100.00%, 860.00%, and 700.00%, respectively. In Aobaotupao, the reed wetland had the largest restoration area, with a value of (0.48, 1.58) × 10 4 hm 2 , which accounted for an increase of (1500.00%, 2538.46%) compared to the initial area. The lower limit of the fish pond area remained unchanged compared to the original scheme, while the lower limit of the crab pond and marsh wetland areas was reduced; whereas, the upper limits of all three functional areas were increased significantly by 100.00%, 479.64%, and 710.67%, respectively. Overall, the marsh wetlands that present the highest ecological service value were observed to have a decrease in the lower limit of their area in all ecosystem types, with the exception of the Boluo Lake region, which was mainly because soak pond water and reed wetland water were prioritized for the ecological water replenishment to achieve the best restoration of the ecological functions. In comparison to the original scheme, the overall functional area optimized using the ITSP method was greatly increased, indicating that the ecological water replenishment was effective in restoring the wetland area.  Green water is mainly rainwater stored in the soil at the roots of plants, i.e., water vapour that originates from precipitation, is stored in the soil and enters the atmosphere through evaporation and transpiration; blue water is water in rivers, lakes and underground aquifers [33][34][35]. The Boluo Lake Wetland National Nature Reserve has a unique terrain. The reserve is located at the edge of the Songliao Sag in the southeast of the Songliao Plain. The terrain is similar to that of a basin. Owing to the unique topography of this reserve, there is poor circulation of water within the various lake bubbles in the reserve, which renders it difficult to reduce the salinity in the recharged water resources through circulation or exchange after the ecological water replenishment. As a consequence, after ecological water replenishment, the salinity in the water supply could easily interact with the groundwater and thereby affect the groundwater quality, and this salty groundwater would enter into the soil, causing soil salinization and ultimately affecting plant growth. This means that the amount of blue water is increased by ecological rehydration and that the different types of blue water circulate with each other and the salts interact. The blue water enters the soil and is also converted into green water to be stored by the plant roots, while at the same time the salts are absorbed by the plant roots. It is reported that reeds may effectively prevent the salinization of the land [36][37][38]. The salt in the water resources within the Boluo Lake Wetland National Nature Reserve could be absorbed by the reeds and removed through harvesting every year. In the optimization scheme proposed in the present report, the total amount of ecological water replenishment was reduced, thereby reducing the total salinity introduced. In addition, for the best ecological function restoration effect, the ecological water replenishment was prioritized for the water used in the pond and the reed wetland, and the area of the reed wetland was restored. As a consequence, the effectiveness of removal of water salinity through reed adsorption in the reserve was also improved. Figure 3 depicts the change in the functional area of each land type in the different study areas compared to the area in the recommended scheme for the project. In comparison to the functional area in the recommended scheme for the project, the lower limits of the respective areas of the fish and crab pond areas in the Toudaogang Reservoir, the marsh wetlands in the Boluo Lake, the fish ponds in Mobopao, the crab ponds in Yuanbaowapao, and the fish ponds and reed wetlands in Aobaotupao were increased in the optimized simulation scheme. In addition, the lower limits of the areas of the fish ponds in the Boluo Lake, the crab ponds in Mobopao, and the crab ponds in Aobaotupao were reduced, while the lower limits of the areas of the other functional areas remained unchanged. In the optimization scheme, the upper limit of the fish pond area of the Toudaogang Reservoir was increased, the upper limit of the crab ponds area in the Toudaogang Reservoir, the crab ponds area in the Boluo Lake, the fish ponds area in Mobopao, the fish and crab ponds areas in Yuanbaowapao, and the fish ponds area in Aobaotupao remain unchanged, and the upper limit of the areas of the other functional areas was reduced. It could be observed that the restoration areas were different for the different lakes and different functional areas, and the function of the Boluo Lake Wetland National Nature Reserve had been adjusted because the model was optimized and simulated using the value of ecological services as a constraint.

Analysis of the Change in the Ecological Service Value of the Boluo Lake Wetland Water Replenishment Project Based on the ITSP Method
The Boluo Lake wetland water replenishment project considered the project benefits while reconstructing the wetland ecosystem and restoring its functions. The project benefits are mainly reflected in the level of restoration of the wetland area and its functions and the value of ecosystem services [39]. The value of ecosystem services refers to the benefits obtained as a result of maintaining the ecosystem functions of the earth's life support system [40]. The recommended scheme for the project regarding the value of ecological services under different indicator systems and the results of optimization using the ITSP method in the Boluo Lake Wetland National Nature Reserve are presented in Figure 4. In the present study, 15 types of ecological service functions were analyzed and the main ecological service value sequence generated was: plant adsorption (150.07, 172.88) × 10 6 CNY) > microclimate adjustment (141.24, 162.71) × 10 6 CNY) > fish (125.56, 147.72) × 10 6 CNY) > pollution absorption capacity (80.25, 92.45) × 10 6 CNY) > crab (46.33, 54.50) × 10 6 CNY) > tourist development (29.69, 34.21) × 10 6 CNY) > volume of flood water mobilized (28.89, 33.28) × 10 6 CNY) > oxygen release (18.50, 21.36) × 10 6 CNY) > biodiversity maintenance (10.43, 12.02) × 10 6 CNY) > urban landscape (8.03, 9.24) × 10 6 CNY) > reed (5.47, 6.32) × 10 6 CNY) > natural landscape (3.21, 3.70) × 10 6 CNY) > scientific culture (2.01, 2.31) × 10 6 CNY) > carbon sequestration (1.90, 2.19) × 10 6 CNY) > water supply (0.00, 0.00) × 10 6 CNY).In this sequence, the top two indicators, namely plant adsorption, and microclimate adjustment, as well as the fourth indicator, i.e., pollution absorption capacity, are all indicators of the environmental regulation of the ecosystem, which suggests that the ecological water replenishment project in the present study focused on environmental benefits. The two indicators capable of generating direct economic value, namely fish and crab, are in the third and fifth place, respectively, indicating that the ecological water replenishment project in the present study does not just focus on the ecological benefits, and rather completely utilizes the local economic benefits while focusing on the ecological benefits alongside. With the exception of the three ecological service functions of fish, crab, and water supply, which remained unchanged compared to the recommended scheme for the project, the upper and lower limits of the ecological service values of the other ecological service functions were improved. In comparison to the ecological service value in the recommended scheme for the project, the largest improvement was observed in the carbon sequestration, with a value of (23.90%, 21.58%), while the upper and lower limits of the ecological service value of the remaining ecological service functions were increased by over 15%. This indicates that compared to the recommended scheme for the project, the ITSP method optimization scheme obtained a higher ecological service value while ensuring the restoration of the ecosystem functions. The ITSP-based model constructed in the present study could be used in a wide range of applications, including the ecological water replenishment projects in different regions across the world. It is noteworthy that during the application of the constructed model, the relevant influencing factors should be adjusted according to the local characteristics, such as climate, topography, and land type of the concerned region. In addition, when applying the ITSP model, the planning area should be redefined for the size of the concerned study area, suitable scenarios should be designed, the water diversion sequence should be established according to the local water resource structure, the proportion of water supply through various modes should be adjusted according to the local water supply capacity, and the ITSP model should be customized for the actual local situation, to provide the decisionmakers with a decision plan that is most consistent with the actual local situation. For instance, when applying the ITSP model to the Momoge Reserve, which has a large area, it would be suitable to study 18 lake bubbles in the reserve and to establish three scenarios for analysis, namely a low year, a medium year, and a high water year. Moreover, when applying the ITSP model, a suitable ecological service indicator system should be selected for the assessment of the ecological service values based on the main ecological functions of the Momoge Reserve, and relevant adjustments should be incorporated in the constraints according to the actual situation of this reserve.

Conclusions
The main objective of the present study was to alleviate the water shortage situation in the lake bubble wetlands within the Boluo Lake Wetland National Nature Reserve, which was achieved by constructing an ecological water replenishment optimization model for the Boluo Lake Wetland replenishment project, which is based on the ITSP method and considered the minimum ecological water replenishment amount as the target. The ITSP method optimized the proposed scheme by scientifically allocating the ecological water resources and increasing the amount of flood diversion to achieve the restoration of the lake bubble wetland area and the reconstruction of the wetland ecosystem in the reserve with a minimum possible amount of ecological water replenishment. A lower ecological water replenishment implies that the amount of salt introduced due to the ecological water replenishment pathway is lower and the restoration of the reed wetland area would remove some of the salt introduced by the ecological water replenishment, thereby effectively preventing the salinization of the land in the reserve. In addition, the value of ecological services in the reserve was considered to ensure both the restoration of the ecosystem functions and the economic benefits. The ITSP method-based ecological water replenishment optimization model for the Boluo Lake wetland water replenishment project constructed in the present study may be widely applied to various ecological water replenishment projects, thereby providing the decision-makers with a feasible decision plan.

Data Availability Statement:
The data presented in this study is available on request from the corresponding author.

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