A Sustainable Industry-Environment Model for the Identification of Urban Environmental Risk to Confront Air Pollution in Beijing, China

In this study, an optimized industry-environment model is proposed for identifying environmental risk under uncertainties. The strategy associated with an emission-permit trading mechanism has been introduced into the industrial-environment regulation (model) for remitting the pressures of frequent/severe haze events in Beijing City. A dual stochastic mixed fuzzy risk analysis method with Laplace’s criterion (DSFRL) can be embedded into industry-environment issues with a trading emission-permit trading mechanism (IEST) for handling uncertainties regarded as possibility and probability distributions. Meanwhile, this can also reflect the environmental risks and corresponding system benefits due to the occurrence of a random event (such as random wind velocity). Based on the application of the proposed IEST with DSFRL, the numbers of the obtained results associated with production reduction, adjustment of industrial layout pattern, emission-permit transactions, pollutant mitigation and system benefits under various Laplace criterion cases can be analyzed. A tradeoff between production development and pollution mitigation based on the preference of policymakers can be used for rectifying current strategies with a sustainable mode, which can prompt an effort to confront air pollution


Introduction
Against the background of rapid urbanization/industrialization, high-speed economic development has greatly improved human lives across the world, but has also brought about detrimental pressures on the air environment. Particularly in a number of fast-growing cities such as Beijing, China, the accelerated exploitation of industry and disorderly industrial layouts can generate excessive pollutant emissions, which exceed what the environmental load can afford, generating environmental risks. From 2013, environmental risks associated with haze events of atmospheric pollution have been frequent. Given the hazed concentration of fine particulate matter (e.g., PM 2.5 ) and traditional pollutants (including SO 2 , NO 2 , and PM 10 ), the incidence of respiratory system diseases due to air pollution has presented an increasing trend. Meanwhile, other damage to cardiovascular/cerebrovascular systems and nervous systems in terms of public health have been reported frequently in response to air pollution. Thus, the air pollution issue can give rise to a health risk in Beijing, which is today challenging policymakers [1][2][3]. Against the stresses of air pollution, a number of anthropogenic engineering techniques (including engineering control measures, clean production techniques, environmental regulation policy and an industrial adjustment strategy) have been advocated to confront such hazardous air pollution issues. Among them, pollutant emission-permit trading is an effective method which can facilitate pollutant permits from a lower to higher value in order to achieve gains in economic efficiency, but also can provide incentives to adopt pollution abatement measures [4]. Meanwhile, a market approach can drive industrial companies to improve the efficiency of pollutant treatments, with the aim of avoiding a higher environmental penalty or cost for buying the emission permit. Moreover, a trading mechanism-based (TM) can adjust the scale and structure of emissions to improve the efficiency of production at the extreme, leading to greater system benefit. Therefore, an optimized industry-environment strategy with an emission-permit trading mechanism (IEST) is desirable for coordinating the relationship between the development of the economy and the protection of the environment in Beijing, reducing the environmental risk and damage to public health [5,6].
However, an IEST issue is complicated with a number of objective and subjective uncertainties, leading to complexities. Above all, uncertainties and their corresponding complexities can increase the risk for policymakers to generate the desired decision alternatives [7]. Previously, many research works such as genetic algorithms, potential economic impact assessments, the asymmetric fuzzy method and the inexact stochastic programming model have been developed for handling inherent uncertainties in the decision process of an IEST issue, which can reduce the difficulties and risk levels of decision-making [8][9][10][11][12][13]. Among them, two-stage stochastic programming (TSP) is a useful method to deal with random problems through rectifying actions (i.e., recourse action with probabilistic event). For example, meteorological conditions in an IEST issue can influence the diffusion and migration processes of pollutants, which can result in random concentration levels. If these levels exceed standards regulated by a government, this implies the random events of air pollution may occur at several times or last a few days. These random hazed air pollution incidents can bring about environmental loss to the first-stage decisions (e.g., economic development targets) periodically over time through TSP [14]. However, this cannot reflect the variability of the second-stage decisions, where a high-variability level of risk violation would challenge policymakers [15]. Therefore, a risk control analysis method (RAM) is proposed to capture variable risk in a TSP contest, which can control the variability of the recourse values (e.g., pollutant penalty) in a rational range [16]. Nevertheless, in an IEST issue, the economic and environmental strategy with scenario assumptions can be also expressed as probabilistic distributions, which result in dual stochastic situations. A stochastic Laplace criterion (SLC) is introduced to handle the probability of each scenario occurrence under the supposition of no data being available on the probabilities of the various outcomes [17,18]. Meanwhile, a number of fuzzy occurrences in an IEST issue (such as missing economic data, an imprecise diffusion process and fuzzy trading efficiency) can intensify the risk of the decision-making process. For example, in an IEST issue, the diffusion or migration process of pollution is regarded as fuzziness since it is often difficult to calculate accurately. Therefore, a fuzzy credibility constraint programming (FCP) can be included while handling of this precise information in goals or constraints with a high satisfaction degree [19,20]. Previously, few works have focused on multiple uncertainties in hybrid formats in an IEST planning issue.
Thus, the objective of this study is to propose an optimized industry-environment strategy with an emission-permit trading mechanism (IEST) for alleviating production-emission conflict in Beijing. A dual stochastic mixed fuzzy risk analysis method with Laplace's criterion (DSFRL) can be embedded into an IEST issue, which can deal with fuzzy and dual stochastic uncertainties in a two-stage context. Meanwhile, it can also reflect and control the risk of a two-stage decision in a rational range according to risk exposure and aspiration benefit levels. The application in Beijing has shown that market trading is an effective mechanism to confront severely diminished environmental capacity and aggravating air pollutant emissions. Results of emission-permit transactions, production reductions, pollution-abatement schemes, adjustment of industrial structures and system benefits under various cases are obtained and analyzed. The results show that a trading mechanism is a more effective method to adjust industrial layout and reduce air pollutant emission in the study region, as well as manage the air-environment crisis. Meanwhile, the risk-averse attitude and robustness coefficient can improve the reliability of decision-making in an IEST issue. These findings can support policymakers in identifying optimized industry-environment policies for coordinating the relationship between economic development and environmental protection, as well as controlling the air pollution crisis at an urban level.

