Human development associated with environmental quality in China

This paper aims to investigate the connection between overall environmental quality and human development. Based on China’s provincial panel data from 2004 to 2017, this study constructed the Environment Degradation Index (EDI) and Human Development Index (HDI) to measure environmental pollution and human development, respectively, and it used the Simultaneous Equations Model (SEM) to assess the relationship between them. The results showed that there was an inverted U-shaped relationship found between EDI and HDI, and the coefficients of the first and second power of HDI were 5.2781 and -2.3476, respectively. Meanwhile, the results also confirmed that environmental pollution, in turn, delayed regional economic growth, and every 0.01 unit increase in EDI was correlated with a 3.15% decrease in GDP per capita. It is recommended that the government should speed up human development to surpass the turning point of the inverted U-shaped curve soonest possible.


Introduction
The traditional environmental Kuznets curve (EKC) depicts the relationship between economic development and a specific type of pollutant, which does not reflect how environmental quality overall changes with human development. On the one hand, the development of human society not only requires economic growth, but it also demands progress in other aspects of life like education and medical treatment. On the other hand, environmental quality is affected by many kinds of pollutants, such as air pollution, water pollution, and solid pollution, and a single pollution index cannot reflect environmental status fully [1]. Therefore, to understand the entire set of characteristics of environmental quality at different stages of human development, it is vital to building composite indicators for both environmental pollution and human development to explore how emission pollution changes with human development.
Although China's economic and social development has made great progress in recent decades, it has also led to the deterioration of environmental quality, which poses a threat to sustainable development. Environmental problems that appeared in developed countries during the industrialization process in the past hundred years are now prevalent in China, which exhibit structural, compound, and compression characteristics. Both Chinese people and the government have realized the harm caused by pollution and have already begun to pay attention to environmental protection. The government has taken a series of measures to control emissions, which included making stronger environmental regulations, adjusting production and consumption structures, and introducing advanced technologies. These changes may eventually reverse the trend of environmental deterioration. Thus, this paper aims to determine whether there is a turning point for environmental pollution in the process of human development. By using the Environment Degradation Index (EDI) and the Human Development Index (HDI) to measure comprehensively environmental quality and human development, respectively, this study utilized the Simultaneous Equations Model (SEM) to examine the relationship between EDI and HDI empirically. The results suggested that there was an inverted U-shaped relationship between human development and environmental pollution, and the worsening of environmental pollution, in turn, delayed regional economic growth. This paper contributes to the literature in two aspects. First, unlike the traditional EKC curve that studies how specific environmental pollution changes with economic development, this paper explored how general environmental quality varied with human development, which can better explain the relationship between human activity and the natural environment. Second, it uses the Simultaneous Equations Model (SEM) to estimate the parameters, which not only dealt with the problem of endogeneity but also provided insights into the relationship among human development, environmental quality, and economic development.
The remainder of the paper proceeds as follows. The second part reviews the existing literature. The third part constructs indicators for human development and environmental pollution and introduces the framework of the research method. The fourth part is the results and sensitivity analysis, and the last part gives conclusions and policy implications.

The relationship between economic growth and environmental quality
Among existing studies on the relationship between economic growth and environmental quality, the most typical model used is the environmental Kuznets curve (EKC). The EKC hypothesis states that at the beginning of economic development, the government and market participants only focus on the growth of output and income, which leads to the overexploitation of natural resources and an increase in environmental pollution. However, when the economy surpasses a certain turning point, people's demands for environmental quality rise and they would urge the government to take measures to deal with pollution, such as changing the mode of economic growth and adjusting industrial structure. Thus, there exists an inverted U-shaped relationship between economic growth and environmental pollution [2,3]. Recently, scholars have also added other variables to the EKC model, which include international trade [4,5], industrial structure [6,7], education [8], political democracy [9], and social corruption [10,11]. In contrast, other scholars use different pollutants and different quantitative analytical methods to further test the EKC hypothesis [12,13]. However, due to different modes in economic development and industrial structures, there does not always exist a single form of relationship within EKC that applies to all countries and all pollutants.
Meanwhile, the deterioration of the environment and excessive consumption of resources also influence economic growth. On the one hand, the exhaustibility of resources restricts the sustainability of an economy [14][15][16]. On the other hand, people's demand for environmental quality increases with economic development, which, in turn, affects economic growth [17].
Besides, the worsening of environmental pollution also has an impact on output and consumption, which often exhibits negative marginal utility and positive marginal output [18,19].

