Journal topic
Atmos. Meas. Tech., 13, 3033–3042, 2020
https://doi.org/10.5194/amt-13-3033-2020
Atmos. Meas. Tech., 13, 3033–3042, 2020
https://doi.org/10.5194/amt-13-3033-2020

Research article 10 Jun 2020

Research article | 10 Jun 2020

A comparison of OH nightglow volume emission rates as measured by SCIAMACHY and SABER

A comparison of OH nightglow volume emission rates as measured by SCIAMACHY and SABER
Yajun Zhu1,2, Martin Kaufmann1,3, Qiuyu Chen1,3, Jiyao Xu2,4, Qiucheng Gong1,3, Jilin Liu1,5, Daikang Wei1,3, and Martin Riese1,3 Yajun Zhu et al.
• 1Institute of Energy and Climate Research, Forschungszentrum Jülich, Jülich, Germany
• 2State Key Laboratory of Space Weather, National Space Science Center, Chinese Academy of Sciences, Beijing, China
• 3Institute for Atmospheric and Environmental Research, University of Wuppertal, Germany
• 4School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing, China
• 5Qian Xuesen Laboratory of Space Technology, China Academy of Space Technology, Beijing, China

Correspondence: Qiuyu Chen (q.chen@fz-juelich.de)

Abstract

Hydroxyl (OH) short-wave infrared emissions arising from OH(4-2, 5-2, 8-5, 9-6) as measured by channel 6 of the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) are used to derive concentrations of OH(v=4, 5, 8, and 9) between 80 and 96 km. Retrieved concentrations are used to simulate OH(5-3, 4-2) integrated radiances at 1.6 µm and OH(9-7, 8-6) at 2.0 µm as measured by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument, which are not fully covered by the spectral range of SCIAMACHY measurements. On average, SABER “unfiltered” data are on the order of 40 % at 1.6 µm and 20 % at 2.0 µm larger than the simulations using SCIAMACHY data. “Unfiltered” SABER data are a product, which accounts for the shape, width, and transmission of the instrument's broadband filters, which do not cover the full ro-vibrational bands of the corresponding OH transitions. It is found that the discrepancy between SCIAMACHY and SABER data can be reduced by up to 50 %, if the filtering process is carried out manually using published SABER interference filter characteristics and the latest Einstein coefficients from the HITRAN database. Remaining differences are discussed with regard to model parameter uncertainties and radiometric calibration.

1 Introduction

Hydroxyl (OH) airglow stems from spontaneous emissions of metastable excited OH molecules which are mainly produced by the exothermic reaction of H and O3 in the upper mesosphere and lower thermosphere (UMLT). Its emission layer peaks at an altitude of approx. 87 km and extends about 8 km . OH airglow covers a broad spectral region from the ultraviolet to near-infrared spectral range and is of importance for studying photochemistry and dynamics in the UMLT region.

Since first confirmed by , OH airglow emissions have been widely observed using various remote spectroscopic techniques . The measurements obtained in such studies have been analyzed for various purposes. For example, rotational temperature can be obtained from OH emissions as a proxy for kinetic temperature under the assumption of rotational local thermodynamic equilibrium (LTE) . Gravity waves passing through the OH airglow layer can be monitored to study the dynamics and energy balance in the UMLT . The understanding of OH relaxation mechanisms with different species can be improved by studying different OH band emissions in the UMLT . Another important application of OH airglow is to derive trace constituents in the UMLT, such as H and O abundances .

OH nightglow has been globally measured, among others, by SABER (Sounding of the Atmosphere using Broadband Emission Radiometry) operating since 2002 and SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) observing from 2002 to 2012. SABER performed observations successfully over a 17-year period, covering one and half solar cycles, and is still measuring, and many outstanding achievements have been accomplished using these data . The OH data obtained by SABER have been used by different investigators to derive atomic oxygen abundance in the UMLT; however, deviations of up to 60 % were found in comparison with atomic oxygen data derived from O(1S) green-line measurements obtained by SCIAMACHY and WINDII (Wind Imaging Interferometer) . This large deviation promoted a discussion on the absolute values of atomic oxygen abundance . derived new atomic oxygen data from SABER OH 2.0 µm absolute radiance measurements in the UMLT under the constraints of the global annual mean energy budget. also retrieved atomic oxygen data from SABER OH 1.6 and 2.0 µm radiance ratios as an alternative approach. Further new atomic oxygen data were recently derived by from SCIAMACHY nighttime OH(9-6) band measurements using rate constants measured in the laboratory by , which agree with atomic oxygen data derived from SCIAMACHY O(1S) green-line and O2 A-band measurements within a range of 10 %–20 % . While the agreement between new atomic oxygen data obtained by SABER and SCIAMACHY has improved, systematic deviations of up to 50 % still persist . This systematic difference needs to be addressed in future studies.

In this study, OH nightglow limb spectra measured by SCIAMACHY were used to derive OH spectrally averaged radiances at 1.6 and 2.0 µm as measured by SABER. The obtained radiances were compared to SABER OH radiometric measurements to investigate whether systematic differences exist between the two datasets.

2 OH nightglow measurements and auxiliary data

