Skip to main content

Dynamic linkages among economic development, environmental pollution and human health in Chinese

Abstract

Background

Research on the relationships between economic development, energy consumption, environmental pollution, and human health has tended to focus on the relationships between economic growth and air pollution, energy and air pollution, or the impact of air pollution on human health. However, there has been little past research focused on all the above associations.

Methods

The few studies that have examined the interconnections between the economy, energy consumption, environmental pollution and health have tended to employ regression analyses, DEA (Data Envelopment Analysis), or DEA efficiency analyses; however, as these are static analysis tools, the analyses did not fully reveal the sustainable economic, energy, environmental or health developments over time, did not consider the regional differences, and most often ignored community health factors. To go some way to filling this gap, this paper developed a modified two stage Undesirable Meta Dynamic Network model to jointly analyze energy consumption, economic growth, air pollution and health treatment data in 31 Chinese high-income and upper-middle income cities from 2013–2016, for which the overall efficiency, production efficiency, healthcare resource utilization efficiency and technology gap ratio (TGR) for all input and output variables were calculated.

Results

It was found that: (1) the annual average overall efficiency in China’s eastern region was the highest; (2) the production stage efficiencies were higher than the healthcare resource utilization stage efficiencies in most cities; (3) the high-income cities had lower TGRs than the upper–middle income cities; (4) the high-income cities had higher average energy consumption efficiencies than the upper-middle income cities; (5) the health expenditure efficiencies were the lowest of all inputs; (6) the high-income cities’ respiratory disease and mortality rate efficiencies were higher than in the upper–middle income cities, which had improving mortality rate efficiencies; and (7) there were significant regional differences in the annual average input and output indicator efficiencies.

Conclusions

First, the high-income cities had higher average efficiencies than the upper-middle income cities. Of the ten eastern region high-income cities, Guangzhou and Shanghai had average efficiencies of 1, with the least efficient being Shijiazhuang. In the other regions, the upper-middle income cities required greater technology and health treatment investments. Second, Guangzhou, Lhasa, Nanning, and Shanghai had production efficiencies of 1, and Guangzhou, Lhasa, Nanning, Shanghai and Fuzhou had healthcare resource utilization efficiencies of 1. As the average production stage efficiencies in most cities were higher than the healthcare resource utilization stage efficiencies, greater efforts are needed to improve the healthcare resource utilization. Third, the technology gap ratios (TGRs) in the high-income cities were slightly higher than in the upper-middle income cities. Therefore, the upper-middle income cities need to learn from the high-income cities to improve their general health treatment TGRs. Fourth, while the high-income cities had higher energy consumption efficiencies than the upper-middle income cities, these were decreasing in most cities. There were few respiratory disease efficiency differences between the high-income and upper-middle income cities, the high-income cities had falling mortality rate efficiencies, and the upper-middle income cities had increasing mortality rate efficiencies. Overall, therefore, most cities needed to strengthen their health governance to balance economic growth and urban expansion. Fifth, the average AQI efficiencies in both the high-income and upper-middle income cities were higher than the average CO2 efficiencies. However, the high-income cities had lower average CO2 emissions and AQI efficiencies than the upper-middle income cities, with the AQI efficiency differences between the two city groups expanding. As most cities were focusing more on air pollution controls than carbon dioxide emissions, greater efforts were needed in coordinating the air pollution and carbon dioxide emissions treatments. Therefore, the following suggestions are given. (1) The government should reform the hospital and medical systems. (2) Local governments need to strengthen their air pollution and disease education. (3) High-income cities need to improve their healthcare governance to reduce the incidence of respiratory diseases and the associated mortality. (4) Healthcare governance efficiency needs to be prioritized in 17 upper-middle income cities, such as Hangzhou, Changchun, Harbin, Chengdu, Guiyang, Kunming and Xi’an, by establishing sound medical management systems and emergency environmental pollution treatments, and by increasing capital asset medical investments. (5) Upper-middle income cities need to adapt their treatment controls to local conditions and design medium to long-term development strategies. (6) Upper-middle income cities need to actively learn from the technological and governance experiences in the more efficient higher-income cities.

Background

Global fossil fuel resources have been rapidly depleted since the beginning of the industrial revolution, and in the past few decades, climate change effects have become more apparent around the world. As a result, there has been an increased focus on energy efficiency and new clean energy technologies to reduce global carbon emissions. Because China is now the largest carbon dioxide emitter at around 28% of all global emissions, the Chinese government launched the “National Ambient Air Quality Standard”, which outlined aims to reduce its PM2.5 national emissions to 35 μg/m3 by 2030 [50]. The 13th Five-Year Plan, which was adopted on 15 March 2016, also included a goal to reduce sulfur dioxide and carbon dioxide emissions by 7% by 2040 [50]. Air pollutants and especially particulate matter have been found to contribute to increases in urban lung and cardiovascular diseases [32], with the World Health Organization (2018) reporting that in 2016, there were 4.2 million premature deaths due to PM2.5; 5 8% due to heart disease and stroke, 18% due to chronic obstructive pulmonary disease and acute lower respiratory infection, and 6% due to lung cancer. Because China has the world’s worst air quality, the Chinese government has begun to heavily invest in air quality improvements.

There has been significant research into the relationships between economic development, energy, and environmental pollution [17, 27, 48, 49, 52]), with some specifically focusing on energy and air pollution factors ([1, 8, 9, 12, 13, 18, 19, 26, 37, 48, 49, 54, 55]). The relationships between air pollution, health, and especially children’s health have also been major research foci [14, 16, 20, 21, 25, 31, 33, 36, 40, 46, 47, 51, 53, 56] and [6, 7, 15, 22, 23, 24, 30, 34, 38, 39]).

However, there has been little research on the interconnections between economic development, energy consumption, environmental pollution and health, or environmental pollution’s impact on social activities due to the health effects. The few studies that have examined these interconnections have tended to employ regression analyses, DEA (Data Envelopment Analysis), or static DEA efficiency tools that have lacked any dynamic considerations, which meant that the results did not fully reveal the sustainable economic, energy, environmental, and health developments over time. Therefore, past research has ignored the impact of social activities on human health, failed to reflect annual changes, and also failed to provide a comprehensive discussion on the connections between the economy, environmental pollution and health issues. Further, past research has only tended to examine the impact of environmental pollution such as CO2 and PM2.5 on first-stage economic production and have rarely considered the important AQI environmental indicators or discussed the relationships between production and health management.

Even though it is well known that different cities and different regions have different development levels, a majority of past studies have failed to provide such regional comparisons. While Feng et al. [14] recently employed a two-stage meta-frontier dynamic network data envelopment analysis (TMDN-DEA) model to explore energy consumption’s environmental pollution effects on child and adult mortality in 28 EU countries and 53 non-EU countries from 2010 to 2014, the research was mainly focused on the EU and non-EU countries. Therefore, to go some way to filling this gap, this study proposes a modified Undesirable Meta Dynamic Network model to explore the economic, energy, environmental, and human health efficiencies in 31 Chinese cities.

This research has two main contributions. First, this study explored the economic, energy and environmental pollution efficiencies across China and the associated government health expenditure and disease efficiencies. Second, to avoid any shortcomings associated with static analyses, carry-over impacts across the periods were included and the regional differences accounted for, for which a modified Undesirable Meta Dynamic Network model was developed to assess data from 2013–2016 in 31 Chinese cities in two stages; a production stage and a health treatment stage. In the production stage, labor and energy consumption were the inputs, GDP was the output, and CO2 (carbon dioxide) and the AQI (air quality index) emissions were the link variables between the production stage and health treatment stage, and in the second health treatment stage, health expenditure was the input, birth rate, respiratory disease rate, and mortality rate were the outputs, with the carryover being fixed assets.

The remainder of this paper is organized as follows. “Literature review” section gives the Literature review,  “Research method” section details the research method, “Empirical study” section examines the empirical results and provides an analytical discussion, and “Discussion” section gives the policy recommendations and managerial implications.

Literature review

Research on the connections between economic development, energy, environmental pollution, and human health has tended to follow three main paths: the relationships between economic growth and air pollution; energy and air pollution efficiencies; and the impact of air pollution on human health.

Relationships between economic growth and air pollution

Because the unfettered global focus on economic growth has led to a significant increase in environmental pollution, there has been a growing body of research focused on the relationships between economic growth and air pollution.

Georgiev and Mihaylo [17] tested the Environmental Kuznets Curve (EKC) hypothesis on four local and two global air pollutants and found that the EKC inverted U-shaped relationship between income and pollution did not hold for all gases, Xie et al. [52] found that an improvement in PM2.5 emissions would result in a growth in China’s GDP, Li et al. [27] conducted an environmental and economic analysis using a willingness to pay model, Wang et al. [48, 49] used panel smooth transition (PSTR) models to study the relationships between China’s economic growth and carbon dioxide emissions.

Energy and air pollution efficiency

