Air Pollutant and Health-Efficiency Evaluation Based on a Dynamic Network Data Envelopment Analysis

Environmental pollution and the associated societal health issues have attracted recent research attention. While most research has focused on the effect of air pollution on human health and local economies, few articles have discussed the environment, health, and economic development in in an integrated analysis. This paper used a Dynamic Network SBM Model to evaluate production and health efficiencies in Chinese cities and found that the production efficiency scores were slightly higher than the health efficiency scores, with the two-stage efficiency scores in most cities having significant fluctuations. Labor, fixed assets, energy, GDP, and lung disease and mortality reduction efficiencies in the first stage were generally high; however, the medical input efficiencies in the second stage were low, indicating that there was there significant room for improvement in many cities.


Introduction
Air pollution is one the most worrying environmental issues in developing countries. Since 2013, the atmospheric pollution; carbon dioxide, sulfur dioxide, and nitrogen oxides; resulting from rapid industrial development has damaged the health of local communities. China's rapid economic development over the past few decades has resulted in serious air pollution and greenhouse gas emissions. Table 1 shows the emissions from six major air pollutants; SO 2 , NO 2 , PM 10 , PM 2.5 , O 3 and CO; in mainland China (not including Hong Kong, Macau and Taiwan) from 2013, when the newest air quality standards were put into practice, to 2017 [1]. As can be seen almost all air pollutants except for O 3 fell, with the reduction in SO 2 being the highest (55.5%). These reductions are the result of multiple government interventions to combat air pollution; the three-year plan for defending the blue sky, the promotion of green development, and enhanced environmental law enforcement. However, air pollution remains a serious problem in China as even after these significant interventions, the PM 10 and PM 2.5 concentration levels have not yet reached the national secondary standard, therefore, more concerted efforts are needed.  10  Data collected by the World Meteorological Organization (WMO) Global Atmosphere Monitoring System show that the CO 2 , methane, and nitrous oxide concentrations today are 145%, 257%, and 122% higher than before industrialization [2]. Excessive carbon dioxide emissions have been a major contributor to climate change, which has resulted in polar ice reductions, sea level rises, and soil salinization in coastal areas. The temperature increases have also resulted in decreased mountain snow cover and reduced water resources. After signing the Paris Agreement, many countries have implemented policies to reduce emissions in an effort to reduce the climate change effects.
Past research on air pollution has focused on association assessments [3,4], risk assessments [5,6], health impact assessments [7], life cycle assessments [8,9], or the effect of air pollution on public health and the social economy in terms of the relationship between environmental pollution, diseases, and social economics. However, there has been little research that has examined the coordinated effect of the environment, health, and economic development from an efficiency perspective. To address this research gap, this paper used a Dynamic Network SBM (DNSBM) Model to examine data from 2013 to 2016 in 31 Chinese cities and assess the combined environmental and health efficiencies. In the first stage (production stage), labor, fixed assets, and energy consumption are the input indicators and GDP is the output indicator, with the links to the health stage being CO 2 and AQI variables. In the second stage (health stage), health expenditure is the input and birth rate is the output variable, with the carry-over variables being the respiratory diseases and the mortality rate reductions.
The remainder of this paper is organized as the follows: Section 2 gives a comprehensive literature review, Section 3 describes the research method, Section 4 gives the comprehensive empirical results and discussion, and Section 5 gives the conclusions and policy proposals.

Air Pollution Intervention, Public Health and the Social Economy
Air pollution has been identified as a significant risk factor for various respiratory and cardiovascular diseases. A meta-analysis [10] found that outdoor air pollution was a major contributor to the increase in chronic obstructive pulmonary diseases (COPD). Clinical evidence from China, the United States and the European Union concluded that a 10 µg/m 3 increase in PM 10 led to a significant rise in COPD exacerbation and deaths. In another study, [11] epidemiological and pathobiological evidence confirmed that both short and long-term exposure to air pollution increased stroke risk, and in other studies, it was found that air pollution could also be associated with other diseases and symptoms, such as still births [12], depression [13] and Alzheimer's [14].
Given the serious air pollution impacts on public health, a rise in pro-environment attitudes in society has led governments to consider a willingness-to-pay (WTP) to improve air quality [15]. Four main strategies have been suggested to improve air quality; general air quality control regulations, traffic related emissions controls, energy related emissions controls, and greenhouse gas emissions controls [16]. However, as air quality interventions can result in social and economic problems such as inconvenient traffic, factory bankruptcies, and unemployment [17], governments need to weigh up the costs and benefits associated with air pollution interventions so as to balance air pollution controls and socio-economic development.
Generally, joint assessments of air pollution interventions, public health, and the socio-economic effects need to include four assessment types; association assessments, risk assessments, health impact assessments, and efficiency assessments.