From 2002 to 2012, OH Meinel-band near-infrared emissions were measured simultaneously by SCIAMACHY on Envisat and by SABER on the TIMED (Thermosphere, Ionosphere, Mesosphere Energetics and Dynamics) satellite. The spectral range of both instruments covers several OH emission bands stemming from different vibrational states . The SCIAMACHY instrument on Envisat operated in a sun-synchronous orbit with an Equator crossing local solar time of 10:00 and 22:00. The limb spectra used here were observed by SCIAMACHY in a dedicated mesosphere–thermosphere mode and the limb observational range covered 24 tangent altitudes from 73 to 148 km with a vertical sampling of 3.3 km. SCIAMACHY was a multi-channel grating spectrometer, and its channel 6 measured OH spectra arising from upper vibrational states in the range of 2 to 9 at a spectral resolution of 1.5 nm. The measurement error of SCIAMACHY channel 6 is about 1.2 % . Channel 6 covers a spectral range from 971 to 1773 nm . In this study, only the spectral range of channel 6 up to 1589 nm was used due to the reduced performance of the detector beyond this wavelength . It should be noted that SCIAMACHY channels 7 and 8 covered spectral ranges of 1934–2044 and 2259–2386 nm, respectively, but unfortunately suffered from ice condensation on their detectors .

SABER is a multi-channel radiometer and observes radiometric OH(9-7, 8-6) ro-vibrational lines with wavelengths around 2.0 µm (channel 8) and OH(5-3, 4-2) band emissions at about 1.6 µm (channel 9) . The altitude range of the observation covers 60–180 km with vertical resolution of approx. 2 km . To the authors' knowledge, there are no publicly available references on the observed accuracy of the SABER OH channels, except for a presentation named “SABER Instrument Performance and Measurement Requirements” published on http://saber.gats-inc.com/overview.php, which is the official source of SABER data products. According to this document, the estimated accuracy of the 1.6 and 2.0 µm channel data is about 3 % at 80–90 km and about 20 % at 90–100 km. Since the SABER instrument is a radiometer, individual OH ro-vibrational emission lines cannot be resolved. Figure 1 shows simulated OH airglow emissions in the spectral range between 1000 and 2400 nm; spectral ranges covered by the instruments are shaded in different colors.

Figure 1Simulated OH airglow emission bands used in this study in a spectral range between 1000 and 2400 nm. Shaded light blue region covers a spectral range observed by SCIAMACHY channel 6; shaded light red and light green regions cover two spectral ranges measured by SABER channel 9 and 8, respectively.

Since the spectral coverage of SCIAMACHY and SABER does not coincide, we can not compare their measurements directly. However, both instruments observed ro-vibrational lines stemming from the same upper vibrational states. This offered us an opportunity to calculate the number densities of the OH upper vibrational states and then simulate the same ro-vibrational emission bands for the purposes of comparison. In our study, OH limb spectra measured by SCIAMACHY at 1078–1100, 1297–1325, 1377–1404, and 1575–1588 nm were used, as shown in Fig. 2. The spectral ranges covered ro-vibrational lines in the OH(5-2), OH(8-5), OH(9-6), and OH(4-2) bands, with low rotational quantum numbers N (N≤3) to reduce the potential uncertainty that can be introduced by overpopulated high-N rotational states ; details are discussed later. From these measurements, number densities of OH(v=4, 5, 8, and 9) are obtained, which are used to simulate corresponding SABER measurements.

For comparison, SABER V2.0 data were used, including OH in-band and “unfiltered” 1.6 and 2.0 µm data . The in-band OH radiance data comprised raw data that did not take into account filter transmission, while the unfiltered OH radiance data consider the interference filter characteristics measured in the laboratory . The unfiltering process depends on the spectral shape of the underlying ro-vibrational distribution of the emission. This shape has to be determined by a model, which depends on the rotational temperature and on the transition probabilities (Einstein coefficients). In addition to using the official in-band and unfiltered data, separate in-band and unfiltered datasets were obtained from the SCIAMACHY measurements, using the bandpass filter transmission of and various Einstein coefficient datasets; for details see Sect. 4. In this procedure, we also considered OH(3-1) and OH(7-5) emission lines observed by SABER 1.6 and 2.0 µm channels, respectively, and their contributions to the two channels were calculated based on SCIAMACHY OH(3-1) and OH(7-4) measurements. In order to enhance the signal-to-noise ratio and to obtain a large number of coincident measurements with both instruments, monthly zonal median data in 5 latitude bins were used. Since the SCIAMACHY instrument can not measure nighttime temperature in the UMLT, co-located SABER measurements were also used here. The coincidence criteria selected were $±{\mathrm{2.5}}^{\circ }$ in latitude and 1 h in local time.

Figure 2Monthly zonal median OH(4-2), OH(5-2), OH(8-5), and OH(9-6) limb spectra at tangent altitude about 86 km for September 2005 in a latitude range 35–40 N. raw: the raw limb spectra measured by SCIAMACHY; fit: simulated limb spectra as measured by SCIAMACHY from retrieval results.

3 Methodology

3.1 OH emission model

The exothermic reaction of H and O3 in the atmosphere was identified by as the major source of vibrationally excited hydroxyl radicals (OH*) near the mesopause region.