Energy consumption and economic development are inevitably linked; however, higher energy consumption generally results in a commensurate increase in air pollution; therefore, there has been a significant research focus on energy efficiencies. For example, Hu and Wang [19] analyzed the energy efficiencies in 29 Chinese administrative regions from 1995 to 2002 and found that Central China had the worst energy efficiency, Fang et al. [12] used DEA to study energy performances in China and the United States, and found that China’s technical efficiency was worse, and Choi et al. [8] explored low carbon dioxide efficiency in China using an SBM model. In other research, Liou and Hu [26] calculated the ecological total-factor energy efficiencies (ETFEE) in 30 Chinese regions from 2005 to 2009 using an SBM model, and found that there was a monotonic increasing relationship between regional ETFEE and China’s per capita GDP. In more recent studies, Zhang and Choi [55] employed an SBM model and found that most provinces had low energy efficiencies, Apergis et al. [1] found that the energy efficiencies in capital-intensive OECD countries were higher than in energy-intensive OECD countries, and Yao et al. [54] employed panel data and a meta-frontier non-radial Malmquist CO2 emissions performance index (MNMCPI) and found that the Chinese industrial sector had an average annual CO2 emissions growth rate of 5.53%, the average industrial sector carbon dioxide emissions in the eastern, central and western regions had declined, and the MCPI has overestimated the carbon dioxide emissions efficiencies.

More recently, Wang et al. [48, 49] examined China’s energy productivity from 1995 to 2012, and found that the capital and energy trend substitutions were similar to the labor and energy trend substitutions and that energy productivity changes were mainly affected by technological progress, Guo et al. [18] used a dynamic SBM DEA to study the energy efficiency in 26 OECD countries, finding that Canada and China had the best, Qin et al. [37] explored the energy efficiency of Chinese coastal areas and found that the energy efficiencies were lower when undesirable output was considered, and Feng et al. [13] explored the effect of industrial structural adjustment market-oriented reforms and the strengthening of environmental protection measures on China’s CO2 emissions efficiency. In other research, Emrouznejad and Yang [9] reviewed energy efficiency DEA research, Li and Lin [28] used bootstrap to survey total factor energy consumption in China, and Li et al. [29] found that high-income Chinese cities had higher technological efficiencies than low-income Chinese cities from 2013 to 2016.

Impact of air pollution on human health

There has been substantial research into the effects of specific air pollutants because of the known adverse effects on community health. For example, Pope [36] found that a local steel mill was emitting 82% of all industrial \(PM_{10}\)(particulate matter) and that in some months from April 1985 to February 1988, the daily \(PM_{10}\) levels exceeded 150 μg/m in the Utah Valley. Using a Poisson regression, Maheswaran et al. [31] examined the effects of NO2 (nitrogen dioxide) PM10 and CO (carbon monoxide) on stroke mortality and hospital admissions in Sheffield, United Kingdom from 1994 to 1998, and found that outdoor air pollution levels were associated with higher stroke mortality and hospital admission risks. Fischer et al. [16] found that long-term exposures to PM10 and NO2 were related to mortality in people over the age of 30 in the Netherlands, Lelieveld et al. [25] found that PM2.5 (particulate matter) caused 3.3% of annual worldwide premature deaths, and Wu et al. [51] concluded that exposure to particulate air pollution affected circulating antioxidant enzymes. In more recent studies, Yang et al. [53] found that all pollutants were positively correlated with prehypertension in northeast China in 2009, Vlaandern et al. [47] found that short-term exposure to air pollution interfered with blood metabolism, Shen et al. [40] found that the current AQI system was unable to accurately estimate the air pollution health risks, and Zhao et al. [56] found that of the 307 cyclists interviewed in Beijing in 2015 on heavily polluted days, most were low-income males over 30 years old or short-distance travelers. In related research, Ngo et al. [33] found that sandstorms were associated with per capita acute respiratory diseases, Torres et al. [46] found that sulfur dioxide and fine particles increased cardiovascular and respiratory diseases, Huang et al. [20] found that family income and education were negatively correlated with ambient air quality in Beijing in 2014, and Lua et al. [30] estimated that in 2017, 3800 out of 124,800 deaths in Hainan were due to PM2.5 exposure.

The impact of air pollution on children’s health has also been widely examined. For example, Feng et al. [14] employed a two-stage meta-frontier dynamic network DEA (TMDN-DEA) model to explore the influence of energy consumption and the associated environmental pollution on the mortality of children and adults in 28 EU countries and 53 non-EU countries from 2010 to 2014, Lee et al. [24] used a generalized additive model (GAM) time series analysis to explore the effect of multiple air pollutants (sulfur dioxide (SO2), nitrogen dioxide (NO2), ozone (O3), carbon monoxide (CO) and aerodynamic diameter (PM10)) on the health of the children under 15 years old in Seoul from 1997 to 1999, and found that nitrogen dioxide and ozone were the main contributors to childhood asthma. In more recent research, Chen et al. [6] found that PM2.5 and PM10 affected the health and lung function of primary school children, Fioravanti et al. [15] found that traffic air pollution had no effect on childhood obesity in children aged 4 and 8 years in Rome, and Knibbs et al. [22] found that exposure to outdoor NO2 affected the respiratory systems of children aged 7–11 in Australia. In other similar studies, Salavati et al. [39] found that air pollutants affected pregnant women’s health, Chen et al. [7] found that a higher air pollution exposure was associated with respiratory disease and impaired lung function prevalence in younger children in China, Nobles et al. [34] found that air pollution caused fetal growth restrictions, Landrigan et al. [23] found that air, water, soil and chemical pollution caused 940,000 child deaths worldwide, and Roberts et al. [38] found no associations between pollution exposure, air quality and mental health issues.

Literature review summary

Most of the above research only focused on one or two aspects, with no studies having analyzed the relationships between economic development, energy, environmental pollution and human health. Therefore, to fill this gap, the method proposed in Feng et al. [14] was used to develop a modified Undesirable Meta Dynamic Network DEA model to examine the economic development, energy, environmental pollution and human health efficiencies in 31 Chinese cities in China, identify the regional differences, and provide suggestions for improvement.

Research method

Methodological framework

Past research on energy and environmental efficiencies ([19]) used labor, fixed assets and energy consumption as the inputs and GDP, CO2 and SO2 emissions as the outputs. This study used panel data from 31 of the most developed high-income and upper-middle income cities in China, the division for which was based on the World Bank’s classification for rich and poor countries, that is, the high-income economies had a GNI per capita of $12,056 or more, and the upper-middle income economies had a GNI per capita of between $3896 and $12,055.

From this calculation, the 31 sample cities were divided into 14 high-income cities: Beijing, Changsha, Fuzhou, Guangzhou, Hangzhou, Huhehot, Jinan, Nanchang, Nanjing, Shanghai, Shenyang, Tianjin, Wuhan, and Zhengzhou: and 17 upper-middle income cities: Chengdu, Changchun, Chongqing, Guiyang, Harbin, Haikou, Hefei, Kunming, Lanzhou, Lhasa, Nanning, Shijiazhuang, Taiyuan, Urumqi, Xian, Xining, Yinchuan. Data from 2013 to 2016 were extracted from the Statistical Yearbooks of China, the Demographics and Employment Statistical Yearbooks of China, and the Statistical yearbooks from each city. Air pollutant data were collected from China Environmental Protection Bureau reports.

In the first stage analysis of the energy and economic efficiencies in each city, labor and energy consumption were the input indicators and GDP was the output indicator, with the air quality index (AQI) and carbon dioxide (CO2) emissions being the link indicators, and in the second stage analysis of the government health expenditure efficiencies in each city, government health expenditure was the input indicator, and births, respiratory disease prevalence, and mortality rates were the output indicators, with the carryover being fixed assets.

The variables are explained in the following.


Input variables

Labor: Employees; the number of employees in each city at the end of each year; Unit = people.

Fixed assets: capital stock in each city calculated as the fixed assets investment in each city; Unit = 100 million CNY;

Energy consumption: total energy consumption in each city; Unit = 100 million tonnes.


Output variables

GDP (Gross Domestic Product): The first stage of the research model in this study used gross domestic product (GDP) as the expected output, which is an important measure of the economic production in a country, a region or a city. Unit = 100 million CNY.

Link Production Stage and health stage variables:

CO2 emissions: the CO2 emissions in each city were estimated from the energy consumption breakdown by fuel category.

Air quality Index (AQI): the measured pollutant concentrations for particulate matter: PM2.5, PM10, sulfur dioxide (SO2), Nitrogen Dioxide (NO2), Ozone(O3)and Carbon monoxide (CO), with the PM2.5 and PM10 being 24-h average concentrations.

Second stage: health treatment stage

Input variables: health expenditure

Output variables: Birth rate; respiratory diseases; mortality rate.

Figure 1 shows the Network Dynamic Model framework from Feng et al. [14] that was employed to explore the inter-period energy consumption, economic growth, air pollution and health efficiencies in 31 Chinese cities.

Fig. 1
figure 1

Network model

Model parameterization

When developing the traditional DEA model, based on Farrell [10] ’s theory for a generalized mathematical linear programming model that could measure multiple inputs and outputs at constant returns to scale, Charnes et al. [5] developed the CCR model, which was extended in 1984 by Banker et al. [2] to the BCC variable returns to scale (VRS) model. However, both the CCR and BCC models measured radial efficiency, which assumed that the input items and/or output items proportionally increased or decreased. However, as this assumption is not applicable in all situations, Tone [41] proposed a Slacks-Based Measure (SBM) model in 2001 that used the slack variables as the basis for measurement, considered the slacks between the input and output items, and used a non-radial estimation method to find the single value (scalar) improvement space. Then, Färe et al. [11] proposed Network Data Envelopment Analysis to determine the optimal solution under the CCR and BCC models, which found that the production process was composed of many sub-production technologies. Tone and Tsutsui [42] then proposed a weighted slacks-based measure in which the links between the departments of the various decision-making units were used as the analysis basis for the network DEA model, with each department being regarded as a sub-DMU, after which the SBM model was used to find the most suitable solution. Unlike traditional DEA models, as these sub-production technologies were seen as “black boxes”, the Network DEA was applied to explore the impacts of the input allocations and intermediate goods on the production process. In a later network DEA model development, Tone and Tsutsui [43] adopted a dynamic method to evaluate the DMU efficiencies at different times, and then introduced a carryover to connect the various DMU stages in the different periods.

