Next Article in Journal
A ‘Control Model’ of Social Media Engagement in Adolescence: A Grounded Theory Analysis
Next Article in Special Issue
Comparing Economics, Environmental Pollution and Health Efficiency in China
Previous Article in Journal
Workplace Bullying as Experienced by Managers and How They Cope: A Qualitative Study of Swedish Managers
Previous Article in Special Issue
Evaluating the Highway Tunnel Construction in Western Sichuan Plateau Considering Vocational Health and Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial-Temporal Effects of PM2.5 on Health Burden: Evidence from China

1
School of Economics, Sichuan University, Chengdu 610065, China
2
School of Public Administration, Sichuan University, Chengdu 610065, China
3
Andalusian Research Institute in Data Science and Computational Intelligence (DaSCI), University of Granada, 18071 Granada, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(23), 4695; https://doi.org/10.3390/ijerph16234695
Submission received: 21 October 2019 / Revised: 20 November 2019 / Accepted: 22 November 2019 / Published: 25 November 2019
(This article belongs to the Special Issue Environmental and Health Efficiency Evaluation)

Abstract

:
By collecting the panel data of 29 regions in China from 2008 to 2017, this study used the spatial Durbin model (SDM) to explore the spatial effect of PM2.5 exposure on the health burden of residents. The most obvious findings to emerge from this study are that: health burden and PM2.5 exposure are not randomly distributed over different regions in China, but have obvious spatial correlation and spatial clustering characteristics. The maximum PM2.5 concentrations have a significant positive effect on outpatient expense and outpatient visits of residents in the current period, and the impact of PM2.5 pollution has a significant temporal lag effect on residents’ health burden. PM2.5 exposure has a spatial spillover effect on the health burden of residents, and the PM2.5 concentrations in the surrounding regions or geographically close regions have a positive influence on the health burden in the particular region. The impact of PM2.5 exposure is divided into the direct effect and the indirect effect (the spatial spillover effect), and the spatial spillover effect is greater than that of the direct effect. Therefore, we conclude that PM2.5 exposure has a spatial spillover effect and temporal lag effect on the health burden of residents, and strict regulatory policies are needed to mitigate the health burden caused by air pollution.

1. Introduction

With the rapid development of China’s economy, the living standard and health of residents have been greatly improved, but ambient air pollution remains a serious problem. The Global Burden of Disease (GBD) and the World Bank ranked air pollution as the fifth and fourth global health risk factor, respectively (Global Burden of Disease (GBD), 2015; World Bank, 2016). Air pollution is associated with increased incidence of diseases (WHO, 2012). As air pollution is one of causes of death and disability, it is increasingly recognized as a worldwide public health concern [1,2,3,4,5,6,7]. Atmospheric particle matter (PM), especially those having an aerodynamic diameter less than 2.5 um (PM2.5), is considered one of the priority pollutants in air [8,9,10]. As the smog worsened, PM2.5 pollution became the main air pollutant in China and played a non-negligible role in affecting Chinese residents’ health [11,12,13]. For example, in 2018, only 35.8% of 338 cities satisfied air-quality standards in China, while the rest of the cities exceeded them. In 2018, all the 338 cities experienced on average 8.05 days with severe air pollution, and the average annual PM2.5 concentrations in 56.20% of the cities exceeded 35.74 μg/m3. Also, of the 169 key cities, 140 cities had PM2.5 as the main pollutant, and 30% cities had an air quality index (AQI) greater than 100 (Air quality index (AQI) is a quantitative description of air quality. The larger the value is, the more serious the air pollution is, and the more harmful it will be to human health. The main pollutants for air quality evaluation are particulate matter (PM2.5), inhalable particulate matter (PM10), sulfur dioxide (SO2), nitrogen dioxide (NO2), ozone (O3) and carbon monoxide (CO)) [14].
PM2.5 often contains heavy metals such as arsenic, chromium, and manganese, and PM2.5 concentrations in the region mainly depend on its energy efficiency [15,16,17,18]. Frequent exposure to PM2.5 can cause great damage to human health, and then leads to the increment of patient visits and health expense, which brings great economic burden to residents. According to China’s National Statistics Bureau (CSB), the total health expenditure increased from 359.39 billion RMB in 2008 to 1639.91 billion RMB in 2018, with an average annual growth rate of 18.73% (CSB, 2018). Share of health expenditure increased from 4.55% in 2008 to 6.57% in 2018 (Share of health expenditure refers to the ratio of total health expenditure to gross domestic product (GDP)). Per capita health expenditure rose from 1094.5 RMB in 2008 to 4237 RMB in 2018, and its average annual growth rate was 15.49% (Per capita health expenditure refers to the ratio of total health expenditure to the total population.). Besides, per capita outpatient visits increased from 3.7 times in 2008 to 6.0 times in 2018 (Per capita outpatient visits refers to the ratio of total number of outpatient visits to the total population) [19]. This suggests that the Chinese government is increasingly investing in health expenditure to improve public health and make or become less the burden of residents.
In recent years, there has been an increasing amount of literature on PM2.5 exposure and health [20,21]. Several studies have explained that PM2.5 exposure can cause many kinds of diseases, such as respiratory diseases [22,23,24,25], cardiopulmonary diseases [26,27], cardiovascular diseases [24,28], lung cancer [25,27,29,30], and brain damage [31,32,33,34], and it is even closely related to mortality [35]. Although studies have directly or indirectly proved the health effects of PM2.5 exposure, conclusions are inconsistent for the differences in the sample selection or data characteristics (such as time series, panel data, experimental data, etc). Some studies found that PM2.5 exposure can harm the residents’ health, reduce their labor capacity and shorten their life expectancy, thus further increasing health expenses, and imposing huge economic burden on the whole society [36,37,38,39]. For example, Yang et al. [37] found that the economic loss caused by PM2.5 pollution in Beijing in 2013 was 1.11 billion RMB. Zeng et al. [40] used the spatial interpolation method to explore the PM2.5 exposure in China in 2007, and found that economic loss was 1,262.5 billion RMB. Although some studies have begun to explore the impact of PM2.5 exposure on economic losses, few studies analyze the PM2.5 exposure on health burden in detail.
Moreover, with the unprecedented economic development and urbanization in recent decades in China, energy consumption has increased significantly and PM2.5 pollution has become a serious problem [41]. Because of the different economic development and urbanization in different regions in China, PM2.5 pollution varies significantly in different regions. More importantly, different regions are not independent with each other, and some phenomena in one region are closely related to the same phenomena in other regions [42,43,44,45]. From here, it is necessary to take the spatial correlation and spatial spillover effects into account when analyzing the impact of PM2.5 on health burden. If the study ignores the spatial effect, the conclusion may not be accurate. However, exiting studies have rarely considered the spatial spillover effect of PM2.5 exposure on health burden [46,47,48].To fill the gaps mentioned above, we devoted this study to explore the impact of PM2.5 exposure on the health burden from the perspective of spatial spillover effect. The main novelties and contributions of this paper were listed as follows:
(1) This study explored the impact of PM2.5 exposure on the residents’ health burden, further enriching the research perspective of economic loss brought by air pollution.
(2) The spatial econometric models were applied to examine the spatial spillover effects of PM2.5 exposure on residents’ health burden, and its spatial dependence and correlation were also discussed.
(3) The temporal lag effect was investigated in this study. Besides, outpatient expense and outpatient visits were presented to measure the health burden, as well as the number of hospitalizations was selected to test robustness.

2. Research Design

2.1. Data

This study focused on the spatial impact of PM2.5 exposure on health burden. We selected 29 regions in China from 2008 to 2017 as samples. These regions including: Beijing, Tianjin, Hebei, Shanxi, Inner Mongolia, Liaoning, Jilin, Heilongjiang, Shanghai, Jiangsu, Zhejiang, Anhui, Fujian, Jiangxi, Shandong, Henan, Hubei, Hunan, Guangdong, Guangxi, Chongqing, Sichuan, Guizhou, Yunnan, Shaanxi, Gansu, Qinghai, Ningxia and Xinjiang. As China only began releasing the ground-level PM2.5 concentrations data in 2013, this study adopted the long-term data published by Columbia University Center for Socio-Economic Data and Applications [49]. The PM2.5 concentration data in different regions of China is aerosol optical depth (AOD) data from satellite retrieval of surface PM2.5 concentration retrievals. The other indicator data were collected from the Chinese Health Statistics Yearbook 2008-2017 [50], including outpatient expense, outpatient visits, the number of hospitalizations, the number of medical institutions, the number of hospital beds and the number of doctors. The data of GDP and the ratio of urban population was from China Statistics Yearbook 2008–2017 [51].

2.2. Variable

2.2.1. Dependent Variable: Health Burden

Most existing studies only adopted the number of outpatient visits to measure the health burden of residents [52,53]. But it cannot accurately measure the health burden of residents since the results may be biased [54,55]. Also, some literature used the number of hospitalization to measure health burden [56,57]. However, Chinese residents have a habit of not seeking medical treatment if they are not seriously ill, and the number of hospitalization used may underestimate the health burden. Therefore, to ensure the reliability of the results, this study measured health burden using outpatient expense (exp_out) and outpatient visits (num_out). Outpatient expense is expressed as the ratio of the total outpatient expense to the total number of outpatient visits, and outpatient visits are expressed as the ratio of the total number of outpatient visits to the total population. Besides, a robustness test was performed by using the number of hospitalizations (num_hos) which is expressed as the ratio of the total number of hospitalizations to the total population.

2.2.2. Independent Variable: PM2.5 Exposure

PM2.5 exposure has a significant impact on the health burden of residents. This study utilized the PM2.5 concentration data to analyze the spatial impact of PM2.5 exposure on health burden. In fact, the accurate estimation of PM2.5 concentrations is one of the most critical preconditions. This study adopted the PM2.5 data in 2007–2017 published by Columbia University (The average PM2.5 concentrations and the maximum PM2.5 concentrations the of 29 regions in China from 2007 to 2017 are listed in Table A1 and Table A2). Some literature used average PM2.5 concentrations to measure PM2.5 exposure [48,58]. However, severe air pollution, such as the maximum concentrations of PM2.5, may cause more harm to residents’ health. Therefore, this study used maximum PM2.5 concentrations (PM2.5_max) to measure PM2.5 exposure. Meanwhile, due to the temporal lag effect of PM2.5 exposure on residents’ health burden [48], this study also used the maximum of PM2.5 concentrations lags by one stage (PM2.5_max(-1)) as the independent variable to verify whether the temporal lag effect exists. In addition, average PM2.5 concentrations (PM2.5_avg) was used as a substitute variable for PM2.5 exposure to test the robustness of the results.

