Resources and environmental systems management under synchronic interval uncertainties

Resources and environmental systems management (RESM) is challenged by the synchronic effects of interval uncertainties in the related practices. The synchronic interval uncertainties are misrepresented as random variables, fuzzy sets, or interval numbers in conventional RESM programming techniques including stochastic programming. This may lead to ineffectiveness of resources allocation, high costs of recourse measures, increased risks of unreasonable decisions, and decreased optimality of system profits. To fill the gap of few corresponding studies, a synchronic interval linear programming (SILP) method is proposed in this study. The proposition of interval sets and interval functions and coupling them with linear programming models lead to development of an SILP model for RESM. This enables incorporation of interval uncertainties in resource constraints and synchronic interval uncertainties in the programming objective into the optimization process. An analysis of the distribution-independent geometric properties of the feasible regions of SILP models results in proposition of constraint violation likelihoods. The tradeoff between system optimality and constraint violation is analyzed. The overall optimality of SILP systems under synchronic intervalness is quantified through proposition of integrally optimal solutions. Integration of these efforts leads to a violation-constrained interval integral method for optimization of RESM systems under synchronic interval uncertainties. Comparisons with selected existing methods reveal the effectiveness of SILP at eliminating negativity of synchronic intervalness, enabling risk management of and achieving overall optimality of RESM systems, and enhancing the reliability of optimization techniques for RESM problems. The exploited framework for analyzing synchronic interval uncertainties in RESM systems is helpful for addressing synchronisms of other uncertainties such as randomness or fuzziness and avoiding the resultant decision mistakes and disasters due to neglecting them.     

Fuzzy chance-constrained linear fractional programming approach for optimal water allocation

A fuzzy chance-constrained linear fractional programming method was developed for agricultural water resources management under multiple uncertainties. This approach improved upon the previous programming methods, and could reflect the ratio objective function and multiple uncertainties expressed as probability distributions, fuzzy sets, and their combinations. The proposed approach is applied to an agricultural water resources management system where many crops are considered under different precipitation years. Through the scenarios analyses, the multiple alternatives are presented. The solutions show that it is applicable to practical problems to address the crop water allocation under the precipitation variation and sustainable development with ratio objective function of the benefit and the irrigation amount. It also provides bases for identifying desired agriculture water resources management plans with reasonable benefit and irrigation schedules under crops.     

Two-stage fuzzy chance-constrained programming: application to water resources management under dual uncertainties

In this study, a two-stage fuzzy chance-constrained programming (TFCCP) approach is developed for water resources management under dual uncertainties. The concept of distribution with fuzzy probability (DFP) is presented as an extended form for expressing uncertainties. It is expressed as dual uncertainties with both stochastic and fuzzy characteristics. As an improvement upon the conventional inexact linear programming for handling uncertainties in the objective function and constraints, TFCCP has advantages in uncertainty reflection and policy analysis, especially when the input parameters are provided as fuzzy sets, probability distributions and DFPs. TFCCP integrates the two-stage stochastic programming (TSP) and fuzzy chance-constrained programming within a general optimization framework. TFCCP incorporates the pre-regulated water resources management policies directly into its optimization process to analyze various policy scenarios; each scenario has different economic penalty when the promised amounts are not delivered. TFCCP is applied to a water resources management system with three users. Solutions from TFCCP provide desired water allocation patterns, which maximize both the system’s benefits and feasibility. The results indicate that reasonable solutions were generated for objective function values and decision variables, thus a number of decision alternatives can be generated under different levels of stream flows, α-cut levels and fuzzy dominance indices.     

Agricultural effluent control under uncertainty: An inexact double-sided fuzzy chance-constrained model