Association Assessments
An association assessment uses mathematical modeling and statistical analysis to characterize the distributional associations between air pollution, public health and social economy. Luo et al. [3] analyzed the relationship between economic conditions and air pollutants in large Chinese cities and found that there was a U-shaped curve for the relationship in which the environmental quality first declined and then improved with income growth. Greenfield et al. [4] found that that while the environmental and socioeconomic factors were both significantly associated with the disease burden, the socio-economic factors were more closely associated with the disease burden than the environmental factors.
A common misuse of association assessments can been to make causal inferences based solely on the associations. Technically, association assessments involve an examination of the statistical distributions of multiple factors; however, causal inferences focus on the treatment effects of certain interventions [18]. This difference between association and causation can be exemplified by differences between sunglasses and sunstroke; that is sunglasses can be associated with sunstroke, but do not cause sunstroke. As air pollution control decision making requires a proof of causation rather than association, it is vital to keep these differences in mind.

Risk Assessment
Risk assessment is a type of causal inference; for example for air pollution, it involves an identification of the air pollution risks, the potential consequences, the occurrence probabilities, the tolerability or acceptability of the air pollution risks, and methods for mitigating or reducing air pollution risk probabilities.
One of the most widely used risk assessment approaches has been a comparative risk assessment framework that estimates the levels and trends in exposure, deaths, and disability-adjusted life-years (DALYs) from risk factors that include but are not limited to air pollution. Smith et al. gave details about the methods used in the comparative risk assessment of household air pollution, such as how the exposure-response model was built to estimate the relative risks.
However, risk assessment models for outdoor air pollution are far less well developed [5,6], primarily because outdoor air pollution exposure is much more difficult to measure. Dreaves claimed that there was less focus on this type of risk assessment because of poor stakeholder involvement and participation [7].

Health Impact Assessment
Health impact assessments (HIA) are a combination of procedures, methods and tools to assess the potential effects of a policy, program, or project on the health of a population and the distribution of those effects within the population [7].
There have been many methodological HIA iterations, and while there has not been a single agreed HIA method, it is generally accepted that an HIA should include [19] screening (establishing whether an HIA is required), scoping (planning what to do and how to do it), appraisal (identifying the health hazards and considering impact evidence), reporting (developing recommendations to reduce hazards and/or improve health), and monitoring (monitoring the implementation of the proposal and evaluating whether the HIA has influenced the decision making process). For instance, the Aphekom project used traditional HIA methods to assess the health impact of air pollution in 25 European cities, finding that life expectancy and monetary benefits increased significantly beyond the achievements of current EU legislation when fine particle levels were reduced in European cities and pollution was regulated near busy roads.
HIAs can also have drawbacks. As HIAs emphasize the outcomes (impact) of a decision process, when confronted with complicated real-world cases, they lack adequate methods to measure and quantify the inputs, making it difficult to holistically evaluate the efficiency of each proposal. Birley, for example, criticized HIAs for failing to validate scientific hypotheses and stop proposals (however inappropriate) [20].

Life Cycle Assessment
As there are many interventions available to reduce the negative effects of air pollution on the environment and public health, it is necessary to select the most efficient. Life-cycle assessments (LCA) have been a key efficiency assessment method as these are able to evaluate the potential impacts associated with air pollution interventions and interpret the results to help make informed decisions. The LCA generally had four main phrases; goal and scope definitions, an inventory analysis, an impact assessment, and interpretation [21].
A recent study examined the environmental effects of PM 2.5 on human health and identified this factor as an important for life cycle impact assessments [8]. The effect of PM on health in the U.S. in 2005 was estimated at 130,000 deaths and 2 million DALY. In another study, Damgaard et al. [9] compared eight different air pollution control technologies using lifecycle-assessment modelling and found that the potential environmental air emissions impacts had decreased over the last 35 years and that these impacts could be partly or fully offset by recovering the energy produced from fossil fuels.
However, as LCA methods are generally only able to assess a single exposure type or intervention, they are unable to comprehensively assess both economic and health efficacy at the same time. Further as LCA allows researchers to choose the parameters and data, there is a risk of conflict when seeking to compare assessments as a favored proposal in one study may be disfavored in another study because of different parameters and data.

