Robust Linear Programming and Its Application to Water and Environmental Decision-Making under Uncertainty

In this study, we introduce a robust linear programming approach for water and environmental decision-making under uncertainty. This approach is of significant practical utility to decision makers for obtaining reliable and robust management decisions that are “immune” to the uncertainty attributable to data perturbations. The immunization guarantees that the chosen robust management plan will be implementable with no violation of the mandatory constraints of the problem being studied—i.e., natural resource supply constraint, environmental carrying capacity constraint, environmental pollution control constraint, etc.—and that the actual value of the objective will be no worse than the given estimation if the perturbations of data fall within the specified uncertainty set. A simplified example in regional water quality management is provided to help water and environmental practitioners to better understand how to implement robust linear programming from the perspective of application, as well as to illustrate the significance and necessity of implementing robust optimization techniques in real-world practices. Robust optimization is a growing research field that requires more interdisciplinary research efforts and engagements from water and environmental practitioners. Both may benefit from the advances of management science.


Introduction
Scarcity of freshwater has become a continuing worldwide problem as population increases, urbanization advances and climate changes.In many regions of the world, freshwater resources are insufficient to meet domestic, economic development and environmental needs.Meanwhile, the lack of adequate clean water has become a constraint on the advancement of economic development and the maintenance of a safe environment and healthy ecosystems [1][2][3][4].To achieve long-term economic prosperity and environmental sustainability, policy-makers and water practitioners must engage with and actively collaborate with one another to explore a safe, sustainable way without compromising water availabilities that are essential to ecosystem functions and services [5].Innovative water and land management approaches are needed for better planning and management of valuable freshwater resources and prevention of water-related risks and vulnerabilities.Multiple considerations in the social, economic and environmental aspects of water also need to be taken into account in decision-making processes in order to develop solutions that promote sustainable development over the long term [6][7][8][9][10].
It is already known that scientific management and contingency planning are the keys for addressing water-related conflicts and problems.A multitude of studies have been conducted to help decision makers to translate the sustainable development concept into operational management plans [11][12][13][14].Mathematical programming, as one of the most useful tools that evaluate the competition scheme between sectors such as agriculture, industry and energy for limited natural resources, receives a continued interest from water, energy and environmental practitioners [15][16][17][18][19]. Mathematical programming has its particular advantages in providing a simplified and effective way for approaching the management problem of interest and in obtaining optimal solutions that reflect the tradeoffs between conflicting interests.A number of academic and industrial practices were conducted on the basis of an assumption that the problem data are known with certainty [20][21][22][23][24].However, this is seldom the situation in practice.For example, the data may be derived from inferior estimation, affected by observational errors or affected by random events.In fact, uncertainty is a certain component of many natural systems.Failure to consider uncertainty in scientific and engineering practices may significantly reduce the reliability and robustness of model results, especially for the case that correction of decisions is not a readily available option for decision makers when the actual data reveals itself.As a result, there is a growing recognition that it is necessary to translate uncertain information into model parameters in order to obtain more safe and reliable solutions under uncertainty.
Previous research efforts have indeed improved our understanding of the applicable methods for dealing with various uncertainties in a water system [25][26][27][28][29].Multiple types of uncertain information, including interval number, fuzzy set and probability distribution, can be translated into model parameters effectively through different mathematical techniques [30][31][32][33][34][35].However, most of these efforts have limited practical values due to the limitation on providing reliable and robust managerial support under uncertainty.This is because these efforts were only focused on uncertainty characterization and propagation rather than the reliability analysis of decisions under uncertainty.For example, the inexact programming technique that has been extensively studied during the past two decades can provide a set of optimal interval solutions for supporting decision-making.However, this set is a collection of optimal solutions under a range of uncertain scenarios and the optimality of each solution from the set can only be guaranteed under specific scenarios rather than all scenarios.Accordingly, the optimal solutions from such a set are as vulnerable as the nominal optimal solutions that are derived without consideration of uncertainty if the risk assessment of constraint violations has not been carried out.Hence the research focus shall be aimed at how to obtain optimal, uncertainty-resistant solutions for supporting decision-making under uncertainty; otherwise, the solutions offered to decision makers are less meaningful in practice or much riskier than expected.
In this study, we introduce a robust linear programming approach that is capable of providing robust optimal solutions that are resistant to the uncertainty attributable to data perturbations.Infeasibility of the nominal optimal solutions can be effectively prevented and decision makers shall be confident of the reliability of their decisions when the actual data is within the given uncertainty set.Robust optimization is a young research field that has been mainly developed in the most recent 10 years.It has been applied in finance, operations research, electrical engineering, chemical engineering and computer science.Only a few exploratory studies have been carried out on the basis of the concept of robust optimization by water and environmental practitioners.It is thus necessary to promote this newly developed approach to water and environmental practitioners for obtaining reliable decision-making support under uncertainty.A semi-hypothetical example in water quality management is provided to illustrate the applicability of the robust linear programming approach and the vulnerability of the nominal optimal solutions in practice.The data of this example is mainly based on government reports and the numerical example has been simplified for the purpose of illustration.This example may help water and environmental practitioners to better understand the significance and necessity of implementing robust optimization techniques in real-world practices.

