Identifying Key Sources of City Air Quality: A Hybrid MCDM Model and Improvement Strategies

Improvements for air quality were prioritized according to gaps between criteria. Existing improvement strategies for air quality often focus on a single criterion, without considering associations among criteria. Moreover, solutions are often temporary, failing to provide long-term improvement. Therefore, this study employed a multiple-criteria decision-making model and a substitution method based on causal relationships to analyze potential improvement strategies for air quality in Kaohsiung, Taiwan. Results revealed that coal-fired power plants and factory emissions are the major sources of pollution in Kaohsiung. This study suggests that environmental authorities in Kaohsiung should facilitate plants to improve energy efficiency with anti-pollution facilities.


Introduction
The greenhouse effect continues to make the environment deteriorate, and air pollution has become a topic of concern for everyone. Air pollution increases social and economic costs and leads to shortened life expectancy [1,2]. According to a World Health Statistics report [3], the number of premature deaths caused by air pollution is twice as high as a previous estimate, and air pollution has caused the death of nearly seven million people worldwide, becoming the greatest health hazard for human beings [4][5][6][7]. Air pollution has been classified as a major environmental carcinogen, with the most resultant deaths in developing countries in the Asia-Pacific region. Even at low concentrations, fine particles can have adverse health effects [8][9][10].
According to the air quality monitoring report of the Environmental Protection Administration (EPA) of Taiwan's Executive Yuan in 2018 [11], the annual average concentration of fine particulate matter (PM 2.5 ) in Taiwan is 2.5 times higher than the standard limit suggested by the World Health