As the Dynamic DEA model measures operational efficiencies across multiple periods and the Network DEA analyzes the efficiency of individual departments, the limitations of traditional DEA models were overcome. Therefore, in 2014, Tone and Tsutsui proposed a weighted slacks-based Dynamic Network DEA, which treated each department as a sub-DMU, and acknowledged the links between the various decision-making unit departments and the carry-over activities between the different periods.

As this study considered both undesirable outputs and regional differences, Tone and Tsutsui’s [44] dynamic network SBM (Slacks-Based Measure) model was modified to an Undesirable Meta-Frontier Dynamic Network Model. Battese and Rao [3] and Battese et al. [4] proposed a meta-frontier model that was able to compare the technical efficiencies of different groups, after which O’Donnell et al. [35] established a meta-frontier model that accurately calculated meta-frontier and group efficiencies. However, most DEAs assume that all DMUs have the same technology levels, but as the DMUs in this study were in different geographical locations and were subject to varying regional policies, they had differing technology levels. Therefore, in reference to Feng et al.’s [14] framework and based on Tone and Tsutsui’s [44] Network Dynamic SBM, Tone and Tsutsui [45] and O’Donnell et al.’s [35] meta-frontier model, and by including the undesirable outputs, this paper developed an Undesirable Meta-Frontier Dynamic Network Model on Chinese data to establish economic, environmental and health models.


Modified Undesirable Meta Dynamic Network model

Meta-frontier (MF)

It is assumed that all units (N) are composed of DMUs in g groups (N = N1 + N2 +…. + NG), where yrj and xij indicate the output item r (r = 1, 2, …, s) for item j (j = 1, 2, …, N) and input item i (i = 1, 2, …, m) for item j (j = 1, 2, …, N) under the meta-frontier. The DMU efficiency meta-frontier k is solved using the following linear programming (LP):

Suppose there are \(n {\text{DMUs }}\left( {j = 1, \ldots ,{\text{n}}} \right)\), with each having \(k\) divisions \(\left( {k = 1, \ldots ,K} \right)\), and \(T\) time periods \(\left( {{\text{t}} = 1, \ldots ,{\text{T}}} \right)\). Each DMU has an input and output at time period \(t\) and a carryover (link) to the next \(t + 1\) time period.

Inputs and outputs

\(X_{ijk}^{t} \in R_{ + } \left( {i = 1, \ldots ,m_{k} ; j = 1, \ldots ,n; k = 1, \ldots ,K; t = 1, \ldots ,T} \right)\) refers to input \(i\) at time period \(t\) for \(DMU_{j}\) division \(k\); \(X_{ijk}^{t}\). In the first stage, the number of employees, fixed assets and energy consumption are used as the input variables, and in the second health expenditure treatment stage, government health expenditure is used as the input variable.

\(Y_{rjk}^{t} \in R_{ + } \left( {r = 1, \ldots ,r_{k} ; j = 1, \ldots ,n; k = 1, \ldots ,K; t = 1, \ldots ,T} \right)\) refers to output r in time period \(t\) for \(DMU_{j}\) division \(k\); \(Y_{rjk}^{t}\). In the first stage, the GDP is the output, and in the second stage, the birth rate is the desirable output and respiratory disease and mortality rates are the undesirable outputs.

Links

\(Z_{{j\left( {kh} \right)t}}^{t} \in R_{ + } (j = 1, \ldots ,n; l = 1, \ldots , L_{hk} ; t = 1, \ldots ,T\) are the period \(t\) links from \(DMU_{j}\) division \(k\) to division \(h\), with \(L_{hk}\) being the number of \(k\) to \(h\) links; \(Z_{{j\left( {kh} \right)t}}^{t}\). CO2 and the AQI are the link indicators in both the first production stage and the second health treatment stage.

Carryovers

\(Z_{jkl}^{{\left( {t,t + 1} \right)}} \in R_{ + } \left( {j = 1, \ldots ,n; l = 1, \ldots , L_{k} ; k = 1, \ldots K; t = 1, \ldots ,T - 1} \right)\) are the carry-overs from \(t\) to the \(t + 1\) from \(DMU_{j}\) division \(k\) to division \(h\), with \(L_{k}\) being the number of carry-over items in division \(k\); \(Z_{jkl}^{{\left( {t,t + 1} \right)}}\). Fixed assets investment is the carry-over indicator in the production stage.

\({\text{Linkin}}_{k}\) is the number of input links for each division k, \({\text{Linkout}}_{k}\) is the number of output links for each division k, \({\text{ngood}}_{k}\) indicates the number of desirable carry-overs for each division \(k\), and \({\text{nbad}}_{k}\) indicates the number of undesirable carry-overs for each division \(k\).

The meta-frontier \(k\) for the DMU efficiency is solved using the following linear programming (LP):

Objective function

Overall efficiency:

$$\theta_{0}^{*} = \hbox{min} \frac{{\sum\nolimits_{{{\text{t}} = 1}}^{T} {W^{t} } \left[ {\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 - \frac{1}{{m_{k} + linkin_{k} + ninput_{k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iokg}^{t} }}} } + \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{(kl)_{l} = 1}}^{{linkin_{l} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} ing}}^{t} }}} } + \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{k_{l} }}^{{ngood_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} input}}^{(t,t + 1)} }}} } } \right)} \right]} } \right]}}{{\sum\limits_{t = 1}^{T} {W^{t} } \left[ {\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 + \frac{1}{{r_{1k} + r_{2k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokggood}^{t} }} + } } \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{r = 1}^{{r_{2k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokgbad}^{t} }}} } } \right)} \right]} } \right]}}$$
(1)

Subject to:

Production stage

\(x_{o1}^{t} = X_{1}^{t} \lambda_{1}^{t} + s_{1o}^{t - } \left( {\forall t} \right)\);

\(y_{o1good}^{t} = Y_{1good}^{t} \lambda_{1}^{t} - s_{1ogood}^{t + } , \left( {\forall t} \right)\);

\(\lambda_{1}^{t} \ge 0, s_{1o}^{t - } \ge 0, s_{1ogood}^{t + } \ge 0, \left( {\forall t} \right)\);

\(Z_{{o\left( {12} \right)in}}^{t} = Z_{{\left( {12} \right)in}}^{t} \lambda_{k}^{t} + S_{{o\left( {12} \right)in}}^{t} \left( {\left( {12} \right)in} \right)\)

Health treatment stage

\(x_{o2}^{t} = X_{2}^{t} \lambda_{2}^{t} + s_{2o}^{t - } \left( {\forall t} \right)\);

\(y_{o2good}^{t} = Y_{2good}^{t} \lambda_{2}^{t} - s_{2ogood}^{t + } \left( {\forall t} \right)\);

\(y_{o2bad}^{t} = Y_{2bad}^{t} \lambda_{2}^{t} + s_{2obad}^{t - } \left( {\forall t} \right)\);

\(\lambda_{2}^{t} \ge 0, s_{2o}^{t - } \ge 0, s_{2ogood}^{t + } \ge 0, s_{2obad}^{t - } \ge 0 \left( {\forall t} \right)\);

\(e\lambda_{k}^{t} = 1\left( {\forall k, \forall t} \right)\);

$$\begin{aligned} \sum\nolimits_{j = 1}^{n} {Z_{{jk_{1} \alpha }}^{{(t,\left( {t + 1} \right)}} \lambda_{jk}^{t} } & & = & & \sum\nolimits_{j = 1}^{n} {Z_{{jk_{1} \alpha }}^{{(t,\left( {t + 1} \right)}} \lambda_{jk}^{t + 1} \left( {\forall k, \forall t} \right)} \\ \sum\nolimits_{j = 1}^{n} {Z_{{jk_{1} \alpha }}^{{(t,\left( {t + 1} \right)}} \lambda_{jk}^{t} } & & & = & \mathop {\sum\nolimits_{j = 1}^{n} {Z_{{jk_{1} \alpha }}^{{(t,\left( {t + 1} \right)}} \lambda_{jk}^{t + 1} + S_{{ok_{l} input}}^{{\left( {t,t + 1} \right)}} (k_{l} = 1, \ldots ,ngood_{k} ;\;\forall k, \forall t);} }\limits_{{}}^{{}} \\ S_{{ok_{l} good}}^{{\left( {t,t + 1} \right)}} & & \ge 0,(\forall k, \forall t) \\ \end{aligned}$$
(2)

(b) Period and division efficiencies

The period and division efficiencies are as follows:

(b1) Period efficiency

$$\partial_{0}^{*} = \hbox{min} \frac{{\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 - \frac{1}{{m_{k} + linkin_{k} + ninput_{k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iokg}^{t} }} + } } \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{(kl)_{l} = 1}}^{{linkin_{l} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} ing}}^{t} }}} } + \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{k_{l} }}^{{ninput_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} ginput}}^{(t,t + 1)} }}} } } \right)} \right]} }}{{\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 + \frac{1}{{r_{1k} + r_{2k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokggood}^{t} }} + } } \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokgbad}^{t} }}} } } \right)} \right]} }}$$
(3)