Two-Stage Stochastic Programming for Risk Control (TSR)
In a decision-making problem, two-stage stochastic programming (TSP) can built a linkage between expected targets and a random event, where the first-stage benefit can be rectified by a second penalty as follows [21]: In Model 1, a recourse action occurs between uw (i.e., first-stage benefit) and second-stage penalty), where δ h is random event after the first-stage decision has been made; and the p h denotes as the probability of realization of event δ h . However, TSP can not control the risk of recourse action against any infeasibilities under various scenarios (h). If an excessive penalty exceeds what the system can afford, it can bring about systemic failure. Therefore, a risk-control method can be joined into a two-stage context [16]: where if the actual benefit (i.e., R h (w)) is smaller than the expected target, the probability (i.e., P h ) that presents the targeted benefit cannot be met under various scenarios h would be 1, otherwise, it would be 0 [18]. Thus, a risk control based TSP can be expressed as follows: where X h is a binary variable defined for each scenario; η is the allowable risk exposure level; ε is target level; RB h is the balancing benefits for scenario h and target level (ε). In fact, the probability of not meeting the targeted benefit in each scenario would be either 0 or 1.

Dual Stochastic Programming for Risk Control with Laplace's Criterion Scenario Analysis (DSLS)
However, in a practical IEST issue, the input of the expected target (i.e., the first-stage decision variable) can be influenced by various factors, which would be tackled by a scenario analysis (SA) [20]. In general, the risk attitudes of policymakers can impact the generation of a scenario due to their experiences and personality traits. Thus, the definition of a stochastic scenario analysis (SSA) can be expressed as follows: where E(A mn ) is the decision outcome; A n is payoff matrix row, (A n ∈ A, n = 1, 2, . . . N); P m is probability of each scenario occurrence; d is the option; D is the options; R mn is the overall performance; PS m 1 n is the progressive scenario, which means that progressive target with optimism attitude; CS m 2 n is the conservative scenario, which means that conservative target with pessimism attitude; m 2 is the numbers of scenarios; among them, m 1 is the number of progressive scenarios; (m 2 − m 1 ) is the number of conservative scenarios. However, since the scenario of each scenario is random, this results in policymakers having difficulty in decision making. Laplace's criterion can handle the probability of each scenario occurrence under the supposition of no data available on the probabilities of the various outcomes, which can be formulated as follows [17,18]: Based on TSR, dual stochastic programming for risk control with Laplace's criterion scenario analysis (DSLS) can be modified as follows: 11 PS 12 · · · PS 1n PS 21 PS 22 · · · PS 2n · · · · · · · · · · · · PS m 1 1 PS m 1 2 · · · PS m 1 n Subject to:

A Dual Stochastic Mixed Fuzzy Risk Analysis Method with Laplace's Criterion (DSFRL)
Nevertheless, a number of uncertainties related to precise estimation and fuzziness cannot handle probabilistic distributions [21]. Therefore, fuzzy credibility constrained programming (FCP) can be advocated as follows: Based on the concept of fuzzy credibility, the credibility measure (Cr) can be expressed as Cr{ς ≤ s} = 1 2 (Pos{ς ≤ s} + Nec{ς ≤ s}) [22]. In general, the credibility level should be greater than 0.5 usually [23]. Thus, the model (6a) can be proven as follows: Thus, a dual stochastic mixed fuzzy risk control method (DFR) can be resolved as follows:

Overview of the Study Area
Beijing is located on the North China Plain, which covers an area of 16,807.8 km 2 . The city has 18 districts, with a population of more than 21.71 × 10 6 people (at the end of 2015). As the capital of China, it has grown fast in recent years. In 2015, the GDP reached 2296.8 billion Yuan, and it has maintained a high growth rate of about 6.90% in recent years. Meanwhile, the growth rate of the population has been maintained at 9.00‰ per year [24,25]. In the context of rapid urbanization/industrialization, industry has developed rapidly. The establishment of a modern agricultural structure, developing a peri-city modernistic industry, and prioritizing tourism have promoted living standards and social material productivity [24]. Due to the expansion of industry, large-scale pollutant discharges have exceeded what the natural environment can afford, which results in frequent haze events of air pollution. Figure 1 displays the framework of a risk-control based pollutant emission management model with Laplace's criterion (RPEL) and its application in the study region. In order to confront the air pollution issue, policymakers have built a linkage between "production" and "emission" for the adjustment of production, which can reduce pollutant emissions. In a traditional urban industry-environment model, enterprises would reduce the scale of production to confront the environmental penalty. However, a single political plan cannot conserve the environment at the extremes due to "governance failure". Thus, a market approach is joined with industrial-environment regulation to reallocate emission permits from lower to higher values. This can increase economic efficiency overall, meanwhile providing incentives to adopt pollution abatement. In an optimized industry-environment strategy with an emission-permit trading mechanism (IEST) issue, the regional policymakers must create a plan to allocate emission permits to every enterprise and maximize the overall system benefit. Based on an initial emission obligation, the enterprise can buy or sell the emission permit in the market, remedying their emission deficit or releasing their excessive permits. Otherwise, excessively polluting enterprises have to pay the higher environmental penalty due to their emission-permit deficit. The trading mechanism can facilitate enterprise to improve the efficiency of pollutant retreatment and productivity, which can be beneficial for confronting the higher environmental penalty. The tradeoff between economic development and emission abatement can generate a comprehensive plan associated with pollution control. However, in an IEST issue, there are a variety of factors such as on-site surveying and monitoring, the analysis of the main affected factors, the determination of pollution-source emission standards, the partition of functional zones, and the design of their respective environmental capacities, as well as the generation of control measures [26]. Thus, a dual stochastic mixed fuzzy risk analysis method with Laplace's criterion (DSFRL) is proposed to handle these objective and subjective uncertainties in an IEST system. Therefore, the steps of the RPEL method applied in an IEST issue in the study region can be formulated as follows: Step 1: Recognition of the regional features of economic development, production scale, pollutant emission and meteorological conditions in Beijing.
Step 2: Construction of a linkage of production-emission to reflect the relationship between production scale/structure and pollutant emission in the study region.
Step 3: Establishment of an emission-permit trading mechanism to reallocate the emission permit from lower to higher value, regulating the regional production layout and emission scale.
Step 4: Identification of various uncertainties and their interactions in an IEST issue in Beijing.
Step 5: Incorporation of two-stage stochastic programming, risk-analysis control method, Laplace scenario and fuzzy credibility programming into a framework to propose a RPEL method to deal with multiple uncertainties in an IEST issue.
Step 6: Formulation of a sustainable industry-environment model for identifying environmental risk under uncertainties and meanwhile calculating optimized results. The model has been calculated by Lingo 12.0, that has the advantage of providing a capacity to deal with unlimited variables and constraints in a friendly way compared to ORS, WinQSB and Matlab. The Lingo 12.0 is an effective tool for policymakers to describe the complex system and interactions among components based on a series of equations and functions.
Step 7: Analysis of obtained results associated with production reduction, adjustment of industrial layout patterns, emission-permit transactions, pollutant mitigation and system benefits under various Laplace criterion cases.
Step 8: Generation of a sustainable industry-environment plan to adjust regional strategies to conform integrity between the regional economy and environment. Complete a plan process.

