# Environmental Efficiency Evaluation of Chinese Industry Systems by Using Non-Cooperative Two-Stage DEA Model

Authors: Xiao Shi
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

In evaluating the environmental efficiency analysis of Chinese industry systems, data envelopment analysis (DEA) has been a popular method. However, the production system is often treated as a black box in conventional DEA models. This study considers the internal structure of the production system to evaluate the environmental efficiency, which is characterized as a two-stage system, i.e., production subsystem and pollutant treatment subsystem. And, in reality, some subsystems in two-stage production systems are not equally important, and this kind of two-stage systems usually has the feature that one subsystem dominates the other. Thus, we consider the leader and follower relationship in the environmental efficiency analysis. A new non-cooperative two-stage DEA model considering undesirable intermediates and undesirable outputs is proposed to calculate the environmental efficiency. The proposed method is then applied to 30 regional industry systems of China in the year 2010. Thus, each DMU’s environmental efficiencies for the overall system as well as both subsystems could be analyzed by the proposed approach. More accurate information could be provided for environmental management.


---

## Body

## 1. Introduction

As the policy of “opening and reforms” in 1978, China has achieved tremendous success in economic development. Its economic growth maintains high speed, but its environment deteriorates as many industrial pollutants are discharged to environment [1]. These environmental problems are mainly due to huge pollutants from industries. The government and people have realized that more efforts should be concentrated on protecting the environment. Thus, the Chinese government has paid much attention to the discharge of waste gas, waste water, and solid waste from industry [2, 3]. However, the environmental issues, such as resource shortage and environmental pollution, are still the problems China faces [4–6]. Environmental issues have attracted many attentions, and environmental efficiency analysis has become more critical.In previous literature, data envelopment analysis (DEA) has been widely applied to evaluate environmental performance. As a nonparametric approach, it has the advantage that it does not require any prior assumptions on the underlying functional relationships between variables of inputs and outputs. It is based on data-driven frontier analysis, and it floats a piecewise linear surface to rest on top of the observed quantities of inputs and outputs (see [7–11]). And, DEA has gained many new developments as a large amount of papers related to variation of models, and DEA applications have been published. DEA has been widely applied in environmental performance evaluation.Many scholars have published papers on the evaluation of environmental performance by DEA. The main issue when using DEA models to evaluate the environmental performance is to deal with undesirable outputs [12]. Many methods are proposed in DEA models to deal with undesirable outputs, for example, the data translation approach (e.g., [13–15]), treating undesirable outputs as inputs (e.g., [1, 16–18]), environmental DEA technology (e.g., [19–22]). However, the existing literature deems each DMU as a black box but neglects that some real production systems may have network structure as a two-stage network.In the real production process, the internal structure of DMU’s production system is often a two-stage production process ([23–29]; Liu et al. 2018). Thus, the two-stage structure of the production system should be considered when evaluating the environmental performance of this kind of production system. In previous literature, Song et al. [29] deemed the Chinese industry system as a two-stage production system but treated the two subsystems as two independent subsystems which ignored that the intermediate should be consistent in the two subsystems. Recently, Wu et al. [35] proposed a two-stage DEA model with shared inputs to evaluate the environmental performance of the Chinese industry, without considering whether there exists a cooperative relationship between two subsystems. Lozano [36] computes the environmental performance of a two-stage production and abatement system by a network slacks-based inefficiency measure. Liu et al. [32] further extended two-stage DEA models with undesirable indexes by considering undesirable input-intermediate-outputs.In reality, the two-stage production system usually has the feature that one subsystem dominates the other, for example, there are usually two participants in a two-stage supply chain, manufacturer and retailer. In such a case, the manufacturer usually holds manipulative power and acts as a leader, and the retailer acts as a follower in non-cooperative supply chains models [35]. Thus, when evaluating the efficiency of such systems, the leader and follower relationship should be considered (e.g., [35, 36]). In the existing literature, some scholars have proposed to model this kind of leader and follower relationship by introducing Stackelberg game theory in the efficiency evaluation. For example, Liang et al. [35] firstly considers the leader-follower relationship in the efficiency evaluation of two-stage supply chain, which deemed the manufacturer and retailer as leader and follower, respectively Zha and Liang [36] considered shared inputs among different stages, which further extended the non-cooperative two-stage approach. Li et al. [32] further extended Liang et al.’s [35] non-cooperative two-stage approach by considering additional inputs to stage two. Tavana and Khalili-Damghani [35] further considered the fuzzy data in the two-stage production process by proposing a new two-stage Stackelberg fuzzy DEA model.The aim of the paper is intended to calculate environmental efficiency for two-stage production systems based on Stackelberg game theory. We use a non-cooperative two-stage approach proposed by Liang et al. [35] to measure the environmental efficiencies of the Chinese regional industry system. Especially, the leader-follower relationship between the subsystem efficiencies will be investigated. And, by considering undesirable intermediates and undesirable outputs, we proposed a new non-cooperative two-stage DEA model to measure the environmental efficiency. Finally, the new proposed method is applied to measure the environmental efficiencies of Chinese regional industry systems.The structure of this paper is organized as follows. In Section2, we propose a new non-cooperative two-stage DEA model to measure the environmental performance of the two-stage production system. In Section 3, the new proposed approach is applied to Chinese regional industry systems in 2010. Conclusions are drawn in Section 4.