The relationship between human development and environmental pollution
Human development can be seen as a process in which people expand the real freedoms they enjoy [20]. Progress in human development promotes people's educational level, and people with higher educational levels not only pursue material wealth, but also require a better living environment. Brasington and Hite [21] found that there was a complementarity between education and environmental quality in the United States. They indicated that when the accumulation of material wealth came at the cost of the environment and resources, the improvement in education caused people to exercise internal restraints from excessive resource demand. One study of an Indian city by Jalan et al. [22] showed that when the educational level of women in a family rose from 0 years to 17 years, the willingness to pay for higher water quality increased significantly from 66 rupees to 144 rupees.
Moreover, based on cross-sectional data on 174 countries, Jha and Bhanu [23] found an inverted N-shaped curve between EDI and HDI. With a similar econometric model and data from 1990 to 2004, Mukherjee [24] proved that there were nonlinear relationships between human development and environmental pollution in India, which varied with different pollutants. Costantini and Monni [25] analyzed cross-sectional data on 179 countries and found that achieving sustainable economic growth by investing in human capital did not harm the environment. Serkan [26] used data from 15 countries in the Mediterranean region for 1970-2006 that showed that human development reduced regional pollution emissions. Sapci and Shogren [27] used county-level panel data for 2007-2009 in the United States to explore the impact of air quality on human capital. Their study found that for every 1% reduction in pollution, human capital stock increased by 0.10%. In conclusion, human development is an effective way to achieve sustainable economic and environmental development under the constraints of physical resources.
In recent years, many scholars have studied China's economic and environmental issues. Fan et al. [28] and Ji et al. [29] discussed the relationship between environmental resources and China's economy from the perspective of resource endowment. Chang et al. [18] and Zheng et al. [30] studied the effect of air pollution on people's happiness and labor productivity in China. In addition, they have also tested China's environmental Kuznets curve to verify whether there is a turning point for China's environmental pollution [31][32][33][34][35]. Some studies also discussed China's environmental issues from the perspective of foreign direct investment [36,37], environmental regulation [38,39], and income distribution [40,41].
In conclusion, although studies were conducted on the relationship between the economy and the environment in China, most of them focused on economic growth. In contrast, our paper emphasizes China's environmental pollution through human development. It analyses the correlation between human development and environmental pollution based on the construction and measurement of a composite HDI and EDI. The targeted aim of this study is to provide a more comprehensive understanding of factors that influence environmental quality in China.