Methodology
Robust optimization is an applied research field that has received a revival of interest by a number of management scientists and operational researchers.The methodology of robust linear programming aims at finding solutions that are resistant to data uncertainty, especially for the case where the actual realization of uncertain parameters is within a given uncertainty set.Here a brief introduction of this methodology is given through a simple linear optimization problem with uncertain parameters.The general model of this type of linear optimization problems is as follows: where x is the decision variable and c ∈ R n , A ∈ R m×n and b ∈ R m are the uncertain parameters.
In robust linear programming, the objective is to find the robust optimal value of t within a set of robust feasible solutions that enable the inequality c T x ≤ t to be active for all realizations of the objective function f (x) from the uncertainty set Z. Hence, Model 1 can be converted into its robust counterpart as follows [36][37][38][39]: Min subject to c T x ≤ t (5) x ≥ 0 (7) Here Constraints 5 and 6 with uncertain parameters c and b can be reformed as the following constraints: where x n+1 = −1.Accordingly, without loss of generality, the above constraints can be expressed as a universal constraint as follows: where b ∈ R m becomes a certain parameter and ∀ε ∈ Z. Since uncertainty is constraint-wise in robust optimization, decision makers may need to focus on the ith constraint of Constraint 10 in a robust linear programming problem: where a ij (ε) = a ij + a ij ε ij ; a ij is the nominal value of a ij (ε); and a ij is the positive length of perturbations.Constraint 11 then can be rewritten as follows: Robust linear programming requires the determination of a robust solution that is feasible for any realization of ε in the given uncertainty set Z. Accordingly, the immunized constraint can be expressed as follows: where Z is the specified uncertainty set and the basic set is the box uncertainty set Z = {ε| ε ∞ ≤ 1 }.

Case Study
A semi-hypothetical regional water quality management problem is provided to demonstrate how to implement robust linear programming in practice.Consider a small watershed system that has a nature reserve in the area downstream of a river, as illustrated in Figure 1.Two dairy and beverage factories are the major polluters in the watershed system and are located at the upstream end of the river.A restrictive water quality regulation has been issued to protect the ecosystem services and functions of the nature reserve.It is estimated that the daily maximum allowable biochemical oxygen demand (BOD) mass load discharged from the dairy and beverage factories is 170 kg, and that this BOD mass limit should be restrictively complied with and enforced so that the downstream water quality in the nature reserve can be maintained pursuant to the local water quality regulation.The dairy factory produces wastewater of 1200 m 3 per day with the BOD concentration of 850 mg/L and the beverage factory produces wastewater of 4000 m 3 per day with the BOD concentration of 400 mg/L.The dairy and beverage factories have their own wastewater treatment systems and the BOD removal rates of the systems are 90% and 95%, respectively.The treated wastewater from both the dairy and beverage factories satisfies the local effluent water quality standards and is permitted to be discharged into the river directly.The untreated wastewater is prohibited to be discharged directly into the river and may be discharged into the municipal sewage system with additional service fees.Different service fees apply for the untreated wastewater discharged into the municipal sewage system depending on the BOD concentration of the untreated wastewater.The watershed manager needs to decide how much wastewater is treated at the dairy and beverage factories by their own capacities and how much wastewater is discharged into the municipal sewage system without treatment.The objective of the watershed manager is to obtain a cost-effective management plan that satisfies the local water quality regulation.The watershed manager cannot tolerate violations of the regulatory water quality requirement.The summary of the problem parameters for the regional water quality management problem is provided in Table 1.
Table 1.Parameter values for the regional water quality management problem illustrated in Figure 1.Table 1.Parameter values for the regional water quality management problem illustrated in Figure 1.This water quality management problem can be formulated as the following optimization model:

