# A Provincial Initial Water Rights Incentive Allocation Model with Total Pollutant Discharge Control

^{1}

^{2}

^{3}

^{*}

Next Article in Journal

Next Article in Special Issue

Next Article in Special Issue

Previous Article in Journal

Previous Article in Special Issue

Previous Article in Special Issue

School of Business Administration, Hohai University, Changzhou 213022, China

School of Business, Jiangsu University of Technology, Changzhou 213001, China

School of Business, Hohai University, Nanjing 211100, China

Author to whom correspondence should be addressed.

Academic Editors: Guangxin Zhang and Y. Jun Xu

Received: 25 August 2016 / Revised: 31 October 2016 / Accepted: 2 November 2016 / Published: 11 November 2016

(This article belongs to the Special Issue Tackling Complex Water Problems in China under Changing Environment)

According to the newly introduced water resources management in China, which is the most stringent thus far, pollution limits of water functional zones should be enhanced. In response to this recent change, we propose a provincial initial water rights incentive allocation model with total pollutant discharge constraint. Firstly, three objective functions are set up, namely the maximum comprehensive economic benefits of basin, the optimal fairness and coordination of the provincial pollutant discharge rights, and the minimum losses of basin ecological environment. Based on the objective functions and constraint of total pollutant discharge, we put forward a multiple-objective provincial initial pollutant discharge rights allocation model with a self-adaptive chaotic optimization algorithm. Secondly, according to the incentive mechanism of rewarding excellence and punishing inferiority, we establish a provincial initial water rights allocation model with an incentive function which includes the amount of provincial initial pollutant discharge. Finally, we select the Taihu Basin as a case study to describe the methodology. Under water frequencies of 50%, 75%, and 90%, the empirical results of provincial initial water rights allocation of the Taihu Basin in 2030 show that Jiangsu Province will obtain the most initial water rights, followed by Zhejiang Province and Shanghai city. The incentive allocation model of the provincial initial water rights is effective and the algorithm is feasible. The model may help to assure the most stringent water resources management strategies that are effectively carried out in the Taihu Basin.

Although China’s water resources are rich, the country’s available water is less than one quarter of the world average, in addition, water use efficiency and benefits are relatively low. At present, along with a rapidly growing population, industrialization, and urbanization, two thirds of China’s cities suffer water shortages, and many farmlands are damaged annually by drought. Water pollution is increasing, and many aquatic ecosystems are in decline. In China, water resources now present the most important obstacle for continued economic growth and sustainable social development.

Clarification of initial water rights is the basis and prerequisite for water rights trade-offs in China. Hence, it is an effective solution to the problem of water shortages and an important approach to enhance the efficiency and benefit of water resources utilization. As an important part of basin initial water rights, provincial initial water rights include provincial initial water quantity rights and provincial initial pollutant discharge rights [1]. The allocation of provincial initial water rights is a complicated problem that involves redistribution of vested interests and needs to give consideration to the requirements of multiple participants. Facing the most stringent water resources management, the theory and practice of provincial initial water rights allocation should meet the new requirements of this institution. In addition, a reasonable allocation of provincial initial water rights, which is tied to the total water use, the water use efficiency and pollutant discharge have become the major elements and core aspects of basin initial water rights allocation during the new period of the most stringent water resources management in China.

The optimal allocation of provincial initial water rights has been an active research field both in foreign and domestic areas in recent years. Preliminary literature were concerned with quantity allocation of provincial initial water rights including hybrid allocation methods [2,3], multi-participation allocation models [3,4,5], multi-objective optimization methods [6,7,8,9], and interactive or harmonious allocation models [10,11,12]. Afterwards, considering the importance of water quality, some research paid attention to provincial initial pollutant discharge rights allocation, which can be divided into the following four aspects:

- (1)
- (2)
- Using a single-objective decision making model, Gao Zhu et al. [15] carried out research on the initial allocation of pollutant discharge based on water function regionalization. Wan Shan et al. [16] established an initial pollutant discharge allocation model based on the simulation study of both economic optimality and fairness of allocation.
- (3)
- With a multi-objective decision making method, Huang Xianfeng et al. [17] took into consideration a two level allocation model of river emission rights. The economical optimization and water quality optimization were regarded as the objective functions. The concentration and total quantity of pollutant were the constraints and the principle of equality and continuity of the production program of the dischargers were taken into account. In addition, Huang Binbin et al. [18] carried out research on an optimal allocation model of river pollutant discharge based on pollution limits.
- (4)
- Using a hybrid allocation model, Liu Gang et al. [19] constructed a cooperation allocating initial discharge permits system for industrial source points in the Taihu Lake basin which was regulated by the policy-oriented government, led by economy-oriented government subjects and included the participation of many stakeholders. In addition, Liu Nianlei et al. [20] studied the pollutant discharge with an entropy weighing method and improved the proportional distribution method.

The literature mentioned above is only concerned with the single perspective of water quantity or water quality. Since water quantity and water quality are two basic properties of water resources, initial water rights allocation should consider water quantity and quality simultaneously. Recently, several studies have focused on this area, but the research is still insufficient. For instance, Wang Zongzhi et al. [21] established a basin initial two-dimensional water rights allocation model in which the maximum water resource system in harmony with a river basin is regarded as an object function. However, they did not take into account the initial pollutant discharge allocation aiming at the most stringent water resources management constraints. They only established a penalty function of a higher-than-standard pollutant discharge, but ignored positive inspiration for the lower-than-standard pollutant discharge. Wu Dan et al. [22] applied the interactive iteration algorithm into the study of basin initial water rights allocation. However, they neglected incentive mechanisms of rewarding excellence and punishing inferiority in the water environmental protection. In addition, Zhang Lina et al. [23] made an analysis of the framework of basin initial water rights allocation from a coupling perspective. However, they did not carry out the specific model of this framework.

