Impacts of Water Resources Allocation on Water Environmental Capacity under Climate Change

: Water environmental capacity (WEC) is an essential indicator for effective environmental management. The designed low water ﬂow condition is a prerequisite to determine WEC and is often based on the stationarity assumption of low water ﬂow series. As the low water ﬂow series has been remarkably disturbed by climate change as well as reservoirs operation and water acquisition, the stationarity assumption might bring risk for WEC planning. As the reservoir operation and water acquisition under climate change can be simulated by a water resources allocation model, the low water ﬂow series outputted from the model are the simulations of the disturbances and often show nonstationary conditions. After estimating the designed low water ﬂow through nonstationary frequency analysis from these low water ﬂow series, the WEC under the nonstationary conditions can be determined. Thus, the impacts of water resources allocation on WEC under climate change can be quantitatively assessed. The mid-lower reaches of the Hanjiang River basin in China were taken as a case study due to the intensive reservoir operation and water acquisition under the climate change. A representative concentration pathway scenario (RCP4.5) was employed to project future climate, and a Soil and Water Assessment Tool (SWAT) model was employed to simulate water availability for driving the Interactive River-Aquifer Simulation (IRAS) model for allocating water. Water demand in 2016 and 2030 were selected as baseline and future planning years, respectively. The results show that water resources allocation can increase the amount of WEC due to amplifying the designed low water ﬂow through reservoir operation. Larger regulating capacities of water projects can result in fewer differences of WEC under varied water availability and water demand conditions. The increasing local water demand will decrease WEC, with less regulating capacity of the water projects. Even the total available water resources will increase over the study area under RCP4.5. More water deﬁcit will be found due to the uneven temporal-spatial distribution as well as the increasing water demand in the future, and low water ﬂow will decrease, which further leads to cut down WEC. Therefore, the proposed method for determining the WEC can quantify the risk of the impacts of water supply and climate change on WEC to help water environmental management.


Introduction
With the increasing population and rapid expansion of industry and agriculture, an enormous amount of sewage has been drained into rivers and aggravated the deterioration of water quality, which has resulted in serious water pollution [1,2]. In order to alleviate water deterioration and safeguard water use, a water quality control scheme is implemented according to the water quality protection target. Water environmental capacity (WEC) indicates the yearly total maximum pollutant load that is permitted to drain into river without water quality deterioration in water quality control scheme. WEC is often taken as an indicator for the sustainable development and the ecological stability [1,3,4]. The amount of WEC is not only determined by the water quality protection target, but also dependent on the designed low water flow. The designed low water flow is often chosen by the minimum monthly average water flow with an annual exceedance probability [5,6] or the seven-day average low flow within a T-year return period [7,8]. Thus, the designed low water flows are determined by frequency analysis that are often based on historical natural low water flow series and its stationarity assumption. As the stationary assumption has been violated by human activities or climate change, a lot of new hydrological frequency analyses for nonstationary conditions have been developed. Cooley eschewed the term return period and instead communicates yearly risk in terms of a probability of exceedance to investigate return level for hydrological extremes under nonstationary climates [9]. Du et al. improved the characterization of the nonstationary return period and risk under the expected number of exceedances interpretation by employing meteorological covariates in the nonstationary frequency analysis [10]. Rootzén and Katz proposed the design life levels for quantifying the probability with exceeding a fixed threshold for nonstationary climates [11]. Obeysekera et al. adopted the expected occurrence frequency of a given hydrological extreme under nonstationary climates for projects designs [12]. However, few of them have simultaneously considered the impacts of climate change and human activities (i.e., reservoirs operation and water acquisition) on the designed low water flow, which probably brings risk for the planning of water environmental management.
As climate change is considered as one of the main driving forces for the variation of water flow [13,14], the temporal-spatial distribution of available water resources will become more uneven under future scenarios [15][16][17], which will bring daunting challenges in water supply and cause water deficit problems, especially on the boundary areas of a basin [18]. The variation of available water resources is prone to lead great impacts on low water flow series and its designed low water flow conditions for WEC [19]. Even the responses of runoff to climate change have been concerned [20,21], the impacts of human activities on low water flow series are rarely discussed. As the amount of water supply is dependent on the water project (e.g., reservoir) operational rules, the temporalspatial distribution of natural water flow for water availability can also be altered by water projects operation [22]. With the operation of reservoirs, the connections between water flows at different locations have been disturbed, resulting in the nonstationary time-varying marginal distributions between different hydrological series [23], especially for the low water flow series. With the rapid expansion of industry and agriculture in the future, the water demand will dramatically increase [24,25] and intensify the alteration of water flow, which alters the stationary assumption of low water flow series frequency analysis for the designed low water flow. Therefore, the impacts of climate change and human activities on low water flow series and the designed low water flow should be assessed.
The water resources allocation model simultaneously incorporates the water projects operation including reservoirs operation and water acquisition for water allocation [26]. The water resources allocation model has become an effective tool for quantitatively assessing the impacts of climate change and human activities on low water flow series and the designed low water flows [27][28][29] and further the impacts on WEC. The aim of our study is to propose a framework for estimating WEC under the nonstationary low water flow conditions that are caused by climate change as well as reservoir operation and water acquisition. The remainder of this paper is organized as follows: Section 2 introduces the framework for assessing the impacts of water resources allocation on WEC under climate change. Section 3 describes the methodologies applied in mid-lower reaches of Hanjiang River basin in China. Section 4 represents the results of projected water availability, water demand, water resources allocation, designed low water flow under stationary and nonstationary conditions and their corresponding WECs. The impacts of Water 2021, 13, 1187 3 of 20 water resources allocation on WEC have also been discussed in this section. Section 5 gives the conclusions.

