Next Article in Journal
Drought-Modulated Boreal Forest Fire Occurrence and Linkage with La Nina Events in Altai Mountains, Northwest China
Next Article in Special Issue
Allergic Rhinitis: Association with Air Pollution and Weather Changes, and Comparison with That of Allergic Conjunctivitis in Taiwan
Previous Article in Journal
Benzo[a]pyrene in the Ambient Air in the Czech Republic: Emission Sources, Current and Long-Term Monitoring Analysis and Human Exposure
Previous Article in Special Issue
Submicron and Ultrafine Particles in Downtown Rome: How the Different Euro Engines Have Influenced Their Behavior for Two Decades
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pollution Characteristics, Chemical Compositions, and Population Health Risks during the 2018 Winter Haze Episode in Jianghan Plain, Central China

1
School of Resource and Environmental Sciences, Wuhan University, Wuhan 430072, China
2
Xiantao Municipal Ecological Environment Bureau, Xiantao 433000, China
*
Authors to whom correspondence should be addressed.
Submission received: 26 August 2020 / Accepted: 4 September 2020 / Published: 7 September 2020
(This article belongs to the Special Issue Exposure Assessment of Air Pollution)

Abstract

:
To determine the pollution characteristics, chemical compositions, and population health risks of PM2.5 at different pollution levels, PM2.5 samples were intensively collected during the long-lasting winter haze episode from 13–23 January 2018 in Xiantao in Jianghan Plain (JHP), central China. The higher PM2.5 levels during the severe pollution period were dominated by the WNW-NNE air-masses, whereas the lower PM2.5 concentrations during other pollution periods were mainly affected by the NE, S, and NW air-masses. The NO3/SO42− and OC/EC ratios indicated a mixed contribution of intensive vehicle exhaust and secondary formation. The enrichment factor and geo-accumulation index for assessing the PM2.5-bound metal(loid)s contamination levels were positively correlated. Ingestion is the dominant exposure pathway of PM2.5-bound metal(loid)s for children and adults, followed by inhalation and dermal contact. As, Cr, and Pb may pose carcinogenic and non-carcinogenic risks, whereas Sb and V may only pose non-carcinogenic risks for children and adults. The population health risks may not depend on the pollution levels but depend on the PM2.5-bound metal(loid)s concentrations. PM2.5-bound metal(loid)s may pose much higher population health risks for adults compared to children. More attentions should be paid to the population health risks of PM2.5-bound metal(loid)s during a long-lasting winter haze episode in JHP.

Graphical Abstract

1. Introduction