According to the existing literature, despite different perspectives, many studies have applied a considerable amount of provincial water rights allocation methods in pursuit of a more scientific and effective allocation plan. However, for the new scenarios of the most stringent water resources management, the existing studies of relevant aspects have the following defects:

- Firstly, the most stringent water resources management set the benchmark of the “three red lines”, which means a “red line” for total water use, a “red line” for water use efficiency, and a “red line” for controlling pollution. The initial water rights allocation models focused on the first two “red lines”, but were less concerned with the third “red line” which illustrates the pollution limits of the water functional zones that should be enhanced. The third “red line” can be considered by the allocation of initial pollutant discharge, but current models of provincial initial pollutant discharge rights allocation are insufficient, and existing research in this area also lacks a comprehensive consideration of utilization benefits including economy benefit, society benefit, and ecological environmental benefit.
- Secondly, the incentive mechanism of rewarding excellence and punishing inferiority is crucial. In the field of water resources management, it is important to punish the over standard pollutant discharge during environmental protection of water resources. However, in reality, it is also necessary to implement positive incentives to the less than standard pollutant discharge. How to embed water pollution problems into the initial water right allocation and establish the incentive function is rarely considered and really needs further discussion.

To address these issues, in this paper we apply a well-developed model that produces new insights regarding the related literature. To be specific, under the new scenario of the most stringent water resources management, according to the situation where a water shortage caused by the poor quality of water has become a bottleneck during the sustainable development of a provincial economy, we propose a provincial initial water right incentive allocation model with a total pollutant discharge constraint. More concretely, the contributions of this paper can be illustrated as follows:

- Firstly, in contrast to the existing models of provincial initial pollutant discharge rights allocation, we focus on the comprehensive consideration of utilization benefits including economic benefits, societal benefits, and ecological environmental benefits. Hence, the first contribution of this paper is that we set up three objective functions separately including the maximum comprehensive economic benefits of basin, the optimal fairness and coordination of provincial pollutant discharge rights, and the minimum losses to the ecological environment of the basin. Furthermore, based on the three objective functions and the constraints of the total pollutant discharge, we put forward a multiple-objective provincial initial pollutant discharge rights allocation model with a self-adaptive chaotic optimization algorithm.
- Secondly, none of the aforementioned papers consider the incentive mechanism for rewarding excellence and punishing inferiority in the water environmental protection; however, awarding the less than standard pollutant discharge is as important as punishing the over-standard pollutant discharge. Hence, the second contribution of this paper is that we establish an incentive function by embedding the amount of provincial initial pollutant discharge and we follow the mechanism of rewarding excellence and punishing inferiority.
- Thirdly, by using this incentive function, we attempt to embed pollutant discharge allocation into the process of initial water quantity rights allocation, so as to internalize externalities of over-standard or under-standard pollutant discharge into the initial water rights allocation. Hence, the third contribution of this paper is that we set up a provincial initial water rights allocation model with the above incentive function considering both water quantity and water quality.

The rest of the paper is organized as follows. Section 2 introduces the framework of our research. Section 3 presents a multi-objective allocation model of provincial initial pollutant discharge rights with a total pollutant discharge constraint. Section 4 provides an incentive function to build a provincial initial water rights allocation model from a coupling perspective with the dual control of water quantity and water quality. Section 5 displays the case study of the Taihu Basin. Section 6 summarizes our work.

The framework of our research will follow these basic ideas: Creation of a provincial pollutant discharge rights allocation model based on three optimal objectives including benefits to the economy, society, and ecological environment under total pollutant discharge control; a coupled provincial initial water rights incentive allocation model with dual control for water quantity and water quality; and empirical research for the Taihu Basin. It is depicted in Figure 1.

As is clearly indicated in Figure 1, the specific framework includes three steps: (1) building the provincial pollutant discharge rights allocation model. According to the requirement of total pollutant control, we adopt three optimization goals for economy, society, and ecological environment to establish a provincial initial pollutant discharge allocation model and obtain an initial pollutant discharge rights solution for different provinces; (2) developing a coupled provincial initial water rights allocation model. With total pollutant discharge control, we set up an incentive function according to the mechanism of rewarding excellence and punishing inferiority. As the incentive allocation model of provincial initial water rights with water quantity and quality is proposed, we can get the provincial initial water rights allocation results; (3) illustrating empirical research for the Taihu Basin.

We divide the comprehensive utilization benefits into those for economy, society, and the ecological environment. This is mainly based on the following considerations: Firstly, economic benefits reflect the utilization level of water resources and the production efficiency. Secondly, social benefits show the equality of social allocation and the harmonious relationship in a region; Thirdly, ecological benefits mainly represent pressure or maintenance because the pollutant discharge into the river affects the ecological system. Based on these, considering the three optimization objectives of economic, social, and ecological environmental benefits, we put forward a provincial initial pollutant discharge rights allocation model to meet the total pollutant discharge control.

${S}_{k}$ denotes the province, $k=1,2,...,q$. ${x}_{{S}_{k}}$ denotes the pollutant discharge amount of the province ${S}_{k}$, $k=1,2,...,q$.

Objective function 1: let economic benefit optimization be ${F}_{1}$ which means the biggest general economic benefit. This can be presented by Equation (1).
where ${F}_{1k}={f}_{1k}\left(X\right)=GD{P}_{k}\left({x}_{{S}_{k}}\right)$, $GD{P}_{k}\left({x}_{{S}_{k}}\right)$ is the revenue function of province ${S}_{k}$ when it obtains the pollutant discharge rights ${x}_{{S}_{k}}$.

$${F}_{1}=\mathrm{max}{f}_{1}(X)=\mathrm{max}{\displaystyle \sum _{k=1}^{q}{F}_{1k}}=\mathrm{max}{\displaystyle \sum _{k=1}^{q}{f}_{1k}\left(X\right)=}\mathrm{max}{\displaystyle \sum _{k=1}^{q}GD{P}_{k}({x}_{{S}_{k}})}$$