Model setting
The traditional EKC hypothesis explains how environmental quality changes with economic development by using GDP per capita to measure economic development. However, this paper aims to explore the non-linear relationship between EDI and HDI, where HDI contains three aspects of human society-life expectancy, average years of education, and per capita income, which can better reflect the status of the economy and human activity [23]. Following the modified EKC model of Jha and Murthy [23], we used HDI to replace GDP and constructed the regression model as: where i and t indicate the subscripts of province and time, respectively. EDI it and HDI it stand for environmental pollution and human development, respectively. Growth it is per capita income. COND it symbolizes other control variables that also affect environmental quality, and α i is the regional fixed effect. Besides, change in environmental quality, in turn, affects economic growth [15,19]. One reason for this is that with the aggravation of environmental pollution, people realize the harm caused by it and turn to higher environmental regulations, which affect economic growth by influencing enterprises' environmental costs, technological innovation, and resource utilization [17,42,43]. For another reason, environmental pollution reduces labor productivity and labor supply [18,19], which would also inhibit economic growth. Therefore, Growth it is an endogenous explanatory variable in the equation above, and the coefficients would be biased if we estimate them directly. Considering the mutual influence between environmental pollution and economic growth, this paper further introduces another equation to build a simultaneous equation model, which is an effective way to solve the problem of endogeneity. Specifically, based on the C-D production function, we took Growth as the explained variable and EDI as one of the explanatory variables to establish the following econometric model: where Growth it is measured by GDP per capita, and k represents capital stock per capita. H it stands for human capital and λ i is the regional fixed effect. Eqs (1) and (2) constitute the simultaneous equation model that we used to evaluate the effect of human development on environmental pollution. The control variables in Eq (1) include industrial structure (Strcu), technological progress (Tech), pollution control intensity (Inst), and foreign trade (Fore). When an economy is in the take-off and acceleration stage, the acceleration of industrialization and the development of secondary industries bring about serious environmental problems. When the driving force of economic growth turns to low pollution and high output sectors (e.g., services and information), the impact of the economy on the environment becomes favorable to environmental quality. Hence, we introduced the industrial structure (Strcu) into the equation, which was measured by the proportion of the output value of the second industry [44,45]. In addition, with the development of the economy, people increase investment in environmental protection technology and R&D, which allows pollution control to become easier. Thus, this study used technological progress (Tech) and pollution control intensity (Inst) as control variables, which were measured by the natural logarithm of R&D expenditure and the proportion of investment to control environmental pollution in GDP. Furthermore, international trade also influences the environment [46]. Therefore, we also took foreign trade (Fore) as one of the control variables and used the ratio of foreign capital to GDP to measure it.
Control variables in Eq (2) include capital per capita (k) and human capital (H). Referring to Zhang et al. [47], we used the perpetual inventory method to estimate the total amount of capital, and we chose 1952 as the base year because the earlier the base year, the smaller impact is the estimation error. Also, we used an average year of education to measure human capital, and the number of years of education for primary school, junior high school, senior high school, and junior college were 6, 9, 12, and 16, respectively. Finally, Growth was measured by GDP per capita. k, H, and Growth all appear in their logarithmic forms.

Estimation method
Due to the existence of endogeneity and heteroscedasticity, our study used the three-stage least square method (3SLS) to estimate simultaneous equations. In the first stage, a simplified form of the simultaneous equation system was estimated. In the second stage, the fitted values of endogenous variables were used to get the 2SLS estimates for all equations. According to the seemingly unrelated regression (SUR) technique, a residual value for each equation was used to estimate the variance and covariance between the equations after the parameters for 2SLS were obtained. In the third stage, a parameter estimator of the generalized least square was obtained. 3SLS produced more effective estimators than 2SLS because it considered the correlation between equations. In a balanced system, the estimators obtained by 3SLS are: Where, ∑ is the covariance matrix of the residuals; operator�represents the Kronecker product; Z is the pre-variable matrix; X i stands for the T×k i order explanatory variable matrix of equation i. The problem of model identification must be considered before model estimation, that is, whether parameters estimation can obtain from the estimated induced coefficients. The structural form of the simultaneous equations of this study is: where the variable matrix and the coefficient matrix are: In the simultaneous equations, Growth and EDI are used as endogenous variables. The others are predetermined variables that are assigned by external conditions. According to the order condition and rank condition of model identification, we have: where k is the number of endogenous variables. For the number of predetermined variables (g) of all equations, we have: It can be seen that the empirical model in this paper is over-identified, which means the coefficients can be estimated. Besides, we also conducted a specification test with the use of the Husman test to further examine the simultaneous problem of the equations. Specifically, Specifically, we took Growth it as the explained variable and performed a regression on the predetermined variables to obtain the estimated value Growtĥ t and estimated residual Then, we took EDI as the explained variable to regress on Growtĥ it andÛ it and performed a T-test on the coefficient ofÛ it . The result showed that the coefficient ofÛ it was significant at the 1% level, which means that there exists an interactive relationship between Growth and EDI, and the estimated parameters would be biased and inconsistent if we use the ordinary least squares as the estimation method. Therefore, the three-stage least square method (3SLS) is suitable for estimating the simultaneous equations in this paper.

