Technical note : Evaluation of the simultaneous measurements of mesospheric OH , HO 2 , and O 3 under a photochemical equilibrium assumption – a statistical approach

This Technical Note presents a statistical approach to evaluating simultaneous measurements of several atmospheric components under the assumption of photochemical equilibrium. We consider simultaneous measurements of OH, HO2, and O3 at the altitudes of the mesosphere as a specific example and their daytime photochemical equilibrium as an evaluating relationship. A simplified algebraic equation relating local concentrations of these components in the 50–100 km altitude range has been derived. The parameters of the equation are temperature, neutral density, local zenith angle, and the rates of eight reactions. We have performed a one-year simulation of the mesosphere and lower thermosphere using a 3-D chemical-transport model. The simulation shows that the discrepancy between the calculated evolution of the components and the equilibrium value given by the equation does not exceed 3–4 % in the full range of altitudes independent of season or latitude. We have developed a statistical Bayesian evaluation technique for simultaneous measurements of OH, HO2, and O3 based on the equilibrium equation taking into account the measurement error. The first results of the application of the technique to MLS/Aura data (Microwave Limb Sounder) are presented in this Technical Note. It has been found that the satellite data of the HO2 distribution regularly demonstrate lower altitudes of this component’s mesospheric maximum. This has also been confirmed by model HO2 distributions and comparison with offline retrieval of HO2 from the daily zonal means MLS radiance.


Introduction
A prominent feature of atmospheric photochemical systems is the presence of a large number of chemical components with short lifetimes and concentrations close to stable photochemical equilibrium at every instant.The condition of balance between their sources and sinks is described by a system of algebraic equations.This system can be used to determine characteristics of hard-to-measure atmospheric species through other measurable components, evaluate results of remote or in situ measurements, estimate reaction rates usually known with significant uncertainty, and to understand processes and chemical reactions that influence the variability in the most important atmospheric components, e.g., ozone, in the geographical region of interest.
This approach has found the following wide applications: 1.In 3-D chemical transport models that include a large set of physical and chemical processes with a broad spectrum of spatio-temporal scales.In particular, the chemical family concept is widely used for simulating gas-phase photochemistry of the lower and middle atmosphere (e.g., Douglass et al., 1989;Kaye and Rood, 1989;Rasch et al., 1995), when transport is taken into account only for the concentration of a chemical family, while relative concentrations of the constituent fast components are calculated from the instantaneous stable equilibrium condition.Complemented with Henry's law (e.g., Djouad et al., 2003;Tulet et al., 2006) in multiphase models, this approach markedly saves calculation time and increases the overall stability of the numerical scheme.Moreover, use of the photochemical Published by Copernicus Publications on behalf of the European Geosciences Union.
2. In investigations of the chemistry of the surface layer and free troposphere in different regions (over megalopolises, in rural areas, in the mountains, over the seas) based on measurements of nitrogen species, peroxy radicals, ozone, aerosols, and other components aimed at understanding processes impacting the surface ozone formation and air quality.The equilibrium condition is most frequently used for nitrogen species.For example, Chameides (1975) proposed a model for determining the vertical distribution of odd nitrogen, in which the HNO 3 profile could be deployed to retrieve profiles of five other components (NO, NO 2 , NO 3 , N 2 O 5 , and HNO 2 ) from their photochemical equilibrium condition.In the paper by Stedman et al. (1975), the equation for NO 2 equilibrium that accounted only for the main source and sink of this component was applied to determine the photodissociation constant J (NO 2 ).A more accurate equation for the NO 2 equilibrium was used by Crawford et al. (1996) and Koike et al. (1996) to determine NO 2 /NO partitioning and NO x , allowing for, in particular, an investigation of the spatial distribution of NO x /NO y over the Pacific.
Nighttime equilibrium in the NO 2 -NO 3 -N 2 O 5 system is used to determine surface layer N 2 O 5 concentration, equilibrium constant of this system, equilibrium partitioning between NO 3 and N 2 O 5 , and loss coefficients of NO 3 , N 2 O 5 , and NO x (Martinez et al., 2000;Brown et al., 2003;Crowley et al., 2010;McLaren et al., 2010;Benton et al., 2010;Sobanski et al., 2016).Ko et al. (2003), Cantrell et al. (2003), and Penkett et al. (1997, 1998) algebraic expressions derived from equilibrium conditions for H 2 O 2 , peroxy radicals, and nitrogen species were used to determine equilibrium values of peroxide concentration, total peroxy radical level, and NO/NO 2 ratio, and to diagnose the ozone production and loss levels in a clean or polluted troposphere.