Stationary Air Pollution Sources
The air pollution generated from existing coal-fired power plants has a significantly negative effect on human health and the environment and necessitates upgrades for existing power plant equipment [16,17]. The emission concentration of pollutants dramatically decreased because of flue gas desulfurization facilities installed at power plants [18].
Through investments in the reduction of pollutant emissions, boiler room parameters have been changed to correspond to environmental regulations [19]. Organic pollutants in machine-dismantling factories have indicated potential harmful effects in human and other in vivo studies [20]. Incidents of allergic symptoms have been significantly higher in children living in high-pollutant-exposure areas than those in the low-exposure areas [21].
Farmers burn crops to prepare biochar to increase soil fertility [22]. Field incineration of crop residues is a common practice in many parts of the world, with potential effects on air quality, the atmosphere, and climate [23]. The main pollutant in agricultural activities is ammonia, which is mainly derived from livestock farming and related manure treatment. Agricultural smoke emissions lead to PM 2.5 , which has a major influence on human health [24].
A large amount of dust is generated during the construction of big projects. Because it is difficult to disperse, it poses a serious threat to human health [24,25]. Determining the source of urban smog is complex; some of the main causes are rapid urbanization and increasing dust emissions [26,27]. Urbanization has led to a significant increase in waste disposal, and the ash produced by incinerators has increased dramatically [28]. The Chinese method of cooking food gives rise to a large number of carcinogenic soot chemicals [29]. Volatile organic compounds (VOCs) are released from the exhaust pipes of restaurants [30].
Joss paper burning is another source of toxic substances [31,32]. Burning joss paper is a common tradition at temples in Asian countries, and this tradition is used for ceremonial purposes in the The purpose of this section is to offer an overview of the whole analytical process, since it was combined with several procedures and methods. This study constructed a hybrid MCDM model. The decision-making trial and evaluation laboratory (DEMATEL) technique is used to build an influential network relation map (INRM); through the analytic network process (ANP), this DEMATEL-based ANP (DANP) was used to obtain influential weights (IWs). Finally, Vlse Kriterijuska Optimizacija Komoromisno Resenje (VIKOR, "Optimization of Compression Resistance"), combined with IWs to determine how to reduce gaps to achieve an ideal level, were used to assess the effectiveness of the enforcement of air pollution improvement by the Environmental Protection Bureau of Kaohsiung City Government (see Section 4.2. for a detailed description of this location and data collection).
The hybrid MCDM model used in this study is the concept of "gap" [49]. This model believes that people will naturally pursue a level of aspiration, which is the concept of a perfect score. As a result, decision makers can clearly understand the gaps that need to be improved for each indicator and the need to achieve continuous improvement at the expected level. The second concept of the model explains how to systematically develop a strategy [50]. After realizing the problem, the decision maker must develop an improvement strategy to achieve the goal. Traditionally, the multivariate analysis model assumes that the factors are independent of each other, so the decision makers will be guided to look at each problem individually rather than looking at the relationship between all the indicators, so they can only address the obvious problems-improve, not solve the root cause of the problem. In the real world, the cause of the problem is not simple and often involves complex systems. In the system, each factor has more or less the influence of the target, and the analyst often deliberately ignores the inconspicuous relationship in order to simplify the analysis process. Therefore, this research model recognizes that all factors are causally related to each other and then finds the root cause of the problem and discusses how to improve to achieve the goal.
In that vein, this study used a hybrid MCDM study model to perform analysis tools, composed of three technologies. The first is the decision-making trial and evaluation laboratory (DEMATEL), which uses an expert questionnaire to assess the relationship between each factor and expresses the interaction relationship through a pairwise comparison matrix to find the total impact relationship matrix [51,52].
The second technique is the analytic network process (ANP), which is designed to loosen the independence assumptions of factors in the analytic hierarchy process (AHP) system to reflect real-world conditions [53,54]. Then, the results of DEMATEL are combined with ANP (also known as DEMATEL-based ANP), and the influential weights (IWs) are derived from the calculation program. IWs are mainly used to express the relative importance of factors between the factors. The third technique is VIKOR. This study used modified VIKOR, which combines IWs with VIKOR. Its purpose is to shorten the gap between factors and expected level and find out where the problem lies, in order to refer back to the INRM drawn by DEMATEL. The INRM is used to identify the root cause of the problem and develop a systematic improvement strategy to prioritize the improvement strategy [55,56].
The MCDM model has fully demonstrated its advantages in various research fields. Some scholars have optimized sustainable energy alternatives in the planning and development of renewable energy [57]. Some scholars assess and build sustainable housing affordability on various complex factors that affect the quality of family life [58]. Some scholars have provided continuous improvement strategies for the establishment of creative cities [59]. In the field of air pollution, some scholars have used MCDM to carry out an assessment of the impact of air pollution on the sustainable economic development of cities [60].
This study used a research method that optimized the analytical procedures of the above three hybrid analysis techniques [61]. Due to the largely reduced number of traditional questionnaire items which allow more time for each question to be answered by experts, the quality of the answers for questions could be ensured. Its operating procedures are divided into 11 steps. The first step to the fifth step can obtain the total influence matrix through the DEMATEL method and obtain the data required to draw the INRM. Steps 6 through 7 obtain the local weights of the facets and criteria through the DANP method. The global weights for all criteria can be obtained through step 8, which can be used to understand the message across the facets. Gain the gap between each factor through modified VIKOR through steps 9 through 11. The operating procedures are as follows.
Step 1: Establish the direct relation matrix D of each expert. The direct relation matrix is established through a pair-wise comparison method to evaluate the degree to which factor i affects factor j of each of the H experts, as shown in Formula (1). It can be expressed as D = d h ij n×n for h = 1, 2, . . . , H, and n is the number of criteria. The evaluation scale is divided into 0, 1, 2, 3 or 4, which represent the magnitude of the different influences, ranging from no effect (0) to extremely high influence (4) [62]. (1) Step 2: Calculate the average direct relation matrix A. The average direct relation matrix A can be obtained by integrating the matrix D of the H-expert with an arithmetic averaging method. a ij = 1 H ∑ H h=1 d h ij indicates the degree to which the criterion i affects the criterion j, as shown in Equation (2).
Step 3: Calculate the normalized relation matrix X.
The average direct relation matrix A is normalized to obtain the normalized relation matrix X, as shown in Equations (3) and (4).
Step 4: Calculate the total relation matrix T.
Through Formula (5), the multiple influences and indirect influence values of the criterion, that is, the total influence value, can be calculated. The degree of influence between the criteria in matrix X is directly or indirectly affected by itself and other criteria, T = t ij n×n , i, j = 1, 2, · · · , n; I is the identity matrix.
Step 5: Draw the facet and criteria of the INRM. Adding each column and each row of matrix T to obtain the value of o and p, o means the total degree of the criterion affecting other criteria, as shown in Formula (6); p means that the criterion is affected by other criteria. The total extent is shown in Equation (7).
(o + p) is called prominence, which indicates the extent to which the criteria affect and the sum of the impacts and represents the importance of this criterion in the system. (o − p) is called relation, which indicates the degree of net impact of the criterion and represents the relative intensity of influence of this criterion in the system. Using the coordinate plot, with (o + p) as the X-axis and (o − p) as the Y-axis, we can get the INRM [63]. The criteria have an influential relationship between the two. The criterion with stronger influence will have another criterion with weaker ray pointing. In the questionnaire input stage, two questionnaires, dimensions and criteria, are used respectively [61].
Step 6: Transpose and normalize the total relation matrix. The normalized total relation matrix T α can be obtained by normalizing the matrix T by Formula (8), and then the transposed total influence relation matrix W α can be obtained by Formula (9).
Step 7: Obtain local weights of dimensions w l d and criteria w l d_c .
W α is multiplied by Formula (10) until it converges, and z represents the power of the self. The local weights of dimension w l d and criteria w l d_c within the dimension are obtained; the local weights in the DANP model is called influential weights (IWs).
Step 8: Get global weights of all criteria W g c . Global weights of all criteria W g c are obtained by integrating the local weights of dimensions with criteria, as shown in Formula (11).
Step 9: Obtain the aspiration levels and worst value. This study defines the criteria for assessing satisfaction by ranging from the worst (0) to the best (4). That is, f aspired is 4, and f worst is 0.
Step 10: Normalize performance of k alternatives and calculate the gap. By normalizing the difference between the performance value f kj of the scheme k and the aspiration level, r kj can be obtained, as shown in Formula (12).
Step 11: Obtain the mean group utility S k for the gap. This study aimed to develop an improvement strategy to achieve the aspiration level, so the modified VIKOR was used instead of the traditional VIKOR [55,64]. Therefore, the total gap of this study is based on the mean group utility S k , as shown in Equation (13).