$\begin{array}{}\text{(R1)}& \mathrm{H}+{\mathrm{O}}_{\mathrm{3}}\to {\mathrm{OH}}^{*}\left(v\le \mathrm{9}\right)+{\mathrm{O}}_{\mathrm{2}}\end{array}$

Metastable excited OH* can be de-excited via radiative, chemical, and collisional relaxation processes. OH(v≤8) is not only initially populated by the reaction of H+O3 but is also produced by the deactivation of higher vibrational states of OH* via radiative relaxation and quenching. The number density of OH(v) can be obtained from its emission line measurements by dividing by the corresponding Einstein coefficients. The volume emission rate ${V}_{v-{v}^{\prime }}\left(i\right)$ of an arbitrary ro-vibrational line within a vibrational band OH($v-{v}^{\prime }$) can be calculated as

$\begin{array}{}\text{(1)}& {V}_{v-{v}^{\prime }}\left(i\right)={n}_{v}\cdot \frac{{g}_{v}\left(i\right)\cdot {\mathrm{e}}^{-{E}_{v}\left(i\right)/\left(k\cdot T\right)}}{{Q}_{v}\left(T\right)}\cdot {A}_{v-{v}^{\prime }}\left(i\right).\end{array}$

nv is the total number density of the corresponding upper vibrational state v, and ${A}_{v-{v}^{\prime }}\left(i\right)$ is the Einstein coefficient of the specific state-to-state transition from vibrational level v to v. Ev(i) and gv(i) are the rotational energy and degeneracy of the upper rotational state of the ith line considered. k is the Boltzmann constant and T is the temperature. Qv(T) is the rotational partition sum of OH(v). This formula is only valid under rotational local thermodynamic equilibrium (LTE) conditions and deviations are discussed later.

3.2 Retrieval model

SCIAMACHY measured integrated OH spectra along the line of sight in the tangent altitude range from approx. 73 km to approx. 149 km. The SCIAMACHY OH limb measurements can be expressed as

$\begin{array}{}\text{(2)}& \mathbit{y}=\mathbf{F}\left(\mathbit{x},\mathbit{b}\right)+\mathbit{ϵ}.\end{array}$

y corresponds to the measured SCIAMACHY OH limb spectra. F is the functional formula of the forward model involving Eq. (1). x represents the number densities (nv) of the corresponding upper vibrational state of the emission lines. b is the parameter vector of the forward model, e.g., Einstein coefficients. ϵ represents stochastic measurement errors. The retrieval can be regarded as an approach to solving an inversion problem in the presence of indirect measurements of the properties of interest. In our setup we assume that each atmospheric layer emits OH airglow homogeneously, and we set the retrieval grid to be identical to the tangent altitude grid of the averaged OH limb measurements. To improve the efficiency of the retrieval, to suppress noise in the solution, and to achieve a smooth transition of the retrieved quantities into model data at the upper boundary, a regularization term is added to the minimization (Rodgers2000):

$\begin{array}{ll}{\mathbit{x}}_{i+\mathrm{1}}=& \phantom{\rule{0.25em}{0ex}}{\mathbit{x}}_{\mathrm{a}}+\left({\mathbf{K}}_{i}^{T}{\mathbf{S}}_{\mathit{ϵ}}^{-\mathrm{1}}{\mathbf{K}}_{i}+{\mathbf{S}}_{\mathrm{a}}^{-\mathrm{1}}{\right)}^{-\mathrm{1}}{\mathbf{K}}_{i}^{T}{\mathbf{S}}_{\mathit{ϵ}}^{-\mathrm{1}}\\ \text{(3)}& & ×\left[\mathbit{y}-\mathbf{F}\left({\mathbit{x}}_{i},\mathbit{b}\right)+{\mathbf{K}}_{i}\left({\mathbit{x}}_{i}-{\mathbit{x}}_{\mathrm{a}}\right)\right].\end{array}$

xi reaches the optimal estimate solution when the retrieval converged. Ki corresponds to the first derivative matrix of the forward model, named the Jacobian matrix. xa represents the a priori knowledge of the total number densities of OH(v), and ${\mathbf{S}}_{\mathrm{a}}^{-\mathrm{1}}$ is the regularization matrix. Sϵ is a diagonal error covariance matrix of y.

4 Results and discussion

4.1 Error analysis

The confidence level of simulated volume emission rates (VERs) can be assessed by considering three main aspects: the uncertainty of the auxiliary atmospheric quantities, i.e., temperature; the uncertainty of rate constants, i.e., Einstein coefficients; and the potential uncertainty introduced by overpopulated high rotational states due to nonlocal thermodynamic equilibrium (non-LTE) effects. The temperature uncertainty in the SABER measurements includes random and systematic errors; summarized SABER temperature uncertainties. We consider only systematic errors in SABER temperatures, because monthly mean data are used. The SABER systematic temperature uncertainty is approx. 1.5 K at 70–80 km, 4 K at 90 km, and 5 K at 100 km. Accordingly, VERs are affected by temperature uncertainties by less than 1 % between 80 and 96 km on average, as obtained by in their investigation of the temperature dependence of the band Einstein coefficients as well.

