Volatile Organic Compounds in a Petrochemical Region in Arid of NW China: Chemical Reactivity and Source Apportionment

We measured volatile organic compounds (VOCs) during the heating, non-heating, and sandstorm periods in the air of the Dushanzi district in NW China and investigated their concentrations, chemical reactivity, and sources. The observed concentrations of total VOCs (TVOCs) were 22.35 ± 17.60, 33.20 ± 34.15, and 17.05 ± 13.61 ppbv in non-heating, heating, and sandstorm periods, respectively. C2-C5 alkanes, C2-C3 alkenes, benzene, and toluene were the most abundant species, contributing more than 60% of the TVOCs. Among these VOCs, alkenes such as propene had the highest chemical reactivity, accounting for more than 60% of total hydroxyl radical loss rate (LOH) and ozone formation potential (OFP). Chemical reactivity was the highest in the heating period. The average reaction rate constant (KOH-avg) and average maximum incremental reactivity coefficient (MIR-avg) of the total observed VOCs were (8.72 ± 1.42) × 10−12 cm3/mol·s and 2.42 ± 0.16 mol/mol, respectively. The results of the source apportionment via the Positive Matrix Factorization (PMF) model showed that coal combustion (43.08%) and industrial processes (38.86%) were the major sources of VOCs in the air of the Dushanzi district. The contribution of coal combustion to VOCs was the highest in the heating period, while that of industrial solvents and oil volatilization was the lowest.


Introduction
Volatile organic compounds (VOCs) are important substances in the atmosphere because they can produce ozone and other oxidants via photochemistry and consequently lead to the increase of atmospheric oxidative capacity and the deterioration of air quality [1][2][3][4].Toxic VOCs (e.g., benzene, toluene, ethylbenzene, and 1,3-butadiene) can stimulate the human respiratory system, damage the nervous system, and increase the risk of cancer [5][6][7][8].
VOCs from petrochemical industries account for 14.5% of the total VOCs emissions in China [9].Generally, VOCs in petrochemical industrial areas mainly come from refining, storage tanks, other processes of leakage, and volatilization.Previous studies have reported the characteristics and chemical reactivity of VOCs in the air of petrochemical industrial parks where aromatics, alkanes, alkenes, and halohydrocarbons were detected [10][11][12][13].Han et al. (2018) [14] found that reducing the emissions of aromatic hydrocarbons and isoprene could effectively reduce the formation of haze in a petrochemical industrial region in the Yangtze River Delta, China.Moghadam et al. (2013) [15] reported that the concentrations of VOCs in the Mahshahr petrochemical area were higher in summer than in winter, with a significant impact of wind speed on the VOCs' distribution.Tiwari et al. (2010) [16] studied the VOCs in the petrochemical area of Yokohama, Japan, and found that aromatics were the largest group among which toluene was the highest in concentrations and m/p-xylene was the largest in the ozone formation potential (OFP) values.Mao et al. (2018) [17] measured the VOCs in a petrochemical industry park in Northwest China and found that alkanes were the largest group and some alkanes also exhibited the highest OFP values.
The Dushanzi district, covering an area of 448 km 2 , is located in Xinjiang province, China, at north latitude 44 • 07 -44 • 23 , east longitude 84 • 43 -85 • 06 .Dushanzi district has a temperate continental climate with less rain and stronger northwest winds during sandstorm (0.96 m/s) and non-heating (0.73 m/s) periods, but with weaker wind during the heating period (0.39 m/s).The temperature was the highest in the non-heating period and the lowest in the heating period.By contrast, relative humidity was the highest during the heating period and lowest during the sandstorm and non-heating periods [18].The large day-night temperature difference and long sunshine duration provide favorable meteorological conditions for the transformation of VOCs [19].Dushanzi is one of the most important petrochemical industrial bases in Northwest China.VOC leakage from industrial production has had an adverse impact on the air quality of Dushanzi.According to data provided by the Dushanzi environmental monitoring station (838600), 41% of the days had maximum 8-h average of O 3 (MDA8) near the surface exceeding the national standard (0.10 mg/m 3 ) in 2013.Our previous study illustrated the chemical characteristics and oxidative damage of polycyclic aromatic hydrocarbons (PAHs, semi-VOCs) in the ambient air of the Dushanzi residential area [20].However, at present, no data for VOCs are available in this petrochemical industrial region.Hence, it is essential to investigate the concentrations, chemical reactivity and sources of VOCs in the Dushanzi district.
In this study, we investigated VOCs' characteristics during the non-heating, heating, and sandstorm periods in the Dushanzi district.We further estimated the chemical reactivity of VOCs by calculating the OFP values and OH• loss rate (L OH ) and attributed the sources of VOCs via the Positive Matrix Factorization (PMF) model.

