Sources and Seasonal Variance of Ambient Volatile Organic Compounds in the Typical Industrial City of Changzhi, Northern China

: Volatile organic compounds (VOCs) emitted from industrial processes, which are major emission sources of air pollutants, could cause signiﬁcant impacts on air quality. However, studies on the comprehensive analysis from sources contributing to the health risk perspective regarding ambient VOCs in industrial cities are limited. In this study, VOC samples were collected from 15 April 2018 to 19 October 2018 in Changzhi, a typical industrial city in northern China, and a total of 57 VOCs were measured for analysis. The average VOC concentrations were 54.4 µ g · m − 3 , with the highest concentrations in autumn (58.4 µ g · m − 3 ). Ambient VOCs in spring, summer and autumn were all dominated by alkanes (66.8%), with contributions of 70.3%, 66.3% and 63.8%, respectively. The top ﬁve concentrations of total VOCs were isopentane (19.0%), ethane (9.5%), n-butane (8.1%), benzene (7.9%) and propane (5.2%), indicating that vehicle exhaust and coal combustion are the main sources of VOCs. Source apportionment by principal component analysis showed that vehicle exhaust (27.5%) and coal combustion (23.5%) were the main sources of VOCs in Changzhi, followed by industrial production (17.4%), solvent evaporation (13.5%), liqueﬁed petroleum gas/natural gas leaking (9.5%), and biogenic emissions (8.7%). Sources of coal combustion and vehicle exhaust contributed more VOCs than industrial production. The carcinogenic risks of benzene (3.4 × 10 − 5 ) and ethylbenzene × 10 − 6 ) were higher than the limit levels (1 × 10 − 6 ). Coal combustion contributed most (25.3%) to the carcinogenic risks because of its high VOC emissions. In an industrial city such as Changzhi, vehicle exhaust and coal combustion have become major sources of ambient air VOCs owing to the increasingly stringent industrial standards. Therefore, VOCs from vehicle exhaust and coal combustion also need to take into account mitigation measures for VOCs from the perspective of source contribution to health risk.