Framework for Assessing the Impacts of Water Resources Allocation on Water Environmental Capacity under Climate Change
The framework for assessing the impacts of water resources allocation on WEC under climate change can be divided into four main modules (shown in Figure 1): (1) water availability estimation module, (2) water demand projection module, (3) water resources allocation module and (4) water environmental capacity determination module.

Water Availability Estimation Module
The hydrological model is considered as the most effective tool to assess the available water resources under different climate scenarios. The SWAT model is developed by the Agricultural Research Service of the United States Department of Agriculture and has become one of the most popular distributed hydrological models around the world [37,38]. Hydrological response unit (HRU) is the basic unit in SWAT model. The routine of water is simulated from the HRUs to the sub-basin level according to the daily water cycles. SWAT model has been proved to be an excellent tool for accessing the impacts of climate change on runoff [39,40]. In order to project the future water resources through SWAT at basin scale, outputs from the GCMs at global scale [41,42] are downscaled by DBC method to derive empirical relationships between the global scale GCM predictors and the basinscale predictors [32]. Therefore, DBC method and SWAT model are the primary tools for water availability estimation. The function of the water availability estimation module is estimating the water availability under historical and future scenarios. As the outputs of the General Circulation Model (GCM) have been widely used across the world for projecting the future climate [30,31], they are inputted into Soil and Water Assessment Tool (SWAT) hydrological model for predicting the water availability in the future. In order to adjust the outputs of large scale GCM for satisfying the inputs of the small-scale SWAT model, an effective statistical downscaling method (i.e., Daily Bias Correction (DBC) method) is adopted in this module [32]. The DBC method takes into account the different changes in frequency distributions of daily precipitation, and corrects the wet-day frequency of precipitation by using a small precipitation threshold [32]. The water demand module is composed of the quota method and the Tennant method. The quota method is used for projecting the water demand in socioeconomic use sectors in the planning year, and the Tennant method can determine the water demand of in-stream ecology [33,34]. The function of the water resources allocation module is to deploy water resource to every water use sector according to their water demands and water availabilities. The Interactive River-Aquifer Simulation (IRAS) model is adopted here due to its flexibility and accuracy in simulation [35]. The disturbed low water flow series outputted from water resources allocation module and the natural low water flow series will be taken as the inputs to determine the designed low water flow for water environmental capacity determination module [36]. Based on the designed low water flows, the impacts of water resources allocation on WEC under climate change will be figured out.

Water Availability Estimation Module
The hydrological model is considered as the most effective tool to assess the available water resources under different climate scenarios. The SWAT model is developed by the Agricultural Research Service of the United States Department of Agriculture and has become one of the most popular distributed hydrological models around the world [37,38]. Hydrological response unit (HRU) is the basic unit in SWAT model. The routine of water is simulated from the HRUs to the sub-basin level according to the daily water cycles. SWAT model has been proved to be an excellent tool for accessing the impacts of climate change on runoff [39,40]. In order to project the future water resources through SWAT at basin scale, outputs from the GCMs at global scale [41,42] are downscaled by DBC method to derive empirical relationships between the global scale GCM predictors and the basin-scale predictors [32]. Therefore, DBC method and SWAT model are the primary tools for water availability estimation.

Water Demand Projection Module
Water uses are often classified into socioeconomic (also called off-stream) use and in-stream ecology use in a water allocation model. Socioeconomic use consists of municipal, rural, industrial, agricultural categories in the water demand projection module.

Socioeconomic Water Demand
Quota method is widely applied to project the annual socioeconomic water demand in a planning year [43]. The amounts of water demand for the four water use sectors can be obtained by Equation (1).
where Wd t i,j is the amount of water demand for the jth use sector in the ith operational zone at the tth time step; Wq t i,j is the water quota unit of water demand use sector; A t i,j is the amount of water units in water use sector; UR t i,j is the utilization rate of water use sector. The water quota units are the amount of water consumption per capita in municipal and rural use sectors, the amount of water consumption per ten thousand Yuan in industrial use sector and the amount of net irrigation water per unit area in agricultural use sector, respectively; the number of water units are the projected population in municipal and rural use sectors, the projected GDP in industrial use sector and the projected irrigated area in agricultural use sector.

In-Stream Ecology Water Demand
The in-stream ecological flow is taken account to protect the local ecological environment satisfying the ecological use and to support the survival of aquatic wildlife in operational zones. Tennant method, a popular method to estimate in-stream ecological water demand [34,44], is applied through the product of the annual average runoff and proportions during the flood season and non-flood season: where Wde t i is the in-stream ecological water demand in the ith operational zone at the tth time step; R i is the average annual runoff; ε i is the proportion coefficient of minimum ecological environment water demand.