Summary and Implications of the Literature Review
The literature review revealed that air pollution is an environmental, health, and socio-economic issue. Previous research has provided sound evidence for the relationship between air pollution, health, and the social economy and given examples of successful policies and regulations that have sought to reduce the negative effects. However, even though there are available references, experiences, and paradigms to assess air pollution control in China, it is necessary to determine which if any are the most appropriate at this time. In this paper, however, we have chosen to analyze effectiveness using efficiency analyses, as these have been insufficient in previous research. Therefore, to examine the association between air pollution, health, and the social economy, this study proposes a new method to holistically evaluate socio-economic and health efficiencies under the same framework.

Research Methods
As in Tone and Tsutsui [22], a two-stage Dynamic Network DEA model was designed for the input and output efficiency indicators. In this section, the Dynamic Network DEA model is analyzed first, after which the two stage input and output efficiency indicators are described and assessed.

Dynamic Network DEA
Data Envelopment Analysis is based on a Pareto optimal solution concept and uses linear programming techniques to evaluate the relative efficiencies of Decision Making Units (DMU). The first use of this type of research method was 1957 in Farrell's [23] study on the "Measurement of Productive Efficiency". However, Farrell's efficiency assessment model was only suitable for a single input and a single output. Because such research generally involves multiple inputs and multiple outputs, Charnes, Cooper and Rhodes [24] proposed the CCR model in 1978, which extended Farrell's model to allow for multiple inputs and multiple outputs, with the optimal solution being solved using linear programming, which is now known as Data Envelopment Analysis. In 1984, Banker, Charnes and Cooper [25] proposed the BCC model that included a variable return to scale (VRS) assumption to replace the constant return to scale (CRS) assumption. Tone [26] then proposed a Slacks-Based Measure (SBM) to allow for any output and input slacks. Using non-radial estimation and a single scalar to present the SBM efficiency, the efficiency values are typically between 0 and 1; when the efficiency of a decision making unit is 1, there is no input or output slack on the production frontier.
Färe et al. [27] proposed Network Data Envelopment Analysis (NDEA) in 2007, in which the production processes were broken up into secondary production technologies called Sub-DMUs, and traditional CCR or BCC models employed to determine the optimal solutions. In 2013, Tone and Tsutsui [22] proposed the SBM (weighted slack-based measures) Dynamic Network DEA, which was based on the links between the decision-making units, included each department as a Sub-DMU, and had carry-over link activities.

Dynamic Network DEA Model and Solution:
It is assumed that the number of DMUs is n (j = 1, ..., n), with each DMU being divided into a number of k, (k = 1, . . . , K), and time periods t, (t = 1, . . . , T). Each DMU has an input and output in period t through a carry over (link) to the next period t + 1.
Set m k and r k as the input and output for each division K, in which (k,h)i indicates division k to h, and L hk denotes the set of k and h.
If part of the output is not good, it is considered an input to division k.

Period and Division Efficiencies
Period efficiency: Division efficiency: s (t,(t+1)) ok l good z (t,(t+1)) ok l good Division Period efficiency: From the above results, the overall efficiency, period efficiency, division efficiency and division period efficiency can be determined. Hu and Wang's [28] total-factor energy efficiency index is employed here to overcome any possible biases in the traditional energy efficiency indicators. There are nine key features used in this present study: fixed assets, labor, energy consumption, GDP, health expenditure, birth rate, respiratory diseases and death rate. In our study, "I" represents area and "t" represents time.

Fixed Asset Efficiency
Fixed asset efficiency is the ratio of target fixed asset input to actual fixed asset input, the model for which is: If the target fixed asset input is equal to the actual input level, then the fixed asset efficiency equals 1, indicating efficiency; however, if the target fixed asset input is less than the actual input level, then the fixed asset efficiency is less than 1, indicating inefficiency.