Figure 3Simulated SABER unfiltered OH 1.6 and 2.0 µm volume emission rates from SCIAMACHY data using Einstein A values of HITRAN, , and at 20–40 N for October 2007 (a); Corresponding retrieved OH number densities of vibrational states 9, 8, 5, and 4 from SCIAMACHY data using HITRAN database (b).

Figure 4SABER 1.6 and 2.0 µm unfiltered (a, b) and in-band (c, d) volume emission rates and corresponding simulations from SCIAMACHY data for September 2005 at latitude bins 0–20 N (a, c) and 20–40 N (b, d) The horizontal lines represent error bars considering the total uncertainties discussed in Sect. 4.1. The grey shaded area represents the observed accuracy of SABER 1.6 and 2.0 µm channels.

Many OH Einstein coefficient datasets can be found in the OH research community (see, e.g., ). We consider the values given in the latest HITRAN molecular spectroscopic database and the OH Einstein A values calculated by and . The uncertainty of the Einstein coefficient affects simulated VERs in two ways. In the retrieval of vibrationally excited OH from SCIAMACHY data and in the simulation of the SABER measurements. Figure 3a shows simulated SABER unfiltered OH 1.6 and 2.0 µm VER profiles obtained from SCIAMACHY measurements using these three Einstein coefficient datasets. Corresponding OH number density profiles as derived using Einstein coefficients from the HITRAN database at vibrational states 9, 8, 5, and 4 are also given. The highest VERs are obtained by using the Einstein coefficients calculated by and the lowest are obtained by using the HITRAN database. The differences between them are approx. 26 % for the simulation of SABER 2.0 µm VERs and approx. 19 % for the simulation of SABER 1.6 µm VERs. Similar values were also obtained if we used data for other latitude bins or time periods. The same procedure was also applied to the simulation of SABER 2.0 and 1.6 µm in-band data, giving similar results. Therefore, we used these results as a proxy to estimate related uncertainties of the Einstein coefficients.

Figure 5Latitude–altitude cross sections of monthly zonal median SABER unfiltered 2.0 µm volume emission rates (b) and corresponding simulations (a) from SCIAMACHY data for the year 2007. The numbers represent the month of the year.

Figure 6Scatter plots of SABER 1.6 µm (a, c) and 2.0 µm (b, d) volume emission rates vs. the corresponding simulations using SCIAMACHY data for the year 2007. The color bar shows the latitude. The plus marker indicates the data at 80–85 km, the x marker represents the data at 86–90 km, and the point marker shows the data at 91–95 km. The solid line shows the linear fit to the data. a and b represent the slope and y intercept of the fitting line, respectively. r represents the correlation coefficient of the fitting.

, , and reported that middle and higher excited rotational states (N≥4) of OH do not meet the LTE hypothesis and that these levels are overpopulated. An estimation of the non-LTE contribution was performed by based on cross-dispersed cryogenic spectrometer measurements in the spectral range of 0.97 to 2.4 µm. A combination of two Boltzmann distribution equations with cold and hot OH rotational temperatures was used to predict the observed intensities of OH emission lines. re-analyzed the data used by to estimate the OH rotational temperatures following the approach taken by and . They found that the thermalization of every OH vibrational level is incomplete.

The low spectral resolution of SCIAMACHY spectra does not allow us to estimate this effect from the measured data. Therefore, we performed model simulations using the same approach and parameter sets as to quantify the effect of incomplete thermalization on the spectral ranges used in this study. We calculated OH 1.6 and 2.0 µm VERs by considering only the cold rotational temperature and then obtained them using cold and hot temperatures together as and did. It was found that differences between them are less than 2 % for both SABER channels. The SABER 1.6 and 2.0 µm channels also observe emission lines from OH(3-1) and OH(7-5), respectively. We estimated their influence on spectrally integrated radiances by the derivation of the corresponding emissions using SCIAMACHY OH(3-1) and OH(7-4) nightglow measurements. These simulations show that the contributions of OH(7-5) and OH(3-1) to the two channels are about 3 % and 1 % on average, respectively.

In summary, the uncertainty of the Einstein coefficient dominates the error budget for the in-band and unfiltered data, which is on the order of 20 % and 26 % for the SABER 1.6 and 2.0 µm VER simulations, respectively.

4.2 Comparison of SABER measurements and simulations

An intercomparison between 1.6 and 2.0 µm in-band and unfiltered VERs as measured by SABER and corresponding simulations using SCIAMACHY data and HITRAN OH Einstein coefficients is given in Fig. 4 for two different latitudes in September 2005. Error bars shown in Fig. 4 represent the root-mean-square value of all uncertainties discussed in the Sect. 4.1. Panels (a) and (b) show a comparison of the unfiltered data, and panels (c) and (d) show the in-band data. SABER measurements are always larger than the simulations using SCIAMACHY data. For the unfiltered data, deviations of SABER OH 1.6 µm measurements with respect to the corresponding simulations increase with altitude from 30 % to 45 % at 83 km to 55 %–80 % at 96 km, depending on latitudes. The difference of SABER OH 2.0 µm measurements with respect to the corresponding simulations is 16 % at 86 km. At 96 km, it reaches 70 % in latitude bins 0–20 N and approx. 90 % in 20–40 N.