Data sources and description
Based on the Human Development Indices and Indicators-2018 Statistical Update released by UNDP, we used the geometric average of the following three dimensions to measure HDI: i) health level, which was measured by life expectancy at birth. The life expectancy data for 1982, 1990, 2000 and 2010 come from the National Bureau of Statistics of China., and the data for other years were estimated by interpolation and extrapolation. ii) educational level, which was measured by years of education, and iii) living standard, which was measured by per capita national income. UNDP used the per capita income to measure the living standard. However, because of the lack of statistics during the sample period, we used the ratio between per capita GDP at the national level and GNP measured by PPP dollar price in 2011 to estimate the per capita GNP of each province. The calculation method for the sub-indicators and HDI are presented by the Eqs (12) and (13).
Among them, HX ij and HDI j represent the ith-dimension of province j and Human Development Index of province j. X ij is the actual value of dimension i province j. X i � and X i �� represent the minimum and maximum value of dimension i, respectively. The higher the value of HDI, the higher the level of human development.
HDI values and the ranking of the provinces' representative years during the sample period of 2004-2017 showed that China's human development level improved, and the average HDI value increased from 0.6322 in 2004 to 0.7374 in 2017 (Table 1). However, the relative ranking among provinces remained almost unchanged. Cities in the eastern region, such as Beijing, Shanghai, Tianjin, and Jiangsu, enjoyed a relatively higher level of human development, but cities in the western region, such as Tibet, Guizhou, Yunnan, Qinghai, and Gansu, had a relatively lower level of human development.
Referring to the construction of HDI and the actual situation in China, we selected 17 indicators to comprehensively measure environmental pollution for six dimensions: wastewater, waste gas, industrial solid waste, air quality, domestic waste treatment, and environmental selfpurification rate ( Table 2). Eqs (16) and (17) shows the specific calculation method. or: EX ij , EDIX ij, and EDI j represent the indicator value of dimension i, the value of dimension i, and the Environment Degradation Index (EDI) of province (city) j, respectively. Guizhou, Yunnan, Qinghai, and Gansu had the lowest level of human development, and their rankings for environmental quality were 9, 5, 25, and 24, respectively (Table 3). Although the human development level for Yunnan and Guizhou was relatively low, their environmental quality was the best in China. This might have been because the industrialization process in the two provinces was relatively undeveloped, and so the damage to the environment was also small. Beijing, Shanghai, Tianjin, Jiangsu, and Liaoning ranked in the top five in terms of HDI with their EDI rankings of 1, 6, 13, 8, and 22, respectively. Those provinces (cities) were highly industrialized, but the rankings of their EDI were not in the top five except for Beijing. Industrialization may be one of the factors that led to the deterioration of the environment.
A scatter diagram between HDI and EDI 7 showed there was a non-linear relationship between them (Fig 1). Besides, scatter diagrams showed there were nonlinear relationships between HDI and the six sub-dimensions of EDI (Figs 2-7).
This study did not include Hainan and Tibet provinces due to a lack of data for the two provinces. Thus, our empirical analysis was based on data for 29   Note: Air quality is measured by the statistical indicators of major cities in each province (city) due to the limit of the availability of data; (#) and (^) represent the use of Eqs (14) and (15) Table 3.

Estimation results
For Eq (1), the estimated values of the first and second powers of HDI were 5.2781 and -2.3476 with a significance level of 1% (Table 4), which proved that there was an inverted U-shaped relationship between China's environmental pollution and human development. The results indicated that increasing investment in human development had a direct and significant effect on China's environmental quality. Meanwhile, the estimation results of the six sub-indicators of EDI, which included EDIX1, EDIX2, EDIX3, EDIX4, EDIX5, and EDIX6 that stood for wastewater, waste gas, industrial solid waste, air quality, domestic waste treatment, and selfpurification rate of the environment, respectively, also showed significant correlations with HDI (Table 4). There was an inverted U-shaped relationship between HDI and exhaust gas, domestic waste treatment, and environmental self-purification rate. However, the relationship between HDI and wastewater was significantly negative, that is, the increase in HDI was favorable for reduction of wastewater discharge. The possible explanation is that in recent years, China has introduced the concept of cleaner production and has been strictly controlling sewage discharge rates. Also, local governments have been actively building sewage treatment systems. In December 2016, the Development and Reform Commission, and the Ministry of

PLOS ONE
Housing and Urban-rural Development jointly issued the 13th Five-Year Plan for the Construction of National Urban Sewage Treatment and Recycling Facilities, which required that by the end of 2020, the urban sewage treatment rate should reach 95%, and this goal had been reached in 2017. Thus, with strict regulations, the relationship between wastewater and human development may have crossed the turning point of the inverted U-shaped curve. Also, HDI was positively correlated with industrial solid waste and air quality; the rise in HDI intensified industrial solid waste pollution and air pollution. With the expansion of industrialization, the treatment of industrial solid waste became more difficult due to the complexity and diversity of wastes. Although China's comprehensive utilization rate of industrial solid waste increased gradually, there is still much room for improvement. For example, the comprehensive utilization rate in 2017 was 54%, with an additional utilization space of 46%. Meanwhile, although people were more concerned about air pollution in recent years, the improvement in air quality still faced great challenges due to the complex and long-term task of controlling air pollution. Therefore, the relationship among industrial solid waste, air quality, and HDI occurred on the left half of the inverted U-shaped curve.
The estimated results for the control variables for Eq (1) showed that data fit well with the model and most of the estimated coefficients were significant at the 10% significance level. For Growth, the estimated coefficients for EDI, EDIX1, EDIX2, EDIX3, EDIX4, EDIX5, and EDIX6 were -0.1832, 0.3287, -0.2097, -0.3062, -0.2822, -0.2814, and -0.3579, respectively, and all of them were significant at the 1% level. Thus, we concluded that China was on the right half of the EKC curve that was composed of GDP per capita and environmental pollution, which meant that environmental quality improved with the increase in GDP per capita.
Industrial structure (Struct) was positively correlated with EDI. That is, with the expansion of industrialization, environmental quality deteriorated. However, within the six sub-indicators, only waste gas, industrial solid waste, and environmental self-purification rates were significantly correlated to industrial structure. Among them, the estimated coefficients of Struct for waste gas and industrial solid waste were 0.0061 and 0.0015, respectively, with P = 10%. This result was consistent with our previous analysis that provinces with a large proportion of the secondary industry, such as Shanxi (0.5335), Henan (0.5332), Shaanxi (0.5280), Shandong (0.5260), and Qinghai (0.5254), also faced relatively higher environmental pressure, with EDI rankings of 28, 18, 20, 14, and 25, respectively. Our results confirmed most of the EKC studies [7] and indicated that industrialization was an important factor that influenced regional environmental pollution. However, there was no significant correlation found between industrial structure and air quality or domestic waste treatment. The reason may be that the impact of industrial structure on environmental pollution mainly suggested that it was industrial pollution, so impacts on other forms of pollution were not significant. Foreign trade (Fore) significantly improved overall environmental quality and its six subindicators, and regression coefficients for EDI, EDIX1, EDIX2, EDIX3, EDIX4, EDIX5, and EDIX6 were -0.0100, -0.0061, -0.0267, -0.0119, -0.0100, -0.0045, and -0.0284, respectively. The outcomes confirmed the analysis of Copeland and Taylor [46], which indicated that as a result of technology spillover effects from foreign investment, foreign trade may not necessarily lead to deterioration of environmental quality in developing countries.
The regression coefficients for technological progress (Tech) in the seven regression equations were -0.0440, -0.0126, -0.1406, -0.0685, -0.0288, -0.0082 and -0.0572, respectively, at the significance level of 10%. Technological progress, without exception, promoted the improvement of regional environmental quality in China, which indicated that increasing scientific research funding was an effective way to improve environmental quality. Meanwhile, we found that the effect of pollution control intensity (Inst) on the improvement in environmental quality was minimal, which indicated that improving environmental quality through investment in R&D beforehand was better than investment in pollution abatement afterward.
The estimated results for Eq (2) suggested that the investment of material capital (k) and the accumulation of human capital (H) significantly promoted regional economic growth in China. The coefficients for both k and H were significant at the level of 5%. In the seven equations (EDI, EDIX1, EDIX2, EDIX3, EDIX4, EDIX5, and EDIX6), the coefficients for k were The rise in educational level and the increase in knowledge-based talents brought by human development were critical to regional economic growth. In addition, comprehensive environmental pollution (EDI), waste gas (EDIX2), industrial solid waste (EDIX3), air quality (EDIX4), domestic waste pollution (EDIX5), and environmental self-purification rate (EDIX6) all had a significant negative impact on economic growth, with their coefficients being -3.1503, -0.6968, -1.9005, -4.9022, -3.4895, and -1.2149, respectively. In contrast, the coefficient of wastewater (EDIX1) was 3.7918, which indicated a positive effect on economic growth. In general, the deterioration in environmental quality significantly delayed China's regional economic growth, and the negative effects showed a certain degree of robustness. However, this study did not analyze the detailed mechanism of action for these impacts, which should be studied further.

Sensitivity analysis
In this section, we changed the calculation method for HDI to test the robustness of the results. Specifically, we measured the human development level with the HDI first proposed by UNDP in the 1990 Human Development Report, which took the arithmetic average of the following three dimensions: i) life expectancy, which was measured by life expectancy at birth, ii)