Influent
Subject to: (Water quality requirement) ) , (Wastewater treatment capacity) (Wastewater discharge capacity) where: This water quality management problem can be formulated as the following optimization model: Subject to: (Water quality requirement) (Wastewater treatment capacity) (Wastewater discharge capacity) where: C i = the cost of BOD removal at factory i ($/ton) B i = the BOD concentration of wastewater at factory i (10 −6 ton/m 3 ) S i = the service fee of wastewater discharged into the municipal sewage system from factory i ($/m 3 ) E i = the BOD removal efficiency of the wastewater treatment system at factory i (%) L = the BOD mass limit (10 −3 ton) W i = the daily volume of wastewater produced at factory i (m 3 /day) x i = the volume of wastewater treated by the wastewater treatment system of factory i (m 3 ) y i = the volume of wastewater discharged into the municipal sewage system from factory i (m 3 ) As the perturbation of BOD removal efficiency at each factory is taken into account, the above optimization model can be converted to the following robust counterpart model: Subject to: (Water quality requirement) (Wastewater treatment capacity) (Wastewater discharge capacity) where E i is the nominal value of E i ; E i is the positive length of the perturbation of is the random variable subject to uncertainty; and Z is the specified uncertainty set.
Here we consider a box uncertainty set and a case that the BOD removal rates drift in the 5% and 2% margins of the nominal values at the dairy and beverage factories, respectively.

