A proxy for atmospheric daytime gaseous sulfuric acid concentration in urban Beijing

Gaseous sulfuric acid (H2SO4) is known as one of the key precursors for atmospheric new particle formation (NPF) processes, but its measurement remains challenging. Therefore, a proxy method that is able to derive gaseous sulfuric acid concentrations from parameters that can be measured relatively easily and accurately is highly desirable for the atmospheric chemistry community. Although such methods are available for clean atmospheric environments, a proxy that works well in a polluted atmosphere, such as that found in Chinese megacities, is yet to be developed. In this study, the gaseous sulfuric acid concentration was measured in February–March 2018, in urban Beijing using a nitrate based – long time-of-flight chemical ionization mass spectrometer (LToF-CIMS). A number of atmospheric parameters were recorded concurrently including the ultraviolet radiation B (UVB) intensity, the concentrations of O3, NOx (sum of NO and NO2), SO2, and HONO, and aerosol particle number size distributions. A proxy for atmospheric daytime gaseous sulfuric acid concentration was derived via a statistical analysis method using the UVB intensity, [SO2], the condensation sink (CS), [O3], and [HONO] (or [NOx]) as the predictor variables, where square brackets denote the concentrations of the corresponding species. In this proxy method, we considered the formation of gaseous sulfuric acid from reactions of SO2 and OH radicals during the daytime, and the loss of gaseous sulfuric acid due to its condensation onto the preexisting particles. In addition, we explored the formation of OH radicals from the conventional gas-phase photochemistry using O3 as a proxy and from the photolysis of HONO using HONO (and subsequently NOx) as a proxy. Our results showed that the UVB intensity and [SO2] are dominant factors in the production of gaseous sulfuric acid, and that the simplest proxy could be constructed with the UVB intensity and [SO2] alone. When the OH radical production from both homogenously and heterogeneously formed precursors were considered, the relative errors were reduced by up to 20 %. Published by Copernicus Publications on behalf of the European Geosciences Union. 1972 Y. Lu et al.: A proxy for atmospheric daytime gaseous sulfuric acid concentration