(b2) Division efficiency

$$\varphi_{0}^{*} = \hbox{min} \frac{{\sum\limits_{{{\text{t}} = 1}}^{T} {W^{t} } \left[ {1 - \frac{1}{{m_{k} + linkin_{k} + ninput_{k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iokg}^{t} }} + } } \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{(kl)_{l} = 1}}^{{linkin_{l} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} ing}}^{t} }} + } } \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{k_{l} }}^{{ninput_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} ginput}}^{(t,t + 1)} }}} } } \right)} \right]}}{{\sum\limits_{t = 1}^{T} {W^{t} } \left[ {1 + \frac{1}{{r_{1k} + r_{2k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokggood}^{t} }} + \sum\nolimits_{{{\text{g}} = 1}}^{\text{G}} {\sum\nolimits_{r = 1}^{{r_{2k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokgbad}^{t} }}} } } } } \right)} \right]}}$$
(4)

(b3) Division period efficiency

$$\rho_{0}^{*} = \hbox{min} \frac{{\left[ {1 - \frac{1}{{m_{k} + linkin_{k} + ninput_{k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iokg}^{t} }} + } } \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{(kl)_{l} = 1}}^{{linkin_{l} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} ing}}^{t} }} + \sum\nolimits_{g = 1}^{G} {\sum\nolimits_{{k_{l} }}^{{ninput_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} ginput}}^{(t,t + 1)} }}} } } } } \right)} \right]}}{{\left[ {1 + \frac{1}{{r_{1k} + r_{2k} }}\left( {\sum\nolimits_{g = 1}^{G} {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokggood}^{t} }} + \sum\nolimits_{{{\text{g}} = 1}}^{\text{G}} {\sum\nolimits_{r = 1}^{{r_{2k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokgbad}^{t} }}} } } } } \right)} \right]}}$$
(5)

From the above, using the meta-frontier model, the overall efficiency, period efficiency, division efficiency, and division period efficiency can be obtained.

Group-frontier (GF)

As each DMU under the group frontier chooses the most favorable final weighted output, the DMU efficiencies under the group frontier are solved using the following equations.

(a) Objective function

Overall efficiency

$$\theta_{0}^{*} = \hbox{min} \frac{{\sum\nolimits_{t = 1}^{T} {W^{t} \left[ {\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 - \frac{1}{{m_{k} + linkin_{k} + ninput_{k} }}\left( {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iok}^{t} }}} + \sum\nolimits_{{(kh)_{l} = 1}}^{{linkin_{k} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} in}}^{t} }} + \sum\limits_{{k_{l} }}^{{ninput_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} input}}^{(t,t + 1)} }}} } } \right)} \right]} } \right]} }}{{\sum\nolimits_{t = 1}^{T} {W^{t} \left[ {\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 + \frac{1}{{r_{1k} + r_{2k} }}\left( {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokgood}^{t} }} + \, \sum\nolimits_{{{\text{r}} = 1}}^{{r_{2k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokbad}^{t} }}} } } \right)} \right]} } \right]} }}$$
(6)

(b) Period and division efficiencies

The period and division efficiencies are as follows

(b1) Period efficiency

$$\partial_{0}^{*} = \hbox{min} \frac{{\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 - \frac{1}{{m_{k} + linkin_{k} + ninput}}\left( {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iok}^{t} }}} + \sum\nolimits_{{(kh)_{l} = 1}}^{{linkin_{k} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} in}}^{t} }} + \sum\limits_{{k_{l} }}^{{ninput_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} input}}^{(t,t + 1)} }}} } } \right)} \right]} }}{{\sum\nolimits_{k = 1}^{K} {W^{k} \left[ {1 + \frac{1}{{r_{1k} + r_{2kk} }}\left( {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokgood}^{t} }} + \, \sum\nolimits_{{{\text{r}} = 1}}^{{r_{2k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokbad}^{t} }}} } } \right)} \right]} }}$$
(7)

(b2) Division efficiency

$$\varphi_{0}^{*} = \hbox{min} \frac{{\sum\nolimits_{t = 1}^{T} {W^{t} \left[ {1 - \frac{1}{{m_{k} + linkin_{k} + ninput_{k} }}\left( {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iok}^{t} }}} + \sum\nolimits_{{(kh)_{l} = 1}}^{{linkin_{k} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} in}}^{t} }} + \sum\nolimits_{{k_{l} }}^{{ninput_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} input}}^{(t,t + 1)} }}} } } \right)} \right]} }}{{\sum\nolimits_{t = 1}^{T} {W^{t} \, \left[ {1 + \frac{1}{{r_{1k} + r_{2k} }}\left( {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokgood}^{t} }} + \, \sum\nolimits_{{{\text{r}} = 1}}^{{r_{2k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokbad}^{t} }}} } } \right)} \right]} }}$$
(8)

(b3) Division period efficiency

$$\rho_{0}^{*} = \hbox{min} \frac{{1 - \frac{1}{{m_{k} + linkin_{k} + ninput_{k} }}\left( {\sum\nolimits_{i = 1}^{{m_{k} }} {\frac{{S_{iok}^{t - } }}{{x_{iok}^{t} }}} + \sum\nolimits_{{(kh)_{l} = 1}}^{{linkin_{k} }} {\frac{{s_{{o(kh)_{l} in}}^{t} }}{{z_{{o(kh)_{l} in}}^{t} }} + \sum\nolimits_{{k_{l} }}^{{ninput_{k} }} {\frac{{s_{{ok_{l} input}}^{(t,t + 1)} }}{{z_{{ok_{l} input}}^{(t,t + 1)} }}} } } \right)}}{{\left( {1 + \frac{1}{{r_{1k} + r_{2k} }}\left( {\sum\nolimits_{r = 1}^{{r_{1k} }} {\frac{{s_{rokgood}^{t + } }}{{y_{rokgood}^{t} }} + \, \sum\nolimits_{r = 1}^{{r_{2k} }} {\frac{{s_{rokbad}^{t - } }}{{y_{rokbad}^{t} }}} } } \right)} \right)}}$$
(9)

From the above results, the overall efficiency, the period efficiency, the division efficiency and division period efficiency are obtained.

Technology gap ratio (TGR)

As the meta-frontier model contains g groups, the technical efficiency of the meta-frontier (MFE) is less than the technical efficiency of the group frontier (GFE); therefore, the ratio value or technology gap ratio (TGR) is:

$${\text{TGR}} = \frac{{\rho^{ *} }}{{\rho_{o}^{ *g} }} = \frac{MFE}{GFE}$$
(10)

Energy consumption, Health expenditure, Respiratory Diseases, Mortality Rate, CO2 and AQI Efficiencies

Hu and Wang’s [19] total-factor energy efficiency index was used to overcome any possible biases in the traditional energy efficiency indicators, for which there were eight key efficiency models; energy, environment, health expenditure, respiratory diseases, mortality rate, CO2 emissions, and the AQI. In this study, “I” represents area and “t” represents time.

The efficiency models are defined in the following:

$${\text{Inputefficiency}} = \frac{\text{Target input}}{\text{Actual input}}$$
(11)
$${\text{Undesirable output efficiency }} = \frac{\text{Target Undesirable output}}{\text{Actual Undesirable output}}$$
(12)

If the target inputs equal the actual inputs, then the efficiencies are 1, which indicates overall efficiency; however, if the target inputs are less than the actual inputs, then the efficiencies are less than 1, which indicates overall inefficiency.

If the target undesirable outputs are equal to the actual undesirable outputs, then the efficiencies are 1, indicating overall efficiency; however, if the target undesirable outputs are less than the actual undesirable outputs, then the efficiencies are less than 1, indicating overall inefficiency.

Empirical study

Statistical analysis of the inputs and outputs

Table 1 shows the average input–output indicator quantities in each Chinese region, from which it can be seen that there are large regional differences.

Table 1 Input and output variables from 2013–2016 by region

The labor force indicators in Table 1 show that the urban labor average in the eastern region was the highest and was rising in all periods to 2016. The urban average inputs in the northeastern region were also rising with some fluctuations in all period to 2016. However, the labor inputs were relatively small in the eastern and western regions, with the average labor input in the western region being the lowest of the four regions in 2013 and 2014 but higher than the central region in the subsequent 2 years and showing a straight upward trend to the its highest level in 2016. Although the labor input in the central region also increased slightly, the input in 2015 and 2016 was the lowest of the four regions.

In 2013, the urban fixed assets average in the northeast was the highest of the four regions; however, after that time, the average fell from 2014 to 2016. The largest fixed assets increases were in the eastern region, which from 2014 to 2016 was the highest of all four regions. While the fixed assets input in the western region also continued to rise, the average urban fixed assets input from 2013 to 2015 was the smallest of the four regions, but in 2016 was higher than in the northeast. In general, the differences between the regions were expanding.

The eastern region had the largest energy consumption inputs in all time periods, the western region had the smallest in 2013, but from 2014 to 2016, the northeast region had the smallest and declining energy consumption inputs, with i s lowest level in all 4 years being in 2016.

The highest urban GDP average in all 4 years was in the eastern region, and even though the urban GDP average was rising in the western region, it was the lowest of the four regions in all 4 years.