Result Analysis
In this example we studied the nominal model at first so that the nominal optimal solution of Model 9 can be used to provide a detailed contrast to the robust optimal solution of the robust counterpart model, i.e., Model 10.The solution of Model 9 can be obtained easily and it suggests such a management plan: treat 1200 m 3 of wastewater (no discharge into the municipal sewage system) at the dairy factory and treat 3400 m 3 of wastewater (600 m 3 of wastewater discharged into the municipal sewage system) at the beverage factory for a total cost of $1238 per day.It is known that not all problem data could be viewed as perfectly reliable in practice.For example, one can hardly believe that the BOD removal rates of the wastewater treatment systems are exactly the nominal data, i.e., 90% at the dairy factory and 95% at the beverage factory.In reality, these rates vary around the nominal values.If the perturbations of the BOD removal rates are taken into account, the nominal optimal solution may make the model infeasible.In fact, this nominal optimal solution is sensitive to the model parameters provided in this case study.In particular, since the water quality requirement constraint (i.e., Constraint 15) is active with this optimal solution (i.e., there is no additional BOD amount allowable to be discharged into the river), a slight perturbation of the parameters that describe this constraint could make this nominal optimal solution infeasible.In other words, the regulatory water quality requirement could be violated when a small change occurs in the BOD removal rates.
Robust linear programming provides a feasible way to prevent such violations of the decision criteria and suggests a more rational and robust optimal solution for the management problem under uncertainty.As the perturbations of the BOD removal rates by 5% and 2% at the dairy and beverage factories are taken into account, the robust optimal solution of Model 10 suggests the following management plan: treat 483.57m 3 of wastewater (716.43 m 3 of wastewater discharged into the municipal sewage system) at the dairy factory and treat 4000 m 3 of wastewater (no discharge into the municipal sewage system) at the beverage factory for a total cost of $1366 per day.Since the nominal solution of Model 9 leads to the violation of the water quality requirement constraint, it is necessary to assume a monetary system penalty with respect to the violation so that the negative effect of the violation caused by the non-robust nominal solution could be evaluated.Hence, the system penalty is assumed to be $5000 per ton for the unauthorized discharge of BOD mass load into the river.Now the negative effect of the violation can be estimated and shows that the nominal optimal solution would result in a much higher total cost, increasing from $1238 to $1597 in the worst scenario.The actual total cost will increase by, at most, 29% compared to the expected value.Compared with the nominal optimal solution, the robust optimal solution suggests a higher total cost to prevent the violation of the regulatory water quality requirement.The estimated robust total cost increases the nominal total cost by 10% (i.e., increasing from $1238 to $1366) but it guarantees that no extra cost will be added to this estimation if the perturbations of parameters are within the given uncertainty set.This provides two decision-making alternatives: a certain 10% increase in the total cost with no penalty or a lower expected total cost with a potential increase by as much as 29%.It is apparent that the robust solution is favorable in practice since most decision makers prefer uncertainty-aversion decisions.Moreover, from a viewpoint of the tradeoff between the cost and the risk, it is still favorable in adding 10% initial cost for insurance to avoid a potential 29% cost overrun.
The comparison between the solutions of Models 9 and 10 is given in Figure 2. It shows two conflicting management plans for this regional water quality management problem.The nominal solution suggests that the dairy factory needs to treat all its wastewater by its own wastewater treatment system and no wastewater is required to be discharged through the municipal sewage system.Indeed, it is a somewhat easy decision to be made since it is consistent with our common sense for a minimization problem.The dairy factory has a lower BOD removal cost and a higher sewage service fee attributable to the higher BOD concentration within its wastewater.It seems to be reasonable by giving the dairy factory a top priority for the purpose of minimizing the total cost.However, this is true in the theoretical condition where no uncertainty is involved.As the perturbations of the BOD removal rates are taken into account, a different priority rule for the dairy and beverage factories may be given.In fact, the robust counterpart model gives a conflicting management plan.The robust solution suggests that the beverage factory needs to treat all its wastewater by its own wastewater treatment system rather than giving the dairy factory such a priority.The dairy factory needs to treat only a portion of its wastewater (i.e., 483.57m 3 ) and the remaining wastewater (i.e., 716.43 m 3 ) needs to be discharged into the municipal sewage system directly with a certain amount of sewage service fee.Although the beverage factory has a higher BOD removal cost per unit, the robust solution indicates that giving the beverage factory a top priority over the dairy factory to treat all its wastewater is a robust and cost-effective decision under uncertainty.This robust solution can be considered "immune" to respective 5% and 2% fluctuations of the BOD removal efficiency at the dairy and beverage factories.The immunization guarantees that the robust water quality management plan will be implementable without violation of the regulatory water quality requirement and that the total cost will be no more than the expected estimation no matter what the BOD removal efficiency ends up being, as long as it falls within the defined uncertainty set.
We also studied four boundary scenarios of uncertain parameter combinations for nominal and robust solutions.As can be seen from Figure 3, the actual cost of the nominal solution could increase by $359 in the worst scenario (i.e., the scenario of "--" where -denotes the lower bound of uncertain parameters) and the BOD over-discharging load due to negligent environmental violation is 72 kg in total.Compared with the worst scenario, the actual cost of the nominal solution in one of the non-pessimistic scenarios (i.e., the scenario of "-+") could increase by about one-third of the portion of cost overrun in the worst scenario and the BOD over-discharging load is about 20 kg.It indicates that the nominal solution could also result in negligent environmental violation in a non-pessimistic scenario.In the other two scenarios, the nominal solution may lead to no environmental violations but this solution is not optimal for the two scenarios.The actual cost could be further reduced by about 2.5% if modification of the nominal decisions is allowable after the uncertain parameters reveal their actual values.Figure 4 shows the risk and optimality analysis of the robust solution under the four boundary scenarios.It shows that the robust solution guarantees no environmental violations in all four scenarios but is less optimal in the last three scenarios.The actual cost could be further reduced by at least 8% in these scenarios if modification of the decisions is allowable after the uncertain parameters reveal their actual values.
The comparison between the nominal and robust optimal solutions with respect to the feasibility of the water quality requirement constraint within the defined uncertainty set is provided in Figure 5.As can be seen from Figure 5a, the nominal solution makes the water quality requirement constraint feasible with a probability of 0.50 within the given box uncertainty set.As a result, the watershed manager has to undertake a 50% chance of violating the regulatory water quality requirement without informed knowledge of such a risk.In contrast to the nominal solution, the robust solution, shown in Figure 5b, makes the water quality requirement constraint feasible within the entire uncertainty set.Hence robust linear programming techniques may give the watershed manager more confidence to implement a reliable and robust water quality management plan under uncertainty.
four boundary scenarios.It shows that the robust solution guarantees no environmental violations in all four scenarios but is less optimal in the last three scenarios.The actual cost could be further reduced by at least 8% in these scenarios if modification of the decisions is allowable after the uncertain parameters reveal their actual values.
The comparison between the nominal and robust optimal solutions with respect to the feasibility of the water quality requirement constraint within the defined uncertainty set is provided in Figure 5.As can be seen from Figure 5a, the nominal solution makes the water quality requirement constraint feasible with a probability of 0.50 within the given box uncertainty set.As a result, the watershed manager has to undertake a 50% chance of violating the regulatory water quality requirement without informed knowledge of such a risk.In contrast to the nominal solution, the robust solution, shown in Figure 5b, makes the water quality requirement constraint feasible within the entire uncertainty set.Hence robust linear programming techniques may give the watershed manager more confidence to implement a reliable and robust water quality management plan under uncertainty.