Introduction
The concentration of particulate matter (PM) in China has been significantly reduced in recent years owing to the implementation of strict clean air policies. For example, PM 2.5 and PM 10 concentrations decreased by 30-33% between 2015 and 2019 in the Beijing-Tianjin-Hebei region, and by 7-13% in the Pearl River Delta [1,2]. Conversely, ozone air pollution has continued to rise, and ozone has become a priority pollutant in many Chinese cities during the summer [3,4]. Ozone pollution is addressed in the ambient air quality standards issued by the Chinese Ministry of Environmental Protection in 2016, and it is the only air pollutant whose concentration has increased over the past few years [5][6][7][8][9][10]. High ozone concentrations deteriorate air quality and cause a range of human health and ecosystem problems [11,12]. Volatile organic compounds (VOCs) are important precursors to ozone production in the near-surface atmosphere, and high concentrations of VOCs often lead to high ozone concentrations. In addition, VOCs also pose both short-and long-term negative impacts to human health [13][14][15]. Recent studies have shown that increased O 3 concentrations also lead to reduced carbon uptake by plants [16]. Therefore, the identification and quantification of potential VOC sources can provide a basis for decision makers to develop control strategies to mitigate VOC emissions and O 3 pollution.
Atmospheric VOCs can be emitted from natural and anthropogenic sources. Plants are the main natural sources of VOCs, including isoprene, monoterpenes, sesquiterpenes, etc. [17]. Studies have shown that 34.27 (±2.06) Tg of VOCs are emitted annually from plant sources in China, including 15.94 (±1.12) Tg of isoprene [18]. Studies for Athens and Patras have shown that VOCs are mainly caused by biogenic and traffic emissions in summer while biomass burning is predominant in winter [19]. Studies in the city of Barykessy, Turkey showed that in a drinking water treatment plant, industrial solvent emissions from pesticide and insecticide applications in agriculture and emissions from traffic were the main sources of VOCs in the atmosphere [20]. Anthropogenic sources primarily include vehicle exhaust, solvent evaporation, industrial production, coal combustion, and liquefied petroleum gas/natural gas (LPG/NG) leaking [21]. The characteristics and sources of VOCs differ according to energy, economy, and industrial structure of a region. In Shanghai, vehicle exhaust is the main source of VOCs, accounting for 34% of ambient VOCs [22]. However, in Nanjing, fossil fuel/biomass/biofuel combustion accounted for 28% of the total VOC concentrations [23]. In rural eastern China, 71% of total VOC emissions are attributed to a combination of vehicle exhaust and biofuel combustion [24]. However, industrial production is also an important source of environmental VOCs, and the high concentration and intensity of emissions from industrial zones can significantly affect surrounding areas [25]. The study of Chang et al. on source identification of environmental VOCs by principal component analysis showed that PU coating, chemical packaging, and lithographic printing industries were the main sources of VOC emissions in industrial parks in Taiwan [26]. Previous studies have also investigated the influence of industrial production on ambient VOCs in industrial cities. For example, studies in Bursa, one of the most industrialized cities in Turkey, have shown that vehicle exhaust (gasoline and diesel engine vehicles) and industrial emissions are the main sources of VOCs [27]. VOC emissions from industrial production have been widely studied in China, and several policies and technological innovations have been developed to address this issue in industries, including coking, chemical, steel, and solvent industries [28,29]. However, these studies mostly focus on industrial parks, so that there is limited information regarding the VOC emissions in industrial cities in China.
Changzhi is a typical industrial city in China, with recent rapid economic and population growth. From 2014 to 2018, the total industrial output value of Changzhi increased from RMB 77.6 billion to RMB 89.2 billion, and the total industrial output accounted for more than 50% of the gross product of Changzhi (http://www.changzhi.gov.cn/gzsj/tjgb/) (accessed on 15 July 2019). Steel, electricity, coal, and coking industries are the main industries in Changzhi. In addition, the number of vehicles in the city increased from 386,000 to 528,000 from 2014 to 2018. With the increasing industrial development and number of vehicles in Changzhi, ozone concentration reached 218 µg·m −3 in the summer of 2018, exceeding the national air quality limit (160 µg·m −3 ), and it has become a serious problem in Changzhi during the summers. Nevertheless, recent studies on the air pollution in Changzhi have mainly focused on PM and heavy metal pollution, with limited information on VOCs [30,31].
In this study, the concentrations of 57 VOCs were measured at an urban site in Changzhi from 15 April 2018 to 19 October 2018 (summer to autumn). The dominant emission sources were identified by principal component analysis-multiple linear regres-sion (PCA-MLR). The ozone formation potential (OFP) of VOC species was investigated based on the maximum incremental reactivity (MIR) method. We also conducted a health risk assessment for the investigated VOCs in Changzhi. The results provide important insights to clarify the source apportionment of VOCs in industrial cities, thereby providing a background for governments to develop more efficient strategies for the control and mitigation of VOC emissions and ozone pollution in Changzhi.

Site Description and Sample Collection
The sampling site was the rooftop of the Changzhi Environmental Monitoring Station (113.12 • N, 36.19 • E), which was 20 m aboveground, with adjacent buildings of almost the same height. The sampling port was set at a height of 1.5-2 m from the roof. The building is located in the center of Changzhi, which is a region with a representative dense population. Air samples were collected using Summa canisters (3 L, Entech, Simi Valley, CA, USA) in stainless steel with a Silonite surface coating on the inside of the canisters. The gas flow rate was controlled via a flow valve (8 h/L, Entech, Simi Valley, CA, USA) by the passive sampling method, which was directly docked to the upper end of the summa canister to obtain a 24 h average air sample. The samples were collected in the spring (April to May), summer (June to August), and autumn (September to October) of 2018 with the interval of 5 or 6 days, which effectively characterized the VOCs in the air during different seasons. Finally, a total of 32 valid samples were collected during the sampling period. The concentrations of NO 2 and O 3 during the sampling period were obtained from China City Air Quality Real-time Release Platform (http://106.37.208.233:20035/) (accessed on 30 October 2018). Meteorological parameters (temperature, relative humidity, wind speed, wind direction, and atmosphere pressure) were obtained from automatic weather stations at the same site of VOC sampling.

Sampling and Analysis
All VOC samples were analyzed in a laboratory in Beijing using a three-stage preconcentration system (7200 concentrator, Entech, Entech, Simi Valley, CA, USA), followed by gas chromatography (GC) analysis (Agilent, 7890B, Santa Clara, CA, USA). In this analytical setup, dual columns were configured and coupled to a mass spectrometer (MS) (Agilent, 5977B, Santa Clara, CA, USA) for the analysis of C4-C12 compounds and to a flame ionization detector (FID) for C2-C3 compounds, including ethane, ethylene, acetylene, propylene, and propane. The FID was used to quantify C2-C3 compounds on a HP-PLOT/Q column (30 m × 0.53 mm × 40 µm, Agilent, Santa Clara, CA, USA). To trap VOCs, 400 mL of air was drawn from each Summa canister through a liquid nitrogen cryogenic trap at −40 • C. The trapped VOCs were then transferred to a secondary trap at −70 • C using pure helium and Tenax-TA (sorbent). Finally, the VOCs were transferred to a tertiary trap at −190 • C. In the above steps, most CO 2 , O 2 , H 2 O, and N 2 were removed. Subsequently, the sample was desorbed at >60 • C and injected into the GC equipment with high purity helium (99.99%) as the carrier gas. The MS was used to quantify VOCs with a C4-C12 carbon backbone on a DB-1 column (60 m × 0.25 mm × 0.25 µm, Agilent, Santa Clara, CA, USA). A multi-ramp temperature program was used as follows: 35 • C (10 min), 5 • C min −1 to 150 • C (5 min hold), and 10 • C min −1 to 210 • C (5 min hold). The total running time was 49 min. The FID was maintained at 250 • C and supplied with high purity hydrogen and compressed air. The MS ionization source temperature was 250 • C, and the scan rate was 4.2 scans per second. The analytical system was calibrated using prepared standards and comprising a minimum of five concentrations and a blank. The system stability was verified before each use. Internal standards (benzene-d6, 2-bromo-1,1,1-trifluoroethane and chlorobenzene-d5) were added to the sample stream prior to the trap.
Quality control was achieved by calibrating the instruments before acquiring and analyzing samples throughout the study period so as to ensure the authenticity and accuracy of the experimental data. The instrument blank was analyzed every 24 h, and its results should be lower than the laboratory detection limit. The C2-C5 hydrocarbons and other VOCs were quantified using external and internal standard methods, respectively. The standard gas was measured once a day before the sample was analyzed, and the relative standard deviation (RSD) for each substance response factor was less than 30%. If the RSD was greater than 30%, the standard curve was recreated.

Air Back Trajectories
HYSPLIT was used to calculate the source and long-range transport paths of VOCs during the study period at the Changzhi Environmental Monitoring Station (113.12 • N, 36.19 • E) with an altitude setting of 500 m and a trajectory duration of 24 h. Trajectory clustering was performed using TrajStat to obtain the major air mass source paths. The detailed calculation method is described in the related article [30].

Principal Component Analysis-Multiple Linear Regression
Principal component analysis-multiple linear regression (PCA-MLR) using the "dimensionality reduction" method, transforms the original variables into a smaller number of new variables (principal components) that are orthogonal to each other and reflect the information provided by the original data to the greatest extent possible. Feng's study compared the FA-NNC, PCA-MLR, and PMF methods in the source apportionment results and showed that PCA-MLR performed best in calculating contributions [32]. PCA can be used for studies with small sample sizes, and Choi's study showed that more than 25 samples were needed to produce stable PCA results [33].
The purpose of the multiple linear regression is to determine the contribution of different pollution sources based on the source identification. The equation is: In this study, y denotes ΣVOCs, p denotes the number of extracted principal components. X i denotes the factor score variables obtained from the principal component analysis. b denotes the information of the remaining variables not explained by the factors. The independent and dependent variables were standardized, and the standardized equations were regressed for analysis.
B i is the coefficient of the requested multiple linear regression, so the source contribution was calculated as follows:

Analysis of Ozone Formation Potential
The level of VOC reactivity and sources are key factors for policy makers to consider when aiming to control local ozone contamination. There are different methods used to estimate O 3 formation from VOCs. Among them, the simple calculation of ozone formation potentials (OFPs) with maximum incremental reactivity (MIR) values of different VOC species has been widely applied in other studies [34,35]. It is calculated by Equation (4), where OFP(i) and Conc(i) are the molar mass of O 3 and individual VOC species, respectively, and MIR(i) is the maximum incremental reactivity coefficient of this VOC species.

Health Risk Estimation
To identify and estimate the health risk of specific VOCs to human health, we used the Integrated Risk Information System (IRIS) by the USEPA. The inhalation risks of noncarcinogenic and carcinogenic VOCs were assessed using Equation (5): where EC (µg·m −3 ) is the concentration of inhalation exposure, CA is the average pollutant concentration (µg·m −3 ), ET is the exposure time (24 h/d), EF is the exposure frequency (365 d/y), ED is the exposure duration (70 y), and AT is the average exposure time (70 × 365 × 24 h). The carcinogenic (R) and non-carcinogenic risks (HQ) of all risk-posing VOC species were estimated using the following equations: R = EC × IUR, where R is the estimated inhalation cancer risk, and IUR (m 3 /µg) is the inhalation unit risk, and data were obtained from USEPA's IRIS; and HQ = EC/(RfC × 1000), where HQ is the non-cancer risk hazard quotient (HQ) of specific VOCs, and RfC (mg·m −3 ) is the chronic reference concentration (mg·m −3 ). Cancer risks lower than 1 × 10 -6 were regarded as acceptable. In addition, we used HI = ∑HQ for the analysis, where HI is the sum of the HQ of several pollutants. If HQ or HI were below 1, there was no chronic toxicity, but if HQ exceeded 1, chronic toxicity after long-term exposure was expected.

Situation of Sampling Point Data
The mean wind speed, temperature, and relative humidity during the study period were 3.1 m·s −1 , 19.  Figure 1). The temperature in Changzhi was higher in the summer (23.1 • C) than in the spring (17.4 • C) and autumn (14.7 • C). Similarly, the humidity was higher in the summer (70%) than in the spring (58%) and autumn (61%). VOC concentrations in autumn were higher than those in the spring and summer, a feature consistent in cities such as Beijing, Nanjing, Wuhan, and Hohhot [12,38,39]. Compared to the spring, the higher temperatures in the summer increased the VOC emissions from sources such as solvent evaporation. In addition, the high temperatures, low relative humidity, and strong solar radiation in the summer are conducive to atmospheric oxidation reactions, which can accelerate the depletion of VOCs. The low wind speeds (3.0 m s −1 ) and low temperatures in autumn reduced the rate of atmospheric photochemical reactions, which can lead to the accumulation of VOCs in the atmosphere and the decrease of O 3 concentrations. The O 3 concentration was higher in the summer (218 µg·m −3 ) than in the autumn (154 µg·m −3 ). However, the mean NO 2 concentration, which is also a precursor of O 3 , was lower in the summer (16 µg·m −3 ) than in the autumn (30 µg·m −3 ). In addition, the photochemical reaction rate decreased at lower temperatures and light intensities under stronger north-westerly wind conditions in the autumn, which decreased the consumption of NO 2 and VOCs, thereby increasing their concentrations in the air.
Atmosphere 2022, 13, x FOR PEER REVIEW 6 of 16 the photochemical reaction rate decreased at lower temperatures and light intensities under stronger north-westerly wind conditions in the autumn, which decreased the consumption of NO2 and VOCs, thereby increasing their concentrations in the air. Alkanes were the main components (66.8%) of VOCs in Changzhi, followed by aromatic compounds (16.4%), alkenes (10.2%), and alkynes (6.6%). Higher alkane concentrations have also been observed in Beijing (58%) and Tianjin (63%) [40,41]. The proportion of different components of atmospheric VOCs in Changzhi presented clear seasonal changes, with higher levels of alkanes in the spring (70.3%) than in the summer (66.3%) and autumn (63.8%). The high proportion of alkanes in the spring was consistent with results obtained for other Chinese cities [12,42]. The contribution of aromatic compounds was higher in the summer (19.5%), followed by autumn (16.9%) and spring (15.6%). The high concentration of aromatics in the summer was also similar to results obtained for Shanghai and Tianjin [8,43]. The intensity of emissions from solvent evaporation and the chemical industry, which are the main sources of aromatics, increased in the summer because of the higher temperatures.
As shown in Figure 2, alkanes in the top 10 VOCs with the highest concentrations accounted for 53.7% of the total VOCs, and included isopentane (19.0%), ethane (9.5%), propane (5.2%), n-butane (8.1%), 2-methylheptane (4.6%), and isobutane (3.3%). Alkanes were followed by aromatics (12.2%), including benzene (7.9%) and m/p-xylene (4.3%), alkenes (5.2%, ethylene), and alkynes (3.3%, acetylene). Isopentane concentrations (11.8 µg·m −3 ) were high throughout the seasons, with the highest concentrations in the summer (20.8 µg·m −3 ), followed by spring (12.0 µg·m −3 ) and autumn (6.1 µg·m −3 ). As isopentane is a VOC commonly emitted by gasoline volatilization, its high concentrations were attributed to the extensive emissions from gasoline vehicles [4]. Ethane (5.9 µg·m −3 ) was the second most emitted VOC, and its concentration was lower in the summer (5.5 µg·m −3 ) than in the spring (6.6 µg·m −3 ) and autumn (6.4 µg·m −3 ). Ethane is a typical VOC emitted from coal combustion [44], and its relatively low concentration in the summer was attributed to the lower contribution of coal combustion to ambient VOCs in the summer. The concentrations of aromatic compounds such as benzene and m-xylene, which were the more reactive VOCs, were higher in the summer (5.6 and 3.3 µg·m −3 , respectively) than in the spring (4.6 and 1.3 µg·m −3 , respectively) and autumn (4.4 and 2.8 µg·m −3 , respectively). This higher concentration of aromatics in the summer inevitably produced more ozone, which increased ozone concentrations. Propane and n-butane were also within the top 10 VOCs identified in Changzhi, with average concentrations of 5.8 and 5.0 µg·m −3 , respectively. They are typical VOCs emitted by LPG leaks [45], which is a major source of VOCs in Changzhi. Alkanes were the main components (66.8%) of VOCs in Changzhi, followed by aromatic compounds (16.4%), alkenes (10.2%), and alkynes (6.6%). Higher alkane concentrations have also been observed in Beijing (58%) and Tianjin (63%) [40,41]. The proportion of different components of atmospheric VOCs in Changzhi presented clear seasonal changes, with higher levels of alkanes in the spring (70.3%) than in the summer (66.3%) and autumn (63.8%). The high proportion of alkanes in the spring was consistent with results obtained for other Chinese cities [12,42]. The contribution of aromatic compounds was higher in the summer (19.5%), followed by autumn (16.9%) and spring (15.6%). The high concentration of aromatics in the summer was also similar to results obtained for Shanghai and Tianjin [8,43]. The intensity of emissions from solvent evaporation and the chemical industry, which are the main sources of aromatics, increased in the summer because of the higher temperatures.

Ratio of Specific Compounds
The ratios between VOC species have been widely used to identify VOC sources and associated chemical reactions [40,45,46]. When two species with low reactivity or similar OH reaction constants (kOH) are emitted from the same source, their ratios should not change significantly by photochemical degradation. However, if they share a common source but have different chemical reactivity, their ratios are expected to change upon chemical reactions [47].
The VOC ratios of ethylbenzene to m-xylene (E/X) and benzene to m-xylene (B/X) are widely used to analyze ambient air aging because of their different photochemical reactivity [21,48]. The E/X and B/X ratios vary with the photochemical reaction of the air mass. When E/X is higher than 0.33 or B/X is higher than 1.7, the given region is more likely to be controlled by aging air masses [49]. The results in Table 1 show an E/X of 0.6 and a B/X of 2.8 for Changzhi, which indicate that atmospheric VOCs were probably controlled by aging air masses. The high E/X and B/X in the spring indicated a greater degree of air mass aging, which could be related to the influence of air mass transport in the atmosphere. The average wind speed in spring (3.6 m·s −1 ) in Changzhi is higher than that in summer (2.9 m·s −1 ) and autumn (3.0 m·s −1 ), which leads to longer air mass transport distances and therefore greater air mass aging in spring. As shown in Figure S1, the results of our analysis of the backward trajectory at the time of sampling in Changzhi indicated that the air masses transported in spring are mainly from the north upwards compared to summer and autumn. Air mass in spring transport was mainly from the northern and northwestern directions, with relatively low emission intensities at the places where it was mainly experienced. Air mass transport in summer and autumn all pass through places with relatively large VOC emissions in Henan Province and Hebei Province, and primary emissions had a significant effect on transport and low oxidation.

Ratio of Specific Compounds
The ratios between VOC species have been widely used to identify VOC sources and associated chemical reactions [40,45,46]. When two species with low reactivity or similar OH reaction constants (k OH ) are emitted from the same source, their ratios should not change significantly by photochemical degradation. However, if they share a common source but have different chemical reactivity, their ratios are expected to change upon chemical reactions [47].
The VOC ratios of ethylbenzene to m-xylene (E/X) and benzene to m-xylene (B/X) are widely used to analyze ambient air aging because of their different photochemical reactivity [21,48]. The E/X and B/X ratios vary with the photochemical reaction of the air mass. When E/X is higher than 0.33 or B/X is higher than 1.7, the given region is more likely to be controlled by aging air masses [49]. The results in Table 1 show an E/X of 0.6 and a B/X of 2.8 for Changzhi, which indicate that atmospheric VOCs were probably controlled by aging air masses. The high E/X and B/X in the spring indicated a greater degree of air mass aging, which could be related to the influence of air mass transport in the atmosphere. The average wind speed in spring (3.6 m·s −1 ) in Changzhi is higher than that in summer (2.9 m·s −1 ) and autumn (3.0 m·s −1 ), which leads to longer air mass transport distances and therefore greater air mass aging in spring. As shown in Figure S1, the results of our analysis of the backward trajectory at the time of sampling in Changzhi indicated that the air masses transported in spring are mainly from the north upwards compared to summer and autumn. Air mass in spring transport was mainly from the northern and northwestern directions, with relatively low emission intensities at the places where it was mainly experienced. Air mass transport in summer and autumn all pass through places with relatively large VOC emissions in Henan Province and Hebei Province, and primary emissions had a significant effect on transport and low oxidation. The greatest sources of pentane are LPG combustion, gasoline vehicle exhaust, and fuel evaporation [15,21,44,50]. The mean i-pentane/n-pentane ratio in this study was 8.1 ± 7.1, indicating that VOC emissions in Changzhi were influenced by vehicle exhaust or associated sources. The lower i-pentane/n-pentane ratio (3.4) in autumn indicated that the sources of high isopentane emissions were significantly reduced in that season. Previous studies have shown that the i-butane/n-butane ratio varies according to source, with values of~0.2-0.3 for vehicle exhaust, 0.46 for LPG emissions, and 0.6-1.0 for NG emissions [51][52][53]. The correlation between isobutane and n-butane in Changzhi was very strong (R 2 = 0.74), and the i-butane/n-butane ratio was 0.8, which indicates that NG emissions had a significant impact on the butane concentration. In addition, the higher i-butane/n-butane ratio (1.0) in the summer indicated that NG emissions contributed more to the VOC concentration in the summer than in other seasons.

Source Apportionment by PCA-MLR
PCA-MLR is often used to assess the contribution of various sources to the concentration of pollutants in the ambient atmosphere. PCA-MLR applies multiple methods to transform the correlation matrix composed of receptor data to extract the source information. The steps of PCA are to first transform the collected data orthogonally, then perform a maximum variance rotation, then select the principal factors (eigenvalues greater than 1 or cumulative contribution of variance value between 85-95%), and finally multiple linear regression based on the calculated source component spectrum to estimate the contribution of each type of pollutant source to the pollutant concentration. In this study, there were 57 VOCs monitored online, and 26 of them were selected for source analysis, considering data integrity and substance types and concentrations. The selected basic statistics of VOCs are shown in Table S1. In this paper, the principal component analysis of VOCs was performed using SPSS 21.0 software, and the explained total variance table and rotation matrix table were obtained as follows. There were six factors with eigenvalues greater than 1 obtained, accounting for 89.34% of the total variance. Based on the VOC species, six sources were identified, namely solvent evaporation, vehicle exhaust, industrial production, coal combustion, biogenic emissions, and LPG/NG leaking. These sources were classified into six factors, as shown in Table 2. Due to the limited data used in this study, there is uncertainty in the results of the sources of VOCs, and more observations could better explain the sources of VOCs in Changzhi. Factor 1 was distinguished by a high proportion of ethane, ethylene, propene, acetylene, and benzene. Benzene, toluene, ethylbenzene, and xylene (BTEX) compounds are usually emitted by coal combustion, solvent evaporation, and industrial processes [54]. Factor 1 was defined as coal combustion considering that coking, iron and steel, and power industries are typical industries in Changzhi that consume large amounts of coal.
Factor 2 encompassed other industrial production. The industrial enterprises in Changzhi were diverse and at large scale, including electricity, coking, iron and steel, and chemical sectors. High contributions of C3-C6 alkanes and aromatics were observed in factor 2. Iron and steel industries emit large amounts of aromatics [55], and the chemical and coking industries produce large amounts of alkanes and aromatics [8,56]. Therefore, factor 2 was identified as industrial production.
Factor 3 consisted mainly of alkanes (such as propane, n-butane, n-pentane and 2,3-dimethylbutane) and aromatics (such as benzene). C3-C6 alkanes and alkenes are typical products of the incomplete combustion of gasoline or LPG [37,57]. In addition, the high proportion of aromatic species is a main characteristic of gasoline. Therefore, factor 3 was defined as vehicle exhaust [58].
Factor 4 comprised biogenic emissions. In 2018, the forest area of Changzhi was 445,800 ha, and the forest coverage rate was 31.9%. Biogenic emissions can release large amounts of isoprene, especially during the summer [59]. Isoprene is the predominant species in factor 4, and therefore factor 4 is considered to be the biogenic emission.
In factor 5, 1-butene and isobutane are the predominant species, and they are usually considered to be from NG and LPG leaking [60]. In recent years, the natural gas sector in Changzhi has rapidly developed, aiming to replace coal for residential fuel. However, natural gas can leak during pipeline transportation and use [41,61]. Therefore, this category was defined as LPG/NG leaking. Factor 6 is distinguished by ethylbenzene, m/p-xylene, isopropylbenzene, n-propylbenzene, and 3-methylhexane. Previous studies have shown that isopropylbenzene and n-propylbenzene can be used as a solvent or as an intermediate in organic synthesis, as well as in textile dyes and printing [21,62]. Therefore, we consider this factor as solvent evaporation. Figure 3 shows the source apportionment of VOCs according to the PCA. Vehicle exhaust was the largest source (27.5%) of atmospheric VOCs, and coal combustion was the second largest one (23.5%). Additionally, the contributions of industrial production, solvent evaporation, LPG/NG, and biogenic emissions to atmospheric VOCs were 17.4%, 13.5%, 9.5%, and 8.7%, respectively.
Changzhi has rapidly developed, aiming to replace coal for residential fuel. However, natural gas can leak during pipeline transportation and use [41,61]. Therefore, this category was defined as LPG/NG leaking. Factor 6 is distinguished by ethylbenzene, m/p-xylene, isopropylbenzene, npropylbenzene, and 3-methylhexane. Previous studies have shown that isopropylbenzene and n-propylbenzene can be used as a solvent or as an intermediate in organic synthesis, as well as in textile dyes and printing [21,62]. Therefore, we consider this factor as solvent evaporation. Figure 3 shows the source apportionment of VOCs according to the PCA. Vehicle exhaust was the largest source (27.5%) of atmospheric VOCs, and coal combustion was the second largest one (23.5%). Additionally, the contributions of industrial production, solvent evaporation, LPG/NG, and biogenic emissions to atmospheric VOCs were 17.4%, 13.5%, 9.5%, and 8.7%, respectively. The contribution of different sources in each season was calculated from the coefficients of each source for different samples and VOC concentration data in the PCA calculation results. The contribution of different VOC sources varied between seasons. The contribution of biogenic emission was significantly higher in the summer (13.2%) than in the spring (6.8%) and autumn (6.2%). The temperatures are higher in summer (23.0 °C) than in spring (17.4 °C) and autumn (15.0 °C), which leads to accelerated plant growth and higher respiration rates, resulting in more emissions from biogenic sources in summer, whereas the contributions of vehicle exhaust were generally consistent across the seasons, at 29.1%, 26.2%, and 27.2% in the spring, summer, and autumn, respectively. The lowest contribution of coal combustion to VOCs was in the summer (18.6%), followed by autumn (24.3%) and spring (27.5%). The highest contribution of industrial production was in the summer (23.3%), which might be related to the increased intensity of industrial emissions in that season. The contribution of LPG/NG was significantly lower in the summer (6.7%) than in the spring (11.1%) and autumn (10.7%). Finally, the contribution of solvent evaporation was mostly consistent across the seasons, at 14.6%, 12.1%, and 13.8% in the spring, summer, and autumn, respectively. Figure S1 showed that the summer air mass transport in Changzhi was mainly from northern Henan Province (51.8%). Henan Province is one of the major grain producing areas in China. Summer crop growth generated a large number of biogenic emissions that reach Changzhi through air mass transport, resulting in a high proportion of biogenic contribution. In addition, both summer and autumn air masses pass through the industrial agglomerations in central Henan province. The transport of industrial emissions of VOCs through the air masses may be one of the reasons for the high proportion of industrial contributions in summer and autumn in Changzhi. The contribution of different sources in each season was calculated from the coefficients of each source for different samples and VOC concentration data in the PCA calculation results. The contribution of different VOC sources varied between seasons. The contribution of biogenic emission was significantly higher in the summer (13.2%) than in the spring (6.8%) and autumn (6.2%). The temperatures are higher in summer (23.0 • C) than in spring (17.4 • C) and autumn (15.0 • C), which leads to accelerated plant growth and higher respiration rates, resulting in more emissions from biogenic sources in summer, whereas the contributions of vehicle exhaust were generally consistent across the seasons, at 29.1%, 26.2%, and 27.2% in the spring, summer, and autumn, respectively. The lowest contribution of coal combustion to VOCs was in the summer (18.6%), followed by autumn (24.3%) and spring (27.5%). The highest contribution of industrial production was in the summer (23.3%), which might be related to the increased intensity of industrial emissions in that season. The contribution of LPG/NG was significantly lower in the summer (6.7%) than in the spring (11.1%) and autumn (10.7%). Finally, the contribution of solvent evaporation was mostly consistent across the seasons, at 14.6%, 12.1%, and 13.8% in the spring, summer, and autumn, respectively. Figure S1 showed that the summer air mass transport in Changzhi was mainly from northern Henan Province (51.8%). Henan Province is one of the major grain producing areas in China. Summer crop growth generated a large number of biogenic emissions that reach Changzhi through air mass transport, resulting in a high proportion of biogenic contribution. In addition, both summer and autumn air masses pass through the industrial agglomerations in central Henan province. The transport of industrial emissions of VOCs through the air masses may be one of the reasons for the high proportion of industrial contributions in summer and autumn in Changzhi. Figure 4 shows the source apportionment of VOCs in seven Chinese cities during a similar period, including Shanghai, Nanjing, Handan, Wuhan, Tianjin, Beijing, Zhengzhou, and Weinan [22,23,40,41,60,[63][64][65]. Vehicle exhaust was the main source of VOCs in several cities. The contribution of vehicle exhaust (27%) in Changzhi was lower than those in densely populated cities such as Beijing (33%) and Shanghai (28%). As a typical industrial city, the contribution of industrial production to VOCs in Changzhi was higher than that of southern cities such as Wuhan (12%) and Handan (11%) but lower than that of the northern Chinese city of Tianjin (20%). In addition, the contribution of solvent evaporation in Changzhi was significantly lower than those in the southern cities of Wuhan (37%) and Nanjing (23%). In recent years, the coal chemical and pharmaceutical industries in Changzhi have grown rapidly. The production processes in these industries consume large amounts of organic materials and produce large amounts of VOCs. Therefore, industrial production and solvent evaporation have become important sources of local VOCs. Nevertheless, steel, electricity and coking are still the most important industries in Changzhi, and they have a high energy demand, which is mostly produced by coal combustion. Therefore, the contribution of coal combustion was significantly higher in Changzhi than in other cities. The contribution of biogenic emissions to VOCs was similar across all investigated cities, ranging from 6 to 9%. Figure 4 shows the source apportionment of VOCs in seven Chinese cities during a similar period, including Shanghai, Nanjing, Handan, Wuhan, Tianjin, Beijing, Zhengzhou, and Weinan [22,23,40,41,60,[63][64][65]. Vehicle exhaust was the main source of VOCs in several cities. The contribution of vehicle exhaust (27%) in Changzhi was lower than those in densely populated cities such as Beijing (33%) and Shanghai (28%). As a typical industrial city, the contribution of industrial production to VOCs in Changzhi was higher than that of southern cities such as Wuhan (12%) and Handan (11%) but lower than that of the northern Chinese city of Tianjin (20%). In addition, the contribution of solvent evaporation in Changzhi was significantly lower than those in the southern cities of Wuhan (37%) and Nanjing (23%). In recent years, the coal chemical and pharmaceutical industries in Changzhi have grown rapidly. The production processes in these industries consume large amounts of organic materials and produce large amounts of VOCs. Therefore, industrial production and solvent evaporation have become important sources of local VOCs. Nevertheless, steel, electricity and coking are still the most important industries in Changzhi, and they have a high energy demand, which is mostly produced by coal combustion. Therefore, the contribution of coal combustion was significantly higher in Changzhi than in other cities. The contribution of biogenic emissions to VOCs was similar across all investigated cities, ranging from 6 to 9%.

Ozone Formation Potential
Because of the large differences in the photochemical reactivities of VOC species, the contributions of different VOC species to the photochemical generation of ozone are different. To better understand the characteristics of VOC species in the atmosphere and determine their role in the photochemical generation of ozone, we used MIR methods to analyze and describe the reactivity and OFP of VOC species.
The contribution of major VOC groups to the OFP is shown in Figure 5. This study assessed the effect of VOCs on ozone formation using the OFP, which was based on the results of the MIR study. The calculated OFP of each VOC species sampled during the observation period is shown in Figure 5. The VOCs contributed to 125.8 µg·m −3 of the atmospheric ozone generation in Changzhi during the study period. Alkenes were the largest contributor to ozone formation in Changzhi (40.9% of OFP), followed by aromatics (28.6%), alkanes (28.2%), and alkynes (2.3%).

Ozone Formation Potential
Because of the large differences in the photochemical reactivities of VOC species, the contributions of different VOC species to the photochemical generation of ozone are different. To better understand the characteristics of VOC species in the atmosphere and determine their role in the photochemical generation of ozone, we used MIR methods to analyze and describe the reactivity and OFP of VOC species.
The contribution of major VOC groups to the OFP is shown in Figure 5. This study assessed the effect of VOCs on ozone formation using the OFP, which was based on the results of the MIR study. The calculated OFP of each VOC species sampled during the observation period is shown in Figure 5. The VOCs contributed to 125.8 µg·m −3 of the atmospheric ozone generation in Changzhi during the study period. Alkenes were the largest contributor to ozone formation in Changzhi (40.9% of OFP), followed by aromatics (28.6%), alkanes (28.2%), and alkynes (2.3%).
Although the concentrations of alkenes and aromatics were significantly lower than those of alkanes, their chemical reactivities were relatively high, and their contributions to the concentration of ozone were greatly increased. The contribution of VOCs to ozone production was higher in June and October. The increased emissions of alkanes, especially isopentane, led to the higher OFP in June, whereas the weaker atmospheric photochemical reactions in autumn can explain the increased OFP in October. Although the concentrations of alkenes and aromatics were significantly lower than those of alkanes, their chemical reactivities were relatively high, and their contributions to the concentration of ozone were greatly increased. The contribution of VOCs to ozone production was higher in June and October. The increased emissions of alkanes, especially isopentane, led to the higher OFP in June, whereas the weaker atmospheric photochemical reactions in autumn can explain the increased OFP in October. As isopentane and O3 concentrations are indicators of vehicle exhaust emissions, their variation suggested that vehicle exhaust in Changzhi had a significant influence on atmospheric O3.
The contribution of the emission sources to ozone production and to VOC pollution is different because the emitted pollutants are subjected to different chemical activities. Coal combustion was the largest source of ozone production, contributing to 28.9% of the OFP. Industrial production was associated with high VOC contribution and contributed to 18.0% of OFP, followed by LPG/NG (15.3%), solvent evaporation (16.7%), vehicle exhaust (12.8%), and biogenic emissions (8.3%). Ethylene and propylene are characteristic species from coal combustion emissions, and their high concentrations and high chemical reactivity led to the greatest contribution to OFP. In addition, solvent evaporation was also a main source to OFP owing to the chemical reactivity and high concentrations of aromatics emitted from solvent evaporation.
The contribution of different sources to OFP also changed according to season. The increased emissions from the local vegetation in the summer resulted in a higher contribution of biogenic emission to OFP. Coal combustion contributed less to the OFP in the summer because of its diminished contribution to VOCs in that season. These changes were attributed to the different VOC concentrations and reactivity of the indicator species. Therefore, source apportionment based on OFP can provide a basis for the development of countermeasures to reduce O3 concentrations.

Human Health Assessment
As a standard to estimate overall risk, six aromatic compounds (benzene, toluene, ethylbenzene, p-xylene, o-xylene, and styrene) were selected for the health risk assessment. The non-carcinogenic risks (HQ) of all compounds were lower than the criterion, ranging from 1.5 × 10 −1 to 1.5 × 10 −5 . The HI (average of 2.0 × 10 −1 ) was higher in the summer (2.4 × 10 −1 ) than in the spring (1.7 × 10 −1 ) and autumn (1.9 × 10 −1 ), and they were all below the criterion during the study period, which indicates that exposure to these benzenes would not pose a non-carcinogenic risk to the public. The R of benzene (average of 3.8 × 10 −5 ) was higher than the criterion (1 × 10 −6 ) in all seasons, at 3.4 × 10 −5 , 4.4 × 10 −5 , and 3.4 × 10 −5 for the spring, summer, and autumn, respectively. The average R of ethylbenzene in the study period was 2.9 × 10 −6 , and it was higher than the criterion in all seasons, at 2.2 × 10 −6 , 3.4 × 10 −6 , and 2.8 × 10 −6 in the spring, summer, and autumn, respectively. The high R indicated that long-term exposure to benzene and ethylbenzene in the study area can The contribution of the emission sources to ozone production and to VOC pollution is different because the emitted pollutants are subjected to different chemical activities. Coal combustion was the largest source of ozone production, contributing to 28.9% of the OFP. Industrial production was associated with high VOC contribution and contributed to 18.0% of OFP, followed by LPG/NG (15.3%), solvent evaporation (16.7%), vehicle exhaust (12.8%), and biogenic emissions (8.3%). Ethylene and propylene are characteristic species from coal combustion emissions, and their high concentrations and high chemical reactivity led to the greatest contribution to OFP. In addition, solvent evaporation was also a main source to OFP owing to the chemical reactivity and high concentrations of aromatics emitted from solvent evaporation.
The contribution of different sources to OFP also changed according to season. The increased emissions from the local vegetation in the summer resulted in a higher contribution of biogenic emission to OFP. Coal combustion contributed less to the OFP in the summer because of its diminished contribution to VOCs in that season. These changes were attributed to the different VOC concentrations and reactivity of the indicator species. Therefore, source apportionment based on OFP can provide a basis for the development of countermeasures to reduce O 3 concentrations.

Human Health Assessment
As a standard to estimate overall risk, six aromatic compounds (benzene, toluene, ethylbenzene, p-xylene, o-xylene, and styrene) were selected for the health risk assessment. The non-carcinogenic risks (HQ) of all compounds were lower than the criterion, ranging from 1.5 × 10 −1 to 1.5 × 10 −5 . The HI (average of 2.0 × 10 −1 ) was higher in the summer (2.4 × 10 −1 ) than in the spring (1.7 × 10 −1 ) and autumn (1.9 × 10 −1 ), and they were all below the criterion during the study period, which indicates that exposure to these benzenes would not pose a non-carcinogenic risk to the public. The R of benzene (average of 3.8 × 10 −5 ) was higher than the criterion (1 × 10 −6 ) in all seasons, at 3.4 × 10 −5 , 4.4 × 10 −5 , and 3.4 × 10 −5 for the spring, summer, and autumn, respectively. The average R of ethylbenzene in the study period was 2.9 × 10 −6 , and it was higher than the criterion in all seasons, at 2.2 × 10 −6 , 3.4 × 10 −6 , and 2.8 × 10 −6 in the spring, summer, and autumn, respectively. The high R indicated that long-term exposure to benzene and ethylbenzene in the study area can cause potential carcinogenic risks. Therefore, their emissions should be controlled to ensure public health.
The relative contribution of each source to health risk was evaluated. The sources containing more toxic pollutants were expected to generate higher health risks. In this study, solvent evaporation contributed to 28.0% of the non-carcinogenic risk, and 19.1% of the carcinogenic risk. The contribution of solvent evaporation to the carcinogenic risk was attributed to its high load of aromatics. Coal combustion contributed to 29.0% and 25.3% of the non-carcinogenic and carcinogenic risks because of its high VOC emissions. Therefore, the health risks of different sources should be considered when developing abatement measures for VOC pollution.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The datasets analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest:
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.