Background and Problem Descriptions
Taiwan's air pollution is severe, especially in the Central and Southern regions. The Environmental Protection Administration [11] compared the pollution situation in various cities and counties in Taiwan in Spring 2018 and Kaohsiung was the most polluted city. See Kaohsiung's location in Figure 1 below. In addition to natural pollution factors, such as its topography and climate, Kaohsiung has many highly polluting heavy industries, factories, and coal-fired power plants. The study used the highly industrialized and commercialized Kaohsiung as an empirical case because it is also a crucial agriculture area and has the largest metropolitan area in the southern Taiwan, with a population of approximately 2.77 million people. This study sought to evaluate air pollution control measures and offer suggestions to improve their effectiveness using a hybrid MCDM model combining DANP with VIKOR.

Mobile air pollution sources (D2) Diesel vehicles (C21)
Diesel vehicle without the installation of filters. Scooters (C22) Two-stroke scooters emitting high PM2.5 concentrations. Ship emissions (C23) Ships emission entering and leaving the port. Railway transportation (C24) Railway transportation equipment emissions.
The location of Kaohsiung is indicated in the following map in Figure 1.

Analysis of Results
In this study, the DEMATEL structure was established and 2 dimensions and 11 criteria of domestic air pollution sources for Kaohsiung were analyzed. Based on the experts' consensus, the total influence matrices and were obtained. The sums of influences given and received for constructing the INRM are shown in Table 2 and the INRM is shown in Figure 2. The direction of the arrow in the INRM represents the direction of influence.

Data Collection
This study reviewed existing research to develop an evaluation model for an empirical case and provided textural refinements and corrections for the pretest questionnaire for these dimensions and criteria undertaken by a committee compromising five experts, including three civil servants from the Environmental Protection Bureau of Kaohsiung City Government and two academic experts in relevant professional fields regarding air pollution sources. A heterogeneous group of at least five experts is favorable for an expert questionnaire survey to achieve group consensus [50]. Therefore, the crucial dimensions and criteria that were adopted for the final questionnaire were suggested by the nine including aforementioned experts. The results are shown in Table 1. Dust from construction sites. Cooking fumes (C 15 ) Oil fumes from restaurants without recycling devices. Customs (C 16 ) Joss paper burning in traditional ceremonies. River dust (C 17 ) Dust from dry river beds.

Mobile air pollution sources (D 2 ) Diesel vehicles (C 21 )
Diesel vehicle without the installation of filters.

Ship emissions (C 23 )
Ships emission entering and leaving the port.

Railway transportation (C 24 )
Railway transportation equipment emissions.
The location of Kaohsiung is indicated in the following map in Figure 1.

Analysis of Results
In this study, the DEMATEL structure was established and 2 dimensions and 11 criteria of domestic air pollution sources for Kaohsiung were analyzed. Based on the experts' consensus, the Appl. Sci. 2019, 9, 1414 8 of 13 total influence matrices T D and T c were obtained. The sums of influences given and received for constructing the INRM are shown in Table 2 and the INRM is shown in Figure 2. The direction of the arrow in the INRM represents the direction of influence.  The matrix is the transposed matrix of . Subsequently, the matrix was obtained. The IWs of the DANP could be obtained by limiting the power of until it reaches a steady state. The global weights and local weights of the dimensions and criteria were obtained based on the IWs from the DANP; subsequently, the DANP was combined with the VIKOR method to assess the effectiveness of existing air pollution control measures in Kaohsiung, as shown in Table 3. Table 3. Performance evaluation of the empirical case using modified Vlse Kriterijuska Optimizacija The matrix W α is the transposed matrix of T α . Subsequently, the matrix W was obtained. The IWs of the DANP could be obtained by limiting the power of W until it reaches a steady state. The global weights and local weights of the dimensions and criteria were obtained based on the IWs from the DANP; subsequently, the DANP was combined with the VIKOR method to assess the effectiveness of existing air pollution control measures in Kaohsiung, as shown in Table 3.