Objective function 2: let social benefits optimization objective function be ${F}_{2}$ which reflects the best equality and coordination of provincial pollutant discharge allocation. This can be described by Equation (2).
where ${F}_{2k}={f}_{2k}(X)=[({x}_{Si}/{x}_{Sj})-{\displaystyle \sum _{l}^{3}{\beta}_{l}{\gamma}_{l}({S}_{i}/{S}_{j})}{]}^{2}$, ${\gamma}_{l}({S}_{i}/{S}_{j})$ separately represents the ratio of population, area, and GDP index value between province ${S}_{i}$ and province ${S}_{j}$, $l=1,2,3$. ${\beta}_{l}$ reflects the relative importance of population, area, and GDP index value between province ${S}_{i}$ and province ${S}_{j}$, and $\sum _{l=1}^{3}{\beta}_{l}}=1$.

$${F}_{2}=\mathrm{min}{f}_{2}\left(X\right)=\mathrm{min}{\displaystyle \sum _{k=1}^{q}{F}_{2k}}=\mathrm{min}{\displaystyle \sum _{k=1}^{q}{f}_{2k}(X)=}\mathrm{min}{\displaystyle \sum _{i=1}^{q}{\displaystyle \sum _{j=1}^{q}[({x}_{Si}/{x}_{Sj})-{\displaystyle \sum _{l}^{3}{\beta}_{l}{\gamma}_{l}({S}_{i}/{S}_{j})}{]}^{2}}}\text{},i\ne j$$

Objective function 3: set ecological environment optimization objective function to ${F}_{3}$ which means the smallest damage to the ecological environment. This can be illustrated by Equation (3).
where ${F}_{3k}={f}_{3k}\left(X\right)={\eta}_{k}{P}_{\mathrm{k}}\left({x}_{{S}_{k}}\right)$, ${\eta}_{k}$ is the environmental damage coefficient of province ${S}_{k}$ after receiving the pollution, $k=1,2,...,q$. When ${\eta}_{k}$ is bigger, the damage situation of province ${S}_{k}$ is worse. ${P}_{\mathrm{k}}\left({x}_{{S}_{k}}\right)$ is the environmental damage function of province ${S}_{k}$ after receiving pollution ${x}_{{S}_{k}}$, $k=1,2,...,q$.

$${F}_{3}=\mathrm{min}{f}_{3}(X)=\mathrm{min}{\displaystyle \sum _{k=1}^{q}{F}_{3k}}=\mathrm{min}{\displaystyle \sum _{k=1}^{q}{f}_{3k}(X)=}\mathrm{min}{\displaystyle \sum _{k=1}^{q}{\eta}_{k}}{P}_{k}({x}_{{S}_{k}})$$

Constraint condition 1: the total amount of every province’s pollutant discharge cannot exceed the amount of the basin’s pollutant discharge ${P}_{total}$, which can be expressed by Equation (4).

$$\sum _{k=1}^{q}\left({x}_{{S}_{k}}\right)\le {P}_{total}$$

Constraint condition 2: non-negative constraints of variables which can be described by Equation (5).
where $k=1,2,...,q$.

$${x}_{{S}_{k}}\ge 0$$

We construct our model using Equations (1)–(5). The model includes the three important objective functions of economic benefits, social benefits and ecological environmental benefits. The detailed procedures of calculation can be described as follows.

Objective function ${F}_{1}$ should be turned into function minimization as follows:

$${F}_{1}=\mathrm{max}{f}_{1}(X)=\mathrm{min}[-{f}_{1}(X)]=\mathrm{min}[-{\displaystyle \sum _{k=1}^{q}GD{P}_{k}({x}_{{S}_{k}})]}$$

However, the objective function ${F}_{1}$, ${F}_{2}$, and ${F}_{3}$ have different scales; therefore, it is necessary to transform them into non-dimensional ones which can be represented by ${F}_{1}^{\prime}$, ${F}_{2}^{\prime}$, and ${F}_{3}^{\prime}$. Here, we use the ideal point method of multiple objective planning to remove dimensions or balance magnitude differences.

Firstly, three single-objective planning models can be expressed by Equations (7)–(9).

