Impacts of Meteorological Factors, VOCs Emissions and Inter-Regional Transport on Summer Ozone Pollution in Yuncheng

: Summer ozone (O 3 ) pollution in China has become increasingly serious in recent years. This study is based on hourly data of near-surface ozone (O 3 ) and nitrogen oxides (NOx) and volatile organic compounds (VOCs) from June to August 2020 in Yuncheng, combined with meteorological data to analyse the characteristics of O 3 pollution in summer and the inﬂuence of meteorological factors, precursors, and long-range transport on O 3 pollution. In this paper, the VOCs/NOx characteristic ratio method was used to explore the sensitivity of O 3 generation. Backward trajectories, cluster analysis, potential source contribution factor (PSCF) analysis and concentration weight trajectory (CWT) analysis were also calculated using Trajstat software. In 2020, Yuncheng had persistent O 3 pollution, with the highest concentrations in June, signiﬁcantly higher than July and August. Conditions of high temperature, low relative humidity and low wind speed contribute to the O 3 accumulation. VOCs are the main precursors to the local production of O 3 . Besides, the long-range transport analysis shows that southeast-oriented air masses are the main direction inﬂuencing summer O 3 pollution. The primary potential source areas of O 3 are in the central and southern part of Henan province, the north-western Anhui province, and the northern Shaanxi. In addition, northern Hubei and southwestern Shandong also inﬂuence O 3 pollution in summer Yuncheng.


Introduction
After implementing the Chinese government's air pollution prevention and control plan, China has achieved remarkable achievement for the frequent PM 2.5 pollution in heating seasons [1,2]. However, tropospheric ozone (O 3 ) pollution in summer has increased year by year [3,4]. According to statistics, O 3 concentrations in key cities in mainland China increased year by year after 2013, with the most pronounced increase in the Beijing-Tianjin-Hebei Region (BTH), Yangtze River Delta (YRD) and Pearl River Delta (PRD) regions [5]. O 3 concentrations across Eastern China increased by about 10% in recent years [6]. Therefore, near-surface O 3 concentrations are a critical factor influencing summer air quality in China.
Near-surface O 3 is a typical secondary pollutant, and its primary source is the photochemical reaction of atmospheric volatile organic compounds (VOCs) and nitrogen oxides (NOx) under solar radiation conditions [7]. Although the stratosphere may transport some O 3 to the troposphere, the concentration is only 1/15 to 1/7 of the O 3 produced by photochemical reactions [8]. At the same time, O 3 pollution tends to occur more readily during Hourly O 3 concentration monitoring data from June to August 2020 are the primary data used in this research. This paper includes four parts: (1) analysed O 3 and its precursors' pollution characteristics, (2) showed the influences of meteorological factors on O 3 concentration, (3) analysed the relationship between precursors and O 3 production by using the characteristic ratio method (VOCs/NOx), (4) using potential source contribution factor (PSCF) analysis and concentration weight trajectory (CWT) analysis provided by Trajstat software to study the regional transport characteristics of O 3 and its potential source areas in different seasons in Yuncheng City. It is expected to provide a scientific reference basis for O 3 pollution control and emission reduction policy formulation in Yuncheng City and other urban clusters in the Fenwei Plain.  35.02 • N) is in the southwestern part of Shanxi province, surrounded by the Luliang mountains, Zhongtiao mountains and Yellow River, respectively. The topography of Yuncheng is dominated by low mountains, hills, plains and Chinese Loess Plateau. The climate is warm-temperate continental monsoon with cold and dry winters, hot and rainy summers, with an average annual temperature and rainfall of 13 • C and 537 mm, respectively. Nearly half of the rain concentrated from July to September [23].