3.
In stratospheric chemistry studies, including determination of a critical parameter in catalytic cycles of ozone destruction in the polar stratosphere.In particular, the equilibrium condition for ClO and Cl 2 O 2 along with the measurement data of daytime and nighttime concentrations of these components in the polar stratosphere are used to evaluate the temperature dependence of the ClO concentration, reaction constants determining the ClO + ClO + M ↔ Cl 2 O 2 + M equilibrium, and the photolysis rate of Cl 2 O 2 (Ghosh et al., 1997;Avallone et al., 2001;Solomon et al., 2002;Stimpfle et al., 2004;von Hobe et al., 2005von Hobe et al., , 2007;;Berthet et al., 2005;Butz et al., 2007;Kremser et al., 2011;Sumińska-Ebersoldt et al., 2012;Wetzel et al., 2012).Pyle et al. (1983) proposed a method for derivation of the OH concentration from satellite infrared measurements of NO 2 and HNO 3 using a simple algebraic relation following from the equilibrium condition for HNO 3 .Algorithms for retrieving distributions of OH and HO 2 from the satellite measurement data of O 3 , NO 2 , H 2 O, HNO 3 by LIMS/Nimbus-7 and UARS with the help of algebraic models following from the photochemical equilibrium of O x , HO x , and HNO 3 components were proposed by Pyle and Zavody (1985) and Pickett and Peterson (1996).It is also noteworthy that similar models are widely used for calculating concentrations of components with a short lifetime (e.g., O( 1 D) and OH) and subsequently evaluating vertical distributions of eddy diffusivity from measurements of trace gas concentration profiles (e.g., see Massie and Hunten, 1981).Kondo et al. (1988) made use of the photochemical equilibrium between NO and NO 2 for understanding diurnal variations in NO concentration measured during aircraft flights.In the paper by Webster et al. (1990) simultaneous in situ balloon-borne measurements of NO, NO 2 , HNO 3 , O 3 , and N 2 O and the photochemical equilibrium condition for various nitrogen components were used to determine OH, N 2 O 5 , and NO y concentrations.A similar approach was employed by Kawa et al. (1990) 4. In investigations of the chemistry of O x -HO x components and atmospheric glows in the mesosphere to lower thermosphere (MLT) area.In particular, Kulikov et al. (2006Kulikov et al. ( , 2009a) ) proposed algorithms for the simultaneous retrieval of O, H, HO 2 , and H 2 O from joint OH and O 3 satellite measurement, in which the assumption of photochemical equilibrium of O 3 , OH, and HO 2 was utilized.For several decades the assumption of the photochemical equilibrium of ozone (PEO) was widely used to determine distributions of atomic oxygen and atomic hydrogen at altitudes of the MLT via satellite and rocket measurements of ozone concentration and airglow emissions (e.g., Evans and Llewellyn, 1973;Good, 1976;Pendleton et al., 1983;McDade et al., 1985;McDade and Llewellyn, 1988;Evans et al., 1988;Thomas, 1990;Llewellyn et al., 1993;Llewellyn and McDade, 1996;Mlynczak et al., 2007Mlynczak et al., , 2013aMlynczak et al., , b, 2014;;Smith et al., 2010;Siskind et al., 2008Siskind et al., , 2015)).Russell and Lowe (2003) applied PEO to infer the seasonal and global climatology of atomic oxygen using WINDII/UARS.PEO was deployed to investigate hydroxyl emission mechanisms, morphology, and variability in the MLT region (Marsh et al., 2006;Xu et al., 2010Xu et al., , 2012;;Kowalewski et al., 2014).Mlynczak andSolomon (1991, 1993) and Mlynczak et al. (2013b) used the equilibrium assumption to derive exothermic chemical heat.The PEO assumption was employed for studying the mesospheric OH* layer response to gravity waves (Swenson and Gardner, 1998).In ultimately theoretical works, e.g., Grygalashvyly et al. (2014) and Grygalashvyly (2015), PEO was used to derive the dependence of excited hydroxyl layer concentration and altitude on atomic oxygen and temperature.In the paper by Sonnemann et al. (2015) it was used to analyze annual variations of the OH* layer.Moreover, PEO is frequently applied implicitly, when authors are equating the nighttime loss of ozone in the reaction with atomic hydrogen and production of ozone by a 3-body reaction of molecular and atomic oxygen (e.g., Nikoukar et al., 2007).
In the present Technical Note we demonstrate how the photochemical equilibrium condition of several atmospheric components may be employed to statistically validate data of their simultaneous measurements, particularly in the case when measurement error is large.
We consider the simultaneous photochemical daytime equilibrium of OH, HO 2 , and O 3 at the altitudes of the mesosphere.We have derived a simplified algebraic equation describing the relationship between local concentrations of the components at the altitudes of 50-100 km.The only parameters of the equation are temperature, neutral density, local zenith angle, and constants of eight reactions.A one-year simulation of the mesosphere and lower thermosphere based on a 3-D chemical-transport model shows that the discrep-ancy between the calculated evolution of the components and the equilibrium value given by the equation does not exceed 3-4 % in the full range of altitudes independent of season or latitude.
We have developed a statistical Bayesian evaluation technique for simultaneous measurement of OH, HO 2 , and O 3 based on the mentioned equilibrium equation taking into account the measurement error.The first results of its application to MLS/Aura data (Wang et al., 2015a, b;Schwartz et al., 2015) are presented.It is found that the satellite data of the HO 2 distribution regularly demonstrate lower altitudes of this component's mesospheric maximum.These results confirm the ones obtained via the offline retrieval of HO 2 from the MLS primary data (Millán et al., 2015).
This Technical Note is structured as follows.A 3-D chemical transport model is briefly described in Sect. 2. In Sect.3, a simplified algebraic relationship between the equilibrium concentrations of OH, HO 2 , and O 3 is derived and verified by 3-D simulations.Section 4 presents the method of statistical evaluation of simultaneous data of OH, HO 2 , and O 3 .The results of applying the method to MLS/Aura data are presented in Sect. 5.The last section contains a discussion of the results followed by concluding remarks.