Sample Collection
The sampling site was on the rooftop of a residential building in Dushanzi district, Xinjiang (north latitude 44 • 19'43", east longitude 84 • 53'0").The building is about 25 m above the ground and is taller than its surroundings.The petrochemical industrial area is to the northwest of the sampling site, and Quitun city is to the north (Figure 1).The sampling site is surrounded by residential buildings, commercial streets, and major roads.Therefore, the sampling site is representative of a mixed-use functional area in the Dushanzi district integrating residential, traffic, and industrial impacts.The campaign was conducted during the non-heating (14 to 19 September 2015, n = 21), heating (10 to 16 December 2015, n = 21), and sandstorm (19 to 24 April 2016, n = 19) periods, respectively.One-hour VOC sampling (based on US EPA TO-14) was conducted using 2-L stainless-steel canisters three times a day (9:30-10:30, 14:30-15:30, and 22:30-23:30).The interior surfaces of the canisters were passivated.Each canister was cleaned five times using a high-purity nitrogen gas and evacuated.After that, one of ten canisters was selected and filled with nitrogen for 24 h, and then analyzed to ensure that the target compounds were not detected.During sampling, a filter was added at the front end of the pipeline to prevent particles from entering into canisters.To avoid contamination, the pressure in each canister was about 30 psi with an oil-free pump at the completion of each sampling.After the sampling, the samples were stored at room temperature and sent to the lab for analysis after a month.completion of each sampling.After the sampling, the samples were stored at room temperature and sent to the lab for analysis after a month.

Chemical Analysis
The VOC samples were analyzed with a Model 7100 preconcentrator (Entech Instrument Inc, USA) coupled with a gas chromatography-mass selective detector/flame ionization detector/electron capture detector (GC-MSD/FID/ECD, 6890/Agilent 5973N and 6890/Agilent 5975).In total, 350 mL of air in each canister was extracted through a three-stage cryogenic trapping system.In the first stage, the VOCs were absorbed on the glass beads at −160 °C and desorbed at 10 °C to enrich the VOC species.In the second stage, the VOCs were trapped at −40 °C and desorbed at 180 °C to remove CO2 and water.Finally, the VOCs were concentrated in a transfer line at −170 °C and heated rapidly to 180 °C for injection into the GC-MSD/FID/ECD system by helium.The VOCs were first separated by a capillary column (HP-1, 60 m × 0.32 mm × 1.0 µm, Agilent Technologies, USA) with helium and then split into two paths: one leading to a PLOT-Q column (30 m × 0.32 mm × 2.0 µm, Agilent Technologies, USA) followed by a FID for the detection of C2-C3 hydrocarbons; and the other was to an MSD for the detection of C4-C11 hydrocarbons.The GC temperature was programmed to −50 °C initially, held for 3 min, increased to 10 °C at 15 °C min −1 and then to 120 °C at 5 °C min −1 .Finally, it was increased to 250 °C at 10 °C min −1 and held for 10 min.The MSD was operated in the selected ion monitoring (SIM) mode.Detailed information about the analytical processes is given elsewhere [21].Quantification was performed by an external calibration method [22].The target compounds were identified based on their retention time and mass spectrum.

Quality Assurance and Quality Control
Every day, high-purity nitrogen gas was used as a blank sample to ensure that the target VOCs were not detected.A mixed standard (100 mL) with a concentration of −1 ppbv was also tested every day to check the changes in the instrument responses.If the difference in the responses between the two days exceeded 20%, another calibration curve was generated.

Chemical Analysis
The VOC samples were analyzed with a Model 7100 preconcentrator (Entech Instrument Inc, USA) coupled with a gas chromatography-mass selective detector/flame ionization detector/electron capture detector (GC-MSD/FID/ECD, 6890/Agilent 5973N and 6890/Agilent 5975).In total, 350 mL of air in each canister was extracted through a three-stage cryogenic trapping system.In the first stage, the VOCs were absorbed on the glass beads at −160 • C and desorbed at 10 • C to enrich the VOC species.In the second stage, the VOCs were trapped at −40 • C and desorbed at 180 • C to remove CO 2 and water.Finally, the VOCs were concentrated in a transfer line at −170 • C and heated rapidly to 180 • C for injection into the GC-MSD/FID/ECD system by helium.The VOCs were first separated by a capillary column (HP-1, 60 m × 0.32 mm × 1.0 µm, Agilent Technologies, USA) with helium and then split into two paths: one leading to a PLOT-Q column (30 m × 0.32 mm × 2.0 µm, Agilent Technologies, USA) followed by a FID for the detection of C 2 -C 3 hydrocarbons; and the other was to an MSD for the detection of C 4 -C 11 hydrocarbons.The GC temperature was programmed to −50 • C initially, held for 3 min, increased to 10 • C at 15 • C min −1 and then to 120 • C at 5 • C min −1 .Finally, it was increased to 250 • C at 10 • C min −1 and held for 10 min.The MSD was operated in the selected ion monitoring (SIM) mode.Detailed information about the analytical processes is given elsewhere [21].Quantification was performed by an external calibration method [22].The target compounds were identified based on their retention time and mass spectrum.

Quality Assurance and Quality Control
Every day, high-purity nitrogen gas was used as a blank sample to ensure that the target VOCs were not detected.A mixed standard (100 mL) with a concentration of −1 ppbv was also tested every day to check the changes in the instrument responses.If the difference in the responses between the two days exceeded 20%, another calibration curve was generated.

Estimation of the Initial Concentration of VOCs
VOCs are a group of reactive substances whose life-time in the air depends highly on chemical reactivity.The concentrations of the hydroxyl radical ([OH]) range from a few hours (e.g., isoprene and α-pinene) to several days (e.g., benzene, propane, and toluene) [23].A previous study by Gouw et al. (2005) [24] demonstrated that the initial concentrations of VOCs in the air could be estimated based on the losses of VOCs in the air with the hydroxyl radical.Since ethylbenzene (E) and m/p-xylene (X) have the same emission sources [25], we selected E and X to calculate the initial concentrations of each VOC as described in Equations ( 1) and ( 2) [26]: where [VOC i ] t and [VOC i ] 0 are the observed and initial concentrations of each VOC, respectively.K i is the reaction rate constant of VOC i with the hydroxyl radical [23].
[OH] is the concentration of the hydroxyl radical.Importantly, results do not depend on the assumed [OH] mixing ratio.∆t is the reaction time, and K E and K X are the reaction rate constants of ethylbenzene and m/p-xylene with the hydroxyl radical, respectively.According to previous studies, the [E]/[X] ratios in the petrochemical industry (i.e., solvent use and coal combustion ([E]/[X] t = t0 )) were reported to be in the ranges of 0.5-0.7,0.4-0.5, and 0.3-0.4,respectively [27,28].We selected an [E]/[X] t = t0 value of 0.5 to calculate ∆t and the initial concentration of each VOC with the assumption that the VOCs were well mixed in the air.

Analysis of Ozone Formation Potential
In this study, we estimated the ozone formation potential of each measured VOC and thereby identified the key species in O 3 formation.The OFP was calculated as the concentration of each VOC and the maximum incremental reactivity (MIR) coefficient according to Equation (3) [23,29]: where OFP i is the ozone formation potential of species i, MIR i stands for maximum incremental reactivity coefficient of VOC i .

OH Loss Rate
VOC oxidation with hydroxyl radicals is one of the most important reactions in the air.The chemical reactivity of VOC species can be evaluated by calculating the hydroxyl radical loss rate (L OH ) according to Equation (4) [4,19]: where L i OH is the hydroxyl radical loss rate of [VOC i ], [VOC i ] is the atmospheric concentration of compound i (molecules/cm 3 ), and K i OH is the reaction rate constant of species i (cm 3 /molecule•s), which was reported by [23].

Source Analysis with Positive Matrix Factorization
The Positive Matrix Factorization (PMF) model was used to investigate the sources of VOCs.Detailed information about PMF has been described in the literature [30,31].PMF 5.0 was used to attribute the sources of VOCs in the air of the Dushanzi district.If the concentrations of VOC species were higher than the method detection limits (MDLs), the uncertainties were calculated via the following equation: Otherwise, the calculation of the uncertainties was u ij = MDL × 5/6.Due to the strong reactivity and short life-time of biogenic terpenoids (e.g., isoprene and monoterpenes), we only apportioned anthropogenic VOCs in this study [32,33].

Characteristics of VOCs
In this study, a total of 57 VOCs were measured including 25 alkanes, 17 alkenes, 14 aromatics, and 1 alkyne.The observed and initial concentrations of the VOCs are shown in Figure 2. The observed concentrations of the total VOCs (TVOCs) were 22.35 ± 17.60, 33.20 ± 34.15, and 17.05 ± 13.61 ppbv in non-heating, heating, and sandstorm periods, respectively.The initial concentrations of the TVOCs were 31.94 ± 12.20, 34.51 ± 22.60, and 21.48 ± 10.66 ppbv in the non-heating, heating, and sandstorm periods, respectively, which are higher than the observed concentrations by ~18%.Both observed and initial levels of TVOCs were the highest in the heating period.This situation was mainly due to meteorological factors, especially wind speed.The average wind speeds were 2.5, 1.4, and 3.7 m/s, respectively, during the non-heating, heating, and sandstorm periods.High wind speed during the non-heating and sandstorm periods lowered the concentrations of pollutants, while low wind speed during the heating period led to pollutant accumulation [21].The concentrations of aromatics were the highest in the non-heating period and the lowest in the heating period, likely due to high humidity and low temperatures in the heating period, facilitating the chemical transformation of aromatics from a gas phase to a particle phase [34], and the low temperature was not conducive to volatile aromatics.Low concentrations of aromatics in the sandstorm period were caused by the low temperatures and high wind speeds, which were beneficial to the diffusion of aromatics.
Alkanes were the major species, with a mean concentration of 14.01 ± 14.33 ppbv, followed by alkenes (4.31 ± 4.38 ppbv) and aromatics (2.04 ± 3.32 ppbv).The VOC compositions were similar to those previously reported in petrochemical industrial areas [35,36].C 2 -C 5 alkanes, C 2 -C 3 alkenes, benzene, and toluene were the most abundant species, accounting for more than 60% of the TVOCs.Acetylene (2.19 ± 2.70 ppbv) contributed 10% of the TVOCs.For the initial VOCs, the fractions of alkenes and aromatics in the initial TVOCs were higher than those in the observed TVOCs, largely due to the high chemical reactivity of alkenes and aromatics.
We compared our measurements with those in other cities of China (Table 1).Dushanzi had lower levels of TVOCs than Nanjing, Guangzhou, Taiyuan, and Jinan, but higher levels than Beijing and Shuozhou.For the major components, alkanes in Dushanzi were lower than those in Nanjing and Guangzhou, and higher than those in Beijing, Taiyuan, Shuozhou, and Jinan.The levels of aromatics in Dushanzi were the lowest among these cities in China.

Chemical Reactivity of VOCs
L OH and OFP were the key parameters describing the chemical reactivity of VOCs.The chemical reactivity of VOCs might be underestimated in this study because the K OH and MIR of some species were not available.Figure 3 shows the L OH and OFP values for each VOC group during the three sampling periods.The chemical reactivity of the initial VOCs was ~20% higher than that of the observed VOCs.Alkenes had the highest chemical reactivity, with contributions of 65% to the total chemical reactivity, followed by alkanes (17%), aromatics (15%), and alkynes (1%).These results are similar to those reported for other cities in China [41][42][43][44], and indicate the priority of alkenes in VOC control.It can be seen that the total L OH and OFP of both the observed and initial VOCs during the heating period were much higher than those in the non-heating and sandstorm periods.
Further examination suggested that the C 2 -C 5 alkenes, toluene, and C 8 aromatics were the largest contributors to the total L OH , accounting for 70% and 50% of the observed and initial VOCs' L OH , respectively.For the OFP, C 2 -C 5 alkenes contributed ~60% to the observed VOC OFP, followed by toluene, m/p-xylene, styrene, and o-xylene (7%).Although the L OH and OFP values of the alkanes were relatively low, their contributions to chemical reactivity could not be ignored due to their high concentrations.
The concentrations of O 3 were 67.71, 30.04, and 66.79 µg/m 3 during the non-heating, heating, and sandstorm periods, respectively, and the consumed VOCs were 9.59, 1.31, and 4.48 ppbv.Thus, the more VOCs that were consumed, the more O 3 was generated.C 2 -C 5 alkenes and C 7 -C 9 aromatics were the most important contributors (~60%) to the consumed VOCs.The contributions of α-pinene and β-pinene were also high during the non-heating period.
Atmosphere 2019, 10, x FOR PEER REVIEW 9 of 16 The concentrations of O3 were 67.71, 30.04, and 66.79 µg/m 3 during the non-heating, heating, and sandstorm periods, respectively, and the consumed VOCs were 9.59, 1.31, and 4.48 ppbv.Thus, the more VOCs that were consumed, the more O3 was generated.C2-C5 alkenes and C7-C9 aromatics were the most important contributors (~60%) to the consumed VOCs.The contributions of α-pinene and β-pinene were also high during the non-heating period.The top ten species with the highest L OH and OFP values in the atmosphere of Dushanzi (in Table 2) accounted for 62.33%, 70.03%, and 68.68% for total L OH and 64.07%, 78.65%, 74.87% for total OFP during the non-heating, heating, and sandstorm periods, respectively.Propene was the dominant contributor to HO• loss and ozone formation, accounting for 12.51%, 26.59%, 18.08% and 18.09%, 31.93%, and 24.00% of the total L OH and OFP during the three sampling periods, respectively.Other large contributors included isoprene, n-hexane, ethylene, and trans-2-butene.

Average Reactivity
Figure 4 shows the significant correlations between the levels of VOCs and the values of L OH and OFP (R 2 > 0.7, p < 0.05 for all).Based on the positive linear correlations, we estimated the average reaction rate constant (k OH -avg) of the total VOCs.The k OH -avg in the air of Dushanzi was (8.72 ± 1.42) × 10 −12 cm 3 /molecules•s for the observed VOCs and (11.62 ± 3.25) × 10 −12 cm 3 /molecules•s for the initial VOCs, which was close to the k OH vale of ethylene (8.52 × 10 −12 cm 3 /molecules•s) but much lower than that of propylene (26.3 × 10 −12 cm 3 /molecules•s) [23].This suggests that the chemical reactivity of VOCs in the air of Dushanzi was relatively low.Additionally, the average maximum incremental reactivity coefficient (MIR-avg) of the observed and initial VOCs was 2.42 ± 0.16 mol/mol and 2.96 ± 0.29 mol/mol, respectively.This also indicates the relatively low chemical reactivity of the VOCs in the air of Dushanzi.
The top ten species with the highest L OH and OFP values in the atmosphere of Dushanzi (in Table 2) accounted for 62.33%, 70.03%, and 68.68% for total L OH and 64.07%, 78.65%, 74.87% for total OFP during the non-heating, heating, and sandstorm periods, respectively.Propene was the dominant contributor to HO• loss and ozone formation, accounting for 12.51%, 26.59%, 18.08% and 18.09%, 31.93%, and 24.00% of the total L OH and OFP during the three sampling periods, respectively.Other large contributors included isoprene, n-hexane, ethylene, and trans-2-butene.

Average Reactivity
Figure 4 shows the significant correlations between the levels of VOCs and the values of L OH and OFP (R 2 > 0.7, p < 0.05 for all).Based on the positive linear correlations, we estimated the average reaction rate constant (k OH -avg) of the total VOCs.The k OH -avg in the air of Dushanzi was (8.72 ± 1.42) × 10 -12 cm 3 /molecules•s for the observed VOCs and (11.62 ± 3.25) × 10 -12 cm 3 /molecules•s for the initial VOCs, which was close to the k OH vale of ethylene (8.52 × 10 -12 cm 3 /molecules•s) but much lower than that of propylene (26.3 × 10 -12 cm 3 /molecules•s) [23].This suggests that the chemical reactivity of VOCs in the air of Dushanzi was relatively low.Additionally, the average maximum incremental reactivity coefficient (MIR-avg) of the observed and initial VOCs was 2.42 ± 0.16 mol/mol and 2.96 ± 0.29 mol/mol, respectively.This also indicates the relatively low chemical reactivity of the VOCs in the air of Dushanzi.

Source Apportionment by PMF
The PMF model was used to attribute the VOC sources in the air of Dushanzi.Five sources were identified in Dushanzi, including oil volatilization, vehicle exhaust, industrial processes, coal combustion, and industrial solvents (Figure 5).
Source 2 was characterized by a high percentage of some heavier alkanes (e.g., n-heptane, 2,3,4-trimethylpentane, 2-methylheptane, and 3-methylheptane) which were abundant in diesel exhaust, especially in vehicle-used diesel [46].N-octane, n-decane, and n-nonane were typical markers of diesel vehicle exhaust [47,48].Aromatics such as toluene, ethylbenzene, o-xylene, and m/p-xylene in urban areas have been related to motor vehicles [21].Thus, source 2 was related to vehicle exhaust.Source 3 was dominated by C 3 -C 6 alkanes (e.g., propane, n-butane, 2,3-dimethylbutane, n-pentane, 2-methylpentane, and n-hexane), benzene, and C 2 -C 4 alkenes (ethylene, propene, and 1,3-butadiene).Industrial process losses can release alkanes and alkenes.For example, n-hexane and ethylene are widely used in the manufacturing industries.Benzene is the basic product in the chemical industries and is also a basic raw material in petrochemical industries [47].The annual production of i-pentane, n-pentane, and cyclopentane in the Dushanzi district were estimated to be 70,000 tons.Thus, source 3 was related to industrial processes.
Source 4 was characterized by high fractions of ethane, propane, i-butane, acetylene, ethylene, 1-butene, trans-2-butene, and benzene.Acetylene is a typical product of incomplete combustion [49].Some studies have noted that 1-butene is an indicator for coal-fired power plants [50][51][52].C 2 species such as ethylene, ethane, and acetylene are abundant in the exhaust of coal combustion [53].As a typical tracer for coal combustion, sulfur dioxide has the same trend as ethane and acetylene (Figure S1).Thus, source 4 was assigned to coal combustion.
Source 5 was dominated by aromatics (ethylbenzene, m-ethyltoluene, p-ethyltoluene, m/p-xylene, styrene, 1,3,5-trimethylbenzene, 1,2,3-trimethylbenzene, and others).Toluene, ethylbenzene, xylene, and 1,3,5-trimethylbenzene are used as organic solvents in various industries and can be regarded as indicator species for industrial solvent volatilization [54].Therefore, source 5 was identified as industrial solvents.Figure 6 shows the source apportionment of the VOCs via the PMF model.Coal combustion was the most important contributor (43.08%), and industrial processes (including industrial solvent volatilization and industrial process leakages) accounted for 38.86%.Coal combustion is the main source of energy in the Dushanzi district, so it accounts for the largest proportion.Since Dushanzi has a small population and urban size, the contributions of oil volatilization and vehicle exhaust are Figure 6 shows the source apportionment of the VOCs via the PMF model.Coal combustion was the most important contributor (43.08%), and industrial processes (including industrial solvent volatilization and industrial process leakages) accounted for 38.86%.Coal combustion is the main source of energy in the Dushanzi district, so it accounts for the largest proportion.Since Dushanzi has a small population and urban size, the contributions of oil volatilization and vehicle exhaust are relatively small.We then analyzed the VOC sources during the non-heating, heating, and sandstorm periods, respectively (Figure S2).We found that the contribution of coal combustion to VOCs during the heating period was much higher than that of during the other sampling periods (Figure 7).This was mainly due to the increase in coal combustion during the winter.At the same time, industrial solvents and oil contributed less to VOCs because of the lower temperatures during the heating period.The contribution of the sources of VOCs during the non-heating period were similar those in the sandstorm period, because the decrease of coal combustion during the non-heating and sandstorm periods and the higher temperatures and wind speeds lead to a significant increase in the volatilization of industrial solvents and oils.Author Contributions: The study was completed with cooperation between all authors.D.T., X.W., X.D., G.W., Y.Z. have conceived and designed the experiments; X.W. provided the laboratory apparatus and accommodation; D.T., X.D., A.A. have provided valuable comments and suggestions for the modification of the manuscript.G.W., X.Z. have performed graphic processing and data analysis; G.W., X.Z.wrote the paper and X.D. contributed grammar correcting of all text.

Figure 1 .
Figure 1.Location of the sampling site and residential area in the Dushanzi district.

Figure 1 .
Figure 1.Location of the sampling site and residential area in the Dushanzi district.

Figure 2 .
Figure 2. The (a) observed concentrations and (b) initial concentrations of 57 volatile organic compounds (VOCs) during the non-heating, heating, and sandstorm periods.The inset illustrations represent the observed and initial concentrations of alkenes, alkanes, aromatics, and alkynes during the three sampling periods.

Figure 2 .
Figure 2. The (a) observed concentrations and (b) initial concentrations of 57 volatile organic compounds (VOCs) during the non-heating, heating, and sandstorm periods.The illustrations represent the observed and initial concentrations of alkenes, alkanes, aromatics, and alkynes during the three sampling periods.

Figure 3 .
Figure 3.The observed and (a-b) initial (c-d) average OH• loss rate (L OH ) and ozone formation potential (OFP) of each group of VOCs of for the non-heating, heating, and sandstorm periods.

Figure 3 .
Figure 3.The observed and (a,b) initial (c,d) average OH• loss rate (L OH ) and ozone formation potential (OFP) of each group of VOCs of for the non-heating, heating, and sandstorm periods.

Figure 4 .
Figure 4.The relationship between the total L OH /OFP and TVOCs of the (a-b) observed and (c-d) initial VOCs for the sampling periods.

Figure 4 .
Figure 4.The relationship between the total L OH /OFP and TVOCs of the (a,b) observed and (c,d) initial VOCs for the sampling periods.

Figure 5 .
Figure 5. Source profiles resolved from Positive Matrix Factorization (PMF) during sampling periods in Dushanzi.

Figure 5 .
Figure 5. Source profiles resolved from Positive Matrix Factorization (PMF) during sampling periods in Dushanzi.

Figure 7 .
Figure 7. Contributions of the five sources to VOCs in Dushanzi during three sampling periods.Supplementary Materials: Figure S1: Relationship between ethane or acetylene and SO2 for sampling period; Figure S2: apportionment of VOCs in the non-heating (a), heating (b) and sandstorm (c) sampling periods.

Table 1 .
Observed VOC concentrations compared among major cities in China.

Table 1 .
Observed VOC concentrations compared among major cities in China.

Table 2 .
Top ten species with the highest L OH and OFP values in the atmosphere of the Dushanzi district.

Table 2 .
Top ten species with the highest L OH and OFP values in the atmosphere of the Dushanzi district.