Next Article in Journal
Evaluating Park Use and Satisfaction: The Case of Trojan Park in St. Louis Missouri
Previous Article in Journal
Comment on Choi, Y.-J.; Choi, H.-B.; O’Donnell, M. Disaster Reintegration Model: A Qualitative Analysis on Developing Korean Disaster Mental Health Support Model. Int. J. Environ. Res. Public Health 2018, 15, 362
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Health Risk and Resilience Assessment with Respect to the Main Air Pollutants in Sichuan

1
School of Civil Engineering and Architecture, Southwest Petroleum University, Chengdu 610500, China
2
State Key Laboratory of Resources and Environmental Information System, Institute of Geographic Sciences and Natural Resources Research, CAS, Beijing 100101, China
3
Synthesis Research Centre of Chinese Ecosystem Research Network, Key Laboratory of Ecosystem Network Observation and Modelling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
4
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(15), 2796; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152796
Submission received: 25 June 2019 / Revised: 31 July 2019 / Accepted: 1 August 2019 / Published: 6 August 2019

Abstract

:
Rapid urbanization and industrialization in developing countries have caused an increase in air pollutant concentrations, and this has attracted public concern due to the resulting harmful effects to health. Here we present, through the spatial-temporal characteristics of six criteria air pollutants (PM2.5, PM10, SO2, NO2, CO, and O3) in Sichuan, a human health risk assessment framework conducted to evaluate the health risk of different age groups caused by ambient air pollutants. Public health resilience was evaluated with respect to the risk resulting from ambient air pollutants, and a spatial inequality analysis between the risk caused by ambient air pollutants and hospital density in Sichuan was performed based on the Lorenz curve and Gini coefficient. The results indicated that high concentrations of PM2.5 (47.7 μg m−3) and PM10 (75.9 μg m−3) were observed in the Sichuan Basin; these two air pollutants posed a high risk to infants. The high risk caused by PM2.5 was mainly distributed in Sichuan Basin (1.14) and that caused by PM10 was principally distributed in Zigong (1.01). Additionally, the infants in Aba and Ganzi had high health resilience to the risk caused by PM2.5 (3.89 and 4.79, respectively) and PM10 (3.28 and 2.77, respectively), which was explained by the low risk in these two regions. These regions and Sichuan had severe spatial inequality between the infant hazard quotient caused by PM2.5 (G = 0.518, G = 0.493, and G = 0.456, respectively) and hospital density. This spatial inequality was also caused by PM10 (G = 0.525, G = 0.526, and G = 0.466, respectively), which is mainly attributed to the imbalance between hospital distribution and risk caused by PM2.5 (PM10) in these two areas. Such research could provide a basis for the formulation of medical construction and future air pollution control measures in Sichuan.

1. Introduction

Air pollution is principally formed by a multifaceted mix of gaseous pollutants and suspended particulates [1], which affects atmospheric emission, absorption, and scattering of light. Recent research has reported that 3.7 million deaths due to exposure to outdoor air pollution in 2012. The number of deaths caused by ambient air pollutants is continuously rising, and will become the main environmental issue causing premature death by 2050 [2]. The continuous increase in human activities, such as urbanization and modernization, industrialization, human population, traffic density, and the number of motor vehicles, are the major contributors to the rising ambient air quality problems [3]. Air pollution has become a new threat to air quality [4], regional sustainable development, global climate change [5,6], and ecological suitability for human settlement [7]. In particular, air pollution continues to be a global environmental problem, and has been recognized as a severe public health risk [8]. Simultaneously, this health risk has drawn media coverage and captured public concern, motivating the government to take aggressive actions designed to control air pollution.
The sensitive parts of the lungs are susceptible to invasion by small, fine, and ultrafine ambient particles; these particles can exacerbate respiratory diseases such as bronchitis and emphysema, as well as existing heart disease, thus further leading to increased hospital admissions [9,10,11]. Additionally, respiratory and cardiovascular illnesses, especially perpetual obstructive aspiratory sicknesses, are some of the main causes of mortality in worldwide [12]. The above-mentioned issues and the enhanced emphasis on public health risk due to ambient air pollutants necessitate the estimation of the association between air pollution exposure and adverse public health effects. At present, public health risk assessments are performed by the human health risk assessment framework originating from the US Environmental Protection Agency (USEPA). This framework is a handy tool that is used to quantitatively calculate human health risk resulting from exposure to ambient air pollutants [13]. Further, inhalation of different atmospheric particulates by children and adults will lead to a different degrees of health risk, as investigated based on the above-mentioned framework [14]. In addition, the current application of this framework mainly considers people as having average exposure or continuous exposure [15].
Resilience is a term with various symbolic meanings and has become a popular theme across many disciplines, including health resilience, organizational resilience, psychological resilience, fund resilience, and urban resilience, etc. At present, urban resilience has been widely used, which has been applied to flash flood risk based on environmental and social indices (population density, building density, economic density, farmland density, road density, and urbanization rate, etc.) [16]. As to the risk caused by ambient air pollutants, one previous study defined the resilience to air pollution (mainly distributed in urban area) as the capacity of local stakeholders to formulate corresponding measures focused on hazard (or vulnerability) [2]. Furthermore, seven air quality resiliency factors for urban systems were defined by previous research, and included natural emissions mitigation, anthropogenic emissions, environmental policy, urban form and land-cover change, access to technology, geography, and socio-cultural dimensions of risk perception [17]. However, there are still limited studies that have used resilience to tackle air pollution, but health resilience can be seen with large and small risk, and the ability of people to recover from this risk [18].
Generally, although many studies have investigated the spatial and temporal patterns, chemical composition, potential sources, and effects of meteorological condition based on six criteria air pollutants, these studies have either ignored public health risk assessment [19,20,21,22,23,24,25], or only focused on a specific city or industrial region [3,15,26]. In the face of increasing risks caused by air pollutants, the public health resilience of different groups of people in Sichuan still needs a systematic and comprehensive analysis. Moreover, the reasonable spatial distribution of hospitals will play a vital role in people’s health, but this region lacks a study on spatial inequality analysis between the risk caused by air pollutants and hospital density. Therefore, the present study attempted to accomplish the following: (1) an assessment of the spatial-temporal patterns of six criteria air pollutants in Sichuan; (2) an evaluation of the risk caused by ambient air pollutants to different age groups on the spatial scale; (3) measurement of the public health resilience to the risk resulting from ambient air pollutants; and (4) analysis of the spatial inequality between the risk caused by air pollutants and hospital density in Sichuan.

2. Data and Methods

2.1. Study Area

Sichuan Province (26°03′–34°19′ N, 97°21′–108°33′ E), located in the transitional zone between the middle and lower reaches of Yangtze River and the Tibetan Plateau in southwest China, contains 21 cities (Figure 1 and Table 1). There are four types of topography: hills, mountains, plains, and plateaus. The annual average temperature is 16–18 °C and the precipitation is 1000–1200 mm in Sichuan Basin, although for northwestern Sichuan these values are 4–12 °C and 500–900 mm, respectively. In addition, Sichuan Province covers 486,000 km2 and has a population of 83.41 million; its gross domestic product (GDP) exceeded 4.06 trillion yuan at the end of 2018. Sichuan is thus an important province in southwest China. With the rapid development of the economy, the pace of urbanization has further increased, and the urbanization rate of the province’s permanent residents reached 50.4% in 2017 [27]. The specific conditions (e.g., geographical position, weather, and social and economic development) may exert a significant effect on air pollution characteristics.