Surprisingly, for the in-band data, the differences for the 1.6 and 2.0 µm channels are significantly smaller at most altitudes. They vary in a range of 8 %–28 % (21 %–50 %) and 8 %–60 % (28 %–100 %) from 83 to 96 km at 0–20 N (20–40 N). It should be noted that SCIAMACHY and SABER have a resolution of about 3.3 and 2 km, respectively. A linear interpolation has been applied to SABER data to make a comparison with SCIAMACHY data. This may underestimate the SABER data at peak altitudes and overestimate the SABER data at two wings besides the peak altitudes.

Figure 5 shows the global spatial distributions of SABER OH 2.0 µm VERs (panel b) and the corresponding simulations (panel a) using SCIAMACHY data for the year 2007. A strong annual variation with a maximum in April and a semiannual oscillation are visible in the radiance data over the Equator region, as was also found by in a study of SCIAMACHY OH(3-1) and OH(6-2) volume emission rates. It is obvious that SABER VERs are significantly larger than corresponding simulated values based on SCIAMACHY observations, as already stated. Comparing the SABER OH 1.6 µm VERs and the corresponding simulations leads to the same conclusion (not shown).

Figure 6 shows two pairs of scatter plots which elucidate the consistency of SABER unfiltered 1.6 µm (Fig. 6a, c) and 2.0 µm (Fig. 6b, d) VERs and corresponding simulated values based on SCIAMACHY observations. Again, SABER data are systematically larger than the SCIAMACHY simulations. The SABER 1.6 µm channel data (Fig. 6a, c) are 44 % larger for the unfiltered data and 23 % larger for the in-band data, if all altitudes and latitudes are considered simultaneously in one fit. For the 2.0 µm data (Fig. 6b, d), the differences are 23 % and 35 % on average.

To illustrate whether this difference changes on long timescales, Fig. 7 shows the ratio of SABER unfiltered and in-band data to the corresponding simulations based on SCIAMACHY data from 2003 to 2011. For the OH 1.6 µm unfiltered (in-band) data, the ratio value varies roughly between 1.2 (1.0) and 1.3 (1.2) for 2003–2009, reaching 1.1 for 2010 and 1.36 for 2011. The ratio varies between 1.0 (1.1) and 1.1 (1.2) for the OH 2.0 µm unfiltered (in-band) data. The data indicate that there are no significant variations in the slope of SABER data vs. SCIAMACHY simulations from 2003 to 2011 and that there is a systematic bias between them in general.

Figure 7Slope of SABER 2.0 and 1.6 µm volume emission rates vs. the corresponding simulations using SCIAMACHY data from 2003 to 2011.

5 Conclusions

Near-infrared OH nightglow emissions measured by SCIAMACHY channel 6 were used in this study to simulate SABER 1.6 and 2.0 µm radiance measurements to assess systematic differences between the two measurements. Two different SABER data products are used for this comparison: so-called in-band data, which are the data directly obtained from the measurements, and unfiltered data. For the latter, the shape, width, and transmission of the instrument's broadband filters have been considered, and the fraction of OH lines passing the interference filter has been “upscaled” to obtain total band intensities of the corresponding vibrational transitions . If, however, in-band data are used, the data user has to apply the broadband filter transmission curve to the model data themselves. This procedure has decisive advantages, because no a priori assumptions have to be made to upscale partial measurements of OH vibrational bands to total band intensities. This allows us to use consistent datasets of Einstein coefficients in all processing steps.

When SABER OH in-band data are compared to model simulations using SCIAMACHY data, the typical differences are 35 % for 2.0 µm and 23 % for 1.6 µm radiances, whereas the differences are 23 % and 44 % for the unfiltered data, respectively. The significance or uncertainty of these differences is affected by uncertainties in the Einstein coefficients used to “map” SCIAMACHY to SABER data. For the in-band and unfiltered data, this uncertainty is estimated to be about 20 % for the OH 1.6 µm channel and 26 % for the OH 2.0 µm channel. Considering the radiometric uncertainty of both instruments, which is estimated to be about 1 % for SCIAMACHY and 3 %–20 % for SABER, OH 2.0 µm in-band and unfiltered data agree within their combined uncertainties; OH 1.6 µm in-band data also agree remarkably well, but not for the unfiltered 1.6 µm data.

The OH 2.0 µm data measured by SABER and O(1S) green-line emission and OH(9-6) nightglow observed by SCIAMACHY were used in the past to obtain atomic oxygen abundances. Significant differences in atomic oxygen absolute values were reported . These differences are of similar magnitude as uncertainties in the Einstein coefficients and other model parameters used in the retrieval of those data.

Data availability
Data availability.

The SCIAMACHY Level 1b version 8 data used in this study are available at ftp://scia-ftp-ds.eo.esa.int (last access: 10 May 2019). SABER version 2.0 data are available at http://saber.gats-inc.com (last access: 2 July 2019). Derived OH volume emission rate data are available on request.

Author contributions
Author contributions.

YZ initiated the topic, processed the data, performed the analysis, and drafted the manuscript. MK and QC provided insight and instructions and discussed the results regularly. All authors contributed to the revision and improvement of the paper.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