Discussion
This regional water quality management example was developed for the purpose of illustrating the application of robust linear programming and the main limitation of previous methods in dealing with water and environmental management problems under uncertainty.A number of simplifications with respect to real-world water quality management problems had been made to help water and environmental practitioners to better understand the significance and necessity of implementing robust optimization techniques in real-world managerial practices.Although this representative example was simplified for the purpose of illustration, its generality and flexibility allow this example to be further improved and extended to real-world applications with some modifications based on decision makers' actual needs.For example, more water quality parameters such as chemical oxygen demand (COD), total suspended solids (TSS), NH4-N, NO3-N, and total phosphorus may need to be included.Detailed cost estimation may need to be performed.The process of environmental pollutant transport in the study area may need to be simulated prior to the development of a regional management model.These improvements will significantly enhance the suitability and practicality of this example in dealing with real-world problems.
Robust optimization is an optimization technique that focuses on the worst scenarios under

Discussion
This regional water quality management example was developed for the purpose of illustrating the application of robust linear programming and the main limitation of previous methods in dealing with water and environmental management problems under uncertainty.A number of simplifications with respect to real-world water quality management problems had been made to help water and environmental practitioners to better understand the significance and necessity of implementing robust optimization techniques in real-world managerial practices.Although this representative example was simplified for the purpose of illustration, its generality and flexibility allow this example to be further improved and extended to real-world applications with some modifications based on decision makers' actual needs.For example, more water quality parameters such as chemical oxygen demand (COD), total suspended solids (TSS), NH4-N, NO3-N, and total phosphorus may need to be included.Detailed cost estimation may need to be performed.The process of environmental pollutant transport in the study area may need to be simulated prior to the development of a regional management model.These improvements will significantly enhance the suitability and practicality of this example in dealing with real-world problems.
Robust optimization is an optimization technique that focuses on the worst scenarios under uncertainty.It guarantees the certainty of the objective value and the feasibility of robust decisions in all scenarios.As a result, the optimality of robust decisions in other scenarios has to be compromised.In this study, we studied the basic uncertainty set in robust optimization-the box uncertainty set.This uncertainty set ensures that no violation of constraints will occur if the perturbations of uncertain parameters fall within the box uncertainty set.However, the robust solution obtained from this set is particularly conservative so that its practicality may also be discounted unless decision makers have no tolerance for the violation of constraints.If the decision makers are willing to accept a small degree of violation of the predetermined constraints, either the ellipsoidal uncertainty set or the polyhedral uncertainty set may be employed for obtaining a relatively conservative solution with competitive robustness under uncertainty.