Labor Efficiency
Labor efficiency is the ratio of target labor input to actual labor input, the model for which is: If the target labor input is equal to the actual input level, then labor efficiency equals 1, indicating efficiency; however, if the target labor input is less than the actual input, then the labor efficiency is less than 1, indicating inefficiency.

Energy Consumption Efficiency
Energy consumption efficiency is the ratio of target energy input to actual energy input, the model for which is: If the target energy input is equal to the actual input level, then energy consumption efficiency equals 1, indicating efficiency; however, if the target energy input is less than the actual input level, then the energy consumption efficiency is less than 1, indicating inefficiency.

GDP Efficiency
GDP efficiency is the ratio of actual desirable GDP output to target desirable GDP output, the model for which is: If the target desirable GDP output is equal to the actual desirable GDP output level, then the GDP efficiency equals 1, indicating efficiency. If the actual desirable GDP output is less than the target desirable GDP output level, then the GDP efficiency is less than 1, indicating inefficiency.

Health Expenditure Efficiency
Health expenditure efficiency is the ratio of target health expenditure input to actual health expenditure input, the model for which is: If the target health expenditure input is equal to the actual health expenditure input level, then the health expenditure efficiency equals 1, indicating efficiency; however, if the target health expenditure input is less than the actual health expenditure input level, then the health expenditure efficiency is less than 1, indicating inefficiency.

Birth Rate Efficiency
Birth rate efficiency is the ratio of actual desirable birth rate output to target desirable birth rate output, the model for which is: If the target desirable birth rate output is equal to the actual desirable birth rate output level, then the birth rate efficiency equals 1, indicating efficiency; however, if the actual desirable birth rate output is less than the target desirable birth rate output level, then the birth rate efficiency is less than 1, indicating inefficiency.

Respiratory Disease Efficiency
Respiratory disease efficiency is the ratio of target undesirable respiratory disease output to actual undesirable respiratory disease output, the model for which is; Respiratory If the target undesirable respiratory disease output is equal to the actual undesirable respiratory disease output, then the respiratory diseases efficiency equals 1, indicating efficiency; however, if the target undesirable respiratory disease output is less than the actual undesirable respiratory disease output, then the respiratory disease efficiency is less than 1, indicating inefficiency.

Death Rate Efficiency
The death rate efficiency is the ratio of target undesirable death rate output to actual undesirable death rate output, the model for which is; Death Rate efficiency = Target Death Rate Undesirable output (i, t) Actual Death Rate Undesirable output (i, t) If the target undesirable death rate output is equal to the actual undesirable death rate output, then the death rate efficiency equals 1, indicating efficiency; however, if the undesirable target death rate output is less than the actual undesirable death rate output, then the death rate efficiency is less than 1, indicating inefficiency.

Data Sources and Description
This study used panel data from 31 of the most developed cities in eastern and western China.    Figure 2 shows the statistical analysis of the input-output indicators. From 2013 to 2016, the growth in employment (em) was relatively slow, so the four years trend was slightly upward. The fixed assets input was large and significantly increased over the years. The maximum employed population increased slightly and the average and minimum employed population fluctuated.