Water Resources Allocation Module
Based on the outputs from water availability estimation module and water demand projection module, water resources allocation module will deploy the available water resources to every water use sector and the in-stream water flow along with river segments. Developed by Cornell University [45], IRAS is a rule-based water system simulation model for water resources allocation and is adopted here due to its flexibility and accuracy in simulation [35].
As the water system is often composed of water transfer, consumption and loss components, it can be sketched by node network topology. Two types of nodes, named the reservoir node and demand node, play vital roles and are directly corresponding to the "supply-side" and "demand-side", respectively. Water release from reservoir node is based on the reservoir operation rules while the water deficit at demand node can be estimated by the following water balance Equations (3) and (4).
where t is the current time step; tst is the total number of the sub-time-step; st is the current sub-time-step; W e st is the amount of water from natural inflow; W dem t is the water demand of the water use sector; f is the demand reduction factor; W d st is the water deficit.

Water Environmental Capacity Determination Module
The disturbed low water flow series can be derived from water resources allocation module while the natural low water flow series can be selected from observed data or from the simulated runoff by SWAT model. The minimum monthly average water flow with an annual exceedance probability of 90% is defined as the designed low water flow for WEC. WEC in every river segment with their corresponding water quality protection target C s can be estimated by Equation (5).
where E is WEC in a river segment; C s is the water quality protection target; C 0 is pollutant background concentration; k is the pollutant attenuation coefficient; x is the river segment length; u is flow velocity related to the designed low water flow; Q 0 is the designed low water flow; q is water flow from the sewage.

Study Area
The mid-lower reaches of Hanjiang River basin is located in the Hubei province of central China (shown in Figure 2), and its area is about 63,800 km 2 . The water from the upper reaches of basin flows into the Danjiangkou reservoir and then releases into the study area. The Danjiangkou reservoir is the water source of the middle route of the south-north water transfer project in China [46], our study area has been impacted by both climate change and water acquisition from the water transfer project. Moreover, due to the fast urbanization and the rapid development of the economy, the local water demand is increasing, and water quality is deteriorating. Therefore, the impacts of climate change and According to the main water use and its water sources, the study area can be classified into three types: type A that is the irrigation area and acquires water from the mainstream of the Hanjiang River; type B that is also the irrigation area which diverts water from the mainstream of Changjiang River through the Three Gorges Reservoir-Hanjiang water transfer project (under planning); type C that is cities while acquires water from the mainstream of Hanjiang River. 40 operational zones are sketched based on the local water systems, and there are 13, 11 and 16 operational zones in type A, B and C, respectively. There are seven medium or large size reservoirs (the total storage volume is 36.7 billion m 3 ) for regulating water flows (shown in Figure 2). In order to illustrate the connections between operational zones and river systems in IRAS model, the water system can be sketched as  According to the main water use and its water sources, the study area can be classified into three types: type A that is the irrigation area and acquires water from the mainstream of the Hanjiang River; type B that is also the irrigation area which diverts water from the mainstream of Changjiang River through the Three Gorges Reservoir-Hanjiang water transfer project (under planning); type C that is cities while acquires water from the mainstream of Hanjiang River. 40 operational zones are sketched based on the local water systems, and there are 13, 11 and 16 operational zones in type A, B and C, respectively. There are seven medium or large size reservoirs (the total storage volume is 36.7 billion m 3 ) for regulating water flows (shown in Figure 2). In order to illustrate the connections between operational zones and river systems in IRAS model, the water system can be sketched as Figure 3.

Description of Water Demand and Availability Schemes
To figure out the potential impacts of water resources allocation on WEC, six schemes are set according to the varied water availability and water demand conditions (shown in Table 1). The schemes I, II, III and IV are set by the combination of two water availability scenarios (i.e., history and future) and water demands in baseline (2016) and planning years (2030). In order to figure out the impacts of human activities on the WEC, schemes V and VI with historical and future natural water flows are set as reference schemes without water resources allocation.

Description of Water Demand and Availability Schemes
To figure out the potential impacts of water resources allocation on WEC, six schemes are set according to the varied water availability and water demand conditions (shown in Table 1). The schemes I, II, III and IV are set by the combination of two water availability scenarios (i.e., history and future) and water demands in baseline (2016) and planning years (2030). In order to figure out the impacts of human activities on the WEC, schemes V and VI with historical and future natural water flows are set as reference schemes without water resources allocation. basin, which are provided by China Meteorological Administration. 261 precipitation gauging stations charged by the Bureau of Hydrology affiliated to the Changjiang Water Resources Commission have also provided daily precipitation data from 1980 to 2000 to drive SWAT model. The daily discharge data from four runoff gauging stations from 1980 to 2000 are also collected. The outputs from the Beijing Climate Center Climate System Model (BCC-CSM) under the representative concentration pathway RCP4.5 scenario are adopted to project future climate change, due to its good performance in reproducing observed climate variations of temperature and precipitation in China [47].

Socioeconomic Water Demand Data
Water consumption ratios, social and economic development projection for each operational zone in the baseline year (2016) and planning year (2030) are provided by the Hubei Provincial Department of Water Resources. Their water demands are projected by water demand projection module. Water uses in operational zones in type C area are only municipal and industrial uses while there are four types of water use in the operational zones in type A and type B areas. The water flow data for estimating in-stream water flow demand are outputted from the SWAT model.

Characteristics of Reservoirs
The characteristics of the seven reservoirs are listed in Table 2 and their operation rules are collected from the dispatching schedules of Hubei provincial large reservoirs compiled by Hubei Provincial Department of Water Resources [48].