Qiuyu Chen, Qiucheng Gong, Jilin Liu, and Daikang Wei were supported in their work by the China Scholarship Council. Martin Kaufmann was supported by Forschungszentrum Jülich. This research was also supported by the National Science Foundation of China (41831073 and 41674152) and in part by the Specialized Research Fund for State Key Laboratories.

Financial support
Financial support.

The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.

Review statement
Review statement.

This paper was edited by William Ward and reviewed by Konstantinos Kalogerakis, Christian von Savigny, and one anonymous referee.

References

Baker, D. J. and Stair, A. T. J.: Rocket measurements of the altitude distributions of the hydroxyl airglow, Phys. Scripta, 37, 611–622, 1988. a

Baker, D. J., Thurgood, B. K., Harrison, W. K., Mlynczak, M. G., and Russell,  J. M.: Equatorial enhancement of the nighttime OH mesospheric infrared airglow, Phys. Scripta., 75, 615–619, 2007. a

Bates, D. R. and Nicolet, M.: The photochemistry of atmospheric water vapor, J. Geophys. Res., 55, 301–327, https://doi.org/10.1029/JZ055i003p00301, 1950. a

Brooke, J. S., Bernath, P. F., Western, C. M., Sneden, C., Afşar, M., Li, G., and Gordon, I. E.: Line strengths of rovibrational and rotational transitions in the X2Π ground state of OH, J. Quant. Spectrosc. Ra., 168, 142–157, https://doi.org/10.1016/j.jqsrt.2015.07.021, 2016. a, b

Cosby, P. C. and Slanger, T. G.: OH spectroscopy and chemistry investigated with astronomical sky spectra, Can. J. Phys., 85, 77–99, https://doi.org/10.1139/p06-088, 2007. a, b, c

Dawkins, E. C. M., Feofilov, A., Rezac, L., Kutepov, A. A., Janches, D., Höffner, J., Chu, X., Lu, X., Mlynczak, M. G., and Russell III, J.: Validation of SABER v2.0 Operational Temperature Data With Ground-Based Lidars in the Mesosphere-Lower Thermosphere Region (75–105 km), J. Geophys. Res.-Atmos., 123, 9916–9934, https://doi.org/10.1029/2018JD028742, 2018. a

Gao, H., Xu, J., and Wu, Q.: Seasonal and QBO variations in the OH nightglow emission observed by TIMED/SABER, J. Geophys. Res.-Space, 115, A06313, https://doi.org/10.1029/2009JA014641, 2010. a, b

Gordon, I., Rothman, L., Hill, C., Kochanov, R., Tan, Y., Bernath, P., Birk,  M., Boudon, V., Campargue, A., Chance, K., Drouin, B., Flaud, J.-M., Gamache,  R., Hodges, J., Jacquemart, D., Perevalov, V., Perrin, A., Shine, K., Smith,  M.-A., Tennyson, J., Toon, G., Tran, H., Tyuterev, V., Barbe, A., Császár,  A., Devi, V., Furtenbacher, T., Harrison, J., Hartmann, J.-M., Jolly, A., Johnson, T., Karman, T., Kleiner, I., Kyuberis, A., Loos, J., Lyulin, O., Massie, S., Mikhailenko, S., Moazzen-Ahmadi, N., Müller, H., Naumenko, O., Nikitin, A., Polyansky, O., Rey, M., Rotger, M., Sharpe, S., Sung, K., Starikova, E., Tashkun, S., Auwera, J. V., Wagner, G., Wilzewski, J., Wcisło, P., Yu, S., and Zak, E.: The HITRAN2016 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 203, 3–69, https://doi.org/10.1016/j.jqsrt.2017.06.038, 2017. a

Kalogerakis, K. S., Matsiev, D., Sharma, R. D., and Wintersteiner, P. P.: Resolving the mesospheric nighttime 4.3 µm emission puzzle: Laboratory demonstration of new mechanism for OH(ν) relaxation, Geophys. Res. Lett., 43, 8835–8843, https://doi.org/10.1002/2016GL069645, 2016. a

Kalogerakis, K. S., Matsiev, D., Cosby, P. C., Dodd, J. A., Falcinelli, S., Hedin, J., Kutepov, A. A., Noll, S., Panka, P. A., Romanescu, C., and Thiebaud, J. E.: New insights for mesospheric OH: multi-quantum vibrational relaxation as a driver for non-local thermodynamic equilibrium, Ann. Geophys., 36, 13–24, https://doi.org/10.5194/angeo-36-13-2018, 2018. a, b

Kaufmann, M., Lehmann, C., Hoffmann, L., Funke, B., López-Puertas, M., Savigny, C., and Riese, M.: Chemical heating rates derived from SCIAMACHY vibrationally excited OH limb emission spectra, Adv. Space Res., 41, 1914–1920, https://doi.org/10.1016/j.asr.2007.07.045, 2008. a, b, c

Kaufmann, M., Ern, M., Lehmann, C., and Riese, M.: The Response of Atomic Hydrogen to Solar Radiation Changes, in: Climate and Weather of the Sun-Earth System (CAWSES), edited by: Lübken, F.-J., Springer Atmospheric Sciences, Springer Netherlands, 171–188, https://doi.org/10.1007/978-94-007-4348-9_10, 2013. a