Discussion
According to the results, the air quality of Kaohsiung needs to be improved by identifying key pollutants, then making effective strategies. In summary, the gap value of mobile air pollution sources (D 2 ) is the largest, with Scooters (C 22 ) and Diesel vehicles (C 21 ) listed in the first and second places. There are many scooters and diesel vehicles in the city, leading to a high degree of historical reliance on such transportation tools. Encouraging electronic vehicles might be a necessary method to achieve the replacement effect. This requires the local government to invest more in subsidy for electronic vehicles for individual and industrial transportation uses.
Because there are coal-fired power plants supporting heavy industries' growth, which is a major strategy of Kaohsiung's economic development., the criteria gap value of coal-fired power plants (C 11 ) and factories (C 12 ) are the first and second largest. According to statistics [11], almost 50% of all monitored (highly pollutant) factories of Taiwan are in Kaohsiung. Moreover, the largest refinery, petrochemical industrial zone and steel plant of Taiwan are also in Kaohsiung. For a long time, policy makers often neglected the upgrading of energy-consuming industries and pollutant processing industries by sacrificing the ecological environment in exchange for economic development. However, scholars have found that energy-consuming industries can cause carbon emissions and pollute the environment, but they do not necessarily lead to economic growth [65,66]. Therefore, in the long run, policymakers should consider implementing more aggressive policies to support energy-consuming industries and power plants to improve their efficiency without affecting economic growth [67].
The improvement strategy of this study could be discussed in dimensions and criteria. In terms of dimensions, through the facet of INRM, it can be seen that mobile air pollution sources (D 2 ) are affected by stationary air pollution sources (D 1 ). Therefore, policymakers need to consider the impact of stationary air pollution sources on mobile air pollution sources. In terms of stationary air pollution sources, experts hope to improve coal-fired power plants (C 11 ), but it is affected by certain factors (C 12 ). In other words, when thinking about how to improve air quality, policy makers should not just prioritize resources to reduce the level of power plants, but they should help the factory to install anti-fouling equipment through counselling mechanisms, requiring plants to gradually improve the power saving efficiency of the machine, reducing the burden on the power plant, increasing the monitoring equipment, and strengthening the elimination of excessive factory emissions.

Conclusions and Limitation
A hybrid modified MCDM model was utilized to evaluate the effectiveness of air pollution control and determine feasible improvement strategies in Kaohsiung City. The main contributions of this study are as follows. The main dimensions and criteria of the evaluation model for domestic air pollution sources were constructed according to the literature review and expert panel consensus. DEMATEL was used to verify the relationships between dimensions and criteria and to construct INRMs. DANP provided the IWs of criteria and modified VIKOR was used to acquire the gaps among dimensions and criteria. We used hybrid research methods to holistically obtain key indicators to target for continuous improvement of air pollution improvement. This empirical case study verified the effectiveness of our research structure and method model. Workplace-generated air pollution is critical for workers' and stakeholders' health [68]. Thus, based on the study result, we recommend that the Environmental Protection Bureau prioritize coal-fired power plants as a direction for improvement, which corresponds to the finding to previous studies well [69].
There are several limitations to this study. First, the study is a case study applied specifically to Kaohsiung City, which means that results might be different due to heterogeneity in different regions and times. Second, the respondents are environmental officials of the Kaohsiung City Government and experts who long live in Kaohsiung City. There may be some bias from their answers. Additionally, we only used Kaohsiung as the sample (which might limit the generalizability of our conclusion), though it is an important and representative port city in the world. Thus, we believe that the case is still worth noting and can be implication stimulating, even every city has its own special conditions. Future studies could follow conceptual and methodological logic to create context-specific investigations of different characterized cities. Also note that there are other possible emission reduction opportunities from coal-fired power plants (not by replacement with other energy sources, but by, for example, modernizing the flue gas treatment system or by changing the coal fuel to more environmentally friendly fuel) and significant impact of mobile air pollution sources on air quality (also resulting from the analysis). Future studies could find it useful to follow such a thought and design a proper research to investigate these possible strategies.