Even though the carbon dioxide emissions in the eastern region were showing a linear downward trend in all 4 years, in 2016, they were still much higher than in the other regions. The carbon dioxide emissions were the lowest in the western region in 2014, but the lowest in the northeastern region in 2013, 2015 and 2016. The carbon dioxide emissions in the northeastern region were also rising slightly and the gap between the regions was narrowing.

The highest air pollution index value was in the central region in 2013, 2014 and 2016; however, in 2015, the average air pollution index in the northeastern region was higher than in the other regions. The lowest average urban air pollution index was in the western region from 2013 to 2015, but in 2016, it was the lowest in the northeastern region. The air pollution index in all four regions had a downward trend, indicating that the air quality in most regions was improving, with only the western region rebounding slightly in 2016.

Results and analysis

Overall efficiency analysis

This study first compared the overall efficiencies in each city from 2013 to 2016, after which the first production stage efficiencies and the second healthcare resource utilization stage efficiencies were compared.

The overall efficiency was the output divided by the input, which was then expressed as a percentage; therefore, a perfect process would have an efficiency of 100%.

Appendix 1 shows the overall efficiency in each city from 2013 to 2016, from which it can be seen that Guangzhou, Lhasa, and Shanghai had overall efficiencies of 1 for all 4 years, Beijing had an overall efficiency of 1 in the first year, which dropped to below 0.6 in 2015 and 2016, Nanning had an overall efficiency of 1 for the first 3 years, which in the last year dropped to around 0.68, and Jinan’s efficiency was poor in the first 3 years, but reached 1 in the final year. Fuzhou, Haikou, and Urumqi had better overall efficiencies than the other cities; however, Urumqi’s efficiency was lower in 2015 at 0.6 but higher than 0.8 in the other years. These results were similar to the empirical results in Choi et al. [8] and Zhang and Choi [55]. In general, the overall efficiencies in most provinces needed improvement.

In most other cities, the overall efficiencies over the 4 years were below 0.6. Shijiazhuang had the lowest, with its highest efficiency in the 4 years being less than 0.3, the best annual efficiency in Taiyuan was less than 0.4 in 2013, Guiyang, Chengdu, Lanzhou, Kunming, Xian, and Xining’s had efficiencies of around 0.4, and Changchun, Harbin, Hangzhou, Hefei, Huhehot, Nanchang, Nanjing, Wuhan and Zhengzhou had efficiencies ranging from 0.4 to 0.6; therefore, most of these cities needed efficiency improvements.

Only six cities: Changsha, Chongqing, Hefei, Jinan, Shenyang, and Taiyuan: had increasing efficiencies, five of which were middle-income cities. Jinan’s efficiency, however, rose from around 0.4 to 1, and Hefei’s rose from below 0.5 in 2013 to close to 0.9 in 2016. Over the period, 23 cities saw decreasing efficiencies, with Beijing, Nanning and Wuhan experiencing the largest declines. Wuhan’s overall efficiency in 2013 was around 0.9, which dropped to only 0.4 in 2016, and Beijing’s overall efficiency in 2013 was 1, which by 2016 had fallen to around 0.5.

In reference to the empirical framework in Feng et al. [14], Table 2 shows the annual efficiencies in China’s eastern, central, northeastern and western regions. The annual average overall efficiency in the eastern region was the highest, but at 0.665, there was a great need for improvements. The average overall efficiency in the western region was significantly higher than in the central and northeastern regions in 2013 and 2015, but fluctuated and declined and by 2016 was lower than the central region for the first time. Of the four regions, the lowest average total efficiency was in the northeastern region at only 0.5, and had a fluctuating downward trend, with the lowest being in 2016 at only 0.381; therefore, significant improvements were needed. The annual average overall efficiency in the central region also had a fluctuating downward trend; however, there was a significant rise in 2016 to be second highest behind the eastern region, which was in line with the empirical results in Yao et al. [54]. Overall, all regions required overall efficiency improvements.

Table 2 Average overall efficiencies from 2013–2016 by region

Production efficiency and healthcare resource utilization efficiency analyses

The significant regional economic growth and social development differences meant that there were large production and health treatment stage efficiencies differences between the cities (see Appendix 2). Guangzhou, Lhasa, Nanning, and Shanghai had production stage efficiencies of 1, and Guangzhou, Lhasa, Nanning, Shanghai and Fuzhou had heath treatment stage efficiencies of 1. Generally, the production stage efficiencies in most cities were higher than the healthcare resource utilization stage efficiencies. A detailed analysis of the two-stage efficiencies is given in the following.

  1. i.

    Production efficiency

Production efficiency is when an economic system is unable to produce any more of one good without sacrificing the production of another good or without improving the production technology. In other words, production efficiency is when a good or a service is produced at the lowest possible cost. In simple terms, production efficiency is illustrated on a production possibility frontier on which all points on the curve are indicators of productive efficiency. However, equilibrium may be productively efficient without being allocative efficient, that is, it may result in a goods distribution in which the social welfare is not maximized.

A city/region’s production efficiency is when all services and enterprises within a city operate using best-practice technological and managerial processes and there are no further reallocations that can result in greater output with the same inputs and production technology. By improving these processes, however, an economy or a business can extend its production possibility frontier outward so that the efficient production yields a greater output than previously.

Twenty-one of the 31 cities had higher production stage efficiencies than healthcare resource utilization stage efficiencies, which indicated that these cities needed to place a greater focus on health treatment governance.

While Nanchang, Nanjing, Jinan and Beijing had production stage efficiencies over 0.8 for 2 years, 3 years and 4 years respectively, except for Lanzhou, Shijiazhuang, Taiyuan, and Xining, which had production stage efficiencies of less than 0.4, most other cities had production stage efficiencies between 0.4 and 0.8.

Hefei, Huhehot, Jinan, Nanjing, Taiyuan, and Xining had increasing production efficiencies, with the largest increase being in Jinan from below 0.6 in 2013 to 1 in 2016. The production stage efficiencies in the other 21 cities fell, with the biggest declines being in Fuzhou, Urumqi, and Zhengzhou to around 0.2. These results were similar to the empirical results in Wang et al. [48, 49], in which it was found that most regions had low production stage efficiencies and large regional differences.

  1. ii.

    Healthcare resource utilization efficiency

The average healthcare resource utilization efficiencies were far lower than the average production stage efficiencies. Fuzhou, Haikou, and Urumqi had higher healthcare resource utilization efficiency stage efficiencies than production stage efficiencies in all 4 years, and Lanzhou, Yinchuan and Xining had higher healthcare resource utilization efficiencies than production stage efficiencies in 3 years. Therefore, the local governments in these ten cities should focus more resources on improving their production efficiencies.

Besides the five cities that had healthcare resource utilization efficiencies of 1 in all 4 years, the efficiencies in the other cities were generally low. Fifteen upper-middle income cities: Changchun, Chengdu, Guiyang, Harbin, Hangzhou, Huhehot, Kunming, Nanchang, Nanjing, Shenyang, Shijiazhuang, Taiyuan, Tianjin, Xian, and Zhengzhou: had healthcare resource utilization efficiencies below 0.4, with eight of these: Changchun, Chengdu, Harbin, Nanchang, Shijiazhuang, Tianjin, Xian and Zhengzhou: having healthcare resource utilization efficiencies below 0.2.

Table 3 shows the average annual production and health treatment stage efficiencies in the four regions. The highest annual average production stage efficiency was in the eastern region, which rose marginally over the 4 years. The central region had an average production stage efficiency just below the eastern region in 2013 and 2014, but it began falling in 2015, when it was slightly lower than in the northeast. The average efficiency in the central region had a fluctuating downward trend, with the lowest efficiency being in 2016. Of the four regions, the western region had the lowest average production stage efficiency, which had a downward fluctuating trend to its lowest in 2016 of only 0.569. Although the production stage efficiency in the northeast region was higher than in the western region and in 2015 was higher than both the western and central regions, in 2016 it fell to only 0.563, which was lowest of the four regions.

Table 3 Average production efficiency and healthcare resource utilization efficiency from 2013–2016 by region

The regional health treatment stage annual average efficiencies varied widely across the four regions. The eastern region had the highest efficiency, followed by the western region, the central region and the northeastern region, which dropped to only 0.085 in 2015 and had a rebound in 2016 to about 0.2, its highest point in the 4 years. The health treatment stage annual average efficiency in the central region was lower than either the eastern or western region and had a U-shaped change, rising to its highest point in 4 years in 2016 at 0.374. However, in 2014 and 2015, the central region’s health treatment stage annual average efficiency was lower than 0.25, which was only slightly better than the northeast region; therefore, there were significant improvements needed.

Average overall efficiency and average TGR comparison between the city groups

Battese and Rao [3] stated that “the technology gap ratio indicates the technology gap for the given group according to currently available technology for firms in that group, relative to the technology available in the whole industry”.

To more clearly analyze the regional differences, the high-income and upper-middle-income city overall efficiencies and average TGRs were compared (Table 4).

Table 4 TGRs and efficiency from 2013–2016 for the high-income and upper-middle income cities

