Influence of Atmospheric Particulate Matter on Ozone in Nanjing,China:Observational Study and Mechanistic Analysis

Yawei QU,Tijian WANG?,Yanfeng CAI,Shekou WANG,Pulong CHEN,Shu LI,Mengmeng LI,Cheng YUAN,Jing WANG,and Shaocai XU

1School of Atmospheric Sciences,CMA-NJU Joint Laboratory for Climate Prediction Studies,Jiangsu Collaborative Innovation Center for Climate Change,Nanjing University,Nanjing 210023,China

2Guangdong Electric Power Design Institution Co.,Ltd.of China Energy Engineering Group,Guangzhou 510663,China

3Jiangsu Academy of Environmental Industry and Technology Crop,Nanjing 210036,China

4Qingdao Environmental Monitoring Center,Qingdao 266003,China

ABSTRACT Particulate matter with diameters of 2.5μm or smaller(PM2.5)and ozone(O3)are major pollutants in the urban atmosphere.PM2.5can affect O3by altering the photolysis rate and heterogeneous reactions.However,these two processes and their relative importance remain uncertain.In this paper,with Nanjing in China as the target city,we investigate the characteristics and mechanism of interactions between particles and O3based on ground observations and numerical modeling.In 2008,the average concentrations of PM2.5and O3at Caochangmen station are 64.6±47.4μg m?3and 24.6±22.8 ppb,respectively,while at Pukou station they are 94.1±63.4μg m?3and 16.9±14.9 ppb.The correlation coefficient between PM2.5and O3is?0.46.In order to understand the reaction between PM2.5and O3,we construct a box model,in which an aerosol optical property model,ultraviolet radiation model,gas phase chemistry model,and heterogeneous chemistry model,are coupled.The model is employed to investigate the relative contribution of the aforementioned two processes,which vary under different particle concentrations,scattering capability and VOCs/NOxratios(VOCs:volatile organic compounds;NOx:nitric oxide and nitrogen dioxide).Generally,photolysis rate effect can cause a greater O3reduction when the particle concentrations are higher,while heterogeneous reactions dominate O3reduction with low-level particle concentrations.Moreover,in typical VOC-sensitive regions,O3can even be increased by heterogeneous reactions.In Nanjing,both processes lead to O3reduction,and photolysis rate effect is dominant.Our study underscores the importance of photolysis rate effect and heterogeneous reactions for O3,and such interaction processes should be fully considered in future atmospheric chemistry modeling.

Key words:PM,ozone,photolysis,heterogeneous reaction,Nanjing,urban atmosphere

1.Introduction

Air pollution has become increasingly significant in Nanjing,China,because of its accelerating urbanization and industrialization(Xie et al.,2016;Chen et al.,2017;Zhu et al.,2017).As the major atmospheric pollutants in the urban area,particulate matter(PM)and ozone(O3)have drawn much attention in recent years for their adverse effects on human health(Kan et al.,2008;Xie et al.,2017).Particles and O3are short-lived and radiatively active chemicals(Liao et al.,2015;Park et al.,2016).The chemical interaction between O3and PM is complicated.O3is believed to have an influence on PM concentrations.O3can act as an oxidant,as well as have an impact on the concentration of hydroxyl(OH)radicals and affect the formation of sulfate(),nitrateammonium(),secondary organic aerosols,etc.(Yang et al.,2014;Shi et al.,2015;Li et al.,2017).PM can strengthen the extinction capacity of the atmosphere directly,and it can also act as the cloud condensation nuclei that influence radiation indirectly through affecting cloud formation.Therefore,PM can influence the photolysis rate of O3and its precursors by changing solar radiation(Cai et al.,2013;Nishanth et al.,2014;Xie et al.,2014).PM can also influence the transportation and diffusion of O3and its precursors through changing the atmospheric stability(Chan et al.,1999;Velasco et al.,2008;Qu et al.,2017).Moreover,the heterogeneous reactions of chemicals,such as NOxand dinitrogen pentoxide(N2O5),on particle surfaces,also exert important influences on the generation and consumption of O3(Jacob,2000;Lou et al.,2014;Akimoto,2016).However,the mechanism and combined effect of particles on both photolysis rates and heterogeneous reactions are not fully understood.

One of the important pathways through which PM can influence O3is by altering the photolysis rate.The scattering and absorbing effect of PM can change the optical properties of the atmosphere and have an impact on the solar radiation that reaches the ground,leading to O3variation through altering the photolysis rate and photochemical reactions(Meng et al.,1997;Nishanth et al.,2014).Existing studies generally agree that a reduction in O3concentration and productivity is closely related to serious atmospheric PM pollution.Liao et al.(1999)used a radiation transfer model,DISORT,to analyze the impact of aerosols on photochemistry.Tie et al.(2001,2005)assessed the effects of aerosols on tropospheric oxidants by using a three-dimensional global chemical/transport aerosol model,MOZART,in which photochemical reactions were implemented.They found that surface absorbent particles reduced the photolysis rates and therefore weakened O3generation and accumulation,while scattering aerosols did the opposite.Deng et al.(2008,2011,2012)compared observations with box model results and analyzed the correlation between PM10(PM with diameters between 2.5 and 10μm),PM2.5(PM with diameters of 2.5μm or smaller),O3,and aerosol optical depth(AOD),in Guangzhou,China.The correlation coefficient between AOD and PM10mass concentration was very high,with a maximum of 0.98,and the AOD and UV radiation/O3was anti-correlated,with a correlation coefficient of 0.90.Li et al.(2011a)used a regional chemical transport model coupled with a radiative transfer model to examine the effects of PM on photochemistry in central East China.They found that summertime O3in the boundary layer(below 1 km and 1–3 km)was reduced by 5%,with a maximum of 9%,in highly polluted regions.

Heterogeneous reactions on the surface of PM2.5can also affect the concentration of O3and secondary particles.Ravishankara(1997)described the role,rates and media of heterogeneous and multiphase reactions in the troposphere,and it was clear from their findings that heterogeneous and multiphase chemistry have the potential to alter the composition of the troposphere.Jacob(2000)proposed that heterogeneous reactions have an impact on O3concentrations in a number of ways,including the generation and consumption of NOx,N2O5,HOx[OH+hydroperoxyl(HO2)]and O3,and the generation of halogen radicals.Moreover,according to previous studies(Remorov et al.,2002;Hoffman et al.,2003;Bauer et al.,2004;McNaughton et al.,2009),heterogeneous reactions on the surface of sea salt and soot have significant effects on trace gases.Heterogeneous reactions can also be accelerated by black carbon(BC),which is mainly composed of PM2.5and has a relatively large surface area,leading to changes in the concentrations of gases and particles(Kleffmann and Wiesen,2005;Kaiser et al.,2011).Deng et al.(2010)applied a box model to study the interactions between aerosols and O3.Their results showed that the impact of heterogeneous reactions strongly depends on the aerosol concentration and the surface uptake coefficients.

Only a few studies have taken into account the impact of PM on both photochemical and heterogeneous reactions,and were limited to modeling and large-scale observation(Tie et al.,2001,2005;Bian et al.,2003;Lou et al.,2014).In contrast,the present paper focuses mainly on the mechanism of the impact of PM on O3generation and consumption in Nanjing,China.By using field observational data and box modeling,we investigate how PM influences O3concentrations through two ways—altering the photolysis rate and affecting heterogeneous reactions on particle surfaces.Section 2 describes the observational data and the box model.The observational data in Nanjing for the whole of 2008 and a typical case are analyzed in section 3.In section 4,we discuss the variations of O3and its precursors and secondary aerosols based on the results from a set of sensitivity experiments,in which we consider the impact of PM on photolysis rates and heterogeneous reactions.Finally,conclusions and perspectives on future work are given in section 5.

2.Methods

2.1.Observation

The measurement data used in this study were collected at two stations.The first of these was Pukou(PK)air monitoring station(32°18′N,118°47′E),on Pukou Campus of Nanjing University,in the northwest suburbs of Nanjing.The Industrial Development Zone along the Yangtze River is about 6 km to the northeast of Pukou,and thus this site can be regarded as an industrially polluted suburb.The second one was Caochangmen(CCM)station(32°18′N,118°26′E),which is located in the urban district and more influenced by human activities.The locations of PK and CCM are shown on the map of Nanjing in Fig.1.Chemical species are monitored continuously by instruments including an O3analyzer(O342M,Environment S.A.,Paris,France)forO3,anammonia(NH3)and NOxanalyzer(AC32M-CNH3,Environment S.A.,Paris,France)for NOx,and a PM analyzer(MP101M,Environment S.A.,Paris,France)for PM.The analyzers were calibrated daily and observations carried out continuously from October 2007 to August 2009.AOD retrievals were collected from sunphotometer AOD measurements at the Nanjing University of Information Science&Technology site(32.2°N,118.72°E)of NASA’s Stennis Aerosol Robotic Network.In this paper,PM2.5,O3,NOxand AOD data collected from 1 January to 31 December 2008 are analyzed.

2.2.Box model

A zero-dimensional box model,which includes the calculation of aerosol optical properties,UV radiation transmission processes,atmospheric photochemical processes,heterogeneous chemical processes,and aerosol thermodynamic equilibrium processes, was constructed and used in this study.

Fig.1.Map of Nanjing,China,and the locations of the observational sites of Pukou(PK,32°18′N,118°47′E)and Caochangmen(CCM,32°18′N,118°26′E).

The OPAC(Optical Properties of Aerosols and Clouds)software package,developed by the University of Munich,Germany,and the Max Planck Institute of Meteorology(Hess et al.,1998),was used as the optical property module.It provides optical properties including extinction,scattering and absorption coefficients,single scattering albedo(SSA),asymmetry parameter,and phase function.These aerosol optical properties are given for up to 61 wavelengths between 0.25 and 40μm and eight values of relative humidity.For further details of the OPAC database,please refer to Hess et al.(1998).The UV radiation transmission module was the Tropospheric Ultraviolet and Visible(TUV)radiation model from the National Center for Atmospheric Research(Madronich et al.,1999).TUV was used for calculating the spectral irradiance,spectral actinic flux,and photo dissociation coefficients,over the wavelength range 121–750 nm.The carbon bond chemical mechanism known as CB05(Yarwood et al.,2005)was applied for the photochemical reactions.This mechanism is a comprehensive upgrade of the CB-IV carbon bond mechanism and contains 51 species and 156 reactions(including 23 photolysis reactions).The aerosol chemistry module was the balancing module known as ISORROPIA,developed by Nenes et al.(1998).This module describes the equilibrium of different phase species,including NH3,hydrogen chloride,nitric acid(HNO3)gas,solid and liquidand chloride.The heterogeneous chemistry module contains ten reactions,four of which occur on the surface of aqueous aerosols(Jacob,2000),and the others on BC surfaces(Deng et al.,2010).The first-order rate constant k for heterogeneous loss of a gas to the aerosol is given by

in which r is the particle radius,Dgis the gas-phase molecular diffusion coefficient of species X in air,m is the mean molecular speed of X in the gas phase,γ is the absorption coefficient on particles,and A is the aerosol surface area per unit volume of air(Jacob,2000).Heterogeneous reactions and uptake coefficients on BC and aqueous aerosols are listed in Table 1.

Table 1.Heterogeneous reactions and uptake coefficients on BC(Deng et al.,2010)and aqueous aerosols(Jacob,2000).

The initial PM2.5concentration was set to 105 μg m?3,based on observations.In addition,according to previous works on PM in Nanjing(Wang et al.,2002;Yang et al.,2005;Yin et al.,2009;Chen et al.,2010;Tan et al.,2010),the proportion of PM components was 16%for soluble particles,25%for15%forand35%for soluble organic particles,and 9%for BC.All these results were applied as our initial proportion settings.As for size distributions,lognormal distributions were used for each chemical species,more details on which can be found in Hess et al.(1998).The initial NOxwas set to 40 ppb,based on observation in Nanjing.The total volume concentration of volatile organic compounds(VOCs)was 70 ppb at the beginning,and the settings for the proportions of different kinds of VOCs are shown in Table 2(Liu et al.,2008;Wang and Zhao,2008).To verify the reasonability of this setup,the SSA was calculated with this box model.The result of 0.76 is close to the SSA observed by Chen et al.(2009),which was 0.75.

3.Observational analysis

3.1.Statistical characteristics

Based on hourly data from 1 January to 31 December 2008,the annual mean,standard deviation,median and maximum of PM2.5and O3concentrations were calculated,and the results are shown in Table 3.The annual mean and standard deviation values for PM2.5and O3at PK station are 94.09±63.42μg m?3and 16.94±14.91 ppb,respectively.Their counterparts at CCM station are 64.63±47.41μg m?3and 24.63±22.81 ppb.The PM2.5concentrations at both sites far exceed the national standard concentration,which is 35μg m?3for the annual mean,according to the China National Ambient Air Quality Standard(China,2012).Clearly,severe PM pollution exists both in the urban and suburban areas of Nanjing,probably originating from industry,vehicle emissions,and dust blown by winds above unpaved roads(Wang et al.,2012).The PK site shows relatively high PM2.5and low O3levels compared with CCM,suggesting severer particulate pollution in the suburbs of Nanjing.One possible explanation is that PK station is situated near an industrial area and has mountains on both its north and southwest side(Shao et al.,2016).The relatively higher levels of industrial emissions and worse diffusion conditions lead to the relatively higher PM2.5concentrations at PK,which may reduce the O3concentration.

Table 2.Composition of VOCs and their proportions in the boxmodel setup.The numbers in the table represent the ratio between the number of functional groups and the number of VOC molecules;therefore,the total percentage is not 100%.

Table 3.Statistics(annual mean,median,standard deviation and maximum)of PM2.5(μg m?3)and O3(ppb)concentrations from CCM and PK stations in 2008.

3.2.Seasonal variation

Figure 2 illustrates the seasonal mean PM2.5and O3concentrations.The seasonal trends of O3are consistent at both sites,with a maximum in spring and minimum in winter.Solar radiation is higher in summer and spring compared with the rest of the year,but a large amount of precipitation in summer—especially during the mei-yu period—reduces the solar radiation and leads to a lower O3concentration in summer than in spring(Tong et al.,2017).Therefore,O3peaks in spring,with a value of 32.41 ppb and 21.77 ppb at CCM and PK,respectively.The PM2.5concentrations at CCM and PK in spring are 63.15 μg m?3and 92.37 μg m?3,respectively.In summer,the PM2.5concentration reduces to 53.33μg m?3at CCM but rises to a peak of 119.41 μg m?3at PK.One possible explanation is that most of the factories in Pukou are located to the southeast of PK station.Southeasterly wind,which is the prevailing wind in summer in Nanjing,will bring pollutants from the factories to PK station,leading to the maximum PM2.5(Xiao et al.,2017).The O3concentrations in summer are 29.07 ppb and 19.03 ppb at CCM and PK,respectively,and the O3concentration in autumn is similar to that in summer.In autumn,the PM2.5concentrations are 73.01μg m?3and 91.27μg m?3at CCM and PK,respectively.The PM2.5concentrations in winter are 77.87μg m?3at CCM and 86.15μg m?3at PK,and the O3concentrations at the two stations are both around 17 ppb.Comparing the two sites,regardless of season,lower O3concentrations will be accompanied by higher PM.O3concentrations are lower at PK than at CCM in all four seasons;whereas,PM2.5at PK is higher than that at CCM in each season.This phenomenon may be linked to the inverse relationship between PM2.5and O3.

3.3.Weekly variation

Fig.2.Seasonal variations of O3and PM2.5at CCM and PK stations from 1 January to 31 December 2008.

Fig.3.Weekly cycles of(a)PM2.5and(b)O3concentration anomalies at CCM and PK stations from 1 January to 31 December 2008.

The weekly variations of PM2.5and O3are demonstrated in Fig.3.PM2.5concentrations are lower at weekends than on weekdays,possibly because PM2.5is tightly linked with human activity,i.e.,the reduced vehicle and factory emissions may be responsible for the low PM2.5levels at weekends.Indeed,this “weekend effect”on PM2.5has been verified by several studies(e.g.,Qin et al.,2004;Jeon,2015).Previous studies on O3have shown that it is higher at weekends than at other times of the week,which is perhaps attributable to the nonlinear response of O3to reduced NOxconcentrations under VOC-sensitive regimes(Heuss et al.,2003;Gao and Niemeier,2007;Pierce et al.,2010).In our observations,the O3“weekend effect”is less significant than that of PM2.5,with higher O3only observed on Sundays.The weekly variation shows an inverse correlation between PM2.5and O3.At PK station,PM2.5rises from Monday to Thursday,while O3concentrations decline from Tuesday;and when the PM2.5concentration declines from Thursday to Sunday,O3shows the opposite trend.Observations at CCM show a similar anticorrelation between PM2.5and O3from Monday to Wednesday,as well as at weekends.

3.4.Daily variation

Figure 4 depicts the diurnal variations of O3and PM2.5in different seasons,which shows an opposite trend between PM2.5and O3.The latter shows a single-peak pattern at both stations.From the figure,the O3concentration is at a minimum in the early morning.Then,due to enhanced atmospheric photochemical activity,O3generates and accumulates rapidly before reaching a peak in the afternoon at around 1400 LST(local standard time).The O3concentration then declines because of destruction by NOxduring the night.Afternoon maxima and early morning minim a have also been found in previous studies(Khoder,2009;Tong et al.,2017).The strongest diurnal variation of O3occursin spring,and the weakest in winter;plus,the diurnal variation is more significant in the urban district compared with the suburban area.

The diurnal variation of PM2.5shows a double peak,with one peak appearing between 0700 and 0900 LST,and the other one between 1900 and 2200 LST.Essentially,this pattern is due to human activities and planetary boundary layer(PBL)mixing.With the high numbers of vehicles in the morning and evening rush hours,emissions of PM2.5grow.Boogaard et al.(2010)and Dos Santos-Juusela et al.(2013)found that the temporal variation of particle concentrations was higher at locations near to streets.The double-peak pattern is more pronounced at CCM,suggesting that this urban site is more influenced by traffic,whereas the suburban site(i.e.,PK)is mainly influenced by industry.The morning PM2.5peak in winter lags that in summer by two hours at PK station.The reason for this phenomenon is that the sun rises around two hours later in the winter(0700 LST)than in summer(0500 LST).Consequently,the morning rush hour and the diffusion of PBL mixing are delayed.During the daytime,stronger solar radiation is beneficial to the development of the PBL.The increased turbulence intensity in the PBL facilitates the diffusion of pollutants,leading to the minimum PM2.5concentration in the afternoon(Qu et al.,2017).The minimum PM2.5appears between 1300 and 1500 LST,during which period the O3concentration is usually at its maximum,due to solar radiation and photochemical production.

3.5.Case study

Negative correlation between PM and O3can be found in all of the variations analyzed above,i.e.,seasonal,weekly and daily.In order to eliminate the effects of meteorological conditions,we analyzed a specific case in which there was a steady meteorological field that lasted from 24–30 March 2008 at PK station.Figure 5 shows the time series of meteorological elements including atmospheric pressure,surface temperature,relative humidity,and surface wind velocity.The pressure remains steady from 24–28 March,except for a slightly lower value 26 March.After 29 March,the atmospheric pressure gradually rises to 1020 hPa because of a high-pressure system.The temperature hardly changes during 24–28 March,with the average temperature at around 12°and a maximum of 20°.Surface temperature declines after 29 March.The relative humidity also remains steady during 24–28 March,and then grows to 90%after 29 March because of rain.The wind velocity is less than 8 m s?1during the whole period,except on 28 March when the weather system changes significantly.According to the meteorological conditions,the whole period can be divided into two parts:a steady period from 24–28 March,and a rainy period from 29–30 March.

Fig.4.Diurnal variations of observed(a,c)O3and(b,d)PM2.5at(a,b)CCM and(c,d)PK stations in different seasons in 2008.

Fig.5.Observed(a)atmospheric pressure,(b)surface temperature,(c)relative humidity and(d)surface wind velocity at PK station from 29–30 March 2008.

Figure 6 shows the time series for PM2.5,O3and NOxconcentrations,as well as AOD,from 24–30 March 2008 at PK station.The daily average PM2.5rises from 47.3μg m?3on 24 March to 124.1 μg m?3on 27 March,and the AOD rises by 0.5(500 nm)and 0.8(380 nm)in the same period.However,O3decreases from 36.7 ppb to 30.7 ppb.On 29 and 30 March,negative correlation between PM and O3is apparent,during which time the PM2.5concentration drops dramatically and reaches a minimum at noon on 30 March.Meanwhile,O3climbs gradually until 1500 LST 30 March.The negative correlation between PM and O3can be observed under the different meteorological conditions of the steady period and rainy period.

Fig.6.Observed(a)AOD and(b–d)concentrations of(b)PM2.5,(c)O3and(d)NOxin Nanjing from 24–30 March 2008.

Fig.7.Relationship between(a)O3and NOx,and(b)O3and PM2.5,from 24–27 March 2008.

To further investigate the negative correlation between PM2.5and O3,observations on 24 March and 27 March were selected and the hourly concentrations of NOx,PM2.5and O3in the daytime(0700–1700 LST)plotted in a scatter diagram(Fig.7).The opposite relationship between O3and NOxindicates that Nanjing is in a VOC-sensitive region,according to the O3isopleths from Empirical Kinetics Modeling Approach(EKMA;Milford et al.,1989;Shiu et al.,2007).An et al.(2015)and Geng et al.(2009)found a similar phenomenon in the Yangtze River Delta region and verified Nanjing as a VOC-sensitive region.The impact of the VOCs/NOxratio on PM and O3interaction will be discussed later.PM2.5and O3concentrations are shown as negatively correlated in Fig.7.The reason is most likely associated with an enhanced atmospheric PM extinction capacity.An increasing concentration of atmospheric PM(especially fine particles)raises the AOD and weakens the photochemical formation rate of O3,therefore depressing the O3concentration(Li et al.,2011b;Deng et al.,2012;Nishanth et al.,2014).In addition,the particle surface area concentration rises correspondingly when the particle number concentration increases,therefore leading to more heterogeneous reactions taking place on particle surfaces(Deng et al.,2010;Kaiser et al.,2011).

4.Mechanistic analysis

Here,we employ the box model to assess the relative contribution of the photochemical and heterogeneous reactions to the interaction between particles and O3under various conditions.The simulation schemes are introduced in section 4.1.Section 4.2 explains the negative correlation between O3and PM in the observations at Nanjing.We also run three groups of sensitivity experiments,to investigate the impact of atmospheric scattering capability(section 4.3),particle concentration(section 4.4),and the VOCs/NOxratio(section 4.5)on photochemical and heterogeneous reactions.

4.1.Simulation schemes

In order to study the impact of particles on photolysis rates and heterogeneous reactions,a set of numerical experiments was designed(Table 4).A control test was carried out without the effect of PM on the photolysis rate or heterogeneous reactions,before the sensitivity runs.A simulation with the effect of aerosols on the photolysis rate but without heterogeneous reactions was set as Case 1.Then,to evaluate the impact of different heterogeneous reactions,four simulations(Case 2-1,2-2,2-3 and 2-4),which ignored the effect of aerosols on photolysis,were performed.Case 2-1 only considered heterogeneous reactions of OH,HO2and hydrogen peroxide(H2O2).Case 2-2 was the same,but with the heterogeneous reaction of O3added.Case 2-3 simply included the heterogeneous reaction of nitrogen dioxide(NO2),whereas Case 2-4 involved the heterogeneous reactions of NO?3,dinitrogen pentoxide(N2O5)and HNO3.Next,Case 3-1 and Case 3-2 were designed to investigate the change in photolysis rate in different AOD and SSA scenarios.Compared with Case 3-1,which represented the clean atmosphere,Case 3-2 was a polluted case with the addition of PM.In Case 3-2,the AOD was set as 1 or 2,and the SSA was 1,0.9 or 0.8.In order to estimate the relative importance of two reactions under different PM concentrations,Case 4,with the effect of PM on photolysis(Case 4-1)and heterogeneous reactions(Case 4-2),was carried out,with the PM2.5concentration varyingfrom 0 to 175 μg m?3.Finally,Case 5 considered photochemical reactions(Case 5-1)and heterogeneous reactions(Case 5-2)with VOCs varying from 0 to 200 ppb and NOxfrom 0 to 100 ppb,to compare the impact of particles on O3by affecting both kinds of reactions in the O3–NOx–VOCs system.

Table 4.Simulation schemes of sensitivity experiments:(a)Control test,Case 1 and Case 2(?1,2,3,4)applied in studying the impact of PM on O3through altering the photolysis rate and heterogeneous reactions;(b)Case 3(?1,2),Case 4(?1,2)and Case 5(?1,2)applied in studying the impact of PM on O3under different scenarios,including different AOD and SSA,different PM2.5concentrations,and different VOCs/NOxregions.

In Eq.(2),a variable dX is defined as the deviation between chemical component concentrations in the sensitivity experiments and the control test:

Here,Xiis the concentration of chemical components in the sensitivity experiments(Case 1–Case 5-2),while Xcontrolis the concentration in the control experiment.

4.2.Relative importance

To estimate the relative importance of photochemical reactions and different kinds of heterogeneous reactions to the impact of particles on O3,with the initial setup based on observational data,the sensitivity experiments including Control,Case 1,and Case 2-1,2-2,2-3 and 2-4,are analyzed.Case 1 represents photochemical reactions,and Case 2-1,2-2,2-3 and 2-4 were designed to represent the effect of different heterogeneous reactions.The time series of NOx,HOx,O3and HNO3deviations from the control experiment in different case studies are shown in Fig.8,where dNOx,dHOx,dO3and dHNO3means the deviation of NOx,HOx,O3and HNO3,respectively.In Case 1,the average NOxincreases by 19.2 ppb,especially after 12–18 h when the deviation reaches 27.2 ppb.This positive deviation indicates a decline in NOxconsumption.One possible explanation is that the enhancing extinction effects of PM weaken the UV radiation,leading to the reduced photolysis rate.The suppression of NOxconsumption affects the photochemical reactions of O3,resulting in a 101-ppb reduction in the average O3concentration.Due to the weakening of atmospheric oxidation,the HOxconcentration in Case 1 is 30.2 ppt less than in the Control test,and the maximum difference(39.2 ppt)appears after 13 model hours.

Different heterogeneous reactions are taken into account in Case 2-1,2-2,2-3 and 2-4.The heterogeneous reactions of HxOy(OH+HO2+H2O2)are added into Case 2-1,and results from the simulation show that the concentration of HOxradicals in Case 2-1 is 31 ppt less than in the Control test,while NOx,O3and HNO3are slightly different in Case 2-1.This indicates that the heterogeneous reaction of HxOyis not the main contributor to NOxphotolysis,leading to negligible changes in HNO3and O3generation.The heterogeneous absorption of O3on BC surfaces is considered in Case 2-2.Due to the uptake of O3,VOC oxidation processes are suppressed and a slight decrease in HOxcan be seen.Case 2-3 simply includes the heterogeneous reaction of NO2,and a signi ficant change can be observed for all four components.NOxconsumption is accelerated compared with Case 1,and thus the average concentration of precursor-NOxis reduced by 2.8 ppb,with the largest decline of 22.2 ppb after 4 h.NO2transforms into NO?3through a heterogeneous reaction,and NO2no longer participates in the generation of O3after the transformation.This leads to the phenomenon that O3and HOxconcentrations reduce by 36.2 ppb and 2.8 ppt,respectively,after 12 h.However,HOxand O3concentrations rise dramatically at the beginning,which can be explained by the generation of nitrous acid(HONO),which is one of the main products in NOxheterogeneous reactions.HONO is a significant source of OH in the troposphere(Aumont et al.,2003);therefore,the generation of HONO leads to more HOxand,subsequently,O3.Case 2-4 simulates the existence of NOz(e.g.,NO3,N2O5,HNO3)heterogeneous reactions.Only the average HNO3concentration is seen to decrease among the four components.Without considering the effect of aerosols on photolysis rate,NO?3,which is generated by the heterogeneous reactions of HNO3,has little influence on the photolysis reactions of NOx,HOxor O3;thus,the heterogeneous reactions of HNO3have negligible impact on NOx,HOxand O3concentrations.

Fig.8.Variation of(a)dNOx,(b)dHOx,(c)dO3and(d)dHNO3in Case 1 to Case 2-4(crosses represent Case 1;dots represent Case 2-1;triangles represent Case 2-2;squares represent Case 2-3;and diamonds represent Case 2-4).

In this case study,the initial concentration of each species is set based on observations in Nanjing.In Fig.8c,dO3in all of the cases is shown as being negative by the end of the reaction,which demonstrates that O3in Nanjing can be reduced by PM through the impact of particles on both photochemical and heterogeneous reactions.Also,this negative correlation is consistent with the observational result.Generally,Case 1 shows a greater reduction in dO3than Case 2-1,Case 2-2,Case 2-3 and Case 2-4 combined.This indicates that the effect on photochemical reactions has a greater impact on O3than that on heterogeneous reactions.The rate of reduction in Case 1 slows at the end of the 48-h reaction in Fig.8c.This is possibly because the emissions replenishment mechanism is not taken into consideration in our zero-dimensional box model;therefore,remaining reactants are not as sufficient as at the beginning.Apart from Case 2-1,it is worth noting that the result from Case 2-2 shows the largest reduction in O3.This demonstrates that the process of O3absorption on particle surfaces plays the major role in O3depletion among all the heterogeneous reactions.

4.3. Scattering capability

In order to understand the impact of atmospheric scattering capability on the photolysis reaction rate,Case 3-1(clean atmosphere)and Case 3-2(AOD=1 or 2 and SSA=1,0.9 or 0.8)were simulated.The photolysis rate(J)was calculated by the UV module of the box model,and the profiles of J[NO2](the photolysis rate of NO2)and J[O3](the photolysis rate of O3)in different AOD and SSA scenarios are shown in Fig.9.According to previous research,surface particles were considered evenly distributed under 2 km in our experiments(Boynard et al.,2014;Pet?j? et al.,2016).When calculating the UV radiation,four schemes were set to represent different atmospheric conditions,among which “clean atmosphere”represented the atmosphere without PM.The J[NO2]and J[O3]affected by the extinction effect of particles are smaller than those in clean atmospheric conditions near the ground.However,photolysis rates change with height because of different extinction effects of PM.The UV radiation that arrives at the surface will be scattered back into the air by scattering aerosols,enhancing the UV radiation above the aerosol layer.Thus,photolysis rate coefficients in high-scattering scheme(SSA=1)elevate correspondingly,and become even larger than in the clean atmosphere scheme.Tie et al.(2005)proposed that absorbent PM(e.g.,BC)can reduce the UV radiation both below and above aerosols layers,and this is consistent with our results showing that the diminution of the photolysis rate coefficient occurs at almost all altitudes when SSA=0.8.In Fig.9,the effects of scattering and absorbing aerosols are more evident in the case of high AOD.The calculated SSA of all pollutants in Nanjing is relatively low at 0.76;thus,the photolysis rate is reduced in association with the existence of PM.

Fig.9.Photolysis rate profile of(a)NO2and(b)O3in different AOD and SSA scenarios(Case 3-1 and Case 3-2)in the box-model simulation.

4.4.PM2.5concentration

Sensitivity experiments were carried out to evaluate the relative importance of photochemical and heterogeneous reactions under different PM concentrations.Photochemical(Case 4-1)and heterogeneous reactions(Case 4-2)were compared when the PM2.5concentration changes from 0 to 175μg m?3,and the variations of NOx,HOx,O3and HNO3in each experiment are shown in Fig.10.Heterogeneous absorption is strengthened when the PM concentration increases due to the rising surface area of particles,which provides a larger reaction surface.Therefore,dNOx,dHOx,dO3and dHNO3in heterogeneous reactions decline with increasing particle concentrations.

The NOxconcentration changes differently under different interaction mechanisms.In Fig.10a,the photochemical reactions enlarge the NOxconcentration by 0.17 ppbμg?1m3on average,while heterogeneous reactions on BC and soluble particle surfaces continue to consume NOxat a rate of?0.01 ppb μg?1m3.The photolysis rate decreases because of the extinction effect of PM,and thus the consumption of NOxslows.The change in HOxradicals shows a similar trend in two different interactions in Fig.10b.When the concentration of PM2.5is less than 90 μg m?3,the rate of decline of dHOx,with consideration of the impact of PM on photochemical reactions,is?0.31 pptμg?1m3.The increasing particle concentrations suppress photolysis and affect the atmospheric oxidation,resulting in the declining concentration of HOxradicals.However,heterogeneous reactions play a more important role in the decline of dHOxthan photolysis under conditions of lower PM2.5levels,with the average rate of change being ?0.43 pptμg?1m3when the PM2.5concentration is less than 35 μg m?3.With a rise in the PM concentration,the effects of the two interaction mechanisms tend to be the same.Heterogeneous and photochemical effects on NOxand HOxradicals finally result in the reduction of dO3.In Fig.10c,the average rates of change of dO3affected by both heterogeneous and photochemical reactions are?0.44 ppbμg?1m3and?0.71 ppbμg?1m3,respectively,when PM2.5is under 175μg m?3.At lower particle concentrations,heterogeneous reactions are the main influence on O3variation;however,when PM2.5＞50μg m?3,photolysis overtakes heterogeneous reactions and plays the major part in O3reduction.In the case study of Nanjing,photochemical reactions are more dominant than heterogeneous reactions for PM2.5=105μg m?3.In Fig.10d,dHNO3rises to a peak of 1 ppb when the PM2.5concentration is 35 μg m?3,after which the change in HNO3concentration is similar to that of HOxradicals.The rates of dHNO3are ?0.03 ppb μg?1m3for photochemistry and ?0.001 ppb μg?1m3for heterogeneous reactions when PM2.5changes from 35μg m?3to 175μg m?3.

4.5.VOCs/NOxratio

Fig.10.Variation of(a)dNOx,(b)dHOx,(c)dO3and(d)dHNO3in two interaction processes under different PM2.5 concentrations.The solid lines represent the photolysis impact process(Case 4-1)and the dashed lines the heterogeneous uptake process(Case 4-2).

Fig.11.With varying VOCs/NOxvalues(VOCs=0–200 ppb and NOx=0–100 ppb),(a,c)show the changes in dO3and dOH when considering photochemical reactions in Case 5-1,and(b,d)show the changes when considering heterogeneous reactions in Case 5-2.

We used a set of sensitivity experiments to verify how O3and OH are affected by photochemical(Case 5-1)and heterogeneous reactions(Case 5-2)in the O3–NOx–VOCs system,with VOCs varying from 0 to 200 ppb and NOxfrom 0 to 100 ppb.Figures 11a and c show dO3and dOH as always negative,meaning PM can reduce O3and OH through photochemical reactions.This negative effect of particles on O3has also been found in other studies.For instance,Li et al.(2011a)reported that surface O3was reduced by 5%and OH concentrations by 51%via changing photolysis frequencies in central eastern China.Lou et al.(2014)also found a reduction in O3when considering the impact of aerosols on O3by altering photolysis rates in eastern China.The impact of particles on heterogeneous reactions is illustrated in Figs.11b and d.Generally,O3will be reduced by heterogeneous reactions;however,when NOxis more than 45 ppb and VOCs less than 60 ppb(VOCs/NOxratio of less than 1.33),which represents a typical VOC-sensitive regime,O3will increase.This is possibly because excess NOxwill consume O3in a VOC-sensitive region(Geng et al.,2009;An et al.,2015).However,when NOxis consumed by heterogeneous reactions on PM surfaces,the reduction of NOxlimits the consumption of O3,leading to a rise in the O3concentration.According to earlier discussion,Nanjing is a VOC-sensitive region and the initial concentrations of NOxand VOCs were 40 ppb and 70 ppb,respectively.The VOCs/NOxratio of Nanjing is shown in Fig.11 as a black circle.Although Nanjing is a VOC-sensitive region,the VOCs/NOxratio is not low enough to be affected by the increasing effect of heterogeneous reactions in such a region.

Generally,PM reduces O3through altering the photolysis rate and heterogeneous reactions.However,O3can also be increased by PM in two situations:with a low concentration of particles(PM2.5＜50μg m?3)and low VOCs/NOxratio(＜ 1.33);and with high PM pollution(PM2.5＞ 50 μg m?3)and high atmospheric scattering capability.In the case study of Nanjing,the average concentration of PM2.5is larger than 50μg m?3,according to the observational analysis.Furthermore,the SSA is 0.76 and VOCs/NOxratio is 1.75.Thus,particles lead to O3reduction through both photochemical and heterogeneous reactions.

5.Conclusions

The interaction between PM and O3in Nanjing,in the western part of the Yangtze River Delta region in China,was investigated in this paper through observational analysis and model simulation results.

One year of observations at two sites in 2008 shows that the average concentrations of PM2.5and O3are 64.6±47.4μg m?3and 24.6±22.8 ppb,respectively,for CCM station,and 94.1±63.4μg m?3and 16.9±14.9 ppb for PK station.A significant diurnal pattern was observed for O3,with a maximum value in the afternoon and minimum at night;whereas,PM2.5shows a double-peak pattern,with one between 0700 and 0900 LST and the other between 1900 and 2100 LST.Negative correlation between PM and O3can be found in the seasonal,weekly and daily variations.In a case study at PK station from 24–27 March 2008,a 6-ppb decrease in O3during daytime was found,with an increase in PM2.5by 76.8μg m?3,along with 0.5(500 nm)and 0.8(380 nm)increases in AOD.The correlation coefficient between PM2.5and O3reached?0.46.

Previous studies show that the O3level is not only affected by the variability of emissions and meteorological conditions,but also chemical processes such as photochemical reactions and heterogeneous reactions(Nishanth et al.,2014;Lou et al.,2014;Akimoto,2016).The box model simulation in the case study suggests that PM can diminish O3through affecting both kinds of reactions,and the effect of photochemical reactions has greater influence than heterogeneous reactions.Of all the heterogeneous reactions considered in this study,the absorption of O3on particle surfaces plays the most important role in O3reduction.

Modeling results also indicate that the scattering capability of particles has an impact on the rate of photolysis.When SSA=1,the photolysis rates of J[NO2]and J[O3]reach a maximum.On the contrary,when SSA=0.8,the scattering capability grows weaker and the absorbing capability grows stronger.Consequently,J[NO2]and J[O3]decline.

The contribution of aerosol-photolysis and heterogeneous reactions varies depending on the concentration of PM.Heterogeneous reactions play a major role in O3destruction at relatively lower particle concentrations(PM2.5＜50μg m?3),whereas aerosol-photolysis can cause a greater O3reduction when the PM concentration is higher(PM2.5＞50μg m?3).When PM2.5reaches 175μg m?3,with the impact of aerosols on photochemical reactions,the rate of variation between O3and PM2.5concentrations,i.e.,dO3/dPM2.5,is?0.71 ppbμg?1m3.For heterogeneous reactions,the rate is?0.44 ppbμg?1m3.

Besides the scattering capability and concentration of PM,this study shows that the ratio of VOCs/NOxcan also affect the interactions between PM and O3.Although particles generally reduce O3by their impact on photolysis and heterogeneous reactions,O3can,however,be increased by particles through heterogeneous reactions in typical VOC-sensitive regions with low VOCs(＜ 60 ppb)and high NOx(＞ 45 ppb).

In the case study of Nanjing,the SSA is relatively low at 0.76.The annual average concentration of PM2.5is larger than 50μg m?3and the ratio of VOCs/NOxis 1.75.Therefore,the effects of both photochemical and heterogeneous reactions lead to O3reduction in Nanjing,and the aerosol photolysis effect plays a major role.

Lastly,we want to emphasize the complexity of the interaction between O3and particles.On the one hand,increasing PM enhances the atmospheric extinction capacity and weakens the photochemical formation of O3.On the other hand,rising particle numbers provide a larger surface area for heterogeneous reactions,thus influencing the O3concentration.Our results are acquired based on a short period of observation in Nanjing and zero-dimensional modeling.Clearly,more studies need to be carried out to verify the robustness of our results.Nonetheless,the mechanism revealed in this study,together with other recent research on the interactions between O3and PM,should be considered in the chemical module of many climate models.In future work,more simulations based on chemistry models could provide more information to qualify the impact of particles on O3generation.In-situ measurements of chemical species and particles are also essential,to evaluate model performance with respect to the interaction mechanisms involved.

Acknowledgements.This work was supported by the National Key Basic Research&Development Program of China(Grant Nos.2016YFC0203303,2016YFC0208504 and 2014CB441203),the National Natural Science Foundation of China(Grant Nos.91544230,41575145 and 41621005),and the EU 7th Framework Marie Curie Actions IRSES project:REQUA(Grant No.PIRSESGA-2013-612671).