2.2. Data

2.2.1. Air Quality Data

In order to understand the characteristics of the six criteria air pollutants and the health risk assessment of main air pollutants in Sichuan, real-time hourly air pollutant data (PM2.5, PM10, SO2, NO2, CO, and O3) from January 2015 to December 2017 were obtained from the China National Environmental Monitoring Center [28], including 126 air-monitoring sites in and around Sichuan (Figure 1). The sites were set in each city according with the China Environmental Protection Standards HJ 655-2013 and HJ 193-2013 [29]. All of the air pollutants were monitored every 1 h except for O3, which was measured as an average over every 8 h. An advisable check was conducted at each site to remove faulty data before publishing these hourly air pollutant data [30].
For the spatial distribution of the six criteria air pollutants in different seasons, the seasonal concentrations of air pollutants were calculated by averaging the real-time hourly data in each site. On this basis, we used the inverse distance weight interpolation method in ArcGIS 10.2 software (Esri, Redlands, CA, USA) to get the air pollution raster data, and the resolution was 1 km × 1 km.

2.2.2. Resilience Assessment Data

The spatial datasets were used to evaluate the public health resilience caused by air pollutants, including the hospital, economy, road density, land use, and night-time light data that related to various social coping capacity. Hospitals are an important indicator for the resilience assessment of public health, but the latitude and longitude data of hospital sites are not directly available at present. In this study, we collected hospital information from the Yao Zhi Data [31], which include information such as name, level, type, and address, and then we obtained the hospital’s longitude and latitude from the Baidu Map using the address information. This paper collected hospital data, but did not include sanitation stations, epidemic prevention and control centers, or community health care centers. In addition, we also selected economy, road density, land use, and night-time light except for hospital data, as these indicators can reflect the coping capacity to deal with the risks caused by air pollutants [32], and the original data sources are shown in Table 2. The data of pertussis vaccination rate was collected from Sichuan Health Statistics Yearbook (2017) [33].
In order to get the raster data of hospital density, the point density analysis was performed in ArcGIS 10.2 software. Previous studies have illuminated that the driving speeds during the trip to the hospital on a freeway, national highway, provincial highway, county highway, and village road are 100, 60, 50, 30, and 20 km/h, respectively, leading to a journey lasting about 1 h on these roads [34]. Therefore, we assumed that the service range of each hospital was 50 km, and took that as the radius for the point density analysis. The road density was calculated in ArcGIS 10.2 software based on the linear density analysis; the roads included freeways, national highways, provincial highways, county highways, and village roads. In addition, six land uses types were assigned values to each cell based on the intensity of human activity as 1, 2, 3, 4, 5, and 6 for unused land, water area, forest land, meadow, cultivated land, and urbanized area, respectively. For the night-time light, we calculated the three-year average as the indicator input. The final resolution of all the data was unified to 1 km.

2.3. Methods

To research the health risk and resilience assessment resulting from main air pollutants, and spatial inequality between the risk caused by ambient air pollutants and hospital density in Sichuan, we first analyzed the spatial characteristics and variation of the six criteria air pollutants. This section is organized as follows. Section 2.3.1 describes the health risk assessment of the main air pollutants. Section 2.3.2 presents the health resilience assessment. The Lorenz curve and Gini coefficient are presented in Section 2.3.3.

2.3.1. Health Risk Assessment of Main Air Pollutants

Health risk assessment is a quantitative method for determining risks to human health caused by exposure to ambient air pollutants [14]. In the present study, it was revealed that PM2.5, PM10, SO2, NO2, CO, and O3 are harmful to public health and associated health risks [35]. Our study mainly explores the chronic (annual) exposure period, whereas the annual mean concentrations published by the ambient air quality standard (GB3095-2012) in China only include PM2.5, PM10, SO2, and NO2 (Table 3) [36]. Therefore, the health risk assessment takes the above-mentioned air pollutants (PM2.5, PM10, SO2, and NO2) as the research objects. In addition, the health risk assessment was performed among different age groups, namely, adults (19–75 years), children (6–12 years), and infants (birth to a year of age) [15].
Risk characterization is performed by the hazard quotient (HQ), which can reflect the probability of an adverse air pollutant occurring among healthy (or sensitive) individuals, thus avoiding the impacts of population density. An HQ ≤ 1.0 indicates a low risk, that is, the concentration of air pollutants is not likely to induce adverse health effects. An HQ > 1.0 indicates that the concentration of air pollutants may induce adverse effects on human health. The HQ was calculated as follows [14]:
H Q = A D D / R E L
where REL is the reference exposure level (REL). Based on the ambient air quality standard (GB3095-2012) [36], this paper takes the average value of Grade I and Grade II as the reference exposure level, as adopted by the previous research [37]. ADD is the daily dose, and is calculated as follows [38]:
A D D = ( C A × I R × E D ) / ( B W × A T )
where CA denotes the concentration of air pollutants (μg m−3), IR denotes the average values of inhalation rate (m3 day−1), ED denotes the exposure duration (day), BW denotes the average weight (kg), and AT denotes the averaging time (day). All parameter standards are derived from the US Environmental Protection Agency [39], which are shown in Table 4.

2.3.2. Health Resilience Assessment

In this paper, the health resilience assessment is defined as the capacity of public health to recover from the risks caused by ambient air pollutants [18]. Based on the result of health risk assessment, we focused health resilience assessment on air pollutants that pose a high risk to public health, and the process is as follows: (1) calculation of the social coping capacity; and (2) calculation of health resilience to the risk caused by air pollutants.
1. Calculation of social coping capacity
In order to quantify social coping capacity, principal component analysis was selected to assist us because this method is the oldest multivariate technique and can reflect important information from the original indicators [40,41]. In this study, the standardized results of the hospital, GDP, road density, land use, and night-time light were transformed into two principal components in ArcGIS 10.2 software. The number of principal components was determined based on previous studies [42], and the contribution rates of these two principal components to the original indicators information were 0.742 and 0.224, respectively. The social coping capacity is calculated as follows [43]:
S = i = 1 n C i P i
where n is the number of principal components, Ci is the contribution rate of principal component i, and Pi is the standardized result of principal component i.
2. Calculation of health resilience
Health resilience is largely determined by the social coping capacity and the risks caused by air pollutants: it is the existence of high social coping capacity and low risks that can lead to high health resilience. The health resilience to the risk caused by air pollutants can be expressed as follows [44]:
R = S s / H Q s
where R is the health resilience to the risk caused by air pollutants, Ss is the standardized result of social coping capacity, and HQs is the standardized result of HQ. An R that is ≤1.0 indicates a low health resilience, that is, the social coping capacity is lower than risks caused by air pollutants. Meanwhile, R > 1.0 indicates a high health resilience, that is, the social coping capacity can deal with risks caused by air pollutants.

2.3.3. Lorenz Curve and Gini Coefficient