## 2. Environmental Efficiency Measurement of Chinese Industry Systems by Using Non-Cooperative Two-Stage DEA Model

Figure1 shows a two-stage production process. Suppose there exist n DMUs. Each DMU is denoted as DMUj1,2,…,n. In the production subsystem, inputs xij1,i=1,…,m are used to produce desirable outputs yrj1,r=1,…,s and undesirable outputs pdk1,d=1,…,D. Here, pdk1,d=1,…,D are denoted as intermediate measures. In the pollutant treatment subsystem, they are also used as the inputs. The exogenous inputs of the pollutant treatment subsystem are denoted as xhj2,h=1,…,H. In the pollutant treatment subsystem, the desirable outputs and undesirable outputs are ygj2,g=1,…,G and pbj2,b=1,…,B, respectively.Figure 1
The two-stage production system.If the DMU is deemed as a black box, and the internal structure is ignored, thenDMUj’s inputs and outputs are as follows: xij1,i=1,…,m and xhk2,h=1,…,H are inputs, yrj1,r=1,…,s and ygj2,g=1,…,G are desirable outputs, and pbj2,b=1,…,B are undesirable outputs. As for the undesirable outputs, pollutants are usually the by-products of the desirable outputs. Thus, in the measurement of environmental efficiency, it is expected as little as possible.To approach the undesirable outputs, a translation method proposed by Seiford and Zhu [36] is adopted to deal with the undesirable outputs. Firstly, each undesirable output is first multiplied by “−1.” Then, to make the undesirable output positive, we add an appropriate translation vector v2 to it. That is, pbj2¯=−pbj2+vb2,b=1,…,B, where vb2 could be obtained by vb2=maxjpbj2+1,b∈B. Thus, the larger the pbj2,b=1,…,B, the smaller the pbj2¯,b=1,…,B. As for evaluating the environmental efficiency of DMUk when deeming it as a black box, we can apply the following CCR model (see [37] to obtain it):(1)ekCCR=max∑r=1suryrk1+∑g=1Gπgygk2+∑b=1Bεbpbk2¯∑i=1mvixik1+∑h=1HQhxhk2,s.t.∑r=1suryrj1+∑g=1Gπgygj2+∑b=1Bεbpbj2¯∑i=1mvixij1+∑h=1HQhxhj2≤1,∀j,ur,vi,πg,εb,Qh≥0,∀i,r,g,b,h.In the production subsystem, the nonnegative weightsvi and ur are attached to its inputs and desirable outputs, respectively. And, in the pollutant treatment subsystem, the nonnegative weights Qh, πg, and εb are attached to inputs, desirable outputs, and undesirable outputs, respectively. ekCCR∗ is the optimal value of model (1), which is the environmental efficiency of DMUk when DMUk is deemed as a black box. Based on this approach, each DMU’s internal structure is ignored. Thus, we could not identify the environmental efficiencies of both subsystems.The study aims to measure the environmental efficiencies of DMUs with a production structure of two-stage network. In this paper, the non-cooperative two-stage DEA approach proposed by Liang et al. [35] is adopted in this study. When the first subsystem is assumed to be more important, it is deemed as a leader, and the second subsystem is deemed as a follower. For a specific DMUk, the environmental efficiency of its first subsystem (the leader) could be calculated by the following traditional CCR model [37]:(2)Ek1∗=Max∑d=1Dwdz¯dk+∑rP=1sPHrPyrPk∑iP=1mPviPxiPk,s.t.∑d=1Dwdz¯dj+∑rp=1sPHrpyrpj∑iP=1mPviPxiPj≤1,∀j,viP,wd,HrP≥0,∀iP,∀d,∀HrP,ur,vi,πg,εb,Qh≥0,∀i,r,g,b,h.In the production subsystem, the nonnegative weightsviP, wd, and HrP are attached to inputs, intermediates, and desirable outputs, respectively. The optimal environmental efficiency of this subsystem is Ek1∗, which is the optimal value of model (2). In model (2), the undesirable outputs are zdj,d=1,…,D. To deal with the undesirable outputs, we firstly multiply “−1” by each undesirable output and then add an appropriate translation vector p to the negative undesirable outputs. Then, the undesirable outputs are positive by this approach. That is, zdj¯=−zdj+pd,d=1,…,D, and pd is obtained by pd=maxjzdj+1,d∈D. In this approach, the larger the zdj,d=1,…,D, the smaller the zdj¯,d=1,…,D.Model (2) is a fractional program, but we can apply Charnes–Cooper (C-C) transformation to convert it into a linear model as follows:(3)Ek1∗=Max∑d=1Dwdz¯dk+∑rP=1sPHrPyrPk,s.t.∑d=1Dwdz¯dj+∑rp=1sPHrpyrpj−∑iP=1mPviPxiPj≤0,∀j,∑iP=1mPviPxiPk=1,viP,wd,HrP≥0,∀iP,∀d,∀HrP.The optimal environmental efficiency of the production subsystem is obtained by model (3). The intermediate measure connects the two subsystems. When the environmental efficiency of the production subsystem is obtained, the pollutant treatment subsystem will consider the variables viP, wd, and HrP which maintain Ek1=Ek1∗. Following the non-cooperative two-stage concept of Liang et al. [35] and Zha and Liang [36]; we propose the following model that maximizes the environmental efficiency of pollutant treatment subsystem in the objective function with the constraint that maintains Ek1=Ek1∗.(4)Ek2∗=Max∑b=1Bεbu¯bk+∑rT=1sTRrTyrTk∑d=1Dwdz¯dk+∑iT=1mTπiTxiTk,s.t.∑d=1Dwdz¯dj+∑rP=1PHrPyrPj∑iP=1mPviPxiPj≤1,∀j,∑b=1Bεbu¯bj+∑rT=1sTRrTyrTj∑d=1Dwdz¯dj+∑iT=1mTπiTxiTj≤1,∀j,∑d=1Dwdz¯dk+∑rP=1sPHrPyrPk∑iP=1mPviPxiPk=Ek1∗.In model (4), the nonnegative weights πiT, RrT, and εb are attached to inputs, desirable outputs, and undesirable outputs in the pollutant treatment subsystem, respectively. The optimal environmental efficiency of this subsystem is Ek2∗, which is the optimal value of model (4). In model (4), the undesirable outputs are ubj,b=1,…,B. To deal with the undesirable outputs, we firstly multiply “−1” by each undesirable output and then add an appropriate translation vector q to the negative undesirable outputs. Then, the undesirable outputs are positive by this approach. That is, ubk¯=−ubk+qb, (b = 1, …, B), and qb is obtained by qb=maxubk+1,b∈B. In this approach, the larger the ubk,b=1,…,B, the smaller the ubk¯. In model (4), the environmental efficiency of DMUk′ pollutant treatment subsystem is maximized with the constraint that the environmental efficiency of the production subsystem remains unchanged. Model (4) is a fractional program, which could be transformed to be a linear program as follows:(5)Ek2∗=Max∑b=1Bεbu¯bk+∑rT=1sTRrTyrTk,s.t.∑d=1Dwdz¯dk+∑iT=1mTπiTxiTk=1,∑d=1Dwdz¯dj+∑rP=1sPHrPyrPj−∑iP=1mPviPxiPj≤0,∀j,∑b=1Bεbu¯bj+∑rT=1sTRrTyrTj−∑d=1Dwdz¯dj−∑iT=1mTπiTxiTj≤0,∀j,∑d=1Dwdz¯dk+∑rP=1sPHrPyrPk=Ek1∗∑iP=1mPviPxiPk.The optimal value of model (5) is denoted as Ek2∗, which is the optimal environmental efficiency of the pollutant treatment subsystem. Then, for the overall system, DMUk′ environmental efficiency is Enon,1=Ek1∗∗Ek2∗.In a similar approach, if the pollutant treatment subsystem is assumed to be the leader,DMUk′ environmental efficiency of the pollutant treatment subsystem πk2∗ could be calculated by traditional CCR model with inputs (zdk and xipk) and desirable outputs (yrpk), and undesirable outputs (ubk). Then, we could calculate the environmental efficiency for the production subsystem πk1∗ by solving a DEA model that maintains πk2=πk2∗. Then, DMUk′ efficiency for the overall system is Enon,2=πk1∗∗πk2∗.