Input-Output Index Statistical Analyses
While there were some fluctuations in energy consumption, the average growth was not significant and decreased in 2016, with the minimum decreasing significantly after 2014; however, maximum energy consumption had an upward trend, particular in 2016 and 2017. The gap between the maximum and minimum energy consumption was therefore increasing.
The government health expenditure input indicators also fluctuated, with the average increasing from 2013 to 2015 and falling in 2016. The maximum government Health Expenditure input in 2014 was slightly lower than in 2013, but increased significantly in 2015 and continued to increase in 2016. The birth rate maximum had a continuous upward trend with a significant overall increase; however, while the mean birth rate had a fluctuating upward trend, the minimum birth rate had a downward trend; therefore, the gap between the maximum and minimum birth rate was increasing.
The total GDP rose significantly over the years. While there was a significant increase in the maximum GDP, the average GDP was relatively flat, and the minimum GDP tended to fluctuate. was slightly lower than in 2013, but increased significantly in 2015 and continued to increase in 2016. The birth rate maximum had a continuous upward trend with a significant overall increase; however, while the mean birth rate had a fluctuating upward trend, the minimum birth rate had a downward trend; therefore, the gap between the maximum and minimum birth rate was increasing.
The total GDP rose significantly over the years. While there was a significant increase in the maximum GDP, the average GDP was relatively flat, and the minimum GDP tended to fluctuate.   Table 2 shows the total efficiency scores and ranks for 31 Chinese cities from 2013 to 2016. Five cities had total efficiency scores of 1; Beijing, Fuzhou, Guangzhou, Lhasa, Urumqi, and Yinchuan; Changsha, Nanning, and Wuhan had total efficiency scores around 0.95, and the total efficiency of the other cities was below 0.7. The lowest five ranked cities were Shijiazhuang, Haikou, Harbin, Chengdu, and Tianjin, all of which had total efficiency scores below 0.   Table 3 shows the total efficiency scores and total ranks, the average efficiency scores, and the ranks in the two stages. Cities with efficiency below 0.5 in the first stage; Guiyang, Kunming, Shijiazhuang, Taiyuan, Xi'an, and Xining; were classified into one category, and cities with efficiency at 0.5 or greater were classified into another category. The cities with efficiencies greater than or equal to 0.5 in the first stage were embedded into the second stage; cities with efficiencies below 0.5 in the second stage were classified as "a city with an efficiency greater than 0.5 in the first stage but an efficiency less than 0.5 in the second stage" (In Table 3, these cities were Chengdu, Chongqing, Harbin, Haikou and Tianjin). "There were 20 (64.52% of our observations) cities with efficiencies greater than 0.5 in both stages". After classifying the 31 provincial capital cities based on the above, the results were as follows:

Total City Efficiency Scores for Each Year
(1) Characteristics of cities with efficiencies less than 0.5 in the first stage From Table 3, there were six cities with efficiencies below 0.5 in the first stage; Guiyang, Kunming, Shijiazhuang, Taiyuan, Xi'an, and Xining. The main features of these cities were single industrial structures, a lack of characteristic industries, and with the main industries being labor intensive with high energy consumption. Therefore, these cities lack innovation and competitiveness, need to adjust their industrial structures away from traditional industries, and need to introduce new energy sources.
(2) Characteristics of cities with efficiencies greater than 0.5 in the first stage but below 0.5 in the second stage.
From Table 3, there were five cities in this category; Chengdu, Chongqing, Harbin, Haikou and Tianjin. The main characteristics were serious air pollution (Chengdu, Chongqing, and Tianjin), cities with large immigrant populations (higher immigrants from other cities increase local medical burdens), and were comparatively less developed than other municipalities (Chongqing and Tianjin). Chengdu for example has high prices, high house prices, and a large population, and due to its geographic location, it does not have abundant available resources for development, all of which places a significant burden on the economy, the environment, and the population, which in turn reduces the economic efficiency and the medical input efficiency.
(3) Characteristics of cities with efficiencies above 0.5 in both stages Twenty cities had efficiencies greater than 0.5 in both stages, with an average efficiency in the first stage of 0.8339, a standard deviation of 0.1738, and a median of 0.9350, and average efficiency in the second stage of 0.8131, a standard deviation of 0.1872, and a median of 0.8894. However, there were two main issues: 1.
There were large differences between the mean and median efficiencies in the two stages, and there were outliers (maximum or minimum). With the outliers, the median should describe central trends.

2.
The median efficiencies in the first stage were about 5% higher than in the second stage, and the total first stage efficiencies were slightly higher than in the second stage.
The two-stage efficiency analysis indicated that out of the 31 cities, 20 cities were performing well as they had efficiencies over 0.5 in both stages; however, the performance of the average production efficiency in the first stage was slightly higher than the average health efficiency in the second stage.