Measurements of gaseous sulfuric acid in the lower troposphere are challenging because its ambient concentration is typically quite low (10 6 -10 7 molecule cm −3 ) (Kerminen et al., 2010;Mikkonen et al., 2011).Reported real-time measurements of gaseous sulfuric acid are currently based on chemical ionization mass spectrometry with NO − 3 and its ligands as reagent ions (nitrate CIMS) because nitrate CIMS with an atmospheric pressure interface (API) has a low detection limit for the atmospheric concentration range of gaseous sulfuric acid (Jokinen et al., 2012); furthermore, a constant fraction of sulfuric acid present in the air sample will be ionized by excessive nitrate ions in CIMS under constant instrumental conditions (Kürten et al., 2012;Zheng et al., 2010), which makes the quantification of gaseous sulfuric acid feasible.Arnold and Fabian (1980) measured the negative ions in the stratosphere and derived the concentration of stratospheric gaseous sulfuric acid from the fractional abundances of a series of stratospheric negative ions as well as the rate constants.Later, real-time measurement of sulfuric acid in the lower troposphere was performed using nitrate CIMS (Eisele and Tanner, 1993), with laboratory calibrations performed via the production of known concentrations of OH radicals that were titrated into gaseous sulfuric acid.Since this work, measurements of sulfuric acid using CIMS have been performed around the world (e.g., Berresheim et al., 2000;Bianchi et al., 2016;Chen et al., 2012;Jokinen et al., 2012;Kuang et al., 2008;Kürten et al., 2014;Kurtén et al., 2011;Petäjä et al., 2009;Weber et al., 1997;Zheng et al., 2011), and CIMS has been proven to be a robust tool for gaseous sulfuric acid detection.However, sulfuric acid measurements are still rather sparse due to the high cost of the CIMS instrument and the extensive demand for special-ized expertise with respect to instrument calibration, maintenance, and data processing, etc.Therefore, a proxy for gaseous sulfuric acid concentration is highly desirable.
Proxies for the estimation of atmospheric gaseous sulfuric acid concentrations were previously developed to approximate measurement results of sulfuric acid in Hyytiälä, southern Finland (Petäjä et al., 2009), assuming that gaseous sulfuric acid is formed from reactions between SO 2 and OH radicals, and lost due to its condensation onto preexisting particles.The derived simplest proxy can be written as Eq.(1) below, where the square brackets denote the concentration, and the authors recognize that the proxies might be site-specific and should be verified prior to utilization in other environments.
[ Mikkonen et al. (2011) later developed a couple of statistical proxies based on measurements of sulfuric acid at six European and North American sites, including urban, rural, and forest areas.Their results showed that the radiation intensity and [SO 2 ] are the most important factors to determine the concentration of sulfuric acid, and that the impact of the condensation sink (CS) for gaseous sulfuric acid, is generally negligible.In several proxies developed by Mikkonen et al. (2011), the correlation between the gaseous sulfuric acid concentration and the CS is positive, which is contrary to what one would expect because a larger CS normally leads to a faster loss for gaseous sulfuric acid.In addition, the performance of a proxy equation is site-specific due to varying atmospheric conditions from one site to another, which implies that the proxy suggested by Mikkonen et al. (2011) might not work well in locations with atmospheric environments different to those at the six sites in that study.
Beijing is a location with typical CS values (e.g., 0.01-0.24s −1 in the 5 %-95 % percentiles in this study) that are 10-100 times higher (Herrmann et al., 2014;Wu et al., 2007;Xiao et al., 2015;Yue et al., 2009;Zhang et al., 2011) and typical SO 2 concentrations that are 1-10 times higher (M.Wang et al., 2011;Wu et al., 2018) than those in Europe and North America (Dunn et al., 2004;Mikkonen et al., 2011), yet measured gaseous sulfuric acid concentrations are relatively similar in these environments (Chen et al., 2012;Smith et al., 2008;Z. B. Wang et al., 2011;Zheng et al., 2011).Whether previous proxies developed for European and North American sites work in Beijing remains to be tested.Furthermore, in addition to the gas-phase reaction between O( 1 D) and water molecules (Crutzen and Zimmermann, 1991;Logan et al., 1981), photolysis of HONO could potentially be an important source of OH radicals in the atmosphere in the early morning (Alicke et al., 2002(Alicke et al., , 2003;;Elshorbany et al., 2009;Li et al., 2012) and during the daytime (Acker et al., 2005;Aumont et al., 2003;Kleffmann, 2007).An experimental study measuring HONO near the surface layer estimated that HONO was a main contributor to OH production in Beijing, with HONO's contribution being larger than 70 % at around 12:00-13:00 UTC/GMT+08:00, except for in summer when the contribution of O 3 dominated (Hendrick et al., 2014).Given the distinct characteristics of these two OH radical formation pathways, they should both be included and evaluated separately when a proxy for atmospheric gaseous sulfuric acid concentration is being built.The reactions between SO 2 and Criegee intermediates formed from the ozonolysis of atmospheric alkenes could be a potential source of sulfuric acid only in the absence of solar radiation (Boy et al., 2013;Mauldin et al., 2012); therefore, these reactions are expected to provide a minor contribution to the formation of gaseous sulfuric acid during the daytime in urban Beijing.
In this study, the gaseous sulfuric acid concentration was measured using a long time-of-flight chemical ionization mass spectrometer (LToF-CIMS) in February-March 2018, in urban Beijing.A number of atmospheric parameters were recorded concurrently, including the ultraviolet radiation B (UVB) intensity, the concentrations of O 3 , NO x , SO 2 , and HONO, and particle number size distributions.The objective of this study is to develop a robust daytime gaseous sulfuric acid concentration proxy for Beijing, a representative Chinese megacity with an urban atmospheric environment.

Ambient measurements
An intensive campaign was carried out from 9 February to 14 March 2018 on the fifth floor of a teaching building on the west campus of Beijing University of Chemical Technology (39 • 94 N, 116 • 30 E).This monitoring site is 2 km to the west of the West Third Ring Road and is surrounded by commercial properties and residential dwellings.Hence, this station can be regarded as a representative urban site.
The sulfuric acid concentration was measured using a LToF-CIMS (Aerodyne Research, Inc.) equipped with a nitrate chemical ionization source.Ambient air was drawn into the ionization source through a 1.6 m stainless-steel tube with a 3/4 inch diameter.A mixture of a 3 standard cubic centimeter per minute (sccm) ultrahigh purity nitrogen flow containing nitric acid and a 20 standard liter per minute (slpm) pure air flow supplied by a zero-air generator (Aadco 737, USA), together as a sheath flow, was introduced into a PhotoIonizer (Model L9491, Hamamatsu, Japan) to produce nitrate reagent ions.This sheath flow was then introduced into a coaxial laminar flow reactor concentric to the sample flow.Nitrate ions were pushed to the middle of the sample flow under an electric field and subsequently charged sample molecules.For example, the atmospheric H 2 SO 4 molecules would be charged by nitrate reagent ions NO − 3 (HNO 3 ) 0−2 and would mainly produce HSO − 4 ions (m/z = 96.9601Th), HSO − 4 q HNO 3 ions (m/z = 159.9557Th), and HSO − 4 q (HNO 3 ) 2 ions (m/z = 222.9514Th).In addition, HSO − 4 q H 2 SO 4 ions (m/z = 194.9275Th) were formed from ion-induced clustering of neutral sulfuric acid and bisulfate ions within the LToF-CIMS ion reaction zone, and also from the evaporation of DMA and the replacement of one molecule of H 2 SO 4 with one bisulfate ion, HSO − 4 , during the NO − 3 reagent ion charging of a stabilized neutral sulfuric acid dimer in the real atmosphere in the presence of DMA or a molecule that works in the same way as DMA.During the campaign, the sample flow rate was kept at 8.8 slpm, as mass flow controllers fixed the sheath flow rate and the excess flow rate, and the flow into the mass spectrometer (around 0.8 slpm) was fixed by the size of a pinhole between the ionization source and the mass spectrometer.The concentration of gaseous sulfuric acid was then determined by Eq. ( 2): where C is the calibration coefficient, and NO − 3 (HNO 3 ) 0−2 , HSO − 4 q (HNO 3 ) 0−2 and HSO − 4 q H 2 SO 4 represent the signals of corresponding ions in units of counts per second (cps).The unit of the resulting [H 2 SO 4 ] is molecule cm −3 .The CIMS was calibrated during the campaign using a homemade calibration box that can produce adjustable concentrations of gaseous sulfuric acid from SO 2 and OH radicals following the protocols in previous literatures (Kürten et al., 2012;Zheng et al., 2015).We obtain a calibration coefficient of 3.79×10 9 molecule cm −3 for our instrument and use 1.1 × 10 10 molecule cm −3 as the effective calibration coefficient, after taking the diffusion losses in the stainless-steel tube and the nitrate chemical ionization source into account.The obtained mass spectra were analyzed with a tofTools package based on MATLAB software (Junninen et al., 2010).
Ambient particle number size distributions down to about 1 nm were measured using a combination of a scanning mobility particle sizer spectrometer (SMPS) equipped with a diethylene glycol-based condensation particle counter (DEG-CPC, ∼ 1-10 nm) and a conventional particle size distribution system (PSD, ∼ 3 nm-10 µm) consisting of a pair of aerosol mobility spectrometers developed by Tsinghua University (Cai et al., 2017;Jiang et al., 2011;Liu et al., 2016).The values of the CS were calculated following Eq.(3) (Dal Maso et al., 2002): where D pi is the geometric mean diameter of particles in the size bin i and N i is the particle number concentration in the corresponding size bin; D is the diffusion coefficient of gaseous sulfuric acid, and β m represents a transition-regime correction factor dependent on the Knudsen number (Fuchs and Sutugin, 1971;Gopalakrishnan and Hogan Jr., 2011).
Y. Lu et al.: A proxy for atmospheric daytime gaseous sulfuric acid concentration SO 2 , O 3 , and NO x concentrations were measured using a SO 2 analyzer (Model 43i, Thermo, USA), an O 3 analyzer (Model 49i, Thermo, USA), and a NO x analyzer (Model 42i, Thermo, USA) with detection limits of 0.1, 0.5, and 0.4 ppbv, respectively.The above instruments were pre-calibrated before the campaign.The UVB (280-315 nm) intensity (UV-S-B-T, KIPP&ZONEN, the Netherlands) was measured on the rooftop of the building.Atmospheric HONO concentrations were measured using a homemade HONO analyzer with a detection limit of 0.01 ppbv (Tong et al., 2016).
Particle number size distributions and concentrations of gaseous sulfuric acid, SO 2 , O 3 , NO x , and HONO were recorded with a time resolution of 5 min, and the UVB intensity with a time resolution of 1 min.A linear interpolation method was used for deriving the variables with the same time intervals, i.e., 5 min.Only data between local sunrise and sunset were used in the subsequent analysis.

Development of a proxy for atmospheric gaseous sulfuric acid
We derived the gaseous sulfuric acid concentration proxy on the basis of currently accepted formation pathways of sulfuric acid in the atmosphere (Reactions R1-R3; Finlayson-Pitts and Pitts, 2000; Stockwell and Calvert, 1983): Reaction (R1) is the rate-limiting step of this formation pathway (Finlayson-Pitts and Pitts, 2000), so our proxy will consider the two major processes that determine the abundance of gaseous sulfuric acid: the formation of gaseous sulfuric acid from reactions between SO 2 and OH radicals, and the loss of gaseous sulfuric acid due to its condensation onto preexisting particles (Dal Maso et al., 2002;Kulmala et al., 2012;Pirjola et al., 1999).
The rate of change of sulfuric acid concentration can be written as Eq. ( 4) (Mikkonen et al., 2011): where k is a temperature-dependent reaction constant given by Eq. ( 5) (DeMore et al., 1997;Mikkonen et al., 2011). where represents the density of the air in molecule cm −3 as calculated by 0.101 To simplify the calculation, the production and loss of sulfuric acid can be assumed to be in pseudo-steady-state (Mikkonen et al., 2011;Petäjä et al., 2009).Then the sulfuric acid concentration can be written as follows: Atmospheric OH radical measurements also represent a major challenge.As previous studies suggest that the OH radical concentration is strongly correlated with the intensity of UVB, [OH] could be replaced with UVB intensity in the proxy equation (Petäjä et al., 2009;Rohrer and Berresheim, 2006).Although the photolysis of O 3 (λ < 320 nm) and subsequent reactions with H 2 O are considered to be the dominant source of OH radicals in the atmosphere (Logan et al., 1981), recent studies argue that photolysis of HONO (λ < 400 nm) is a potentially important OH radical formation pathway (Hendrick et al., 2014;Kleffmann, 2007;Su et al., 2011;Villena et al., 2011).Thus, we attempt to introduce both O 3 and HONO into the proxy equation and evaluate their effects on the concentration of OH radicals.
In practice, the exponents for variables in nonlinear fitting procedures are rarely equal to one (Mikkonen et al., 2011), so we replaced the factors x i with x w i i in the proxy, where x i can be an atmospheric variable such as UVB and [SO 2 ], and w i defines x i weight in the proxy.As k is a temperaturedependent reaction constant and varies within a 10 % scope in the atmospheric temperature range of 267.6-292.6K, i.e., the actual atmospheric temperature variation in this study, we approximately regard k as a constant and use a new scaling factor k 0 .This methodology has been used previously in the proxies of gaseous sulfuric acid in Hyytiälä, southern Finland (Petäjä et al., 2009).As a result, the general proxy equation can be written as Eq. ( 7), with the UVB intensity, [SO 2 ], the condensation sink (CS), [O 3 ], and [HONO] (or [NO x ]) as predictor variables: The nonlinear curve-fitting procedures using iterative least square estimation for the proxies of gaseous sulfuric acid concentration based on Eq. ( 7) were performed using a custom-made MATLAB software.In addition to the correlation coefficient (R), the relative error (RE) is used to evaluate the performance of proxies in the statistical analysis and can be written as follows:  (Hendrick et al., 2014;Spataro et al., 2013;Wang et al., 2017), which is a factor of 4-10 higher than HONO concentrations measured in Europe (Alicke et al., 2002(Alicke et al., , 2003)).In addition, Beijing is dry in winter and had a mean ambient relative humidity of 28 % during the campaign.

Correlations between [H 2 SO 4 ] and atmospheric variables
Table 2 summarizes the correlation coefficients between [H 2 SO 4 ] and atmospheric variables using a Spearman-type correlation analysis.Clearly, the UVB intensity is an isolated variable that is independent of all the other variables but that imposes a positive influence on O 3 due to the photochemical formation of ozone, and a negative influence on HONO because of HONO's photochemical degradation.(Calvert et al., 1994) or homogeneously formed by the gas-phase NO + OH reaction, the former of which is likely dominate for the daytime HONO production in urban Beijing (Liu et al., 2014).
As the UVB intensity and [SO 2 ] have been reported as the dominating factors for the formation of sulfuric acid (Mikkonen et al., 2011;Petäjä et al., 2009), we further explored the relationship of the measured sulfuric acid concentrations with the UVB intensity and [SO 2 ] using a nonlinear curvefitting method with a single variable.Figure 1a presents a scatterplot of [H 2 SO 4 ] against the UVB intensity, colorcoded by [SO 2 ].A good correlation between [H 2 SO 4 ] and the UVB intensity, laying with [SO 2 ], is evident.A similar scatterplot (Fig. 1b) of [H 2 SO 4 ] against [SO 2 ], color-coded by the UVB intensity, leads to a similar conclusion.

Proxy construction
Similar to the nonlinear proxies suggested by Mikkonen et al. (2011), we tested a number of proxies for gaseous sulfuric acid; these proxies are listed in Table 3 with their respective fitting parameters, and their performance summarized in Table 4.The scatterplots of observed [H 2 SO 4 ] versus predicted values given by proxies are presented in Fig. S1 in the Supplement.In these proxies, the concentration of a gaseous species is in molecule cm −3 , the UVB intensity is in W m −2 , the CS is in s −1 , and k 0 is a scaling factor.
Proxy N1 was built using the UVB intensity and [SO 2 ] as the source terms and the CS as the sink term, which follows the conventional idea of H 2 SO 4 formation and loss in the atmosphere.The CS was then removed from this proxy to examine the performance of proxy N2, which then had the UVB intensity and [SO 2 ] as the only predictor variables.As the formation of OH radicals in the atmosphere depends on precursors in addition to UVB, we further attempted to introduce the OH precursor term into the H 2 SO 4 proxy.Proxies N3 and N4 were built by introducing O 3 as the only OH precursor to evaluate its influence on the formation of sulfuric acid.Furthermore, we added HONO as another potential precursor for OH radicals, resulting in proxies N5 and N6.Lastly, proxy N7 was built by replacing [HONO] with [NO x ] because (1) HONO is not regularly measured, and (2) a good linear correlation between [HONO] and [NO x ] was generally observed during the daytime in this campaign, although higher [HONO]/[NO x ] ratios were noted in the morning due to the accumulation of HONO during the night (Fig. 2).The RH was not considered in the current study because a test that introduced RH into the proxies did not result in significantly better performance, which is consistent with conclusions from Mikkonen et al. (2011).As shown in Table 4, the correlation coefficients are in the range of 0.83-0.86 and REs are in the range of 19.1 %-20.0 %.The exponents for the UVB intensity range from 0.13 to 0.15, and those for [SO 2 ] generally range from 0.38 to 0.41, except in case of proxy N6 (b = 0.33).The obtained exponent b for [SO 2 ] is significantly smaller than 1 unlike the assumption in Eq. ( 6), mainly because [SO 2 ] is also an indicator of air pollution that usually influences the sinks of both OH radicals and sulfuric acid.The exponent for [SO 2 ] ranged from 0.5 to 1.04 in the previous proxy study for European and North American sites (Mikkonen et al., 2011), including values from 0.48 to 0.69 in Atlanta, GA, USA, which was probably quite a polluted site as the measurements were conducted only 9 km from a coal-fired power plant.The obtained value range of the exponent b for [SO 2 ] in our study is probably related to the urban nature of Beijing.The value of exponent c for the CS in proxy N1 is as low as 0.03, which might be due to the covariance of the CS and certain H 2 SO 4 sources that cancels the dependence on the CS, or it might indicate that the CS is actually insufficient in regulating the H 2 SO 4 concentration, as recently suggested by Kulmala et al. (2017).By comparing proxies N1 and N2, we can see that the CS plays a minor role because the exponents of [SO 2 ] and UVB, the overall correlation coefficient, and the REs are almost identical with and without the CS.We can also see the negligible role of the CS when comparing the results of proxies N3 and N4 where O 3 is considered.However, the role of the CS becomes evident between proxies N5 and N6 when HONO is introduced: the exponents of [SO 2 ], [O 3 ], and [HONO] significantly increased when taking the CS into account, suggesting that the covariance between HONO and the CS can explain, at least partially, the close-to-zero exponent of the CS in proxies N1-N4.In addition, when [O 3 ] is introduced as the only precursor for OH radicals, minor Table 3. Proxy functions for the nonlinear fitting procedure.

Proxy Function
is the concentration of sulfur dioxide in molecule cm −3 ; CS is the condensation sink in s −1 ; [O 3 ] is the concentration of ozone in molecule cm −3 ; [HONO] is the concentration of nitrous acid in molecule cm −3 ; [NO x ] is the concentration of nitrogen oxides in molecule cm −3 ; and k 0 is a scaling factor.The superscript numbers define the variables' weight in the proxy, and their values for the best fit are shown in Table 4. improvements in the correlation coefficient and RE were obtained, as suggested by comparing proxies N3 and N1.When both [O 3 ] and [HONO] were introduced as OH precursors in proxies N5-N7, the REs show noticeable improvements, and correlation coefficients improved slightly.Altogether, these observations suggest that it is crucial to introduce HONO into the proxy, both in our study and also likely for the previous work where the exponent of the CS is close-to-zero (Mikkonen et al., 2011).
Although so far proxy N5 had the best fitting quality, it is impractical to explicitly include [HONO] because HONO measurements are very challenging.As shown in Fig. 2, [HONO] and [NO x ] tended to correlate linearly with each other during the daytime in this campaign, with a linearly fitted [HONO]/[NO x ] ratio of around 0.03 and a relative error   of 0.42.Occasionally higher [HONO]/[NO x ] ratios could be seen in the morning, which may have been due to the fact that the HONO concentration could have undergone an accumulation process during the nighttime and lead to a deviation from the steady state.Therefore, due to the good correlation, proxy N7 replaces [HONO] by [NO x ], a more easily measured variable, and performs as well as proxy N5.
Clearly, proxy N2 provides the simplest parameterization, but proxies N5 and N7 result in the best fitting quality due to the introduction of [HONO].Figure 3 presents the RE values for proxies N2 and N7, respectively, as a function of linear bins of measured sulfuric acid concentrations.The performance of proxy N7 is considerably better than that of proxy N2 in the sulfuric acid concentration range of (2.2-10) × 10 6 molecule cm −3 , which covers the 5 %-95 % percentiles of the sulfuric acid concentration in this study.In the worst scenario, the RE of proxy N2 is 1.2 times as high as that of proxy N7, e.g., the REs are 16.75 % and 13.99 %, respectively, in the sulfuric acid concentration bin of (4-5) × 10 6 molecule cm −3 , and 16.71 % and 14.42 %, respectively, in the bin of (7-8) × 10 6 molecule cm −3 .

Comparison of measured and predicted [H 2 SO 4 ]
A comparison between measured and predicted [H 2 SO 4 ] was performed.Figure 4 includes calculated results from proxies N2 and N7 as well as from a proxy constructed according to measurement in a boreal forest site, Finland, i.e., Eq. (1) (Petäjä et al., 2009).The measured daytime [H 2 SO 4 ] on 10 March 2018, was above 4×10 6 molecules cm −3 when averaged to a time resolution of 5 min.The predicted [H 2 SO 4 ] using proxies N2 and N7 both track the measured [H 2 SO 4 ] quite well, even when an unexpected dip in the sulfuric acid concentration was observed at around 10:00-11:00.The performance of proxy N7 is better than that of proxy N2 during the entire day, consistent with our results in Fig. 3.The proxy by Petäjä et al. (2009) underestimated the concentrations of sulfuric acid by a factor of 20 or so, which can be attributed to the very different values of the CS between Beijing and the boreal forest.The fact that [H 2 SO 4 ] Petäjä et al. does not track the measured [H 2 SO 4 ] even after including a scaling factor indicates that proxies are site-specific and do not necessarily work well in locations other than the location that they were originally developed for.In addition, the direct performance comparison between proxy N2 and the proxy by Petäjä et al. (2009) indicates the importance of assigning exponential weights to variables in the nonlinear fitting procedures, which is consistent with results from Mikkonen et al. (2011).We consider this last proxy more reasonable than the others due to the following reasons: it makes the equation physically meaningful as the CS is included as a sink term, and the RE is reduced considerably compared with the other proxies.
Overall, this suggests that the photolysis of O 3 and HONO are both important OH sources in urban Beijing.In summary, we recommend using the simplest proxy (proxy N2 as shown in Eq. 9) and a more accurate proxy (proxy N7 as shown in Eq. 11) for calculating daytime gaseous sulfuric acid concentrations in the urban Beijing atmosphere.It is clear that the current proxies are only based on a month-long campaign of sulfuric acid measurements in urban Beijing during winter.Given the dramatic reduction in the concentration of SO 2 in recent years (Wang et al., 2018) and the strong dependence of calculated [H 2 SO 4 ] on [SO 2 ], the performance of the proxies in past and future years remains to be evaluated.Furthermore, the proxies might be site-specific and season-specific.As the proxies were derived with atmospheric parameters in winter, in urban Beijing, the exponents for atmospheric variables in the proxy could have different values for other cities or other seasons.Thus, the proxies in this study should be further tested before their application to other Chinese megacities or other seasons.
Data availability.The data in this article are available from the corresponding author upon request (lin_wang@fudan.edu.cn).
Supplement.The supplement related to this article is available online at: https://doi.org/10.5194/acp-19-1971-2019-supplement.Author contributions.LW designed this study.YLu, CY, YF, YC, YLiu, GY, YW, YZ, RY, RB, and CD conducted the field campaign.YLu analyzed data with contributions from LW and all of the other co-authors.YLu and LW wrote the paper with contributions from all of the other co-authors.
Competing interests.The authors declare that they have no conflict of interest.

Figure 1 .
Figure 1.Correlations (a) between [H 2 SO 4 ] and UVB intensity, and (b) between [H 2 SO 4 ] and [SO 2 ] during the 9 February to 14 March 2018 campaign.k m is a constant term.

Figure 2 .
Figure 2. Correlation between [HONO] and [NO x ] during the 9 February to 14 March 2018 campaign.The black line represents a linear fit with a zero intercept.

Figure 3 .
Figure 3. Performance assessments of proxy N2 and proxy N7.The REs are used to evaluate the performances of proxy N2 and N7, respectively, as a function of linear bins of measured sulfuric acid concentrations.
proxy,i − [H 2 SO 4 ] meas.,i| , is usually accompanied by a high concentration of both SO 2 and NO x in urban China, indicating co-emissions.A strong correlation between [HONO] and [NO x ] (r = 0.88) in our measurement is supported by the fact that HONO can either be heterogeneously formed by reactions of NO 2 on various surfaces 4 ] and the CS, which is in contrast to the conventional thought that the CS describes the loss of gaseous sulfuric acid molecules onto preexisting particles and thus should show a negative correlation.The CS correlates well with [SO 2 ] (r = 0.83) and [NO x ] (r = 0.77); thus, a high CS value, as an indicator of atmospheric particle pollu-tion

Table 2 .
Correlation coefficients (Spearman type) between [H 2 SO 4 ] and atmospheric variables during the daytime.

Table 4 .
Results of the nonlinear fitting procedure for different proxy functions, in addition to the correlation coefficient (R, Pearson type) and the relative error (RE).
As HONO measurements are not a regular practice, we can further replace [HONO] with [NO x ], shown in Eq. (11), which can be justified by the strong linear correlation between [HONO] and [NO x ] observed in this study.The unit for [NO x ] is molecule cm −3 .