PLOS ONE
knowledge, which was measured by adult literacy rate (with a weight of 2/3) and the comprehensive enrollment rate for primary, intermediate, and advanced education (with a weight 1/3), and iii) resource acquisition, which was measured by GDP per capita. To avoid the multicollinearity problem that may occur in the econometric model, we eliminated the resource acquisition dimension following Costantini and Monni [25] and only utilized the longevity and knowledge dimension to measure human development level of the regions. According to the Human Development Report, life expectancy at birth ranges from 25 to 85 years; the adult literacy rate and the comprehensive enrollment rate of primary school, secondary school, and university range from 0% to 100%, with the adult literacy rate being the proportion of literate persons aged 15 and over in the total population aged 15 and over. In addition, because the data for the total number of people aged 6 to 22 is difficult to obtain during the sample period, we used the proportion of enrolled students (ordinary primary schools, middle schools, high schools, secondary vocational schools, and colleges) to represent the comprehensive enrollment rate.
The relationship between HDI and EDI still had an inverted U-shape after we changed the measurement method for HDI (Table 5). Also, the positive effects of production factors and the negative effect of EDI on economic growth were also consistent with the previous analysis. Thus, we concluded that the main results of this paper were robust.

Conclusions
This study verified the existence of an inverted U-shaped curve between environmental pollution and human development, and we discussed the influence of pollution on regional economic growth. Meanwhile, there were also significant correlations between the six subdimensions of EDI and HDI. For waste gas, domestic waste treatment, environmental selfpurification rate, the shape of the curve was still an inverted U, and for industrial solid waste and air quality, the relationship was positive. For wastewater, the correlation was negative, which indicated that wastewater may have been on the right side of the inverted U-shaped curve. Although the relationships were slightly different, it is indisputable that the improvement in China's human development played an important role in promoting pollution control and in improving environmental quality in the long run. As a whole, China's economic growth helped to improve regional environmental quality, and environmental deterioration significantly delayed regional economic growth. At the same time, physical capital investment and human capital accumulation were essential to regional economic growth; this emphasized the importance of education and human development.
We believe that the key to improving environmental quality is to speed up human development and surpass the turning point of the inverted U-shaped curve as soon as possible, and increasing investment in human development is vital to solving the dilemma among economic growth, resource depletion, and environmental degradation in China.

Acknowledgments
We would like to thank Thomas A. Gavin, Professor Emeritus, Cornell University, for help with editing this paper.