From Table 4, it can be seen that the high-income cities had higher average efficiencies than the middle-income cities, but there were some variances. The average efficiency in the high-income cities increased from 0.63 in 2013 to around 0.66 in 2014, declined in 2015 to around 0.57, and rebounded slightly in 2016 to 0.60. However, the upper-middle-income cities as a whole had lower average overall efficiencies. The average overall efficiency decreased slightly from 0.52 in 2013 to 0.51 in 2014, fell to 0.47 in 2015 and rose slightly to 0.49 in 2016; however, the overall efficiency decline was slightly less than in the high-income cities. The Wilcoxon Test indicated that the differences in the high-income and upper-middle income city efficiencies were significant at 5% from 2013 to 2016, which indicated that the efficiencies in the high-income cities were significantly better than in the upper-middle income cities.

The average TGRs were similar in the high-income and upper- middle income cities, with the average TGR in the high-income cities being only slightly higher at 0.8 in 2014 than the upper-middle income cities. However, in 2015 and 2016, the average TGR in the upper-middle income cities was slightly higher than in the high-income cities, indicating that the TGR was falling in the high-income cities and rising in the upper-middle income cities, and that the TGR gap between the two city types was narrowing. The Wilcoxon Test found that the TGR differences between the high-income and upper-middle income cities was insignificant at a 5% level from 2013 to 2016; that is, the average TGR in the upper-middle income cities was slightly higher than in the high-income cities from 2013 to 2016.

Average efficiency indicators in the high income and upper middle-income cities

Table 5 compares the average CO2, AQI, energy consumption, the average GDP, respiratory disease, mortality rate and health expenditure efficiencies in the high-income and upper-middle income cities.

Table 5 Average indicator efficiencies in the high-income and upper-middle income cities

Except for 2015, the average CO2 emissions and AQI were lower in the high-income cities than in the upper-middle income cities. The average CO2 emissions efficiency in the high-income cities declined from 0.81 in 2013 to 0.67 in 2016, and in the upper-middle income cities, first rose from 0.66 in 2013 to around 0.67 in 2015 and then fell to 0.52 in 2016; therefore, the average CO2 emissions efficiencies between the city groups narrowed slightly in 2016. In general, however, the high-income cities needed to strengthen their CO2 emissions management and consider coordinating their CO2 and air pollutant emissions management.

The AQI efficiency in the high-income cities increased from 0.84 in 2013 to 0.97 in 2014, fell to 0.71 in 2015, then rebounded to 0.89 in 2016, and in the upper-middle income cities increased from 0.80 in 2013 to 0.94 in 2014, declined in 2015 to 0.67, and rebounded slightly to 0.80 in 2016; therefore, the average AQI efficiency gap widened between the two city groups.

While the high-income cities had higher energy consumption efficiency than the upper-middle income cities, this fluctuated down from around 0.89 in 2013 to around 0.67 in 2016. The energy consumption efficiency in the upper-middle income cities also fluctuated down from 0.65 in 2013 to around 0.52 in 2016. Therefore, improvements were needed in both city groups. These results were similar to the empirical results in Li et al. [29] in which it was found that the high-income cities had higher efficiencies than the upper-middle income cities.

The average health expenditure efficiencies in both city groups were low, with the average in the high-income cities being 0.49 and the average in the upper- middle income cities being 0.461. In 2014, the average health expenditure efficiency in the high-income cities rose to around 0.53, but in the upper-middle income cities fell slightly to below 0.46. In 2015, the health expenditure efficiency in the high-income cities fell sharply to 0.36 and in the upper-middle income city average health expenditure efficiency fell to 0.37. In 2016, both city groups had significant increases to around 0.45, with the upper-middle income cities being slightly higher than the high-income cities.

The average mortality rate efficiency gap between the two city groups also narrowed. The average mortality rate efficiency rose in the high-income cities from 0.81 in 2013 to 0.84 in 2016, but dropped in the upper-middle income cities from around 0.88 in 2013 to 0.84 in 2016.

While the average healthcare expenditure efficiency in the high-income cities was slightly higher than in the upper-middle income cities, it was declining in both city groups, reaching its lowest in 2015. The average healthcare expenditure efficiencies in the high-income cities declined faster over the 4 years, and by 2015 and 2016, was lower than in the upper-middle income cities; therefore, significant improvements were needed across the board.

Average efficiency indicators by region

Table 6 shows the main annual average input and output indicator efficiencies in each region.

Table 6 Average efficiencies from 2013–2016 by region

The highest regional average carbon dioxide emissions efficiency in 2013 and 2014, was in the northeast region and the lowest was in the western region, and in 2015 and 2016, the highest carbon dioxide emissions efficiency was in the eastern region, and the lowest efficiencies were in the western and northeastern regions. From 2013 to 2015, the annual average carbon dioxide emissions efficiency in the western region was the lowest of all four regions, but in 2016, was only lower than the eastern region. Regardless, the average carbon dioxide emissions efficiency in the western region was fluctuating and declining, and there was a great need for improvement.

The regional average AQI index efficiency in the eastern region was higher than in the other four regions, and in 2013, 2014 and 2015, the AQI index efficiency was higher than the central and northeastern regions. Except for 2014, the average AQI index efficiency in the northeastern region was the lowest of the four regions and in 2015 was only 0.387, and although there was a significant increase to 0.78 in 2016, there was a great need for improvement. In the western region, the average AQI index efficiency in 2013 was the highest of the four regions. It increased significantly in 2014 but was still slightly lower than the eastern region, and then began to decline significantly. Although it rebounded in 2016, it was lower than in the eastern and central regions. The average AQI index efficiency in the central region was low in 2013 and 2015, was the lowest of the four regions in 2014, but was still high at 0.934, but in 2016 was higher than in the western and northeastern regions. Overall, the average AQI index efficiencies in the western and northeastern regions had a slight fluctuating upward trend.

The average annual respiratory diseases efficiency was the lowest in the northeastern region, and was only 0.524 in 2016. The average respiratory diseases efficiencies in the eastern, central, and western regions, however, were relatively close.

The northeast region had the lowest average mortality rate efficiency, which except for 2014, was only around 0.52 in the other 3 years. The average mortality rate efficiency in the eastern region was the highest of the four regions in 2015 and 2016 at above 0.9, the western region had the highest average mortality rate efficiency in 2013 and 2014 at above 0.85, and the average mortality rate efficiency in the central region ranked second in 3 years, but was ranked third in 2014.

The energy consumption efficiency was high in all regions. The energy consumption efficiency in the northeast region was the highest in 2013 and 2014, was slightly lower than the eastern region in 2015 and was lower than both the eastern and western regions in 2016. The average energy consumption efficiency in the western region from 2013 to 2015 was the worst at below 0.65. In the eastern region, the energy consumption efficiency was slightly lower than in the northeast in 2013 and 2014 but was the highest in 2015 and 2016. The energy consumption efficiency in the central region continued to decline and by 2016 was only 0.453; therefore, there was significant need for improvements.

The health expenditure efficiencies were low in all regions. The health expenditure efficiency in the eastern region was higher than in other regions but in 2016, was only 0.56. The health expenditure efficiency in the northeast region was lower than 0.3 in all 4 years and only 0.27 in 2016, and in the central region was only slightly higher. The health expenditure efficiency in the western region was slightly higher than that in the northeast and central regions, but in 2016, was only about 0.45.

Of the six above-mentioned indicators, the average health expenditure efficiencies in each region were the lowest, followed by the energy consumption and carbon dioxide emissions efficiencies. Therefore, future regional governance attention needs to be paid to health expenditure, energy consumption and carbon dioxide efficiency improvements.

The AQI index, respiratory disease, mortality, and energy consumption efficiencies in the northeastern region were the lowest of the four regions for most of the years; therefore, governance in these areas needs attention.

While the eastern region had higher AQI index, mortality rate and health expenditure efficiencies than the other regions, the health expenditure efficiency needed significant improvements.

The average energy consumption efficiency in the western region was low and while the health expenditure ranked second to the eastern region, governance improvements are still needed.

Discussion

Using a modified two stage Undesirable Meta Dynamic Network model, this study jointly analyzed energy consumption, economic growth, air pollution and health treatment expenditure data from 31 Chinese high-income and upper-middle income cities from 2013 to 2016, from which the following conclusions were made found:

(i) The high-income cities had higher average efficiencies than the upper-middle income cities. Of the 10 eastern region high-income cities, Guangzhou and Shanghai had average efficiencies of 1, with the least efficient being Shijiazhuang. In the other regions, the upper-middle income cities required greater technology and health treatment investments. These results indicated that the environmental and health expenditure efficiencies of the more economically developed regions were significantly higher than in the less developed regions. Economically developed areas are mainly located in the eastern coastal areas of China, most of which have relatively good weather conditions and the economic capacity to support healthy investment and development.

(ii) Guangzhou, Lhasa, Nanning, and Shanghai had production stage efficiencies of 1, and Guangzhou, Lhasa, Nanning, Shanghai and Fuzhou had healthcare resource utilization stage efficiencies of 1. The average production stage efficiencies in most cities were higher than the healthcare resource utilization stage efficiencies; therefore, greater efforts are needed to improve the healthcare resource utilization efficiencies. Guangzhou and Shanghai were the most economically developed cities in China in terms of economic strength and social development and achieved the best production and environmental efficiencies. Lhasa, which is is dominated by tourism, has a sparse population, and few high-polluting industries, had good environmental production stage and health efficiencies. Fuzhou, which is located in the coastal area of southeastern China, has good diffusion conditions and a high economic level, had the most efficient production, environmental and health input efficiencies.

(iii) The technology gap ratios (TGRs) in the high-income cities were slightly higher than in the upper-middle income cities, indicating that the high-income cities had slightly higher technological levels. Therefore, the upper-middle income cities need to learn from the high-income cities to improve their general TGR.

