A Recourse-Based Type-2 Fuzzy Programming Method for Water Pollution Control under Uncertainty

In this study, a recourse-based type-2 fuzzy programming (RTFP) method is developed for supporting water pollution control of basin systems under uncertainty. The RTFP method incorporates type-2 fuzzy programming (TFP) within a two-stage stochastic programming with recourse (TSP) framework to handle uncertainties expressed as type-2 fuzzy sets (i.e., a fuzzy set in which the membership function is also fuzzy) and probability distributions, as well as to reflect the trade-offs between conflicting economic benefits and penalties due to violated policies. The RTFP method is then applied to a real case of water pollution control in the Heshui River Basin (a rural area of China), where chemical oxygen demand (COD), total nitrogen (TN), total phosphorus (TP), and soil loss are selected as major indicators to identify the water pollution control strategies. Solutions of optimal production plans of economic activities under each probabilistic pollutant discharge allowance level and membership grades are obtained. The results are helpful for the authorities in exploring the trade-off between economic objective and pollutant discharge decision-making based on river water pollution control.


Introduction
The trade-off between water pollution control and economic development is of great concern in many basins since it is essential to local sustainable development [1,2].It is difficult to keep the economy growing under the utilization of water resources and the deterioration of environmental conditions [3,4].Meanwhile, it is hard to promote human society improvement if the authorities excessively restrict economic development.Under such a contradictory situation, optimization techniques are proper to detect the economic and environmental impacts of alternative pollution control actions from a system point of view, and thus aid the authorities in formulating and adopting cost-effective water pollution plans and policies.However, water pollution control planning is governed by significant sources of uncertainty associated with different variables, and uncertainty is a non-negligible constituent of such a procedure [5].There are significant uncertainties in not only how the system might develop, but also in how the system is expected to adjust when many system components are altered (e.g., pollutant discharge amount, cost/benefit coefficient and economic activity scale).
In water pollution control problems, values of associated parameters (e.g., cost/benefit coefficients) are usually determined via tests, experiences and expertises, while these methods may fail in determining accurate values, resulting in the parameters being described by fuzzy membership functions.Such deviations in subjective estimations can lead to fuzziness being inherent in the real-world decision problems (e.g., vagueness and/or impreciseness in the outcomes of a water pollution control sample), neglect of which can cause the solutions of problems deviating greatly from their true values.Fuzzy mathematic programming (FMP), based on fuzzy sets theory, can facilitate the analysis of system associated with uncertainties being derived from vagueness and imprecision [6,7].FMP is capable of handling decision problems under fuzzy goal and constraints and tackling ambiguous coefficients in the objective function and constraints.Previously, a wide range of FMP methods were developed for water pollution control [8][9][10][11][12][13][14][15].For example, Liu et al. [16] improved a two-stage fuzzy robust programming model for water pollution control to address fuzzy parameters, which were represented by possibility distributions on the left-and right-hand sides of the constraints.Tavakoli et al. [17] developed an interactive two-stage stochastic fuzzy programming method to handle uncertainties expressed as fuzzy boundary intervals (i.e., the lowerand upper-bounds of intervals are presented as possibility distributions).Ji et al. [18] enhanced an inexact left-hand-side chance-constrained fuzzy multi-objective programming approach to cope with fuzziness in the constraints and objectives.Generally, the conventional fuzzy programming methods could only tackle fuzzy uncertainty with precise membership grades, which may encounter difficulty when the membership grades are also obtained as fuzzy sets.
In many real-world situations, related data such as unit net benefits of economic activities and pollutant discharge allowances are often highly uncertain, which could not be handled by the conventional fuzzy programming methods.When it is challenging to identify the membership grade of a fuzzy set as crisp values (e.g., unit benefit), type-2 fuzzy sets (T2FS) can effectively determine the membership function through defying membership grades of T2FS are fuzzy sets within [0,1] [ [19][20][21][22].In addition, membership functions cannot express uncertainties featured with randomness.Two-stage stochastic programming with recourse is effective for handling decision-making problems in which an examination of policy levels is desired and the system data is characterized by uncertainty [23].Therefore, as an extension to the existing approaches, a recourse-based type-2 fuzzy programming (RTFP) method incorporating the concepts of type-2 fuzzy programming (TFP) within a two-stage stochastic programming with recourse (TSP) framework will be developed to address the above deficiencies.Then, the RTFP method will be applied to water pollution control in the Heshui River Basin in China.Results of optimal agriculture, industry, forestry, fishery, and livestock husbandry activities will be generated, which will be used for providing insight into the trade-off among system benefit, water pollution control, and sustainability.

Methodology
In water pollution control decision making problems, uncertainties may arise due to subjective estimation.For instance, decision makers may estimate the unit benefit from planting fruit/vegetable being [424.1, 575.6] × 10 3 $/km 2 with additional information as the possibility of "the most possible unit benefit of 499.8 × 10 3 $/km 2 " is 0.8, and the possibilities of "there is no possibility that the unit benefit is lower than 424.1 × 10 3 $/km 2 or higher than 575.6 × 10 3 $/km 2 " are 0.2 and 0.3, respectively.Under such a situation, unit benefit from planting fruit/vegetable should be described as type-2 fuzzy sets (T2FS).Thus, type-2 fuzzy programming (TFP) can be adopted to tackle such uncertainties, which can be presented as [24]: subject to: where C ∈ {R} 1×n and A ∈ {R} m×n are vectors of T2FS.A type-2 fuzzy set C in X is a fuzzy set in which the membership function is also fuzzy (i.e., type-2 membership function).The C defined on the universe of discourse X is represented as [25]: where 0 ≤ µ C (x, u) ≤ 1 is the type-2 membership function, J x ⊆ [0, 1] is the primary membership of x ∈ X, which is the domain of the secondary membership function µ C (x) so that all u ∈ J x are the primary membership grades of the point x.The secondary membership function of a triangular C = (c 1 , c 2 , c 3 , θ 1 , θ r ) can be defined as: where c 1 , c 2 , c 3 are real numbers and θ 1 ; θ r ∈ [0, 1] are two parameters representing the spreads of primary membership grades of C. One of the main limitations of the TFP method remains in its difficulty in coping with uncertainties described as probability distributions when the available historical data is sufficient (e.g., pollutant discharge allowances) [26][27][28].Such a problem can be formulated as a two-stage stochastic programming with a recourse (TSP) model [23].Through incorporating the TFP method within the TSP framework, a recourse-based type-2 fuzzy programming (RTFP) method can be formulated as follows: subject to: y js ≥ 0, j = 1, 2, . . ., n 2 ; s = 1, 2, . . ., v.
Decision variables are divided into two subsets: those that must be determined before the realizations of random variables are known, and those (recourse variables) that are determined after the realized random variables are disclosed.x j is the first-stage decision made before the random variable is observed, ω s is the random variable with a probability level p s (i.e., the probability of realization of ω s , with p s > 0 and ∑ S s=1 p s = 1) [29], and y js is the second-stage adaptive decision, which depends on the realization of the random variable.An extra type reduction is needed to convert the output of T2FS into conventional fuzzy sets so that they can be defuzzified to give crisp outputs.Suppose that h 1 , h 2 , . . ., h i are the value of c ij (for at least one pair of i = j, define h i = h j ) evaluated by t different experts.The relative distances of h i are used to approximate the center.Values lying closer to the center are considered more important.Generally, the fuzzification of c ij can be represented as: (1) The conventional fuzzy number F = (a, m, b) of c ij can thus be acquired.Then, the defuzzification of T2FS can be conducted according to the critical value (CV)-based reduction method [30].
Suppose that c is a triangular type-2 fuzzy variable with secondary possibility distribution function µ c (x) (which represents a regular fuzzy variable).The method is introducing the CVs as representing values of the regular fuzzy variable µ c (x), i.e., CV * [µ c (x)] (optimistic CV), CV * [µ c (x)] (pessimistic CV) and CV[µ c (x)].Then, the corresponding fuzzy variables are derived using these CVs of the secondary possibilities: where α 1 , α 2 , and α 3 are primary membership grade of c, and m, 1, and n are corresponding secondary membership grades.
Then, CV * [µ c (x)], CV * [µ c (x)] and CV[µ c (x)] are defined as follows: Finally, to obtain crisp values, centroid method (∑ x xµ c (x)/∑ x µ c (x)) is used for these reduced conventional fuzzy variables.Generally, the detailed computational processes for solving the RTFP method can be summarized as follows: Step 1. Formulate the RTFP model.
Step 2. Discrete probability distribution into several values with each corresponds to one probability.
Step 3. Convert the output of T2FS into conventional fuzzy sets.
Step 4. Conduct defuzzification of T2FS according to the critical value (CV)-based reduction method.
Step 6. Obtain the optimal solutions of the objective function ( f ), first-stage decision variable (x j ), and second-stage decision variable (y js ).

Case Study
Yongxin, a county of Jiangxi Province, is located in the southeast part of China (as shown in Figure 1).It ranges in longitude from 113.83 • E to 114.31 • E, and in latitude from 26.78 • N to 27.23 • N. It occupies a total area of around 2194.57km 2 , and the majority of the county (i.e., approximately 1800 km 2 ) lies within the middle reaches of the Heshui River Basin [31].The basin features subtropical monsoon humid climatic conditions (abundant rainfall and sunlight, and long frost-free periods) with an average annual precipitation of 1530.7 mm and an average annual temperature of 18.2 • C. The Heshui River has a total length of 225 km, with 77 km flowing within the borders of Yongxin County from west to east.The county mainly relies on the Heshui River to support its agriculture, industry, livestock husbandry, forestry, and fishery [32].
In the county, farmland, orchard, and woodland occupy15.2%,0.6%, and 70.1% of the total land, receptively.Grassland accounts for 6.1% of the total area, most of which is in good condition and suitable for feeding a large number of livestock.Land for urban construction, traffic, and water conservancy facilities cover 3.4%, 0.4%, and 0.6% of the total land, respectively.Agriculture is traditionally the primary sector of the County.Paddy soil is the main cultivation soil type, and it is not only suitable for rice planting but can also be used for rapeseeds and other cash crops.Fruit orchards produce a variety of fruits such as pear, peach, orange and plum.Agriculture includes paddy, dry and vegetable/fruit farms.The main types of livestock raised in the County are hogs, cattle and poultry.Pork and beef account for around 75.5 and 13.2% of the total meat production, respectively.The second industrial sector in the county is mainly comprised of mining, manufacturing, construction, transportation and other industries.The area possesses more than 20 types of mineral deposits, and the mineral production was 315.0 thousand tonnes [32,33].
Symmetry 2017, 9, x FOR PEER REVIEW 6 of 18 manufacturing, construction, transportation and other industries.The area possesses more than 20 types of mineral deposits, and the mineral production was 315.0 thousand tonnes [32,33].In the past twenty years, driven by poverty, authorities have given top priority to the booming economy in the local strategic plans [33].According to the report of local government in 2016, the total population was approximately 53,139, and the net income per capita was 8520 Chinese yuan (approximately $1290).One major environmental issue in the study area is water pollution, which is mainly caused by excessive pollutant loadings from agriculture, forestry, fishery, and livestock husbandry, and industry.Water pollution problems pose great obstacles to sustainable development in the county; thus, it is essential to optimize economic structure from a systematic point of view.Based on field investigations and related literature, chemical oxygen demand (COD), total nitrogen (TN), total phosphorus (TP), and soil loss have been selected as the pollutant types to control water pollution.Five zones (zones 1 to 5) are chosen to control water pollution.After the first-stage decisions (i.e., planning targets of economic activities) are made, associated pollution, the amount of which is in proportion to the economic production, is discharged into the water body.Based on the calculation of discharged wastewater and the measurement of incoming water quality, the mitigation schemes can be determined as the second-stage decisions in order to meet the environmental standards.If the production targets of them are made too high, the additionally generated pollutants will have to be mitigated in a more expensive way or discharged/drained into the stream (leading to penalties from the government).Conversely, if the production targets are made too low, the system will encounter opportunity losses of economic income, leading to a reduced system benefit.The In the past twenty years, driven by poverty, authorities have given top priority to the booming economy in the local strategic plans [33].According to the report of local government in 2016, the total population was approximately 53,139, and the net income per capita was 8520 Chinese yuan (approximately $1290).One major environmental issue in the study area is water pollution, which is mainly caused by excessive pollutant loadings from agriculture, forestry, fishery, and livestock husbandry, and industry.Water pollution problems pose great obstacles to sustainable development in the county; thus, it is essential to optimize economic structure from a systematic point of view.Based on field investigations and related literature, chemical oxygen demand (COD), total nitrogen (TN), total phosphorus (TP), and soil loss have been selected as the pollutant types to control water pollution.Five zones (zones 1 to 5) are chosen to control water pollution.After the first-stage decisions (i.e., planning targets of economic activities) are made, associated pollution, the amount of which is in proportion to the economic production, is discharged into the water body.Based on the calculation of discharged wastewater and the measurement of incoming water quality, the mitigation schemes can be determined as the second-stage decisions in order to meet the environmental standards.If the production targets of them are made too high, the additionally generated pollutants will have to be mitigated in a more expensive way or discharged/drained into the stream (leading to penalties from the government).Conversely, if the production targets are made too low, the system will encounter opportunity losses of economic income, leading to a reduced system benefit.The pollutant discharge allowances are highly uncertain, which are presented as probability distributions.Seven pollutant-discharge allowance levels are generated in association with different probabilities, which are identified as very-low, low, low-medium, medium, medium-high, high, and very-high, respectively.More uncertainties associated with cost/benefit coefficients may come from measurement errors and/or subjective judgment, and they can be expressed as T2FS.Therefore, the developed RTFP method can be adopted to plan water pollution control in the study area.The formulation of the RTFP model is presented as follows: (1) COD discharge constraints: (2) Phosphorus discharge constraints: (3) Nitrogen discharge constraints: (4) Soil loss constraints: (5) Water supply balance constraint: (6) Product demand constraints: Technical constraints: The nomenclature of variables and parameters is provided at the end of this ariticle.The imprecise input parameters are investigated according to field surveys, statistical data [31], government reports [34], and related literature [32,33].Table 1 presents unit benefit of agriculture and livestock husbandry, which are expressed as type-2 fuzzy sets.Table 2 provides pollutant allowances and the associated probabilities of occurrence.

Results Analysis and Discussion
Figure 2 represents the benefits under various pollutant discharge allowances and the system benefit, indicating that a higher allowance would correspond to a higher benefit.For instance, when the allowable pollutant discharge changes from very low to very high, the total benefit (i.e., sum of benefit of each activity) would be raise from $172.1 × 10 6 to $195.7 × 10 6 .Decisions at a higher allowable pollutant discharge would lead to a higher system benefit, but the reliability in fulfilling the environmental requirements would decrease; on the contrary, decisions at a lower allowable pollutant discharge would lead to a decreasing of risk for violating the pollutant discharge constraints, but with a lower system benefit.It demonstrates a trade-off between environmental requirement violation risk and benefit due to the uncertainties existing in various system components.In practice, when the plan aims to a higher system benefit, the environmental requirements may not be adequately satisfied; contrarily, planning with a lower system benefit may guarantee that the requirements be met.Additionally, the benefit of agriculture activity would take the largest proportion in total benefit and would increase slightly with the raising of pollutant discharge allowances.Moreover, the benefit of fishery would be stable at a low level, approximately occupying 2.1% of the total benefit.

Results Analysis and Discussion
Figure 2 represents the benefits under various pollutant discharge allowances and the system benefit, indicating that a higher allowance would correspond to a higher benefit.For instance, when the allowable pollutant discharge changes from very low to very high, the total benefit (i.e., sum of benefit of each activity) would be raise from $172.1 × 10 6 to $195.7 × 10 6 .Decisions at a higher allowable pollutant discharge would lead to a higher system benefit, but the reliability in fulfilling the environmental requirements would decrease; on the contrary, decisions at a lower allowable pollutant discharge would lead to a decreasing of risk for violating the pollutant discharge constraints, but with a lower system benefit.It demonstrates a trade-off between environmental requirement violation risk and benefit due to the uncertainties existing in various system components.In practice, when the plan aims to a higher system benefit, the environmental requirements may not be adequately satisfied; contrarily, planning with a lower system benefit may guarantee that the requirements be met.Additionally, the benefit of agriculture activity would take the largest proportion in total benefit and would increase slightly with the raising of pollutant discharge allowances.Moreover, the benefit of fishery would be stable at a low level, approximately occupying 2.1% of the total benefit.The water pollution control problem can simplified as a recourse-based fuzzy programming (RFP) problem by transforming the membership grades of T2FS into deterministic values.System benefit obtained from RFP ($180.8× 10 6 ) would be lower than that from RTFP ($191.1 × 10 6 ).This is due to the fact that cost/benefit coefficients are handled by RFP, resulting in higher loss of uncertain information than that handled by RTFP. Figure 3 presents the optimal target and standard production scale of each activity under RFP and RTFP.Target of each activity discharge pollutant exceeds standards can be calculated through multiplying excess scale of economic activity (i.e., targetstandard) by pollutant discharge rate.Results indicate that the excess planning scale of agricultural activity would be high.This may be attributed to their high crop yields and great selling prices.It is also depicted that the excess feeding size of livestock husbandry activity would also maintain high levels due to its high annual incomes.The excess outputs of industrial activity would also be significant since industry is promoted by the authority to push up the local income.Excess fishery The water pollution control problem can simplified as a recourse-based fuzzy programming (RFP) problem by transforming the membership grades of T2FS into deterministic values.System benefit obtained from RFP ($180.8× 10 6 ) would be lower than that from RTFP ($191.1 × 10 6 ).This is due to the fact that cost/benefit coefficients are handled by RFP, resulting in higher loss of uncertain information than that handled by RTFP. Figure 3 presents the optimal target and standard production scale of each activity under RFP and RTFP.Target of each activity discharge pollutant exceeds standards can be calculated through multiplying excess scale of economic activity (i.e., target-standard) by pollutant discharge rate.Results indicate that the excess planning scale of agricultural activity would be high.This may be attributed to their high crop yields and great selling prices.It is also depicted that the excess feeding size of livestock husbandry activity would also maintain high levels due to its high annual incomes.The excess outputs of industrial activity would also be significant since industry is promoted by the authority to push up the local income.Excess fishery and forestry activities would be low due to their limited planning lands.Furthermore, excess economic activities corresponding to RFP would be higher than that corresponding to RTFP.For instance, target of agricultural activity would be 180 ha in zone 1, while the standard production scale under RTFP and RFP would, respectively, be 156.4ha and 148.5 ha corresponding to very-low level.Thus, excess agricultural activity would be 23.6 ha and 31.5 ha, respectively.It is revealed that the varied uncertain information would affect the water pollution control plans.Any simplifications may result in unreliable or misleading plans.and forestry activities would be low due to their limited planning lands.Furthermore, excess economic activities corresponding to RFP would be higher than that corresponding to RTFP.For instance, target of agricultural activity would be 180 ha in zone 1, while the standard production scale under RTFP and RFP would, respectively, be 156.4ha and 148.5 ha corresponding to very-low level.Thus, excess agricultural activity would be 23.6 ha and 31.5 ha, respectively.It is revealed that the varied uncertain information would affect the water pollution control plans.Any simplifications may result in unreliable or misleading plans.Figure 4 displays the excess pollutant discharges under different pollutant discharge allowance levels corresponding to RFP and RTFP.It is indicated that amounts of excess pollutant discharges would reduce with the increased levels.For instance, under RTFP, the amounts of excess soil loss, COD, TP, and TN discharges would be 52.3 × 10 3 , 152.1, 22.4, and 99.2 t under low pollutant discharge allowance level; in comparison, under high pollutant discharge allowance level, they would respectively decrease to 1.3 × 10 3 , 8.2, 7.2, 5.4 t.In general, results discover that a more restrictive pollution control would result in a higher excess pollutant discharge while a looser pollution control would bring on a lower excess pollutant discharge.Figure 4 displays the excess pollutant discharges under different pollutant discharge allowance levels corresponding to RFP and RTFP.It is indicated that amounts of excess pollutant discharges would reduce with the increased levels.For instance, under RTFP, the amounts of excess soil loss, COD, TP, and TN discharges would be 52.3 × 10 3 , 152.1, 22.4, and 99.2 t under low pollutant discharge allowance level; in comparison, under high pollutant discharge allowance level, they would respectively decrease to 1.3 × 10 3 , 8.2, 7.2, 5.4 t.In general, results discover that a more restrictive pollution control would result in a higher excess pollutant discharge while a looser pollution control would bring on a lower excess pollutant discharge.Figure 4 displays the excess pollutant discharges under different pollutant discharge allowance levels corresponding to RFP and RTFP.It is indicated that amounts of excess pollutant discharges would reduce with the increased levels.For instance, under RTFP, the amounts of excess soil loss, COD, TP, and TN discharges would be 52.3 × 10 3 , 152.1, 22.4, and 99.2 t under low pollutant discharge allowance level; in comparison, under high pollutant discharge allowance level, they would respectively decrease to 1.3 × 10 3 , 8.2, 7.2, 5.4 t.In general, results discover that a more restrictive pollution control would result in a higher excess pollutant discharge while a looser pollution control would bring on a lower excess pollutant discharge.Figure 5 displays proportion of excess pollutant discharges under medium-high level.It is indicated that the excess soil loss, TN, and TP discharges of agriculture would account for about 97%, 88%, and 87% of the of the total discharge amounts, respectively.The high excess pollutants discharges are mainly related to its high targeted planning scale and discharge rates.The excess COD discharged from fishery would be high (accounting for 58% of the total excess discharge) due to its high COD discharge rate.In summary, the trade-off between agricultural income and pollution control (i.e., soil loss, TP, and TN discharges) would be of great concern for the local authority; fishery would be ceased due to its low benefit and high pollutant discharges.
Figure 5 displays proportion of excess pollutant discharges under medium-high level.It is indicated that the excess soil loss, TN, and TP discharges of agriculture would account for about 97%, 88%, and 87% of the of the total discharge amounts, respectively.The high excess pollutants discharges are mainly related to its high targeted planning scale and discharge rates.The excess COD discharged from fishery would be high (accounting for 58% of the total excess discharge) due to its high COD discharge rate.In summary, the trade-off between agricultural income and pollution control (i.e., soil loss, TP, and TN discharges) would be of great concern for the local authority; fishery would be ceased due to its low benefit and high pollutant discharges.Generally, compared with RTFP, RFP that tackles uncertainty in a single fuzzy set is not satisfactory to sufficiently reflect enough uncertain information in the decision-making process.Such a simplification may lead to unreliable or even misleading solutions.For instance, due to loss of the additional design degree of freedom, which is useful in water pollution control systems where input data are highly uncertain, economic coefficients (unit cost and unit benefit) treated by RFP are higher than that treated by RTFP.The higher economic coefficients lead to a lower system benefit, which may present unreliable decision support for local authorities.RTFP can tackle more complex uncertainty in terms of T2FS due to the extra degree of freedom.Moreover, solutions of economic activities scale and excess pollutant discharge also display that RTFP is more capable of balancing economic development and pollutant discharge control in the study area.Thus, RTFP is more enhanced in uncertainty reflection, as well as water pollution control.

Conclusions
In this study, a recourse-based type-2 fuzzy programming (RTFP) approach is developed for water pollution control planning.RTFP has incorporated the techniques of type-2 fuzzy programming and two-stage stochastic programming with recourse (TSP) within a general framework.RTFP can handle uncertainties presented as type-2 fuzzy sets (i.e., a fuzzy set in which the membership function is also fuzzy) and probability distributions.RTFP could provide benefit assessment for random pollutant discharge allowance levels through the two-stage framework.The  Generally, compared with RTFP, RFP that tackles uncertainty in a single fuzzy set is not satisfactory to sufficiently reflect enough uncertain information in the decision-making process.Such a simplification may lead to unreliable or even misleading solutions.For instance, due to loss of the additional design degree of freedom, which is useful in water pollution control systems where input data are highly uncertain, economic coefficients (unit cost and unit benefit) treated by RFP are higher than that treated by RTFP.The higher economic coefficients lead to a lower system benefit, which may present unreliable decision support for local authorities.RTFP can tackle more complex uncertainty in terms of T2FS due to the extra degree of freedom.Moreover, solutions of economic activities scale and excess pollutant discharge also display that RTFP is more capable of balancing economic development and pollutant discharge control in the study area.Thus, RTFP is more enhanced in uncertainty reflection, as well as water pollution control.

Conclusions
In this study, a recourse-based type-2 fuzzy programming (RTFP) approach is developed for water pollution control planning.RTFP has incorporated the techniques of type-2 fuzzy programming and two-stage stochastic programming with recourse (TSP) within a general framework.RTFP can handle uncertainties presented as type-2 fuzzy sets (i.e., a fuzzy set in which the membership function is also fuzzy) and probability distributions.RTFP could provide benefit assessment for random pollutant discharge allowance levels through the two-stage framework.The results are helpful for the authorities in exploring the trade-off between economic objective and pollutant discharge decision-making based on river water pollution control.
The developed RTFP method has been demonstrated through its application to water pollution control of Heshui River Basin.Different levels that pollutant discharge allowances are assumed to be random and benefit/cost coefficients are specified as type-2 fuzzy sets have been investigated.Solutions for production scale of agriculture, livestock husbandry, industry, forestry, and fishery are generated.Several findings can be concluded: (i) results reveal that uncertainties in pollutant discharge allowance level have significant effects on the city's future economic structure; (ii) constrained with low pollutant discharge allowance, the excess pollutant discharge would be high; (iii) the increased restriction of discharge would stimulate the basin taming the pace of economic growth; (iv) results also disclose that the benefit of agriculture activity would take the largest proportion in total benefit and would increase slightly with the raising of pollutant discharge allowances; and (v) system benefit is powerfully impressed by the pollutant discharge allowance level.
It is the first attempt to apply the RTFP method to water pollution control planning, and results indicate that (1) RTFP cannot only handle uncertainty expressed as T2FS, but also effectively cope with uncertainty described as probability distribution; (2) TFCP can help authorities to make trade-offs among environmental violation risk and system benefit.Nevertheless, there are also potential extensions of the proposed method.For example, RTFP has difficulty in addressing uncertainties that cannot be expressed as type-2 fuzzy sets or probability distributions; interval is a proper type to present such a kind of uncertainty.Furthermore, construction of T2FS membership function is usually based on experts' subjective evaluation, which may also lead to high uncertainty.At present, there is no appropriate evaluation method except experts' decisions.It is worthwhile to search for more effective methods to improve the method of constructing membership functions for T2FS, which can be helpful to reduce uncertainty.
js denotes the second-stage cost function ( e j is cost coefficient of y js ).Inequality (4e) presents the relationship among x j , y js , and ω s .

Figure 2 .
Figure 2. Benefit of each activity under each pollutant discharge allowance.

Figure 2 .
Figure 2. Benefit of each activity under each pollutant discharge allowance.
calculate the relative distance matrix D = d ij t×t , where d ij = h i − h j ; (2) calculate the average of relative distancesd i = ∑ t i=1 d ij /(t − 1); (3) introduce a pair-wise comparison p ij (p ij = d j /d i ), and the pair-wise matrix P = p ij t×t ; (4) obtain the true-importance degree w i of h i (w j = 1/∑ t j=1 p ij ); (5) assess the mode of m (m = ∑ t i=1 w i h i ) of the fuzzy number; (6) choose s = ∑ t i=1 w i |h i − m| to approximate the unknown mean deviation σ; (7) acquire the ratio of left spread to the right spread η

Table 1 .
Unit benefits of agriculture and livestock husbandry.

Table 2 .
Pollutant allowances under different levels.
PEF ikreduction of net benefit from fishery activity i for excess discharge of pollutant k (RMB¥/kg) DPF ijk discharge rate of pollutant k from fishery activity i in zone j (kg/km 2 )XF ijhk decision variables representing amount by which target of fishery activity i discharge pollutant k exceeds standards in zone j when level is h (km 2 ) BL ± ij unit benefit from livestock husbandry activity i in zone j (RMB¥/head) TL ij target for livestock husbandry activity i in zone j (head) PEL i reduction of net benefit from livestock husbandry activity i for excess discharge of pollutant (i.e., COD) (RMB¥/kg) DPL ij discharge rate of pollutant (i.e., COD) from livestock husbandry activity i in zone j (kg/head) XL ijh decision variables representing amount by which target of livestock husbandry activity i discharge pollutant (i.e., COD) exceeds standards in zone j when level is h (head) TI ij output target for industrial activity i in j (RMB¥) PEI i reduction of net benefit from industrial activity i for excess discharge of pollutant (i.e., COD) (RMB¥/kg) DPI ij discharge rate of pollutant (i.e., COD) from industrial activity i in zone j (kg/RMB¥) XI ijh decision variables representing amount by which target of industrial activity i discharge pollutant (i.e., COD) exceeds standards in zone j when level is h (RMB¥) BW ij unit benefit from forestry activity i in zone j (RMB¥/head) TW ij land area target for forestry activity i in zone j (unit) PEW i reduction of net benefit from forestry activity i for excess discharge of pollutant (i.e., soil loss) (RMB¥/tonne) DPW ij discharge rate of pollutant (i.e., soil loss) from forestry activity i in zone j (tonne/km 2 ) XW ijh decision variables representing amount by which target of forestry activity i discharge pollutant (i.e., soil loss) exceeds standards in zone j when level is h (unit) COF ij COD discharge from fishery farming activity i in zone j (kg/km 2 ) PCF jh maximum allowable COD discharge for fishery farming activities in zone j with probability p h of occurrence under level h (kg) PCL jh maximum allowable COD discharge for livestock husbandry activities in zone j with probability p h of occurrence under level h (kg) PCI jh maximum allowable COD discharge for industrial activity i in zone j with probability p h of occurrence under level h (kg) MCL h maximum allowable COD discharge from economic activities with probability p h of occurrence under level h (kg) SL ij soil loss from agricultural activity i in zone j (tonne/km 2 ) AP ij phosphorous content of soil corresponding to agricultural activity i in zone j (kg/tonne) RA ij runoff from agricultural activity i in zone j (kg/km 2 ) RP ij dissolved phosphorous content of runoff corresponding to agricultural activity i in zone j (%) PAP jh maximum allowable phosphorous loss from agricultural activities in zone j with probability p h of occurrence under level h (kg) FP ij dissolved phosphorous loss from fishery farming activity i in zone j (kg/km 2 ) PFP jh maximum allowable phosphorous loss from fishery farming activities in zone j with probability p h of occurrence under level h (kg) MPL h maximum allowable phosphorous loss from economic activities with probability p h of occurrence under level h (kg) AN ij nitrogen content of soil corresponding to agricultural activity i in zone j (kg/tonne) RN ij dissolved nitrogen content of runoff corresponding to agricultural activity i in zone j (%) PAN jh maximum allowable nitrogen loss from agricultural activities in zone j with probability p h of occurrence under level h (kg) FN ij dissolved nitrogen loss from fishery activity i in zone j (kg/km 2 ) PFN jh maximum allowable nitrogen loss from fishery farming activities in zone j with probability p h of occurrence under level h (kg)