Haze is associated with the explosive growth of airborne fine particulate matter (PM2.5) in ambient air [1]. PM2.5 has attracted worldwide concerns over the several years because of its adverse effects on atmospheric visibility and human health [2]. For example, PM2.5 can deposit in lungs through the inhalation exposure pathway and then result in lung cancer and acute respiratory infections [3]. In addition, it was classified as belonging to cancerogenic substance Group I by the International Agency for Research on Cancer [4]. However, the impacts of PM2.5 on population health rely on its chemical compositions, such as metal(loid)s, which represent all kinds of PM2.5-bound metals and metalloids [5], and organic components [6]. It was reported that PM2.5-bound metal(loid)s have influences on acute changes in cardiorespiratory physiology [7] and allergic airways disease [8], whereas PM2.5-bound elemental carbon is the main contributor to epidemiological disease [9]. Hence, necessary knowledge of PM2.5 chemical compositions is beneficial for the public and government to evaluate the population health risks.
Several studies have investigated the pollution characteristics and chemical compositions of PM2.5 in China, including Xi’an [10], Beijing-Tianjin-Hebei Region [11,12], Nanjing [13], Shanghai [14], Chengdu [15], Zhengzhou [16], Wuhan [17], and Xiangyang [18]. However, information regarding PM2.5 dynamic variations in Jianghan Plain (JHP), central China has not yet been reported so far. The study related to ambient air in JHP has focused on the spatiotemporal variation of atmospheric organochlorine pesticides [19] and ozone [20]. To our best knowledge, this is first the comprehensive study of the pollution characteristics, chemical compositions, and population health risks of PM2.5 in JHP, especially during a long-lasting winter haze episode.
JHP is located in the central and southern Hubei Province and is named after the alluvial deposits of the Yangtze River and Han River, where rice, cotton, and rape are mainly planted. With an average altitude of only 27 m, it is an important part of the plains of the middle and lower reaches of the Yangtze River in central China and is identified as a pollution hotspot [20]. Xiantao (30°04′~30°32′ N, 112°55′~113°49′ E) is a traditional agriculture city in JHP, central China, with an area and population of 2538 km2 and 1.14 million, respectively. According to the Xiantao Yearbook (http://tjj.hubei.gov.cn/tjsj/sjkscx/tjnj/gsztj/xts/), the gross domestic product (GDP) composed of primary industries (10.9%), secondary industries (51.6%), and tertiary industries (37.5%) of Xiantao reached 11.4 billion USD in 2018. These main industries in Xiantao include non-woven, electronic and electrical, textile and clothing, food processing, new materials and pharmaceutical manufacturing, and machinery and automotive parts, which accounted for 26.3%, 5.38%, 9.28%, 17.8%, 11.9%, and 11.6% of the GDP of industries, respectively. Xiantao’s car ownership increased from 171,264 in 2017 to 193,614 in 2018 by 13.1%. There are approximately 143 Gg of coal consumption for power plants and 13.5 Gg of coal consumption for residential biofuel, which are potential emission sources of PM2.5. According to the statistics of the Environmental Air quality Report of key cities in Hubei Province (http://sthjt.hubei.gov.cn/fbjd/zwgk/jcsjfb/hjkq/), the average PM2.5 concentration in January from 2012–2017 was 123 μg m−3, which was 1.64 times than the Grade II of National Ambient Air Quality Standards of China (NAAQS, 75 μg m−3) (GB 3095-2012, https://www.mee.gov.cn/ywgz/fgbz/bz/bzwb/dqhjbh/dqhjzlbz/201203/t20120302_224165.shtml), indicating a heavy PM2.5 pollution during wintertime in JHP, central China. As reported in Environmental Quality Bulletin of Xiantao City in 2018 (http://hbj.xiantao.gov.cn/hjzl/hjzlzggg/201908/P020200102465288147981.pdf), the frequency of haze events in this region is 26.8% in 2018. In January 2018, the frequency of haze events in this region is 50.2%. This makes it urgent for the government and public to take insight into the pollution characteristics, chemical compositions, and population health risks of PM2.5 for mitigating the haze episode in JHP.
In this study, a total of 43 pairs of PM2.5 samples were effectively collected during the long-lasting winter haze episode from 13–23 January 2018 in Xiantao in JHP, central China. The objectives of this study were to: (1) analyze the pollution characteristics, chemical compositions of PM2.5, including water-soluble inorganic ions (WSIIs), carbonaceous species (CAs), and trace elements (TE); (2) evaluate the contamination levels of PM2.5-bound metal(loid)s; (3) assess the population health carcinogenic and non-carcinogenic risks for children and adults due to exposure to PM2.5-bound metal(loid)s through ingestion, inhalation, and dermal contact pathways. This study provides important information for population health risks posed by PM2.5-bound metal(loid)s at four different pollution levels of Xiantao, which has a great praxis function as a reference for other cities in JHP, central China.

2. Materials and Methods

2.1. Field Sampling Campaign

According to the forecast results of Regional Air Quality Forecast System (http://aqf.climblue.com/) based on the Weather Research and Forecasting - Community Multiscale Air Quality (WRF-CMAQ) model provided by the International Institute for Earth System Science, Nanjing University, there would be a long-lasting winter haze episode from 13 January 2018 in Hubei Province, especially in Xiangyang, Jingmen, Wuhan, and Xiantao. Therefore, we conducted the intensive sampling campaign in 13 January 2018. PM2.5 field sampling campaign was carried out at an urban site from 13–23 January 2018 in Xiantao in JHP, central China. The geographical location of this sampling site is plotted in Figure 1. There is no obvious pollution source which may have effect on the sampling work around this urban sampling site, because the scale of industries near the sampling site, such as Ruiyang Automotive Parts (Xiantao) Co. Ltd. and Xiangyuan Electromechanical Equipment Co. Ltd. (Xiantao, China), are very small or do not have an independent production process. To track the haze bloom-decay process, samples were simultaneously collected on the rooftop of the library of Xiantao Vocational College (~20 m height, 30°20′27″ N, 113°25′16″ E) four times a day (at local times 06:00–11:00, 11:30–16:30, 17:00–22:00, and 22:30–05:30 the following day) with two medium-volume samplers (TH-150F, Wuhan Tianhong Instrument Co., Ltd., Wuhan, China) at an air flow rate of 100 L min−1 to increase the number of PM2.5 samples classified at different pollution levels as much as possible. This prevented the stoppage of sampling due to excessive filter resistance caused by a haze episode. However, heavy rains and snowy weather were observed from the evening of 23 January till 27 January 2018. So, we terminated this intensive sampling process on 24 January 2018 after the fourth sample collection was completed on 22:30 23 January–05:30 24 January 2018. To subtract possible contamination occurring during or after sampling, blank samples were collected for about 20 min by mounting blank filters onto the samplers without pumping into any air separately before and after the sampling [21]. In sum, a total of 43 pairs of PM2.5 samples as well as three pairs of blank samples were effectively collected on quartz fiber filters (Whatman, UK) and Teflon filters (Munktell, Sweden), respectively. The sample data on 06:00–11:00 21 January was missing because the power supply line was cut off by accident, which meant that the two samplers could not work normally during this period. The quartz fiber filters were baked at 500 °C for 6 h in a muffle furnace and the Teflon filters were prepared at a constant temperature and relative humidity (25 ± 1 °C, 50 ± 5%) for 48 h before use, respectively. An electronic microbalance with a resolution of 1 μg (Sartorius SECURA 125-1S, Sartorius Lab Instruments GmbH & Co. KG, Göttingen, Germany) was used to determine the PM2.5 mass on Teflon filters. After sampling, all the quartz fiber filters folded, wrapped with aluminum foil, then were sealed in plastic bags, and finally were stored in a refrigerator at −18 °C as well as Teflon filters for chemical analysis. According to the Technical Regulation on Air Quality Index of China (HJ 633-2012, http://www.mee.gov.cn/gkml/hbb/bgg/201203/t20120302_224146.htm), the sampling periods were divided into four pollution levels: mild, moderate, heavy, and severe pollution when the PM2.5 concentrations were between 75 and 115, 115 and 150, 150 and 250, and greater than 250 μg m−3, respectively. The sample numbers for mild pollution, moderate pollution, heavy pollution, and severe pollution levels are 17, 10, 10, and 3, respectively.

2.2. Chemical Analysis

The analysis method and procedures of WSIIs, CAs, and TE have been described in our previous studies [22,23,24,25]. To determine the WSIIs, a punch (3.14 cm2) of quartz fiber filter was extracted using 20 mL ultrapure Milli-Q water (18.25 MΩ cm). After 45 min ultrasonic ice-water bath extraction, the solution was filtered through a 0.22 μm pore size hydrophilic filter and then stored in a clean HDPE bottle. Nine WSIIs in the extraction solution, including four anions (F, Cl, NO3, and SO42−) and five cations (Na+, NH4+, K+, Mg2+, and Ca2+), were measured using an Ion Chromatography (883 Basic IC plus, Metrohm, Switzerland). The detection limits of F, Cl, NO3, SO42−, Na+, NH4+, K+, Mg2+, and Ca2+ were 0.010, 0.012, 0.027, 0.030, 0.019, 0.020, 0.025, 0.020, and 0.037 μg m−3, respectively. The OC and EC were analyzed using a thermal/optical reflectance carbon analyzer (DRI 2001 A, Desert Research Institute, USA). Briefly, a punch (0.518 cm2) of quartz fiber filter was stepwise heated at 140 °C (OC1), 280 °C (OC2), 480 °C (OC3), and 580 °C (OC4) in a pure He gas atmosphere for OC volatilization and at 580 °C (E1), 740 °C (E2), and 840 °C (E3) in a 2% O2/98% He atmosphere for EC oxidation for 150 s at each stage in the oven. The differences from replicate analyses on 10% of the total samples were within 3.95% for OC and 2.53% for EC. Standard concentrations of CH4/CO2 mixed gases were used to calibrate the analyzer each day. TE was measured using an Inductively Coupled Plasma-Mass Spectrometer (ICP-MS, PerkinElmer 350D, Boston, USA). Briefly, a Teflon filter was digested with a 10 mL acid mixture of 68% HNO3 and 40% HF (1:1, v/v) in a microwave system for 45 min. The digested solution was filtered through a 0.22 μm pore size hydrophilic filter after it cooled to room temperature and then was measured to determine the mass concentrations of twenty-five TE, including Li, Be, Al, Si, V, Cr, Mn, Co, Ni, Cu, Zn, As, Se, Sr, Mo, Ag, Cd, Sn, Sb, Ba, Tl, Pb, Bi, Th, and U. The relative standard deviations of soil national standard materials real values (Sigma-Aldrich, Saint Louis, USA) ranged from 0.10 to 18.00%. The detection limits of TE ranged from 0.00001 to 0.0005 μg L−1. Analyses of filed blank samples were performed using the above same methods. All the reported data of WSIIs, CAs, and TE were corrected by the field blanks.

2.3. Data Analysis Methods

2.3.1. Meteorological Parameters, Trace Gases, and Air Mass Back Trajectory Analysis

Hourly meteorological parameters in January 2018 (Supplementary Figure S1), including wind direction (WD), wind speed (WS), temperature (Temp.) relative humidity (RH), precipitation (Prec.), and visibility (Vis.) were collected from the nearest meteorological station to the sampling site (3 km, Figure 1), which can represent the ambient meteorology for the Xiantao site. Trace gases (Supplementary Figure S2), including SO2 and NO2, during the sampling time at Xiantao site were collected from the same place of Xiantao Industrial Park Station, which is a provincial controlled air quality monitoring station. To determine the long-range transport of air masses from different potential regions, 3 days (72 h) air mass back trajectories were calculated every 4 h (00:00, 06:00, 12:00, and 18:00 UTC) at 500 m A.G.L and were cross-checked at 1000 m A.G.L using the hybrid single particle Lagrangian integrated trajectory (HYSPLIT) model (https://www.arl.noaa.gov/hysplit/), provided by the US NOAA. All the air mass back trajectories were clustered into four typical types, including NW (20%), S (24%), WNW-NNE (50%), and NE (7%), by HYSPLIT 4.8 and GIS 10.2 software using a hierarchical clustering method.

2.3.2. Geo-Accumulation Index

To compare the concentration of PM2.5-bound metal(loid)s in ambient air with the concentration in the earth’s crust, the geo-accumulation index (IGeo) was used to evaluate contamination levels of metal(loid)s in PM2.5 aerosol samples [26] as follows:
I G e o = log 2 C i Aerosol ( 1.5 × C i Crust ) ,
where C i Aerosol and C i Crust represent the concentration of measured metal(loid) ‘i’ in the PM2.5 aerosol samples and the earth’ crust, respectively. The background metal(loid)s concentrations ( C i Crust ) of Li, Be, V, Cr, Mn, Co, Ni, Cu, Zn, As, Se, Sr, Mo, Ag, Cd, Sn, Sb, Ba, Tl, Pb, Bi, Th, U, and Al of Hubei Province in the earth’s crust were 37.1 ± 10.2, 2.09 ± 0.50, 110 ± 41.9, 86.0 ± 36.9, 712 ± 309, 15.4 ± 5.24, 37.3 ± 15.0, 30.7 ± 14.1, 83.6 ± 36.2, 12.3 ± 7.31, 0.28 ± 0.35, 116 ± 80.7, 1.70 ± 2.52, 0.18 ± 0.11, 0.17 ± 0.19, 2.20 ± 0.71, 1.65 ± 0.61, 542 ± 229, 0.58 ± 0.10, 26.7 ± 7.86, 0.34 ± 0.11, 14.5 ± 1.60, 3.20 ± 1.00 mg kg−1, and 7.19 ± 1.80%, respectively, provided by the Background Values of Soil Elements in China [27]. The constant of 1.5 was to verify the natural fluctuations of a specific metal(loid) in the environment. The IGeo values of PM2.5-bound metal(loid)s were classified as uncontaminated, uncontaminated to moderately contaminated, moderately contaminated, moderately to heavily contaminated, heavily contaminated, heavily to extremely contaminated and extremely contaminated, when IGeo values were less than 0, between 0 and 1, 1 and 2, 2 and 3, 3 and 4, 4 and 5, and greater than 5 [28]. Generally, it may indicate the influence of anthropogenic emissions when IGeo values were higher than 1.

2.3.3. Population Exposure Assessment Model

Population Exposure Dose

Due to long-term expose to PM2.5-bound metal(loid)s in ambient air through ingestion, inhalation, and dermal contact exposure pathways, it may cause potentially adverse population health risks, including non-carcinogenic and carcinogenic risks. A population exposure assessment model was used to evaluate the health risks of PM2.5-bound metal(loid)s in ambient air for humans, including children and adults. According to Technical Guideline for Population Exposure Assessment of Environmental Pollutant (HJ 875-2017, http://www.mee.gov.cn/gkml/hbb/bgg/201711/t20171129_427128.htm), population exposure was defined as average daily population exposure dose (ADPED) of each metal(loid) and was then calculated individually for each metal and exposure pathway, including ingestion ( ADPED Ing ), inhalation ( ADPED Inh ), and dermal contact ( ADPED Der ) as follows:
ADPED Ing = C i Aerosol × R Ing × EF × ED BW × AT × CF ,
ADPED Inh = C i Aerosol × R Inh × EF × ED BW × AT × PEF ,
ADPED Der = C i Aerosol × SA × AF × ABS × EF × ED BW × AT × CF ,
where R Ing and R Inh represent the ingestion and inhalation rate, respectively; ABS represents the dermal absorption factor; EF represents the exposure frequency; ED represents the exposure duration; BW represents the body weight; AT represents the averaging time; PEF represents the particle emission factor; SA represents the skin surface area in contact with air; AF represents the adherence factor for fine particulates to skin; CF represents the conversion factor. The Cr concentration was typically equal to the one-seventh of the total Cr concentration because only Cr(VI) was carcinogenic, while Cr(III) was non-carcinogenic [29]. These variables have been described in previous studies [5,30,31,32] and are listed in Supplementary Table S1.

Population Health Carcinogenic and Non-Carcinogenic Risk

The carcinogenic risks (CRs) due to exposure to As, Cd, Co, Cr(VI), Ni, and Pb were equal to ADPED multiplied by a specific slope factor (SF). The non-carcinogenic risks due to exposure to As, Cd, Cr(III), Co, Cu, Mn, Ni, Zn, Pb, Ag, Al, Ba, Mo, Sb, Sr, U, and V were evaluated using the hazard quotient (HQ), which were calculated by dividing ADPED into a specific reference dose (RfD) [29,30,33]. The total carcinogenic risk (TCR) and total hazard index (THI) indicated the mixed carcinogenic and non-carcinogenic risk due to exposure to an individual metal(loid) and multiple metal(loid)s through three exposure pathways in ambient air, respectively, which were estimated as follows:
CR i = ADED i × SF i ,
HQ i = ADED i RfD i ,
TCR = CR i ,
THI = H Q i ,
Generally, acceptable CR and TCR values ranging from 1 × 10−6 to 1 × 10−4 suggest that the PM2.5-bound metal(loid)s in ambient air does not cause a carcinogenic risk for population health. The carcinogenic risk is categorized as very low, low, moderate, high, and very high for population health when the CR values were less than 1 × 10−6, between 1 × 10−6 and 1 × 10−4, 1 × 10−4 and 1 × 10−3, 1 × 10−3 and 1 × 10−1, and greater than 1 × 10−1 [34]. There was no significant population health non-carcinogenic risk when the HI values were less than 1. Otherwise, there may be non-carcinogenic risks.

3. Results and Discussion

3.1. PM2.5 Mass Concentrations and Pollution Characteristics

Figure 2 presents the time series of mass concentration of PM2.5, WSIIs, CAs, and TE during 13–23 January 2018 at Xiantao site. There was low atmospheric visibility (2.34 km) for 13–23 January and Xiantao experienced an eleven day-long haze episode with an average daily PM2.5 mass concentration of 132 ± 56.0 μg m−3. This level was higher than the levels (123 μg m−3) in January from 2012–2017, and approximately 1.76 and 2.49 times higher than the Grade II of NAAQS and the “normal” levels (53 μg m−3) without pollution events there in January 2018 (http://sthjt.hubei.gov.cn/fbjd/zwgk/jcsjfb/hjkq/), respectively, indicating a heavy PM2.5 pollution during wintertime in JHP, central China. The average daily WSIIs, OC, EC, and TE concentrations were 83.6 ± 46.7, 16.4 ± 5.75, 5.56 ± 1.78, and 7.97 ± 2.36 μg m−3, respectively. The PM2.5 mass concentrations during severe pollution period (277 ± 20.2 μg m−3) was approximately 2–3 times higher than the levels during mild, moderate, and heavy pollution periods (95 ± 10.5, 131 ± 8.52, and 175 ± 27.3 μg m−3, respectively).
When compared to other cities around the world (Table 1), the PM2.5 concentration in Xiantao was comparable to the level in Tianjin, China (124 μg m−3) [11]. The level was much higher compared to most cities in China, including Hefei (81.0 μg m−3) [35], Shanghai (92.9 ± 44.4 μg m−3) [14], Huangshi (98.3 μg m−3) [36], Nanjing (112 ± 26.4 μg m−3) [13], Chengdu (113 ± 60.3 μg m−3) [15], and Beijing (117 μg m−3) [11], and in other countries, such as Iasi (Romania, 23.4 ± 11.7 μg m−3) [37], Zonguldak (Turkey, 37.3 μg m−3) [38], and Riyadh (Saudi Arabia, 71.9 μg m−3) [39]. The level was significantly lower compared to some cities in China, such as Guilin (144 ± 28.5 μg m−3) [40], Wuhan (160 ± 66.3 μg m−3) [17], Xiangyang (169 ± 57.0 μg m−3) [18], Zhengzhou (188 ± 52.4 μg m−3) [16], Shijiazhuang (234 ± 140 μg m−3) [12], and Xi’an (537 ± 146 μg m−3) [10], and in other developing countries, such as Delhi (India, 293 ± 36.7 μg m−3) [41].
The PM2.5 concentration reached the highest value on 20 January during severe pollution period. The phenomenon could have been caused by the air masses from the WNW-NNE direction (first comes from the WNW direction and then turn around from the NNE direction to approach the Xiantao site), accounting for 50% of all the trajectories (Figure 3), which could be found by the highest average concentrations of PM2.5 (153 ± 69.5 μg m−3) and WSIIs (104 ± 55.6 μg m−3) (Supplementary Table S2). These air-masses originated from Xiangyang, passed through Suizhou, Xiaogan and Wuhan, and then approached the Xiantao site with a much lower transport speed during severe pollution period (Supplementary Figure S1), which was not conducive to the horizontal transport for ambient air pollutants. In contrast, the PM2.5 concentration reached the lowest value on 21 January during mild pollution period. This phenomenon could have been partially caused by the air masses from the NE direction, accounting for 7% of all the trajectories (Figure 3), which could be found by the lowest average concentration of PM2.5 (96.7 ± 19.8 μg m−3) (Supplementary Table S2). These air-masses initiated in Anhui Province, then passed through Jiangsu province, and finally entered Hubei Province to approach Xiantao site with a higher transport speed during mild pollution period, which was conducive to the horizontal transport for ambient air pollutants. In addition, rainy and snowy weather (0.30 mm prec.) was observed from midnight of 20 January until the afternoon of 21 January, and thus significantly reduced the PM2.5 concentration through the wet deposition (Supplementary Figure S1).

3.2. Chemical Compositions

3.2.1. WSIIs

Figure 4a presents the WSIIs concentration, Secondary inorganic aerosol (SNA, including SO42−, NO3, and NH4+)/WSIIs and WSIIs/PM2.5 ratios at different pollution levels. SNA accounted for 91.4 ± 4.89% of total WSIIs, ranging from 76.9% to 96.3%. The average WSIIs/PM2.5 ratios were 60.8 ± 9.82%, indicating that WSIIs accounted for the most percentage of PM2.5 mass concentration. It should be noted that the WSIIs concentration increased as the pollution levels increased: mild pollution (54.7 ± 11.0 μg m−3) < moderate pollution (76.3 ± 14.7 μg m−3) < heavy pollution (119 ± 23.5 μg m−3) < severe pollution (208 ± 10.1 μg m−3). Moreover, the SNA/WSIIs and WSIIs/PM2.5 ratios increased as pollution levels increased, which indicated the intensification of secondary formation during the winter haze episode in Xiantao, consistent with our previous simultaneously field sampling campaign conducted in Xiangyang [18].
The acidity of PM2.5 (Supplementary Text S1) could be roughly evaluated through the balance of anions and cations [42]. The average anions/cations ratios during mild, moderate, heavy, and severe pollution periods were 1.04 ± 0.06, 1.00 ± 0.09, 1.05 ± 0.02, and 1.09 ± 0.03, respectively. All these ratios approached 1 and strong positive correlations (r = 0.93, 0.86, 0.99, and 0.82, respectively, p < 0.05) between anions and cations were observed at different pollution levels in Xiantao (Figure 4b), which indicated that almost all the WSIIs were identified and these WSIIs were important alkaline and acidic species in the PM2.5 aerosol [43]. Moreover, all these ratios were also greater than 1, which may be due to a deficiency of H+ in the calculation and/or to NH4+ being converted into the gaseous phase [22]. The anions/cations ratios is also a good indicator for studying the acidity of PM2.5 aerosols [22], which indicated that the PM2.5 aerosols were acidic during the winter haze episode in Xiantao.
The NO3/SO42− ratio has been generally used as an indicator to evaluate the relative contribution of local stationary sources (e.g., coal-fired power plant) and local mobile sources (e.g., vehicle exhaust) to sulfur and nitrogen in the ambient air [44,45]. In this study, the average NO3/SO42− ratios during mild, moderate, heavy, and severe pollution periods were 2.35 ± 0.78, 2.39 ± 1.08, 2.31 ± 0.62, and 1.70 ± 0.15, respectively, which indicated that local mobile sources contributed more to PM2.5 compared to local stationary sources during the winter haze episode. Generally, the sulfur oxidation ratio (SOR) and nitrogen oxidation ratio (NOR) values (Supplementary Text S1) were less than 0.10 for primary pollutants, whereas they were greater than 0.10 when gaseous precursors (SO2 and NO2) were photochemically oxidized in the ambient air [46,47]. The SOR and NOR values during mild, moderate, heavy, and severe pollution periods were 0.38 ± 0.25 and 0.28 ± 0.09, 0.45 ± 0.27 and 0.30 ± 0.09, 0.65 ± 0.11 and 0.47 ± 0.05, and 0.66 ± 0.04 and 0.47 ± 0.05, respectively (Supplementary Figure S2), which indicated that the significantly secondary transformations of SO2 to SO42− and NO2 to NO3 were observed during the winter haze episode.

3.2.2. CAs

The average OC and EC concentrations during the winter haze episode were 16.4 ± 5.74 and 5.56 ± 1.78 μg m−3, respectively, which accounted for 74.3 ± 4.60% and 25.7 ± 4.60% of TC. The OC, EC, and TC concentrations basically increased as the pollution levels increased, except for OC and TC during heavy pollution period. The OC concentration reached the maximum and minimum values at 17:00–22:00 on 15 January (32.5 μg m−3) during the moderate pollution period and at 06:00–11:00 on 22 January (7.21 μg m−3) during mild pollution period, respectively. The EC concentration reached the maximum and minimum values at 17:00–22:00 on 14 January (9.05 μg m−3) during moderate pollution period and at 06:00–11:00 on 22 January (2.66 μg m−3) during mild pollution period, respectively. In addition, the OC/EC ratios ranged from 1.89 to 4.83 with an average OC/EC ratio of 3.01 ± 0.76.
The OC/EC ratios increased as the pollution levels increased from mild to moderate pollution, whereas the ratios decreased as the pollution levels increased from moderate to severe pollution, thus reached the highest levels during moderate pollution period (Figure 5, which indicated that more SOA were formed during moderate pollution period compared to the other three pollution levels. Moreover, the average OC/EC ratios at four pollution levels (3.06 ± 0.85, 3.42 ± 0.87, 2.81 ± 0.39, and 2.42 ± 0.46, respectively) were higher than 2, indicating that secondary organic aerosol (SOA) were formed during the winter haze episode. OC was positively correlated with EC (r = 0.71, p < 0.05), which indicated that there may be a certain same emission source during the winter haze episode. 32.6% of OC/EC ratios ranged from 1.5 to 2.5 (r = 0.95, p < 0.05) (Supplementary Figure S3), it may indicate the contribution of intensive vehicle exhaust to CAs, whereas 67.4% of the OC/EC ratios were higher than 2.5 (r = 0.83, p < 0.05), which might suggest a rapid transformation of secondary organic carbon (SOC) precursors such as VOCs [48].

3.2.3. TE

The TE concentrations during mild, moderate, heavy, and severe pollution periods were 6.64 ± 2.18, 9.26 ± 1.89, 9.59 ± 1.70, and 6.52 ± 2.30 μg m−3, respectively. The proportion of TE to PM2.5 was 6.04 ± 1.79%. The sum of Sb, Al, Si, Zn, Pb, Co, Cr, Mn, Sn, and Cu concentrations accounted for 95.7% of total TE. Figure S4a shows the dynamic variations of TE mass concentrations at the four ambient air quality levels. The Si, Mn, Ni, Cu, Zn, As, Se, Sr, Ag, Cd, Sn, Ba, Pb, and Bi concentrations increased as the pollution levels increased from mild to moderate pollution, whereas these concentrations decreased as the pollution levels from moderate to severe pollution. The Al, Cr, Co, and Sb concentrations increased as the pollution levels increased from mild to heavy pollution, whereas these concentrations decreased as the pollution levels from heavy to severe pollution. The other TE concentrations fluctuated at four pollution levels. The average TE concentration reached the higher value during 18 January to 23 January. The phenomenon could have been partially caused by the air masses from the NE and WNW-NNE directions, accounting for 7% (8.93 ± 1.91 μg m−3) and 50% (8.54 ± 2.31 μg m−3) of all the trajectories (Figure 3 and Supplementary Table S2), respectively. The former air-masses originated from Anhui Province, passed through Shandong Province, and then turn around to approach the Xiantao site with a much higher transport speed during mild pollution period (Supplementary Figure S1) due to long-range transport, which may carry lots of ambient air pollutants to the region. The latter air-masses originated from originated from Xiangyang, passed through Suizhou, Xiaogan, and Wuhan, and then approached the Xiantao site with a much lower transport speed due to combined effect of long-range transport and local emissions (Supplementary Figure S1). Moreover, the rains and snowy weather occurred in the evening of 23 January may increase the hygroscopic growth. In contrast, the average TE concentration reached the lower value during 13 January to 17 January. This phenomenon could have been partially caused by the air masses from the S and NW directions, accounting for 24% (7.59 ± 2.76 μg m−3) and 20% (6.68 ± 2.10 μg m−3) of all the trajectories (Figure 3 and Supplementary Table S2). The former air-masses initiated in Shiyan, then passed through Jingmen, and finally approached the Xiantao site with a lower transport speed (Supplementary Figure S1). The latter airmasses originated from originated from Xiantao, passed through Yueyang and Changsha, and then turn around to approach the Xiantao site with a much higher transport speed (Supplementary Figure S1). Both may be due to the combined effect of long-range transport and local emissions.
The Pb/Cd ratio is generally considered as an indicator to identify the sources of TE [49]. In our study, the Pb/Cd ratios during mild, moderate, heavy, and severe pollution periods were 35.9 ± 33.7, 34.7 ± 9.14, 40.5 ± 7.94, and 44.4 ± 2.89, respectively. These ratios were close to the ratios for anthropogenic emissions (46), indicating that Pb and Cd may be emitted from anthropogenic emissions. Moreover, the V/Ni ratio is commonly used to distinguish industrial emissions (0.7–1.9) and shipping emissions (2.1–3.1) [49,50]. The V/Ni ratios were found to be within the range of 0.7–1.9 during mild and severe pollution periods (1.80 and 1.59, respectively), which indicated that V and Ni may come from industrial emissions during these two periods. This might be mainly affected by the airmasses from the NE and WNW-NNE directions due to combined effect of long-range transport and local emissions as mentioned above. In addition, Pb was positively correlated with Cd (r = 0.96, p < 0.05) during the moderate pollution period, which indicated that the Pb and Cd may have same emission sources during this period. Similar results for V and Ni (r = 0.69, p < 0.05) were observed during heavy pollution period.
Figure 6a and Supplementary Figure S5 present EF and IGeo values for PM2.5-bound metal(loid)s at different pollution levels. Except for Al and Th, the average EF values (Supplementary Text S2) for PM2.5-bound metal(loid)s widely ranged from 3 to 90,717 during the winter haze episode in Xiantao (Figure 6a), which indicated that PM2.5-bound metal(loid)s were affected by natural and anthropogenic activities. The EF values of U, Ba, Sr, Mn, Li, and V were within the range of 2–10, which suggested that these metal(loid)s were mixed affected by crustal sources and anthropogenic emissions. In contrast, the EF values of most metal(loid)s, including Be, Cr, Ni, Cu, Mo, Tl, Zn, Ag, As, Pb, Co, Bi, Cd, Sn, Se, and Sb were significantly higher than 10, indicating the predominate influence from anthropogenic emissions on these metal(loid)s (Supplementary Figure S5). As shown in Figure 6a, the IGeo values of Li, Al, Th, U, Ba, Sr, V, and Mn were lower than 0, indicating that these metal(loid)s were uncontaminated. The IGeo values of Be ranged from 1 to 2, which indicated that there was moderate contamination for Be. The metals of Cr and Ni were moderately to heavily contaminated because the IGeo values of them ranged from 2 to 3. Mo was heavily contaminated because the IGeo values ranged from 3 to 4. Cu and Tl were heavily to extremely contaminated because the IGeo values of them ranged from 4 to 5. What calls for special attention is the extreme contamination for Zn, As, Ag, Pb, Co, Bi, Cd, Sn, Se, and Sb because the IGeo values of these metals were higher than 5, which may have caused population health carcinogenic or non-carcinogenic risk to some extent. It should be noted that EF values were positively correlated with IGeo values at the four pollution levels (r = 0.99, p < 0.01) (Supplementary Figure S4b), which indicated that the dynamic variation trends of EF and IGeo values remained consistent during the winter haze episode.

3.3. Chemical Mass Closure

The chemical mass closure method (Supplementary Text S3) was used to better understand PM2.5 chemical compositions, including mineral dust (MD), trace element oxides (TEO), OM, EC, SNA, Cl and unidentified matters (UM) [51,52,53]. Figure 7a–e presents the chemical mass closure for PM2.5 at different pollution levels at Xiantao site. As shown in Figure 7a, SNA occupied the largest proportion (58.6%) of PM2.5, followed by OM (19.9%), MD (4.31%), EC (4.21%), TEO (4.21%), and Cl (2.76%). The contributions of MD, TEO, OM, EC, and Cl to PM2.5 decreased as the pollution levels increased, whereas the contribution of SNA to PM2.5 increased as pollution levels increased (Figure 7b–e). The results may indicate the contribution variations of emission sources to PM2.5 chemical compositions. On average, only a small proportion of PM2.5 (UM, 6.02%) cannot be identified in this study as well as each pollution level, which indicated that the chemical masses of PM2.5 aerosol samples were balanced during the winter haze episode in Xiantao.

3.4. Population Exposure Assessment

3.4.1. Population Exposure Dose

Table 2 presents ADPED values of PM2.5-bound metal(loid)s through ingestion, inhalation, and dermal contact pathways for children and adults during the winter haze episode. Sb showed the maximum ADPED values of PM2.5-bound metal(loid)s through these three different exposure pathways for children (3.66 × 101) and adults (5.96 × 10−2), whereas U showed the minimum ADPED values of PM2.5-bound metal(loid)s through these three different exposure pathways for children (3.64 × 10−5) and adults (5.92 × 10−6), respectively. Ingestion is found to be the dominant exposure pathway of PM2.5-bound metal(loid)s for children and adults in Xiantao, followed by inhalation and dermal contact, which agreed with previous studies in Kanpur, India [54] and Xiangyang, central China [18]. As shown in Supplementary Table S3, the ADPED values of PM2.5-bound metal(loid)s through ingestion pathway were 1–4 orders of magnitude higher for children and adults at the four pollution levels compared to inhalation and dermal contact pathways. It should be noted that ADPED values of PM2.5-bound metal(loid)s through these three exposure pathways for children were 5.74-fold higher than the values for adults, which indicated that children were more likely to be exposed to PM2.5-bound metal(loid)s than adults during the winter haze episode in Xiantao.

3.4.2. Population Health Risks

Figure 6b and Supplementary Figure S6a–d presents the population health carcinogenic risk and non-carcinogenic risk due to personal exposure to PM2.5-bound metal(loid)s at different pollution levels for children and adults, respectively. The highest carcinogenic risk for children and adults during the winter haze episode was As (1.49 × 10−3 and 2.59 × 10−3), followed by Pb (7.40 × 10−4 and 3.17 × 10−4), Cr(VI) (6.37 × 10−4 and 2.72 × 10−4), Co (1.91 × 10−6 and 6.22 × 10−7), Cd (4.35 × 10−8 and 1.42 × 10−8), and Ni (4.22 × 10−8 and 1.37 × 10−8). The CR values through ingestion, inhalation, and dermal contact exposure pathways widely ranged from 4.22 × 10−8 to 1.49 × 10−3 for children and 1.37 × 10−8 to 2.59 × 10−3 for adults. Moreover, the TCR values were approximately 1.12-fold higher for adults (3.22 × 10−3) than for children (2.87 × 10−3), which indicated that PM2.5-bound metal(loid)s may pose much higher carcinogenic risk for adults compared to children during the winter haze episode in Xiantao (Figure 6b). The CRs value for As is 1~5 orders of magnitude higher than the values for other PM2.5-bound metal(loid)s for children and adults. The TCR value is the sum of CR values for As, Cd, Co, Cr(VI), Ni, and Pb, which mainly depended on the CR value for As. The CR value for As for children (1.49 × 10−3) is lower than the value for adults (2.63 × 10−3), although the CR values for other PM2.5-bound metal(loid)s for children were higher than the values for adults. Therefore, the TCR values were approximately 1.12-fold higher for adults than for children. The TCR values through the three exposure pathways for children and adults during moderate pollution period (1.86 × 10−3 and 3.28 × 10−3) were the highest, followed by mild pollution (1.55 × 10−3 and 2.73 × 10−3), heavy pollution (1.20 × 10−3 and 2.11 × 10−3), and severe pollution periods (3.62 × 10−4 and 6.38 × 10−4) (Supplementary Figure S6a–d). The TCR values through the three exposure pathways for children and adults in Xiantao were lower than the values during the four pollution periods (mild pollution: 1.64 × 10−1, 2.21 × 10−1; moderate pollution: 2.71 × 10−3, 3.77 × 10−3; heavy pollution: 2.27 × 10−3, 2.30 × 10−3; severe pollution: 1.48 × 10−3, 1.65 × 10−3) in Xiangyang [18]. Among all the PM2.5-bound metal(loid)s, the TCR of As, Cr(VI), and Pb for children and adults were at the moderate or high levels, whereas the TCR of Cd, Ni and Co for children and adults were below moderate levels, except for the low carcinogenic risk level of Co for children.
The highest non-carcinogenic risk during the winter haze episode for children and adults was Sb (2.43 × 103 and 2.80 × 103), followed by As (38.7 and 17.1), Cr(III) (13.6 and 15.6), Pb (10.9 and 4.54), and V (1.34 and 1.85). In addition, Co (2.01) and Cd (1.48) may only pose non-carcinogenic risks for children during the winter haze episode. The HQ values through the three exposure pathways widely ranged from 4.13 × 10−3 to 2.43 × 103 for children and 1.45 × 10−3 to 2.80 × 103 for adults. Moreover, the THI values were approximately 1.14-fold higher for adults (2.84 × 103) than for children (2.50 × 103), which indicated that PM2.5-bound metal(loid)s may pose much higher non-carcinogenic risk for adults compared to children during the winter haze episode in Xiantao (Figure 6b). The THI values through the three exposure pathways for children and adults decreased as pollution levels increased: mild pollution (2.42 × 103 and 2.74 × 103) > moderate pollution period (2.07 × 103 and 2.32 × 103) > heavy pollution (1.87 × 103 and 2.12 × 103) > severe pollution (1.13 × 103 and 1.29 × 103) (Supplementary Figure S6a–d). The THI values through the three exposure pathways for children and adults in Xiantao were lower than the values during mild (4.22 × 103, 4.84 × 103) and moderate (2.49 × 103, 2.81 × 103) pollution periods in Xiangyang, whereas the THI values for children and adults were higher than the values during heavy (1.65 × 103, 1.86 × 103) and severe (1.05 × 103, 1.19 × 103) pollution periods in Xiangyang [18]. The total non-carcinogenic risk of Sb for children and adults through the three exposure pathways were 3–6 orders of magnitude higher at the four pollution levels compared to other PM2.5-bound metal(loid)s. It should be noted that As, Cr, and Pb may pose carcinogenic and non-carcinogenic risks for children and adults, whereas Sb and V may only pose non-carcinogenic risks for children and adults. Moreover, the population health risks of PM2.5-bound metal(loid)s may not depend on the pollution levels but depend on the PM2.5-bound metal(loid)s concentrations. Other CAs, such as polycyclic aromatic hydrocarbons (PAHs) could also pose population health risks and needed to be addressed in the future. The government and the public should pay more attention to the population health risks posed by PM2.5-bound metal(loid)s during the winter haze episode in JHP, central China.

4. Conclusions

In this study, a total of 43 pairs of PM2.5 samples were effectively collected to determine the pollution characteristics, chemical compositions, and population health risks during 13–23 January 2018 in Xiantao in JHP, central China. During the sampling period, Xiantao experienced an eleven-day long-lasting haze episode with an average daily PM2.5 mass concentration of 132 ± 56.0 μg m−3. The level was higher than the levels (123 μg m−3) in January from 2012 to 2017, and approximately 1.76 and 2.49 times than the Grade II of NAAQS the “normal” levels (53 μg m−3) without pollution events there in January 2018, respectively, indicating a heavy PM2.5 pollution during wintertime in JHP, central China. The PM2.5 typical chemical compositions for the average WSIIs, OC, EC, and TE concentrations were 83.6 ± 46.7, 16.4 ± 5.75, 5.56 ± 1.78, and 7.97 ± 2.36 μg m−3, respectively. The higher PM2.5 levels during severe pollution period were dominated by the air masses from the WNW-NNE direction, whereas the lower PM2.5 concentrations during other pollution periods were mainly affected by the air masses from the NE, S, and NW directions.
The anions/cations ratios indicated that the PM2.5 aerosols were acidic during the winter haze episode in Xiantao. The NO3/SO42− ratio indicated that the mobile sources contributed more to PM2.5 compared to the stationary sources. The OC/EC ratios indicated a mixed contribution of intensive vehicle exhaust and secondary formations. The Pb/Cd and ratios indicated that Pb and Cd may be emitted from anthropogenic emissions. The V/Ni ratios indicated that V and Ni may come from the industrial emissions during mild and severe pollution periods. The EF values were positively correlated with the IGeo values at the four pollution levels. Ingestion is found to be the dominant exposure pathway of PM2.5-bound metal(loid)s for children and adults in Xiantao, followed by inhalation and dermal contact. As, Cr, and Pb may pose carcinogenic and non-carcinogenic risks for children and adults, whereas Sb and V may only pose non-carcinogenic risks for children and adults. PM2.5-bound metal(loid)s may pose much higher population health risks for adults compared to children during the winter haze episode in Xiantao. Moreover, the population health risks of PM2.5-bound metal(loid)s may not depend on the pollution levels but depend on the PM2.5-bound metal(loid)s concentrations. Other CAs, such as polycyclic aromatic hydrocarbons (PAHs) could also pose population health risks and need to be addressed in the future. The government and the public should pay more attention to the population health risks posed by PM2.5-bound metal(loid)s during a winter haze episode in JHP, central China.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4433/11/9/954/s1, Supplementary Materials: Text S1: PM2.5 Acidity, Sulfur, and Nitrogen Oxidation Ratios, Text S2: Enrichment Factor, Text S3: Chemical Mass Closure Method, Table S1: Summary of parameter values for population health risk assessment, Table S2: The mass concentrations (Mean ± Standard deviation, μg m−3) of PM2.5 and the corresponding species for the air mass back trajectory clusters during 13–23 January 2018 in Xiantao, Table S3: The ADPED (mg kg−1 day−1) of PM2.5-bound metal(loid)s through ingestion, inhalation, and dermal contact pathways for children and adults at different pollution levels, Figure S1: Time series of hourly meteorological parameter, including wind direction (WD), wind speed (WS), temperature (Temp.) relative humidity (RH), precipitation (Prec.), and visibility (Vis.) at Xiantao site in January 2018, Figure S2: Time series of trace gases (SO2 and NO2), SOR, and NOR during the sampling time at Xiantao site, Figure S3: The OC/EC ratios during the winter haze episode at Xiantao site, Figure S4: The (a) TE concentration and (b) correlations between EF and IGeo at different pollution levels, Figure S3: The contamination levels for PM2.5-bound metal(loid)s at different pollution levels: (a) mild pollution, (b) moderate pollution, (c) heavy pollution, and (d) severe pollution, Figure S4: The population health carcinogenic risk and non-carcinogenic risk due to personal exposure to PM2.5-bound metal(loid)s for children and adults, respectively, at different pollution levels: (a) mild pollution, (b) moderate pollution, (c) heavy pollution, and (d) severe pollution.

Author Contributions

Conceptualization, data curation, and formal analysis, H.J.; methodology, H.J. and H.X.; investigation, X.L., H.X., D.C., P.Z., H.J., and X.Y.; writing—original draft, H.J.; validation and writing—review and editing, H.C.; supervision, resources, project administration, and funding acquisition, H.C. and Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R & D Program of China (Nos. 2017YFC0212603 and 2019YFB2102900) and the Natural Science Foundation of China (Grant No. 41673102).

Acknowledgments

We thank the engineer X.G. (Xinli Guo) from Wuhan Tianhong Instrument Co., Ltd. for the maintenance of sampling instrument. We are also grateful to the staff from Xiantao Municipal Ecological Environment Bureau for their support to the sampling work and the NOAA’s Air Resources Laboratory for providing the HYSPLIT model and the READY website (https://www.arl.noaa.gov/hysplit/).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sun, Y.; Zhuang, G.; Tang, A.; Wang, Y.; An, Z. Chemical characteristics of PM2.5 and PM10 in haze-fog episodes in Beijing. Environ. Sci. Technol. 2006, 40, 3148–3155. [Google Scholar] [CrossRef] [PubMed]
  2. Fiore, A.M.; Naik, V.; Spracklen, D.V.; Steiner, A.; Unger, N.; Prather, M.; Bergmann, D.; Cameron-Smith, P.J.; Cionni, I.; Collins, W.J.; et al. Cationic polymers and their therapeutic potential. Chem. Soc. Rev. 2012, 41, 6663–6683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Hamra, G.B.; Guha, N.; Cohen, A.; Laden, F.; Raaschou-Nielsen, O.; Samet, J.M.; Vineis, P.; Forastiere, F.; Saldiva, P.; Yorifuji, T.; et al. Outdoor particulate matter exposure and lung cancer: A systematic review and meta-analysis. Environ. Health Perspect. 2014, 122, 906–911. [Google Scholar] [CrossRef] [Green Version]
  4. Loomis, D.; Grosse, Y.; Lauby-Secretan, B.; El Ghissassi, F.; Bouvard, V.; Benbrahim-Tallaa, L.; Guha, N.; Baan, R.; Mattock, H.; Straif, K. The carcinogenicity of outdoor air pollution. Lancet Oncol. 2013, 14, 1262–1263. [Google Scholar] [CrossRef]
  5. Hu, X.; Zhang, Y.; Ding, Z.; Wang, T.; Lian, H.; Sun, Y.; Wu, J. Bioaccessibility and health risk of arsenic and heavy metals (Cd, Co, Cr, Cu, Ni, Pb, Zn and Mn) in TSP and PM2.5 in Nanjing, China. Atmos. Environ. 2012, 57, 146–152. [Google Scholar] [CrossRef]
  6. Zhang, R.; Wang, G.; Guo, S.; Zamora, M.L.; Ying, Q.; Lin, Y.; Wang, W.; Hu, M.; Wang, Y. Formation of Urban Fine Particulate Matter. Chem. Rev. 2015, 115, 3803–3855. [Google Scholar] [CrossRef]
  7. Cakmak, S.; Dales, R.; Kauri, L.M.; Mahmud, M.; Van Ryswyk, K.; Vanos, J.; Liu, L.; Kumarathasan, P.; Thomson, E.; Vincent, R.; et al. Metal composition of fine particulate air pollution and acute changes in cardiorespiratory physiology. Environ. Pollut. 2014, 189, 208–214. [Google Scholar] [CrossRef] [Green Version]
  8. Gavett, S.H.; Haykal-Coates, N.; Copeland, L.B.; Heinrich, J.; Gilmour, M.I. Metal composition of ambient PM2.5 influences severity of allergic airways disease in mice. Environ. Health Perspect. 2003, 111, 1471–1477. [Google Scholar] [CrossRef] [Green Version]
  9. Delfino, R.J.; Staimer, N.; Tjoa, T.; Gillen, D.L.; Polidori, A.; Arhami, M.; Kleinman, M.T.; Vaziri, N.D.; Longhurst, J.; Sioutas, C. Air pollution exposures and circulating biomarkers of effect in a susceptible population: Clues to potential causal component mixtures and mechanisms. Environ. Health Perspect. 2009, 117, 1232–1238. [Google Scholar] [CrossRef]
  10. Elser, M.; Huang, R.; Wolf, R.; Slowik, J.G.; Wang, Q.; Canonaco, F.; Li, G.; Bozzetti, C.; Daellenbach, K.R.; Huang, Y.; et al. New insights into PM2.5 chemical composition and sources in two major cities in China during extreme haze events using aerosol mass spectrometry. Atmos. Chem. Phys. 2016, 3207–3225. [Google Scholar] [CrossRef] [Green Version]
  11. Shen, R.; Liu, Z.; Chen, X.; Wang, Y.; Wang, L.; Liu, Y.; Li, X. Atmospheric levels, variations, sources and health risk of PM2.5-bound polycyclic aromatic hydrocarbons during winter over the North China Plain. Sci. Total Environ. 2019, 655, 581–590. [Google Scholar] [CrossRef] [PubMed]
  12. Xie, Y.; Liu, Z.; Wen, T.; Huang, X.; Liu, J.; Tang, G.; Yang, Y.; Li, X.; Shen, R.; Hu, B.; et al. Characteristics of chemical composition and seasonal variations of PM2.5 in Shijiazhuang, China: Impact of primary emissions and secondary formation. Sci. Total Environ. 2019, 677, 215–229. [Google Scholar] [CrossRef] [PubMed]
  13. Li, S.; Chang, M.; Li, H.; Cui, X.; Ma, L.Q. Chemosphere Chemical compositions and source apportionment of PM2.5 during clear and hazy days: Seasonal changes and impacts of Youth Olympic Games. Chemosphere 2020, 256, 127163. [Google Scholar] [CrossRef] [PubMed]
  14. Wei, N.; Xu, Z.; Wang, G.; Liu, W.; Zhouga, D.; Xiao, D.; Yao, J. Source apportionment of carbonaceous aerosols during haze days in Shanghai based on dual carbon isotopes. J. Radioanal. Nucl. Chem. 2019. [Google Scholar] [CrossRef]
  15. Qu, Y.; Gao, T.; Yang, C. Elemental characterization and source identification of the near-road PM2.5 using EDXRF in Chengdu, China. X-Ray Spectrom. 2019, 48, 232–241. [Google Scholar] [CrossRef]
  16. Wang, S.; Yin, S.; Zhang, R.; Yang, L.; Zhao, Q.; Zhang, L.; Yan, Q.; Jiang, N.; Tang, X. Insight into the formation of secondary inorganic aerosol based on high-time-resolution data during haze episodes and snowfall periods in Zhengzhou, China. Sci. Total Environ. 2019, 660, 47–56. [Google Scholar] [CrossRef]
  17. Liu, J.; Li, J.; Vonwiller, M.; Liu, D.; Cheng, H.; Shen, K.; Salazar, G.; Agrios, K.; Zhang, Y.; He, Q.; et al. The importance of non-fossil sources in carbonaceous aerosols in a megacity of central China during the 2013 winter haze episode: A source apportionment constrained by radiocarbon and organic tracers. Atmos. Environ. 2016, 144, 60–68. [Google Scholar] [CrossRef]
  18. Jiang, H.; Xiao, H.; Song, H.; Liu, J.; Wang, T.; Cheng, H.; Wang, Z. A long-lasting winter haze episode in Xiangyang, Central China: Pollution characteristics, chemical composition, and health risk assessment. Aerosol Air Qual. Res. 2020, 20. [Google Scholar] [CrossRef]
  19. Qu, C.; Xing, X.; Albanese, S.; Doherty, A.; Huang, H.; Lima, A.; Qi, S.; De Vivo, B. Spatial and seasonal variations of atmospheric organochlorine pesticides along the plain-mountain transect in central China: Regional source vs. long-range transport and air-soil exchange. Atmos. Environ. 2015, 122, 31–40. [Google Scholar] [CrossRef]
  20. Liu, R.; Ma, Z.; Liu, Y.; Shao, Y.; Zhao, W.; Bi, J. Spatiotemporal distributions of surface ozone levels in China from 2005 to 2017: A machine learning approach. Environ. Int. 2020, 142, 105823. [Google Scholar] [CrossRef]
  21. Liu, Q.; Baumgartner, J.; Zhang, Y.; Schauer, J.J. Source apportionment of Beijing air pollution during a severe winter haze event and associated pro-inflammatory responses in lung epithelial cells. Atmos. Environ. 2016, 126, 28–35. [Google Scholar] [CrossRef]
  22. Cheng, H.; Gong, W.; Wang, Z.; Zhang, F.; Wang, X.; Lvy, X.; Liu, J.; Fu, X.; Zhang, G. Ionic composition of submicron particles (PM1.0) during the long-lasting haze period in January 2013 in Wuhan, central China.pdf. J. Environ. Sci. 2014, 26, 810–817. [Google Scholar] [CrossRef]
  23. Zhang, F.; Cheng, H.R.; Wang, Z.W.; Lv, X.P.; Zhu, Z.M.; Zhang, G.; Wang, X.M. Fine particles (PM2.5) at a CAWNET background site in Central China: Chemical compositions, seasonal variations and regional pollution events. Atmos. Environ. 2014, 86, 193–202. [Google Scholar] [CrossRef]
  24. Zhang, F.; Wang, Z.W.; Cheng, H.R.; Lv, X.P.; Gong, W.; Wang, X.M.; Zhang, G. Seasonal variations and chemical characteristics of PM2.5 in Wuhan, Central China. Sci. Total Environ. 2015, 518, 97–105. [Google Scholar] [CrossRef] [PubMed]
  25. Lyu, X.-P.; Wang, Z.-W.; Cheng, H.-R.; Zhang, F.; Zhang, G.; Wang, X.-M.; Ling, Z.-H.; Wang, N. Chemical characteristics of submicron particulates (PM1.0) in Wuhan, Central China. Atmos. Res. 2015, 161, 169–178. [Google Scholar] [CrossRef]
  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] [Green Version]
  27. State Environmental Protection Bureau. Background Values of Soil Elements in China, 1st ed.; China Environmental Monitoring Station, Ed.; China Environmental Science Press: Beijing, China, 1990. [Google Scholar]
  28. Li, N.; Han, W.; Wei, X.; Shen, M.; Sun, S. Chemical characteristics and human health assessment of PM1 during the Chinese Spring Festival in Changchun, Northeast China. Atmos. Pollut. Res. 2019, 10, 1823–1831. [Google Scholar] [CrossRef]
  29. Massey, D.D.; Kulshrestha, A.; Taneja, A. Particulate matter concentrations and their related metal toxicity in rural residential environment of semi-arid region of India. Atmos. Environ. 2013, 67, 278–286. [Google Scholar] [CrossRef]
  30. US Environmental Protection Agency. Exposure Factors Handbook: 2011 Edition; US Environmental Protection Agency: Washington, DC, USA, 2011; EPA/600/R-090/052F. [Google Scholar]
  31. Wang, N.; Zhao, X.; Wang, J.; Yin, B.; Geng, C.; Niu, D. Chemical Composition of PM2.5 and Its Impact on Inhalation Health Risk Evaluation in a City with Light Industry in Central China. Atmosphere 2020, 11, 340. [Google Scholar] [CrossRef] [Green Version]
  32. Griselda, J.; Mar, R.; Antonio, A.; Guzman, E.; Guarnaccia, C.; Mart, S.; Lara, C.; Rangel, M.; Hern, G.; Lozada, C.; et al. Trace Metal Content and Health Risk Assessment of PM10 in an Urban Environment of León, Mexico. Atmosphere 2019, 10, 573. [Google Scholar] [CrossRef] [Green Version]
  33. International Agency for Research on Cancer. Agents Classified by the IARC Monographs; World Health Organization: Geneva, Switzerland, 2020; Volume 1-116. [Google Scholar]
  34. Roy, D.; Singh, G.; Seo, Y.C. Carcinogenic and non-carcinogenic risks from PM10- and PM2.5-Bound metals in a critically polluted coal mining area. Atmos. Pollut. Res. 2019, 10, 1964–1975. [Google Scholar] [CrossRef]
  35. Xue, H.; Liu, G.; Zhang, H.; Hu, R.; Wang, X. Similarities and differences in PM10 and PM2.5 concentrations, chemical compositions and sources in Hefei City, China. Chemosphere 2019, 220, 760–765. [Google Scholar] [CrossRef]
  36. Mao, Y.; Hu, T.; Shi, M.; Cheng, C.; Liu, W.; Zhang, J. PM2.5-bound PAHs during a winter haze episode in a typical mining city, central China: Characteristics, influencing parameters, and sources. Atmos. Pollut. Res. 2020, 11, 131–140. [Google Scholar] [CrossRef]
  37. Giorgiana Galon-Negru, A.; Iulian Olariu, R.; Arsene, C. Chemical characteristics of size-resolved atmospheric aerosols in Iasi, north-eastern Romania: Nitrogen-containing inorganic compounds control aerosol chemistry in the area. Atmos. Chem. Phys. 2018, 18, 5879–5904. [Google Scholar] [CrossRef] [Green Version]
  38. Akyüz, M.; Çabuk, H. Meteorological variations of PM2.5/PM10 concentrations and particle-associated polycyclic aromatic hydrocarbons in the atmospheric environment of Zonguldak, Turkey. J. Hazard. Mater. 2009, 170, 13–21. [Google Scholar] [CrossRef] [PubMed]
  39. Modaihsh, A.S.; Al-Barakah, F.N.; Nadeem, M.E.A.; Mahjoub, M.O. Spatial and Temporal Variations of the Particulate Matter in Riyadh City, Saudi Arabia. J. Environ. Prot. 2015, 6, 1293–1307. [Google Scholar] [CrossRef] [Green Version]
  40. Zhong, S.; Zhang, L.; Jiang, X.; Gao, P. Comparison of chemical composition and airborne bacterial community structure in PM2.5 during haze and non-haze days in the winter in Guilin, China. Sci. Total Environ. 2019, 655, 202–210. [Google Scholar] [CrossRef] [PubMed]
  41. Khanna, I.; Khare, M.; Gargava, P.; Khan, A.A. Effect of PM2.5 chemical constituents on atmospheric visibility impairment.pdf. J. Air Waste Manag. Assoc. 2018, 68, 430–437. [Google Scholar] [CrossRef] [Green Version]
  42. Chow, J.C.; Fujita, E.M.; Watson, J.G.; Lu, Z.; Lawson, D.R. Evaluation of Filter-based Aerosol Measurements During the 1987 Southern California Air Quality Study. Environ. Monit. Assess. 1994, 87, 49–80. [Google Scholar] [CrossRef]
  43. Wang, Y.; Zhuang, G.; Tang, A.; Yuan, H.; Sun, Y.; Chen, S.; Zheng, A. The ion chemistry and the source of PM2.5 aerosol in Beijing. Atmos. Environ. 2005, 39, 3771–3784. [Google Scholar] [CrossRef]
  44. Arimoto, R.; Duce, R.A.; Savoie, D.L.; Prospero, J.M.; Talbot, R.; Cullen, J.D.; Tomza, U.; Lewis, N.F.; Ray, B.J. Relationships among aerosol constituents from Asia and the North pacific during PEM-West A. J. Geophys. Res. 1996, 101, 2011–2023. [Google Scholar] [CrossRef]
  45. Yin, L.; Niu, Z.; Chen, X.; Chen, J.; Zhang, F.; Xu, L. Characteristics of water-soluble inorganic ions in PM2.5 and PM2.5–10 in the coastal urban agglomeration along the Western Taiwan Strait Region, China. Environ. Sci. Pollut. Res. 2014, 21, 5141–5156. [Google Scholar] [CrossRef] [PubMed]
  46. Yang, S.; Ma, Y.L.; Duan, F.K.; He, K.B.; Wang, L.T.; Wei, Z.; Zhu, L.D.; Ma, T.; Li, H.; Ye, S.Q. Characteristics and formation of typical winter haze in Handan, one of the most polluted cities in China. Sci. Total Environ. 2018, 613–614, 1367–1375. [Google Scholar] [CrossRef] [PubMed]
  47. Feng, J.; Yu, H.; Mi, K.; Su, X.; Li, Y.; Li, Q.; Sun, J. One year study of PM2.5 in Xinxiang city, North China: Seasonal characteristics, climate impact and source. Ecotoxicol. Environ. Saf. 2018, 154, 75–83. [Google Scholar] [CrossRef]
  48. Pachauri, T.; Satsangi, A.; Singla, V.; Lakhani, A.; Maharaj Kumari, K. Characteristics and sources of carbonaceous aerosols in PM2.5 during wintertime in Agra, India. Aerosol Air Qual. Res. 2013, 13, 977–991. [Google Scholar] [CrossRef] [Green Version]
  49. Bressi, M.; Sciare, J.; Ghersi, V.; Mihalopoulos, N.; Petit, J.-E.; Nicolas, J.B.; Moukhtar, S.; Rosso, A.; Féron, A.; Bonnaire, N.; et al. Sources and geographical origins of fine aerosols in Paris (France). Atmos. Chem. Phys. 2014, 14, 8813–8839. [Google Scholar] [CrossRef] [Green Version]
  50. Mamane, Y.; Perrino, C.; Yossef, O.; Catrambone, M. Source characterization of fine and coarse particles at the East Mediterranean coast. Atmos. Environ. 2008, 42, 6114–6130. [Google Scholar] [CrossRef]
  51. Zhang, R.; Jing, J.; Tao, J.; Hsu, S.-C.; Wang, G.; Cao, J.; Lee, C.S.L.; Zhu, L.; Chen, Z.; Zhao, Y.; et al. Chemical characterization and source apportionment of PM2.5 in Beijing: Seasonal perspective. Atmos. Chem. Phys. 2013, 13, 7053–7074. [Google Scholar] [CrossRef] [Green Version]
  52. Kong, S.F.; Li, L.; Li, X.X.; Yin, Y.; Chen, K.; Liu, D.T.; Yuan, L.; Zhang, Y.J.; Shan, Y.P.; Ji, Y.Q. The impacts of firework burning at the Chinese Spring Festival on air quality: Insights of tracers, source evolution and aging processes. Atmos. Chem. Phys. 2015, 15, 2167–2184. [Google Scholar] [CrossRef] [Green Version]
  53. Zheng, H.; Kong, S.; Wu, F.; Cheng, Y.; Niu, Z.; Zheng, S.; Yang, G.; Yao, L.; Yan, Q.; Wu, J.; et al. Intra-regional transport of black carbon between the south edge of the North China Plain and central China during winter haze episodes. Atmos. Chem. Phys. 2019, 19, 4499–4516. [Google Scholar] [CrossRef] [Green Version]
  54. Izhar, S.; Goel, A.; Chakraborty, A.; Gupta, T. Annual trends in occurrence of submicron particles in ambient air and health risk posed by particle bound metals. Chemosphere 2016, 146, 582–590. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The geographical location of Xiantao sampling site (© Google Earth).
Figure 1. The geographical location of Xiantao sampling site (© Google Earth).
Atmosphere 11 00954 g001
Figure 2. Time series of mass concentration of (a) PM2.5 and WSIIs, (b) CAs, and (c) TE during 13–23 January 2018 at Xiantao site.
Figure 2. Time series of mass concentration of (a) PM2.5 and WSIIs, (b) CAs, and (c) TE during 13–23 January 2018 at Xiantao site.
Atmosphere 11 00954 g002
Figure 3. Four simulated clusters for 3 days (72 h) air masses back trajectories every four times a day at 500 m A.G.L using HYSPLIT model (https://www.arl.noaa.gov/hysplit/) during 13–23 January 2018 at Xiantao site.
Figure 3. Four simulated clusters for 3 days (72 h) air masses back trajectories every four times a day at 500 m A.G.L using HYSPLIT model (https://www.arl.noaa.gov/hysplit/) during 13–23 January 2018 at Xiantao site.
Atmosphere 11 00954 g003
Figure 4. Dynamic variations of (a) WSIIs, SNA/WSIIs and WSIIs/PM2.5 ratios, and (b) acidity of PM2.5 at different pollution levels at Xiantao site.
Figure 4. Dynamic variations of (a) WSIIs, SNA/WSIIs and WSIIs/PM2.5 ratios, and (b) acidity of PM2.5 at different pollution levels at Xiantao site.
Atmosphere 11 00954 g004
Figure 5. Dynamic variations of CAs (OC, EC, and TC) and OC/EC ratios at different pollution levels at Xiantao site.
Figure 5. Dynamic variations of CAs (OC, EC, and TC) and OC/EC ratios at different pollution levels at Xiantao site.
Atmosphere 11 00954 g005
Figure 6. The (a) EF and IGeo, and (b) population health risks for PM2.5-bound metal(loid)s during the winter haze episode.
Figure 6. The (a) EF and IGeo, and (b) population health risks for PM2.5-bound metal(loid)s during the winter haze episode.
Atmosphere 11 00954 g006
Figure 7. Chemical mass closure for PM2.5 at Xiantao site at different pollution levels: (a) winter haze episode, (b) mild pollution, (c) moderate pollution, (d) heavy pollution, and (e) severe pollution.
Figure 7. Chemical mass closure for PM2.5 at Xiantao site at different pollution levels: (a) winter haze episode, (b) mild pollution, (c) moderate pollution, (d) heavy pollution, and (e) severe pollution.
Atmosphere 11 00954 g007
Table 1. Reported PM2.5 mass concentrations (Mean ± Standard deviation, μg m−3) in Xiantao compared to other cities around the world.
Table 1. Reported PM2.5 mass concentrations (Mean ± Standard deviation, μg m−3) in Xiantao compared to other cities around the world.
Station (Country)Site TypeField Sampling CampaignMean ± Standard DeviationReferences
Xiantao (China)Urban13–23 January 2018132 ± 56.0This Study
Zonguldak (Turkey)UrbanWinter, 200737.3[38]
Xi’an (China)Urban17–26 December 2013537 ± 146[10]
Riyadh (Saudi Arabia)UrbanJanuary 201371.9[39]
Wuhan (China)Urban9 January 2013–6 February 2013160 ± 66.3[17]
Delhi (India)Near-roadDecember 2013–January 2014293 ± 36.7[41]
Nanjing (China)UrbanWinter, 2014112 ± 26.4[13]
Beijing (China)Urban29 December 2014–27 January 2015117[11]
Tianjin (China)Urban124
Iasi (Romania)UrbanJanuary 201623.4 ± 11.7[37]
Shijiazhuang (China)Urban15 December 2015–14 January 2016234 ± 140[12]
Hefei (China)IndustrialOctober 2016–January201781.0[35]
Guilin (China)Urban10–15 December 2016–4–7 January 2017144 ± 28.5[40]
Shanghai (China)Suburban2018 haze days92.9 ± 44.4[14]
Chengdu (China)Near-roadWinter, 2018113 ± 60.3[15]
Huangshi (China)Industrial13–25 January 201898.3[36]
Zhengzhou (China)Urban11–23 January 2018188 ± 52.4[16]
Xiangyang (China)Urban13–24 January 2018169 ± 57.0[18]
Table 2. The ADPED (mg kg−1 day−1) of PM2.5-bound metal(loid)s through ingestion, inhalation, and dermal contact pathways for children and adults during the winter haze episode.
Table 2. The ADPED (mg kg−1 day−1) of PM2.5-bound metal(loid)s through ingestion, inhalation, and dermal contact pathways for children and adults during the winter haze episode.
PM2.5-Bound Metal(loid)sChildrenAdults
ADPEDIng *ADPEDInh *ADPEDDer *ADPEDADPEDIng *ADPEDInh *ADPEDDer *ADPED
Sb3.54 × 10−12.60 × 10−51.24 × 10−23.66 × 10−13.79 × 10−22.12 × 10−62.16 × 10−25.96 × 10−2
Al2.24 × 10−11.65 × 10−57.85 × 10−32.32 × 10−12.40 × 10−21.34 × 10−61.37 × 10−23.78 × 10−2
Zn6.45 × 10−24.74 × 10−62.26 × 10−36.67 × 10−22.29 × 10−31.28 × 10−71.30 × 10−33.59 × 10−3
Pb3.35 × 10−22.46 × 10−61.17 × 10−33.47 × 10−24.44 × 10−32.48 × 10−72.53 × 10−36.97 × 10−3
Co3.22 × 10−22.37 × 10−61.13 × 10−33.33 × 10−23.86 × 10−32.15 × 10−72.20 × 10−36.05 × 10−3
Cr1.74 × 10−21.28 × 10−66.09 × 10−41.80 × 10−22.72 × 10−31.52 × 10−71.55 × 10−34.27 × 10−3
Mn1.41 × 10−21.03 × 10−64.92 × 10−41.46 × 10−21.51 × 10−38.42 × 10−88.59 × 10−42.36 × 10−3
Cu1.16 × 10−28.53 × 10−74.06 × 10−41.20 × 10−21.24 × 10−36.95 × 10−87.09 × 10−41.95 × 10−3
As1.00 × 10−27.38 × 10−71.05 × 10−31.11 × 10−21.33 × 10−37.43 × 10−82.27 × 10−33.60 × 10−3
Ni8.66 × 10−36.37 × 10−73.03 × 10−48.96 × 10−31.15 × 10−36.41 × 10−86.54 × 10−41.80 × 10−3
Ba6.95 × 10−35.11 × 10−72.43 × 10−47.20 × 10−33.03 × 10−41.69 × 10−81.72 × 10−44.75 × 10−4
Sr2.11 × 10−31.55 × 10−77.38 × 10−52.18 × 10−32.26 × 10−41.26 × 10−81.29 × 10−43.55 × 10−4
V2.08 × 10−31.53 × 10−77.29 × 10−52.16 × 10−32.23 × 10−41.25 × 10−81.27 × 10−43.51 × 10−4
Cd1.19 × 10−38.75 × 10−84.17 × 10−61.19 × 10−31.58 × 10−48.82 × 10−98.99 × 10−61.67 × 10−4
Mo9.30 × 10−46.84 × 10−83.26 × 10−59.63 × 10−49.97 × 10−55.57 × 10−95.68 × 10−51.56 × 10−4
Ag1.47 × 10−41.08 × 10−85.13 × 10−61.52 × 10−41.57 × 10−58.78 × 10−108.96 × 10−62.47 × 10−5
U3.52 × 10−52.59 × 10−91.23 × 10−63.64 × 10−53.77 × 10−62.11 × 10−102.15 × 10−65.92 × 10−6
* ADPEDIng, ADPEDInh, and ADPEDDer represent the average daily population exposure dose due to exposure to PM2.5-bound metal(loid)s through ingestion, inhalation, and dermal contact pathways, respectively.

Share and Cite

MDPI and ACS Style

Jiang, H.; Lei, X.; Xiao, H.; Chen, D.; Zeng, P.; Yang, X.; Wang, Z.; Cheng, H. Pollution Characteristics, Chemical Compositions, and Population Health Risks during the 2018 Winter Haze Episode in Jianghan Plain, Central China. Atmosphere 2020, 11, 954. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos11090954

AMA Style

Jiang H, Lei X, Xiao H, Chen D, Zeng P, Yang X, Wang Z, Cheng H. Pollution Characteristics, Chemical Compositions, and Population Health Risks during the 2018 Winter Haze Episode in Jianghan Plain, Central China. Atmosphere. 2020; 11(9):954. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos11090954

Chicago/Turabian Style

Jiang, Huimeng, Xiaowei Lei, Han Xiao, Danhong Chen, Pei Zeng, Xingyun Yang, Zuwu Wang, and Hairong Cheng. 2020. "Pollution Characteristics, Chemical Compositions, and Population Health Risks during the 2018 Winter Haze Episode in Jianghan Plain, Central China" Atmosphere 11, no. 9: 954. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos11090954

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