## 3. Application to Regional Industry Systems in China

In previous literature of evaluating the environmental efficiency, most studies treat the production process of Chinese regional industry systems as a black box and neglect its network structure. To overcome the problem, we consider the internal structure of Chinese regional industry systems by dividing the process into two subsystems: production subsystem and pollutant treatment subsystem. As the status of these two subsystems may not be equal, we consider the non-cooperative relationship between these two subsystems. In this section, the proposed non-cooperative two-stage DEA model will be applied to measure the environmental efficiency of Chinese regional industry systems. And, we provide some practical suggestions for governments to govern environmental protection in industrial sectors.In this section, the environmental efficiency of Chinese regional industry systems in the year 2010 will be analyzed by applying our non-cooperative two-stage DEA approach. We do not use the data after 2010 as the indicators of the statistical system and statistical technologies were revised by the Ministry of Environmental Protection after 2010 (some indicators are absent). As shown in Figure1, each regional industry system is divided into two subsystems: production subsystem and pollutant treatment subsystem.As for selecting the variables, we follow the existing work of Song et al. [29]. In the production subsystem, inputs- resources are utilized to generate desirable output-products and undesirable outputs-pollutants. The pollutants from the production subsystem are then disposed in the pollutant treatment subsystem by using pollutant investments. In the pollutant treatment subsystem, the desirable output is the value of comprehensive utilization of the three wastes (VCU), and the undesirable outputs are the emitted pollutants. In the production subsystem, we select the total amount of employees, fixed assets, and electricity as inputs. We select gross domestic product (GDP) as the desirable output as it can reflect the production subsystem intuitively. We select COD, SO2, and solid waste generated from the production subsystem as undesirable outputs as they are harmful. And, they are selected as re-inputs as they are disposed in the pollutant treatment subsystem. In the pollutant treatment subsystem, we select the value of comprehensive utilization of the three wastes (VCU) as desirable outputs, and COD, SO2, and solid waste generated as undesirable outputs. We obtain the data from China Statistical Yearbook 2010 published by the China National Bureau of Statistics in 2010. Some statistical data of Tibet are not available, thus Tibet is excluded in our study. The characteristics of the data are documented in Table 1.Table 1
Descriptive statistics of raw data.
Variables
Units
Max
Min
Mean
Standard deviation
Employees
10000 persons
2141.907
53.568
734.312
636.557
Fixed assets
1 billion yuan
12372.821
177.744
3810.847
2940.517
Electricity consumption
(100 million kW·h)
4060.130
159.020
1399.280
1012.942
COD generation
10000 tons
261.256
4.237
61.668
59.218
SO2 generation
10000 tons
453.535
11.598
172.277
106.494
Solid waste generation
10000 tons
31688.000
212.000
8031.000
6635.777
GDP of the second industry
1 billion yuan
21462.720
385.210
6442.604
5475.785
Pollutant investment
1 billion yuan
456759.300
4353.500
132325.607
105065.660
VCU
10000 yuan
49.268
0.488
14.488
10.878
COD emission
10000 tons
138.287
2.817
62.144
36.129
SO2 emission
10000 tons
12011.451
0.001
1925.202
2644.794
Solid waste emission
10000 tons
2863867.000
31623.000
592826.533
650121.434In Table2, we list the results of overall efficiency, production efficiency, and pollutant treatment efficiency based on our proposed approach. In order to analyze the difference that existed among the regions in terms of the overall efficiency scores when the production subsystem is a leader and the pollutant treatment subsystem is a leader, we present the cluster analysis map of the overall efficiency level of Chinese 30 regional industry systems in Figures 2 and 3, respectively.Table 2
Environmental efficiency scores of all regions.
Subsystem 1 as a leader
Subsystem 2 as a leader
Black-box
DMU
Region
E
non
,
1
E
k
1
∗
E
k
2
∗
E
non
,
2
π
k
1
∗
π
k
2
∗
CCR
Eastern area
1
Beijing
0.9949
1.0000
0.9949
0.9949
1.0000
0.9949
0.3566
2
Tianjin
0.9368
1.0000
0.9368
0.9368
1.0000
0.9368
0.6395
3
Hebei
0.2557
0.6047
0.4228
0.6047
0.6047
1.0000
0.5653
4
Liaoning
0.0787
0.7558
0.1041
0.5668
0.7558
0.7499
0.3431
5
Shanghai
0.9769
1.0000
0.9769
0.9769
1.0000
0.9769
0.3756
6
Jiangsu
0.3954
0.7981
0.4954
0.7981
0.7981
1.0000
0.6724
7
Zhejiang
0.8115
0.8115
1.0000
0.8115
0.8115
1.0000
1.0000
8
Fujian
0.5358
0.8517
0.6290
0.8479
0.8479
1.0000
0.3574
9
Shandong
0.1582
0.9132
0.1732
0.9132
0.9132
1.0000
0.6540
10
Guangdong
0.7690
1.0000
0.7690
0.7690
1.0000
0.7690
0.2271
11
Hainan
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
Average
0.6284
0.8850
0.6820
0.8382
0.8847
0.9480
0.5628
Central area
12
Shanxi
0.3314
0.5955
0.5564
0.5651
0.5810
0.9725
0.5649
13
Inner Mongolia
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
0.7123
14
Jilin
0.9602
1.0000
0.9602
0.9658
0.9989
0.9669
1.0000
15
Heilongjiang
0.7600
0.9522
0.7981
0.9301
0.9396
0.9899
0.6739
16
Anhui
0.4378
0.7967
0.5495
0.7344
0.7344
1.0000
0.6055
17
Jiangxi
0.5053
0.9413
0.5368
0.8959
0.8959
1.0000
0.9578
18
Henan
0.2086
0.7935
0.2629
0.7662
0.7905
0.9692
0.3441
19
Hubei
0.1943
0.8173
0.2378
0.7382
0.8157
0.9050
0.6699
20
Hunan
0.2438
0.8489
0.2872
0.7334
0.8380
0.8751
0.8293
21
Guangxi
0.5955
0.7289
0.8169
0.6168
0.6630
0.9304
0.5963
Average
0.5237
0.8474
0.6006
0.7946
0.8257
0.9609
0.6954
Western area
22
Chongqing
0.7845
0.9489
0.8267
0.8346
0.9389
0.8889
0.6388
23
Sichuan
0.2240
0.7736
0.2895
0.6235
0.7736
0.8059
0.3353
24
Guizhou
0.3837
0.4294
0.8935
0.4212
0.4212
1.0000
0.4755
25
Yunnan
0.4184
0.5153
0.8121
0.5054
0.5054
1.0000
1.0000
26
Shaanxi
0.3557
0.8670
0.4102
0.6584
0.8639
0.7622
0.4484
27
Gansu
0.2474
0.4123
0.6001
0.4462
0.4462
1.0000
0.6937
28
Qinghai
0.8805
0.9227
0.9542
0.9061
0.9061
1.0000
0.9728
29
Ningxia
0.6207
0.6872
0.9031
0.6381
0.6799
0.9384
0.9091
30
Xinjiang
0.6995
0.9217
0.7589
0.8244
0.9110
0.9050
1.0000
Average
0.5127
0.7198
0.7165
0.6509
0.7163
0.9223
0.7193
National average
0.5588
0.8229
0.6652
0.7675
0.8145
0.9446
0.6540Figure 2
Cluster analysis map of the overall efficiency level of 30 Chinese regional industry systems when subsystem 1 is a leader.Figure 3
Cluster analysis map of overall efficiency level of 30 Chinese regional industry systems when subsystem 2 is a leader.According to Figure2, when the production subsystem is a leader, the overall efficiency scores of 30 regional industry systems are significantly different. In the eastern area, 5 out of 11 regions (Beijing (P1), Tianjin (P2), Shanghai (P5), Zhejiang (P7), and Hainan (P11)) are highly efficient as their overall efficiency scores are above 0.8, 1 out of 11 regions (Guangdong (P10)) has overall efficiency scores between 0.6 and 0.8, 1 out of 11 regions (Fujian (P8)) has overall efficiency scores between 0.4 and 0.6, and 3 out of 11 regions (Hebei (P3), Liaoning (P4), Jiangsu (P6), and Shandong (P9)) are worst with overall efficiency scores below 0.3. In the central area, 2 out of 10 regions (Inner Mongolia (P13) and Jilin (P14)) obtain the highest overall efficiency scores, 1 out of 10 regions (Heilongjiang (P15)) has overall efficiency scores between 0.6 and 0.8, and 1 out of 10 regions (Guangxi (P21)) has an overall efficiency score between 0.4 and 0.6. Shanxi (P12), Anhui (P16), Jiangxi (P17), Henan (P18), Hubei (P19), and Hunan (P20) are the six worst DMUs with overall efficiency scores below 0.4. In the western area, Qinghai (P28) obtains the highest overall efficiency scores as it has an overall efficiency score above 0.8, 3 out of 9 regions (Chongqing (P22), Ningxia (P29), and Xinjiang (P30)) have overall efficiency scores between 0.6 and 0.8, 2 out of 9 regions (Yunnan (P25) and Shaanxi (P26)) have overall efficiency scores between 0.4 and 0.6, and 3 out of 9 regions (Sichuan (P23), Guizhou (P24), and Gansu (P27)) have overall efficiency scores below 0.4.According to Figure3, when the pollutant treatment subsystem is a leader, the overall efficiency scores of 30 regional industry systems are also significantly different. In the eastern area, 7 out of 11 regions (Beijing (P1), Tianjin (P2), Shanghai (P5), Zhejiang (P7), Fujian (P8), and Shandong (P9)) are highly efficient as their overall efficiency scores are above 0.8, and 3 out of 11 regions (Hebei (P3), Jiangsu (P6), and Guangdong (P10)) have overall efficiency scores between 0.6 and 0.8. In the central area, 6 out of 10 regions (Inner Mongolia (P13), Jilin (P14), Heilongjiang (P15), and Jiangxi (P17)) obtain the highest overall efficiency scores, 5 out of 10 regions (Anhui (P16), Henan (P18), Hubei (P19), Hunan (P20), and Guangxi (P21)) have overall efficiency scores between 0.6 and 0.8. Shanxi (P12) is the worst DMU with overall efficiency scores below 0.6. In the western area, 3 out of 9 regions (Chongqing (P22), Qinghai (P28), and Xinjiang (P30)) are highly efficient as their overall efficiency scores are above 0.8, 3 out of 9 regions (Sichuan (P23), Shaanxi (P26), and Ningxia (P29)) have overall efficiency scores between 0.6 and 0.8. Guizhou (P24), Yunnan (P25), and Gansu (P27) are the three worst DMUs with overall efficiency scores below 0.6.The results of the environmental efficiencies when deeming the DMUs as “black box” and the environmental efficiencies based on our new proposed approach are documented in Table2. When the production subsystem is assumed to be a leader, the corresponding maximum environmental efficiency of this subsystem Ek1∗, the minimum environmental efficiency of the pollutant treatment subsystem Ek2∗, and the environmental efficiency of the whole system Enon,1 are documented in column 3 to 5. When the pollutant treatment subsystem is assumed to be a leader, the corresponding maximum environmental efficiency of the production subsystem πk1∗, the minimum environmental efficiency of the pollutant treatment subsystem πk2∗, and the environmental efficiency of the whole system Enon,2 are documented in column 6 to 8. When the internal structure of regional industry system is neglected, the corresponding environmental efficiency ekCCR∗ is documented in column 9. It could be found that 5 regions are efficient by using Black-box DEA, while only 2 regions are efficient in the whole system when the production subsystem is considered as a leader in the two-stage production process (it is similar when the pollutant treatment subsystem is considered as a leader). It can be seen that the proposed non-cooperative two-stage DEA model identifies more sources of inefficiency than the traditional Black-Box DEA model. The reason is that the discriminating power of it is better than that of the traditional Black-Box DEA model. Based on our proposed approach, Hainan (P10) and Inner Mongolia (P13) are efficient in the overall system. It could be seen that any DMU is efficient in the overall system if and only if it is efficient in both subsystems.In order to compare the difference of the environmental efficiencies in three areas (i.e., eastern area, central area, and western area), we compare the efficiency results in Table2. It could be found that the average environmental efficiency for the overall system in the eastern area at 0.6248 is higher than that in the central area at 0.5237 and western area at 0.5127. As for the average environmental efficiency for the production subsystem and pollutant treatment subsystem, it could be found that there exists significant efficiency difference between these two subsystems.Based on the results, suggestions for improving the environmental efficiencies of 30 regional industry systems can be given. In the case that the production subsystem is assumed as a leader and in the case that the pollutant treatment subsystem is assumed as a leader, the environmental efficiency scores of both subsystems are used to plot the efficiency matrix as shown in Figures4 and 5, respectively. In these two figures, the horizontal axis of the efficiency matrix represents the environmental efficiency for the production subsystem and the vertical axis represents the environmental efficiency for the pollutant treatment subsystem. Each region is located in the matrix. When the production subsystem is a leader, the average environmental efficiency for the production subsystem of the 30 regions is 0.8229, and that of the pollutant treatment subsystem is 0.6652; when the pollutant treatment subsystem is a leader, the average environmental efficiency for the production subsystem of the 30 regions is 0.8145, and that of the pollutant treatment subsystem is 0.9446. Thus, the efficiency matrix could be divided into four quadrants using both average environmental efficiencies as shown in Figures 4 and 5. In the case of the production subsystem acting as a leader and that of the pollutant treatment subsystem acting as a leader, policy and strategic suggestions for improving environmental efficiency could be given based on the location of the regions in the four quadrants in Figures 4 and 5, respectively. For example, as seen in Figure 4, when the production subsystem is a leader, Jiangxi (17) is located in the fourth quadrant; it means Jiangxi (17) has large environmental efficiency in the production subsystem while it has small environmental efficiency in the pollutant treatment subsystem. Thus, if Jiangxi (17) wants to improve its environmental efficiency in the overall system, then it should put more emphasis on the pollutant treatment subsystem.Figure 4
Scatter plot of efficiency for two subsystems when the production subsystem is assumed to be a leader.Figure 5
Scatter plot of efficiency for two subsystems when the pollutant treatment subsystem is assumed to be a leader.In order to study the impact of pollutant intensities on the overall efficiency of a DMU (e.g., a region or a country), carbon intensity was used to monitor a DMU’s environmental performance in previous studies on environmental performance analysis [38]. The previous approach is adopted to define the intensities in both subsystems as follows: Electricity intensity (EI) is defined as the ratio of electricity consumption and GDP, COD-generation intensity (CGI) is defined as the ratio of COD generation and GDP, SO2-generation intensity (SGI) is defined as the ratio of SO2 generation and GDP, solid waste-generation intensity (SWGI) is defined as the ratio of solid waste generation and GDP, COD-discharge intensity (CDI) is defined as the ratio of COD discharge and VCU, SO2-discharge intensity (SDI) is defined as the ratio of SO2 discharge and VCU, and solid waste-discharge intensity (SWDI) is defined as the ratio of solid waste discharge and VCU. To describe the relationship between the overall efficiency score and its determinant factors, we formulate the following multiple regression equation:(6)OE=β0+β1CGI+β2SGI+β3SWGI+β4EI+β5CDI+β6SDI+β7SWDI+ε,where ε is a random disturbance term such that ε∼N0,oε2.Tables3 and 4 summarize the main regression results when the production subsystem is a leader and when the pollutant treatment subsystem is a leader, respectively. The coefficient of determination is denoted as R2 or R2adj. It shows that the multiple regression model fits the data well.Table 3
Determinants of the overall efficiency scores when subsystem 1 is a leader.
Independent variables
Coefficients
Standard error
Constant
0.211∗∗∗
0.156
CGI
0.148
0.259
SGI
−0.591∗∗
0.323
SWGI
−0.430
0.302
EI
0.474∗∗
0.257
CDI
−0.280
0.292
SDI
−0.835∗∗
0.307
SWDI
−0.192
0.232
Note:R2=0.667 and R2adj=0.603. Dependent variables are the overall efficiency scores of each province. ∗∗∗, ∗∗, and ∗ indicate the parameter estimate is significantly different from zero at 0.001, 0.01, and 0.05, respectively.Table 4
Determinants of the overall efficiency scores when subsystem 2 is a leader.
Independent variables
Coefficients
Standard error
Constant
0.137∗∗
0.109
CGI
0.259
0.181
SGI
−0.990∗∗∗
0.226
SWGI
−0.351
0.211
EI
0.294
0.180
CDI
−0.413
0.205
SDI
−0.633∗∗
0.215
SWDI
−0.111
0.162
Note:R2=0.853 and R2adj=0.641. Dependent variables are overall efficiency scores of each province. ∗∗∗, ∗∗, and ∗ indicate the parameter estimate is significantly different from zero at 0.001, 0.01, and 0.05, respectively.From Tables3 and 4, we can see that there are three statistically significant determinant factors when the production subsystem is assumed to be a leader, and there are two statistically significant determinant factors when the pollutant treatment subsystem is assumed to be a leader. When the production subsystem is a leader, SGI and SDI have significantly negative relationships with the overall efficiency, EI has significantly positive relationship with the overall efficiency, while CGI, SWGI, CDI, and SWDI have no significant relationship with the overall efficiency. Thus, for the government, to improve the overall environmental performance, more efforts should be made to reduce the emission of SO2 in the production subsystem and SO2 in the pollutant treatment subsystem and make full use of the electricity in the production subsystem. If the pollutant treatment subsystem is assumed to be a leader, two determinant factors (SGI and SDI) have significantly negative relationships with the overall efficiency, while CGI, SWGI, EI, CDI, and SWDI have no significant relationship with the overall efficiency. Thus, for the government, to improve the overall environmental performance, more efforts should be made to reduce the emission of SO2 in the production subsystem and SO2 in the pollutant treatment subsystem.


## 4. Conclusion

The Chinese regional industry system has structure of two-stage network. However, most published papers on environmental performance analysis ignored the internal structure of the regional industry system. In our study, the internal structure of the Chinese regional industry system is considered. In real production system, some subsystems in two-stage systems are not of equal importance, and this kind of two-stage systems usually has the feature that one subsystem dominates the other. Thus, in the efficiency evaluation of such systems, the leader and follower relationship should be considered. This is the first paper which introduced Stackelberg game theory to the environmental efficiency evaluation of two-stage production systems. Thus, we evaluate the environmental performance of the Chinese regional industry system not only considering the internal structure of production system but also considering the non-cooperative relationship between subsystems. These two subsystems are production subsystem and pollutant treatment subsystem, respectively. The newly proposed non-cooperative two-stage DEA approach is applied to the case of the Chinese regional industry system. Some interesting conclusions are drawn. To make our conclusions more accurate, lengthening the time span of the empirical study may be a direction of future studies.



---

*Source: 9208367-2019-10-09.xml*