Model and calculations
For our calculations we used the global 3-D chemical transport model (CTM) of the middle atmosphere developed by the Leibniz Institute of Atmospheric Physics (IAP; e.g., Berger, 1994;Ebel et al., 1995;Kremp et al., 1999;Berger and von Zahn, 1999;Hartogh et al., 2004Hartogh et al., , 2011;;Sonnemann et al., 1998Sonnemann et al., , 2006Sonnemann et al., , 2007)).It was designed particularly for investigation of the spatio-temporal structure of phenomena in the MLT region and specifically in the extended mesopause region.The grid-point model extends from the ground up to the middle thermosphere (0-150 km; 118 pressure-height levels).The horizontal resolution amounts to 5.625 • latitudinally and 5.625 • longitudinally.The chemical module described in numerous papers (e.g., Sonnemann et al., 1998;Körner and Sonnemann, 2001;Grygalashvyly et al., 2009Grygalashvyly et al., , 2011Grygalashvyly et al., , 2012) ) consists of 19 constituents, 49 chemical reactions, and 14 photo-dissociation reactions (see Table 1).The reaction rates used in the model are taken from Burkholder et al. (2015).The temperature-dependent reaction rates are calculated online; thus, they are sensitive to small temperature fluctuations.We make use of the pre-calculated dissociation rates (Kremp et al., 1999).
The evolution of the components of the HO x (H, OH, HO 2 , H 2 O 2 ) and NO x (N, NO, NO 2 , NO 3 ) families is calculated using the chemical family concept proposed by Shimazaki (Shimazaki, 1985).This is done because of the presence of short-lived components among these families, with lifetimes much shorter than those of the families themselves, which imposes significant restrictions on the value of the CTM's For example, the daytime lifetimes of OH and HO 2 above 70 km are about 1 s or less, while the lifetime of the HO x family is about 10 4 s or more.Therefore, when calculating these components individually it is necessary to set the CTM's integration step to be much less than 1 s.In our work, the Shimazaki technique is applied for calculating the evolution of each component of the HO x and NO x families.We emphasize that this technique does not explicitly use the steady-state approximation for the components, instead it utilizes the approach based on an implicit Euler scheme (see Shimazaki, 1985).This allows for increasing the integration step of the CTM significantly without loss of accuracy in calculating the short-lived components.In our work the integration time is chosen to be 9 s.The model includes 3-D advective and vertical diffusive transport (turbulent and molecular).Three-dimensional fields of temperature and winds are taken from the Canadian Middle Atmosphere Model (CMAM) for the year 2000 (de Grandpre et al., 2000;Scinocca et al., 2008).We use the Walcek scheme (Walcek and Aleksic, 1998;Walcek, 2000) for advective transport and the implicit Thomas algorithm as described in Morton and Mayers (1994) for diffusive transport.The vertical eddy diffusion coefficient is based on the results by Lübken (1997).
We calculate the annual variation in spatio-temporal distributions of OH, HO 2 , and O 3 and constructed distributions of the F (OH, HO 2 , O 3 ) function introduced in Sect. 1.To remove transitional regions that correspond to sunset and sunrise, we only take into account periods of local time with the solar zenith angle χ < 85 • .The obtained results are pre-sented in the model coordinates, so the pressure-height levels are used for the vertical axes.In addition, the approximate altitudes are shown in the figures of Sect. 1, calculated for a given month utilizing averaged temperature profiles of the model and hydrostatic equilibrium.
The daytime balance of OH concentration at mesospheric altitudes is determined by the following primary reactions (Brasseur and Solomon, 2005): The daytime balance of the HO 2 concentration is Atmos.Chem.Phys., 18, 7453-7471, 2018 www.atmos-chem-phys.net/18/7453/2018/The daytime balance of the O 3 concentration is Expressions for local concentrations of OH, HO 2 , and O 3 in the photochemical equilibrium are written in the form where k i are the corresponding reaction constants from Burkholder et al. (2015).
We eliminate O and H from Eqs.
Almost everywhere in the mesosphere and lower thermosphere (with the exception of 85-95 km, see Kulikov et al., 2017) the photodissociation is the main ozone sink, i.e., k 52 + k 53 k 21 • H. Therefore, in the zero order approximation, Eq. ( 3) can be simplified and the concentration of atomic oxygen can be defined in terms of ozone concentra-tion: Making use of Eq. ( 4) we can derive, from Eq. ( 2), an expression for the concentration of H in terms of concentrations of OH, HO 2 , and O 3 : By substituting this equation and Eq.(4) into Eq.( 1) we obtain an expression relating OH, HO 2 , and O 3 : Figure 1 shows height-latitude cross-sections of F (OH, HO 2 , O 3 ) for each month (in this section an-gle brackets denote monthly averaged zonal mean values).
The gray area corresponds to χ > 85 • .One can see that Eq. ( 6) is most accurate within the 50-76 km range and above 86 km, where |< F > −1| ≤ 1 %.The difference reaches 3-4 % in the region between 76 and 86 km.The altitude of this region has an annual variation with a maximum deviation in the winter hemisphere.Below 50 km the value of < F > increases up to 1.2 at 40 km, thus below the stratopause Eq. ( 6) no longer describes the simultaneous photochemical equilibrium of OH, HO 2 , and O 3 .Note that these components remain short-lived below 50 km (with the lifetimes of about 10 2 -10 3 s; Brasseur and Solomon, 2005) depending on height and duration of daylight.However, for a quantitative description of their daytime equilibrium, it is necessary to include additional reactions involving, in particular, the components of the NO x family.Note also that Eqs. ( 1) and ( 6) take into account only the main daytime source of OH (P OH ) specified by reactions (R18), (R14), and (R21): These reactions run "inside" the HO x (H, OH, HO 2 , H 2 O 2 ) family and do not perturb its total concentration.The heightlatitude cross-sections of < P OH > for each month are presented in Fig. 2. The next important daytime source of OH is specified by reactions (R59) and (R7) involving H 2 O, the main source for the HO x family: Figure 3 shows height-latitude cross-sections of < P H 2 O OH /P OH > for each month.Comparing Figs. 1 and 3, we conclude that the previously indicated 3-4 % deviation of < F > from 1 in the region between 76 and 86 km is largely due to the neglect of these reactions.
Another source of OH is sporadically activated during charged particle precipitation events and exists for a relatively short time (several days).Solar proton events (SPEs) perturb the ionic composition in the mesosphere and the upper stratosphere considerably and trigger a whole cascade of reactions involving ions, neutral components, and their clusters (e.g., O + 2 • H 2 O).This leads to an additional (to reactions R59 and R7) conversion of H 2 O molecules into OH and H (Solomon et al., 1981).The maximum of the OH production rate (P SPE OH ) induced by SPE is located in the polar latitudes in the region of 60-80 km and, as a rule, does not exceed 2×10 3 cm −3 s −1 (Jackman et al., 2011(Jackman et al., , 2014)).It can be seen from Fig. 2 that at these latitudes and altitudes the P SPE OH /P OH ratio does not exceed 1-2 %, even for the maximum values of P SPE OH .This means that the impact of P SPE OH on Eq. ( 6) is of the same order of smallness as in the case of reactions (R59) and (R7); hence, it may be neglected.A similar conclusion can be made for other reactions from Table 1, not accounted for by Eq. ( 6), including the ones involving NO x in both quiet and perturbed conditions in the mesosphere.
4 Method of statistical evaluation of simultaneous measurement of OH, HO 2 , and O 3 The proposed method is based on the statistical Bayesian procedure described in the works by Kulikov et al. (2009b) and Nechaev et al. (2016).It was originally developed for retrieving trace gas concentrations in the mesosphere from ground-based and satellite measurements of other mesospheric components.With respect to the considered evaluation problem, this procedure consists of three steps: (i) constructing a conditional probability density function (PDF) of OH, HO 2 , and O 3 concentration values at each altitude z in the selected interval assuming that there is certain measurement data of these components and the algebraic relationship (6) is valid; ii) calculating the first moments of this distribution, i.e., expected value and dispersion of each component using the Metropolis-Hastings algorithm (Chib and Greenberg, 1995) for multidimensional integration; (iii) comparing the obtained results with the initial measurement data.
For constructing posterior PDF it is convenient to introduce vector u HO ret 2 , O ret 3 , OH ret , whose components are the retrieved values of chemical species concentrations at a certain altitude z, and vector x HO m 2 , O m 3 , OH m composed of experimentally measured values of the components of vec- tor u, x j = u j + ξ j , j = 1., where ξ j is a random error of measuring the j th component of vector u at the altitude z.It is assumed that 1. random variables ξ j are distributed normally with densities 2. ξ j are mutually independent: where W ξ (ξ ) is the total PDF of all ξ j ; 3. dispersions σ j in Eq. ( 7), that are expected error values, are assumed to be known a priori (in our case they are provided by the MLS retrieval algorithm along with measured data).
Then the probability to observe vector x is given by the conditional PDF where δ (...) is the delta function.
www.atmos-chem-phys.net/18/7453/2018/Atmos.Chem.Phys., 18, 7453-7471, 2018 The prior relationship of HO ret 2 , O ret 3 , and OH ret concentrations (Eq.6) can be written as u 3 = G (u 1 , u 2 ).Integrating the left-hand side of Eq. ( 9) with the conditional PDF of the variable u 3 : yields a likelihood function of the model According to Bayes' theorem, the posterior function, i.e., the probability density of latent variables u 1 and u 2 , under the condition that x is observed, is defined by the expression in which P apr (u 1 , u 2 ) defines the prior PDF of u 1 and u 2 .The retrieved value of the latent variable u 1,2,3 is hereinafter understood as the mean value of the function in Eq. ( 11): Its dispersion defines the uncertainty in the retrieval: where the angle brackets denote averaging in the sense of Eq. ( 12).