(iv) While the high-income cities had higher energy consumption efficiencies than the upper-middle income cities, this was decreasing in most cities. There were insignificant respiratory disease efficiency differences between the high-income and upper-middle income cities; however, the high-income cities had decreasing mortality rate efficiencies while the upper-middle income cities had increasing mortality rate efficiencies. Overall, most cities needed to strengthen their health governance to ensure balanced economic growth and urban expansion.

(v) The average AQI efficiencies in both the high-income and upper-middle income cities were higher than the average CO2 efficiencies. However, the high-income cities had lower average CO2 emissions and AQI index efficiencies than the upper-middle income cities, and the AQI efficiency differences between the two city groups was expanding. As most cities were focusing on air pollution controls rather than carbon dioxide emissions controls, greater efforts are needed to coordinate the air pollution and carbon dioxide emissions treatments.

In summary, climate change, economic growth and social development all pose great challenges to development. Therefore, to effectively respond to these challenges and problems, governments need to actively adapt measures to local conditions, develop scientific governance systems, and formulate short, medium- and long-term dynamic strategic management directions. Therefore, the following policy suggestions are given.

(1) Hospital and medical system reform is needed and a wide-ranging medical management system developed that has local medical systems based on regional characteristics and local disease characteristics to ensure a complete governance system from prevention to emergency response to supervision to security. For example, for areas that need to deal with the air pollutant discharges such as Zhengzhou, the focus should be on strengthening investment in the prevention and treatment of the diseases caused by air pollution and especially by PM2.5. Investment should also be strengthened in labor and the materials to ensure active prevention, active response and effective supervision and control.

(2) Local governments need to strengthen awareness in adolescents and residents about the health problems and diseases caused by air pollution. There is a large local education gap in China, with the high-income cities generally having better education than the upper-middle income cities; therefore, the education and policy promotions in the upper middle-income cities needs to be improved. Due to regional differences, different cities need to take different measures.

(3) High-income cities need to take action based on their specific problems and challenges. Overall, healthcare governance requires medium and long-term response measures and coordinated governance responses to reduce the incidence of respiratory diseases and the associated mortality.

(4) Of the 17 upper-middle income cities, Hangzhou, Changchun, Harbin, Chengdu, Guiyang, Kunming and Xi’an need to prioritize healthcare governance efficiency by establishing sound medical management systems and emergency environmental pollution treatments. For example, while Hangzhou is a coastal city that is conducive to the discharge of harmful air pollutants, industrial environmental pollution has adversely affected resident health. Therefore, the city needs to put more emphasis on healthcare treatment and air pollution treatment investments. Chengdu is in a basin surrounded by mountains, which is not conducive to the discharge of air pollutants; therefore, medical monitoring should be prioritized, and in addition to strengthening medical monitoring governance in Kunming, the treatment of carbon dioxide emissions should take precedence over the treatment of air pollutants.

(5) As high-income cities have more advanced technologies and fewer polluting industries, the CO2 and AQI emissions treatments have achieved some good results. However, upper-middle income cities such as Lanzhou, Taiyuan and Shijiazhuang, all of which have high energy consumption and high polluting industry structures, need to focus on reducing their CO2 and AQI emissions. Therefore, the local governments need to adapt the treatment controls to local conditions and design medium to long-term development strategies.

(6) Upper-middle income cities need to actively learn from the technology governance experiences in the more efficient higher-income cities. Based on their city’s economic development stage, economic resources and industry structures, local and regional governments need to develop health governance models that meet their own characteristics, and when reforming their current medical systems, high-income cities could learn from countries that have good medical systems such as Japan and the United Kingdom. For example, special attention should be given to developing and employing competent medical workers. Education programs for adolescents and residents about the health problems linked to air pollution could help them take environmental protection action

Conclusion

This study used DEA to determine the input and output improvements needed in 31 Chinese cities and to propose corresponding policy recommendations. Climate change, economic growth and social development all pose great challenges to human development. Therefore, to effectively respond to these challenges and problems, the government needs to actively adapt measures to local conditions, develop scientific governance systems, and formulate short, medium- and long-term dynamic strategic management directions. All regions should take active measures based on local meteorological conditions, geological characteristics, resource endowments and population characteristics to improve their energy consumption, economic, environmental and health efficiencies.

As the sample data were obtained from publicly released Chinese government data, the data used in this research were the latest available. Therefore, future research could continue to track the air pollution, economic growth and health governance efficiencies based on updated data and provide the central and local governments with more timely and effective policy recommendations and management decision-making references.

Availability of data and materials

Not applicable.