$$\{\begin{array}{c}{F}_{1}=\mathrm{min}[-{f}_{1}(X\left)\right]\\ S.T.\{\begin{array}{c}{\displaystyle \sum _{k=1}^{q}\left({x}_{{S}_{k}}\right)\le {P}_{\mathrm{t}otal}}\\ {x}_{{S}_{k}}\ge 0\end{array}\end{array}$$

$$\{\begin{array}{c}{F}_{2}=\mathrm{min}\left[{f}_{2}\right(X\left)\right]\\ S.T.\{\begin{array}{c}{\displaystyle \sum _{k=1}^{q}\left({x}_{{S}_{k}}\right)\le {P}_{\mathrm{t}otal}}\\ {x}_{{S}_{k}}\ge 0\end{array}\end{array}$$

$$\{\begin{array}{c}{F}_{3}=\mathrm{min}\left[{f}_{3}\right(X\left)\right]\\ S.T.\{\begin{array}{c}{\displaystyle \sum _{k=1}^{q}\left({x}_{{S}_{k}}\right)\le {P}_{\mathrm{t}otal}}\\ {x}_{{S}_{k}}\ge 0\end{array}\end{array}$$

Secondly, the optimal solutions of Equations (7)–(9) can be expressed by $-{f}_{1}\left({X}^{*}\right)$ ${f}_{2}({X}^{*})$ and ${f}_{3}({X}^{*})$. In order to remove the dimensions of the objective function ${F}_{1}$, ${F}_{2}$, and ${F}_{3}$, we define ${F}_{1}^{\prime}$ as the non-dimensional objective function of ${F}_{1}$, define ${F}_{2}^{\prime}$ as the non-dimensional objective function of ${F}_{2}$, and define ${F}_{3}^{\prime}$ as the non-dimensional objective function of ${F}_{3}$. Then, according to the ideal point method of multiple objective planning to remove dimensions, ${F}_{1}^{\prime}$, ${F}_{2}^{\prime}$ and ${F}_{3}^{\prime}$ can be illustrated as follows:

$${F}_{1}^{\prime}=\left[{f}_{1}\right(X)-{f}_{1}({X}^{*}\left)\right]/{f}_{1}\left({X}^{*}\right)$$

$${F}_{2}^{\prime}=\left[{f}_{2}\right(X)-{f}_{2}({X}^{*}\left)\right]/{f}_{2}\left({X}^{*}\right)$$

$${F}_{3}^{\prime}=\left[{f}_{3}\right(X)-{f}_{3}({X}^{*}\left)\right]/{f}_{3}\left({X}^{*}\right)$$

Simplify the vector optimization problem of the three objective functions ${F}_{1}$, ${F}_{2}$, and ${F}_{3}$ into scalar optimization and it can be written as follows:

$${F}^{\prime}={\left({F}_{1}^{\prime}\right)}^{2}+{\left({F}_{2}^{\prime}\right)}^{2}+{\left({F}_{3}^{\prime}\right)}^{2}$$

Transfer the former multi-objective planning model into a nonlinear problem which can be expressed by Equation (14).

$$\{\begin{array}{c}\mathrm{min}{F}^{\prime}={({F}_{1}^{\prime})}^{2}+{({F}_{2}^{\prime})}^{2}+{({F}_{3}^{\prime})}^{2}\\ S.T.\{\begin{array}{c}{\displaystyle \sum _{k=1}^{q}({x}_{{S}_{k}})\le {P}_{total}}\\ {x}_{{S}_{k}}\ge 0\end{array}\end{array}$$

In order to obtain the amount of initial pollutant discharge rights ${x}_{{S}_{k}}$ for the different province ${S}_{k}$, $k=1,2,...,q$, we need to solve Equation (14) which is a single-objective nonlinear planning equation. The regular methods to solve Equation (14) can be adopted, including the chaotic optimization algorithm (COA), genetic algorithm (GA), and simulated annealing algorithm. As a new optimization technique, COA is developed with the characteristics of sensibility, pseudo-randomness, ergodicity, and self-similar fractals, but it needs to calculate many objective functions and its disadvantage it its slow convergence [24]. Although GA is an easy solution, it still has many disadvantages including precocious, inferior local search ability. The simulated annealing algorithm has a better local search capability and a shorter running time, but its global search capability is poor, and it is easily affected by parameters.

To overcome the problems of the above algorithms, we adopt the following method: (1) improve the quality of individuals of the initial group through chaotic ergodicity; (2) according to the adaptability of crossing operators and mutation operators, the search capability and the convergence rate of the genetic algorism can be improved. Furthermore, the algorithm makes it almost impossible to fall into local optimization; (3) the whole annealing selection model with the parent involved competition is used to increase diversity. In order to avoid premature convergence, we choose a hybrid parent body Therefore, we adopt a self-adoptive chaotic optimization algorithm [25] to solve Equation (14). The detailed steps of the self-adoptive chaotic optimization algorithm are as follows:

Step 1: Parameter initializing. Set $B$ to the population size, set $M$ to the chaotic iteration, set ${T}_{0}$ to the initial temperature, set ${T}_{END}$ to the final temperature, set ${P}_{c1}$, ${P}_{c2}$, ${P}_{m1}$, and ${P}_{m2}$ to the self-adaptive parameters.

Step 2: Initial population generation. Chaotic equation is produced by the logistic map formula in Equation (15).
where $\vartheta $ is the control parameter, then when $\vartheta =4$ the system is in chaos. According to Equation (15), we can randomly generate several initial populations formed by chaotic series.

$${X}_{k+1,j}=\vartheta \cdot {X}_{k,j}\cdot (1-{X}_{k,j})$$

Step 3: Evaluation of the individual fitness function. Make the objective function in Equation (15) be the self-adaption evaluate function, and punish the individuals who break the constraint conditions. Record the best individual of every generation, after solving it, replace the individual whose fitness is the lowest. Then the best individual is ensured to enter into the next generation.

Step 4: Calculation. According to the whole annealing mechanism [26], allow the parent to take part in the competition. According to the fitness of the individual, the formula of self-adaption, and the formula of selection probability [25], we can determine the crossover probability and the mutation rate, then process the crossover calculation and mutation calculation. Finally, we can get the next generation by this formula of selection probability.

Step 5: Judgment of the termination of algorithm. The descent rate of temperature ${T}_{k}$ can be presented by Equation (16).
where $k$ is iteration time. When temperature ${T}_{k}$ is bigger than the final temperature ${T}_{END}$, go back to Step 3, and evaluate fitness of the individual. Or, we can eventually obtain the amount of initial pollutant discharge rights ${x}_{{S}_{k}}$ for the different province ${S}_{k}$, where $k=1,2,...,q$.

$${T}_{k}=1/\mathrm{ln}(k/{T}_{0}+1)$$

Assume that the total amount of initial allocation water rights is W, then the pre-configuration allocation plans of provincial initial water quantity rights ${w}_{{S}_{k}}$ can be presented by Equation (17).
where ${\omega}_{{S}_{k}}$ is the pre-configuration allocation proportion of the initial water quantity rights of province ${S}_{k}$; in this paper we assume this is already known. Then water quantity rights ${w}_{{S}_{k}}$ need to be coupled with provincial initial pollutant discharge rights ${x}_{{S}_{k}}$ which were calculated by Equation (14). In order to achieve this goal, we absorb the idea from “rewarding excellence and punishing inferiority” to design the incentive mechanism which can be described as follows:

$${w}_{{S}_{k}}=W\times {\omega}_{{S}_{k}}$$

- Set the real amount of pollutant discharge of province ${S}_{k}$ to ${{x}_{{S}_{k}}}^{R}$. When the real amount of pollutant discharge ${{x}_{{S}_{k}}}^{R}$ exceeds the initial pollutant discharge rights ${x}_{{S}_{k}}$, the negative incentive should be implemented to reduce its initial water quantity rights.
- When the real provincial amount of pollutant discharge ${{x}_{{S}_{k}}}^{R}$ is lower than the initial pollutant discharge rights ${x}_{{S}_{k}}$, the positive incentive should be implemented to increase its initial water quantity rights.

According the above incentive mechanism, different provinces in the same basin can get their own initial water rights through the coupling control of total water use and total pollutant discharge. The specific processes can be presented as follows:

Step 1: Set the incentive function to $\mu ({{x}_{{S}_{k}}}^{R}/{x}_{{S}_{k}})$, and then the incentive function can be designed by Equation (18).
where $\partial $ is an adjustment coefficient which reflects the degree of incentive $\partial \in [2,10]$. When the value of $\partial $ is higher, the changing trend of incentive degree is smoother.

$$\mu ({{x}_{{S}_{k}}}^{R}/{x}_{{S}_{k}})=\{\begin{array}{cc}{\left|{x}_{{S}_{k}}/{{x}_{{S}_{k}}}^{R}-1\right|}^{\partial}& {x}_{{S}_{k}}<{{x}_{{S}_{k}}}^{R}\\ 1& {x}_{{S}_{k}}={{x}_{{S}_{k}}}^{R}\\ {\left|{{x}_{{S}_{k}}}^{R}/{x}_{{S}_{k}}-1\right|}^{\partial}& {x}_{{S}_{k}}>{{x}_{{S}_{k}}}^{R}\end{array}$$

Step 2: According to the incentive function, the adjustment of the proportion of the water rights allocation of province ${S}_{k}$ can be described by Equation (19).
where ${{\omega}_{{S}_{k}}}^{\prime}$ is the adjusted allocation proportion of the initial water rights of province ${S}_{k}$.

$${{\omega}_{{S}_{k}}}^{\prime}=\{\begin{array}{cc}{\omega}_{{S}_{k}}(1-\mu ({{x}_{{S}_{k}}}^{R}/{x}_{{S}_{k}}\left)\right)& {x}_{{S}_{k}}<{{x}_{{S}_{k}}}^{R}\\ {\omega}_{{S}_{k}}& {x}_{{S}_{k}}={{x}_{{S}_{k}}}^{R}\\ {\omega}_{{S}_{k}}(1+\mu ({{x}_{{S}_{k}}}^{R}/{x}_{{S}_{k}}\left)\right)& {x}_{{S}_{k}}>{{x}_{{S}_{k}}}^{R}\end{array}$$

Step 3: ${{\omega}_{{S}_{k}}}^{\prime}$ can be normalized by Equation (20).
where $k=1,2,...,q$ and ${{\omega}_{{S}_{k}}}^{\u2033}$ is the normalized allocation proportion of initial water rights of province ${S}_{k}$.

$${{\omega}_{{S}_{k}}}^{\u2033}={{\omega}_{{S}_{k}}}^{\prime}/{\displaystyle \sum _{k=1}^{q}{{\omega}_{{S}_{k}}}^{\prime}}$$

Step 4: Finally, we can figure out the final quantity of initial water rights of province ${S}_{k}$ by Equation (21).
where $k=1,2,...,q$.

$${{w}_{{S}_{k}}}^{\prime}=W\times {{\omega}_{{S}_{k}}}^{\prime \prime}$$

The administration areas in the Taihu Basin mainly include Jiangsu province, Zhejiang province, and Shanghai city. The map of the main administrative divisions in Taihu Basin are illustrated in Figure 2. The Taihu Basin is one of the most developed and most modernized regions in China. The amount of annual water use has reached 36.5 billion cubic meters, while the amount of the perennial water resources of the basin is only 19.6 billion cubic meters, so the demand-supply gap is huge. The problems of water shortage, worsening water pollution, and increasing demand for water resources have imposed undue pressure on the carrying capacities of the water resources and water environment in the Taihu Basin.

According to the materials of the Taihu Basin and Southeast Rivers Water Resources Bulletin and the Water Quality Bulletin of Taihu Basin and Important Water Function Areas of Southeast Rivers, and others, the pre-configuration allocation plan of the initial water rights for the Taihu Basin can be set as follows: the ratios of the initial water quantity rights of Jiangsu province, Zhejiang province, and Shanghai city are, respectively, 35.62%, 31.94%, and 32.45%. Furthermore, the area of those three provinces are 19,399 square kilometers, 12,093 square kilometers, and 5178 square kilometers, respectively. The specific relevant data of the Taihu Basin from 2011 to 2014 are shown in Table 1.

According to the materials of the Taihu Basin and Southeast Rivers Water Resources Bulletin in planning year 2030, the total amount of pollutant discharge rights will be 54.46 hundred million cubic meters. Using the data in Table 1, according to a multi-objective provincial initial pollutant discharge rights allocation model which is expressed by Equation (14), we can obtain relevant solutions which are shown in Table 2.

In planning year 2030, using the amount of provincial initial pollutant discharge ${x}_{{S}_{k}}$ in Table 2, according to the coupled provincial initial water rights incentive allocation model which is expressed by Equations (18)–(21), we can obtain the normalized allocation ratios of the initial water rights of the provinces which are shown in Table 3.

According to the materials of Water Resources Planning for Taihu Basin (2012–2030), under three different water frequencies of 50%, 75%, and 90%, in the planning year 2030, the amount of initial allocation water rights of Taihu Basin are, respectively, 329.2 hundred million cubic meters, 363.3 hundred million cubic meters, and 392.6 hundred million cubic meters.

Using the Equation (21), in the planning year 2030, the provincial initial water rights allocation plans in the Taihu Basin under water frequencies of 50%, 75%, and 90% can be calculated and are shown in Table 4.

For comparison with the “water use efficiency” allocation method provided by [28], the provincial initial water rights allocation plans are also presented in Table 4. The contrasted result can be shown as follows: (1) Under the different water frequencies of 50%, 75%, and 90%, in planning year 2030, our solutions are approximately consistent with the results of the “water use efficiency” allocation method, namely, Jiangsu Province obtains the most initial water rights, followed by Zhejiang Province and Shanghai City; (2) According to the pollutant discharge, due to our incentive function of rewarding excellence and punishing inferiority, when water frequency is the same, compared with the result of the “water use efficiency” allocation method, our result for Shanghai City is lower, while our results for Jiangsu province and Zhejiang province are higher. Since water resources have natural coupling characteristics of quantity and quality, our incentive allocation model emphasizes that considering both water quantity and quality is important to the provincial initial water rights allocation.

In the context of the newly introduced water resources management in China, which is the most stringent thus far, the major elements and the core aspects of the basin initial water rights allocation are: considering a reasonable allocation of provincial initial water rights tied to the total water use, the water use efficiency and pollutant discharge. This research is focused upon the provincial initial water rights allocation model based on the dual control of water quantity and quality, which supports policy making decisions for the government with respect to sustainable water resources development. It also contributes to the policy implementation of the most stringent water resources management. In this paper, according to the provincial initial water rights incentive allocation model with total pollutant discharge constraint, the empirical result of the Taihu Basin shows that, under the different water frequencies of 50%, 75%, and 90%, in planning year 2030, Jiangsu Province obtains the most initial water rights, followed by Zhejiang Province and Shanghai City. Through the comparison with the “water use efficiency” allocation method, this paper illustrates that the province initial water rights incentive allocation model with the mechanism of rewarding excellence and punishing inferiority is effective and the algorithm is feasible. The incentive allocation model in this paper provides an important reference for similar research in the other basins. Of course, this study still has limitations. How to design a more reasonable incentive function will be our future research focus and requires further discussion.

The authors would like to thank the National Natural Science Foundation of China (Grant No. 71301063; Grant No. 41271537), sponsorship provided by “Jiangsu Overseas Research & Training Program for University Prominent Young & Middle-aged Teachers and Presidents”, and “333” high-caliber talents cultivation project of Jiangsu Province (No. (2013)III-2110) for their financial support.

Min Ge designed the research and drafted the manuscript; Feng-Ping Wu conducted the model simulation, reviewed, commented on, and revised the manuscript; Min You collected the data and participated in the model simulation. All authors have read and approved the final manuscript.

The authors declare no conflict of interest.

- Wang, H.; Dang, L.; Xie, X.; Qi, Y. Theory and Practice on Allocation of Initial Water Rights in a Watershed; China Water Power Press: Beijing, China, 2008. (In Chinese) [Google Scholar]
- Wang, H.; Wang, J.; Qin, D. Research advances and direction on the theory and practice of reasonable water resources allocation. Adv. Water Sci.
**2004**, 15, 123–128. [Google Scholar] - Wu, D.; Wu, F.; Chen, Y. The bi-level optimization model of the compound system for basin initial water right allocation. Syst. Eng. Theory Pract.
**2012**, 32, 196–202. (In Chinese) [Google Scholar] - Ralph, W.A. Modeling river-reservoir system management, water allocation, and supply reliability. J. Hydrol.
**2005**, 300, 100–113. [Google Scholar] - Chen, Y.; Wu, F.; Zhou, Y. Analysis of evolutionary game between strong group and vulnerable group in initial water rights allocation. Soft Sci.
**2011**, 25, 11–15. (In Chinese) [Google Scholar] - Read, L.; Madani, K.; Inanloo, B. Optimality versus stability in water resource allocation. J. Environ. Manag.
**2014**, 133, 343–354. [Google Scholar] [CrossRef] [PubMed] - Wu, F.; Ge, M. Initial allocation model for water right of the first hierarchy. J. Hohai Univ. (Nat. Sci.)
**2005**, 33, 216–219. (In Chinese) [Google Scholar] - Huang, X.; Shao, D.; Gu, W.; Dai, T. Optimal water resources deployment based on multi-objective chaotic optimization algorithm. J. Hydraul. Eng.
**2008**, 39, 183–188. (In Chinese) [Google Scholar] - Condon, L.E.; Maxwell, R.M. Implementation of a linear optimization water allocation algorithm into a fully integrated physical hydrology model. Adv. Water Resour.
**2013**, 60, 135–147. [Google Scholar] [CrossRef] - Wang, Z.J.; Zheng, H.; Wang, X.F. A harmonious water rights allocation model for Shiyang river basin, Gansu province, China. Int. J. Water Resour. Dev.
**2009**, 25, 355–371. (In Chinese) [Google Scholar] - Wu, F.; Ge, M. Method for interactive water right initial allocation based on harmoniousness judgment. J. Hohai Univ. (Nat. Sci.)
**2006**, 34, 104–107. (In Chinese) [Google Scholar] - Wu, F.; Chen, Y. Research on the Harmonious Allocation Method of Initial Water Rights in a Basin; China Water Power Press: Beijing, China, 2010. (In Chinese) [Google Scholar]
- Mostafavi, S.A.; Afshar, A. Waste load allocation using non-dominated archiving multi-colony ant algorithm. Procedia Comput. Sci.
**2011**, 3, 64–69. [Google Scholar] [CrossRef] - Sun, T.; Zhang, H.; Wang, Y. The application of information entropy in basin level water waste permits allocation in China. Resour. Conserv. Recycl.
**2013**, 70, 50–54. (In Chinese) [Google Scholar] [CrossRef] - Gao, Z.; Li, S. Research on the initial allocation of emission permits based on water function regionalization. Shanghai Manag. Sci.
**2010**, 5, 36–38. (In Chinese) [Google Scholar] - Wan, S.; Li, S.; Lin-Jie, M.A. The initial allocation of emission permits in river basin. J. Syst. Manag.
**2013**, 2, 278–281. (In Chinese) [Google Scholar] - Huang, X.; Shao, D.; Gu, W. Multi objective optimal allocation model of river emission rights. J. Hydraul. Eng.
**2008**, 39, 73–78. (In Chinese) [Google Scholar] - Huang, B.; Wang, X.; Hu, Z.; Gui, F.L. Optimal allocation model of river emission rights based on pollution Limits. Resour. Environ. Yangtze Basin
**2011**, 20, 1508–1513. (In Chinese) [Google Scholar] - Liu, G.; Wang, H.; Qiu, L. Construction of a cooperation allocating initial discharge permits system for industrial source points in a lake basin. Resour. Environ. Yangtze Basin
**2012**, 21, 1223–1229. (In Chinese) [Google Scholar] - Liu, N.; Jiang, H.; Lu, Y.; Zhang, J. Study on allocation for total amount controlling Objectives of Water Pollutants—Considering the Constraint of Environmental Goals of National Main Function regions. China Popul. Resour. Environ.
**2014**, 24, 82–89. (In Chinese) [Google Scholar] - Wang, Z.; Hu, S.; Wang, Y. Initial two-dimensional water rights allocation modeling based on water quantity and water quality in the river basin. J. Hydraul. Eng.
**2010**, 41, 524–530. (In Chinese) [Google Scholar] - Wu, D.; Wu, F. The coupling allocation of initial two-dimensional water rights in basin based on the bi-hierarchy optimal model. China Popul. Resour. Environ.
**2012**, 22, 26–34. (In Chinese) [Google Scholar] - Zhang, L.; Wu, F.; Jia, P. A preliminary theoretical framework of basin initial water rights allocation from a coupling perspective and when constrained by the strictest water resources management system. Resour. Sci.
**2014**, 36, 2240–2247. (In Chinese) [Google Scholar] - Liu, Z.; Yang, D. Research advances of chaos optimization algorithms for engineering global optimization. Chin. J. Comput. Mech.
**2016**, 33, 269–286. (In Chinese) [Google Scholar] - Sen, W.; Chun, C.; Xin, W.; Bao, L. Application of self-adaptive chaos whole annealing genetic algorithm to optimal operation of hydropower station groups. J. Hydroelectr. Eng.
**2014**, 33, 63–71. (In Chinese) [Google Scholar] - Zhang, J.S.; Xu, Z.B.; Liang, Y. The whole annealing genetic algorithms and their sufficient and necessary conditions of convergence. Sci. China (Ser. E)
**1997**, 27, 154–164. (In Chinese) [Google Scholar] - Xin, G. Research of Water Environment Ecological Risk Based on Pollution Footprints Theory. Master’s Thesis, Nanjing University of Finance and Economics, Nanjing, China, 2013. [Google Scholar]
- Li, Z.; Feng, W.; Chen, Z. Initial Water Rights Differentiated Allocation among Provinces in Taihu Basin Constrained by Water use Efficiency in Multiple Scenarios. China Popul. Resour. Environ.
**2015**, 25, 124–132. (In Chinese) [Google Scholar]

Province | Year | Population (Ten Thousand Persons) | Per Capita GDP (Yuan) | Pollutant Discharge (One Hundred Million Cubic Meters) |
---|---|---|---|---|

Jiangsu | 2011 | 2939 | 8.21 | 29.0 |

2012 | 2960 | 8.97 | 29.1 | |

2013 | 2985 | 9.84 | 29.2 | |

2014 | 2996 | 10.68 | 28.3 | |

Zhejiang | 2011 | 1763 | 6.54 | 12.1 |

2012 | 1775 | 6.97 | 12.1 | |

2013 | 1791 | 7.53 | 12.6 | |

2014 | 1797 | 7.98 | 12.9 | |

Shanghai | 2011 | 1175 | 12.01 | 22.6 |

2012 | 1184 | 12.53 | 23.1 | |

2013 | 1194 | 13.44 | 22.9 | |

2014 | 1198 | 14.47 | 22.9 |

Notes: Data source: Taihu Basin and Southeast Rivers Water Resources Bulletin (2011–2014) and Water Quality Bulletin of Taihu Basin and Important Water Function Areas of Southeast Rivers (2011–2014).

Provincial Initial Pollutant Discharge Rights Allocation | Calculation Results |
---|---|

Objective function F_{1} | Revenue function of Jiangsu’s pollutant discharge rights ^{1}: ${F}_{11}=GD{P}_{1}\left({x}_{1}\right)=1410.8{x}_{1}-32200$ |

Revenue function of Zhejiang’s pollutant discharge rights: ${F}_{12}=GD{P}_{2}\left({x}_{2}\right)=1200.7{x}_{2}-10586$ | |

Revenue function of Shanghai’s pollutant discharge rights: ${F}_{13}=GD{P}_{3}\left({x}_{3}\right)=4248.8{x}_{3}-80597$ | |

${F}_{1}=\mathrm{max}{f}_{1}(X)=\mathrm{min}[-{f}_{1}(X)]=\mathrm{min}[(32200-1410.8{x}_{1})+(10586-1200.7{x}_{2})+(80597-4248.8{x}_{3})]$ Solving the Equation (7), ${f}_{1}({X}^{*})=16359$; Using Equation (10), ${F}_{1}^{\prime}=\mathrm{min}(0.086{x}_{1}+0.073{x}_{2}+0.26{x}_{3}-8.542)$ | |

Objective function F_{2} | Parameters selected ^{2}: relative importance of province population, area, and GDP index: ${\beta}_{\text{}1}=1/3$, ${\beta}_{\text{}2}=1/3$, ${\beta}_{\text{}3}=1/3$ |

Using Equation (2) to calculate: ${F}_{21}={({x}_{1}/{x}_{2}-{\beta}_{1}\times 1.667-{\beta}_{1}\times 1.604-{\beta}_{1}\times 2.500)}^{2}$; ${F}_{22}={({x}_{1}/{x}_{3}-{\beta}_{1}\times 2.500-{\beta}_{1}\times 3.746-{\beta}_{1}\times 0.625)}^{2}$; ${F}_{23}={({x}_{2}/{x}_{3}-{\beta}_{1}\times 1.500-{\beta}_{1}\times 2.335-{\beta}_{1}\times 0.250)}^{2}$; ${F}_{2}={F}_{21}+{F}_{22}+{F}_{23}=\mathrm{min}[{f}_{2}(X)]=\mathrm{min}[{({x}_{1}/{x}_{2}-1.92)}^{2}+{({x}_{1}/{x}_{3}-2.29)}^{2}+{({x}_{2}/{x}_{3}-1.36)}^{2}]$; Solving the Equation (8), ${f}_{1}({X}^{*})=0.0172$; Using Equation (11), ${F}_{2}^{\prime}=\mathrm{min}[{({x}_{1}/{x}_{2}-1.92)}^{2}+{({x}_{1}/{x}_{3}-2.29)}^{2}+{({x}_{2}/{x}_{3}-1.36)}^{2}-0.0172]/0.0172$ | |

Objective function F_{3} | Parameters selected: environment damage coefficient ^{3}: Jiangsu ${\eta}_{1}$ = 0.6; Zhejiang ${\eta}_{2}$ = 0.5; Shanghai ${\eta}_{3}$ = 0.5 |

Environmental damage function of Jiangsu province ^{4}: ${P}_{1}\left({x}_{{S}_{1}}\right)=3075.6{x}_{1}-47695$; Environmental damage function of Zhejiang province: ${P}_{2}\left({x}_{{S}_{2}}\right)=2890.2{x}_{2}-9226.2$; Environmental damage function of Zhejiang province: ${P}_{3}\left({x}_{{S}_{3}}\right)=1062.9{x}_{3}-13426$; ${F}_{3}=\mathrm{min}[{f}_{3}(X)]=\mathrm{min}[0.6\times (3075.6{x}_{1}-47695)+0.5\times (2890.2{x}_{2}-9226.2)+0.5\times (1062.9{x}_{3}-13426)]$; Solving the Equation (7), ${f}_{3}\left({X}^{*}\right)=12292$; Using Equation (11), ${F}_{3}^{\prime}=\mathrm{min}(0.150{x}_{1}+0.118{x}_{2}+0.043{x}_{3}-4.25)$ | |

Constraint condition | ${x}_{1}+{x}_{2}+{x}_{3}\le 54.46;\text{}{x}_{1}\ge 0$; ${x}_{2}\ge 0$; ${x}_{3}\ge 0$ |

Algorism parameter selecting | Initial population size B = 200, Chaotic iteration M = 100, Initial temperature T_{0} = 100 |

Initial pollutant discharge of province (a hundred million cubic meters) | Using the Equation (14) to obtain the results as follows: Jiangsu ${x}_{{S}_{1}}=27.55$; Zhejiang ${x}_{{S}_{2}}=14.93$; Shanghai ${x}_{{S}_{3}}=11.98$ |

Notes: ^{1} Revenue function of pollutant discharge can be obtained mainly by fitting analysis between the GDP of the province and the total pollutant discharge of the province; ^{2} Considering the population, area, and GDP as equally important to the initial water rights allocation, we assume ${\beta}_{\text{}1}=1/3$, ${\beta}_{\text{}2}=1/3$, ${\beta}_{\text{}3}=1/3$; ^{3} The environmental damage coefficients are obtained by the combination of actual provincial situation and expert opinion approach; ^{4} Environmental damage function can be obtained by fitting analysis between the pollution footprint and the pollutant discharge of the province according to the evaluation model of the water environmental ecological damage [27].

Variables for Initial Water Rights Incentive Allocation | Calculation Results |
---|---|

Parameters Selected | According to the Ratio between Real Pollutant Discharge and Initial Pollutant Discharge and the Expert Opinion Approach, We Select the Adjustment Coefficient $\partial =3$ |

Adjusted value of incentive function for provinces | Using Equation (18) to calculate: Jiangsu Province: $(1-27.55/28.97{)}^{2}=0.002$; Zhejiang Province: $(1-11.97/14.93{)}^{2}=0.039$; Shanghai City: $(1-11.98/22.66{)}^{2}=0.222$ |

Pre-configuration allocation ratio of the initial water quantity rights of province | According to [28], the pre-configuration allocation ratio of the initial water quantity rights of province can be expressed as follows: Jiangsu Province ${\omega}_{{S}_{1}}=35.63\%$; Zhejiang Province ${\omega}_{{S}_{2}}=31.94\%$; Shanghai City ${\omega}_{{S}_{3}}=32.45\%$ |

Adjusted allocation ratio of initial water rights of province | Using Equation (19) to calculate: Jiangsu Province ${{\omega}_{{S}_{1}}}^{\prime}=35.53\%$; Zhejiang Province ${{\omega}_{{S}_{2}}}^{\prime}=33.20\%$; Shanghai City ${{\omega}_{{S}_{3}}}^{\prime}=25.24\%$ |

The normalized allocation ratio of initial water rights of province | Using Equation (20) to calculate: Jiangsu Province ${{\omega}_{{S}_{1}}}^{\u2033}=37.81\%$; Zhejiang Province ${{\omega}_{{S}_{2}}}^{\u2033}=35.33\%$; Shanghai City ${{\omega}_{{S}_{3}}}^{\u2033}=26.86\%$ |

Different Scenarios | Jiangsu Province | Zhejiang Province | Shanghai City | ||||||
---|---|---|---|---|---|---|---|---|---|

Different Water Frequency | Different Water Frequency | Different Water Frequency | |||||||

50% | 75% | 90% | 50% | 75% | 90% | 50% | 75% | 90% | |

Our incentive allocation plans (a hundred million cubic meters) | 124.47 | 137.36 | 148.44 | 116.31 | 128.35 | 138.71 | 88.42 | 97.58 | 105.45 |

“Water use efficiency” allocation model (a hundred million cubic meters) | 117.25 | 129.39 | 139.83 | 105.13 | 116.02 | 125.38 | 106.82 | 117.89 | 127.39 |

© 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).