An inexact double-sided fuzzy chance-constrained programming (IDFCCP) method was developed in this study and applied to an agricultural effluent control management problem. IDFCCP was formulated through incorporating interval linear programming (ILP) into a double-sided fuzzy chance-constrained programming (DFCCP) framework, and could be used to deal with uncertainties expressed as not only possibility distributions associated with both left- and right-hand-side components of constraints but also discrete intervals in the objective function. The study results indicated that IDFCCP allowed violation of system constraints at specified confidence levels, where each confidence level consisted of two reliability scenarios. This could lead to model solutions with high system benefits under acceptable risk magnitudes. Furthermore, the introduction of ILP allowed uncertain information presented as discrete intervals to be communicated into the optimization process, such that a variety of decision alternatives can be generated by adjusting the decision-variable values within their intervals. The proposed model could help decision makers establish various production patterns with cost-effective water quality management schemes under complex uncertainties, and gain in-depth insights into the trade-offs between system economy and reliability.     

Bivariate frequency analysis of flood and extreme precipitation under changing environment: case study in catchments of the Loess Plateau,China

Floods have changed in a complex manner, triggered by the changing environment (i.e., intensified human activities and global warming). Hence, for better flood control and mitigation in the future, bivariate frequency analysis of flood and extreme precipitation events is of great necessity to be performed within the context of changing environment. Given this, in this paper, the Pettitt test and wavelet coherence transform analysis are used in combination to identify the period with transformed flood-generating mechanism. Subsequently, the primary and secondary return periods of annual maximum flood (AMF) discharge and extreme precipitation (Pr) during the identified period are derived based on the copula. Meanwhile, the conditional probability of occurring different flood discharge magnitudes under various extreme precipitation scenarios are estimated using the joint dependence structure between AMF and Pr. Moreover, Monte Carlo-based algorithm is performed to evaluate the uncertainties of the above copula-based analyses robustly. Two catchments located on the Loess plateau are selected as study regions, which are Weihe River Basin (WRB) and Jinghe River Basin (JRB). Results indicate that: (1) the 1994–2014 and 1981–2014 are identified as periods with transformed flood-generating mechanism in the WRB and JRB, respectively; (2) the primary and secondary return periods for AMF and Pr are examined. Furthermore, chance of occurring different AMF under varying Pr scenarios also be elucidated according to the joint distribution of AMF and Pr. Despite these, one thing to notice is that the associate uncertainties are considerable, thus greatly challenges measures of future flood mitigation. Results of this study offer technical reference for copula-based frequency analysis under changing environment at regional and global scales.     

Total consumption controlled water allocation management for multiple sources and users with inexact fuzzy chance-constrained programming: a case study of Tianjin,China

Environmental and ecological issues caused by water resources crisis have brought enormous challenges to the sustainable development of water-deficient area. Water resources allocation management balancing the relationship between the social-economic development and the ecological environment has become a hot topic in recent years. In this paper, an inexact fuzzy chance-constrained programming (IFCCP) approach is proposed for regional water resource allocation optimization with the aim of promoting the harmonious development of the social economic and the ecological environment, improving water utilization efficiency, and realizing water resources consumption control under uncertainties. The method is incorporated with interval parameter programming, fuzzy programming, and chance-constrained programming, for handling system uncertainties and balancing the optimal objectives with the risk of violating system constraints. Under this framework, an IFCCP model for water resources allocation management was successfully formulated and applied to a typical water-deficit area, Tianjin, China, for obtaining a better water resources plan among multiple users under resources and environmental limitation. Different total water consumption control policies are designed for assessing regional water allocation schemes. The results indicated that the gap of supply and demand will only be solved by foreign water, the transferred water from Luan River and Changjiang River would still be the main supplier in planning horizon. Moreover, the strict total water consumption control policy would guarantee the water requirement of ecological environment, lead to changes in the structure of water supply, actively guide on water conservation, and promote the large-scale utilization of desalted water and recycle water.     

An inexact fuzzy-random-chance-constrained air quality management model