References

  1. Apergis N, Aye GC, Barros CP, Gupta R, Wanke P. Energy efficiency of selected OECD countries: a slacks based model with undesirable outputs. Energy Economics. 2015;51:45–53.

    Google Scholar 

  2. Banker D, Charnes AC, Cooper WW. Some models for estimating technical and scale inefficiencies in data envelopment analysis. Manage Sci. 1984;30:1078–92.

    Google Scholar 

  3. Battese GE, Rao DSP. Technology gap, efficiency, and a Stochastic Metafrontier Function. Int J Bus Econ. 2002;1:87–93.

    Google Scholar 

  4. Battese GE, Rao DSP, O’Donnell CJ. A meta frontier production function for estimation of technical efficiencies and technology gaps for firms operating under different technologies. J Prod Anal. 2004;21:91–103.

    Google Scholar 

  5. Charnes AC, Cooper WW, Rhodes E. Measuring the efficiency of decision making units. Eur J Oper Res. 1978;2:429–44.

    Google Scholar 

  6. Chen C, Li C, Li Y, Liua J, Meng C, Han J, Zhang Y, Xua D. Short-term effects of ambient air pollution exposure on lung function: a longitudinal study among healthy primary school children in China. Sci Total Environ. 2018;645:1014–20.

    CAS  PubMed  Google Scholar 

  7. Chen Z, Cui L, Cui X, Li X, Yu K, Yue K, Dai Z, Zhou J, Jia G, Zhang J. The association between high ambient air pollution exposure and respiratory health of young children: a cross sectional study in Jinan, China. Sci Total Environ. 2019;656:740–9.

    CAS  PubMed  Google Scholar 

  8. Choi Y, Zhang N, Zhou P. Efficiency and abatement costs of energy-related CO2 emissions in China: a slacks-based efficiency measure. Appl Energy. 2012;98:198–208.

    Google Scholar 

  9. Emrouznejad A, Yang GL. A survey and analysis of the first 40 years of scholarly literature in DEA: 1978–2016. Socio-Econ Plan Sci. 2018;61:4–8.

    Google Scholar 

  10. Farrell MJ. The measurement of productive efficiency. Journal of Royal Statistical Society. 1957;120(3):253–90.

    Google Scholar 

  11. Färe R, Grosskopf S, Whittaker G. Network DEA modeling data irregularities and structural complexities in data envelopment analysis. Boston: Springer; 2007. p. 209–40.

    Google Scholar 

  12. Fang H, Wu J, Zeng C. Comparative study on efficiency performance of listed coal mining companies in China and the US. Energy Policy. 2009;37(12):5140–8.

    Google Scholar 

  13. Feng C, Zhang H, Huang JB. The approach to realizing the potential of emissions reduction in China: an implication from data envelopment analysis. Renew Sustain Energy Rev. 2017;71:859–72.

    CAS  Google Scholar 

  14. Feng Y, Yu X, Chiu YH, Chang TH. Dynamic linkages among economic development, energy consumption, environment and health sustainable in EU and Non-EU Countries. Healthcare. 2019;7:138.

    PubMed Central  Google Scholar 

  15. Fioravanti S, Cesaroni G, Badaloni C, Michelozzi P, Forastiere F, Porta D. Traffic-related air pollution and childhood obesity in an Italian birth cohort. Environ Res. 2018;160:479–86.

    CAS  PubMed  Google Scholar 

  16. Fischer PH, Marra M, Ameling CB, Hoek G, Beelen R, de Hoogh K, Breugelmans O, Kruize H, Janssen NAH, Houthuijs D. Air pollution and mortality in seven million adults: the Dutch Environmental Longitudinal Study (DUELS). Environ Health Perspect. 2015;123(7):697.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Georgiev E, Mihaylo E. Economic growth and the environment: reassessing the environmental Kuznets Curve for air pollution emissions in OECD countries. Lett Spat Resour Sci. 2015;8(1):29–47.

    Google Scholar 

  18. Guo X, Zhu Q, Lv L, Chu J, Wu J. Efficiency evaluation of regional energy saving and emission reduction in China: a modified slacks-based measure approach. J Cleaner Prod. 2017;140:1313–21.

    Google Scholar 

  19. Hu JL, Wang SC. Total-factor energy efficiency of regions in China. Energy Policy. 2006;34:3206–17.

    Google Scholar 

  20. Huang G, Zhou W, Qian Y, Fisher B. Breathing the same air? Socioeconomic disparities in PM25 exposure and the potential benefits from air filtration. Sci Total Environ. 2019;657:619–26.

    CAS  PubMed  Google Scholar 

  21. Khaniabadi YO, Sicard P, Takdastan A, Hopke PK, Taiwo AM, Khaniabadi F, Omidi MAD, Daryanoosh M. Mortality and morbidity due to ambient air pollution in Iran. Clin Epidemiol Global Health. 2018;7(2):222–7.

    Google Scholar 

  22. Knibbs LD, Cortés de Waterman AM, Toelle BG, Guo Y, Denison L, Jalaludin B, Marks GB, Williams GM. The Australian Child Health and Air Pollution Study (ACHAPS): a national population-based cross-sectional study of long-term exposure to outdoor air pollution, asthma, and lung function. Environ Int. 2018;120:394–403.

    CAS  PubMed  Google Scholar 

  23. Landrigan PJ, Fuller R, Fisher S, Suk WA, Sly P, Chiles TC, Bose-O’Reilly S. Pollution and children’s health. Sci Total Environ. 2019;650(2):2389–94.

    CAS  PubMed  Google Scholar 

  24. Lee J, Kim H, Song H, Hong Y, Cho Y, Shin S. Air pollution and asthma among children in Seoul, South Korea. Epidemiology. 2002;13:481–4.

    PubMed  Google Scholar 

  25. Lelieveld J, Evans JS, Fnais M, Giannadaki D, Pozzer A. The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature. 2015;525:367–71.

    CAS  PubMed  Google Scholar 

  26. Liou LB, Hu JL. Ecological total-factor energy efficiency of regions in China. Energy Policy. 2012;46:216–24.

    Google Scholar 

  27. Li L, Lei Y, Email P, Chen DY, Si C. Economic evaluation of the air pollution effect on public health in China’s 74 cities. Springer Plus. 2016;5:402.

    PubMed  Google Scholar 

  28. Li K, Lin B. An application of a double bootstrap to investigate the effects of technological progress on total-factor energy consumption performance in China. Energy. 2017;128(1):575–85.

    Google Scholar 

  29. Li Y, Chiu YH, Lin TY. Energy and environmental efficiency in different Chinese regions. Sustainability. 2019;11:1216–42.

    CAS  Google Scholar 

  30. Lua X, Lin C, Li W, Chen Y, Huang Y, Funga JCH, Lau AKH. Analysis of the adverse health effects of PM2.5 from 2001 to 2017 in China and the role of urbanization in aggravating the health burden. Sci Total Environ. 2019;652(20):683–95.

    Google Scholar 

  31. Maheswaran R, Haining RP, Brindley P, Law J, Pearson T, Fryers PR. Outdoor air pollution and stroke in Sheffield, United Kingdom. A small-area level geographic study. Stroke. 2005;36:239–43.

    Google Scholar 

  32. Morabet REl. Effects of Outdoor Air Pollution on Human Health. Reference Module in Earth Systems and Environmental Sciences. 2018.

  33. Ngo NS, Bao X, Zhong N. Local pollutants go global: the impacts of intercontinental air pollution from China on air quality and morbidity in California. Environ Res. 2018;165:473–83.

    CAS  PubMed  Google Scholar 

  34. Nobles CJ, Grantz KL, Liu D, Williams A, Ouidir M, Seeni I, Sherman S, Mendola P. Ambient air pollution and fetal growth restriction: physician diagnosis of fetal growth restriction versus population-based. Sci Total Environ. 2019;650:2641–7.

    CAS  PubMed  Google Scholar 

  35. O’Donnell CJ, Rao DSP, Battese GE. Metafrontier frameworks for the study of firm-level efficiencies and technology ratios. Empiri Econ. 2008;34:231–55.

    Google Scholar 

  36. Pope CA. Respiratory diseases associated with community air pollution and a steel mill, Utah Valley. Am J Public Health. 1989;79(5):623–8.

    PubMed  PubMed Central  Google Scholar 

  37. Qin Q, Li X, Li L, Zhen W, Wei YMi. Air emissions perspective on energy efficiency: an empirical analysis of China’s coastal areas. Appl Energy. 2017;185:604–14.

    CAS  Google Scholar 

  38. Roberts S, Arseneault L, Barratt B, Beevers S, Danese A, Odgers CL, Moffitt TE, Reuben A, Kelly FJ, Fisher HL. Exploration of NO2 and PM2.5 air pollution and mental health problems using high-resolution data in London-based children from a UK longitudinal cohort study. Psychiatry Res. 2019;272:8–17.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Salavati N, Strak M, Burgerhof JGM, de Walle HEK, Erwich JJHM. Bakker MK Air pollution removal by urban forests in Canada and its effect on air quality and human health. Int J Hyg Environ Health. 2018;221:1061–7.

    CAS  PubMed  Google Scholar 

  40. Shen F, Ge X, Hu J, Nie D, Tian L, Chen M. Air pollution characteristics and health risks in Henan Province, China. Environ Res. 2017;156:625–34.

    CAS  PubMed  Google Scholar 

  41. Tone K. A slacks-based measure of efficiency in data envelopment analysis. Eur J Oper Res. 2001;130(3):498–509.

    Google Scholar 

  42. Tone K, Tsutsui M. Network DEA: a slacks based measurement approach. Eur J Oper Res. 2009;197:243–52.

    Google Scholar 

  43. Tone K, Tsutsui M. Dynamic DEA: a Slacks-based measure approach. Omega. 2010;38:145–56.

    Google Scholar 

  44. Tone K, Tsutsui M. Dynamic DEA with network structure: a slacks-based measure approach. Omega. 2014;42(1):124–31.

    Google Scholar 

  45. Tone K, Tsutsui M. Malmquist index in the dynamic & network DEA model. Oper Res Soc Japan. 2013;31:206–7.

    Google Scholar 

  46. Torres P, Ferreira J, Monteiro A, Costa S, Pereira MC, Madureira J, Mendes A, Teixeira JP. Air pollution: a public health approach for Portugal. Sci Total Environ. 2018;643:1041–53.

    CAS  PubMed  Google Scholar 

  47. Vlaandern JJ, Janssen NA, Hoek G, Keski-Rahkonen P, Barupal DK, Cassee FR, Gosens I, Strak M, Steenhof M, Lan Q, Brunekreef B, Scalbert A, Vermeulen RCH. The impact of ambient air pollution on the human blood metabolome. Environ Res. 2017;156:341–8.

    Google Scholar 

  48. Wang Q, Zhang C, Cai W. Factor substitution and energy productivity fluctuation in China: a parametric decomposition analysis. Energy Policy. 2017;109:181–90.

    Google Scholar 

  49. Wang ZX, Hao P, Yao PY. Non-linear relationship between economic growth and CO emissions in China: an empirical study based on panel smooth transition regression models. Int J Environ Res Public Health. 2017;14(12):1568.

    PubMed Central  Google Scholar 

  50. World Energy Statistics. International Energy Agency. 2017 edition.

  51. Wu S, Wang B, Yang D, Wei H, Li H, Pan L, Huang J, Wang X, Qin Y, Zheng C, Shima M, Deng F, Guo X. Ambient particulate air pollution and circulating antioxidant enzymes: a repeated-measure study in healthy adults in Beijing, China. Environ Pollut. 2016;208:16–24.

    CAS  PubMed  Google Scholar 

  52. Xie Y, Dai H, Dong H, Hanaoka T, Masui T. Economic impacts from PM2.5 pollution-related health effects in China: a provincial-level analysis. Environ Sci Technol. 2016;50(9):4836–43.

    CAS  PubMed  Google Scholar 

  53. Yang BY, Qian ZM, Vaughn MG, Nelson EJ, Dharmage SC, Heinrich J, Lin S, Lawrence WR, Ma H, Chen DH, Hu LW, Zeng XW, Xu SL, Zhang C. Dong GH is prehypertension more strongly associated with long-term ambient air pollution exposure than hypertension? Findings from the 33 Communities Chinese Health Study. Environ Pollut. 2017;229:696–704.

    CAS  PubMed  Google Scholar 

  54. Yao X, Guo C, Shao S, Jiang Z. Total-factor CO2 emission performance of China’s provincial industrial sector: a meta-frontier non-radial Malmquist index approach. Appl Energy. 2016;184:1142–53.

    CAS  Google Scholar 

  55. Zhang N, Choi Y. Environmental energy efficiency of China’s regional economies: a non-oriented slacks-based measure analysis. Soc Sci J. 2013;50(2):225.

    Google Scholar 

  56. Zhao P, Li S, Li P, Liu J, Long K. How does air pollution influence cycling behavior? Evidence from Beijing. Transp Res Part D Transp Environ. 2018;63:826–38.

    Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by the National Natural Science Foundation of China (No. 71773082), Sichuan Science project, No. 2020JDR0079 and the Fundamental Research Funds for the central Universities (Grants No. SCU-BS-PY201016)

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, YL and YHC; Methodology, YL; Software, YHC; Validation, YL, YHC, and T-YL; Formal Analysis, YL; Resources, T-YL; Data Curation, YL; Writing—review and editing, YL; Visualization, T-YL; Supervision, Y-HC; Funding Acquisition, YL. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yung-Ho Chiu.

Ethics declarations

Ethics approval and consent to participate

This manuscript is a commentary article which has no human data collection.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendices

Appendices

Appendix 1

See Table 7.

Table 7 Overall efficiency by city and region from 2013–2016

Appendix 2

See Table 8.

Table 8 Production efficiency and healthcare resource utilization efficiencies from 2013–2016

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Lin, TY. & Chiu, YH. Dynamic linkages among economic development, environmental pollution and human health in Chinese. Cost Eff Resour Alloc 18, 32 (2020). https://doi.org/10.1186/s12962-020-00228-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12962-020-00228-6

Keywords