The Lorenz curve is a model used in economics to measure the equitability of the distribution of social resources across the population, and the Gini coefficient can quantify the extent of inequality between two variables [45]. Previous studies identified the regional disparity of emission distribution of air pollutants and socio-economic condition by the Lorenz curve and Gini coefficient [46]. In this study, to calculate the spatial inequality between the risk caused by air pollutants and hospital density (risk is represented by HQ, hospital density is quantified by density analysis result of hospital point), this relationship is performed as a curve with a concave slope (Figure 2). The perfect equality line is the upper limit of the Lorenz curve, and a greater gap between the perfect equality line and Lorenz curve that leads to greater distribution inequality. The Gini coefficient evaluates the extent of inequality by comparing the area between the Lorenz curve and the perfect equality line (A) and the area under the equality line (A + B) (Figure 2) [47]:
G = A A + B
where G ranges from 0 to 1; G = 0 indicates that the risk and hospital density are in balance. The larger Gini coefficient leads to greater spatial imbalance between the risk and hospital density. The Lorenz curve and Gini coefficient were obtained using the package “lawstat” in R software (R Core Development Team, R Foundation for Statistical Computing, Vienna, Austria).

3. Results

3.1. Spatial and Temporal Patterns of Six Criteria Air Pollutants

According to the spatial distribution of six criteria air pollutants in Sichuan from spring to winter (Figure 3), we discovered that the PM2.5 and PM10 have the same pattern, with the maximum distribution in Sichuan Basin, and the minimum distribution in the northwest part of Sichuan (Figure 3a–h). Furthermore, the highest concentrations of PM2.5 and PM10 occurred in winter, and ranged from 17.3 to 112.9 μg m−3 and from 32.3 to 157.3 μg m−3, respectively. According to Figure 3i–l,q–t, Panzhihua had the highest concentrations of SO2 and CO in spring (29.9 μg m−3 and 1.3 mg m−3, respectively), summer (28.2 μg m−3 and 1.3 mg m−3, respectively), fall (36.6 μg m−3 and 1.6 mg m−3, respectively), and winter (40.6 μg m−3 and 4.7 mg m−3, respectively). The concentration ranges of NO2 were 7.3–61.4 μg m−3, 5.1–51.4 μg m−3, 6.2–59.4 μg m−3, and 11.0–70.4 μg m−3, in spring, summer, fall, and winter, respectively (Figure 3m–p); the maximum values were primarily distributed in Chengdu. However, the highest O3 concentration occurred in summer (25.7–101.4 μg m−3), and the lowest concentration occurred in winter (11.1–66.9 μg m−3) (Figure 3v,x). We also found that in summer the maximum O3 concentration was chiefly located in Sichuan Basin, while the minimum concentration was in Bazhong.

3.2. Health Risk Assessment of Chronic Exposure Across the Population

To better understand the health risk assessment of chronic (annual) exposure in Sichuan, the risk caused by PM2.5, PM10, SO2, and NO2 are shown in Figure 4. For infants, according to Figure 4a, the high-risk values of PM2.5 are mainly distributed in Leshan (1.25), Meishan (1.18), Mianyang (1.05), Nanchong (1.09), Neijiang (1.35), Chengdu (1.23), Dazhou (1.15), Deyang (1.19), Suining (1.08), Yibin (1.34), Ziyang (1.06), Zigong (1.62), and Luzhou (1.21). The high-risk values of PM10 occurred mainly in Zigong (1.01) (Figure 4d). However, no high-risk values of SO2 and NO2 were observed in Sichuan (Figure 4g,j). For children and adults, there were no high-risk areas for PM2.5, PM10, SO2, and NO2 (Figure 4b,c,e,f,h,i,k,l).

3.3. Health Resilience Assessment of Infants for the Risk Caused by PM2.5 and PM10

On the basis of the risk caused by PM2.5 and PM10 to the infants, the health resilience assessment is shown in Figure 5. For the risk caused by PM2.5, the high values of health resilience were mainly distributed in the southwest of Aba and the south of Ganzi, and the low values were around the borders of Aba and Chengdu, Deyang, and Mianyang, as well as in Leshan, Yibin, Zigong, Dazhou and Luzhou (Figure 5a). As to the average health resilience of infants to the risk caused by PM2.5, the cities of Ganzi (4.79), Liangshan (3.91), Aba (3.89), and Guangyuan (3.42) have high health resilience. The health resilience of infants to the risk caused by PM10 had the same pattern as PM2.5, but with Yaan having a large area of low health resilience (Figure 5b).
Since research assessing the health resilience to the risk caused by ambient air pollutants is in progress, there is a current lack of studies that verify the validity of health resilience assessment. Hence, a correlation analysis was conducted on the pertussis vaccination rate and average health resilience of infants to the risk caused by PM2.5 in 21 cities in Sichuan. There is a correlation coefficient of 0.79, which indicates that a better health resilience to the risk caused by air pollutants will lead to weaker vaccination awareness (Figure 5c).

3.4. Spatial Inequality Analysis between Risk Caused by Air Pollutants and Hospital Density

As shown in Figure 6 and Figure 7, we conducted a spatial inequality analysis between risk caused by PM2.5 (PM10) and hospital density. High spatial inequalities level between risk caused by PM2.5 and hospital density mainly occurred in Aba (G = 0.518) and Ganzi (G = 0.493), whereas low spatial inequalities were chiefly concentrated in Meishan (G = 0.06), Chengdu (G = 0.04), Deyang (G = 0.072), Panzhihua (G = 0.059), Suining (G = 0.035), Ziyang (G = 0.066), and Zigong (G = 0.087) (Figure 6a,c,e,f,l,p,q,t,u). In addition, we found that the entire province of Sichuan also remained high spatial inequality (G = 0.456) (Figure 6v). The spatial inequality between risk caused by PM10 and hospital density had the same trend as PM2.5: Aba (G = 0.525) and Ganzi (G = 0.526) had high spatial inequalities, while low spatial inequalities primarily occurred in Meishan (G = 0.037), Neijiang (G = 0.064), Chengdu (G = 0.036), Deyang (G = 0.058), Panzhihua (G = 0.068), Suining (G = 0.021), Ziyang (G = 0.028), and Zigong (G = 0.084) (Figure 7a,c,e,f,l,o–q,t,u). Simultaneously, Sichuan also had high spatial inequality between risk caused by PM10 and hospital density had the same trend as PM2.5 (G = 0.466) (Figure 7v).

4. Discussion

4.1. Spatial and Temporal Patterns of Air Pollutants

Our study identified that the Sichuan Basin had high concentrations of PM2.5 and PM10, and low concentrations were mainly distributed in the northwestern Sichuan region (Figure 3a–h). A few studies have revealed that the Sichuan Basin had high populations, dense urban traffic, rapid urban expansion, and large exhaust emissions, which are the direct causes of the high concentration of suspended particles in this region [29,48]. Furthermore, a previous study reported that industrial emissions contributed 27% to PM2.5 concentration in Wuhan during 2011–2012 [49]. However, the northwestern Sichuan region is a high-altitude area, and the local residents mainly rely on grazing [50], and thus low industrial activity might be the main reason for the PM2.5 and PM10 patterns. We also found that Panzhihua has the highest concentrations of SO2 and CO (Figure 3i–l,q–t). Usually, secondary industry plays an important role in the high concentration of SO2 and CO. For example, a previous study showed that steel manufacturing contributed 24% and 30% to SO2 and CO concentrations, respectively, in Tianjin [51]. Therefore, the high concentrations of SO2 and CO in Panzhihua are attributed to the more developed steel industry in this region [52]. As to NO2, Chengdu had the highest concentration (Figure 3m–p), this region had 3128 thousand vehicles and a GDP of 1005.7 billion yuan in 2015 [29], which we consider to have contributed to the high emissions of NO2 that was supported by the previous research [53].
The spatial characteristics of PM2.5, PM10, SO2, NO2, and CO were basically the same in each season, and the highest concentrations occurred in winter (Figure 3a–t). This pattern is a typical characteristic in East Asia [54], the concentration of air pollutants in winter are chiefly affected by emission sources and meteorological conditions. Usually, Sichuan has frequent stagnant weather conditions in winter, characterized by slow winds, little precipitation, and shallow mixing layers; this special condition can capture the air pollutants emitted by local sources (or air pollutants transported from other regions) and elevates winter pollutant concentrations [55]. Coal combustion is another important factor influencing the air pollutant patterns in winter. Previous studies have reported that emissions of CO (32%), NO2 (70%), SO2 (90%), and PM2.5 (97%) were attributed to coal combustion [56]. Regarding O3, which showed different seasonal distributions, the lowest and highest concentrations occurred in winter and summer, respectively. In summer, the high concentration of O3 was mainly distributed in Sichuan Basin (Figure 3u–x); research based on fossil fuels and shortwave radiation provided direct evidence that the high fossil fuel consumption and strong solar radiation in the summer significantly influenced the O3 concentration in this region [57,58]. In winter, the highest concentration of O3 was concentrated in Liangshan; a previous study confirmed that the increase in O3 concentration in this area is mainly caused by the secondary conversion of inorganic ions [59].