MLS/Aura data evaluation and results
We used the latest version (v4.2) of the MLS "standard" product (Livesey et al., 2017) for trace gas concentrations and temperature T within the 1-0.046 mbar pressure interval, where all data are suitable for scientific applications (Wang et al., 2015a, b;Schwartz et al., 2015).We took the daytime data when the solar zenith angle χ < 80 • for January, May, and September 2005.All data were appropriately screened."Pressure", "estimated precision", "status flag", "quality", "convergence", and "clouds" fields were taken into account.HO 2 data were seen as the day-minus-night difference as prescribed by the MLS data guidelines (Livesey et al., 2017).Following Pickett et al. (2008), each daytime profile of this component measured on a given day at a latitude Lat, a profile resulting from averaging the nighttime profiles of HO 2 , measured on the same day in the latitude range of Lat ± 5 • , was subtracted.This operation eliminates systematic biases affecting HO 2 retrievals, but limits the studied latitude range to the one where MLS observes both daytime and nighttime data.
The integrals in Eqs. ( 12)-( 13) were calculated at every pressure level p for each set of simultaneously measured vertical profiles OH MLS (p), HO MLS distributed with normalized probability density given by Eq. ( 11) with P apr (u 1 , u 2 ) ≡ 1 was generated with the help of the Metropolis-Hastings algorithm (Chib and Greenberg, 1995).In this case, the statistical moments in Eqs. ( 12)-( 13) were determined by summation over the sample.
A typical example of retrieved profiles HO ret 2 , O ret 3 , and OH ret (black curves) in comparison with the measured HO MLS .Second, the results of a single measurement of all three components and their retrieved values have considerable uncertainties relative to their means within the whole interval of altitudes.Therefore, the observed and retrieved data should be compared using the commonly accepted approach (e.g., Pickett et al., 2008) of averaging large ensembles of profiles within certain latitude and time ranges, or zones.It is supposed that the noise of satellite measurement instruments is delta-correlated, so that random values corresponding to each single measured or retrieved profile are statistically independent.In this case the dispersion of a measured or retrieved zonal mean profile is determined by summation where N is the number of measured or retrieved profiles within the zone and σ 2 k is the dispersion of the kth measured or retrieved profile.
The range of latitudes covered by the satellite trajectory was divided into 17 bins 10 • each.About 3000 single profiles of each chemical component fall into one bin during a month of MLS/Aura observations.Therefore, the resulting uncertainties due to measurement noise in OH, HO 2 , and O 3 concentration profiles (both measured and retrieved) averaged over such ensembles are significantly lower (about one and a half orders of magnitude) than the uncertainties in individual profiles.Examples of such profiles for January, May, and September 2005 are presented in Fig. 5.One can see that the indicated uncertainties are now small enough to make clear conclusions about the extent to which the observed and retrieved profiles agree by comparing their averaged values only, i.e., < OH MLS >, < HO MLS > as well as similar characteristics for OH and O 3 concentration profiles for three months of the year 2005.First, clearly, the distributions of < OH ret > and < O ret 3 > are in good qualitative and quantitative agreement with the initial MLS/Aura measurement data at lower altitudes, below ∼ 0.07 and 0.1 mbar, correspondingly.At higher altitudes, the distributions of < OH ret > reproduce all the main structural features of < OH MLS >, but the retrieved OH concentration has lower values than the observed one with a relative difference OH reaching ∼ 15 % at the top.The distribution of < O ret 3 > above 0.1 mbar, in turn, differs considerably from < O MLS 3 >, both in quantity and quality, and O 3 locally reaches 50-60 % and more.Second, for all months there are significant qualitative and quantitative differences between < HO ret 2 > and < HO MLS 2 >, the most noticeable one being location of the mesospheric maximum of this component's concentration.According to the observations it is close to 0.1 mbar, while the retrieved data demonstrate the altitudes of about ∼ 0.046 mbar or higher.Our analysis of the applied method of statistical evaluation demonstrates that the higher position of this maximum in the distributions of < HO ret 2 > is influenced by the OH MLS data in which the mesospheric maximum (see Figs. 6-8) is also located notably higher than 0.1 mbar.