Two stage Relative Change Rate
A relative change rate was used to assess the efficiency score change trends in the two stages; relative change rate = (second stage efficiency score − first stage efficiency score)/first stage efficiency score. A negative relative change rate indicated that the city's second stage efficiency was lower than the first stage efficiency. Conversely, a positive relative change rate indicated that the city's second stage efficiency was higher than the first stage efficiency. Table 4 shows the relative change rates for the 20 cities with two stage efficiencies greater than 0.5. As can be seen in Table 4 there were are eight cities with negative relative change rates, with a median relative change rate value of −0.1184. There were six cities with a positive relative change rate, with a median relative change rate value of 0.0626, which indicated that the efficiency improvements from the first stage to the second stage were less than the reduction in efficiency. Of the cities with negative relative change rates, Shanghai (−0.4472) and Hangzhou (−0.20431) had outliers. Of the cities with a positive relative change rate, the city with the largest change was Lanzhou (0.5673). As Lanzhou is a western city, this indicated that efficiency improvements are not necessarily reliant on the amount of capital. In fact, without the necessary adjustment measures, if the city were inefficient, capital accumulation could result in a waste of resources, as can be seen in the total efficiency scores in Table 5. The total efficiency scores in the western region cities of Lhasa, Urumqi, and Yinchuan were 1, while the total efficiency scores in the eastern region cities of Haikou, Hangzhou, and Jinan were very low (0.4167, 0.6063, and 0.6018), which indicated that the more developed eastern regions were not performing better than the western regions in terms of their rational use of resources. Therefore, there is an urgent need in China to optimize energy, economic, and medical input efficiencies.  Figure 3 shows that from 2013 to 2016, the cities in which the efficiencies continued to rise in the first stage were Guiyang, Hangzhou, Kunming, Nanchang, Nanjing, Xian, and Zhengzhou. However, while the other cities experienced fluctuations, the overall efficiency in the first stage in 2016 was higher than in the previous three years, indicating that in general, the production input and output efficiencies in the first stage improved.
The  Figure 4 shows each city's efficiency for the production input and output indicators. Of the firststage inputs (employment, fixed assets, and energy consumption), the average employment efficiency was the highest at more than 0.9, the average fixed asset efficiency was the lowest at around 0.77, the average energy consumption efficiency was around 0.85, and the average GDP efficiency was around 0.78. Overall, employment efficiency was the highest; except for Chongqing and Harbin in which the efficiencies were lower than 0.7, most cities had efficiencies higher than 0.8, and 18 cities has efficiencies of 1. However, there were only 9 cities with overall efficiencies of 1. Generally, in most cities, the fixed asset and energy consumption input needed small adjustments.  Figure 4 shows each city's efficiency for the production input and output indicators. Of the first-stage inputs (employment, fixed assets, and energy consumption), the average employment efficiency was the highest at more than 0.9, the average fixed asset efficiency was the lowest at around 0.77, the average energy consumption efficiency was around 0.85, and the average GDP efficiency was around 0.78. Overall, employment efficiency was the highest; except for Chongqing and Harbin in which the efficiencies were lower than 0.7, most cities had efficiencies higher than 0.8, and 18 cities has efficiencies of 1. However, there were only 9 cities with overall efficiencies of 1. Generally, in most cities, the fixed asset and energy consumption input needed small adjustments.

Efficiency Scores and Rankings for Birth Rate and Medical Input from 2013 to 2016
Table 5 and Figure 5 shows the birth rate and medical input efficiency scores in each city from 2013 to 2016. As can be seen, there were significant differences in the health input efficiencies, with most cities being very low; for example, Chengdu, Harbin, Xian, and Zhengzhou's medical input efficiencies were less than 0.3. The differences in health input efficiencies were greater than the differences in the birth rate efficiencies, which tended to be high in most cities, with only Harbin, Huhehot, and Tianjin having an efficiency lower than 0.8. In 2016, the medical input and birth rate efficiencies rose significantly compared to the previous three years. However, as the average medical input efficiency was low, the differences between the cities needs to be addressed.  Table 5 and Figure 5 shows the birth rate and medical input efficiency scores in each city from 2013 to 2016. As can be seen, there were significant differences in the health input efficiencies, with most cities being very low; for example, Chengdu, Harbin, Xian, and Zhengzhou's medical input efficiencies were less than 0.3. Huhehot, Lanzhou, Nanning, Taiyuan, Xining The differences in health input efficiencies were greater than the differences in the birth rate efficiencies, which tended to be high in most cities, with only Harbin, Huhehot, and Tianjin having an efficiency lower than 0.8. In 2016, the medical input and birth rate efficiencies rose significantly compared to the previous three years. However, as the average medical input efficiency was low, the differences between the cities needs to be addressed.   Table 6 shows the respiratory disease and death rate reduction efficiency scores in each city. The respiratory disease reduction efficiencies in Beijing, Changsha, Fuzhou, Guangzhou, Lhasa, Wuhan, Urumqi, and Yinchuan were 1, indicating high efficiency. Chengdu had the lowest respiratory disease reduction efficiency at below 0.4 in all four years, Chongqing was slightly better at below 0.6, and Hangzhou was slightly higher at close to 0.7. The death rate reduction efficiency in Changsha in 2016 dropped slightly; however in all other years as in nine other cities, it was 1; therefore, there is no room for improvement. Although the death rate reduction efficiency was slightly higher than the respiratory disease reduction efficiency in Chengdu and Chongqing, it was still low at around 0.6. Xian had a death rate reduction efficiency of around 0.7 and Guiyang, Hangzhou, Kunming, Nanchang, Nanjing, and Shijiazhuang had efficiencies of around 0.8. Therefore, generally, the respiratory disease and death rate reduction efficiencies were high (more than 0.7); however, Chongqing and Chengdu need improvements.