Water Quality Data
Water quality data from the nine gauge stations along the mainstream are obtained from Hubei Provincial Department of Water Resources to force the water environmental capacity determination module. These nine stations divide the main stream into eight river segments (shown in Figure 4). Chemical oxygen demand (COD Mn ) and ammonium nitrogen (NH 3 -N) are taken as water quality indicators according to the water quality control scheme. The total actual pollution discharges of COD Mn and NH 3 -N in eight segments in 2016 were obtained from the environmental pollution source investigation report provided by Hubei Environmental Protection Bureau [49]. The total discharge of COD Mn and NH 3 -N in 2016 are 307,209 t and 8961 t, respectively. According to the technical guideline of environmental capacity of surface water in China complied by the Chinese Academy for Environmental Planning [50], the parameter k for COD Mn and NH 3 -N are set to be 0.25 and 0.10, respectively. trol scheme. The total actual pollution discharges of CODMn and NH3-N in eight segments in 2016 were obtained from the environmental pollution source investigation report provided by Hubei Environmental Protection Bureau [49]. The total discharge of CODMn and NH3-N in 2016 are 307,209 t and 8961 t, respectively. According to the technical guideline of environmental capacity of surface water in China complied by the Chinese Academy for Environmental Planning [50], the parameter k for CODMn and NH3-N are set to be 0.25 and 0.10, respectively.

Water Availability under Historical and Future RCP4.5 Scenarios
In order to project the future water availability, the DBC method is calibrated firstly by historical data. The relative errors of precipitation are less than 15% and the absolute errors for maximum and minimum temperature are less than 1.5 °C, so the calibrated DBC can be applied in downscaling future precipitation and surface temperature. The SWAT model is calibrated and validated by the daily discharge data under historical scenario from 1980 to 2000. The period  is split into calibration period (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993) and validation period (1994)(1995)(1996)(1997)(1998)(1999)(2000), respectively. The Nash-Sutcliffe efficiency coefficient (NSE) values for hydrological simulation at the four runoff gauging stations are listed in Table 3. The NSE values during the calibration and validation periods range from 0.66 to 0.93, suggesting that the calibrated SWAT model can apply to project the future runoff. When the calibrated SWAT model is applied in both the historical period (1956-2016) and the future period (2020-2080), the amounts of local available water resources during a year, flood season (May-September) and non-flood season (January-April and October-December) in type A and type B areas can be estimated (listed in Table 4). The total annual available water resources will increase from 7.03 billion m 3 in historical period to 8.01 billion m 3 under RCP4.5 scenario due to the increasing precipitation. However, the available water resources in non-flood season will decrease from 1.76 billion m 3 to 1.57 billion m 3 . Thus, the temporal distribution of available water resources will become more uneven.

Water Availability under Historical and Future RCP4.5 Scenarios
In order to project the future water availability, the DBC method is calibrated firstly by historical data. The relative errors of precipitation are less than 15% and the absolute errors for maximum and minimum temperature are less than 1.5 • C, so the calibrated DBC can be applied in downscaling future precipitation and surface temperature. The SWAT model is calibrated and validated by the daily discharge data under historical scenario from 1980 to 2000. The period  is split into calibration period (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993) and validation period (1994)(1995)(1996)(1997)(1998)(1999)(2000), respectively. The Nash-Sutcliffe efficiency coefficient (NSE) values for hydrological simulation at the four runoff gauging stations are listed in Table 3. The NSE values during the calibration and validation periods range from 0.66 to 0.93, suggesting that the calibrated SWAT model can apply to project the future runoff. When the calibrated SWAT model is applied in both the historical period (1956-2016) and the future period (2020-2080), the amounts of local available water resources during a year, flood season (May-September) and non-flood season (January-April and October-December) in type A and type B areas can be estimated (listed in Table 4). The total annual available water resources will increase from 7.03 billion m 3 in historical period to 8.01 billion m 3 under RCP4.5 scenario due to the increasing precipitation. However, the available water resources in non-flood season will decrease from 1.76 billion m 3 to 1.57 billion m 3 . Thus, the temporal distribution of available water resources will become more uneven. The available water resources in the type C area are from the main stream of the Hanjiang River and are dependent on the inflow and the operational rule of the Danjiangkou reservoir. Figure 5 shows the percentage changes of the amount of available water resources under RCP4.5 scenario comparing to that under historical scenario. The available water resources in type A area will significantly increase with more than 10%, while decreasing trends are found in type B area with a range from 5 to 10%. As the geographic topography serves as a barrier to block the east monsoon winds going west, the precipitation is going to decrease in the west of mountain areas (i.e., type B area). The uneven spatial distribution of available water resources will be exacerbated by climate change according to the results of runoff projection in the future scenarios. under RCP4.5 scenario comparing to that under historical scenario. The available water resources in type A area will significantly increase with more than 10%, while decreasing trends are found in type B area with a range from 5 to 10%. As the geographic topography serves as a barrier to block the east monsoon winds going west, the precipitation is going to decrease in the west of mountain areas (i.e., type B area). The uneven spatial distribution of available water resources will be exacerbated by climate change according to the results of runoff projection in the future scenarios.