Modeling Formulation
Based on the DSFRL method, an optimized industry-environment strategy with emission-permit trading mechanism (IEST) is formulated. The policymakers should consider the issues as follows: (a) how to balance the relationship between economy and environment to confront urban diseases and the air crisis; (b) how to introduce a trading mechanism to improve the production mode and efficiency of pollutant mitigation; (c) how to adjust the current industrial production scale and pollutant emission permit based on a trading mechanism to meet air-quality standards in the study region; (d) how to generate a reliable plan against environmental risk. Thus, an optimized industry-environment with emission-permit trading mechanism (IEST) through DSFRL method can be formulated as follows: (1) Benefits from various industrial sectors based on expected economic development (IIP): (2) Environmental penalty for excessive emission from industrial sectors (PEP): (3) Loss for reduced production according to adjustment of industrial structure (LRP): (4) Environmental retreatment cost of various industrial sectors (CER): System benefit from emission-permit transaction (BET): Cost for emission-permit transaction (CET): The detailed nomenclatures for the variables and parameters is provided in the notation and subscript. Where f Laplace is the green gross domestic product (GGDP) in Beijing, which equals the total system benefit of economy-added return from emission-permit trading minus the loss from reduced production, the environmental penalty for excessive emissions, retreatment costs and trading costs. In this study, industry is deemed to be the main source of pollutant emissions, which can be divided into three sources whereby high, medium and low pollution industries can be considered in the modelling formulation. Meanwhile, industrial productiveness, pollutant generation and retreatment efficiency can be included in the model to reflect various pollution patterns. Under these situations, varied constraints such as mitigation requirements, capacity of pollutant-mitigation, emission allowance, air quality, industrial structure and scale and non-negativity can be displayed as follows: (1) Pollutant-mitigation requirement: (2) Pollutant-mitigation capacity: (3) Available emission permit for transaction: (4) Ambient air-quality requirement: (5) Industrial production scale: W IL tij · l p tij ≤ IPL tijmax (9g) (6) Economic risk based on DSFRL method: Constraint (1) displays the demand of pollutant-mitigation in an IEST issue, where the pollutant emission from the industrial sector should be restricted by the regional environmental load. Constraint (2) shows the capacity of pollution mitigation according to current clean production techniques, where the discharged pollution can be processed with meteorological data (e.g., wind velocity). Constraint (3) presents the emission-permit transaction, which equals the expected pollutant emission minus maximum emission-permit allowance. This means that after the emission-permit is allocated to enterprises, if they discharge more pollutants that exceed their initial permit obligations they have to buy a permit from the market; otherwise, they have to afford the environmental penalty. Constraint (4) displays the ambient air quality based on the Gaussian dispersion model (GDM), where the wind velocity (i.e., s) expressed as fuzziness can be handled by credible measures (as shown in the "Methodology" section). Constraint (5) shows the scale of industrial development, which can be restricted by the minimum/maximum development scale in the study region. Constraint (6) presents economic risk based on the DSFRL method, which can reflect the risk between the expected benefit and the targeted value (as shown in the "Method" section). Constraint (7) is non-negativity restrictions.