An inexact fuzzy-random-chance-constrained programming model (IFRCCMM) was developed for supporting regional air quality management under uncertainty. IFRCCMM was formulated through integrating interval linear programming within fuzzy-random-chance-constrained programming framework. It could deal with parameter uncertainties expressed as not only fuzzy random variables but also discrete intervals. Based on the stochastic and fuzzy chance-constrained programming algorithms, IFRCCMM was solved when constraints was satisfied under different satisfaction and violation levels of constraints, leading to interval solutions with different risk and cost implications. The proposed model was applied to a regional air quality management problem for demonstration. The obtained results indicated that the proposed model could effectively reflect uncertain components within air quality management system through employing multiple uncertainty-characterization techniques (in random, fuzzy and interval forms), and help decision makers analyze trade-offs between system economy and reliability. In fact, many types of solutions (i.e. conservative solutions with lower risks and optimistic solutions with higher risks) provided by IFRCCMM were suitable for local decision makers to make more applicable decision schemes according to their understanding and preference about the risk and economy. In addition, the modeling philosophy is general and applicable to many other environmental problems that may be complicated with multiple forms of uncertainties.     

An inexact inventory-theory-based chance-constrained programming model for solid waste management

In this study, an inexact inventory-theory-based chance-constrained programming (IICP) model is proposed for planning waste management systems. The IICP model is derived through introducing inventory theory model into a general inexact chance-constrained programming framework. It can not only tackle uncertainties presented as both probability distributions and discrete intervals, but also reflect the influence of inventory problem in decision-making problems. The developed method is applied to a case study of long-term municipal solid waste (MSW) management planning. Solutions of total waste allocation, waste allocation batch and waste transferring period associated different risk levels of constraint violation are obtained. The results can be used to identify inventory-based MSW management planning with minimum system cost under various constraint-violation risks. Compared with the ICP model, the developed IICP model can more actually reflect the complexity of MSW management systems and provide more useful information for decision makers.     

Inexact quadratic joint-probabilistic programming for water quality management under uncertainty in the Xiangxi River, China

An inexact quadratic joint-probabilistic programming model for water quality management (IQJWQ) is developed and applied to supporting multiple-point-source waste reduction in the Xiangxi River, China. The IQJWQ is a hybrid of interval quadratic programming, joint probabilistic programming and multi-segment water quality simulation. It has advantages in reflecting uncertainties expressed as joint probabilities of system risk, probability distributions of water quality standards, interval parameters and nonlinearities in the objective function. An interactive and derivative algorithm is employed for solving the IQJWQ model. The results indicate that the Pingyikou chemical plant and Liucaopo chemical plant contribute more to pollution of the main stream in the Xiangxi River, which should be the prior plants to reduce the wastewater discharge and enhance the wastewater treatment efficiencies. Meanwhile, the environmental agencies should choose the joint probability carefully to balance the tradeoff between production development and pollution control. Compared with the conventional chance-constrained programming method, the IQJWQ exhibits an increased robustness in handling the overall system risk in the optimization process. Although this study is the first application of the IQJWQ to water quality management, the proposed methods in the IQJWQ can also be applicable to many other environmental management problems under uncertainty.     

Many-year variations of water pollution and the state of river ecosystems in different latitudinal zones in European Russia

Many-year variations of water pollution and the state of more than two hundred of rivers (or their segments) in different latitudinal zones in European Russia are estimated. In accordance with latitudinal zonality, space and time trends in variations of water quality and the state of rivers are identified. By the level of water pollution, most river ecosystems (or their segments) were classified as polluted and very polluted, very polluted and dirty, dirty and very dirty. A tendency toward water quality improvement could be seen in the recent decade. More than 50% of the examined river ecosystems are at the transitional state of passing from natural (and/or equilibrium) to crisis state.     

Stochastic models of severe weather watches and warnings: transition probabilities