Water Demand in 2016 and 2030
The water demands in 2016 and 2030 are estimated by the water demand projection module, with the results are listed in Table 5. The percentages of water demand among socioeconomic water use are shown in Figure 6. The in-stream ecological water demand in every node has been estimated on basis of the runoff series outputted from the SWAT model. Figure 7a shows the in-stream ecological water demand in type A and type B areas. The in-stream ecological water demand in type C is set as the same as the value in type A at their sharing nodes. The Huangjiagang station (shown in Figure 7b) is set as a typical node at main stream of mid-lower reaches of the Hanjiang River for estimating the in-stream ecological water demand.

Water Demand in 2016 and 2030
The water demands in 2016 and 2030 are estimated by the water demand projection module, with the results are listed in Table 5. The percentages of water demand among socioeconomic water use are shown in Figure 6. The in-stream ecological water demand in every node has been estimated on basis of the runoff series outputted from the SWAT model. Figure 7a shows the in-stream ecological water demand in type A and type B areas. The in-stream ecological water demand in type C is set as the same as the value in type A at their sharing nodes. The Huangjiagang station (shown in Figure 7b) is set as a typical node at main stream of mid-lower reaches of the Hanjiang River for estimating the in-stream ecological water demand.    The annual water demand will gradually increase from 11,734 million m 3 in 2016 to 13,897 million m 3 in 2030 (listed in Table 5). The municipal water demand increases from 946 million m 3 (the percentage is 8.06%, hereinafter the same) to 1178 million m 3 (8.48%), while the rural water demand declines from 184 million m 3 (1.56%) to 170 million m 3 (1.22%). The industrial water demand increases from 3467 million m 3 (29.55%) to 5973 million m 3 (42.98%). Although the agricultural water demand slightly declines, it still occupies the largest percentage among the water uses (i.e., 60.82% in 2016 and 47.32% in 2030). As the Hanjiang River basin is a pilot area for implementing the strictest water resources control system for water resources management policy in China, the water consumption will be restricted by this water use efficiency control policy. The water demand    The annual water demand will gradually increase from 11,734 million m 3 in 2016 to 13,897 million m 3 in 2030 (listed in Table 5). The municipal water demand increases from 946 million m 3 (the percentage is 8.06%, hereinafter the same) to 1178 million m 3 (8.48%), while the rural water demand declines from 184 million m 3 (1.56%) to 170 million m 3 (1.22%). The industrial water demand increases from 3467 million m 3 (29.55%) to 5973 million m 3 (42.98%). Although the agricultural water demand slightly declines, it still occupies the largest percentage among the water uses (i.e., 60.82% in 2016 and 47.32% in 2030). As the Hanjiang River basin is a pilot area for implementing the strictest water resources control system for water resources management policy in China, the water consumption will be restricted by this water use efficiency control policy. The water demand The annual water demand will gradually increase from 11,734 million m 3 in 2016 to 13,897 million m 3 in 2030 (listed in Table 5). The municipal water demand increases from 946 million m 3 (the percentage is 8.06%, hereinafter the same) to 1178 million m 3 (8.48%), while the rural water demand declines from 184 million m 3 (1.56%) to 170 million m 3 (1.22%). The industrial water demand increases from 3467 million m 3 (29.55%) to 5973 million m 3 (42.98%). Although the agricultural water demand slightly declines, it still occupies the largest percentage among the water uses (i.e., 60.82% in 2016 and 47.32% in 2030). As the Hanjiang River basin is a pilot area for implementing the strictest water resources control system for water resources management policy in China, the water consumption will be restricted by this water use efficiency control policy. The water demand in agriculture thus decreases. However, water demand in industry still increases enormously, even taking the water-saving technologies into consideration due to the rapid social and economic development. Therefore, the water supply security for social economic sustainable development will be seriously challenged.

Results of Water Resources Allocation Model under Climate Change
According to the Water Resources Management Plan over Hanjiang River Basin [51], the water demand for domesticity and ecology should be firstly ensured. Priorities for the water use from high to low are assigned as municipal use, rural use, in-stream ecology use, industrial use and agricultural use, respectively. Taking ten days as the time step, the allocated water resources for water use sectors in operational zones are obtained from the outputs of the IRAS model. The results under the schemes I, II, III and IV are listed in Table 6. The total water deficits under the four schemes are 977, 1635, 1160, 1780 million m 3 , respectively. The amount of increasing available water resources in the future (0.98 billion m 3 , from 7.03 billion m 3 to 8.01 billion m 3 ) is much less than the amount of dramatic increase of water demand (2.17 billion m 3 , from 11.73 to 13.90 billion m 3 ). Thus, the water deficit under future water demand scenario is obviously higher than that under historical scenario. Moreover, with the more uneven temporal-spatial distribution of available water resources under future water availability scenario (shown in Table 4 and Figure 5), more water stress will be directly exerted on reservoirs operation for increasing water supply. The risk of water deficit will increase (scheme III (9.40%) > scheme I (7.91%) and scheme IV (12.27%) > scheme II (11.27%)). The water deficits are mainly found in industry and agriculture while the other water uses can be satisfied. The water deficit rates of agriculture are even less than those of industry (scheme III and scheme IV), although the industry has a higher priority. Because the water consumption in agriculture mainly happens in flood season (according to water demand projection results, the water demands from May to September occupy 80.49% and 78.87% in 2016 and 2030, respectively). Water demand in agriculture matches the inflow process fitly (shown in Table 4) and is effectively satisfied by natural water flow.
In order to analyze the spatial distribution of water deficit, the water deficit rates of every operational zone in type A and type B areas are presented in Figure 8. The water deficit rates are less than 1% in type C area and less than 5% in type A area. As the operational zones in type A and type C areas are located in the downstream of Danjiangkou reservoir, they can extract the water released from Danjiangkou reservoir. Even the water demand will dramatically increase in the scheme II and IV, there are few water deficits in these two type areas. However, water deficits in type B area are serious, as its water supply are only from the local available water resources and the local reservoirs without enough storages, which cannot meet their increasing demand. ational zones in type A and type C areas are located in the downstream of Danjiangkou reservoir, they can extract the water released from Danjiangkou reservoir. Even the water demand will dramatically increase in the scheme II and IV, there are few water deficits in these two type areas. However, water deficits in type B area are serious, as its water supply are only from the local available water resources and the local reservoirs without enough storages, which cannot meet their increasing demand.