Data Sources
O 3 and NOx concentration data and meteorological data (temperature, wind speed, wind direction, relative humidity) for Yuncheng from 1 January 2018 to 31 December 2020 can be obtained from the National Air Quality Networked Monitoring Management Platform at hourly resolution.
The monitoring station for this study was the roof of Yuncheng Middle School, Yanhu District, Yuncheng City, Shanxi Province (35.038 • N, 111.029 • E). The monitoring site is surrounded by mainly residential and commercial areas and is close to the main traffic roads, a typical urban area. The operational heights of the monitoring points for O 3 and NOx are selected according to the technical regulation for selection of ambient air quality monitoring stations (HJ 664-2013) for automatic monitoring equipment. The specification states that if the average height of buildings within a radius of 300-500 m around the selected point is above 20 m, the height of the sampling port or monitoring beam above the ground should be set within the range of 15-25 m [24]. Since the building near this sampling point is seven storeys (>20 m) in height, the selected point in this study is at the top of Yuncheng Middle School, with a height of 23 m from the ground. The sampling ports of the monitoring instruments are all greater than 1 m away from the surface of the building walls, roofs, and other supports.
The analysers used for the O 3 and NOx observations were manufactured by Thermo Fisher, a US-based thermoelectric environmental instrument company. O 3 concentrations were measured using a Model 49i O 3 analyser with a UV absorption method. This instrument first uses a pump to draw ambient air into the absorption cell. The reaction cell is then irradiated with UV light at a wavelength of 254 nm. The detector will monitor the unabsorbed light passing through the unit, and the Beer-Lambert equation and output will calculate the O 3 concentration through the screen. NOx concentrations are measured using the chemiluminescence-based Model 42i NOx analyser. Ambient air enters the system via a pump and dryer unit and is split into two streams. The NO in the first stream reacts with O 3 to form active NO 2 , which quickly changes from an excited state to a ground state and emits a photon. The released photon results from the chemiluminescent radiation, which depends on the amount of NO in the air. Finally, radiated photons are detected using a photomultiplier or photodiode, which outputs voltages of different intensities. The NO 2 in the second gas stream can be reduced to NO, and then the sum of the NO and NO 2 concentrations are measured according to the above principle to obtain the NOx concentration. According to G. Villena et al., the concentration of NOx measured by chemiluminescence in urban areas is higher than other methods [25]. However, this study focuses on the effect of NOx trends on O 3 formation and pays no attention to the NOx's concentration changes. Therefore, this error has little influence on the results of the study.
This study uses a Wuhan Tianhong TH-300B atmospheric VOC online monitoring system to observe VOCs in the environment, mainly consisting of an ultra-low temperature online pre-concentrator and a GC-MS/FID detector. After the atmospheric samples have been collected and entered the ultra-low temperature online pre-concentrator, the VOCs in the environment are frozen and trapped in a capillary trap column. When rapid thermal desorption by the system, they enter the GC-MS analysis system. The components from C 2 to C 5 are detected in the FID detector, while the MS detects the remaining ingredients. In order to ensure the observation accuracy, aldehydes and ketones, PAMS, TO-15 and internal standards were used in the monitoring process. The three standard gases, aldehydes and ketones, PAMS and TO-15, are certified as BW with an uncertainty of 5%. In contrast, the typical internal gas is licensed as GBWE and contains four substances, namely chlorobromomethane, chlorobenzene-d5, 4-bromofluorobenzene and 1,4-difluorobenzene. The standard gases are diluted and calibrated using zero gas (de-hydrocarbonized air) and the responding standard gas without humidification during dilution. It is worth noting that the standard gas used in the VOCs online monitoring process is non-humid. However, the air compressor draws the carrier gas directly from the atmosphere, so the mixed standard and carrier gases will be humid. This measurement process will essentially have no significant effect on the VOCs monitoring values. Six concentration points (0.5, 2, 4, 6, 8 and 10 µg·L −1 ) were set up for multi-point calibration. At least five concentration points were fitted to verify that the correlation coefficient of the target compounds at the time of calibration was 0.98 or higher. Besides, single-point calibration was performed daily using standard gases so that more than 80% of the target compounds had a concentration deviation of less than 30% and no more than five key compounds failed. A weekly blank check with zero gas ensured that the monitoring system was not contaminated and that the blank met the target compound concentration of less than 0.1 µg·L −1 . During the observation period, 68 VOCs were detected, including 29 alkanes, 10 olefins, 1 alkyne, 16 aromatic hydrocarbons and 12 oxygenated volatile organic compounds (OVOCs).

Hysplit Model
The Hysplit model, developed jointly by the National Oceanic Atmospheric Administration (NOAA) and the Australian Bureau of Meteorology, can be used to calculate and analyse airflow movement, deposition and dispersion [26,27]. It is also often used by researchers in air pollution to analyse long-range emission sources and transport analysis of different types of pollutants.
The backward trajectory of Yuncheng City, cluster analysis, PSCF and CWT analysis were all calculated based on Trajstat software [27]. The researcher set the starting point and the height of the simulation as the centre of Yuncheng city (110.98 • E, 35.02 • N) and 500 m, respectively, because this height can reflect the inter-regional transportation flow characteristics and reduces the influence of near-surface friction [28]. The duration of the simulation for this study was 48 h. Trajectories began at four times (00:00, 06:00, 12:00, 18:00, UTC). The clustering is based on the angle distance algorithm [29] and the mutation points total spatial variance (TSV) [26].

Potential Source Contribution Factor Analysis (PSCF)
The PSCF is based on the HYSPLIT backward trajectory endpoint of the receptor model, which can identify potential sources of atmospheric pollutants using spatial grid probability calculations [30,31]. This study divided the research area into 0.5 • × 0.5 • grids, and the pollution threshold of O 3 is set as 144 µg·m −3 (the 75th percentile of the daily average concentration from June to August in 2020) [32]. Pollution thresholds can help define pollution trajectories. When a trajectory corresponding concentration is higher than this value, it is a contaminated transport trajectory.
where, "ij" represents the specific location of each grid, "m ij " shows the number of contaminated trajectory endpoints passing through grid "ij", "n ij " means the endpoints amounts of all trajectories passing through grid "ij". According to the calculation process of PSCF, when the total number of endpoints (n ij ) of the trajectory in the grid is relatively small, the calculation result may have considerable uncertainty. Therefore, this paper uses a weighting factor (W ij ) to enhance its statistical balance and reduce the error [33].
where, "avg" represents the average number of endpoints of all grids in the study area, "W ij " is the weighting factor, "WPSCF ij " is the calculation result of the PSCF by "W ij " equilibrium.

Concentration Weighted Trajectory Analysis (CWT)
Since the PSCF calculation is based on the ratio of the number of contaminated trajectories passing through a grid to the number of endpoints in a whole way, it is possible that when the concentration is a medium to high strength source, the cells that pass through will have the same calculation result [34,35]. CWT is a gridded result of pollution sources obtained by back-projecting using local pollution concentration data combined with air mass trajectories. CWT can be weighted by trajectory residence time and pollutant concentration, taking into account the concentration of pollutants crossing that grid cell and weighting the concentration by residence time [36]. Therefore, CWT can calculate relatively accurate weights for each grid cell traversed at medium and high concentrations, improving the accuracy of the analysis results. Because CWT and PSCF can compensate for shortcomings, researchers always combine the two techniques to identify the source area [34,37].
where, "ij" represents the specific location of each grid, "l " means a particular transportation trajectory, "τ ijl " is the residence time of trajectory "l " in grid "ij" (or the number of trajectory endpoints), "C l " is the concentration of O 3 as trajectory "l" passes through grid "ij".To reduce errors, CWT also applies the weighting factor "W ij " quoted by PSCF:

Ozone Pollution Occurred in Summer 2020
As shown in Figure 2, the time series of O 3 hourly pollution levels in Yuncheng from 2018 to 2020 are divided into five classes according to the O 3 monitoring concentration values, reflecting the degree of exceedance of the hourly O 3 standard (>160 µg·m −3 ) in each month [38]: excellent (0-100 µg·m −3 ), good (100-160 µg·m −3 ), mild pollution (160-215 µg·m −3 ), moderate pollution (215-265 µg·m −3 ) and heavily pollution (265-500 µg·m −3 ) [11]. Yuncheng has experienced varying O 3 exceedances from March to September over three consecutive years, predominating mild and moderate pollution. The higher levels of pollution are usually in the months of June to August each year, during which there are several consecutive hours of monthly average concentrations reaching moderate pollution. Although the O 3 concentration exceedances proportion in 2020 is smaller than in 2019, the pollution characteristics of O 3 are more pronounced in the summer of 2020. Not only did hourly concentrations of heavy pollution occur during the night of June, but moderate pollution also occurred more frequently and continuously. Therefore, 2020 was chosen as the base year to identify and analyse local O 3 generation factors and potential source areas.

Opposite Changes in Ozone and Its Precursors
Both O 3 and its precursors show an apparent seasonal or monthly variation. As can be seen from Figure 3, O 3 concentrations in 2020 show a clear seasonal pattern of summer > autumn > spring > winter, which is very similar to its pollution situation in many countries and regions in the northern hemisphere [39,40]. Because O 3 is a photochemical product, intense solar radiation in the northern hemisphere summer is favourable for O 3 production [41]. Therefore, ozone concentrations tend to reach higher values more easily from June to August. At the same time, O 3 in Yuncheng showed two distinct pollution peaks in 2020, in June and September, with higher pollution levels in June than in September. The decrease in O 3 concentrations in July and August may be due to the increase in air humidity, leading to a weakening of solar radiation intensity and thus photochemical reactions. The relationship between O 3 concentration and RH will be discussed in Section 3.2.
The NOx and VOCs concentrations variation are the opposite of O 3 , showing a high winter and low summer pattern. During the Chinese heating period in winter, continuous heating requires consuming large amounts of solid fuels, which emit more NOx, toluene, and benzene. At the same time, the increased hours of indoor activities of people and the high frequency of private car use also emit more NOx and VOCs [42,43]. Ethane, propane, n-butane, iso-butane, n-pentane, iso-pentane and a proportion of olefins and benzene can be emitted from transportation. From winter to summer, as the weather warms and solar radiation increases, some of the NOx and VOCs are consumed in the photochemical reactions, resulting in a gradual decrease in their concentrations. Therefore, the variation in NOx and VOCs concentrations is closely related to climatic conditions and emission patterns. In addition to the effect of VOCs emissions from anthropogenic sources on O 3 concentrations, emissions from plant sources in summer may also influence O 3 production. The main VOCs that can represent natural sources are isoprene [44]. The contribution of isoprene to VOC concentrations in June, July and August was calculated to be 2.38%, 2.32% and 1.11%, respectively. Because the monitoring station represents the Yuncheng downtown's pollutants concentration, the influence of biogenic VOCs is much smaller than the emissions from the anthropogenic sources. Therefore, this study mainly focuses on the impact of VOCs generated during anthropogenic activities on O 3 in summer in Yuncheng City. Figure 4 uses monthly average hourly data to show the daily trend between O 3 concentrations and temperature (Temp, • C), relative humidity (RH, %) and wind speed (WS, m/s) for June to August 2020. O 3 concentrations show a positive correlation with temperature throughout June, which means that O 3 concentrations increase with climbing air temperature. Solar radiation is the main factor that reflects the intensity of O 3 photochemical reactions. In hotter weather, stronger solar radiation tends to occur, so the air temperature level can also show whether the photochemical reactions that generate O 3 arise at a given time of year [45]. Therefore, high temperatures in summer promote the occurrence of photochemical reactions of O 3 . O 3 contamination showed a continuous peak of approximately 150 µg·m −3 from 14:00 to 17:00. In addition, the RH can also determine the strength of the photochemical reaction of O 3 [46]. Studies have shown that the intensity of the photochemical reaction reaches a critical value when the RH is 60%. When RH is above 60%, there is a decrease in O 3 concentration because cloudy weather tends to occur when RH is high, declining the solar radiation intensity and thus leading to lower O 3 concentrations [45]. The RH can reach around 59% in the early morning of June in Yuncheng when O 3 concentration was less than 100 µg·m −3 , and at the moment of peak O 3 concentrations in the afternoon, the RH was less than 53%.

This phenomenon indicates that the effect of RH on O 3 concentrations is significant.
Besides, the short-range transport degree of O 3 is influenced by the magnitude of WS, which can also lead to O 3 concentrations fluctuations. If under static WS conditions, poorer dispersion conditions for O 3 can lead to a concentration accumulating [47]. At 4:00 in the early morning, there was a significant peak in WS in June, which was about 0.3 m/s higher than other times, so that O 3 concentrations decreased rapidly from 105 µg·m −3 to 82 µg·m −3 . O 3 would diffuse downwind with the near-surface wind, thus reducing the local O 3 pollution level to some extent. However, it is worth noting that, although the increase in wind speed helped to improve diffusion conditions, the eventual decrease in O 3 concentrations was due to the depletion of O 3 by NO titration at night.
The temperature, RH and WS change trends in July were very different in Yuncheng. Firstly, the peak temperature occurred between 8:00 and 9:00, after which it gradually dropped to around 29 • C and fluctuated up and down. The reason that temperature cannot reach higher values in the afternoon may be due to the influence of RH. The peak RH (70%) appeared around 10:00, with a high temperature. Later, the RH dropped rapidly to around 58% and remains stable. This is probably because rainy and humid weather tends to occur in the morning when the overall summer temperatures are high, and evaporation absorbs heat, resulting in a rapid drop in RH and a decrease in temperature. Therefore, the O 3 concentrations were low in the early morning (<100 µg·m −3 ) due to high humidity and low temperatures and gradually increased to 140 µg·m −3 as the RH decreased rapidly. The overall O 3 concentrations are lower in July than in June. The main reason is that, during the time between 14:00 and 15:00, when solar radiation is most robust, temperatures remain relatively low, and dispersion conditions (2.2 m/s) are favourable for downwind dispersion.
Similar to the meteorological conditions in July, the RH in August continues to show a "single valley" pattern. The RH peak occurred at 8:00 and can reach around 78%. The troughs appeared after 16:00, but were still above 67%. The temperature in August peaks at midday and then drops as the RH drops and eventually remains stable (28 • C). Compared to July, the RH and WS were higher in August, while the temperature did not reach the same level as in July, resulting in a decrease in peak O 3 concentrations to around 120 µg·m −3 .
The meteorological factors' influences on O 3 concentration analysis are not a simple linear relationship. There are common effects between different meteorological conditions. Thus, this study used multiple linear regression to analyse the association between meteorological factors and O 3 concentrations. According to Table 1, the temperature in June has a very significant effect on O 3 concentration, while the process of increasing RH in July and August reduces the O 3 concentration considerably. The results are similar to those previously discused. Under high temperature, low humidity conditions will positively affect the photochemical reaction and accumulation of O 3 . The effect of WS on O 3 in Yuncheng shows a positive correlation. When other factors remain constant and WS increases, O 3 concentrations increase to some extent, which indicates that long-range transport is one of the key factors affecting O 3 concentrations. When continuous hot, dry, and soft WS weather occurs during the summer months, more attention should be paid to controlling O 3 precursor emissions and long-distance transportation.

VOCs Are Critical Precursors for Ozone Formation
The NOx concentrations from June to August 2020 were averaged into four intervals according to their percentile to ensure that the number of hourly NOx concentrations within each concentration range was approximately the same [48]. These are the 0-25% (4-11 µg·m −3 ), 25-50% (11-14 µg·m −3 ), 50-75% (14-20 µg·m −3 ) and 75-100% (20-60 µg·m −3 ) quartiles, respectively. Figure 5 reflects the effects of different wind direction NOx concentrations on O 3 in different percentile intervals. When NOx concentrations are below the 50% percentile, the main sources are to the north and northeast of the monitoring site. Whereas, at higher levels of NOx concentrations, the contribution from the northeast and east of the city of Yuncheng gradually dominates. The main reason for this may be that the northeastern part of the monitoring site is the main road of Yuncheng (Hongqi East Street), and traffic sources mainly influence NOx emissions. In addition, Figure 5 also reveals the relationship between NOx and O 3 concentrations. When NOx concentrations are low, most O 3 concentrations are at high levels, with the most significant points lying in the 100-150 µg·m −3 interval. When comparing the O 3 concentrations corresponding to NOx above the 75th percentile, they are less than 100 µg·m −3 . Thus, NOx concentrations show a negative correlation with O 3 concentrations. When the NOx concentration gradually decreased, the O 3 concentration increased rather than decreased, suggesting that O 3 pollution is likely in the VOCs-limited regime [49]. O 3 production is mainly closely related to NOx and VOCs. In Figure 6, the O 3 concentration was divided into three ranges of <100 µg·m −3 , 100-160 µg·m −3 and >160 µg·m −3 to plot a three-dimensional scatter plot of NOx-VOCs-O 3 . According to the characteristic ratio method, 4 and 8 are the criteria for classifying VOCs-limited and NOx-limited regimes. When O 3 concentrations were between 100 and 160 µg·m −3 , the number of points in the NOx-limited regime was significantly lower. Comparing the range of NOx concentrations on clean days shows that NOx concentrations range was more extensive on clean days (0-40 µg·m −3 ), which further confirms the negative correlation between NOx concentration and O 3 concentration. In other words, when NOx concentrations are high, O 3 pollution is reduced.
According to the results on polluted days (>160 µg·m −3 ), almost all of the points fall in the VOCs-limited regime and the transition zone. The dispersion degree of these points decreases. When controlling the NOx concentration as constant, the O 3 concentration decreases with the reduction in VOCs. In contrast, maintaining constant VOCs concentration, as NOx concentration decreases, O 3 concentration increases significantly. So, the main precursors of O 3 production on polluted days are VOCs.
Overall, in terms of the current O 3 pollution situation in Yuncheng, VOCs emission control is more critical to reduce summer O 3 concentration. However, it is worth noting that, as air quality is optimised and O 3 pollution will gradually decline, the precursor control situation will move progressively towards the NOx-limited regime. In the current O 3 control process, the control of VOC concentrations will have a meaningful impact. However, in the future, it will become more and more significant when discussing O 3 pollution control.   According to Figure 6 and Table 1, track 3 s probability is the highest (31.79%), starting from the south of Hebei, then through the north of Henan to Yuncheng. However, the average ρ (O 3 ) and contamination trajectory share are not the highest, around 104.93 µg·m −3 and 23.66%, respectively. The highest probability of pollution track (track 4) from north-western Anhui, arriving at Yuncheng via most of Henan, can be 44.09%. The mean and contaminated trajectory concentrations were 134.89 µg·m −3 and 173.20 µg·m −3 , which are the highest of all tracks. Track 2 from southern Inner Mongolia, across northern Shaanxi to Yuncheng, has the second probability of occurrence after track 3. However, its pollution concentration is low, probably due to its long transmission distance and fast-moving speed, which is not conducive to the accumulation of pollutants. Track 1 from western Hubei has the lowest probability for all trajectories and contaminated paths, at 13.04% and 6.45%, and the lowest pollution concentrations, at 97.43 µg·m −3 and 168.50 µg·m −3 .

Long-Range Transport Sources of Ozone Pollution
Based on the Section 3.4.1 analysis, this research used Trajstat software to do a potential source analysis (WPSCF) of O 3 in Yuncheng (Figure 8a). Areas with a WPSCF value greater than 0.5 were the most significant potential source areas subject to point contamination. In Figure 8a, the primary potential sources of O 3 pollution in summer are in the southeast. The WPSCF results greater than 0.5 are mainly in central and southern Henan, northwestern Anhui and northern Hubei. The regions calculated 0.4 to 0.5 are mainly Shaanxi north, central and west Henan, and Hubei north. As most of the areas with high WPSCF values are in the southeast of Yuncheng, this may be related to the prevailing southeasterly winds aloft during the summer months in Yuncheng. To further explore the degree of pollution contribution of source areas to O 3 in Yuncheng City and to verify the results of the WPSCF analysis, the WCWT method was carried out (Figure 8b). Larger WCWT values represent a more significant pollution contribution to the receiving point, and this study defined 100 µg·m −3 as major contributing areas. The results of the WCWT analysis are similar to the WPSCF. Regions with WCWT values greater than 120 µg·m −3 are mainly in central Henan and northern Hubei, and northwestern Anhui's border. Besides, in north Shaanxi, southwestern Shanxi and southwestern Shandong, which also have sporadic areas greater than 120 µg·m −3 . The WCWT values of 100-120 µg·m −3 were mainly in northern Shaanxi, southwestern Shanxi, and northern Hebei. Some high WCWT also exist in north Henan, southern Inner Mongolia, southern Hebei, southwestern Shandong, north of Hubei and northwestern Anhui.
The main potential source areas for summer O 3 in Yuncheng are in the southeast, mainly in the central and southern part of Henan province and the northwestern part of Anhui province. In addition, northern Shaanxi, northern Hubei and southwestern Shandong also influence summer O 3 pollution.

Conclusions
Compared to 2018 and 2019, the summer of 2020 showed more persistent O 3 pollution in Yuncheng. In particular, the O 3 concentration in June shows a continuity of multi-hour pollution and severity of contamination, even higher than 265 µg·m −3 . Due to photochem-ical reaction conditions, the seasonal characteristics of O 3 pollution in Yuncheng in 2020 are summer > autumn > spring > winter. In summer, the peak pollution concentrations are highest in June (155 µg·m −3 ), significantly higher than in July (145 µg·m −3 ) and August (120 µg·m −3 ). The NOx and VOC concentrations vary opposite to O 3 because their concentrations are influenced by the emission conditions and photochemical reactions degree. The O 3 pollution scenario for 2020 cannot occur without meteorological factors, local generation, and transmission from air masses.
According to the influence of meteorological factors on O 3 concentrations, there is a positive correlation between temperature and O 3 concentrations. A decreasing trend in O 3 concentration occurs gradually when the RH is more significant than 60%. Afternoon peak temperatures in June 2020 were above 30 • C, significantly higher than in July (29 • C) and August (28 • C). At the same time, the RH was higher than 65% in July and August, which contributes to afternoon temperatures staying at lower levels. The low temperatures and high humidity (cloudy weather) in July and August indicate that the intensity of solar radiation is weaker, so O 3 concentrations are significantly lower than in June. According to the multiple linear regression results, the wind speed was positively correlated with O 3 concentrations in Yuncheng from June to August with correlation coefficients of 22.951, 17.046 and 9.050, respectively. Higher O 3 concentrations may occur when higher wind speeds arise, suggesting that O 3 concentrations from regional transport need to be noted.
The local formation of O 3 in Yuncheng is inseparable from the emission of VOCs. O 3 concentration increased when NOx concentration decreased suggests that the VOCs concentration control may be more effective in controlling local O 3 pollution. When O 3 concentrations are high, the VOCs/NOx ratio is less than 8, indicating that O 3 generation is sensitive to VOCs. When O 3 concentrations are low, VOCs/NOx ratio move towards the transition zone and the NOx-limited regime, meaning that when the O 3 pollution level is reduced, the pollution of O 3 by NOx emissions should not be underestimated. According to the current pollution situation in Yuncheng, VOCs are the main control factor for controlling O 3 pollution. The most significant OFP contributing species were acetaldehyde, propionaldehyde, ethylene, isoprene, toluene, m/p-xylene, propylene, isopentane, o-xylene and n-butane. Therefore, the emission control of the top ten OFP VOCs mentioned above will be crucial to managing O 3 pollution.
The results of the transport trajectories analysis show that southeast-oriented air masses are the main track affecting O 3 pollution in Yuncheng in summer 2020, with a seasonal share of 44.09% and a pollution concentration as high as 173.20 µg·m −3 . The WPSCF and WCWT analyses indicate that the main potential source areas of O 3 are mainly in central and southern Henan and northwestern Anhui. In addition, northern Shaanxi, northern Hubei, and southwestern Shandong also influence O 3 pollution in the summer Yuncheng.