To alert the public to the possibility of tornado (T), hail (H), or convective wind (C), the National Weather Service (NWS) issues watches (V) and warnings (W). There are severe thunderstorm watches (SV), tornado watches (TV), and particularly dangerous situation watches (PV); and there are severe thunderstorm warnings (SW), and tornado warnings (TW). Two stochastic models are formulated that quantify uncertainty in severe weather alarms for the purpose of making decisions: a one-stage model for deciders who respond to warnings, and a two-stage model for deciders who respond to watches and warnings. The models identify all possible sequences of watches, warnings, and events, and characterize the associated uncertainties in terms of transition probabilities. The modeling approach is demonstrated on data from the NWS Norman, Oklahoma, warning area, years 2000–2007. The major findings are these. (i) Irrespective of its official designation, every warning type {SW, TW} predicts with a significant probability every event type {T, H, C}. (ii) An ordered intersection of SW and TW, defined as reinforced warning (RW), provides additional predictive information and outperforms SW and TW. (iii) A watch rarely leads directly to an event, and most frequently is false. But a watch that precedes a warning does matter. The watch type \(\{SV\), TV, \(PV\}\) is a predictor of the warning type \(\{SW\), RW, \(TW\}\) and of the warning performance: It sharpens the false alarm rate of the warning and the predictive probability of an event, and it increases the average lead time of the warning.     

Simulating and mapping the risk of surge floods in multiple typhoon scenarios: a case study of Yuhuan County,Zhejiang Province,China

We propose a scenario-based method for simulating and mapping the risk of surge floods for use by local authorities concerned with public safety and urban planning in coastal areas. Focusing on the triad of hazard, vulnerability and adaptation capability, we estimate the comprehensive risk and display its spatial distribution using the raster calculation tool in ArcGIS. The detailed methodology is introduced via a case study of Yuhuan, an island county in Zhejiang Province, China, which is frequently affected by typhoon storm surges. First, we designed 24 typhoon scenarios and modeled the flood process in each scenario using the hydrodynamic module of MIKE 21. Second, flood depth and area were used for hazard assessment; an authorized indicator system of land use categories and a survey of emergency shelters were used for vulnerability and adaptation capability assessment, respectively; and a quantified model was used for assessment of the comprehensive risk. Lastly, we used the GIS raster calculation tool for mapping the risk of storm surges in multiple typhoon scenarios. Our principal findings are as follows: (1) Seawalls are more likely to be overtopped or destroyed by more severe storm surges with increasing typhoon intensity. (2) Most of the residential areas with inadequate emergency shelters are highly vulnerable to flood events. (3) As projected in the risk mapping, if an exceptional typhoon with a central pressure of 915 or 925 hPa made a landfall in Yuhuan, a wide range of areas would be flooded and at high risk. (4) Determining optimal strategies based on identification of risk-inducing factors is the most effective way of promoting safe and sustainable development in coastal cities.     

An interval-parameter multi-stage stochastic programming model for water resources management under uncertainty

In this study, an interval-parameter multi-stage stochastic linear programming (IMSLP) method has been developed for water resources decision making under uncertainty. The IMSLP is a hybrid methodology of inexact optimization and multi-stage stochastic programming. It has three major advantages in comparison to the other optimization techniques. Firstly, it extends upon the existing multi-stage stochastic programming method by allowing uncertainties expressed as probability density functions and discrete intervals to be effectively incorporated within the optimization framework. Secondly, penalties are exercised with recourse against any infeasibility, which permits in-depth analyses of various policy scenarios that are associated with different levels of economic consequences when the promised water-allocation targets are violated. Thirdly, it cannot only handle uncertainties through constructing a set of scenarios that is representative for the universe of possible outcomes, but also reflect dynamic features of the system conditions through transactions at discrete points in time over the planning horizon. The developed IMSLP method is applied to a hypothetical case study of water resources management. The results are helpful for water resources managers in not only making decisions of water allocation but also gaining insight into the tradeoffs between environmental and economic objectives.     

Achieving the objective of ecological planning for arid inland river basin under uncertainty based on ecological risk assessment

To study the impact of the uncertainties existed in nature and management on the ecological planning and the achievement of the planning objective under uncertainty, the General restoration planning of Shiyanghe river basin, which is for Shiyanghe river basin located in an arid inland area in China, is selected as the research object. The general thought of the essay is to get the probability of the objective achievement of ecological planning by analysing the floating region and possibility distribution of the uncertainties which is of the critical factors in achieving the planning objective. Firstly, the uncertainties of the critical factors are identified and analysed separately. Secondly, the two uncertainties are combined to study its impacts on the achievement of planning objective in each scenario. Lastly, the probability and the possible deviation of ecological objective achievement are analysed in each scenario. The methods used in the study include the scenario analysis, BPF based on BP neural network, PPCC analysis and interval analysis. The results show that the achievement of ecological objective is greatly influenced by the uncertainties, the objective achievement could be get only in the scenario of water resources utilization and management meet the planning requirement and high flow year in planning year, the results also show that the greater volume of local water resources and the higher level of local water resources utilization and management lead to the higher probability of achieving ecological objective and the lower possible deviation.     