2.2.3. Control Variable

As all know, health burden is affected not only by PM2.5 exposure, but also by many others. Referring to the existing studies [30,48], this study controlled the following variables: per capita gross domestic product (PGDP), the ratio of urban population (urban), the number of medical institutions (num_inst), the number of hospital beds (num_bed) and the number of doctors (num_doctor).
The variables and their definitions were presented in Table 1, and the descriptive statistics of all were shown in Table 2.

2.3. Method

2.3.1. Spatial Autocorrelation Test

To comprehensively explore the spatial spillover effect of PM2.5 exposure, we used the global and the local spatial correlation indices in Exploratory Spatial Data Analysis (ESDA) [59] to test the spatial correlation. The global and the local spatial correlation indices were measured by Moran’s index (Moran’s I), and their calculation formulas were shown as follows:
I = i = 1 n j = 1 n w i j ( x i x ¯ ) ( x j x ¯ ) S 2 i = 1 n j = 1 n w i j , x ¯ = 1 n i = 1 n x i , S 2 = 1 n i = 1 n ( x i x ¯ ) 2
I i = ( x i x ¯ ) S 2 j = 1 n w i j ( x j x ¯ ) , x ¯ = 1 n i = 1 n x i , S 2 = 1 n i = 1 n ( x i x ¯ ) 2
where, xi and xj represent the observed values of region i and region j, respectively; n represents the number of all regions; wij is the element in the spatial weight matri x; x ¯ is the mean value of the sample, and S2 is the variance of the sample.
For Moran’s I, its value range is [−1,1]. If its value is greater than 0, it indicates that there is a positive spatial correlation among variables. If its value is smaller than 0, it indicates that there is a negative spatial correlation among variables. Meanwhile, the values of Moran’s I in different regions can be plotted as scatter plots to view the degree of spatial agglomeration in the regions. The slope of the regression line of the scatter plot is equal to the value of global Moran’s I.
In the spatial econometric analysis, it is necessary to introduce spatial weighting matrices to describe the relationship among different regions. To systematically explore the spatial correlation characteristics among different regions in China, this study modeled the following three spatial weight matrices: spatial contiguity matrix W1, spatial distance matrix W2, and spatial economy matrix W3.
The spatial contiguity matrix is most widely used in spatial econometric analysis, but sometimes the relationship among regions is so simplified that research conclusions are biased. The elements of spatial contiguity matrix W1 were defined as follows [47]:
w i j = { 1 , i j 0 , i = j         i , j = 1 , 2 , , n
To enhance the robustness of the analysis results, we also constructed the spatial distance matrix. The elements of spatial distance matrix W2 were defined as follows [47]:
w i j = { 1 / d i j , i j 0 , i j         i , j = 1 , 2 , , n
where dij represents road distance between region i and region j.
The spatial contiguity matrix and spatial distance matrix only reflect the influence of geographical location, but do not reflect the economic correlation among regions and their influence. For example, the influence of Hebei province on Beijing municipal is much smaller than that of Beijing municipal on Hebei province. Referring to the related studies [47], this study defined the spatial economy matrix W3 as follows:
W 3 = W 1 * 1 Y ¯ d i a g ( Y 1 ¯ , Y 2 ¯ , , Y n ¯ ) Y i ¯ = 1 t 1 t 0 + 1 t = t 0 t 1 Y i t , Y ¯ = 1 n i = 1 n Y i ¯   i , j = 1 , 2 , , n
where Yit is the per capita real GDP of region i in year t; Y i ¯ represents the annual average of per capita real GDP of region i. Y is the average of Y i ¯ for all the regions. As can be seen from the formula (5), if the per capita real GDP of region i is greater than that of other regions, the region also has more influence on the other regions.

2.3.2. Spatial Econometric Model

This study used spatial econometric models to analyze the impact of PM2.5 exposure on residents’ health burden in China, and to measure the direct effect and the spatial spillover effect. Three widely used spatial econometric models were adopted to examine the spatial effects, including the spatial autoregression model (SAR), spatial errors model (SEM) and spatial Durbin model (SDM) [60]. SAR only includes the lag term of the spatial dependent variable, and SEM only includes spatial spillover effects of independent variables, while SDM includes both the lag term of the spatial dependent variable and spatial spillover effects of independent variables. Based on this, these three spatial econometric models were constructed as follows:
SAR:
Y i t = α   +   ρ W Y i t + β 1 P M 2.5 _ m a x i t + β 2 P M 2.5 _ m a x ( 1 ) i t                 + β 3 P G D P i t + β 4 u r b a n i t + β 5 n u m _ i n s t i t                 + β 6 n u m _ b e d i t + β 7 n u m _ d o c t o r i t + ε i t
SEM:
Y i t   =   α + β 1 P M 2.5 _ m a x i t + β 2 P M 2.5 _ m a x ( 1 ) i t                 + β 3 P G D P i t + β 4 u r b a n i t + β 5 n u m _ i n s t i t                 + β 6 n u m _ b e d i t + β 7 n u m _ d o c t o r i t + u i t u i t = λ W u i t + ε i t ,   ε ~ N ( 0 , σ 2 I n )
SDM:
Y i t   =   α + ρ W Y i t + β 1 P M 2.5 _ m a x i t + β 1 P M 2.5 _ m a x ( 1 ) i t                 + β 3 P G D P i t + β 4 u r b a n i t + β 5 n u m _ i n s t i t + β 6 n u m _ b e d i t                 + β 7 n u m _ d o c t o r i t + σ W X k i t + ε i t
where Y is the dependent variable; PM2.5_max and PM2.5_max(−1) are the core independent variable; PGDP, urban, num_inst, num_bed, and num_doctor are the control variables; X represents all of the above core independent variables and control variables; W is the spatial weighting matrix; εit and μit are normally distributed random error vector; α denotes the intercept item; β denotes the influence coefficient of independent variables on dependent variable; ρ denotes the spatial autoregressive coefficient; λ denotes the spatial error coefficient; θ denotes the space lag coefficient of the independent variables; i represents regions, and t represents year.

2.3.3. Model Test

In general, the Lagrange multiplier tests (i.e., LM-lag and LM-err) were used to determine a proper spatial econometric model [61], but these methods were only suitable for sectional data. For panel data, referring to the study of Belotti et al. [62], this study tested the conditions given in Table 3 to select the spatial econometric model.
In Table 3, the test results are all significant at the 1% level (p < 1), and we should reject these null assumptions that λ = 0 and λ = −ρβ. In other words, SDM cannot be simplified into SAR or SEM, and should be adopted to analyze the effect of PM2.5 exposure on residents’ health burden in China. Also, through the Hausman test (i.e., space fixed effect or time fixed effect), the result shows that the spatial econometric model should adopt fixed effect (χ2 = 11.35, p = 0.0782 < 10%). From here, this study should use the SDM model with fixed effect to analyze.

3. Spatial Distribution and Spatial Autocorrelation Analysis

3.1. Spatial Distribution

Figure 1 shows the spatial distribution of PM2.5 concentrations (Figure 1(a1,a2)), outpatient expense (Figure 1(b1,b2)), outpatient visits (Figure 1(c1,c2)) and the number of hospitalization (Figure 1(d1,d2)) in all selected regions of China in 2008 and 2017. In absolute terms, compared with 2008, PM2.5 concentrations, outpatient visits and the number of hospitalization in different regions of China increased in 2017, whereas outpatient expense varied slightly. In relative terms, Chinese regions displayed similarities during the periods. For PM2.5 concentrations, outpatient visits and outpatient expense, regions with high-values were concentrated in the eastern districts, while regions with low-values were concentrated in the western and central districts. For the number of hospitalization, regions with high-values were concentrated in the western and central districts, and regions with low-values were concentrated in the eastern districts. Due to space limitation, we only chose the spatial distribution of PM2.5 concentrations, outpatient expense, outpatient visits and the number of hospitalizations in 2008 and 2017 to analyze.
From Figure 1(a2), we can see that PM2.5 pollution is the worst in Beijing, Hebei, Tianjin, and Henan, followed by Shaanxi, Hubei, Anhui, and Xinjiang. While the regions of southwest, southeast, and northeast in China, such as Yunnan, Qinghai, Guizhou, and Fujian, have the least pollution. What stands out in this figure is that PM2.5 concentration distribution in China has obvious spatial clustering characteristics. That is, the regions with the most serious PM2.5 pollution are concentrated together, and the regions with the least pollution are also concentrated together.
Figure 1(b2,c2,d2) present the distribution of residents’ health burden in all the selected regions of China. Similar to PM2.5 pollution, the distribution of residents’ health burden in China also shows obvious spatial clustering characteristics. Figure 1(b2,c2) show that outpatient expense and outpatient visits are positively correlated with the economic level in all different regions. For example, outpatient expense and outpatient visits in developed regions are significantly greater than those in backward regions. Besides, outpatient expense and outpatient visits are also significantly higher in regions with serious PM2.5 pollution than in regions with less pollution, such as Beijing, Shandong, Jiangsu, etc. However, for Figure 1(d2), the number of hospitalization per capita in western China is higher than that in eastern China. Since the medical treatment insurance system of the western regions in China is not perfect, some patients are hospitalized to get medical insurance compensation, even if they do not meet the standards of hospitalization, resulting in false hospitalization [63]. Because of the above shortages, this study chose outpatient expense and outpatient visits as the main independent variables.

3.2. Spatial Autocorrelation Analysis

Global Moran’s I of exp_out and PM2.5_max based on the three spatial matrices from 2008 to 2017 were shown in Table 4. The results indicate that all the Moran’s I of exp_out and PM2.5_max are significantly positive at the 5% level for the three spatial matrices, except that Moran’s I of PM2.5_max are not all significant in the spatial economic matrix. These suggest that exp_out and PM2.5_max are not randomly distributed in different regions of China, but have obvious spatial correlation and spatial clustering characteristics during the study period. In a word, the residents’ health expense and PM2.5 pollution in different regions of China are characterized by high-high (H-H) value aggregation and low-low (L-L) value aggregation. Therefore, the spatial correlation must be considered when studying the impact of PM2.5 exposure on health burden of residents in China, otherwise, the research results may be biased.
To analyze the local agglomeration characteristics of all the selected regions in China, Figure 2 presents the local Moran’s I scatter plots of exp_out and PM2.5_max based on spatial contiguity matrix W1 for all the regions in 2008 and 2017. In the local Moran’s I scatter plot, the horizontal axis represents the observations of the local region, and the vertical axis represents observations of the adjacent regions. In other words, the first quadrant and the third quadrant point out the existence of positive spatial correlation, representing high-high (H-H) value clustering and low-low (L-L) value clustering, respectively. While the second quadrant and the fourth quadrant mark the existence of negative spatial correlation, representing low-high (L-H) value clustering and high-low (H-L) value clustering, respectively.
As can be seen from Figure 2, for both 2008 and 2017, most of the regions are located in the first quadrant or the third quadrant, indicating that there is a positive spatial correlation between residents’ health expense and PM2.5 concentrations in most regions. This means that the residents’ health burden and PM2.5 concentrations of different regions in China are not random but show significantly positive spatial autocorrelation. For example, for the local Moran’s I scatter plots of exp_out in 2017, there are 9 regions located in the first quadrant, such as Beijing, Tianjin, Shanghai, Jiangsu, etc., and 10 regions are in the third quadrant, including Yunnan, Sichuan, Guizhou, Gansu, etc. The others are in the second or fourth quadrants. In the scatter plots of PM2.5_max in 2017, there are 23 regions located in the first or third quadrants, such as Beijing, Tianjin, Hebei, Shanghai, etc. In a word, there is a positive spatial correlation between PM2.5 concentrations and residents’ health burden, which is consistent with the previous conclusion.
We also used hot spot analysis [64] to study the local agglomeration characteristics, and Figure 3 presents the analysis results of exp_out and PM2.5_max for all the regions in China in 2008 and 2017. By calculating the Z value (namely Getis-Ord Gi*), we can identify hot spots and cold spots with statistical significance. If the absolute value of Z of the region is smaller than 1.65, it indicates that the region is less likely to be related to its neighboring regions. While if the absolute value is greater than 1.65, it indicates that there is a close connection among the regions [64]. Furthermore, the positive Z value indicates that the region is a hot spot, that is, the values of exp_out or PM2.5_max of the region and its neighboring regions are all high. While the negative Z value indicates that the region is a cold spot, that is, the values of exp_out or PM2.5_max of the region and its neighboring regions are all low. Meanwhile, 1.65, 1.96 and 2.58 are the threshold of 10%, 5%, and 1% significance level, respectively. As can be seen from Figure 3, for residents’ health expenses, Liaoning, Jilin, and Zhejiang are hot spots in 2008, and Liaoning, Jilin, Shangdong, and Jiangsu are hot spots in 2017. Gansu, Sichuan, Yunnan, and Chongqing are cold spots in 2008, and Sichuan, Yunnan, and Chongqing are cold spots in 2017. It can be seen that hot spots are distributed in the eastern districts, and cold spots in the western districts. For PM2.5 concentrations, there are a great number of hot spots and cold spots. Hot spots are mainly distributed in the Beijing–Tianjin–Hebei regions, central districts, southeast districts, Liaoning and Jilin, whereas cold spots are mainly distributed in northwest districts and Sichuan.

4. Empirical Analysis and Discussion

This study used SDM to explain the impact of PM2.5 exposure on residents’ health burden in China. In this section, we discussed the influence of PM2.5 exposure on outpatient expense and outpatient visits, respectively. According to the test results of the models, the fitting degree of the SDM under the space and time fixed effect is superior to other models, and the following tests are based on it (The results of SDM model with space-fixed effect and time-fixed effect were also tested. The SDM model based on the space-and-time fixed effect was superior to the above two models. Therefore, this study adopted the SDM model with space-and-time fixed effect for empirical testing).

4.1. Impact of PM2.5 Exposureon Outpatient Expense

The results based on the three spatial weight matrices for the whole sample were presented in Table 5. Among them, the column (1), column (3) and column (5) in Table 5 are the results of the impact of PM2.5 concentrations (PM2.5_max) on outpatient expense (exp_out) based on spatial contiguity matrix W1, spatial distance matrix W2 and spatial economy matrix W3, respectively. The column (2), column (4) and column (6) in Table 5 are the results of the influence of one stage lag of PM2.5 concentrations (PM2.5_max(−1)) on outpatient expense (exp_out) when spatial contiguity matrix W1, spatial distance matrix W2 and spatial economy matrix W3 are adopted, respectively.
First, when the three different spatial weight matrices are used, all these coefficients of PM2.5_max are significantly positive at the 1% level. This shows that PM2.5 exposure will increase the outpatient expense of residents during the research period. Since exposure to air pollutants, such as PM2.5, causes health problems of residents, including increased incidence of respiratory and cardiovascular diseases, and then increases outpatient expense. Also, when the spatial contiguity matrix W1 and the spatial distance matrix W2 are used, the spatial lag coefficients of PM2.5_max are also significantly positive at the 5% level. While the spatial lag coefficient of PM2.5_max is positive, but it does not exceed the 10% significance level when the spatial economic matrix W3 is used. This suggests that the increase of PM2.5 concentrations in geographically close regions will cause an increase of outpatient expense, and increase the health burden of residents in a particular region. But the increase of PM2.5 concentrations in economically related regions has almost no impact on outpatient expense in a particular region. Therefore, the spatial effect of PM2.5 exposure on health burden is more influenced by neighboring or geographically close regions.
Second, the coefficients of PM2.5_max(−1) are positive at the 1% significance level when the three spatial weight matrices are used, respectively. This expresses indirectly that PM2.5 exposure has a long-term impact on residents’ health burden during the research period. Because long-term PM2.5 exposure can increase the incidence of chronic diseases such as cardiovascular disease, cancer, and diabetes [65,66,67]. As for the spatial lag of PM2.5_max(−1), the coefficients are positive at the 1% level and 5% level, respectively, when we use the spatial contiguity matrix W1 and the spatial distance matrix W2. While the coefficient is positive based on the spatial economic matrix W3, it does not pass the 10% significance level test. This points out that the increase of PM2.5 concentrations in neighboring or geographically close regions will promote the increase of outpatient expense in this region, while the increase of PM2.5 concentrations in economically similar regions will not have a significant impact on the outpatient expense in the local region.
Third, the spatial autoregressive coefficients ρ are significantly negative at the 1% and 10% levels, respectively, when the spatial distance matrix W2 and spatial economy matrix W3 are adopted. However, when we apply the spatial contiguity matrix W1, the spatial autoregressive coefficients are negative but do not pass the significant 10% level test. This suggests that there is a negative correlation of outpatient expense between a particular region and geographically or economically similar regions, and PM2.5 exposure has a spatial spillover effect on outpatient expense. That is to say, health expense in a particular region is not only affected by PM2.5 exposure in the region, but also influenced by outpatient expense in neighboring or economically connected regions. Therefore, if the spatial spillover effect and temporal lag effect are ignored, the research conclusion will be biased, and the impact of PM2.5 exposure on residents’ health burden will be greatly underestimated.
Finally, for all the control variables, the direct coefficients and the spatial lag coefficients of PGDP are significant negatives at the 10% level, suggesting that there is clear evidence for a negative correlation between PGDP and residents’ health burden. Moreover, with the improvement of economic level, the medical conditions will be improved and the possibility of getting sick is reduced to a certain extent. For the urbanization rate, the direct coefficients and the spatial lag coefficients are positive at the significant 10% level, indicating that urbanization in the region and its adjacent regions can significantly promote the increase of outpatient expense. As for the number of medical institutions on the impacts of outpatient expense, the direct coefficients are significantly negative at the 10% level based on the matrices W1 and W2, while the spatial lag coefficients are not robust. These results indicate that the increase in the number of medical institutions will increase the degree of competition, promote the institutions to improve the medical level and lower the medical price, and then make or become less the outpatient expense. For the number of hospital beds, the coefficients are all positive at significant 1% level, but the spatial lag coefficients are positive at significant 10% level when the matrix W1 is adopted. However, as for the number of doctors, the direct coefficients and the spatial lag coefficients have not passed the significant tests except for the spatial contiguity matrix W1.
Given the particularity of the spatial econometric model, we cannot directly find out the spatial spillover effect from the above models. This study calculated the direct effect, the spatial spillover effect and the total effect of PM2.5 concentrations on outpatient expense, respectively. The results were given in Table 6.
The direct effect, the spatial spillover effect and the total effect of PM2.5_max and PM2.5_max(−1) are all significantly positive at the significance 5% level. Specifically, the direct effect value of PM2.5_max is 0.0987, the spatial spillover effect value of PM2.5_max is 0.1245, the total effect value of PM2.5_max is 0.2232, and the spatial spillover effect accounts for about 55.78% of the total effect. This suggests that for every 1% increase in PM2.5 concentrations in the local region and adjacent regions, the outpatient expense will increase by approximately 9.87% and 12.45%, respectively. As for the PM2.5_max(−1), the direct effect value is 0.0942, the spatial spillover effect value of PM2.5_max is 0.1283, the total effect value of PM2.5_max is 0.2225, and the spatial spillover effect accounts for about 57.66% of the total effect. This points out that every 1% increase in PM2.5 concentrations in a particular region and adjacent regions, outpatient expense will increase about 9.42% and 12.83%, respectively. In other words, the influence of PM2.5 exposure on outpatient expense in China is mainly due to spatial spillover effect. This further proves that the spatial spillover effects of PM2.5 exposure on health burden cannot be ignored.

4.2. Impact of PM2.5 Exposure on Outpatient Visits

As the health burden of residents is a comprehensive concept, this study adopted another variable, outpatient visits, to measure the health burden of residents. The results of the impact of PM2.5 exposure on outpatient visits were presented in Table 7. Among them, the column (1), column (3) and column (5) in Table 7 are the results of PM2.5 exposure (PM2.5_max) on outpatient visits (num_out) based on the three spatial matrices from W1 to W3. The column (2), column (4) and column (6) in Table 7 are the results of the influence of one stage lag of PM2.5 exposure (PM2.5_max(−1)) on outpatient visits (num_out) when the spatial weight matrices change from W1 to W3.
First, the direct coefficients of PM2.5_max and PM2.5_max(−1) are significantly positive at the 1% level when the spatial contiguity matrix W1 and spatial distance matrix W2 are used. While the coefficients are positive but not significant at 10% level based on the spatial economy matrix W3. This indicates that PM2.5 exposure has a positive impact on outpatient visits, and the influence of PM2.5 exposure on outpatient visits has a temporal lag effect. The main reason is that, with the increase of PM2.5 concentration, the incidence of acute and chronic diseases such as respiratory system diseases and cardiovascular diseases will increase, which leads to an increase in outpatient visits.
Second, as for the spatial lag of PM2.5_max and PM2.5_max(−1), the coefficients are all negative at the 1% significance level when the spatial weight matrices change from W1 to W3, showing that there is a negative spatial spillover effect of PM2.5_max and PM2.5_max(−1) on outpatient visits in China. In other words, whether in the short or long term, the increase of PM2.5 concentrations in geographically close or economically similar regions will promote the decrease of outpatient visits in a particular region. These results are different from the spatial impact of PM2.5 exposure on outpatient expense discussed above.
Third, the spatial autoregressive coefficients ρ are significantly positive at the 1% level when the spatial matrix is W1, and significantly negative at the 1% level based on the spatial economy matrix W3. However, they do not pass the significant test when the spatial matrix is W2. This proves that the outpatient visits of adjacent regions have a positive influence on that of the particular region, while there is a negative correlation of outpatient visits between the particular region and economically similar regions. Therefore, spatial spillover effect should be considered when analyzing the impact of PM2.5 exposure on health burden of residents.
We also calculated the direct effect, the spatial spillover effect and the total effect of PM2.5 concentrations on outpatient visits of residents, respectively, and the results are shown in Table 8. The results suggest that the direct effect of PM2.5_max and PM2.5_max(−1) are significantly positive at the 1% level, while both the spatial spillover effect and the total effect are significantly negative at the significance level of 5%. These results are consistent with the above conclusions.
For the coefficient value of PM2.5_max, the direct effect value is 0.1944, the spatial spillover effect value is −0.3516, and the total effect value is −0.1572. The spatial spillover effect of PM2.5_max is much larger than the direct effect. This indicates that for every 1% increase in PM2.5 concentrations in a particular region and adjacent regions, it will lead to an increase of about 19.44% and a decrease of 35.16% in outpatient visits in the particular region, respectively.
As for the PM2.5_max(−1), the direct effect value is 0.1984, the spatial spillover effect value is −0.3497, and the total effect value is −0.1513. The value of the spatial spillover effect is also larger than that of the direct effect. The conclusions further confirm that PM2.5 exposure has an important spatial spillover effect on residents’ health burden. It also points out that if the spatial models are not adopted, the direct impact of PM2.5 exposure may be overestimated and the spatial spillover effect may be underestimated.

4.3. Robustness Tests

4.3.1. Alternative Independent Variable Estimation

This study took the maximum value of PM2.5 concentrations as the core independent variable to analyze the spatial impact of PM2.5 exposure on residents’ health burden. To avoid the selection bias of the independent variable, the mean value of PM2.5 concentrations were selected as the substitution variable of PM2.5 exposure for the robustness test. The results were given in column (1) and column (2) of Table 9.
The results show that the direct coefficient and the spatial lag coefficient of PM2.5_avg are significantly positive at the 1% level when the spatial matrix is W1. As for the PM2.5_avg(−1), the direct coefficient and the spatial lag coefficient are also significantly positive at 5% level based on the spatial contiguity matrix W1. This suggests that the average of PM2.5 concentrations in a particular region or the adjacent regions has a positive impact on outpatient expense in the particular region, with spatial spillover effect and temporal lag effect. The results are consistent with the above conclusions, indicating that the results are stable and reliable. In other words, PM2.5 exposure has temporal lag effect and spatial spillover effect on residents’ health burden.

4.3.2. Alternative Dependent Variable Estimation

In this study, outpatient visits and outpatient expenses were selected as the substitution variables of health burden. To eliminate bias in the selection of indicators, this study used the number of hospitalizations (num_hos) as the dependent variable for robustness tests. The results were given in column (3) and column (4) of Table 9. The direct coefficients of PM2.5_max and PM2.5_max(−1) are all significant positive at 1% level when the spatial matrix is W1, while the spatial lag coefficients are significantly negative at 1% level. The results are consistent with the above conclusions and are robust, suggesting that PM2.5 exposure has spatial spillover effect and temporal lag effect on residents’ health burden.

4.3.3. Endogenous Test

Due to the two-way influence between the dependent variable and the independent variable or some important omitted variables, there may be endogenous problems between PM2.5 exposure and health burden. To reduce the possibility of estimation errors caused by endogenous problems, this study used the spatial Generalized Method of Moments (GMM) [60] to verify the reliability of the main empirical results. Kelejian et al. [68] proposed that W ( 1 λ W ) 1 X β was the relatively ideal instrumental variable, but the value of λ cannot be obtained in advance. Referring to the study of Yu and Liu [69], W*PM2.5_max and W*PM2.5_max(−1) were selected as the instrumental variables of the spatial GMM method, and Hansen J test was used to verify the rationality of the selected instrumental variables.
The results of the spatial GMM method were presented in column (5) and column (6) of Table 9. As can be seen from the results, the p values of Hansen J test are 0.2952 and 0.2276, respectively, indicating that W*PM2.5_max and W*PM2.5_max(−1) are suitable as instrumental variables. The coefficients are all significant positive at the 1% level, indicating that the results of the spatial GMM estimation are consistent with the previous results. Therefore, the results of this paper are robust and reliable.

5. Conclusions

Using the panel data of 29 Chinese regions from 2007 to 2017, this study used the spatial Durbin model (SDM) under space and time fixed effect to estimate the direct and the spatial lag effects of PM2.5 exposure on residents’ health burden in China based on three representative spatial weight matrices. The main conclusions drawn were as follows:
(1) Residents’ health burden and PM2.5 exposure are not randomly distributed among different regions in China, but have obvious spatial correlation and spatial clustering characteristics. There is a positive spatial correlation between PM2.5 concentrations and residents’ health burden. Also, the health burden in a particular region is not only affected by PM2.5 exposure in this region, but also influenced by the health burden in neighboring or economically similar regions.
(2) PM2.5 exposure has a significant positive impact on the health burden of residents in China. For example, PM2.5 pollution will increase outpatient expense and outpatient visits of residents. The possible reason is that exposure to air pollutants, such as PM2.5, makes the incidence of acute and chronic diseases such as respiratory diseases and cardiovascular diseases increase, and then the outpatient expense and outpatient visits will also increase.
(3) PM2.5 exposure has a spatial spillover effect on health burden. The increase of PM2.5 concentrations in surrounding regions or geographically close regions will lead to an increase of outpatient expense, but reduce outpatient visits in a particular region. Moreover, the value of the spatial spillover effect is larger than that of the direct effect, and the influence of PM2.5 exposure on health burden in China is mainly due to the spatial spillover effect.
(4) PM2.5 exposure has a long-term impact on residents’ health burden, that is to say, it has a temporal lag effect. Health burden in a particular region is not only affected by PM2.5 exposure in this region for a long time, but also affected by PM2.5 exposure in adjacent or geographically close regions for a long time. Long-term PM2.5 exposure can increase the incidence of chronic diseases such as cardiovascular disease, cancer, and diabetes.
According to the conclusions above, we have drawn some policy implications as follows:
(1) The impact of PM2.5 exposure on the health burden of residents in China must be given attention in the future. Frequent exposure to PM2.5 pollution can cause many acute and chronic diseases, which, eventually, brings great economic burden to residents. Hence, the government should strengthen environmental inspection and punishment, and shut down the enterprises with high pollution and low efficiency.
(2) Governments at all levels should break the administrative monopolies and achieve cross-regional cooperation in the field of environmental protection. The influence of PM2.5 exposure on health burden in China is mainly due to the spatial spillover effect, and this non-negligible fact urgently requires governments at all levels to strengthen the sense of cooperation, share information and technology of pollution control, and establish a common environmental protection system.
(3) Governments should establish a warning system and a long-term governance mechanism for environmental pollution. The impact of PM2.5 exposure on health burden is a long-term process, which can be easily overlooked. Therefore, the government should increase environmental protection publicity and improve residents’ awareness of environmental protection, so as to reduce residents health damage caused by environmental pollution.

Author Contributions

The research was designed and performed by M.Z. and W.Z. The data was collected and analyzed by M.Z., J.D. and W.Z. The paper was written by M.Z. and W.Z., and finally checked and revised by M.Z., J.D. and W.Z. All authors read and approved the final manuscript.

Funding

The authors would like to acknowledge the financial support from the Fundamental Research Funds for the Central Universities and National Social Science Foundation (14BJ176). Additionally, Weike Zhang gratefully acknowledges the financial support of the China Scholarship Council (Nos. 201806240231).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Average PM2.5 concentrations of 29 regions in China from 2007 to 2017 (μg/m3).
Table A1. Average PM2.5 concentrations of 29 regions in China from 2007 to 2017 (μg/m3).
Region20072008200920102011201220132014201520162017
Beijing50.8148.2147.3548.5943.0248.0641.5550.0044.0648.4645.00
Tianjin81.9374.8474.9578.1470.8171.8862.0381.7971.4774.8270.79
Hebei62.9560.3254.6957.1452.6453.1349.8461.0053.5255.3555.99
Shanxi33.9833.9326.7127.7626.5627.1124.8230.0825.3725.9224.97
Inner Mongolia11.0712.2811.9211.7712.0410.5210.1221.0510.8813.5212.49
Liaoning33.0234.4937.1938.3435.8533.1328.4935.8134.4747.6634.61
Jilin28.8229.7932.6534.6232.5929.4226.2433.2432.1247.5334.47
Heilongjiang19.2318.3620.3921.8321.5918.4216.8322.0922.4932.6826.69
Shanghai52.0756.9556.7858.0151.6149.9644.7054.1547.3161.0850.85
Jiangsu61.2361.7459.5860.0059.9758.0650.2060.6557.2265.3958.31
Zhejiang33.3837.8638.3534.2933.8131.6231.7034.9034.2733.2128.58
Anhui49.5758.1155.1952.2453.3849.6745.4653.1353.8157.0246.15
Fujian23.7324.6523.2321.7420.6819.9619.5620.3721.2919.9120.00
Jiangxi37.6341.0439.7437.4636.7233.6734.5634.9337.9934.8631.36
Shandong64.4469.3160.9558.2464.1257.3655.3564.7757.8161.6562.53
Henan60.3165.4450.6650.8754.5152.1048.7461.3351.5652.5648.91
Hubei45.8249.1846.8845.5849.4045.4740.3546.2948.1447.2937.68
Hunan41.6346.7945.0243.0540.5837.9939.4737.9340.8836.5531.43
Guangdong31.3334.2035.2834.3230.7429.1328.6028.9333.4926.7525.49
Guangxi35.2538.7638.2037.7133.9234.5136.1735.0836.9729.9528.67
Chongqing39.0136.1832.1332.3035.4330.3730.7730.9428.9825.9023.28
Sichuan37.1629.4829.7128.3934.6030.0029.6831.1128.5323.1422.85
Guizhou29.9329.1929.7129.9828.5528.8128.7726.4128.9323.1420.74
Yunnan16.2716.0916.3716.6116.5717.6215.7518.0717.2614.7714.25
Shaanxi32.3132.7425.6827.3628.2128.0726.3231.8225.7626.2824.23
Gansu21.3922.1119.3918.2218.4417.7116.3221.0518.3415.1515.38
Qinghai9.719.9310.049.1010.928.708.1610.659.706.947.85
Ningxia24.0220.8520.5119.9221.0317.3416.9121.9419.6017.3217.30
Xinjiang9.047.789.068.398.697.707.759.778.7210.4311.50
Table A2. Maximum PM2.5 concentrations of 29 regions in China from 2007 to 2017 (μg/m3).
Table A2. Maximum PM2.5 concentrations of 29 regions in China from 2007 to 2017 (μg/m3).
Region20072008200920102011201220132014201520162017
Beijing102.1090.8088.7089.7085.6092.0078.6097.0084.3093.5088.90
Tianjin96.1089.0088.5090.8083.8087.0075.0096.2086.1084.6084.10
Hebei90.2486.9180.8484.2678.5480.2974.9488.2581.1681.9784.65
Shanxi62.6360.2551.4953.2651.2951.7749.2655.4350.4550.3551.29
Inner Mongolia26.6827.5727.8727.5627.6025.0824.5330.0626.6232.6729.68
Liaoning48.0849.3452.6554.3651.1548.9542.8051.1849.3763.4850.55
Jilin38.9740.8342.8346.2743.7139.9036.4344.8443.8367.1947.13
Heilongjiang34.2433.1834.2238.8137.2933.0530.9238.7538.9556.5243.22
Shanghai60.1065.8064.4065.2056.6057.5052.7062.7057.9073.9059.50
Jiangsu67.9268.5865.7665.9266.4564.0555.9467.0463.9372.2564.88
Zhejiang51.6656.8557.2753.2951.4649.5350.0454.3552.7352.2547.15
Anhui59.1568.3164.8862.3963.0058.5154.2963.0865.0166.8855.85
Fujian47.6248.2446.4944.9343.8743.2242.3743.8845.0643.8743.70
Jiangxi50.3354.0152.2250.6149.7745.9547.2147.9151.8748.4444.33
Shandong80.2884.7475.9972.7479.4372.1071.0880.4573.6377.1279.54
Henan78.5683.9367.7468.1672.0469.8565.7080.0770.3470.0967.53
Hubei57.6361.8558.8057.4762.5857.4451.2259.2961.5860.6249.24
Hunan53.3958.9958.1555.7653.3450.0751.7650.8655.2949.6243.48
Guangdong41.0844.3444.8644.2340.2538.0437.7937.9643.8335.5634.21
Guangxi44.7448.7147.5447.4742.9443.3345.3244.4545.8938.8437.56
Chongqing72.5072.2060.5058.4063.6057.4057.7059.1055.6051.7044.70
Sichuan57.3247.4647.1945.5653.4847.5346.9149.9846.1939.7739.54
Guizhou40.9740.3241.1641.8140.1440.3640.7737.4240.6234.5231.49
Yunnan29.5329.7329.5129.9930.5131.4329.3232.1031.0528.1427.66
Shaanxi50.0250.3841.1242.9943.9644.2041.2049.6842.0942.6440.87
Gansu30.5932.3428.4626.2226.3125.8323.4930.0626.7622.3623.14
Qinghai17.7318.7518.7516.7319.2415.7915.2819.3117.6613.5415.66
Ningxia29.6026.0625.8624.8226.0622.0621.9827.3224.5822.6823.08
Xinjiang25.4422.8623.9123.3123.6621.2821.3625.5623.6726.5229.57

References

  1. Sarkodie, S.A.; Strezov, V.; Jiang, Y.J.; Evans, T. Proximate Determinants of Particulate Matter (PM2.5) Emission, Mortality and Life Expectancy in Europe, Central Asia, Australia, Canada and the US. Sci. Total Environ. 2019, 683, 489–497. [Google Scholar] [CrossRef] [PubMed]
  2. Cohen, G.; Gerber, Y. Air Pollution and Successful Aging: Recent Evidence and New Perspectives. Curr. Environ. Health Rep. 2017, 4, 1–11. [Google Scholar] [CrossRef] [PubMed]
  3. Apte, J.S.; Brauer, M.; Cohen, A.J.; Ezzati, M.; Pope, C.A. Ambient PM2.5 Reduces Global and Regional Life Expectancy. Environ. Sci. Technol. Lett. 2018, 5, 546–551. [Google Scholar] [CrossRef]
  4. Liu, X.; Xu, Y.J.; Ge, Y.; Zhang, W.K.; Herrera, F. A Group Decision Making Approach Considering Self-Confidence Behaviors and Its Application in Environmental Pollution Emergency Management. Int. J. Environ. Res. Public Health 2019, 16, 385. [Google Scholar] [CrossRef]
  5. Zhang, W.K.; Du, J.; Tian, X.L. Finding a Promising Venture Capital Project with TODIM under Probabilistic Hesitant Fuzzy Circumstance. Technol. Econ. Dev. Econ. 2018, 24, 2026–2044. [Google Scholar] [CrossRef]
  6. Li, Y.; Chiu, Y.H.; Wang, L.H.; Liu, Y.C.; Chiu, C.R. A Comparative Study of Different Energy Efficiency of OECD and Non-OECD Countries. Trop. Conserv. Sci. 2019, 12, 1940082919837441. [Google Scholar] [CrossRef]
  7. Gu, H.J.; Cao, Y.X.; Ehan, E.; Sunil, K.J. Human Health Damages Related to Air Pollution in China. Environ. Sci. Pollut. Res. 2019, 26, 13115–13125. [Google Scholar] [CrossRef]
  8. Lin, Y.L.; Zou, J.L.; Yang, W.; Li, C.Q. A Review of Recent Advances in Research on PM2.5 in China. Int. J. Environ. Res. Public Health 2018, 15, 438. [Google Scholar] [CrossRef]
  9. Munir, S.; Habeebullah, T.M. Vehicular Emissions on Main Roads in Makkah, Saudi Arabiaa Dispersion Modelling Study. Arabian J. Geosci. 2018, 11, 531. [Google Scholar] [CrossRef]
  10. Zhang, S.; Li, Y.; Hao, Y.; Zhang, Y.P. Does Public Opinion Affect Air Quality? Evidence based on the Monthly Data of 109 Prefecture-Level Cities in China. Energy Policy 2018, 116, 299–311. [Google Scholar] [CrossRef]
  11. Gao, J.H.; Woodward, A.; Vardoulakis, S. Haze, Public Health and Mitigation Measures in China: A Review of the Current Evidence for Further Policy Response. Sci. Total Environ. 2017, 578, 148–157. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, Y.; Li, J.; Zhu, G.; Yuan, Q. Spatio-Temporal Relationship and Evolvement of Socioeconomic Factors and PM2.5 in China During 1998–2016. Int. J. Environ. Res. Public Health 2019, 16, 1149. [Google Scholar] [CrossRef] [PubMed]
  13. Li, Y.; Zheng, C.; Ma, Z.; Quan, W. Acute and Cumulative Effects of Haze Fine Particles on Mortality and the Seasonal Characteristics in Beijing, China, 2005–2013: A Time-Stratified Case-Crossover Study. Int. J. Environ. Res. Public Health 2019, 16, 2383. [Google Scholar] [CrossRef] [PubMed]
  14. The State of China’s Ecological Environment in 2018. Available online: http://www.mee.gov.cn/home/jrtt1/201905/t20190529_704841.shtml (accessed on 17 August 2019).
  15. Li, Y.; Chiu, Y.H.; Lu, L.C. New Energy Development and Pollution Emissions in China. Int. J. Environ. Res. Public Health 2019, 16, 1764. [Google Scholar] [CrossRef]
  16. Chen, H.M.; Liu, J.; Li, Y.; Chiu, Y.H.; Lin, T.Y. A Two-stage Dynamic Undesirable Data Envelopment Analysis Model Focused on Media Reports and the Impact on Energy and Health Efficiency. Int. J. Environ. Res. Public Health 2019, 16, 1535. [Google Scholar] [CrossRef]
  17. Li, Y.; Chiu, Y.H.; Lu, L.C. Energy, CO2, AQI and Economic Performance in 31 Cities in China: A Slacks-Based Dynamic Data Envelopment Analysis. Carbon Manag. 2019, 10, 269–286. [Google Scholar] [CrossRef]
  18. Li, Y.; Chiu, Y.H.; Lin, T.Y. Research on New and Traditional Energy Sources in OECD Countries. Int. J. Environ. Res. Public Health 2019, 16, 1122. [Google Scholar] [CrossRef]
  19. National Bureau of Statistics of China. Available online: http://data.stats.gov.cn/easyquery.htm?cn = C01 (accessed on 19 August 2019).
  20. Lelieveld, J.; Evans, J.S.; Fnais, M.; Giannadaki, D.; Pozzer, A. The Contribution of Outdoor Air Pollution Sources to Premature Mortality on a Global Scale. Nature 2015, 525, 367–371. [Google Scholar] [CrossRef]
  21. Wu, Y.; Lu, B.; Zhu, X.L. Seasonal Variations, Source Apportionment, and Health Risk Assessment of Heavy Metals in PM2.5 in Ningbo, China. Aerosol Air Qual. Res. 2019, 19, 2083–2092. [Google Scholar] [CrossRef]
  22. Polezer, G.; Tadano, Y.S.; Siqueira, H.V. Assessing the Impact of PM2.5 on Respiratory Disease using Artificial Neural Networks. Environ. Pollut. 2018, 235, 394–403. [Google Scholar] [CrossRef]
  23. Andrea, R.V.L.; Yezid, R.R.N.; Camilo, B.B.L. Short-Term Effects of Air Pollution on Respiratory and Circulatory Morbidity in Colombia 2011–2014: A Multi-City, Time-Series Analysis. Int. J. Environ. Res. Public Health 2018, 15, 1610. [Google Scholar]
  24. Weber, S.A.; Insaf, T.Z.; Hall, E.S. Assessing the Impact of Fine Particulate Matter (PM2.5) on Respiratory-Cardiovascular Chronic Diseases in the New York City Metropolitan Area Using Hierarchical Bayesian Model Estimates. Environ. Res. 2016, 151, 399–409. [Google Scholar] [CrossRef] [PubMed]
  25. Cao, Q.L.; Rui, G.Q.; Liang, Y. Study on PM2.5 Pollution and the Mortality Due to Lung Cancer in China Based on Geographic Weighted Regression Model. BMC Public Health 2018, 18, 925. [Google Scholar] [CrossRef] [PubMed]
  26. De, M.A.; Amoatey, P.; Khaniabadi, Y.O. Mortality and Morbidity for Cardiopulmonary Diseases Attributed to PM2.5 Exposure in the Metropolis of Rome, Italy. Eur. J. Intern. Med. 2018, 57, 49–57. [Google Scholar]
  27. Badyda, A.J.; Grellier, J.D.P. Ambient PM2.5 Exposure and Mortality Due to Lung Cancer and Cardiopulmonary Diseases in Polish Cities. Adv. Exp. Med. Biol. 2017, 944, 9–17. [Google Scholar]
  28. Dabass, A.; Talbott, E.O.; Venkat, A.; Rager, J.; Marsh, G.M. Association of Exposure to Particulate Matter (PM2.5) Air Pollution and Biomarkers of Cardiovascular Disease Risk in Adult NHANES participants (2001–2008). Int. J. Hyg. Environ. Health 2016, 219, 301–310. [Google Scholar] [CrossRef]
  29. Chen, X.Z.; Zhang, W.K.; Yang, K.; Wang, L.L.; Liu, J.; Wang, L.; Hu, J.K.; Zhang, B.; Chen, Z.X.; Chen, J.P.; et al. Correlation between Serum CA724 and Gastric Cancer: Multiple Analyses based on Chinese Population. Mol. Biol. Rep. 2012, 39, 9031–9039. [Google Scholar] [CrossRef]
  30. De Grubb, M.C.M.; Kilbourne, B.; Kilbourne, K. Socioeconomic, Environmental, and Geographic Factors and US Lung Cancer Mortality, 1999–2009. Fam. Med. Community Health 2017, 5, 3–12. [Google Scholar] [CrossRef]
  31. Chen, J.C.; Wang, X.H.; Wellenius, G.A. Ambient Air Pollution and Neurotoxicity on Brain Structure: Evidence from Women’s Health Initiative Memory Study. Ann. Neurol. 2015, 78, 466–476. [Google Scholar] [CrossRef]
  32. Calderon, G.L.; Calderon, G.A.; Torres, J.R. Air Pollution and Your Brain: What do You Need to Know Right Now. Prim. Health Care Res. Dev. 2015, 16, 329–345. [Google Scholar] [CrossRef]
  33. Shin, J.; Han, S.H.; Choi, J. Exposure to Ambient Air Pollution and Cognitive Impairment in Community-Dwelling Older Adults: The Korean Frailty and Aging Cohort Study. Int. J. Environ. Res. Public Health 2019, 16, 3767. [Google Scholar] [CrossRef] [PubMed]
  34. Li, R.L.; Ho, Y.C.; Luo, C.W.; Lee, S.S.; Kuan, Y.H. Influence of PM2.5 Exposure Level on the Association between Alzheimer’s Disease and Allergic Rhinitis: A National Population-Based Cohort Study. Int. J. Environ. Res. Public Health 2019, 16, 3357. [Google Scholar] [CrossRef] [PubMed]
  35. Tran, B.L.; Chang, C.C.; Hsu, C.S.; Chen, C.C.; Tseng, W.C.; Hsu, S.H. Threshold Effects of PM2.5 Exposure on Particle-Related Mortality in China. Int. J. Environ. Res. Public Health 2019, 16, 3549. [Google Scholar] [CrossRef] [PubMed]
  36. Xia, Y.; Guan, D.B.; Jiang, X.J. Assessment of Socioeconomic Costs to China’s Air Pollution. Atmos. Environ. 2016, 139, 147–156. [Google Scholar] [CrossRef]
  37. Yang, Y.; Luo, L.W.; Song, C. Spatiotemporal Assessment of PM2.5-Related Economic Losses from Health Impacts during 2014–2016 in China. Int. J. Environ. Res. Public Health 2018, 15, 1278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Miao, W.J.; Huang, X.; Song, Y. An Economic Assessment of the Health Effects and Crop Yield Losses Caused by Air Pollution in Mainland China. J. Environ. Sci. 2017, 56, 102–113. [Google Scholar] [CrossRef]
  39. Liu, Y.W.; Xie, S.G.; Yu, Q. Short-Term Effects of Ambient Air Pollution on Pediatric Outpatient Visits for Respiratory Diseases in Yichang City, China. Environ. Pollut. 2017, 227, 116–124. [Google Scholar] [CrossRef]
  40. Zeng, X.G.; Ruan, F.F.; Peng, Y.Y. Health Effects’ Spatial Distribution Analysis of PM2.5 Pollution in China Based on Spatial Grid Scale. China Environ. Sci. 2019, 39, 2624–2632. [Google Scholar]
  41. Sun, W.; Sun, J.Y. Daily PM2.5 Concentration Prediction based on Principal Component Analysis and LSSVM Optimized by Cuckoo Search Algorithm. J. Environ. Manag. 2017, 188, 144–152. [Google Scholar] [CrossRef]
  42. Tobler, W.R. A Computer Movie Simulating Urban Growth in the Detroit Region. Econ. Geogr. 1970, 46, 234–240. [Google Scholar] [CrossRef]
  43. Anselin, L.; Griffith, D.A. Do Spatial Effects Really Matter in Regression Analysis? Pap. Reg. Sci. 1988, 65, 1–34. [Google Scholar]
  44. Chen, X.Y.; Shao, S.; Tian, Z.H.; Xie, Z.; Yin, P. Impacts of Air Pollution and its Spatial Spillover Effect on Public Health Bbased on China’s Big Data Sample. J. Clean. Prod. 2017, 142, 915–925. [Google Scholar] [CrossRef]
  45. Shao, S.; Liu, A.J.; Yu, C.J.; Yang, H.J.; Li, H.; Li, B. Spatial Time-Frequency Distribution of Cross Term-Based Direction-of-Arrival Estimation for Weak Non-Stationary Signal. Eurasip J. Wirel. Commun. Netw. 2019, 2019, 1–12. [Google Scholar] [CrossRef]
  46. Cheng, Z.H.; Li, L.S.; Liu, J. Identifying the Spatial Effects and Driving Factors of Urban PM2.5 Pollution in China. Ecol. Indic. 2017, 82, 61–75. [Google Scholar] [CrossRef]
  47. Feng, Y.; Cheng, J.H.; Shen, J. Spatial Effects of Air Pollution on Public Health in China. Environ. Resour. Econ. 2019, 73, 229–250. [Google Scholar] [CrossRef]
  48. Li, G.Q.; He, R.W. PM2.5 Pollution and Health Spending: Temporal Lag Effects and Spatial Spillover Effects. J. Saf. Environ. 2019, 19, 326–336. (In Chinese) [Google Scholar]
  49. Columbia University Center for Socio-Economic Data and Applications. Available online: https://sedac.ciesin.columbia.edu/ (accessed on 20 July 2019).
  50. Chinese Health Statistics Yearbook 2008–2017. Available online: http://data.cnki.net/yearbook/Single/N2019030282/ (accessed on 23 July 2019).
  51. Chinese Health Statistics Yearbook 2008–2017. Available online: http://www.stats.gov.cn/tjsj/ndsj/ (accessed on 23 July 2019).
  52. Xie, Y.B.; Chen, J.; Li, W. An Assessment of PM2.5 Related Health Risks and Impaired Values of Beijing Residents in a Consecutive High-Level Exposure During Heavy Haze Days. Environ. Sci. 2014, 35, 1–8. [Google Scholar]
  53. Mu, Q.; Zhang, S.Q. An Evaluation of the Economic Loss Due to the Heavy Haze during January 2013 in China. China Environ. Sci. 2013, 33, 2087–2094. (In Chinese) [Google Scholar]
  54. Wei, G.R.; Shi, X.M. Evaluation the Extent of Health Damage Caused by PM2.5 Particulate in Xi’an City. Environ. Sci. 2018, 39, 3014–3021. (In Chinese) [Google Scholar]
  55. Yang, D.R.; Wang, H.X.; Lan, L.; Tao, Y.; Mi, S.Q. Evaluation of Economic Losses Due to Human Health Damage from Air Pollution before and after Energy Reform in Lanzhou. Environ. Eng. 2013, 31, 112–116. (In Chinese) [Google Scholar]
  56. Dong, Y.; Xu, G.Z.; Wang, A.H.; Lu, B.B. Evaluation of PM2.5 Air Pollution Impacts on the Public Health and Related Economic Loss in a City from 2014–2016. China Prev. Med. 2018, 19, 579–582. [Google Scholar]
  57. Cao, C.H.; Han, L.Y. The Assessment on the Social Health Costs Caused by Fog and Haze. Stat. Res. 2015, 32, 19–23. [Google Scholar]
  58. Shao, S.; Li, X.; Cao, J.H.; Yang, L.L. China’s Economic Policy Choices for Governing Smog Pollution Based on Spatial Spillover Effects. Econ. Res. 2016, 51, 73–88. [Google Scholar]
  59. Moran, P. Notes on Continuous Stochastic Phenomena. Biometrika 1950, 37, 17–23. [Google Scholar] [CrossRef] [PubMed]
  60. Elhorst, J.P. Spatial Econometrics: From Cross-Sectional Data to Spatial Panels; Springer: Heidelberg, Germany, 2014. [Google Scholar]
  61. Anselin, L.; Arribas, B.D. Spatial Fixed Effects and Spatial Dependence in A Single Cross-Section. Pap. Reg. Sci. 2013, 37, 17–23. [Google Scholar] [CrossRef]
  62. Belotti, F.; Hughes, G.; Mortari, A.P. Spatial Panel-Data Models Using Stata. Stata J. 2017, 17, 139–180. [Google Scholar] [CrossRef] [Green Version]
  63. Qiu, Y.L.; Wang, Z.Q. The Forty Years Development of Medical Insurance System in China: Process, Experience and Prospect. J. Cent. China Norm. Univ. (Human. Soc. Sci.) 2019, 58, 23–30. (In Chinese) [Google Scholar]
  64. Getis, A.; Ord, J. The Analysis of Spatial Association by Use of Distance Statistics. Geogr. Anal. 1992, 24, 189–206. [Google Scholar] [CrossRef]
  65. Gan, W.Q.; Davies, H.W.; Koehoorn, M.; Brauer, M. Association of Long-Term Exposure to Community Noise and Traffic-Related Air Pollution with Coronary Heart Disease Mortality. Am. J. Epidemiol. 2012, 175, 898–906. [Google Scholar] [CrossRef] [Green Version]
  66. Heinrich, J.; Thiering, E.; Rzehak, P.; Krämer, U.; Hochadel, M.; Rauchfuss, K.M.; Gehring, U.; Wichmann, H. Long-Term Exposure to NO2 and PM10 and All-Cause and Cause−Specific Mortality in a Prospective Cohort of Women. Occup. Environ. Med. 2013, 70, 179–186. [Google Scholar] [CrossRef] [Green Version]
  67. Lanzinger, S.; Rosenbauer, J.; Sugiri, D. Impact of Long–Term Air Pollution Exposure on Metabolic Control in Children and Adolescents with Type 1 Diabetes: Results from the DPV Registry. Diabetologia 2018, 61, 1354–1361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Kelejian, H.H.; Prucha, I.R.; Yuzefovich, Y. Instrumental Variable Estimation of a Spatial Autoregressive Model with Autoregressive Disturbances: Large and Small Sample Results. Spat. Spatiotemporal Econom. 2004, 18, 163–198. [Google Scholar]
  69. Yu, Y.Z.; Liu, D.Y. The Effect of the Space Outflow of China’s Regional Innovation and the Effect of the Outflow of Value Chanins: A Study, from the Perspective of the Innovative Value Chain, on the Model of the Panel of Multidimentsional Space. Manag. World 2013, 7, 6–20, 70, 187. (In Chinese) [Google Scholar]
Figure 1. The spatial distribution of the core variables in 2008 and 2017.
Figure 1. The spatial distribution of the core variables in 2008 and 2017.
Ijerph 16 04695 g001
Figure 2. Local Moran’s I scatter plot in 2008 and 2017. Note: Numbers 1 to 29 represent Beijing, Tianjin, Hebei, Shanxi, Inner Mongolia, Liaoning, Jilin, Heilongjiang, Shanghai, Jiangsu, Zhejiang, Anhui, Fujian, Jiangxi, Shandong, Henan, Hubei, Hunan, Guangdong, Guangxi, Chongqing, Sichuan, Guizhou, Yunnan, Shaanxi, Gansu, Qinghai, Ningxia and Xinjiang, respectively.
Figure 2. Local Moran’s I scatter plot in 2008 and 2017. Note: Numbers 1 to 29 represent Beijing, Tianjin, Hebei, Shanxi, Inner Mongolia, Liaoning, Jilin, Heilongjiang, Shanghai, Jiangsu, Zhejiang, Anhui, Fujian, Jiangxi, Shandong, Henan, Hubei, Hunan, Guangdong, Guangxi, Chongqing, Sichuan, Guizhou, Yunnan, Shaanxi, Gansu, Qinghai, Ningxia and Xinjiang, respectively.
Ijerph 16 04695 g002
Figure 3. Hot spot analysis results in 2008 and 2017.
Figure 3. Hot spot analysis results in 2008 and 2017.
Ijerph 16 04695 g003
Table 1. Description of the variable.
Table 1. Description of the variable.
TypeVariableSymbolDefinition
Dependent variableOutpatient expenseexp_outThe ratio of the total outpatient expense to the total number of outpatient visits in the form of the natural logarithm
Outpatient visitsnum_outThe ratio of the total number of outpatient visits to the total population in the form of the natural logarithm
The number of hospitalizationnum_hosThe ratio of the total number of hospitalization to the total population
Independent variableMaximum PM2.5 concentrationsPM2.5_maxThe maximum values of PM2.5 concentrations in the form of natural logarithm
Maximum PM2.5 concentrations lag by one stagePM2.5_max(−1)The maximum values of the last year’s PM2.5 concentrations in the form of the natural logarithm
Average PM2.5 concentrationsPM2.5_avgThe average values of PM2.5 concentrations in the form of the natural logarithm
Average PM2.5 concentrations lag by one stagePM2.5_avg(−1)The average values of the last year’s PM2.5 concentrations in the form of natural logarithm
Control variablePer capita GDPPGDPThe ratio of gross domestic product to the total population in the form of the natural logarithm
The ratio of urban populationurbanThe ratio of the urban population to the total population
The number of medical institutionsnum_instThe ratio of the total number of medical institutions to the total population in the form of the natural logarithm
The number of hospital bedsnum_bedThe ratio of the total number of hospital beds to the total population in the form of the natural logarithm
The number of doctorsnum_doctorThe ratio of the total number of doctors to the total population in the form of the natural logarithm
Table 2. Descriptive statistics.
Table 2. Descriptive statistics.
VariableObsMeanS.D.MinMedianMax
exp_out2905.2300.2974.3855.2426.248
num_out2901.5480.3240.8321.5012.397
num_hos2900.1250.0430.0390.1260.224
PM2.5_max2903.8470.4132.6053.9034.575
PM2.5_max(−1)2903.8410.4142.6053.8974.575
PM2.5_avg2903.4460.5341.9383.5194.404
PM2.5_avg(−1)2903.4220.5491.9383.4884.404
PGDP2901.3680.514−0.0101.3692.557
urban2900.5480.1340.2910.5300.896
num_inst2901.7860.5100.2081.9492.455
num_bed2903.7730.2363.1403.8024.227
num_doctor2904.2930.1983.6894.3104.978
Table 3. The selection results of spatial autoregression model (SAR), spatial errors model (SEM), and spatial Durbin model (SDM).
Table 3. The selection results of spatial autoregression model (SAR), spatial errors model (SEM), and spatial Durbin model (SDM).
NameModelSelection CriteriaChi-Square Valuep-Value
SAR y = ρ W y + X β + ε λ = 0 32.320.0000
SEM y = X β + u , u = λ W u + ε λ = ρ β 31.370.0000
SDM y = ρ W y + X β + λ W X + ε λ 0 & λ ρ β
Hausman test: The Chi-square value is 11.35, and the p-value is 0.0782.
Table 4. Global Moran’s I values of exp_out and PM2.5_max (2008–2017).
Table 4. Global Moran’s I values of exp_out and PM2.5_max (2008–2017).
Yearexp_outPM2.5_max
W1W2W3W1W2W3
20080.201 **0.058 **0.079 *0.527 ***0.238 ***0.097 *
20090.278 ***0.167 ***0.307 ***0.519 ***0.243 ***0.091 *
20100.270 ***0.168 ***0.296 ***0.514 ***0.233 ***0.080
20110.256 ***0.163 ***0.339 ***0.504 ***0.236 ***0.067
20120.227 ***0.149 ***0.310 ***0.511 ***0.225 ***0.038
20130.215 ***0.136 ***0.267 ***0.512 ***0.256 ***0.062
20140.197 ***0.122 ***0.254 ***0.542 ***0.242 ***0.064
20150.170 ***0.101 ***0.234 ***0.525 ***0.261 ***0.112 *
20160.167 ***0.102 ***0.245 ***0.545 ***0.289 ***0.081
20170.163 **0.091 ***0.245 ***0.468 ***0.239 ***0.099 *
Note: ***, **, and * denote significance at the 1%, 5%, and 10% levels, respectively; The Moran’s I values of other variables are not reported for space limitation.
Table 5. Estimation results of the impact of PM2.5 exposure on outpatient expense.
Table 5. Estimation results of the impact of PM2.5 exposure on outpatient expense.
VariableSpatial Contiguity Matrix W1Spatial Distance Matrix W2Spatial Economy Matrix W3
(1)(2)(3)(4)(5)(6)
PM2.5_max0.1017 ***0.1282 ***0.1773 ***
(2.79)(4.41)(8.58)
PM2.5_max(−1)0.0971 ***0.1186 ***0.1713 ***
(2.65)(4.16)(8.34)
PGDP−0.2317 ***−0.2323 ***−0.3029 ***−0.3097 ***−0.1916 ***−0.1975 ***
(−4.76)(−4.76)(−5.39)(−5.51)(−2.70)(−2.77)
urban1.3876 ***1.3840 ***1.7333 ***1.7686 ***1.2335 ***1.2656 ***
(7.32)(7.30)(10.04)(10.22)(6.86)(7.03)
num_inst−0.0710 **−0.0717 **−0.0469 *−0.0449 *−0.0226−0.0191
(−2.47)(−2.49)(−1.76)(−1.68)(−0.78)(−0.66)
num_bed0.4853 ***0.4929 ***0.3549 ***0.3659***0.3868 ***0.3915 ***
(5.76)(5.85)(5.29)(5.47)(5.57)(5.58)
num_doctor0.01780.01180.10520.0887−0.0487−0.0556
(0.23)(0.15)(1.26)(1.06)(−0.57)(−0.64)
W*PM2.5_max0.1531 ***0.4066 **0.0106
(2.76)(2.12)(0.17)
W*PM2.5_max(−1)0.1551 ***0.4771 **0.0012
(2.78)(2.55)(0.02)
W*PGDP−0.4213 ***−0.4178 ***−1.4444 ***−1.4888 ***−0.0383−0.0382
(−4.62)(−4.58)(−4.70)(−4.84)(−0.31)(−0.31)
W*urban2.7245 ***2.7882 ***5.4106 ***5.6841 ***0.79620.8438 *
(6.62)(6.77)(4.43)(4.66)(1.56)(1.65)
W*num_inst0.1294 *0.1458 **−0.3039−0.2645−0.3431 ***−0.3369 ***
(1.87)(2.08)(−1.23)(−1.07)(−4.05)(−3.95)
W*num_bed0.3319 *0.3426 **0.22440.3220−0.0646−0.0803
(1.91)(1.96)(0.47)(0.68)(−0.30)(−0.37)
W*num_doctor−0.7153 ***−0.7542 ***0.0433−0.0844−0.0392−0.0449
(−3.55)(−3.71)(0.08)(−0.15)(−0.15)(−0.17)
ρ−0.1312−0.1230−0.7074 ***−0.7217 ***−0.2088 *−0.2034 *
−1.49)(−1.40)(−2.95)(−3.01)(−1.78)(−1.73)
sigma2_e0.0122 ***0.0123 ***0.0118 ***0.0118 ***0.0123 ***0.0124 ***
(11.94)(11.95)(12.03)(12.03)(12.30)(12.29)
N290290290290290290
Notes: ***, **, and * represent significance at the 1%, 5% and 10% levels, respectively; The numbers in brackets are t statistic values.
Table 6. The direct effects, the spatial spillover effects and the total effects of SDM (the dependent variable is exp_out).
Table 6. The direct effects, the spatial spillover effects and the total effects of SDM (the dependent variable is exp_out).
Type VariableCoefficientt−Valuep−Value
Direct effectsPM2.5_max0.0987 **2.550.011
PM2.5_max(−1)0.0942 **2.420.015
Spatial Spillover EffectsPM2.5_max0.1245 **2.370.018
PM2.5_max(−1)0.1283 **2.420.016
Total EffectsPM2.5_max0.2232 ***7.110.000
PM2.5_max(−1)0.2225 ***7.020.000
Note: *** and ** represent significance at the 1% and 5% levels, respectively.
Table 7. The results of the spatial impact of PM2.5 exposure on outpatient visits.
Table 7. The results of the spatial impact of PM2.5 exposure on outpatient visits.
VariableSpatial Contiguity Matrix W1Spatial Distance Matrix W2Spatial Economy Matrix W3
(1)(2)(3)(4)(5)(6)
PM2.5_max0.2114 ***0.3311 ***0.0070
(4.10)(8.56)(0.24)
PM2.5_max(−1)0.2154 ***0.3169 ***0.0178
(4.16)(8.22)(0.61)
PGDP0.1840 ***0.1874 ***0.3900 ***0.3944 ***0.3556 ***0.3569 ***
(2.74)(2.79)(5.25)(5.23)(3.50)(3.52)
urban0.6517 **0.6314 **−0.4861 **−0.5139 **0.23400.2304
(2.54)(2.46)(−2.13)(−2.22)(0.90)(0.90)
num_inst−0.0963 **−0.1007 **−0.2164 ***−0.2206 ***−0.1941 ***−0.1933 ***
(−2.40)(−2.52)(−6.12)(−6.16)(−4.68)(−4.68)
num_bed−0.7171 ***−0.7150 ***−0.6721 ***−0.6841 ***−0.4667 ***−0.4541 ***
(−6.10)(−6.09)(−7.50)(−7.54)(−4.70)(−4.57)
num_doctor0.8027 ***0.8038 ***0.9740 ***0.9828 ***0.4070 ***0.3994 ***
(7.65)(7.65)(8.75)(8.71)(3.34)(3.28)
W*PM2.5_max−0.3241 ***−2.3216 ***−0.2375 ***
(−4.33)(−9.33)(−2.66)
W*PM2.5_max(−1)−0.3232 ***−2.1647 ***−0.2460 ***
(−4.29)(−8.76)(−2.77)
W*PGDP−0.1371−0.14160.62400.6221−0.3989 **−0.3803 **
(−1.11)(−1.15)(1.53)(1.50)(−2.21)(−2.11)
W*urban−1.7539 ***−1.7123 ***−2.6430−2.9903 *1.9654 ***1.8389 ***
(−3.23)(−3.16)(−1.62)(−1.83)(2.78)(2.61)
W*num_inst−0.5490 ***−0.5367 ***−0.8454 **−0.8419 **−0.3412 ***−0.3607 ***
(−5.05)(−4.90)(−2.47)(−2.42)(−2.83)(−2.99)
W*num_bed−0.1504−0.1298−3.0001 ***−2.7721 ***0.35990.3428
(−0.59)(−0.51)(−4.66)(−4.27)(1.17)(1.12)
W*num_doctor0.6217 **0.5846 *2.2185 ***2.1068 ***−1.2269 ***−1.2073 ***
(1.96)(1.83)(2.89)(2.69)(−3.29)(−3.26)
ρ0.2720 ***0.2758 ***0.10900.1156−0.2840 **−0.2862 **
(3.51)(3.56)(0.56)(0.59)(−2.36)(−2.38)
sigma2_e0.0238 ***0.0238 ***0.0208 ***0.0213 ***0.0252 ***0.0251 ***
(11.95)(11.93)(12.08)(12.09)(12.13)(12.14)
N290290290290290290
Note: ***, **, and * represent significance at the 1%, 5% and 10% levels, respectively; The numbers in brackets are t statistic values.
Table 8. The direct effects, the spatial spillover effects and the total effects of SDM (the dependent variable is num_out).
Table 8. The direct effects, the spatial spillover effects and the total effects of SDM (the dependent variable is num_out).
Type VariableCoefficientt-Valuep-Value
Direct EffectsPM2.5_max0.1944 ***3.920.000
PM2.5_max(−1)0.1984 ***3.990.000
Spatial Spillover EffectsPM2.5_max−0.3516 ***−4.090.000
PM2.5_max(−1)−0.3497 ***−4.030.000
Total EffectsPM2.5_max−0.1572 **−2.260.024
PM2.5_max(−1)−0.1513 **−2.160.031
Note: *** and ** represent significance at the 1% and 5% level, respectively.
Table 9. Results of the robustness tests.
Table 9. Results of the robustness tests.
Variableexp_outexp_outnum_hosnum_hosGMMGMM
(1)(2)(3)(4)(5)(6)
PM2.5_max0.0335 ***
(8.60)
PM2.5_max(−1)0.0329 ***
(8.33)
PM2.5_avg0.0603 ***
(3.16)
PM2.5_avg(−1)0.0607 ***
(3.35)
PGDP−0.3068 ***−0.3404 ***0.0141 ***0.0143 ***1.3684 ***1.3684 ***
(−5.21)(−5.81)(2.73)(2.76)(45.35)(45.35)
urban1.7402 ***1.8307 ***−0.0923 ***−0.0956 ***0.5480 ***0.5480 ***
(9.64)(10.21)(−4.61)(−4.74)(69.55)(69.55)
num_inst−0.0421−0.0398−0.0112 ***−0.0119 ***1.7864 ***1.7864 ***
(−1.48)(−1.42)(−3.66)(−3.86)(59.62)(59.62)
num_bed0.3559 ***0.3656 ***0.0985 ***0.0997 ***3.7732 ***3.7732 ***
(5.36)(5.23)(11.01)(11.07)(271.75)(271.75)
num_doctor0.13630.1422 *0.00510.00534.2932 ***4.2932 ***
(1.55)(1.64)(0.63)(0.65)(370.05)(370.05)
W*PM2.5_max−0.0484 ***3.8471 ***
(−8.44)(158.76)
W*PM2.5_max(−1)−0.0472 ***3.8407 ***
(−8.11)(157.95)
W*PM2.5_avg0.2953 **
(2.18)
W*PM2.5_avg(−1)0.4607 ***
(3.31)
W*PGDP−1.3705 ***−1.6526 ***−0.0333 ***−0.0341 ***
(−4.18)(−4.94)(−3.49)(−3.55)
W*urban5.3744 ***5.8229 ***0.1284 ***0.1368 ***
(4.27)(4.68)(3.09)(3.28)
W*num_inst−0.3455−0.3123−0.00040.0010
(−1.34)(−1.23)(−0.05)(0.12)
W*num_bed−0.3486−0.0944−0.0570 ***−0.0556 ***
(−0.75)(−0.20)(−3.04)(−2.92)
W*num_doctor0.48730.5663−0.0430 **−0.0472 **
(0.83)(0.97)(−2.00)(−2.17)
ρ−0.5548 **−0.6443 ***0.6150 ***0.6131 ***
(−2.38)(−2.72)(11.86)(11.75)
sigma2_e0.0132 ***0.0128 ***0.0001 ***0.0001 ***
(12.04)(12.03)(11.57)(11.58)
N290290290290290290
Notes: ***, **, and * represent significance at the 1%, 5% and 10% level, respectively; The numbers in brackets are t statistic values.

Share and Cite

MDPI and ACS Style

Zeng, M.; Du, J.; Zhang, W. Spatial-Temporal Effects of PM2.5 on Health Burden: Evidence from China. Int. J. Environ. Res. Public Health 2019, 16, 4695. https://doi.org/10.3390/ijerph16234695

AMA Style

Zeng M, Du J, Zhang W. Spatial-Temporal Effects of PM2.5 on Health Burden: Evidence from China. International Journal of Environmental Research and Public Health. 2019; 16(23):4695. https://doi.org/10.3390/ijerph16234695

Chicago/Turabian Style

Zeng, Ming, Jiang Du, and Weike Zhang. 2019. "Spatial-Temporal Effects of PM2.5 on Health Burden: Evidence from China" International Journal of Environmental Research and Public Health 16, no. 23: 4695. https://doi.org/10.3390/ijerph16234695

APA Style

Zeng, M., Du, J., & Zhang, W. (2019). Spatial-Temporal Effects of PM2.5 on Health Burden: Evidence from China. International Journal of Environmental Research and Public Health, 16(23), 4695. https://doi.org/10.3390/ijerph16234695

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop