Assessment of the Characteristics and Influencing Factors of Ozone in Fuzhou, China, Using Wavelet Analysis

In recent years, ozone (O3) air pollution has become a major problem—which is still increasing—in China. In this study, data on the O3 concentration as well as on two of its influencing factors, meteorological parameters and traffic conditions, were collected in Fuzhou and then analyzed via wavelet transforms to evaluate the variation in O3 concentration and its relationships with the factors. High concentrations appeared during spring and summer, and notable fluctuations occurred primarily from March till April, from June till July and during September. The level of O3 exhibited significant positive correlations with the temperature and sunshine duration but negative ones with the relative humidity, precipitation and air pressure; additionally, it increased when the wind speed was low and fluctuated when the wind speed was high. The concentration was also significantly negatively correlated with the concentrations of O3 precursors (NOx and CO) and hence closely related to the traffic conditions, as the reduced vehicle speed during rush hour due to increased vehicle flow and traffic density led to higher precursor emissions in the vehicular exhaust. A prominent “weekend effect” was observed with the precursor levels, which displayed greater fluctuations during the weekdays than the weekends; the vehicle flow and the O3 concentration.


INTRODUCTION
With the acceleration of industrialization and urbanization, the air pollution problems have drawn unprecedented attention of the society in China. In recent years, the surface ozone (O 3 ) has replaced inhalable particles (PM 10 ) and fine particulate matter (PM 2.5 ) as the primary pollutant in many Chinese cities.
Tropospheric O 3 is a secondary air pollutant, which is produced by the photochemical reaction involving NO x (NO and NO 2 ) and volatile organic compounds (VOCs) in the presence of sunlight. At a height of 10-50 km, stratospheric O 3 forms a gas shield around the earth which can absorb solar ultraviolet radiation. However, the tropospheric O 3 can stimulate and damage mucosal tissues, such as the eyes and respiratory system, and hinder blood oxygen transmission, which will do harm to human health (Jung et al., 2016).
To date, O 3 pollution has been found in most of Chinese areas, such as Beijing-Tianjin-Hebei region (BTH) (Li et al., 2015;Li et al., 2020), the Pearl River Delta (PRD) Zheng et al., 2010;Lin et al., 2019) and the Yangtze River Delta (YRD) (Wu et al., 2017;Liu et al., 2018). To reduce O 3 pollution, many measures have been taken by the Chinese government. During the 12 th Five-Year Plan period (2011)(2012)(2013)(2014)(2015), the nitrogen oxides were included as indicators of air quality, and VOCs were included in the total control range during the 13 th Five-Year Plan period (2016)(2017)(2018)(2019)(2020). Moreover, the VOC Pollution Prevention and Control Plan in the 13 th Five-Year Plan was issued by the Ministry of Environmental Protection of the People's Republic of China (MEP) on September 19, 2017.
A number of O 3 -related pollution studies have been conducted in recent years, including pollution process and source identification of O 3 (Tong et al., 2018;Wang et al., 2019), the effect of ozone on crops and human (Feng et al., 2015;Pun et al., 2019), the spatial and temporal variations of O 3 concentrations Ercan et al., 2019;Zhang et al., 2020), the relationship between surface O 3 concentrations and meteorological conditions Qiao et al., 2019), and the precursors and their contributions to O 3 (Zhang et al., 2008). However, most of these studies were restricted to the overall change characteristics of O 3 , there are few researches about mutation characteristic. Therefore, the wavelet analysis method was introduced. Compared with traditional methods, wavelet analysis can quickly find the mutation characteristics of time series. Moreover, the wavelet coherence spectrum can reflect both the correlation and time lag between variables which makes up the shortage of the traditional correlation analysis method (such as Pearson, Spearman and Kendall). Besides, traffic is a significant source of O 3 (Colvile et al., 2001), but very limited studies have analyzed the relationship between O 3 concentrations and traffic parameters (such as vehicle flow, vehicle speed and traffic flow density). Finally, researches on O 3 have given priority to the developed areas such as BTH, PRD, and YRD, but there are relatively few studies on O 3 in the second tier and third tier cities, whose O 3 pollution is prominent day by day.
Fuzhou is the capital of Fujian Province, a second tier city, one of the Chinese open coastal port cities and one of the starting cities of the Maritime Silk Road. With the fast growth of the economy and society, the O 3 pollution has attracted the attention of the Fuzhou municipal government and the public. From 2015 to 2018, the average annual concentration of O 3 has increased from 79 to 96 µg m -3 , approaching Grade I of the national standard. To cope with this, the environmental protection departments issued many policies, such as the O 3 Pollution Prevention and Control Plan (2018), the Blue Sky Action Plan of Fuzhou and "Accelerate elimination of vehicles with yellow license plates" . However, there are very limited studies on O 3 pollution in Fuzhou.
In this study, the time series of O 3 concentrations was decomposed and reconstructed by wavelet transform, and the variation characteristics of O 3 concentrations in Fuzhou were analyzed, to account for the fluctuations of O 3 concentrations. In this paper the correlation between O 3 concentrations and main influence factors such as meteorological conditions, gaseous precursors and traffic condition based on wavelet coherence were further analyzed. Consequently, the results are helpful for government departments to formulate pollution prevention strategies based on understanding the O 3 variation characteristics and influence factors. In addition, the analysis on traffic parameters and O 3 concentrations can provide a basis for the formulation of transportation management policy, including transportation demand management (TDM) and transportation system management (TSM), to achieve emission reduction goals.

Study Sites
Fuzhou is located in the eastern coast of Fujian Province, at 26°4ʹ48ʺN latitude, 119°18ʹE longitude with an average elevation of 80 m. It covers an area of 11,968 km 2 , and the urban area is 1,219.37 km 2 . The landform of Fuzhou is estuary basin (roughly similar to Los Angeles, CA, USA), and the existing forest cover percentage has reached 57.26% (up to March 2019), which occupies second place in provincial capital cities.
In recent years, with the development of economy and acceleration of urbanization, Fuzhou has become a modern city and the central city of West Taiwan Strait City Belt. The population of Fuzhou has increased continuously from 6.74 million in 2014 to 7.74 million in 2018. The car ownership had reached 1.55 million at the end of 2018, increasing by 31.2% over the end of 2017. The industrial energy consumption was 15.36 million tons of standard coal and the electricity consumption increased by 21.7% in 2018 (Bureau of Statistics of Fuzhou, 2018). Due to the rapid urbanization, the resources and environmental problems become increasingly prominent.
As to the climate conditions, the climate in Fuzhou is subtropical monsoon climate. The annual average sunshine duration varies from 1,700 hours to 1,980 hours, the annual average rainfall varies from 900 mm to 2,100 mm. The annual average temperature ranges from 20°C to 25°C, the predominant wind direction is the northeaster and partial south wind is boosted up during the period from June to August .
Currently, Fuzhou has six national air quality monitoring stations (Fig. 1), the spatial distribution, monitored frequency, and so on are illustrated in our previous study .

Data Sources
For this study, the daily mean data of air pollutant (includes NO 2 , CO and O 3 ) corresponding to the period 2017 to 2018 have been collected from the website of the China National Environmental Monitoring Centre (CNEMC) (http://www.cnemc.cn). According to the requirements of the National Ambient Air Quality Standards (GB 3095-2012), the maximum daily 8-hour average (MDA8) O 3 concentration was used as a criterion for evaluating the O 3 pollution level, and the value of primary standards (100 µg m -3 ) was used to judge whether the O 3 concentration exceeded limits. Moreover, the VOCs had a certain influence on the process of photochemical reaction, but they were not monitored by CNEMC, so VOC data were not used in this study.
The meteorological parameters data, including temperature, sunshine duration, relative humidity, precipitation, air pressure and wind speed, were derived from the website of National Meteorological Center (NMC) (http://www.nmc.cn).
In addition, to analyze the relationship between O 3 concentration and traffic parameters, the vehicle flow data were acquired from video monitoring, and the vehicle speed data were derived from laser Doppler anemometry (LDA) devices installed on the monitoring road segment.

Wavelet Transform The Continuous Wavelet Transform
To analyze the mutation characteristic, the time series of O 3 concentration was decomposed and reconstructed by the continuous wavelet transform (CWT). The CWT is defined as: where WT f (a,τ) is wavelet coefficients, a is scale parameter, τ is translation parameter, t is time and ψ a,τ (t) is the wavelet basis function.
Considering the characteristics of data and the analysis request, the Morlet and Daubechies were used in this paper.

Morlet Wavelet
The Morlet wavelet is defined as: where ω 0 is dimensionless frequency and η is dimensionless time.

Daubechies Wavelet
The series of Daubechies wavelet is abbreviated as "dbN," where N is order. The Daubechies wavelet has no analytic expression, in addition to db1 (which equal to Haar) and it is defined as:

The Cross Wavelet Transform
The relationship between two time series (O 3 concentration and meteorological parameters) can be found by analyzing the time-frequency domain over multiple time scales.
Defining W x (s) and W y (s) as the wavelet transform of two time series (X and Y, s means time), and the cross wavelet transform (XWT) can be defined as: The density of cross wavelet spectrum is defined as |W n XY (S)| (n means the length of time series), which means that the two time series have a common wavelet energy, and the correlation becomes more significant as this value increases.
The cross wavelet spectrum can be verified by comparison with the red noise standard spectrum. If the red noise significance test is passed, the two sequences are considered to have a good correlation, otherwise the correlation is not significant. The red noise standard spectrum of two time series with P k X and P k Y is given, and the theoretical distribution of the cross wavelet power can be defined as: where σ X and σ Y are the standard deviation of two time series, Z v (p) is the confidence level, v is the degree of freedom using Morlet wavelet. For example, the Z V (0.95) = 3.99 is calculated in the 5% significance level.

Cross Wavelet Phase Angle
The O 3 concentration is not completely synchronized with the variation of meteorological parameters, so the cross wavelet phase angle is introduced to describe the time lag of time series. Within the cone of influence (abbreviated as "COI," and the confidence is 95%), the phase angle can quantitatively describe the phase relationship. The mean phase angle of a set of angles a i (i = 1, …, n) is defined as:

Wavelet Coherence
The wavelet coherence (WTC) spectrum is used to analyze the correlation between O 3 concentration and meteorological parameters. The wavelet coherence spectrum can be illustrated as: where S(W) is a smoothing operator, S scale denotes smoothing along the wavelet scale axis and S time denotes smoothing in time. A suitable smoothing operator for the Morlet wavelet is given as follows: where Π is the rectangle function, c 1 and c 2 are normalization constants, and the factor of 0.6 is the empirical scale for the Morlet wavelet. In addition, the significance test of wavelet coherence spectrum was performed by Monte Carlo method in this study. And the mean phase angle and 95% confidence interval drawn in this paper are for wavelet coherence spectrum, and only the arrows with R n 2 ≥ 0.5 are marked in wavelet coherence spectrum.

Traffic Parameters
To explore the relationship between traffic and O 3 pollution, the traffic parameters were introduced in this study. In traffic engineering, vehicle flow, vehicle speed and traffic flow density are the basic parameters of traffic flow characteristics. The basic relation between parameters is: where Q is vehicle flow, V is vehicle speed, and K denotes traffic flow density.

Characteristic Analysis Annual Variation Characteristic
According to the trend between O 3 concentration and meteorological parameters, one year is divided into four seasons: spring (March-May), summer (June-August), autumn (September-November), and winter (December-February). The MDA8 O 3 concentrations in 2018 are shown in Fig. 2, ranging from 16 to 201 µg m -3 . Fig. 2 shows O 3 mass concentration variation has the characteristics of instability, nonlinearity and multiple time scales, which are higher in spring and summer, lower in autumn and winter, ranging from 20 to 201 µg m -3 and from 16 to 176 µg m -3 , respectively. The O 3 concentration exhibits multiple-peak distributions and averages 96 µg m -3 (slightly lower than the Grade I national standards for air quality). More specifically, the mean O 3 concentrations during spring, summer, autumn and winter are 118 µg m -3 , 104 µg m -3 , 92 µg m -3 and 69 µg m -3 , respectively, and the numbers of days exceeding the value of primary standards are 65 d (accounted for 41%), 51 d (accounted for 31%), 37 d (accounted for 24%) and 7 d (accounted for 4%), respectively.
Such an asymmetric distribution may be explained by various reasons. Firstly, Fuzhou was affected by the monsoon climate, the average daily temperature and light intensity reached higher value, which has proved to be advantageous for photochemical reaction. However, both temperature and light intensity were lower in winter, so that the photochemical reaction rate maintained at a low level. Secondly, there was a development process of atmospheric convection, increased gradually in spring and early summer, reached the peak in summer, slowed down after summer, the tropopause folding often induces the exchange of stratospheric and tropospheric when atmospheric convection reached its peak, and O 3 was transported vertically downward to near the ground which caused an increase in near-surface ozone concentration. In addition, lower wind speed, less rainfall and weaker sunshine duration in autumn and winter can leads to the poor atmospheric diffusion conditions. Furthermore, Fig. 2 indicates significant fluctuations of O 3 concentration on following months: March-April, June-July and September. This phenomenon may be caused by the following reasons. To mention first, the rainy season in Fuzhou was mainly concentrated March-April and June-July, and the fluctuations in September were mainly affected by typhoon weather. Less solar radiation and lower temperature during the rainy or cloudy days were not conducive to the process of photochemical reaction, and lower concentration of NO x , VOCs and other precursors results in the decrease of O 3 generation.

Mutation Characteristic
The variation characteristic in mutation days may contain more information related to influencing factors, so the wavelet transform was used to analyze further. The db1 wavelet has good regularity and it can detect the mutation points of time series accurately, so db1 wavelet was selected to conduct decomposition (decomposition level: 4 levels) for the concentration time series (Weng et al., 2011).  Table 1.
The mutation of O 3 concentration was the role of comprehensive factors, mainly including human activities and meteorological factors (Steiner et al., 2002). The ozone precursors (NO x , CO and VOCs) from human activities can be converted into O 3 by photochemical reactions. Furthermore, meteorological factors also have an important impact on the mutation of O 3 concentration. The mutation points shown in Table 1 are the drop points of O 3 concentration, and it can be inferred that lower temperature, larger precipitation, shorter sunshine duration and higher relative humidity may reduce O 3 concentration.

Correlation Analysis Meteorological Parameters
In order to explore the influence of meteorological factors on O 3 variation, the wavelet correlation spectrum was used based on O 3 pollutant concentration and meteorological data in 2018. And the results are shown in Fig. 4.

Temperature
Fig. 4(a) presents the wavelet correlation spectrum between O 3 concentration and daily mean temperature (DMT). As shown in figure, the disorder of phase characteristics between O 3 concentration and DMT was increased as the period becomes higher, and the wavelet energy mainly concentrated in the period range of 0 d to 12 d. With more particular reference to Fig. 4(a), the value of phase difference in May is 45 ± 15° at 0 d to 6 d scale which indicates the changes of O 3 concentration lags behind DMT 1/8 ± 1/24 period (about 0.5 d), but the phase difference is small from August to September at 2 d to 6 d scale and the O 3 concentration has a positive correlation with DMT (correlation coefficient > 0.7), additionally, there is no obvious time lag relationship from late February to April at 4 d to 8 d scale.
The relationship between O 3 concentration and DMT may be further verified by analyzing Fig. 4(b). The mean value and median of O 3 concentration increased gradually with the growth of the temperature, lower concentration mainly concentrated in the range of T ≤ 15°C and the highest concentration appeared when 20°C < T ≤ 0°C. Moreover, the over-standard rate of O 3 increased obviously as temperature growing and reached a peak at T > 30°C, and there was a dramatic growth when 10°C < T ≤ 25°C.
The results revealed that temperature has a positive correlation with O 3 , higher temperature accelerates the speed of photochemical reactions, causing an increase in the amount of O 3 .  992.1 0 WD = wind direction; RH = relative humidity; WS = wind speed; P = Precipitation; T = Temperature; AP = air pressure; SD = sunshine duration. Fig. 4(c) shows O 3 concentration has positive significantly correlated relationship to the sunshine duration (SD; Fig. 4(c)). Looking more closely, wavelet energy distributed in most of the period, a large amount of wavelet energy appeared in June to December and the phase difference was 0 ± 25°, O 3 concentration and sunshine duration existed significant synergistic changes.

Solar Radiation
Like temperature, the O 3 concentration presented synchronous increase trend with sunshine duration (Fig. 4(d)). In particular, as the sunshine duration grew, the mean value, the first quartile and the third quartile of O 3 concentration had experienced process of continuous increased (0 h < SD ≤ 8 h) and changed smoothly (8 h < SD). Furthermore, the over-standard rate of O 3 had a sharp rise from 21.30% (0 h ≤ SD ≤ 2 h) to 84.21% (10 h < SD).
The change of temperature was closely correlated with the solar radiation, and the increase of sunshine duration is conducive to the accumulation of solar radiation energy. The rise of ambient temperature and solar energy can promote the photochemical reaction, accelerate the rate of O 3 generation. Besides, the difference in sunshine duration was an important reason leading to the seasonal fluctuation of O 3 variation. Fig. 4(e) shows O 3 concentration is significantly and negatively correlated with the daily mean relative humidity (DMRH). The energy regions with a correlation coefficient is greater than 0.8 occupies about 2/3 of the figure. In particular, the distribution of wavelet energy was scattered and discontinuous at 0 d to 8 d scale, but energy distributed comparatively centralized at 0d to 8d and > 64 d scale and the phase difference was 0 ± 25°, which indicated that the O 3 concentration varied in a phase opposite to DMRH and there was a significantly negative correlation.

Relative Humidity
As can be seen in Fig. 4(f), the negative correlation between O 3 concentration and relative humidity is further confirmed. Particularly, as relative humidity (RH; Fig. 4(f)), the lower mean value of O 3 concentration occurred, and RH = 70% was a dividing point of the large fluctuations, and higher O 3 concentration was mainly concentrated in the range of RH = 40-70%. Moreover, the exceeding rate of O 3 was in a tendency of decrease (RH ≥ 70%) and the rate reached zero when RH ≥ 90%.
The influence of relative humidity on O 3 concentration is mainly represented in the following aspects: First, solar radiation will be attenuated by extinction mechanism under the action of water vapor, specifically, a large amount of water vapor tends to form clouds, thereby inhibiting ultraviolet rays reaching near the ground and delaying the occurrence of photochemical reactions. Second, in the case of high humidity, the free radicals H, OH, etc. contained in the water vapor rapidly decompose O 3 into oxygen molecules, resulting in the reduction of O 3 concentration (Wennberg et al., 1998). Third, high humidity can decrease O 3 concentration by promoting the dry deposition (a process of atmospheric aerosol particles and trace gas components settling in the absence of precipitation).

Precipitation
In general, the wavelet coherence spectrum showed that a negative correlation existed between O 3 concentration and daily mean precipitation (DMP) (Fig. 4(g)). It is noteworthy that the negative correlation was mainly reflected in spring and summer at 0 d to 16 d scale (correlation coefficient > 0.7). Specifically, wavelet energy appeared in February to June and the phase difference was 180 ± 20° at 0 d to 16 d scale, which indicated the changes of O 3 concentration lagged behind DMP 1/2 ± 1/12 period (about 1-6 d), and the phase difference was 180 ± 25° at 48 d to 64 d scale, which showed that relations between O 3 concentration and DMP were obvious antiphase changed.
As shown in Fig. 4(h), along with the precipitation (PR; Fig. 4(h)) increases, O 3 concentration shows decreasing trend. The highest value of O 3 concentration appears in the range of PR ≤ 3 mm, and the mean value of O 3 concentration decreases continuously (PR ≤ 9 mm) first and then declines slightly (9 mm < PR). The exceeding rate of O 3 decreased dramatically firstly (PR ≤ 9 mm) and then increases slightly (9 mm < PR ≤ 12 mm), finally, dropped to low levels (12 mm < PR). It is obvious that precipitation has a significant elimination effect on ozone at PR ≤ 9 mm.
The precipitation has a certain influence on the O 3 production process. In one aspect, the cloud layer is thicker during precipitation, which reduces the sunshine duration and temperature, and cloud layer scatters and absorbs solar radiation reaching the surface layer. In addition, precipitation has the effect of wet deposition, so that the O 3 concentration on precipitation day is lower than that without precipitation. Fig. 4. The wavelet coherence spectrum (left) and the distribution of meteorological parameters (right). In the spectrum, the 5% significance level against red noise is shown as a thick contour, the relative phase relationship is shown as arrows (with in-phase pointing right, antiphase pointing left), and the orange represents a larger correlation coefficient and blue represents a smaller correlation coefficient.

Air Pressure
Overall, Fig. 4(i) shows a low correlation between O 3 concentration and daily mean air pressure (DMAP). The distribution of the wavelet energy region with correlation coefficient greater than 0.7 is scattered and the proportion of remarkable wavelet energy is small, therefore, no obvious correlation can be seen from the wavelet coherence spectrum. However, it can be found in Fig. 4(j) that there is a negative correlation existed between O 3 concentration and air pressure (AP; Fig. 4(j)). With the increase of air pressure, the mean value, the quartile and the exceeding rate of O 3 concentration decreases at first (AP ≤ 995 hPa), then increases (995 hPa < AP ≤ 1000 hPa), and decreased finally (1000 hPa < AP). It is noteworthy that the over-standard rate of O 3 decreases from 80% to 18.75%, which shows a significantly negative correlation between O 3 and air pressure.

Wind Speed
As can be seen in Fig. 4(k), significant correlation is not observed between O 3 concentration and wind speed, the result is similar to the discussion of O 3 concentration and DMAP.
Further studies were performed to explore the correlation between O 3 concentration and wind speed (WS; Fig. 4(l)). Fig. 4(l) shows that the highest value of O 3 concentration grows gradually first (0.5 m s -1 ≤ WS ≤ 1.9 m s -1 ) and decreased continually (1.9 m s -1 < WS). Likewise, the relationship between the exceeding rate of O 3 concentration and wind speed is not linear, decrease of the exceeding rate is obvious at 2.6 m s -1 < WS ≤ 4.0 m s -1 , but the exceeding rate increases constantly in the low wind speed interval (0.5 m s -1 ≤ WS ≤ 2.6 m s -1 ). In addition, the exceeding rate is maintained higher level at above 30% in all wind speed intervals.
The wind speed can not only reflect the efficiency of pollutant transport and atmospheric diffusion, but also reflect the intensity of stability in the atmospheric boundary layer to some extent. The change of O 3 concentration with the increase of wind speed is mainly reflected in two aspects: First, as wind speed increases, the height of atmospheric boundary layer is elevated, the vertical momentum transport is strengthened, and then the O 3 in the troposphere is transmitted to the ground layer. Second, the horizontal diffusion of O 3 is enhanced. Moreover, when the wind speed is low (0.5 m s -1 ≤ WS ≤ 2.6 m s -1 ), the horizontal diffusion of O 3 is weaker than its downward transport effect, resulting in an increase of O 3 concentration near the ground, but when the wind speed exceeds a certain value (2.6 m s -1 < WS ≤ 4.0 m s -1 ), the horizontal diffusion gradually dominates, and the O 3 concentration will decrease significantly with the increase of wind speed. Besides, the O 3 concentration fluctuation interval (4 m s -1 < WS) may be attributable to the stronger downward transport effect than the horizontal diffusion.

Gaseous Precursor and Traffic Parameters
The NO x and CO contained in motor vehicle exhaust are precursors of O 3 , and exhaust emissions are closely related to real-time traffic conditions. Therefore, the precursors and traffic parameters were combined in this study. In addition, with the consideration of data acquisition, the pollutants and traffic hourly data (2017) were used, based on fitting and verifying the annual data of six monitoring stations, the data of typical working day (Wednesday, April 19, 2017) and rest day (Sunday, April 23, 2017) from Fujian Normal University monitoring station were taken as examples for analysis.
As can be seen from Fig. 5, the concentration of NO x and CO changes with the similar tendency, and there is a significant negative correlation between the precursors and O 3 on working day and rest day. On the working day, the initial concentrations of NO x and CO are higher, both NO x and CO drop gradually from 00:00 to 14:00, the valley values appear at 14:00 to 17:00, the fluctuation time periods of NO x and CO are mainly concentrated at 05:00 to 10:00 and 17:00 to 19:00, moreover, the O 3 concentrations remains steady at a low level before 07:00, increases sharply after 07:00, peaks in the time period from 12:00 to 13:00, and continues to decline after 16:00. In contrast, on the rest day, the initial concentration of NO x and CO are lower than working day, and there is a significant fluctuation from 00:00 to 07:00. Additionally, unlike the working day, the O 3 concentration on rest day gradually decreases from 00:00 to 04:00 and reaches the lowest value (53 µg m -3 ) at 04:00, increases continuously (04:00-12:00), and then maintains at a higher level (12:00-18:00), finally decreases gradually after 17:00.
As shown in Fig. 5, traffic parameters are closely related to O 3 concentration, which also reflects the influence of real-time traffic conditions on the O 3 concentration variation. On the working day, the vehicle flow (Q) and traffic flow density (K) are at a low value before 07:00, the morning peak is from 07:00 to 09:00 and the vehicle flow in this period are greater than 1,000 veh h -1 . According to the hydrodynamics simulation in traffic flow theory, when the traffic density increases due to the change of traffic conditions, the traffic flow wave will be generated in the traffic flow, and the traffic flow will experience the process of congestion to dissipation, therefore, the dissipating process of traffic flow occurred in both morning peak and evening peak (17:00-19:00). The vehicle speed (V) decreases as vehicle flow and traffic flow density increases, reaching a lower value in traffic peak period. The trend of O 3 concentration shows a unimodal shape, and the higher value appears in 12:00 to 16:00. In contrast, on rest day, the morning peak and evening peak also exists, but the vehicle flow and traffic flow density are relatively small, vehicle flow between 07:00 to 19:00 maintains above 1,000 veh h -1 and there are multiple fluctuation intervals of vehicle speed and traffic flow density, this phenomenon may be attributable to the decrease of commuting and the increase of leisure travel on rest day, which results in a discrete traffic flow distribution. And the trend of O 3 concentration on rest day displays roughly the same tendency with that on working day.
Under the influence of human activities, compared with working day, the concentration of precursors decreased at rest day, while the O 3 concentration increased, and this phenomenon is called "weekend effect" (Cai et al., 2010). In this study, the weekend effect was obvious, in particular, the hourly mean value of NO x concentration, CO concentration and vehicle flow decreased on rest day, 91 to 30 µg m -3 , 1.10 to 0.71 mg m -3 and 824 to 778 veh h -1 , respectively, while O 3 concentration increased from 49 to 123 µg m -3 . This result was related to resident trip characteristics on the weekend, different travel patterns caused the change of traffic condition.
According to the traffic flow theory, the traffic flow is in a uninterrupted traffic state and the vehicle speed is high when the road is unobstructed, as the vehicle flow increases, the traffic flow density increases accordingly, the traffic flow turns to interrupted traffic when the vehicle flow reaches a certain value, and the phenomenon of car following occurs (Verhoef, 2005). It is known from the formation mechanism of vehicle exhaust that frequent acceleration, deceleration and idling state will result in incomplete combustion of fuel and the engine speed change (Chan et al., 1991), which leads to an increase in NO x and CO emissions, therefore, the concentration of NO x and CO fluctuated during morning and evening peak. Moreover, combined with the discussion in the previous section, the conditions of strong solar radiation, high temperature and low relative humidity can promote the photochemical reaction, which also explains that NO x and CO were consumed from 09:00 to 17:00 and the precursor concentrations increased at night. Besides, the change of precursors at night was closely related to freight vehicles. Due to policy restrictions, the freight vehicles can only drive at a fixed time during the night, and the emission factor of freight vehicles (such as high tonnage diesel vehicles) is much higher than ordinary motor vehicles.

CONCLUSIONS
Our results, which confirm the advantages of analyzing time series of pollutant concentrations via wavelet transforms, indicate that meteorological factors, gaseous precursor levels and traffic conditions all contribute to the concentration of O 3 . The Ministry of Ecology and Environment should focus on the changes in O 3 concentration during spring and summer and treat the control of precursors as the core of air pollution management, including vehicular emission reduction. However, the accuracy of our assessment may be hampered by our limited data. Hence, additional cities must be sampled and analyzed in the future, and a vehicular emission database that monitors the entire urban road network should be established in order to determine the relationship between traffic conditions and O 3 pollution.