Kaufmann, M., Zhu, Y., Ern, M., and Riese, M.: Global distribution of atomic oxygen in the mesopause region as derived from SCIAMACHY O(1S) green line measurements, Geophys. Res. Lett., 41, 6274–6280, https://doi.org/10.1002/2014GL060574, 2014. a, b

Lichtenberg, G., Kleipool, Q., Krijger, J. M., van Soest, G., van Hees, R., Tilstra, L. G., Acarreta, J. R., Aben, I., Ahlers, B., Bovensmann, H., Chance, K., Gloudemans, A. M. S., Hoogeveen, R. W. M., Jongma, R. T. N., Noël, S., Piters, A., Schrijver, H., Schrijvers, C., Sioris, C. E., Skupin, J., Slijkhuis, S., Stammes, P., and Wuttke, M.: SCIAMACHY Level 1 data: calibration concept and in-flight calibration, Atmos. Chem. Phys., 6, 5347–5367, https://doi.org/10.5194/acp-6-5347-2006, 2006. a, b, c

Liu, W., Xu, J., Smith, A. K., and Yuan, W.: Comparison of rotational temperature derived from ground-based OH airglow observations with TIMED/SABER to evaluate the Einstein coefficients, J. Geophys. Res.-Space, 120, 10069–10082, https://doi.org/10.1002/2015JA021886, 2015. a, b

Meinel, I. A. B.: OH Emission Bands in the Spectrum of the Night Sky, Astrophys. J., 111, 555–564, 1950. a

Mertens, C. J., III, J. M. R., Mlynczak, M. G., She, C.-Y., Schmidlin, F. J., Goldberg, R. A., López-Puertas, M., Wintersteiner, P. P., Picard, R. H., Winick, J. R., and Xu, X.: Kinetic temperature and carbon dioxide from broadband infrared limb emission measurements taken from the TIMED/SABER instrument, Adv. Space Res., 43, 15–27, https://doi.org/10.1016/j.asr.2008.04.017, 2009. a

Mlynczak, M. G., Hunt, L. A., Thomas Marshall, B., Martin-Torres, F. J., Mertens, C. J., Russell, J. M., Remsberg, E. E., López-Puertas, M., Picard,  R., Winick, J., Wintersteiner, P., Thompson, R. E., and Gordley, L. L.: Observations of infrared radiative cooling in the thermosphere on daily to multiyear timescales from the TIMED/SABER instrument, J. Geophys. Res.-Space, 115, a03309, https://doi.org/10.1029/2009JA014713, 2010. a

Mlynczak, M. G., Hunt, L. A., Mast, J. C., Thomas Marshall, B., Russell, J. M., Smith, A. K., Siskind, D. E., Yee, J.-H., Mertens, C. J., Javier Martin-Torres, F., Earl Thompson, R., Drob, D. P., and Gordley, L. L.: Atomic oxygen in the mesosphere and lower thermosphere derived from SABER: Algorithm theoretical basis and measurement uncertainty, J. Geophys. Res.-Atmos., 118, 5724–5735, https://doi.org/10.1002/jgrd.50401, 2013. a, b, c

Mlynczak, M. G., A., H. L., Russell III, J. M., and Marshall, B. T.: Updated SABER Night Atomic Oxygen and Implications for SABER Ozone and Atomic Hydrogen, Geophys. Res. Lett., 45, 5735–5741, 2018. a, b, c, d, e

Noll, S., Kausch, W., Kimeswenger, S., Unterguggenberger, S., and Jones, A. M.: OH populations and temperatures from simultaneous spectroscopic observations of 25 bands, Atmos. Chem. Phys., 15, 3647–3669, https://doi.org/10.5194/acp-15-3647-2015, 2015. a, b

Offermann, D. and Gerndt, R.: Upper mesosphere temperatures from OH*-emissions, Adv. Space Res., 10, 217–221, https://doi.org/10.1016/0273-1177(90)90399-K, 1990. a

Offermann, D., Hoffmann, P., Knieling, P., Koppmann, R., Oberheide, J., and Steinbrecht, W.: Long-term trends and solar cycle variations of mesospheric temperature and dynamics, J. Geophys. Res.-Atmos., 115, d18127, https://doi.org/10.1029/2009JD013363, 2010. a

Oliva, E., Origlia, L., Scuderi, S., Benatti, S., Carleo, I., Lapenna,  E., Mucciarelli, A., Baffa, C., Biliotti, V., Carbonaro, L., Falcini, G., Giani, E., Iuzzolino, M., Massi, F., Sanna, N., Sozzi, M., Tozzi, A., Ghedina, A., Ghinassi, F., Lodi, M., Harutyunyan, A., and Pedani, M.: Lines and continuum sky emission in the near infrared: observational constraints from deep high spectral resolution spectra with GIANO-TNG, Astron. Astrophys., 581, A47, https://doi.org/10.1051/0004-6361/201526291, 2015. a, b, c, d, e, f, g