Discussion and conclusion
On the basis of the data presented in Sect. 5 we can conclude that, upon the whole, simultaneous OH, HO 2 and O 3 satellite measurements poorly satisfy the photochemical equilibrium condition.This bias is most prominent for the HO 2 component.We can conjecture that a possible explanation for the bias is the significant systematic error in HO 2 measurements, particularly at the height of the mesospheric maximum.This assumption is supported by the calculation of the HO 2 distributions with the use of our 3-D chemical transport model (see Fig. 9).It can be seen that the mesospheric maximum of HO 2 in these months, as well as of the < HO ret 2 > distributions lies above 0.046 mbar.
Moreover, new data on the HO 2 distributions were recently obtained from the MLS measurements.Millán et al. (2015) performed the offline retrieval of daily zonal means of HO 2 profiles using averaged MLS radiances measured in 10 • latitude bins.Averaged spectra have a better signalto-noise ratio, which removes many of the limitations of the MLS standard product for HO 2 .In particular, the upper boundary of the altitude region in which daytime data are suitable for scientific use has reached 0.0032 mbar, and the "day-minus-night" correction is not needed at altitudes above 1 mbar.Comparison with various experimental and model data has shown that the offline retrieval reproduces the basic properties of the HO 2 distribution in the mesosphere relatively well (at least qualitatively; Millán et al., 2015).
The offline retrieval product, the alternative data set of daytime HO 2 , has recently become publicly available at https://mls.jpl.nasa.gov(last access: 26 May 2018).Figure 10 shows the monthly averaged zonal means of the offline retrieval data (< HO MLS 2 offline >) and relative differences with retrieved and MLS standard product data (< HO MLS 2 > − < HO MLS 2 offline >)/ < HO MLS 2 offline > and (< HO ret 2 > − < HO MLS 2 offline >)/ < HO MLS 2 offline >, correspondingly.Figure 10 represents the same time periods as Figs.6-8.It is worth noting that the distributions < HO MLS 2 offline > depicted in Fig. 10 represent significantly different amounts of data.The data sets for May and September include 31 and 27 days of measurements, respectively, whereas the January data set encompasses only 4 days.The latter makes the graphs in the first row in Fig. 10 noisier than the others.One can see that the results of the offline HO 2 retrieval show the same features as the results of our evaluation technique in comparison to the standard MLS retrieval, i.e., the height of mesospheric HO 2 maximum is notably higher.We can conclude that the distributions of < HO ret 2 > better match < HO MLS 2 offline > than < HO MLS 2 >, although some quantitative discrepancy between < HO ret 2 > and < HO MLS 2 offline > also exists.Note that this may be due to systematic errors in the HO MLS 2 distributions, which cannot be excluded within the framework of the introduced technique.For a detailed qualitative and quantitative comparison of < HO ret 2 > and < HO MLS 2 offline >, one should modify the method so that a statistical evaluation of the OH MLS and O MLS 3 standard products, and the data of the offline HO 2 retrieval, could be conducted within the framework of a single procedure with no account for the HO MLS 2 distributions.This modification is under way and will be presented elsewhere.
The proposed method for statistical evaluation of mesospheric species measurements can be readily generalized to other atmospheric photochemical systems that contain shortlived components (see Introduction).It may also be modified for assessing hard-to-measure chemical components, characteristics of atmospheric processes (like wind speed or turbulent diffusion rate), or poorly known reaction rates.

Figure 3 .
Figure 3. Daytime monthly averaged zonal mean P H 2 O OH /P OH distributions.

Figure 5 .
Figure 5. Examples of monthly averaged zonal mean vertical profiles of OH, HO 2 , and O 3 measured (red curves) in January, May, and March 2005 and corresponding retrieved profiles (black curves).Solid curves: boundaries of the 65 % confidence intervals, dashed curves: medians.

Figure 6 .
Figure 6.Daytime monthly averaged zonal mean retrieved (a) and measured (b) distributions of HO 2 , OH, and O 3 and their relative difference (c) in January 2005.

Figure 7 .
Figure 7. Daytime monthly averaged zonal mean retrieved (a) and measured (b) distributions of HO 2 , OH, and O 3 and their relative difference (c) for May 2005.

Figure 8 .
Figure 8. Daytime monthly averaged zonal mean retrieved (a) and measured (b) distributions of HO 2 , OH, and O 3 and their relative difference (c) for September 2005.

Figure 9 .
Figure 9. Daytime monthly averaged zonal mean model distributions of HO 2 for January, May, and September.