Nitrate formation from heterogeneous uptake of dinitrogen pentoxide during a severe winter haze in southern China

Nitrate (NO−3 ) has become a major component of fine particulate matter (PM2.5) during hazy days in China. However, the role of the heterogeneous reactions of dinitrogen pentoxide (N2O5) in nitrate formation is not well constrained. In January 2017, a severe haze event occurred in the Pearl River Delta (PRD) of southern China during which high levels of PM2.5 (∼ 400 μg m−3) and O3 (∼ 160 ppbv) were observed at a semi-rural site (Heshan) in the western PRD. Nitrate concentrations reached 108 μg m−3 (1 h time resolution), and the contribution of nitrate to PM2.5 was nearly 40 %. Concurrent increases in NO−3 and ClNO2 (with a maximum value of 8.3 ppbv at a 1 min time resolution) were observed in the first several hours after sunset, indicating an intense N2O5 heterogeneous uptake by aerosols. The formation potential of NO−3 via N2O5 heterogeneous reactions was estimated to be between 29.0 and 77.3 μg m−3 in the early hours (2 to 6 h) after sunset based on the measurement data, which could completely explain the measured increase in the NO−3 concentration during the same time period. Daytime production of nitric acid from the gasphase reaction of OH+NO2 was calculated with a chemical box model built using the Master Chemical Mechanism (MCM v3.3.1) and constrained by the measurement data. The integrated nocturnal nitrate formed via N2O5 chemistry was comparable to or even higher than the nitric acid formed during the day. This study confirms that N2O5 heterogeneous chemistry was a significant source of aerosol nitrate during hazy days in southern China.