4.2. Health Risk Assessment of the Population

Air pollutants remain a global environmental threat and cause a public health risk. Some researchers assumed that health risk from exposure to ambient air pollutants can happen below or at the concentrations allowed by the international and national air quality standards. Findings in this study revealed that, for infants, Leshan (HQ = 1.25), Meishan (HQ = 1.18), Mianyang (HQ = 1.05), Nanchong (HQ = 1.09), Neijiang (HQ = 1.35), Chengdu (HQ = 1.23), Dazhou (HQ = 1.15), Deyang (HQ = 1.19), Suining (HQ = 1.08), Yibin (HQ = 1.34), Ziyang (HQ = 1.06), Zigong (HQ = 1.62), and Luzhou (HQ = 1.21) had high health risks caused by PM2.5 (Figure 4a). Evidence from the PM2.5 concentration shows the ambient concentration (47.7 μg m−3) in the Sichuan Basin exceeded the Grade II standard (35 μg m−3). Furthermore, previous research revealed that respiratory infection disease in infants (in the absence of congenital diseases) are positively correlated with PM2.5 concentration, and they also reported that the Hengshui People’s Hospital received a total of 160 infants with diseases related to air pollutants from June 2014 to January 2015 [60]. This study further revealed that the high health risk of infants caused by PM10 occurred mainly in Zigong (HQ = 1.01) (Figure 4d). Zigong had the highest average annual concentration of PM10 (87.0 μg m−3) in Sichuan and was higher than the Grade II standard (70 μg m−3). Moreover, it was confirmed that the concentration of air pollutants on the ground tends to be higher [61]; hence, infants are more likely to inhale PM, which leads to higher health risks. A previous study also reported that a 10 μg m−3 increase in PM10 concentration can lead to an increase in hospital admissions for ischemic heart disease and congestive heart failure [62]. Our study showed that SO2 and NO2 did not pose high health risks for infants (Figure 4g,j); however, another researcher revealed that Chengdu had high health risks caused by NO2 in 2014 and 2015 [38]. The low risk after 2015 might be supported by Chengdu’s efforts to improve air quality in recent years [63]. Regarding children and adults, the PM2.5, PM10, SO2, and NO2 did not cause high health risk in Sichuan (Figure 4b,c,e,f,h,i,k,l). It has been confirmed that infants have a higher susceptibility to ambient air pollutants than children and adults. Furthermore, infants are considered a health risk group for numerous reasons, chiefly due to the relatively higher amount of air inhalation (as their air intake per weight unit in a resting state is twice that of an adult), and the fact that the lungs and immune system are not yet fully developed [13].

4.3. The Health Resilience Assessment and Spatial Inequality Analysis

The health resilience assessment clearly presented the health resilience of infants to the risk caused by PM2.5 and PM10, and our spatial inequality analysis revealed the spatial inequality between risk caused by PM2.5 (PM10) and hospital density in Sichuan. Our study demonstrated that a high health resilience to the risk caused by PM2.5 was mainly distributed in Aba and Ganzi (Figure 5a). These cities have approximately 53.7% and 61.7% natural grassland, respectively [64,65]; furthermore, local residents mainly rely on grazing, and there are lower levels of industrial emissions and automobile exhaust in these two regions [50]. In this condition, air pollution did not pose a high risk, and therefore the people living there have a high level of health resilience to the low risk caused by PM2.5. Additionally, the Ten-Year Action Plan for Health Development in Ethnic Regions of Sichuan Province, a significant plan for the medical development of ethnic areas launched in 2011, improved the medical conditions and promoted health resilience in Aba and Ganzi [66]. However, there is still a high spatial inequality between risk caused by PM2.5 and hospital density in Aba (G = 0.518) and Ganzi (G = 0.493) (Figure 6a,f), which indicates that these two regions need to continue to strengthen their medical systems with a reasonable spatial distribution. In addition, we further revealed that low health resilience to the risk caused by PM2.5 was principally concentrated in Leshan, Yibin, Zigong, Dazhou, and Luzhou (Figure 5a). A previous study reported that Sichuan Basin had a high level of PM2.5 concentration, and a cumulative concentration of PM2.5 with a three-day lag increased by 1 inter-quartile range (IQR) will led to the number of asthma patients increasing by 3.27% [67]. Additionally, Leshan (47.31%), Yibin (45.10%), Zigong (47.88%), Dazhou (40.87%), and Luzhou (46.08%) have a lower urbanization rates than Chengdu (71.47%) [68], and their spatial inequalities between risk caused by PM2.5 and hospital density were greater than that of Chengdu (Figure 6c,d,i,k,s,u), which together determined that the low health resilience to the risk caused by PM2.5 in the above-mentioned cities. Furthermore, we also found that the health resilience of infants to the risk caused by PM10 was the same as that for PM2.5, which is supported by previous researches based on the correlation analyses between PM10 and PM2.5 [35,69].

4.4. Uncertainties and Limitations

Our results showed that PM2.5 and PM10 pose a high risk to infants. The high risk caused by PM2.5 is mainly distributed in Sichuan Basin, and that caused by PM10 principally occurs in Zigong. Additionally, the high health resilience of infant to the risks caused by PM2.5 and PM10 were chiefly concentrated in Aba and Ganzi. Although uncertainties exist in human health risk and resilience assessments, their application have found that the usefulness in providing a consistent and quantitative method to systematically evaluate public health risk and resilience caused by air pollutants, and inform decisions for air pollution control. Human health resilience assessment used in this study is conservative as it includes many social factors that are built into the social coping capacity. To address these uncertainties, we verified the validity of the health resilience assessment by analyzing the correlation between the pertussis vaccination rate and average health resilience of infants to the risk caused by PM2.5 in 21 cities in Sichuan.
The findings in this study should be considered in light of the following limitations. The public health risk and resilience assessment in this study adopted different age groups of people as the unit of research rather than population density. The AT used in this study for adult group is the average value (30 × 365), which may bring errors to the risk assessment results of young people and the elderly. As to the validity of the health resilience assessment, we measured the correlation between the pertussis vaccination rate and average health resilience of infants to the risk caused by PM2.5 in 21 cities in Sichuan, but the local economic and quality of public health service play an important role in controlling the pertussis vaccination rate. We assumed that the service capabilities of all hospitals are the same. Also, public health resilience refers to the ability of the public health to recover from the risk caused by air pollutants. The high social coping capacity and low risks caused by air pollutants can determine the high public health resilience.
The strengths of our study are worth mentioning. We described the public health risk and resilience associated with human exposure to PM2.5, PM10, SO2, and NO2 on the spatial scale in Sichuan, which can inform the decisions for air pollutant control. In addition, the spatial inequality analysis between risk caused by air pollutants and hospital density can reflect the spatial distribution of medical care in a region, thus providing a certain reference for the government to build medical services.