Water Environmental Capacity Response to Water Resources Allocation under Climate Change
There are two types of low water flow series. One is outputted from water resources allocation model. The other one is natural low water flow series. The designed low water flows are determined by frequency analysis on these two low water flow series. As the mainstream can be divided into eight river segments (shown in Figure 4), there are eight designed low water flows (listed in Table 7). The background concentrations for pollutants are set as the maximum value to meet their water protection targets (i.e., the class II in our study area according to the water environment function area in Hanjiang River basin). WEC in every segment is subsequently determined by the water environmental capacity determination module (listed in Table 8). To achieve water quality protection

Water Environmental Capacity Response to Water Resources Allocation under Climate Change
There are two types of low water flow series. One is outputted from water resources allocation model. The other one is natural low water flow series. The designed low water flows are determined by frequency analysis on these two low water flow series. As the mainstream can be divided into eight river segments (shown in Figure 4), there are eight designed low water flows (listed in Table 7). The background concentrations for pollutants are set as the maximum value to meet their water protection targets (i.e., the class II in our study area according to the water environment function area in Hanjiang River basin). WEC in every segment is subsequently determined by the water environmental capacity determination module (listed in Table 8). To achieve water quality protection target, the permissible discharges of pollutants in these eight river segments can be estimated by deducting discharged pollutants from WECs as listed in Table 9.   Table 9. c The permissible COD Mn and NH 3 -N in eight river segments (t/a). The designed low water flows along the main streams under the six schemes are quite different. Specifically, the designed low water flows under the schemes I, III, V and VI have notably increased in the river segment S2 and then have fluctuated slightly in the downstream river segments. The designed low water flows under the schemes II and IV are found to significantly decline in the river segment S2. As the industrial city, Xiangyang, locates near the river segment S2 (shown in Figure 3) and its industrial water demand is going to enormously increase in 2030, large amount of water is extracted from the main stream to satisfy the demand all the year. Thus, the low water flow is reduced and then WECs decrease. The maximum WEC is also found in the river segment S2. Specifically, the WECs for COD Mn under schemes I, II, III, IV, V and VI are 158,948, 138,636, 157,390, 135,565, 100,412 and 95,520 t/a, respectively. For NH 3 -N, the WECs are 2904, 2519, 2875, 2461, 1768 and 1663 t/a, respectively. However, as the discharged pollutants have exceeded its WEC in 2016, there are no more permissible COD Mn and NH 3 -N. Similarly, no more permissible NH 3 -N have also been found in the river segment S7. Therefore, the pollution discharge should be reduced in the river segments S2 and S7.

Pollutant
To further investigate the impacts of the varied water availability and water demand conditions on the WEC, the percentage changes of WECs are calculated. Schemes I and III, II and IV as well as V and VI are paired to analyze the response to variation of water availability, while the schemes I and II as well as III and IV are paired to investigate the response to different water demand scenarios (shown in Figure 9). Three types of river segments are classified based on their WECs (shown in Figure 10). The river segment where the WEC is slightly responsive to the variations of water availability and water demand is defined as type (i). The river segment S1 is located in the reach between Huangjiagang and Jiangjiaying stations. The river segments S2, S3, S4, S5, S6 are taken as type (ii) where their WECs are sensitive to the changing water demand scenarios and not to the variation of water availability. These river segments are located in the reach between Jiangjiaying and Yuekou stations. The river segments S7 and S8, are defined as type (iii) where their WECs are sensitive to both the variations of water availability and water demand, and these two river segments are located between Yuekou and Wuhan stations. The designed low water flows along the main streams under the six schemes are quite different. Specifically, the designed low water flows under the schemes I, III, V and VI have notably increased in the river segment S2 and then have fluctuated slightly in the downstream river segments. The designed low water flows under the schemes II and IV are found to significantly decline in the river segment S2. As the industrial city, Xiangyang, locates near the river segment S2 (shown in Figure 3) and its industrial water demand is going to enormously increase in 2030, large amount of water is extracted from the main stream to satisfy the demand all the year. Thus, the low water flow is reduced and then WECs decrease. The maximum WEC is also found in the river segment S2. Specifically, the WECs for CODMn under schemes I, II, III, IV, V and VI are 158,948, 138,636, 157,390, 135,565, 100,412 and 95,520 t/a, respectively. For NH3-N, the WECs are 2904, 2519, 2875, 2461, 1768 and 1663 t/a, respectively. However, as the discharged pollutants have exceeded its WEC in 2016, there are no more permissible CODMn and NH3-N. Similarly, no more permissible NH3-N have also been found in the river segment S7. Therefore, the pollution discharge should be reduced in the river segments S2 and S7.
To further investigate the impacts of the varied water availability and water demand conditions on the WEC, the percentage changes of WECs are calculated. Schemes I and III, II and IV as well as V and VI are paired to analyze the response to variation of water availability, while the schemes I and II as well as III and IV are paired to investigate the response to different water demand scenarios (shown in Figure 9). Three types of river segments are classified based on their WECs (shown in Figure 10). The river segment where the WEC is slightly responsive to the variations of water availability and water demand is defined as type (i). The river segment S1 is located in the reach between Huangjiagang and Jiangjiaying stations. The river segments S2, S3, S4, S5, S6 are taken as type (ii) where their WECs are sensitive to the changing water demand scenarios and not to the variation of water availability. These river segments are located in the reach between Jiangjiaying and Yuekou stations. The river segments S7 and S8, are defined as type (iii) where their WECs are sensitive to both the variations of water availability and water demand, and these two river segments are located between Yuekou and Wuhan stations.