Introduction
Severe haze in China has been a major concern of the regulatory and scientific communities in recent years.Nitrate has also been identified as an important component of PM 2.5 during hazy days in both summer and winter (e.g., Huang et al., 2014;Li et al., 2018;Pathak et al., 2009;Zhang et al., 2015).Moreover, the proportion of nitrate in PM 2.5 has steadily increased over the last decade due to the lagged control of NO x emissions compared to SO 2 (Fu et al., 2014;Geng et al., 2017;Qu et al., 2017;Reuter et al., 2014;X. Wang et al., 2016).As a result, the concentrations of nitrate in PM 2.5 / PM 1.0 have been even higher than those of sulfate during some haze events (Ge et al., 2017;Li et al., 2017;Liu et al., 2015;Yang et al., 2017;Yue et al., 2015).
Nitrate is formed from NO x during both the daytime and nighttime.During the day, nitric acid (HNO 3 ) is produced through the gas-phase reaction between OH and NO 2 (Reaction R1), while this pathway is insignificant at night due to very low OH concentrations (e.g., Seinfeld and Pandis, 2016).Nitric acid can react with ammonia (NH 3 ) to form ammonium nitrate (NH 4 NO 3 ), and an equilibrium can be reached for these three compounds between the gas phase and the particle phase (Reactions R2, R3).During nighttime, the heterogeneous uptake of N 2 O 5 , which is formed from reactions involving O 3 , NO 2 , and NO 3 , becomes a source of nitrate and also produces gaseous ClNO 2 when chloridecontaining aerosol is present (Reactions R4-R7; Finlayson-Pitts et al., 1989).This nitrate formation pathway is only important at night due to the fast photolysis of NO 3 during the day.Compared to the relatively well-understood formation of aerosol nitrate via the OH+NO 2 reaction, the contribution from N 2 O 5 heterogeneous reactions has been poorly quantified due to the limited knowledge of key factors controlling the heterogeneous processes, such as the N 2 O 5 uptake coefficient (γ N 2 O 5 ) and ClNO 2 yield (φ ClNO 2 ) (Brown and Stutz, 2012;Chang et al., 2011).
In addition to the modeling approach, field measurements of trace gases and aerosol composition have been used to infer the contribution of N 2 O 5 heterogeneous chemistry to nitrate formation.Pathak et al. (2009) postulated the importance of N 2 O 5 heterogeneous reactions with respect to the high aerosol nitrate observed in Beijing and Shanghai in summertime by examining the variation of nitrate with the change in relative humidity (RH) and the equilibrium between anions and cations in PM 2.5 .Pathak et al. (2011) further investigated nitrate formation using a coupled aqueous phase radical mechanism (CAPRAM) and a gas-phase chemistry mechanism (RACM, without ClNO 2 chemistry).By constraining the uptake coefficient of N 2 O 5 in the range from 0.001 to 0.1, they reproduced the observed enhancement of nitrate and suggested that N 2 O 5 uptake by aerosols contributed to 50 % to 100 % of the nighttime increase in nitrate.A similar method was recently used by Wen et al. (2018) to simulate the summertime nitrate formation in the North China Plain (NCP), which demonstrated the dominant contribution of N 2 O 5 heterogeneous reactions to night-time nitrate formation.Based on the observed covariation of nitrate and RH, Wang et al. (2009) speculated that N 2 O 5 reactions dominated the nitrate formation on polluted days with high NO 2 and O 3 concentrations in Shanghai.Neither N 2 O 5 nor ClNO 2 was measured during these early observationbased studies.A recent study (H.Wang et al., 2017) inferred γ N 2 O 5 from the measured N 2 O 5 on 4 days in urban Beijing and estimated the lower limit of the formation potential of aerosol nitrate assuming a unity φ ClNO2 because ClNO 2 was not measured.Their result showed a comparable contribution to nitrate formation from the N 2 O 5 heterogeneous chemistry and from the daytime pathway of the OH + NO 2 reaction.
In the present study, N 2 O 5 , ClNO 2 , and the related chemical and meteorological parameters were measured at a semirural site in the Pearl River Delta of southern China from 2 to 15 January 2017.A severe haze event was observed during the field study with PM 2.5 reaching 400 µg m −3 and O 3 peaking at 160 ppbv.ClNO 2 , which is only known to be produced from N 2 O 5 heterogeneous uptake, reached 8.3 ppbv, which is the largest reported value to date and revealed extremely active N 2 O 5 chemistry during the abovementioned episode.The concurrent measurements of N 2 O 5 , ClNO 2 , and aerosol nitrate provide better constraints for elucidating nighttime NO 3 / N 2 O 5 chemistry and aerosol nitrate formation.An overview of the measurement data was first presented.The nighttime processes that led to the formation of nitrate (e.g., the production of NO 3 and N 2 O 5 , N 2 O 5 uptake coefficient, and ClNO 2 yield) were analyzed.The nighttime formation potential of nitrate was estimated based on these data and compared to the measured increase in nitrate.The daytime production of nitric acid via the OH + NO 2 reaction was calculated based on a box model using the Master Chemical Mechanism (MCM v3.3.1) and compared to the nighttime formation potential of nitrate.

Site description
The field observation was conducted at the Guangdong atmospheric supersite, a semi-rural site located at Hua Guo Shan (HGS, 22.728 • N, 112.929 • E) southwest of the city of Heshan from 2 to 15 January 2017.As shown in Fig. 1, HGS is a hill with a height of 60 m a.s.l.All measurement instruments were located on the fourth floor of a four-story building on the top of the hill.The observation site was located in the western PRD where the economic activity and population density are much lower than the central PRD.There are five main roads near the HGS site, including three national roads (G325, G94, and G15), and two provincial roads (S272 and S270).The hill is covered by subtropical trees and surrounded by similar hills within close proximity; a few residents live at the foot of the hill, and there is some farmland in the area.