Data Collection
A number of input parameters are calculated by government reports, statistical yearbooks, and related research work [24]. Table 1 displays economic data, where the net benefit of unit production and corresponding environmental penalty of excessive pollutant emission have been listed. They can be estimated according to the statistical yearbook of Beijing from 2000 to 2013, with consideration of the development of the economy [24,25]. Meanwhile, since the meteorological conditions (e.g., wind velocity) can influence the capacity of pollution diffusion, varied meteorological data (i.e., wind velocity) are provided in Table 2. The data can be obtained from the weather stations of Beijing. Through simulation and calculation, the levels of wind velocity and corresponding probabilities have been listed in Table 2. In the study region, the meteorological condition can impact the capacity of the diffusion of pollutant, which can be calculated by the Gaussian dispersion model (GDM) (as shown in "Modelling Formulation"). The coefficient of the Gaussian dispersion model can be estimated according to the Pasquill-Gifford curve [26]. In this study, the regulation of the concentrations of pollutant emission from the industrial section can be obtained from the ambient air-quality standard (GB 3095 2012) [27] in China, which can be denoted "AAQS". Since Beijing confronts great air pollution stress, the air-quality requirement can be advanced to regional policymakers. From 2013, the Grade-II standard of the ambient air quality had been regulated, where the average maximum concentration limits of SO 2 , NO 2 , PM 10 and PM 2.5 were 60, 40, 70 and 35 µg/m 3 per year, respectively. Based on these regulations of AAQS in Beijing, we have obtained the right-hand side of the constraint (4), where the future air-quality requirement can be calculated in the fuzzy credibility manner: (i.e.,  Table 3 presents various cases associated with different risk attitudes of policymakers for developing the economy and protecting the environment, where 10 scenarios are designed to reflect different risk attitudes (Laplace scenarios).  Figure 2 shows total excessive pollutant emissions under various cases when ε and α levels are 0.95 and 0.99, respectively. In this study, excessive pollutant emissions can bring about an increment of concentration of pollution in the air, which would influence public health. The results indicate that lower wind velocities would generate higher concentrations of pollution. For example, when wind velocity is at a high level, the ultra-concentration ratios of PM 10 and PM 2.5 were 3.23 and 0.52 times than when wind velocity is at a low level under case 1 in period 1. Meanwhile, the results show the highest ultra-concentration ratio of SO 2 and NO 2 would be 7.34 and 6.23 times compared to normal conditions due to bad meteorological conditions for diffusion. This demonstrate that the current industrial production mode (including scale of production, emission pattern and retreatment) is backward and cannot accommodate regional environmental regulation. Thus, adjustment of the pollution mitigation scheme, industrial structure and clean production technique should be required.    Figure 2, a number of companies have obtained permits. However, their pollutant emissions have exceeded the allowance. Therefore, they have to reduce the scale of production. The results obtained from Figure 3 presenting the scales of reduced production would vary with different pollutant diffusion capacity such as wind velocity. In Beijing, higher wind velocity occurs in summer and autumn, which would be suitable for pollutant diffusion. Under these situations, the scales of reduced production would be relatively lower due to its higher capacity of diffusion, and vice-versa. Meanwhile, various production scales in different industrial sectors can generate varied reductions. In comparison, the industrial sector with the feature of higher pollutant emission and lower economic returns would be restricted first. For instance, iron processing (denoted as "IP") with features of higher pollutant emission and lower economic returns can generate the highest reduction. This implies that current scale of IP and corresponding production mode would bring about and exorbitant environmental penalty, leading to reduced benefits. Thus, adjustment of the IP scale or the improvement of its production and emission modes would be required. Similar situations occurred in petroleum refining (denoted as "PR") and power generation (denoted as "PG") in response to their higher production scale, pollutant emissions and lower economic benefits. For instance, when wind velocity is at a low level in period 1, the reduced scale of production in PR, PG and IP would be 445.2 × 10 6 tons, 244.3 × 10 9 kwh and 594.4 × 10 6 tons under case 1 (α = 0.95). By contrast, the reduced scale of vehicle manufacturing (denoted as "VM") and chemical industry (denoted as "CI") are relatively lower than in PG, PR and IP due to their lower scales of pollutant emission. This implies that their scales of production are relatively rational, since their scales of production were within the limits of environmental capacity.    Figure 6 presents the system benefits under cases 1, 5, 7, 9, 10 and an optimal state when α level is varied (ε = 0.95). The results indicate that a lower α level can result in a lower benefit and vice versa. Meanwhile, varied cases can bring about various system benefits as follows: (a) the conservative case (i.e., case 5) can generate lower environmental penalties and retreatment costs, since it require higher environmental protection than economic development; (b) the progressive cases (i.e., cases 7 and 9) can generate the opposite results compared to the conservative cases. Their higher target requirements for the economy would bring about higher environmental penalties; (c) in comparison, the basic case (i.e., case 1) can generate higher benefits than cases 5, 7 and 9. For instance, when α level is 0.9, system benefit would be RMB ¥3.87 × 10 12 under case 1, while system benefits are RMB ¥3.45 × 10 9 and $1.05 × 10 9 under case 5 and case 7; (d) the Laplace scenario (case 10) can generate the highest system benefit since it can generate a reliable result according to an element of risk. This implies that the current plan (i.e., case 1) should be adjusted to approach the Laplace scenario (case 10), achieving an optimal state. Figure 6. System benefits under cases 1, 5, 7, 9, 10 and optimal state when α level is varied (ε = 0.95). Figure 7 displays the solutions for excessive NO 2 emissions under case 1 with trading and non-trading mechanisms when ε and α levels are 0.6 and 0.5, respectively. The results show that a trading-oriented mechanism is superior to non-trading over the planning horizon. In a trading mechanism, an emission permit can be encouraged from low-emission to high-emission industrial sectors. Since a buyer should pay too much money for excessive emissions, the buyer should consider improving the industrial productiveness or retreatment efficiency instead, which can prompt emission mitigation. Under this situation, the pollutant emission under the trading mechanism is lower than that under the non-trading mechanism. For example, the NO 2 emissions with the trading scheme would be 157 mg/m 3 , which is lower than that with the non-trading scheme (197 mg/m 3 ). Under same policy scenario, the trading-oriented mechanism can improve emission mitigation, which can reduce the environmental penalty of excessive pollutant emission, leading to higher system benefits.

Discussion
Although a trading mechanism can prompt emission mitigation, reduced production can cut down the economic benefit from industrial sectors. Therefore, the optimized system benefit can be calculated based on the trade-off between an industrial economy and an environmental penalty. Comparing system benefits with trading and non-trading mechanisms, the results show that the trading approach is a more effective method than a non-trading mechanism for an industry-environment plan, which can increase and improve economic returns on the whole (as shown in Figure 8). Emission permits can be encouraged to move from a low value to a high value in a market, which can generate higher benefits than with a non-trading approach. For example, under case 1, when η is 0.9, system benefits would be RMB ¥5.42 × 10 12 with a trading mechanism, while system benefits are RMB ¥4.58 × 10 9 without a trading approach. Meanwhile, the system benefit can vary with α level.  The results obtained can bring about the following findings: (a) unreasonable industrial structure and excessive production has enhanced the occurrences of hazed air pollution. A laissez-faire policy leading to the overdevelopment of industry would lead to excessive pollutant emissions exceeding the atmospheric load. Therefore, improving the productivities of production and outdated retreatment techniques would aggravate pollution issues. (b) A trading mechanism is a more effective manner of improving the efficiency of pollutant mitigation on the whole. Therefore, improving market construction/market behavior can improve the efficiency of emission permit trading to coordinate industrial production and environmental protection. (c) An unscientific risk choice (including confidence levels and scenario planning) in a decision process with uncertainties can affect an IEST plan, which adversely affects the policymaker leading to neither adventurous nor conservative decisions.

Conclusions
In this study, a dual stochastic mixed fuzzy risk analysis method embedded into an optimized industry-environment strategy with an emission-permit trading mechanism (IEST) is proposed for remitting the pressures of frequent/severe haze events in Beijing in conditions of uncertainty. The novelty of the developed method is incorporating TSP, SLC, and FCP into a general framework to deal with multiple formats of uncertainty in a practical IEST issue, which can prompt a capacity to handle complexities and uncertainties. The developed method couples risk analysis, optimization and strategy adjustment into a framework to reflect objective uncertainties (e.g., meteorological conditions, economic development and environmental pollution) regarded as possibility and probability distributions. Meanwhile, it can reflect subjective uncertainties such as the risk attitudes of policymakers in a practical IEST issue. This has assessed the risk based on desired risk exposure level and aspiration benefit level within a two-stage context, with the aim of controlling the environmental penalty in a rational range. Moreover, it can optimize the industry-environment plan with a trading-mechanism, which can prompt pollutant emission from low to high values, encouraging pollutant mitigation.
Although the developed DSFRL method applied to an IEST issue can effectively deal with multiple uncertainties and identify environmental risk under uncertainties. It can help authorities to assess and manage urban environmental risk in order to obtain sustainable industry-environment solutions. However, in an industry-environment system, the environmental loading for various pollutant emissions is difficult to calculate accurately over an extensive area with sporadic time frames. It should take account of various diffusion mechanisms and characterize the spatial and temporal variations of meteorological conditions for emission permit trading. Coupling simulation and a genetic algorithm can reflect, above all, the uncertainties and their interactions in an effective manner in further research.