5. Conclusions

Ambient air pollution is harmful to public health. Our study focused on the spatial-temporal characteristics of air pollutants and the public health risk and resilience assessment with respect to the main air pollutants. We also analyzed the spatial inequality between risk caused by air pollutants and hospital density in Sichuan. High concentrations of PM2.5 and PM10 were mainly concentrated in Sichuan Basin, which posed a high risk to infants. The high risk caused by PM2.5 was mainly distributed in Sichuan Basin, and that caused by PM10 principally occurred in Zigong. Additionally, Aba and Ganzi had high health resilience of infants to the risks caused by PM2.5 and PM10, while low health resilience was chiefly concentrated in Leshan, Yibin, Zigong, Dazhou, and Luzhou. Furthermore, Sichuan, Aba, and Ganzi had serious spatial inequality between risk caused by PM2.5 (PM10) and hospital density.
This study attempted to evaluate the health resilience to the main air pollutants in Sichuan, but it is conservative as it includes limited social factors (hospital, GDP, road density, land use, and night-time light). However, the identification of the potential of air pollution to pose public health hazards and the measured public health resilience to the risks caused by air pollutants provide a valuable contribution for environmental departments, and will aid in the creation of more concrete policies to protect and prolong human lives. On the basis of the results of this study, policymakers should strengthen existing legislation to limit the release of PM2.5 and PM10 in Sichuan Basin, especially in Zigong; moreover, Aba and Ganzi need to continue to strengthen their medical systems, with a reasonable spatial distribution.

Author Contributions

Conceptualization, J.X. and C.Y.; formal analysis, C.Y.; data and resources, C.Y.; writing—original draft preparation, J.X. and C.Y.; writing—review and editing, T.Z.; supervision, J.X.; funding acquisition, J.X. and W.C.

Funding