Conclusions
This paper used a DNSBM model to analyze two stages; production efficiency and health efficiency; in 31 provincial capital cities in China, and identified the areas that needed improvement as a reference for reform, health promotion, and socio-economic development. The following conclusions were made from the overall analysis: (1) As only five cities had overall efficiencies of 1, 26 cities need improvement. The total efficiencies and the efficiencies in the two stages varied widely across the cities, which could not be explained by the eastern and western regional differences. For example, some western cities (such as Lhasa, Urumqi and Yinchuan) had higher efficiency in both stages than some eastern cities (such as Hangzhou and Haikou), which indicated that cities with better economic development do not necessarily have higher efficiency. These results clearly indicated that efficiency improvements depend on the overall development of all aspects of society such as politics, culture, education, and medical care. The analysis in this paper demonstrated that even though the eastern cities are known to have higher economic development, this does not mean that they were efficient, and pointed to the fact that China still needs to focus on economic, energy, and medical input efficiency improvements. (2) The two stage efficiency scores in most cities had significant fluctuations, with only a few cities having a continuous rise (Guiyang, Kunming and Nanjing). This indicated again that overall efficiency is influenced by many factors. From 2013 to 2016, China experienced many major events, such as the "Internet +" action plan, industrial restructuring, supply structure reform, and "Health China 2020". The effect of these major events and the many practical problems such as limited industrial structures, a lack of particular industries, and a lack of innovation and competitiveness, indicated the many facets of efficiency. Further, even though the central and local governments are seeking measures to improve efficiency, it is not an easy task. Therefore, it is necessary to develop effective measures to fit China's national and regional conditions. (3) From the analysis of the individual stage input and output efficiencies, it was evident that there needs to be improvements in government medical expenditures, with special focus needed on respiratory disease and death rate decreases, especially in Chongqing and Chengdu. Although the Chinese government's medical input has gradually increased in recent years, the medical input efficiency had not significantly improved. Therefore, in further health care reforms, there needs to be a more rational allocation of medical resources. According to Yip and Hsiao [28], China's primary health care is very weak and ineffective in terms of disease prevention, consultation and management, patient referrals, and medical coordination (especially in the prevention and control of non-contagious diseases). Yip and Hsiao [29] used diabetes as an example, pointing out that the hospitalization rate for diabetes complications in China was five times higher than the OECD average.

Policy Recommendations
Based on the results of this study, the following policy recommendations are proposed: First, supply side structural reform is necessary. Local governments need to optimize their industrial structures rather than only pursuing GDP growth; therefore. The government should support innovation and enhance industrial competitiveness.
Second, local governments need to ensure the clean use of traditional energy as well as promoting the use of new energy. Further, when adjusting the energy structure, the government needs to focus on energy-conservation, emissions reductions, ecological protection, and pollution control to reduce the negative environmental impacts of economic development.
Third, China should improve medical investments, regulations, guidance, and supervision. Disease prevention needs to be promoted and national medical funding control policies implemented to ensure precise medical investment management and improved medical input efficiency.

Conflicts of Interest:
The authors declare no conflict of interest.