Conclusions
In this study, we introduced a robust linear programming approach for water and environmental decision-making under uncertainty.This approach is of significant practical utility to decision makers for obtaining reliable and robust management decisions that are "immune" to the uncertainty attributable to data perturbations.The immunization guarantees that the chosen robust management plan will be implementable with no violation of the mandatory constraints of the problem being studied-i.e., natural resource supply constraint, environmental carrying capacity constraint, environmental pollution control constraint, etc.-and that the actual value of the objective will be no worse than the given estimation if the perturbations of data fall within the specified uncertainty set.Prevention or reduction of the violation of the mandatory constraints is particularly important to a successful water or environmental management plan, especially when it is known that the nominal solution may have a 50% chance of violating the mandatory constraints under uncertainty.
Nowadays the common management strategy that has been widely applied in water and environmental management, i.e., hedging the risk of environmental violations by economic compensation, is not recommended in practice.The natural world can no longer function as it was and has become less tolerant of negligent environmental violations than us-the debtors.So-called "economic compensation" is collected by people rather than the natural world itself as the creditor.Furthermore, environmental and ecological restoration requires a long-term effort and is far more difficult and costly than our expectation.It thus requires decision makers to be more thoughtful about developing an implementable management plan with resistance to uncertainty.A conservative attitude from decision makers with respect to environmental violations may provide more flexibility in developing a sustainable future for our society.
A simplified numerical example in water quality management was provided to help water and environmental practitioners to better understand how to implement robust linear programming from the perspective of application.This representative example is extendable and can be further improved based on actual needs in real-world practices.Overall, robust optimization is a growing research field that requires more interdisciplinary research efforts and engagements from water and environmental practitioners.Its main practical value lies in potentially providing water and environmental decision makers with a more rational and reliable way to address management and planning problems under uncertainty.In fact, the risk of negligent environmental violations attributable to non-robust decisions shall not be underestimated under uncertain conditions.The risk level associated with the optimal solution needs to be realized by decision makers and be appropriately considered in the decision-making process.Expected optimality shall no longer be the dominant decision criterion and a sustainable management plan shall always be a result of the tradeoff between the optimality and robustness of the management decisions.It is expected that both water and environmental decision makers may benefit from the advances of management science.In future research, complicated real-world cases will be studied, and the ellipsoidal and polyhedral uncertainty sets will be employed to evaluate the practicality of robust optimization.

Figure 1 .
Figure 1.A schematic diagram of the regional water quality management problem.

iCFigure 1 .
Figure 1.A schematic diagram of the regional water quality management problem.

Figure 2 .
Figure 2. The comparison of nominal optimal and robust optimal solutions.Figure 2. The comparison of nominal optimal and robust optimal solutions.

Figure 2 . 13 Figure 3 .
Figure 2. The comparison of nominal optimal and robust optimal solutions.Figure 2. The comparison of nominal optimal and robust optimal solutions.Sustainability 2018, 10, x FOR PEER REVIEW 9 of 13

Figure 3 .
Figure 3. Risk and optimality analysis of the nominal solution under four boundary scenarios of uncertain parameters combinations (where − and + denote the lower and upper bounds).

Figure 3 .
Figure 3. Risk and optimality analysis of the nominal solution under four boundary scenarios of uncertain parameters combinations (where -and + denote the lower and upper bounds).

Figure 4 .
Figure 4. Risk and optimality analysis of the robust solution under four boundary scenarios of uncertain parameters combinations (where -and + denote the lower and upper bounds).

Figure 4 . 13 Figure 5 .
Figure 4. Risk and optimality analysis of the robust solution under four boundary scenarios of uncertain parameters combinations (where − and + denote the lower and upper bounds).Sustainability 2018, 10, x FOR PEER REVIEW 10 of 13

Figure 5 .
Figure 5.The comparison of the feasible region between nominal optimal solution and robust optimal solution with respect to the water quality requirement constraint within the specified uncertainty set.