This research was supported by the Strategic Priority Research Program of Chinese Academy of Sciences (XDA20030302), Science and Technology Project of Xizang Autonomous Region (Grant No. XZ201901-GA-07), Open Subject of Big Data Institute of Digital Natural Disaster Monitoring in Fujian (NDMBD2018003), Southwest Petroleum University of Science and Technology Innovation Team Projects (2017CXTD09), and National Flash Flood Investigation and Evaluation Project (SHZH-IWHR-57).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stockman, J.A. Main air pollutants and myocardial infarction: A systematic review and meta-analysis. JAMA 2013, 307, 713–721. [Google Scholar] [CrossRef]
  2. Cariolet, J.M.; Colombert, M.; Vuillet, M.; Diab, Y. Assessing the resilience of urban areas to traffic-related air pollution: Application in Greater Paris. Sci. Total Environ. 2018, 615, 588–596. [Google Scholar] [CrossRef] [PubMed]
  3. Kumar, A.; Mishra, R.K. Human health risk assessment of major air pollutants at transport corridors of Delhi, India. J. Trans. Health 2018, 10, 132–143. [Google Scholar] [CrossRef]
  4. Huang, R.J.; Zhang, Y.; Bozzetti, C.; Ho, K.F.; Cao, J.J.; Han, Y.; Daellenbach, K.R.; Slowik, J.G.; Platt, S.M.; Canonaco, F.; et al. High secondary aerosol contribution to particulate pollution during haze events in China. Nature 2014, 514, 218–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Betha, R.; Behera, S.N.; Balasubramanian, R. 2013 Southeast Asian smoke haze: Fractionation of particulate-bound elements and associated health risk. Environ. Sci. Technol. 2014, 48, 4327–4335. [Google Scholar] [CrossRef]
  6. Yu, G.; Wang, F.; Hu, J.; Liao, Y.; Liu, X. Value assessment of health losses caused by PM2.5 in Changsha City, China. Int. J. Environ. Res. Public Health 2019, 16, 2063. [Google Scholar] [CrossRef] [PubMed]
  7. Yun, G.; He, Y.; Jiang, Y.; Dou, P.; Dai, S. PM2.5 spatiotemporal evolution and drivers in the Yangtze River Delta between 2005 and 2015. Atmosphere 2019, 10, 55. [Google Scholar] [CrossRef]
  8. Lim, S.S.; Vos, T.; Flaxman, A.D.; Danaei, G.; Shibuya, K.; Adair-Rohani, H.; AlMazroa, M.A.; Amann, M.; Anderson, H.R.; Andrews, K.G.; et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: A systematic analysis for the Global Burden of Disease Study 2010. Lancet 2012, 380, 2224–2260. [Google Scholar] [CrossRef]
  9. Zhang, Z.; Wang, J.; Chen, L.; Chen, X.; Sun, G.; Zhong, N.; Kan, H.; Lu, W. Impact of haze and air pollution-related hazards on hospital admissions in Guangzhou, China. Environ. Sci. Pollut. Res. Int. 2014, 21, 4236–4244. [Google Scholar] [CrossRef]
  10. Liu, T.; Zhang, Y.H.; Xu, Y.J.; Lin, H.L.; Xu, X.J.; Luo, Y.; Xiao, J.; Zeng, W.L.; Zhang, W.F.; Chu, C.; et al. The effects of dust-haze on mortality are modified by seasons and individual characteristics in Guangzhou, China. Environ. Pollut. 2014, 187, 116–123. [Google Scholar] [CrossRef]
  11. Ghozikali, M.G.; Mosaferi, M.; Safari, G.H.; Jaafari, J. Effect of exposure to O3, NO2, and SO2 on chronic obstructive pulmonary disease hospitalizations in Tabriz, Iran. Environ. Sci. Pollut. Res. 2014, 22, 2817–2823. [Google Scholar] [CrossRef] [PubMed]
  12. Zanobetti, A.; Schwartz, J.; Dockery, D.W. Airborne particles are a risk factor for hospital admissions for heart and lung disease. Environ. Health Perspect. 2000, 108, 1071–1077. [Google Scholar] [CrossRef] [PubMed]
  13. Thabethe, N.D.; Engelbrecht, J.C.; Wright, C.Y.; Oosthuizen, M.A. Human health risks posed by exposure to PM10 for four life stages in a low socio-economic community in South Africa. Pan Afr. Med. J. 2014, 18, 206. [Google Scholar] [PubMed]
  14. Matooane, M.; Diab, R. Health risk assessment for sulfur dioxide pollution in South Durban, South Africa. Arch. Environ. Health Int. J. 2003, 58, 763–770. [Google Scholar] [CrossRef] [PubMed]
  15. Morakinyo, O.M.; Adebowale, A.S.; Mokgobu, M.I.; Mukhola, M.S. Health risk of inhalation exposure to sub-10 μm particulate matter and gaseous pollutants in an urban-industrial area in South Africa: An ecological study. BMJ Open 2017, 7, e013941. [Google Scholar] [CrossRef] [PubMed]
  16. Ding, M.; Heiser, M.; Hübl, J.; Fuchs, S. Regional vulnerability assessment for debris flows in China—A CWS approach. Landslides 2015, 13, 537–550. [Google Scholar] [CrossRef]
  17. Duh, J.D.; Shandas, V.; Chang, H.; George, L.A. Rates of urbanisation and the resiliency of air and water quality. Sci. Total Environ. 2008, 400, 238–256. [Google Scholar] [CrossRef] [PubMed]
  18. Wulff, K.; Donato, D.; Lurie, N. What is health resilience and how can we build it? Annu. Rev. Public Health 2015, 36, 361–374. [Google Scholar] [CrossRef]
  19. Chai, F.; Gao, J.; Chen, Z.; Wang, S.; Zhang, Y.; Zhang, J.; Zhang, H.; Yun, Y.; Ren, C. Spatial and temporal variation of particulate matter and gaseous pollutants in 26 cities in China. J. Environ. Sci. 2014, 26, 75–82. [Google Scholar] [CrossRef]
  20. Chen, Y.; Xie, S. Characteristics and formation mechanism of a heavy air pollution episode caused by biomass burning in Chengdu, Southwest China. Sci. Total Environ. 2014, 473–474, 507–517. [Google Scholar] [CrossRef]
  21. Coates, J.; Mar, K.A.; Ojha, N.; Butler, T.M. The influence of temperature on ozone production under varying NOx conditions–A modelling study. Atmos. Chem. Phys. 2016, 16, 11601–11615. [Google Scholar] [CrossRef]
  22. Dimitriou, K.; Kassomenos, P. Three year study of tropospheric ozone with back trajectories at a metropolitan and a medium scale urban area in Greece. Sci. Total Environ. 2015, 502, 493–501. [Google Scholar] [CrossRef] [PubMed]
  23. Gao, J.; Zhu, B.; Xiao, H.; Kang, H.; Hou, X.; Shao, P. A case study of surface ozone source apportionment during a high concentration episode, under frequent shifting wind conditions over the Yangtze River Delta, China. Sci. Total Environ. 2016, 544, 853–863. [Google Scholar] [CrossRef] [PubMed]
  24. Guan, Q.; Cai, A.; Wang, F.; Yang, L.; Xu, C.; Liu, Z. Spatio-temporal variability of particulate matter in the key part of Gansu Province, Western China. Environ. Pollut. 2017, 230, 189–198. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, H.; Tian, C.; Wang, W.; Luo, X. Temporal cross-correlations between ambient air pollutants and seasonality of tuberculosis: A time-series Analysis. Int. J. Environ. Res. Public Health 2019, 16, 1585. [Google Scholar] [CrossRef] [PubMed]
  26. Cai, K.; Li, C.; Na, S. Spatial distribution, pollution source, and health risk assessment of heavy metals in atmospheric depositions: A case study from the sustainable city of Shijiazhuang, China. Atmosphere 2019, 10, 222. [Google Scholar] [CrossRef]
  27. Sichuan Bureau of Statistics. Available online: http://tjj.sc.gov.cn/ (accessed on 5 May 2019).
  28. China National Environmental Monitoring Center. Available online: http://113.108.142.147:20035/emcpublish/ (accessed on 5 October 2018).
  29. Zhao, S.; Yu, Y.; Yin, D.; Qin, D.; He, J.; Dong, L. Spatial patterns and temporal variations of six criteria air pollutants during 2015 to 2017 in the city clusters of Sichuan Basin, China. Sci. Total Environ. 2018, 624, 540–557. [Google Scholar] [CrossRef]
  30. Hu, J.; Wang, Y.; Ying, Q.; Zhang, H. Spatial and temporal variability of PM2.5 and PM10 over the North China Plain and the Yangtze River Delta, China. Atmos. Environ. 2014, 95, 598–609. [Google Scholar] [CrossRef]
  31. Yao Zhi Data. Available online: https://db.yaozh.com/hmap (accessed on 16 December 2018).
  32. Bae, H.-J.; Kang, J.E.; Lim, Y.-R. Assessing the health vulnerability caused by climate and air pollution in Korea using the fuzzy TOPSIS. Sustainability 2019, 11, 2894. [Google Scholar] [CrossRef]
  33. Sichuan Health Statistics Yearbook. Available online: https: //www.xnjdcbs.com (accessed on 15 March 2019).
  34. Liu, Z.; Yang, H.; Xiong, W.; Chen, G. Spatial accessibilities of medical services at county level based on optimized two-step floating catchment area method. Sci. Geogr. Sin. 2017, 37, 728–737. [Google Scholar]
  35. Zhou, T.; Sun, J.; Yu, H. Temporal and spatial patterns of China’s main air pollutants: Years 2014 and 2015. Atmosphere 2017, 8, 137. [Google Scholar] [CrossRef]
  36. Ambient Air Quality Standard (GB3095-2012). Available online: http://www.jingbian.gov.cn/zdzl/jbxzwgksdzl/zwwgk/jcgk/hjbh/62929.htm (accessed on 26 December 2018).
  37. Zhong, M. Health Risk and Health Impact Assessment of Exposure to Air Pollutants in Xi’an; Xi’an University of Architecture and Technology: Xian, China, 2017. [Google Scholar]
  38. Sun, J.; Zhou, T. Health risk assessment of China’s main air pollutants. BMC Public Health 2017, 17, 212. [Google Scholar] [CrossRef] [PubMed]
  39. US Environmental Protection Agency. Available online: http://www.epa.gov/ncea/expofac.htm (accessed on 5 January 2019).
  40. Wang, H.; Zhang, B.; Liu, Y.; Liu, Y.; Xu, S.; Deng, Y.; Zhao, Y.; Chen, Y.; Hong, S. Multi-dimensional analysis of urban expansion patterns and their driving forces based on the center of gravity-GTWR model: Acase study of the Beijing-Tianjin-Hebei urban agglomeration. Acta Geogr. Sin. 2018, 73, 1076–1092. [Google Scholar]
  41. Xiong, J.; Ye, C.; Cheng, W.; Guo, L.; Zhou, C.; Zhang, X. The spatiotemporal distribution of flash floods and analysis of partition driving forces in Yunnan Province. Sustainability 2019, 11, 2926. [Google Scholar] [CrossRef]
  42. Abdi, H.; Williams, L.J. Principal component analysis. Wiley Interdiscip. Rev. Comput. Stat. 2010, 2, 433–459. [Google Scholar] [CrossRef]
  43. He, C.; Li, M.; Li, T.; Peng, X.; Li, J.; Zhao, J. Comparison and analysis of the four methods of determining weights in multi-objective comprehensive evaluation. J. Hubei Univ. (Nat. Sci.) 2016, 38, 172–178. [Google Scholar]
  44. Xiu, C.; Wei, Y.; Wang, Q. Evaluation of urban resilience of Dalian city based on the perspective of “Size-Density-Morphology”. Acta Geogr. Sin. 2018, 73, 2315–2328. [Google Scholar]
  45. Safi, M.N.; Zobitz, J.M. Quantifying distribution in carbon uptake and environmental measurements with the Gini coefficient. Lett. Biomath. 2016, 3, 1–12. [Google Scholar] [CrossRef]
  46. Dong, L.; Liang, H. Spatial analysis on China’s regional air pollutants and CO2 emissions: Emission pattern and regional disparity. Atmos. Environ. 2014, 92, 280–291. [Google Scholar] [CrossRef]
  47. Saha, D.; Kemanian, A.R.; Montes, F.; Gall, H.; Adler, P.R.; Rau, B.M. Lorenz curve and Gini coefficient reveal hot spots and hot moments for nitrous oxide emissions. J. Geophys. Res. Biogeosci. 2018, 123, 193–206. [Google Scholar] [CrossRef]
  48. Ke, B. Research on the Emission Inventory of Major Air Pollutants in 2012 for the Sichuan City Cluster in China; Southwest Jiaotong University: Chengdu, China, 2014. [Google Scholar]
  49. Cheng, H.; Wang, Z.; Feng, J.; Chen, H.; Zhang, F.; Liu, J. Carbonaceous species composition and source apportionment of PM2.5 in urban atmosphere of Wuhan. Ecol. Environ. Sci. 2012, 21, 1574–1579. [Google Scholar]
  50. Zhou, X. Tibetan area economic development of Sichuan Province profile in Ganzi, Aba, and Muli County. J. Southwest Univ. Natl. (Humanit. Soc. Sci.) 2008, 29, 27–31. [Google Scholar]
  51. Zhao, B.; Ma, J. Development of an air pollutant emission inventory for Tianjin. Acta Sci. Circumst. 2008, 28, 368–375. [Google Scholar]
  52. Zhong, X. The only way to upgrade the industrial structure of Panzhihua is to extend and develop the vanadium titanium industrial chain. Sci. Technol. Inf. 2013, 25, 165–168. [Google Scholar]
  53. Yan, J.; Wang, W.; Yang, N.; Jiang, Y. Effect of increasing motor vehicles on atmospheric nitrogen dioxide concentration in Chengdu. Sichuan Environ. 2012, 31, 34–36. [Google Scholar]
  54. Wang, Y.; Ying, Q.; Hu, J.; Zhang, H. Spatial and temporal variations of six criteria air pollutants in 31 provincial capital cities in China during 2013–2014. Environ. Int. 2014, 73, 413–422. [Google Scholar] [CrossRef] [PubMed]
  55. Tai, A.P.K.; Mickley, L.J.; Jacob, D.J. Correlations between fine particulate matter (PM2.5) and meteorological variables in the United States: Implications for the sensitivity of PM2.5 to climate change. Atmos. Environ. 2010, 44, 3976–3984. [Google Scholar] [CrossRef]
  56. Xie, Y.; Zhao, B.; Zhang, L.; Luo, R. Spatiotemporal variations of PM2.5 and PM10 concentrations between 31 Chinese cities and their relationships with SO2, NO2, CO and O3. Particuology 2015, 20, 141–149. [Google Scholar] [CrossRef]
  57. Li, G.; Bei, N.; Cao, J.; Wu, J.; Long, X.; Feng, T.; Dai, W.; Liu, S.; Zhang, Q.; Tie, X. Widespread and persistent ozone pollution in eastern China during the non-winter season of 2015: Observations and source attributions. Atmos. Chem. Phys. 2017, 17, 2759–2774. [Google Scholar] [CrossRef]
  58. Millán, M.M.; Mantilla, E.; Salvador, R.; Carratalá, A.; Sanz, M.J.; Alonso, L.; Gangoiti, G.; Navazo, M. Ozone Cycles in the Western Mediterranean Basin: Interpretation of Monitoring Data in Complex Coastal Terrain. J. Appl. Meteorol. 2000, 39, 487–508. [Google Scholar] [CrossRef]
  59. Zeng, S.; Zheng, Y. Analysis of a Severe Pollution Episode in December 2017 in Sichuan Province. Atmosphere 2019, 10, 156. [Google Scholar] [CrossRef]
  60. Liu, N. Relevance of Air Pollutants with Respiratory Tract Infection in Children in Hengshui; Hebei Medical University: Shijiazhuang, China, 2015. [Google Scholar]
  61. Shi, J. Analytical solution for the dispersion of pollutants from elevated source particles under gravity sedimentation. Environ. Res. Monit. 1985, 1, 1–5. [Google Scholar]
  62. Morris, R.D. Airborne particulates and hospital admissions for cardiovascular disease: A quantitative review of the evidence. Environ. Health Perspect. 2001, 109, 495–500. [Google Scholar] [PubMed]
  63. Qiao, X.; Jaffe, D.; Tang, Y.; Bresnahan, M.; Song, J. Evaluation of air quality in Chengdu, Sichuan Basin, China: Are China’s air quality standards sufficient yet? Environ. Monit. Assess. 2015, 187, 250. [Google Scholar] [CrossRef] [PubMed]
  64. Mou, D. Prevention and control of animal diseases or an important breakthrough for industrial poverty alleviation—Thinking and suggestions on the poverty alleviation in Sichuan Tibetan areas. Sichuan Anim. Vet. Sci. 2018, 45, 20–21. [Google Scholar]
  65. Yang, T.; Wang, D. Governance of prairie rat wasteland in Ganzi. Chin. J. Anim. Sci. 2016, 54–56. [Google Scholar] [CrossRef]
  66. Mou, J. Economic scholar Du Shouhu: Paying attention to people’s livelihood is the key to development. Sichuan Prov. Cond. 2013, 27–29. [Google Scholar] [CrossRef]
  67. Zhang, Z. Spatiotemporal Distributions of PM2.5 Using Ground Monitor Data and Its Relationship with Hospital Visiting of Respiratory Diseases; Zhejiang University: Hangzhou, China, 2016. [Google Scholar]
  68. Ning, G.; Wang, S.; Ma, M.; Ni, C.; Shang, Z.; Wang, J.; Li, J. Characteristics of air pollution in different zones of Sichuan Basin, China. Sci. Total Environ. 2018, 612, 975–984. [Google Scholar] [CrossRef]
  69. Chen, P.; Yang, J.; Pu, T.; Li, C.; Guo, J.; Tripathee, L.; Kang, S. Spatial and temporal variations of gaseous and particulate pollutants in six sites in Tibet, China, during 2016–2017. Aerosol Air Qual. Res. 2019, 19, 516–527. [Google Scholar] [CrossRef]