Multi-objective optimization for sustainable groundwater resource management in a semiarid catchment

Abstract

Groundwater is an important water resource and its management is vital for integrated water resources development in semiarid catchments. The River Shiyang catchment in the semiarid area of northwestern China was studied to determine a sustainable multi-objective management plan of water resources. A multi-objective optimization model was developed which incorporated water supplies, groundwater quality, ecology, environment and economics on spatial and temporal scales under various detailed constraints. A calibrated groundwater flow model was supplemented by grey simulation of groundwater quality, thus providing two lines of evidence to use in the multi-objective water management. The response matrix method was used to link the groundwater simulation models and the optimization model. Multi-phase linear programming was used to minimize and compromise the objectives for the multi-period, conjunctive water use optimization model. Based on current water demands, this water use optimization management plan was able to meet ecological, environmental and economic objectives, but did not find a final solution to reduce the overall water deficit within the catchment.     

Effect of zonal E × B plasma drift on electron density in the low-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at a solar activity maximum near vernal equinox

The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N _e and N _i ) and temperatures (T _e and T _i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N _e , N _i , T _e , and T _i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.     

Detailed Study of Time Variations in the Gutenberg–Richter <Emphasis Type="Italic">b</Emphasis>-value Based on Highly Accurate Seismic Observations at the Garm Prognostic Site,Tajikistan

The time variations in the Gutenberg–Richter b-value are minutely studied based on the data of highly accurate seismological observations at the Garm prognostic site, Tajikistan, where a stationary network of seismic stations of the Complex Seismological Expedition (CSE) of Schmidt Institute of Physics of the Earth (IPE) of the USSR (Russian) Academy of Sciences was in operation from 1955 to 1992. A total of 93035 local earthquakes ranging from 0.0 to 6.3 in the Ml magnitudes are considered. The spatiotemporal fluctuations in the minimal magnitude of completeness of the earthquakes, Mc, are analyzed. The study considers a 25-year interval of the observations at the center of the observation system within which Mc = 0.9. It is shown that in most cases, the b-value and log₁₀E^2/3 experience characteristic time variations before the earthquakes with magnitudes higher than the minimal magnitude of the predicted earthquake (MPE). The 6-year anomaly in the parameters’ b-value, log₁₀E^2/3, and log₁₀N associated with the single strongest earthquake with M = 6.3 that occurred in the observation region on October 26, 1984 is revealed. The inversely proportional relationship is established between the time variations in the b-value and the time variations in the velocities of seismic waves Vp and Vp/Vs. It is shown that the exponent p in the power function which links the time variations of the b-value and log₁₀E^2/3 is higher in the zones of crustal compression than in the zones of extension. It is simultaneously confirmed that the average b-value in the zones of compression is lower than in the zones of extension. It is established that in the case of earthquakes with M ≥ 2.6, the time series of seismic activity log₁₀N_i and the time series of the b-value are highly cross correlated with a coefficient of r ≈ 0.75, whereas in the case of earthquakes with M ≥ 0.9, the coefficient of cross correlation between these time series is close to zero (r ≈ 0.06). The law of variations in the slope of the lines approximating the relationship between the log₁₀N_i time series in the different magnitude ranges (M ≥ Mc_i) and b-value time series is obtained. It is hypothesized that the seismic activity of the earthquakes with high magnitudes can be estimated provided that the parameters of the time series of the b-value and time series of the number of earthquakes logN(ΔM_i) in the range of low magnitudes are known. It is concluded that using the parameter log₁₀N for prognostic estimates of the strong earthquakes only makes sense for earthquakes having moderate and large magnitudes. It is inferred that the time variations in the b-value are predominantly contributed by the time variations of the earthquakes with relatively large magnitudes.     

Dual inexact fuzzy chance-constrained programming for planning waste management systems

This study develops a dual inexact fuzzy chance-constrained programming (DIFCCP) method for planning municipal solid waste (MSW) management systems. The concept of random boundary interval (RBI) is introduced to address the high uncertain parameters in the studied system. Fuzzy flexible programming and chance-constrained programming are also introduced to take into account the uncertainties of RBIs and various uncertainties in MSW management system. Compared with the existing methods, the developed method could deal with the uncertainty without simplification and thus is more robust. Moreover, the potential system-failure risks in MSW management system due to the existing uncertainties could be quantified by means of violation levels and satisfaction levels in DIFCCP. The developed method then is applied to a MSW management system. The obtained solutions could be used for generating efficient management schemes. The values of violation and satisfaction levels could help decision makers understand the tradeoffs between system cost and system-failure risk, and identify desired strategy according to the practical economic and environmental situation.     

Effects of temporally correlated infiltration on water flow in an unsaturated–saturated system

Effects of temporally correlated infiltration on water flow in an unsaturated–saturated system were investigated. Both white noise and exponentially correlated infiltration processes were considered. The moment equations of the pressure head (ψ) were solved numerically to obtain the variance and autocorrelation functions of ψ at 14 observation points. Monte Carlo simulations were conducted to verify the numerical results and to estimate the power spectrum of ψ (S _ψψ ). It was found that as the water flows through the system, the variance of the ψ (\( \sigma_{\psi }^{2} \)) were damped by the system: the deeper in the system, the smaller the \( \sigma_{\psi }^{2} \), and the larger the correlation timescale of the infiltration process (λ _I ), the larger the \( \sigma_{\psi }^{2} \). The unsaturated–saturated system gradually filters out the short-term fluctuations of ψ and the damping effect is most significant in the upper part of the system. The fluctuations of ψ is non-stationary at early time and becomes stationary as time progresses: the larger the value of λ _I , the longer the non-stationary period. The correlation timescale of the ψ (λ _ψ ) increases with depth and approaches a constant value at depth: the larger the value of λ _I , the larger the value of λ _ψ . The results of the estimated S _ψψ is consistent with those of the variance and autocorrelation function.     

Identification of optimal strategies for agricultural nonpoint source management in Ulansuhai Nur watershed of Inner Mongolia,China

In this study, an environmental-friendly modeling system was developed and applied to an agriculture nonpoint source (AGNPS) management in Ulansuhai Nur watershed. In this system, water environmental capacity, credibility-based chance-constrained programming (CCCP), and AGNPS optimization models were integrated into a general modeling framework. It could be used to calculate water environmental capacity of total nitrogen and total phosphorus in Ulansuhai Nur watershed, which could consequentially provide input data for the developed AGNPS optimization model. Also, the inherent uncertainties in estimating water environmental capacities that can be expressed as possibilistic distributions were reflected and addressed based on computational results of three widely used methods. Such uncertainties were consequentially transferred to the proposed CCCP model based on the adoption of multiple credibility satisfactory levels, significantly facilitating objectivity reflection of decision alternatives. The developed modeling system was then applied to Ulansuhai Nur watershed of Inner Mongolia, a semi-arid river basin in northwestern China. Optimal strategies for AGNPS management in Ulansuhai Nur watershed were generated with consideration of the maximum total agricultural income under multiple policy scenarios. The results showed that the total agricultural income would increase with point source pollution being cut down, and would decrease with rising credibility levels, representing decreasing system violation risks. It was indicated that the higher of total nitrogen/phosphorus discharge being less than water environmental capacity of Ulansuhai Nur, the lower the total agriculture incomes. The proposed methods could help decision makers establish various production patterns with cost-effective agriculture nonpoint source management schemes in the basin of Ulansuhai Nur, and gain in-depth insights into the trade-offs between total agricultural incomes and system reliabilities.