Chemical ionization mass spectrometer
N 2 O 5 and ClNO 2 were simultaneously observed using a quadrupole chemical ionization mass spectrometer (CIMS, THS Instruments, Atlanta).The same instrument had been used in several previous studies in southern and northern China (Tham et al., 2016;T. Wang et al., 2016;Z. Wang et al., 2017).The reader can refer to these earlier papers for a detailed description of the measurement principle, calibration, and maintenance procedures.Briefly, ambient N 2 O 5 and ClNO 2 are converted to ion clusters of I(N 2 O 5 ) − and I(ClNO 2 ) − by Iodide ions (I − ) that are produced by exposing CH 3 I / N 2 (0.3 % v/v) to an alpha radioactive source, 210-Po, and are subsequently detected at 235 and 208 m z −1 , respectively.Activated carbon packed in a filter was used to determine the instrument background, which were 10.2±2.2 and 8.9 ± 2.0 Hz on average for N 2 O 5 and ClNO 2 , respectively.In situ offline calibration was carried out every day for N 2 O 5 and every 2 days for ClNO 2 by mixing the respective synthetic standard into humidified zero air (with RH controlled at 60 % in the present study).The N 2 O 5 standard was generated by reacting excess NO 2 with O 3 and determined from the decrease of NO 2 , and the ClNO 2 was synthesized by the uptake of a known concentration of N 2 O 5 on a NaCl slurry (see T. Wang et al., 2016 andTham et al., 2016 for details).The average sensitivity of N 2 O 5 and ClNO 2 was 0.9 ± 0.3 and 0.7 ± 0.2 Hz pptv −1 , respectively.The dependence of the sensitivity on the relative humidity was measured during the field study (see Fig. S1) and was used to correct for the RH effect based on the measured ambient RH values.The detection limits of N 2 O 5 and ClNO 2 were both 6 pptv (2σ , 1 min-averaged data).
The inlet of the CIMS instrument was set approximately 1.5 m above the roof with 6 m long PFA-Teflon tubing as the sampling line.The total sampling flow was set as 11 stan-dard liters per minute (SLPM).Four SLPM were diverted into the CIMS, the O 3 analyzer, and NO x analyzer, and the remaining part was evacuated directly from the system.The total residence time was less than 1 s in the sampling system.Following our previous practice, the inlet tubing and fittings were replaced every afternoon and washed in an ultrasonic bath to reduce the influence of particle deposition on the tubing wall.The loss of N 2 O 5 on the tubing wall was checked on site by injecting N 2 O 5 into the ambient air before and after the tubing replacement; the loss was found to be around 10 % in the"clean" tubing and nearly 40 % by the next afternoon.Because our analysis mainly focused on data from the first few hours of evening, the loss was insignificant and thus was not corrected in our final data.However, this bias could be important in the later period before tube replacement.The uncertainty of the measurement was estimated to be ±25 % for both N 2 O 5 and ClNO 2 (T.Wang et al., 2016).The time resolution for the measurement was approximately 10 s, and the derived data were later averaged to a time resolution of 1 min for further analysis.

Other measurements
Trace gases of CO, SO 2 , O 3 , NO x , total reactive nitrogen (NO y ), nitrous acid (HONO), C2 to C10 non-methane hydrocarbons (NMHCs), oxygenated hydrocarbons (OVOCs), and aerosol chemical composition and number concentrations were also measured.Table 1 summarizes the principle, detection limit and uncertainty of the measuring instruments.
CO was observed using a gas filter correlation analyzer (Thermo Scientific, Model 48i).SO 2 was measured using a pulsed fluorescence analyzer (Thermo Scientific, Model 43i).O 3 was determined using a UV photometric analyzer (Thermo Scientific, Model 49i).NO and NO 2 were detected with a chemiluminescence instrument (Thermo Sci- entific, Model 42i) with a photolytic converter to convert NO 2 to NO (Xu et al., 2013).NO y was determined using a chemiluminescence analyzer which was equipped with a molybdenum oxide (MoO) catalytic converter (Thermo Scientific, Model 42i-Y).HONO was detected using a long path absorption photometer (QUMA, Model LOPAP-03) (Xu et al., 2015).NMHCs were determined using an online gas chromatograph (GC) coupled with a flame ionization detector (FID) and a mass spectrometer (MS) (Wang et al., 2014).NMHCs were only measured from 2 to 8 January 2017 due to the maintenance of the GCMS after 8 January.OVOCs (e.g., formaldehyde, acetaldehyde, acetone, and methyl ethyl ketone) were sampled using 2,4dinitrophenylhydrazine (DNPH) cartridges every 3 h and were later analyzed with a high-performance liquid chromatography (HPLC) system (Cui et al., 2016).Concentrations of PM 2.5 were detected with a multi-angle absorption photometer (MAAP; Thermo Scientific, Model 5012).The ionic compositions of PM 2.5 were measured with an ion chromatography (GAC-IC) system equipped with a gas and aerosol collector at a time resolution of 30 min (Dong et al., 2012), and the data were also averaged every 1 h to meet the time resolution of other components of PM 2.5 .Organic carbon (OC) and elemental carbon (EC) were measured using an online OC / EC analyzer (RT-4, SUNSET) with a time resolution of 1 h (Bauer et al., 2009).A scanning mobility particle sizer (SMPS Model 3936L75, TSI) was used to determine the dry-state particle number size distribution, covering the size range from 16.5 to 1000 nm.The ambient (wet) particle number size distributions were estimated based on a size-resolved kappa-Köhler function considering the variation with relative humidity (Hennig et al., 2005;Liu et al., 2014).In the present study, data with a RH greater than 90 % were excluded due to the large uncertainty of the growth factor at very high RH.The aerosol surface area density was then derived using the ambient particle number size distribution (wet) and an assumption of spherical particles with an estimated uncertainty of around 30 % (Tham et al., 2016;Z. Wang et al., 2017).
Meteorological parameters were measured with a portable weather station (WXT520, Vaisala, Finland), including temperature, relative humidity (RH), wind direction, wind speed, and pressure.A pyranometer (CMP22, Kipp & Zonen B.V., Holland) was used to measure the solar radiation, and the data were then utilized to derive the photolysis frequency of NO 2 based on the method of Trebs et al. (2009).

Chemical box model
To estimate the daytime formation of nitric acid via the reaction of OH + NO 2 , an observation-based chemical box model developed with the latest version of the MCM (v3.3.1)(Jenkin et al., 2003(Jenkin et al., , 2015;;Saunders et al., 2003) and an updated chlorine (Cl) radical chemistry module (Xue et al., 2015) was utilized to estimate the mixing ratio of OH radicals and the reaction rate of OH + NO 2 .The integrated production of nitric acid during daytime was then calculated based on the simulation results.The box model was constrained with the observation data every 10 min, including the data relating to N 2 O 5 , ClNO 2 , HONO, O 3 , NO, NO 2 , SO 2 , CO, C2 to C10 NMHCs, OVOCs (formaldehyde, acetaldehyde, acetone, and MEK), temperature, aerosol surface area density, and J (NO 2 ), which were first averaged or interpolated.Average concentrations of NMHC species during the day (07:00 to 17:00) and at night (17:00 to 07:00 of the next day) are shown in Table S1.A function considering the variation of the solar zenith angle (Saunders et al., 2003) was used to calculate the photolysis frequencies of HONO, O 3 , and other species in clear sky, which were then corrected with the J (NO 2 ) values in the real environment.The J (ClNO 2 ) was treated the same as in Tham et al. (2016).The lifetime of unconstrained species with respect to the physical loss was set as 8 h in a boundary layer with a depth of 1000 m (equivalent to 3.47 × 10 −5 s −1 ) in order to avoid their accumulation.The model was run from 00:00 on 3 January to 11:00 on 8 January 2017.To stabilize the intermediate species, the simulation for the first 24 h was repeated six times.Sensitivity tests were carried out by reducing the input concentrations by 10 % to check the deviation of the average daytime (07:00-17:00) rate of the OH + NO 2 reaction.The simulated rate of the OH + NO 2 reaction was most sensitive to HONO, followed by NO x and OVOCs (see Sect.S1 and Fig. S2 in the Supplement).
3 Results and discussion

Overview of the observation
Figure 2 shows the time series of N 2 O 5 , ClNO 2 , components of PM 2.5 , and related trace gases and meteorological parameters from 18:40 on 2 January to 11:00 on 15 January 2017.The average temperature and RH values during the measurement period were 17 ± 4 • C and 86 ± 14 %, respectively.A severe pollution episode occurred from 3 to 7 January due to stagnant meteorological conditions (Fig. 3a), and the concentrations of most pollutants decreased to very low levels on 9 January and from 12 to 15 January, which corresponded to the change in weather conditions.The most polluted days were 5 and 6 January with the highest PM 2.5 value of 400 µg m −3 and the highest O 3 concentration of 160 ppbv.The PM 2.5 data from the PRD regional air quality monitoring network revealed that the HGS site was within the most polluted area during this haze event (Fig. 3b).This pollution event was characterized by concurrent high levels of PM 2.5 and O 3 and was in contrast to the winter haze in north China, which experienced high PM 2.5 values but low O 3 concentrations (e.g., Sun et al., 2016;H. Wang et al., 2018a).The mixing ratios of N 2 O 5 and ClNO 2 were up to 3358 and 8324 pptv (1 min time resolution), respectively, indicating active N 2 O 5 heterogeneous chemistry.
Very high concentrations of aerosol nitrate (up to 108 µg m −3 ; 1 h time resolution) were also observed during the multi-day episode.Nitrate contributed to 24 % of the total PM 2.5 mass concentration on average, which was comparable to that of organic matter (OM = 1.7•OC, 28 %) and much higher than that of sulfate (16 %) and ammonium (11 %).The contribution of nitrate to the PM 2.5 increased with an increase in the nitrate concentration, and reached nearly 40 % at its peak, indicating that nitrate was a dominant component of the PM 2.5 on the most polluted days.The concentration of NO − 3 exhibited a concurrent increase with that of ClNO 2 in the early nighttime from 3 to 4, 4 to 5, 5 to 6, and 9 to 10 January (see Fig. 4), suggesting that N 2 O 5 heterogeneous reactions significantly contributed to the formation of nitrate during the night.The measured increases of the NO − 3 concentration during these four nights were 17.1, 50.9, 43.3, and 32.7 µg m −3 , respectively.A similar increase in ClNO 2 was observed from 6 to 7 January, but the composition of the PM 2.5 was not available due to instrument maintenance.Apart from chemical reactions, the evolution of the planetary boundary layer (PBL) also affects the concentrations of trace gas and aerosols.The height of PBL generally decreases after sunset due to the fast decrease in the land temperature, which can lead to the accumulation of primary pollutants (and secondary pollutants) at the surface if significant local sources are present.For example, on the night from 4 to 5 January (see Fig. 5), the CO and NO y levels increased between 18:00 and 19:00 LT with the enhancement of ClNO 2 and nitrate, indicative of an accumulation of primary emissions; however, after this period of time the primary pollutants decreased for 3 h while ClNO 2 and nitrate continued to increase due to the nighttime chemical process.
In the remainder of this paper, we will focus on the detailed analysis of the abovementioned five nights to investigate the role of N 2 O 5 heterogeneous chemistry in nitrate formation.

Production of NO 3 and N 2 O 5
The first step in the nighttime nitrate formation via N 2 O 5 chemistry is the production of NO 3 and N 2 O 5 .To gain insight into the key factors affecting the NO 3 / N 2 O 5 chemistry, the variation of N 2 O 5 and the production rate of NO 3 were examined in relation to some relevant gases and meteorological parameters on the selected five nights.Figure 5 shows the data for the night from 4 to 5 January as an example.Some common features were identified for all five nights.In general, low wind speed (<2.0 m s −1 ) at night facilitated the accumulation of air pollutants, and high RH was favorable for N 2 O 5 heterogeneous uptake.In addition, high aerosol surface area density provided interfaces for N 2 O 5 heterogeneous reactions.
In the first couple of hours after sunset (Fig. 5, red rectangle), N 2 O 5 exhibited a peak and then quickly dropped to hundreds of pptv, while nitrate and ClNO 2 concurrently increased, which was indicative of the local production and loss of N 2 O 5 .NO was below the detection limit during this period.The production rates of NO 3 (P NO 3 = k NO 2 +O 3 [NO 2 ] [O 3 ]) were the fastest just after sunset and decreased gradually due to reduced O 3 levels.There was a period later in the night (22:00 to 01:00) when fresh emissions of NO were observed, and the production of NO 3 was suppressed due to the titration of O 3 by NO.In the later nighttime, NO was below the detection limit (Fig. 5, blue rectangle).During this period, NO 3 and N 2 O 5 were produced at moderate rates, and the very low N 2 O 5 concentrations (below the detection limit) suggested a fast loss of N 2 O 5 prob-

The N 2 O 5 uptake coefficient and ClNO 2 yield
The N 2 O 5 uptake coefficient and ClNO 2 yield, in combination with the reactivity of NO 3 with NO and VOCs, determines the loss pathways of NO 3 and N 2 O 5 .To derive the uptake coefficient of N 2 O 5 , a method suggested by McLaren et al. (2010) was applied which treated NO 3 and N 2 O 5 as a whole ([NO 3 ] + [N 2 O 5 ]) without assuming that the chemical system was in a steady state.This approach considers that the change of NO 3 and N 2 O 5 concentrations is mainly due to NO 3 / N 2 O 5 chemistry; thus, it requires that the air mass has relatively stable chemical conditions and is not subject to fresh NO emissions.It also requires that ClNO 2 is produced from the N 2 O 5 chemistry and has an increasing trend to derive the yield of ClNO 2 .This method is applicable for the early nighttime (Fig. 5, red rectangle; Sect.3.2.1)for these five nights.The loss of ] is via the reaction of NO 3 with VOCs and N 2 O 5 heterogeneous reactions, which can both be expressed as pseudo-first-order losses: where k NO 3 and k N 2 O 5 represent the total first-order rate constants for NO 3 and N 2 O 5 , respectively.The loss rate of N 2 O 5 can then be obtained from Eq. (3): Because NO 3 was not measured, it was calculated by assuming an equilibrium of NO 2 -NO 3 -N 2 O 5 as shown in Eq. ( 4).High levels of NO would break this equilibrium.Thus, the periods with detected NO were excluded.5) by interpolating the data of the NMHCs to a 10 min time resolution.The NO 3 reactivity with VOCs (k NO3 ) in the early nighttime ranged from 0.516 to 1.54 × 10 −3 s −1 (Table 2), which was higher than those derived at Mt. TMS in winter 2013 (0.17 to 1.1 × 10 −3 s −1 ) (Brown et al., 2016), but lower than those in the North China Plain during the summertime (2 to 57 × 10 −3 s −1 ) (Tham et al., 2016;H. Wang et al., 2017H. Wang et al., , 2018b;;Z. Wang et al., 2017).NMHCs were not measured from 9 to 10 January 2017.We used the average k NO 3 in the early nighttime from 3 to 4 January as a replacement because these two periods had similar pollution levels for most pollutants.For the later nighttime (Fig. 5, blue rectangle), low levels of N 2 O 5 and moderate levels of P NO 3 made Eq. ( 3) inapplicable even though NO was not detected.
Finally, the uptake coefficient of N 2 O 5 was derived using Eq. ( 6) for every 10 min and averaged for each selected period.In Eq. ( 6), C N 2 O 5 is the mean molecular speed of N 2 O 5 , and S a is the aerosol surface area density.The yield of ClNO 2 was derived via dividing the integrated production of ClNO 2 ([ClNO 2 ] max ) by the integrated loss of N 2 O 5 since sunset as Eq. ( 7).
The relative importance of the NO 3 reactions with VOCs and N 2 O 5 heterogeneous reactions can be examined by comparing the values of the loss coefficient of the NO 3 reac- •K eq and N 2 O 5 heterogeneous reactions (k N 2 O 5 ) (Tham et al., 2016).Based on the calculations, the values of [NO 2 ]•K eq were 1.40 × 10 −5 to 6.07 × 10 −5 s −1 (see Table 2), while those of k N 2 O 5 were 3.78 × 10 −3 to 9.00 × 10 −3 s −1 , which was 2 orders of magnitude higher than k NO 3 [NO 2 ]•K eq , suggesting that N 2 O 5 heterogeneous reactions were the dominant loss pathway for both NO 3 and N 2 O 5 .
The average γ N 2 O 5 and φ ClNO 2 derived for the early night in the five cases are listed in Table 2.The data show that the uptake coefficient ranged from 0.009 to 0.066, which was comparable with previous values derived at Mt. Tai Mo Shan (TMS) in Hong Kong (0.004 to 0.022) (Brown et al., 2016)  and in the North China Plain (0.006 to 0.102) (Tham et al., 2016(Tham et al., , 2018;;H. Wang et al., 2017H. Wang et al., , 2018b;;X. Wang et al., 2017;Z. Wang et al., 2017;Zhou et al., 2018).It is interesting to see a much higher γ N 2 O 5 value (0.066) on 3 January than those on other four nights (0.009-0.015), which resulted from higher P NO 3 but much lower Sa and relatively low N 2 O 5 concentrations on 3 January.We examined known factors affecting the loss of NO 3 and N 2 O 5 such as the concentrations of NO, NMHCs, and aerosol compositions, but found no obvious difference between 3 January and the other nights.The yield in this study varied from 0.18 to 0.32, which was similar to most studies in China (Tham et al., 2016(Tham et al., , 2018;;Z. Wang et al., 2017;Yun et al., 2018;Zhou et al., 2018).

Nitrate formation potential p(NO −
3 ) through N 2 O 5 chemistry

Nighttime p(NO −
3 ) The formation potential of NO − 3 through N 2 O 5 chemistry is the total amount of NO − 3 accumulated from N 2 O 5 heterogeneous loss.It can be calculated by deducting the integrated production of ClNO 2 from the integrated loss of N 2 O 5 as follows: In the early nighttime, the average loss rate of N 2 O 5 (L N 2 O 5 ) ranged from 1.9 to 4.3 ppbv h −1 (Table 2), which was close to the average P NO 3 due to the dominance of the N 2 O 5 heterogeneous reactions in NO 3 and N 2 O 5 loss.Based on the derived N 2 O 5 loss rate and the maximum ClNO 2 concentration, the formation potential of NO − 3 was derived and ranged from 29.0 to 77.3 µg m −3 as shown in Fig. 6.The measured increase of the NO − 3 concentration in the early nighttime can be completely explained by the integrated production of NO − 3 via the N 2 O 5 heterogeneous reactions during the same period.
In the later nighttime, the method described in Sect.3.2.2 was not valid for calculating the N 2 O 5 heterogeneous loss rate as mentioned above.We attempted to estimate the formation potential of nitrate by assuming that the N 2 O 5 heterogeneous reactions continued to dominate the loss of NO 3 + N 2 O 5 in the later nighttime.The k NO 3 values in the later nighttime were comparable to those in the early nighttime, and the high RH (close to 100 %) in the later nighttime was favorable for N 2 O 5 heterogeneous reactions.We assumed that all NO 3 was quickly consumed by the N 2 O 5 heterogeneous reactions, which means that the loss rate of N 2 O 5 approximated to the production rate of NO 3 (L N 2 O 5 ≈ P NO 3 ).As listed in Table 3, the N 2 O 5 loss rates ranged from 0.82 to 1.26 ppbv h −1 , which were significantly lower than those derived in the early nighttime.The derived N 2 O 5 loss rate here and the yield of ClNO 2 in the early nighttime were used to eswww.atmos-chem-phys.net/18/17515/2018/Atmos.Chem.Phys., 18, 17515-17527, 2018 Figure 7.Comparison between the daytime (07:00 to 17:00 LT; assuming that all gas-phase HNO 3 is partitioned into particle phase) and nighttime (17:00 to 07:00 LT of the next day) NO − 3 formation potential.The early nighttime of each day represents the periods in Table 2, including 3 January 17:40-19:00, 4 January 17:00-22:00, 5 January 17:00-22:00, 6 January 17:00-22:40, and 9 January 19:00-00:20.The later nighttime of each day represents the periods in Table 3, including 3-4 January 21:00-05:00, 5 January 01:30-06:50, 5-6 January 23:40-01:10, 6-7 January 23:00-06:00, and 10 January 01:50-03:30.The intercomparison of the NO − 3 formation potential during the day and night on 9 and 10 January was not conducted due to the lack of NMHC data after 8 January which made the model simulation of OH infeasible during the day on 9 January.timate the formation potential of NO − 3 in the later nighttime.As shown in Fig. 6, the nitrate produced during these later periods ranged from 7.3 to 40.3 µg m −3 , which was lower than values in the early nighttime for four nights, indicating that the nighttime nitrate from N 2 O 5 chemistry was mainly produced in the early nighttime.

Comparison with the daytime production of HNO 3
During the daytime, the formation of NO − 3 is mainly from the gas-particle partitioning of the gas-phase HNO 3 formed via the OH + NO 2 reaction.Hence, the daytime formation potential of HNO 3 (p(HNO 3 )) can be treated as the upper limit for the locally produced daytime aerosol nitrate.To calculate the daytime p(HNO 3 ), a box model based on MCM v3.3.1 was used to derive the mixing ratio of OH and the rates of OH + NO 2 as described in Sect.2.4.This model was previously used in our study at Wangdu in northern China (Tham et al., 2016).The calculated mixing ratios of OH at Wangdu with this model compared well with those observed by the laser-induced fluorescence (LIF) technique (Tan et al., 2017).In the present study, the average daytime OH (07:00 to 17:00) mixing ratios were 1.71 to 3.82 × 10 6 molec cm −3 from 3 to 7 January, as listed in Table 4, with the maximum values reaching 3.24 to 6.71 × 10 6 molec cm −3 .The detailed results for OH can be found in Fig. S3.The average production rates of HNO 3 via the OH + NO 2 reaction were 1.40 to 5.21 ppbv h −1 from 3 to 7 January, and the integrated formation potential of HNO 3 during the daytime was 35.7 to 131.8 µg m −3 , which was comparable to the nighttime p(NO − 3 ) that ranged from 69.3 to 102.9 µg m −3 (Fig. 7).The nighttime production of nitrate via the heterogeneous uptake of N 2 O 5 accounted for 43.8 % to 57.7 % of the total nitrate (NO − 3 + HNO 3 ) produced in a 24 h period at the site.These results underscored the important role of N 2 O 5 heterogeneous chemistry in nitrate formation in this severe winter haze in southern China.

Concluding remarks
With the use of concurrent measurements of nitrate, ClNO 2 , and related pollutants, this study demonstrates the important contribution of N 2 O 5 heterogeneous uptake in nitrate formation.Current chemical transport models have difficulties simulating this nitrate production pathway.Therefore, more research efforts are needed to improve the representations of γ N 2 O 5 and φ ClNO 2 for the better prediction of nitrate in the models.The observation-based approach presented here can be applied to investigate nitrate formation in other areas.
Data availability.The data used in this study are available upon request from Tao Wang (cetwang@polyu.edu.hk) and Dingli Yue (dingliyue@163.com).
Author contributions.TW designed the research.WW conducted the CIMS measurements, and YZ, DY, HY, MX, CY, and PS performed the measurements of the other parameters used in this study.HY, TW, MX, and WW analyzed the data.HY and TW wrote the paper, and all authors contributed to discussion and commented on the paper.

Figure 1 .
Figure 1.(a) Google map images of the Pearl River Delta in Guangdong Province and the measurement site (Hua Guo Shan).(b) The topography and major roads (shown by number) adjacent to the measurement site.

Figure 2 .
Figure 2. Time series of N 2 O 5 , ClNO 2 , components of PM 2.5 , and related trace gases and meteorological parameters from 18:40 on 2 January to 11:00 on 15 January 2017.The inset figure shows the variation of the ratio of nitrate to PM 2.5 with an increasing nitrate concentration.The green rectangles in the figure indicate the five days used for detailed analysis.

Figure 4 .
Figure 4.The covariance of aerosol nitrate and ClNO 2 in the early nighttime (at a 30 min time resolution) for four nights.

Figure 5 .
Figure 5. Variation of N 2 O 5 , ClNO 2 , NO − 3 , and trace gases and meteorological conditions during the night from 4 to 5 January 2017 as an example from the five selected nights.

Table 1 .
Technique, limit of detection, and uncertainty of measuring instruments for trace gases and aerosols.

Table 2 .
Average N 2 O 5 concentration values, N 2 O 5 uptake coefficients, ClNO 2 yields, and other related parameters and maximum ClNO 2 concentration values in the early nighttime for the five selected nights.