Impacts of Variations of Water Availability on Water Environmental Capacity
The results from the reference schemes V and VI show that the average designed water flow has declined from 260 m 3 /s to 226 m 3 /s. Due to more uneven temporal-spatial distribution of water availability in the future, there will be less available water resources in dry season (from November to next January), which leads to decreasing the designed low water flow. Thus, there will be notable declines of WEC and their percentage changes with range from 5.17 to 8.66% for CODMn and from 5.36 to 8.85% for NH3-N (shown in Figure 9), respectively. More amounts of pollutants are required to be restricted from discharging into river system (from 2190 t/a to 21,572 t/a for CODMn and from 3304 t/a to 3678

Impacts of Variations of Water Availability on Water Environmental Capacity
The results from the reference schemes V and VI show that the average designed water flow has declined from 260 m 3 /s to 226 m 3 /s. Due to more uneven temporal-spatial distribution of water availability in the future, there will be less available water resources in dry season (from November to next January), which leads to decreasing the designed low water flow. Thus, there will be notable declines of WEC and their percentage changes with range from 5.17 to 8.66% for COD Mn and from 5.36 to 8.85% for NH 3 -N (shown in Figure 9), respectively. More amounts of pollutants are required to be restricted from discharging into river system (from 2190 t/a to 21,572 t/a for COD Mn and from 3304 t/a to 3678 t/a shown in Table 9). Hence, WECs in every river segment from the schemes V and VI will decrease under future water availability scenario.
Water resources allocation for water demand can alter water availability. The WEC varies with the water availability conditions. In the river segments of type (i) and type (ii), WECs are resilient to the varied water availability conditions. The river segments S1, S2, S3, and S4 are located in the downstream of Danjiangkou reservoir while the river segments S5 and S6 are located in the downstream of Changjiang-Hanjiang Water Diversion Project (shown in Figures 3 and 10). There are abundant water resources and strong regulating capability of water projects in S1-S6 river segments. The negative impacts from the uneven temporal-spatial distribution of available water resources on low water flow series are effectively alleviated by water projects operation. Little differences can be found in the designed low water flows of these six river segments under the history and future water availability scenarios, ranging from 0 m 3 /s to 18 m 3 /s (listed in Table 7). Hence, there are little differences of WEC change percentages between schemes I and III, schemes II and IV. Even the maximum differences are even less than 5% (i.e., 3.44% and 4.32% for COD Mn and 3.41% and 4.28% for NH 3 -N according to Figure 9). As the current pollution loads scenario in 2016 is simultaneously applied to the six schemes, there is no increasing pollutants discharge and little differences (less than 5%) can be found in their corresponding permissible pollutants.
The river segments S7 and S8 of the type (iii) are located in the downstream of the main stream (shown in Figure 10), that is far away from Danjiangkou reservoir. There is only one middle size reservoir (i.e., median) for regulating water flow (shown in Figure 2). The water deficit will become more serious and will decrease low water flow series, especially in dry season, due to the more uneven temporal-spatial distribution of available water resources in the future. The designed low water flows will subsequently decrease. Thus, the WECs in S7 and S8 river segments will decrease under future water availability scenarios. The maximum decreases of WECs have even reached 12.19% for COD Mn and 12.04% for NH 3 -N in the river segments S7 and S8 of the type (iii) area, respectively (shown in Figure 9).

Impacts of Water Demand Scenarios on Water Environmental Capacity
The water sources for the cities Xiangyang (operational zones A1, A2, C2 and C6) and Jingmen (operational zones A3, A4, C7 and C16) are mainly from the river segments S2, S3, S4, S5 and S6 (shown in Figure 3). The river segments S7 and S8 are the primary water sources for the cities Xiantao (operational zones A9, A10, C10 and C11) and Wuhan (operational zones A11, A12, A13, C1 and C9) (shown in Figure 3). Their water demands in 2030 will substantially increase with the developments of industry and increasing population. More water will be extracted from these seven river segments (S2-S8) to satisfy their increasing water demands in the future. The designed low water flow series from the schemes III and IV with water demand in future scenarios are obviously lower than that from the schemes I and II. Even the minimum decreases in these river segments are more than 118 m 3 /s and 123 m 3 /s under schemes III and IV, respectively. Thus, their WECs have remarkably decreased (shown in Figure 9), the percentage decreases of which range from 12.78 to 47.77% for COD Mn and from 13.27 to 47.20% for NH 3 -N. The permissible pollutants discharge is then going to decrease if the amount of pollutants discharge keeps the current level. Therefore, the water demand scenarios in the future will increase water deficits and decrease WECs.
However, there are no significant differences of WEC between the water demand scenarios in 2016 and 2030 in the river segment S1. The water demand along the river segment S1 (shown in Figure 3) occupies a very-small percentage of their available water resources, and the increasing amount of water demand can be satisfied through the Danjiangkou reservoir operation. Therefore, the low water flow series in the river segment S1 decreases slightly and there is no significant difference of the designed low water flows (less than 5 m 3 /s, shown in Table 7). The percentage changes of WECs for both COD Mn and NH 3 -N are less than 1.5%, and the impacts of increasing water demand on the permissible pollutants are slight with the amount of pollution loads at current level. Therefore, the increasing water demand does not guarantee the decrease of WEC, as the negative impacts can be effectively alleviated by reservoirs operation.

Impacts of Water Resources Allocation on Water Environmental Capacity
Due to the uneven temporal distribution of water resources, the available water resources in non-flood season only account for 25% and 20% under the history and future water availability scenarios, respectively (listed in Table 4). The designed low water flows under the reference schemes are found much smaller than those under water resources allocation schemes. The average designed low water flows of eight river segments under schemes V and VI are 260 m 3 /s and 226 m 3 /s, while the average designed low water flows under schemes I, II, III and IV are 634 m 3 /s, 519 m 3 /s, 614 m 3 /s and 505 m 3 /s, respectively. Even the minimum difference of the average designed low water flow is 245 m 3 /s. As water resources allocation model can effectively regulate water flow through reservoirs operation, much water is stored in reservoirs in flood season and then is released to satisfy the water demand in dry season. Water demand are effectively assured. The low water flow series is notably increased and further results in the remarkable increase of designed low water flows. The increasing difference of the designed low water flows ranges from 230 m 3 /s (in the segment S5 under scheme II) to 418 m 3 /s (in the segment S7 under scheme I). Thus, WECs increase significantly. The increases of the total WEC under schemes I and II are 244,123 t/a, 130,905 t/a for COD Mn and 4757 t/a, 2595 t/a for NH 3 -N, comparing with those under the reference scheme V (listed in Table 8). Similar results are also found between schemes III, IV and scheme VI. The increasing total WECs are 241,488 and 134,459 t/a for COD Mn and 4726 and 2680 t/a for NH 3 -N, respectively. As the current pollution discharges under the schemes V and VI exceed the WECs, it is required to reduce by 2910 and 21,572 t/a for COD Mn , and 3304 and 3678 t/a for NH 3 -N. While the reservoirs operation in water resources allocation increases WEC feasibly, the exceeding parts of pollutants are significantly decreased under schemes I, II, III and IV. The reduction of the total amount of pollutants is only required for NH 3 -N under schemes II and IV (i.e., 709 t/a and 998 t/a, respectively). Therefore, water resources allocation can regulate water flows to alleviate the negative impacts of uneven temporal-spatial distributions of water resources and can assure the water supply even in dry season. The designed low water flows will increase significantly and will further result in the remarkable increase of WEC.

Conclusions
As the low water flow series is disturbed by climate change as well as reservoirs operation and water acquisition, the stationarity assumption of low water flow series is disturbed and will probably bring risk to WEC planning. In order to figure out the impacts of climate change and water resources allocation in the form of water acquisition and reservoir operation on WEC, a useful framework incorporating water availability estimation module, water demand projection module, water resources allocation module and water environmental capacity determination module was proposed, which offers an effective tool for estimating WEC under nonstationary water flow conditions. The mid-lower reaches of the Hanjiang River basin in China were taken as a case study.
The more uneven temporal-spatial distribution of available water resources due to climate change will bring considerable challenges to WEC. The total available water resources will increase over the study area. However, its uneven temporal-spatial distribution as well as the increasing water demand in the future will bring more stress on water supply, especially in the area with large amount of water demand and weak local storage capacity of the water projects. The low water flow will decrease, which will result in the decreases of their WECs. However, in the areas with the adequate regulating capacities of water projects, little differences of WEC can be found under the varied water availability and water demand conditions. As the water flow is regulated by the water resources allocation model through the reservoirs operation, the uneven temporal-spatial distribution of the water resources is remarkably relieved. More water is stored to satisfy the water demand in the dry season through the reservoirs operation, significantly increasing the low water flow series to further increase the designed low water flow and also the WEC in these areas.
The uncertainties of the results from the proposed framework are not neglected. The uncertainties are not only from the modules themselves, but also from the processes in the connected modules. Specifically, the uncertainties from the climate change scenarios propagate into the calibrated hydrological model in the water availability estimation module. The uncertainties from the water demand projection module and the water resources allocation module can also impact the water environmental capacity determination. As the available water resources are the main driver to the water resources allocation model for determining low water flow for WEC, more attention should be paid to the water availability estimation module in the uncertainty analysis. However, the proposed framework and our research results will not only offer useful guidelines for local water resources management but also demonstrate the potential for effective application in other basins.