Panka, P. A., Kutepov, A. A., Rezac, L., Kalogerakis, K. S., Feofilov, A. G., Marsh, D., Janches, D., and Yiğit, E.: Atomic Oxygen Retrieved From the SABER 2.0- and 1.6 µm Radiances Using New First-Principles Nighttime OH(v) Model, Geophys. Res. Lett., 45, 5798–5803, https://doi.org/10.1029/2018GL077677, 2018. a, b, c, d

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding: Theory and Practice, vol. 2 of Series on Atmospheric, Oceanic and Planetary Physics, World Scientific, Singapore, 2000. a

Smith, A. K., Marsh, D. R., Russell, J. M., Mlynczak, M. G., Martin-Torres,  F. J., and Kyrölä, E.: Satellite observations of high nighttime ozone at the equatorial mesopause, J. Geophys. Res., 113, D17312, https://doi.org/10.1029/2008JD010066, 2008. a

Smith, A. K., Marsh, D. R., Mlynczak, M. G., and Mast, J. C.: Temporal variations of atomic oxygen in the upper mesosphere from SABER, J. Geophys. Res., 115, D18309, https://doi.org/10.1029/2009JD013434, 2010. a, b

Teiser, G. and von Savigny, C.: Variability of OH(3-1) and OH(6-2) emission altitude and volume emission rate from 2003 to 2011, J. Atmos. Sol.-Terr. Phy., 161, 28–42, https://doi.org/10.1016/j.jastp.2017.04.010, 2017. a

van der Loo, M. P. J. and Groenenboom, G. C.: Theoretical transition probabilities for the OH Meinel system, J. Chem. Phys., 126, 114314, https://doi.org/10.1063/1.2646859, 2007. a, b, c

van der Loo, M. P. J. and Groenenboom, G. C.: Erratum: “Theoretical transition probabilities for the OH Meinel system” [J. Chem. Phys. 126, 114314 (2007)], J. Chem. Phys., 128, 159902, https://doi.org/10.1063/1.2899016, 2008. a, b, c

von Savigny, C., Eichmann, K.-U., Llewellyn, E. J., Bovensmann, H., Burrows,  J. P., Bittner, M., Höppner, K., Offermann, D., Taylor, M. J., Zhao, Y., Steinbrecht, W., and Winkler, P.: First near-global retrievals of OH rotational temperatures from satellite-based Meinel band emission measurements, Geophys. Res. Lett., 31, l15111, https://doi.org/10.1029/2004GL020410, 2004. a

von Savigny, C., McDade, I. C., Eichmann, K.-U., and Burrows, J. P.: On the dependence of the OH* Meinel emission altitude on vibrational level: SCIAMACHY observations and model simulations, Atmos. Chem. Phys., 12, 8813–8828, https://doi.org/10.5194/acp-12-8813-2012, 2012.  a

Xu, J., Smith, A. K., Yuan, W., Liu, H.-L., Wu, Q., Mlynczak, M. G., and Russell, J. M.: Global structure and long-term variations of zonal mean temperature observed by TIMED/SABER, J. Geophys. Res.-Atmos., 112, D24106, https://doi.org/10.1029/2007JD008546, 2007. a

Xu, J., Gao, H., Smith, A. K., and Zhu, Y.: Using TIMED/SABER nightglow observations to investigate hydroxyl emission mechanisms in the mesopause region, J. Geophys. Res.-Atmos., 117, D02301, https://doi.org/10.1029/2011JD016342, 2012. a, b, c, d, e

Xu, J., Li, Q., Yue, J., Hoffmann, L., Straka, W. C., Wang, C., Liu, M., Yuan,  W., Han, S., Miller, S. D., Sun, L., Liu, X., Liu, W., Yang, J., and Ning,  B.: Concentric gravity waves over northern China observed by an airglow imager network and satellites, J. Geophys. Res.-Atmos., 120, 11058–11078, https://doi.org/10.1002/2015JD023786, 2015. a

Zhu, Y. and Kaufmann, M.: Atomic oxygen abundance retrieved from SCIAMACHY hydroxyl nightglow measurements, Geophys. Res. Lett., 45, 9314–9322, https://doi.org/10.1029/2018GL079259, 2018. a, b, c, d, e

Zhu, Y. and Kaufmann, M.: Consistent Nighttime Atomic Oxygen Concentrations from O2 A-band, O(1S) Green-Line, and OH Airglow Measurements as Performed by SCIAMACHY, Geophys. Res. Lett., 46, 8536–8545, https://doi.org/10.1029/2019GL083550, 2019. a

Zhu, Y., Xu, J., Yuan, W., and Liu, X.: First experiment of spectrometric observation of hydroxyl emission and rotational temperature in the mesopause in China, Sci. China Technol. Sc., 55, 1312–1318, https://doi.org/10.1007/s11431-012-4824-7, 2012. a, b

Zhu, Y., Kaufmann, M., Ern, M., and Riese, M.: Nighttime atomic oxygen in the mesopause region retrieved from SCIAMACHY O(1S) green line measurements and its response to solar cycle variation, J. Geophys. Res.-Space, 120, 9057–9073, https://doi.org/10.1002/2015JA021405, 2015. a

Zoutman, E., Derksen, G., Bokhove, H., and Snel, R.: Error budget for SCIAMACHY calibration, Tech. rep., TPD, 2000. a