Figure 1. The study area.
Figure 1. The study area.
Ijerph 16 02796 g001
Figure 2. Lorenz curve diagram.
Figure 2. Lorenz curve diagram.
Ijerph 16 02796 g002
Figure 3. Spatial patterns of the six criteria air pollutants showing the seasonal mean concentrations from 2015 to 2017: (a) PM2.5 *, (e) PM10 *, (i) SO2 *, (m) NO2 *, (q) CO #, and (u) O3 * in spring; (b) PM2.5 *, (f) PM10 *, (j) SO2 *, (n) NO2 *, (r) CO #, and (v) O3 * in summer; (c) PM2.5 *, (g) PM10 *, (k) SO2 *, (o) NO2 *, (s) CO #, and (w) O3 * in fall; (d) PM2.5 *, (h) PM10 *, (l) SO2 *, (p) NO2 *, (t) CO #, and (x) O3 * in winter. Units: * μg m−3; # mg m−3.
Figure 3. Spatial patterns of the six criteria air pollutants showing the seasonal mean concentrations from 2015 to 2017: (a) PM2.5 *, (e) PM10 *, (i) SO2 *, (m) NO2 *, (q) CO #, and (u) O3 * in spring; (b) PM2.5 *, (f) PM10 *, (j) SO2 *, (n) NO2 *, (r) CO #, and (v) O3 * in summer; (c) PM2.5 *, (g) PM10 *, (k) SO2 *, (o) NO2 *, (s) CO #, and (w) O3 * in fall; (d) PM2.5 *, (h) PM10 *, (l) SO2 *, (p) NO2 *, (t) CO #, and (x) O3 * in winter. Units: * μg m−3; # mg m−3.
Ijerph 16 02796 g003
Figure 4. Risk assessment of chronic exposure to main air pollutants in Sichuan. (ac)/(df)/(gi)/(jl) represent hazard quotient for infants, children, and adults for PM2.5 (HQ_Infants_PM2.5, HQ_Children_PM2.5, and HQ_Adults_PM2.5)/PM10 (HQ_Infants_PM10, HQ_Children_PM10, and HQ_Adults_PM10)/SO2 (HQ_Infants_SO2, HQ_Children_SO2, and HQ_Adults_SO2)/NO2 (HQ_Infants_NO2, HQ_Children_NO2, and HQ_Adults_NO2), respectively.
Figure 4. Risk assessment of chronic exposure to main air pollutants in Sichuan. (ac)/(df)/(gi)/(jl) represent hazard quotient for infants, children, and adults for PM2.5 (HQ_Infants_PM2.5, HQ_Children_PM2.5, and HQ_Adults_PM2.5)/PM10 (HQ_Infants_PM10, HQ_Children_PM10, and HQ_Adults_PM10)/SO2 (HQ_Infants_SO2, HQ_Children_SO2, and HQ_Adults_SO2)/NO2 (HQ_Infants_NO2, HQ_Children_NO2, and HQ_Adults_NO2), respectively.
Ijerph 16 02796 g004
Figure 5. The health resilience of infants to the risk caused by (a) PM2.5 (Infants_Resilience_PM2.5) and (b) PM10 (Infants_Resilience_PM10); (c) the correlation between the pertussis vaccination rate and average health resilience of infants to the risk caused by PM2.5 in 21 cities in Sichuan.
Figure 5. The health resilience of infants to the risk caused by (a) PM2.5 (Infants_Resilience_PM2.5) and (b) PM10 (Infants_Resilience_PM10); (c) the correlation between the pertussis vaccination rate and average health resilience of infants to the risk caused by PM2.5 in 21 cities in Sichuan.
Ijerph 16 02796 g005
Figure 6. Spatial inequality analysis between the infant hazard quotient caused by PM2.5 (HQ_PM2.5) and hospital density: (a) Aba, (b) Bazhong, (c) Chengdu, (d) Dazhou, (e) Deyang, (f) Ganzi, (g) Guangan, (h) Guangyuan, (i) Leshan, (j) Liangshan, (k) Luzhou, (l) Meishan, (m) Mianyang, (n) Nanchong, (o) Neijiang, (p) Panzhihua, (q) Suining, (r) Yaan, (s) Yibin, (t) Ziyang, (u) Zigong, and (v) Sichuan.
Figure 6. Spatial inequality analysis between the infant hazard quotient caused by PM2.5 (HQ_PM2.5) and hospital density: (a) Aba, (b) Bazhong, (c) Chengdu, (d) Dazhou, (e) Deyang, (f) Ganzi, (g) Guangan, (h) Guangyuan, (i) Leshan, (j) Liangshan, (k) Luzhou, (l) Meishan, (m) Mianyang, (n) Nanchong, (o) Neijiang, (p) Panzhihua, (q) Suining, (r) Yaan, (s) Yibin, (t) Ziyang, (u) Zigong, and (v) Sichuan.
Ijerph 16 02796 g006
Figure 7. Spatial inequality analysis between the infant hazard quotient caused by PM10 (HQ_PM10) and hospital density: (a) Aba, (b) Bazhong, (c) Chengdu, (d) Dazhou, (e) Deyang, (f) Ganzi, (g) Guangan, (h) Guangyuan, (i) Leshan, (j) Liangshan, (k) Luzhou, (l) Meishan, (m) Mianyang, (n) Nanchong, (o) Neijiang, (p) Panzhihua, (q) Suining, (r) Yaan, (s) Yibin, (t) Ziyang, (u) Zigong, and (v) Sichuan.
Figure 7. Spatial inequality analysis between the infant hazard quotient caused by PM10 (HQ_PM10) and hospital density: (a) Aba, (b) Bazhong, (c) Chengdu, (d) Dazhou, (e) Deyang, (f) Ganzi, (g) Guangan, (h) Guangyuan, (i) Leshan, (j) Liangshan, (k) Luzhou, (l) Meishan, (m) Mianyang, (n) Nanchong, (o) Neijiang, (p) Panzhihua, (q) Suining, (r) Yaan, (s) Yibin, (t) Ziyang, (u) Zigong, and (v) Sichuan.
Ijerph 16 02796 g007
Table 1. The 21 cities in Sichuan Province.
Table 1. The 21 cities in Sichuan Province.
CityAbbreviationCityAbbreviationCityAbbreviation
AbaABGuangyuanGYNeijiangNJ
BazhongBZLeshanLSPanzhihuaPZH
ChengduCDLiangshanLSHSuiningSN
DazhouDZLuzhouLZYaanYA
DeyangDYMeishanMSYibinYB
GanziGZMianyangMYZiyangZY
GuanganGANanchongNCZigongZG
Table 2. Source and resolution of resilience assessment data.
Table 2. Source and resolution of resilience assessment data.
IndicatorsData SourceResolutionData TypeTime
HospitalChina: Yao Zhi Data2015
GDPChina: Resource and Environment Data Cloud Platform1 km × 1 kmRaster data2015
RoadsChina: National Earth System Science Data Sharing Infrastructure1:250,000Vector data2016
Land useChina: Resource and Environment Data Cloud Platform1 km × 1 kmRaster data2015
NPP-VIIRS Night-time lightUSA: National Oceanic and Atmospheric Administration500 × 500 mRaster data2015–2017
Table 3. Ambient air quality standard in China.
Table 3. Ambient air quality standard in China.
Air PollutantsAnnual Mean ConcentrationUnit
Grade IGrade II
PM2.51535μg m−3
PM104070μg m−3
SO22060μg m−3
NO24040μg m−3
Table 4. Parameters for health risk assessment of main air pollutants through the inhalation pathway. IR is the average values of inhalation rate, ED is the exposure duration, BW is the average weight, and AT is the averaging time.
Table 4. Parameters for health risk assessment of main air pollutants through the inhalation pathway. IR is the average values of inhalation rate, ED is the exposure duration, BW is the average weight, and AT is the averaging time.
Exposed GroupIR (m3 day−1)ED (day)BW (kg)AT (day)
Infant6.81 × 35011.31 × 365
Child13.512 × 35045.312 × 365
Adult13.330 × 350 71.830 × 365

Share and Cite

MDPI and ACS Style

Xiong, J.; Ye, C.; Zhou, T.; Cheng, W. Health Risk and Resilience Assessment with Respect to the Main Air Pollutants in Sichuan. Int. J. Environ. Res. Public Health 2019, 16, 2796. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152796

AMA Style

Xiong J, Ye C, Zhou T, Cheng W. Health Risk and Resilience Assessment with Respect to the Main Air Pollutants in Sichuan. International Journal of Environmental Research and Public Health. 2019; 16(15):2796. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152796

Chicago/Turabian Style

Xiong, Junnan, Chongchong Ye, Tiancai Zhou, and Weiming Cheng. 2019. "Health Risk and Resilience Assessment with Respect to the Main Air Pollutants in Sichuan" International Journal of Environmental Research and Public Health 16, no. 15: 2796. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152796

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