OH level populations and accuracies of Einstein-A coefficients from hundreds of measured lines
OH level populations and accuracies of Einstein-A coefficients from hundreds of measured lines
Noll, Stefan;Winkler, Holger;Goussev, Oleg;Proxauf, Bastian
OH level populations and accuracies of Einstein-A coefﬁcients from hundreds of measured lines 1,2 3 2 4 Stefan Noll , Holger Winkler , Oleg Goussev , and Bastian Proxauf Institut für Physik, Universität Augsburg, Augsburg, Germany Deutsches Fernerkundungsdatenzentrum, Deutsches Zentrum für Luft- und Raumfahrt, Weßling-Oberpfaffenhofen, Germany Institut für Umweltphysik, Universität Bremen, Bremen, Germany Max-Planck-Institut für Sonnensystemforschung, Göttingen, Germany Correspondence: S. Noll (email@example.com) Abstract. OH airglow is an important nocturnal emission of the Earth’s mesopause region. As it is chemiluminescent radiation in a thin medium, the population distribution over the various roto-vibrational OH energy levels of the electronic ground state is not in local thermodynamic equilibrium (LTE). In order to better understand these non-LTE effects, we studied hundreds of OH lines in a high-quality mean spectrum based on observations with the high-resolution Ultraviolet and Visual Echelle Spectrograph at Cerro Paranal in Chile. Our derived populations cover vibrational levels between v = 3 and 9, rotational levels up to N = 24, and individual -doublet components when resolved. As the reliability of these results critically depends on the Einstein-A coefﬁcients used, we tested six different sets and found clear systematic errors in all of them, especially for Q-branch lines and individual -doublet components. In order to minimise the deviations in the populations for the same upper level, we used the most promising coefﬁcients from Brooke et al. (2016, JQSRT 168, 142) and further improved them with an empirical correction approach. The resulting rotational level populations show a clear bimodality for each v, which is characterised by a probably fully thermalised cold component and a hot population where the rotational temperature increases between v = 9 and 4 from about 700 to about 7,000 K and the corresponding contribution to the total population at the lowest N decreases by an order of magnitude. The presence of the hot populations causes non-LTE contributions to rotational temperatures at low N , which can be estimated quite robustly based on the two-temperature model. The bimodality is also clearly indicated by the dependence of the populations on changes in the effective emission height of the OH emission layer. The degree of thermalisation decreases with increasing layer height due to a higher fraction of the hot component. Our high-quality population data are promising with respect to a better understanding of the OH thermalisation process. 1 Introduction The nighttime emission of the Earth’s atmosphere in the near-infrared is dominated by hydroxyl (OH) airglow (Meinel, 1950; Rousselot et al., 2000; Hanuschik, 2003; Noll et al., 2012, 2015), which originates in the mesopause region in a layer with a width of about 8 km and a typical peak height of 87 km (Baker and Stair, 1988). The various bright roto-vibrational bands of the OH electronic ground state X represent an important tracer for atmospheric dynamics (especially wave propagation), ambient temperatures, and chemical composition (especially atomic oxygen) at these high altitudes, which are mostly probed arXiv:2005.06603v2 [physics.ao-ph] 15 May 2020 by ground- and satellite-based remote sensing (e.g. Taylor et al., 1997; Beig et al., 2003; von Savigny et al., 2012; Mlynczak et al., 2013; Reisin et al., 2014; Sedlak et al., 2016; Noll et al., 2017). For these applications, it is crucial to understand the physical mechanisms that lead to the observed line emission. In the mesopause region, OH is mostly formed by the reaction of hydrogen and ozone (Bates and Nicolet, 1950; Xu et al., 2012), which excites the electronic ground state up to the ninth vibrational level v (Charters et al., 1971; Llewellyn and Long, 1978; Adler-Golden, 1997). The nascent population distribution over the roto-vibrational levels is far from local thermody- namic equilibrium (LTE). As the subsequent relaxation processes by collisions with other atmospheric species are relatively slow compared to the radiative lifetimes of the excited states (e.g., Adler-Golden, 1997; Xu et al., 2012; Kalogerakis et al., 2018; Noll et al., 2018b), the OH emission bands (which contribute to the vibrational relaxation) reveal strong non-LTE effects. The vibrational level populations can be ﬁtted as a function of energy by an exponentially decreasing (i.e. Boltzmann-like) dis- tribution with a pseudo-temperature of around 10,000 K (Khomich et al., 2008; Noll et al., 2015; Hart, 2019a). Hence, OH bands with upper state vibrational levels v up to the highest nascent state can easily be measured. Moreover, the rotational level populations for the different v reveal high overpopulations for high rotational states N compared to the lowest three or four levels under the assumption of a thermal distribution (Pendleton et al., 1989, 1993; Dodd et al., 1994; Cosby and Slanger, 2007; Oliva et al., 2015; Noll et al., 2018b). The pseudo-temperatures for the high-N populations achieve values up to those found for the v levels (Oliva et al., 2015). The theoretical explanation of these populations especially for low v is still uncertain (Dodd et al., 1994; Kalogerakis et al., 2018; Noll et al., 2018b) as their modelling suffers from limitations in the data sets and uncertain input parameters (especially rate coefﬁcients for collisional transitions). It is usually assumed that the ratios of lines related to the lowest N of a ﬁxed v are sufﬁciently close to LTE for a reliable estimate of the ambient temperature (e.g. Beig et al., 2003). However, this assumption appears to be insufﬁcient at least for the highest v, where deviations of several kelvins were found (Noll et al., 2016, 2018b). In addition, small modiﬁcations in the set of considered levels in terms of N can already signiﬁcantly change the corresponding population temperature (Noll et al., 2015). A successful study of OH level populations requires accurate molecular parameters, i.e. line wavelengths, level energies, and Einstein-A coefﬁcients. In particular, the latter suffer from relatively high uncertainties despite numerous dedicated studies for their calculation (e.g., Mies, 1974; Langhoff et al., 1986; Turnbull and Lowe, 1989; Nelson et al., 1990; Goldman et al., 1998; van der Loo and Groenenboom, 2007; Brooke et al., 2016) and evaluation (e.g., French et al., 2000; Pendleton and Taylor, 2002; Cosby and Slanger, 2007; Liu et al., 2015; Hart, 2019b). Apart from the derivation of absolute OH level populations or densities (Noll et al., 2018b; Hart, 2019b), the quality of these transition probabilities especially affects OH-based temperature estimates (Liu et al., 2015; Noll et al., 2015; Parihar et al., 2017; Hart, 2019b) and abundance retrievals for species like atomic oxygen (Mlynczak et al., 2013; Noll et al., 2018b). The persistent uncertainties in the Einstein-A coefﬁcients are obviously related to the molecular structure of OH and the lack of adequate data for the calculation of the molecular parameters (Nelson et al., 1990; Pendleton and Taylor, 2002; Cosby and Slanger, 2007; van der Loo and Groenenboom, 2007; Brooke et al., 2016). In order to improve our knowledge on OH level populations and Einstein-A coefﬁcients, high-quality measurements of a large number of OH lines and a detailed analysis are required. We could perform such a study based on high-resolution spectroscopic data taken with the Ultraviolet and Visual Echelle Spectrograph (UVES; Dekker et al., 2000) at the Very Large Telescope at Cerro Paranal in Chile (24.6 S, 70.4 W). A mean spectrum of the highest-quality spectra (totalling 536 hours of exposure time) allowed us to investigate 723 lines with upper vibrational levels v between 3 and 9 in the optical and near- infrared regime in detail. In many cases, the small doubling effect due to rotational–electronic perturbations between the ground and excited electronic states (Pendleton and Taylor, 2002) was resolved. In Sect. 2, we describe the UVES data set. Then, we discuss the data analysis involving the calculation of the mean spectrum, the measurement of line intensities, and a check of the line positions (Sect. 3). Section 4 discusses the differences in the derived OH level populations for the Einstein-A coefﬁcients of Mies (1974), Langhoff et al. (1986), Turnbull and Lowe (1989), van der Loo and Groenenboom (2008), Rothman et al. (2013), and Brooke et al. (2016). Moreover, the Brooke et al. (2016) reference is used as the basis for an empirical improvement of the coefﬁcients. The corresponding OH level populations are then investigated in detail (Sect. 5). This involves population ﬁtting, a study of the non-LTE contributions to rotational temperatures, and the investigation of population differences caused by a change in the OH emission altitude. Finally, we draw our conclusions (Sect. 6). 2 Data set This study is based on so-called Phase 3 products of the astronomical echelle spectrograph UVES (Dekker et al., 2000) provided by the European Southern Observatory. Noll et al. (2017) selected about 10,400 archived spectra taken between April 2000 and March 2015, extracted the night-sky emission, and performed a complex ﬂux calibration procedure in order to investigate long- term variations in the mesopause region based on OH emission. The studied spectra comprise the wavelength range between 570 and 1040 nm covered by two set-ups centred on 760 and 860 nm. Depending on the width of the entrance slit, the spectral resolving power varied between 20,000 and 110,000. Hence, these data are well suited for OH level population studies as they allow one to measure numerous resolved emission lines. As the exposure time (between 1 and 125 min) and the contamination of the night-sky emission by the astronomical target (the slit length is only between 8 and 12 ) are also strongly varying, it is important to focus on spectra of sufﬁcient quality, especially if very weak lines are studied. The ﬁnal sample of Noll et al. (2017), who studied relatively bright P-branch lines related to low rotational levels, included 3,113 suitable spectra. We use an even smaller subsample of 2,299 spectra as the basis for this study. It is related to the investigation of the faint K(D ) potassium line at 769.9 nm with a mean intensity of about 1 R (rayleigh) by Noll et al. (2019). For that sample, the selected spectra were carefully checked around the K(D ) line. In order to be able to measure even fainter lines in the entire wavelength regime and to have a homogeneous data set for the calculation of a mean spectrum, we further reduced the sample. 45 spectra of the set-up centred on 860 nm were rejected as they showed severe ﬂaws (wrong continuum levels) below 730 nm. This was not a problem for the potassium study. Moreover, we increased the minimum exposure time from 10 to 45 min and reduced the maximum continuum limit around K(D ) from 1 00 100 to 40 R nm . Finally, we only considered spectra that were taken with the standard slit width of 1 , which corresponds to a resolving power of 42,000. 63% of the sample of Noll et al. (2019) was taken with this slit width. 3 The resulting sample consists of 533 high-quality spectra with a total exposure time of 536 hours at a telescope with a diameter of the primary mirror of 8 m. 3 Analysis 3.1 Mean spectrum In order to calculate the probably best high-resolution airglow mean spectrum in the covered wavelength regime so far, we ﬁrst mapped the 533 selected spectra (Sect. 2) with the ﬂux calibration of Noll et al. (2017) applied to a common wavelength grid from 560 to 1061 nm with a step size of 1 pm that well samples airglow lines, which have a full width at half maximum of about 20 pm close to 800 nm. The mapping is necessary since each UVES spectrum has its own wavelength grid. The set-up positioning appears to have an uncertainty of the order of 1 nm. Moreover, the original step sizes vary from 1.8 to 5.2 pm depending on the central wavelength of the set-up (760 or 860 nm), the chip (two chips with spectra separated by a small gap at the central wavelength), and the pixel binning. Pixel pairs in dispersion direction on the chips were merged for 85% of the sample. The rest of the data are unbinned. Before the mean calculation, the spectra were also scaled to be representative of the zenith by using the van Rhijn correction (van Rhijn, 1921) for a thin layer at an altitude of 90 km. As the zenith angles at the mid-exposure times vary from 3 to 64 , this is a crucial correction with factors between 0.46 and 1.00. These factors do not signiﬁcantly change across the entire OH emission layer, i.e. the choice of the reference altitude is not critical. The mean spectrum was calculated by means of a pixel-dependent -clipping approach in order to avoid the contribution of strong sporadic outliers due to technical issues or the contamination by an astronomical target. As the threshold was set to 10 standard deviations, statistical noise and natural variations of the airglow emission do not cause the rejection of a spectrum at a certain pixel. The ﬁnal number of considered spectra at each wavelength after three iterations of the -clipping is displayed in the lower panel of the upper plot in Fig. 1. The clipping only reduced the numbers by a few spectra. There is a trend towards more rejections at longer wavelengths. The plot also reveals the impact of the combination of the two set-ups centred on 760 and 860 nm with 231 and 302 spectra, respectively. The gap between the spectra of the two chips of each set-up and very narrow gaps between the spectral orders at long wavelengths can also be seen. The rounded edges of the sample-related steps in the histogram reﬂect the variation in the wavelength positioning of a certain set-up. The mean spectrum in the upper plot in Fig. 1 shows 15 OH bands marked by the upper and lower vibrational levels v and 00 0 00 v . Bands with v = v v between 3 and 6 are covered. The band strength strongly increases from OH(5-0) to OH(4-1) as bands with higher v and lower v tend to be stronger in the covered wavelength regime. Each band is split into the three branches R, Q, and P, which are characterised by changes of the rotational quantum number N of 1, 0, and +1. While the R and Q branches at the short-wavelength side and central part of the band are relatively compact, the P branch shows relatively wide spaces between the lines. P-branch lines with high upper rotational quantum numbers N are located at distinctly longer wavelengths than those with low N , i.e. they are found in regions that are dominated by other OH bands. For this reason, high-N P-branch lines of the faint OH(7-1) band can also be detected in the UVES mean spectrum. 4 4-1 UVES 760 and 860 nm mean spectrum 8-4 9-5 7-3 10 3-0 6-2 5-1 9-4 8-3 4-0 7-2 9-3 6-1 8-2 5-0 0.6 0.7 0.8 0.9 1.0 UVES 760 nm: N = 231 UVES 860 nm: N = 302 0.6 0.7 0.8 0.9 1.0 Wa0eleng.h [μm] UVES 760 and 860 nm mean spectrum: OH(8-3) 535 0.730 0.735 0.740 0.730 0.735 0.740 Wa0ele)g.h [μm] Figure 1. UVES mean spectrum in rayleighs per picometre. The full spectral range with labelled OH bands is shown at the top. At the bottom, a narrower wavelength range focusing on the Q and P branch of OH(8-3) is plotted. The lines are labelled. The wavelength-dependent number of spectra involved in the calculation of the mean spectrum is indicated in additional subpanels. This number varies due to the use of two set-ups centred on 760 and 860 nm (with some variation) and the -clipping procedure for the mean calculation. Number of spectra Number of spectra −1 Radiance [R pm ] −1 Radiance [R pm ] Q2(1) Q1(1) Q2(2) Q1(2) Q1(3) P2(1) Q1(4) P1(1) P2(2) P1(2) P2(3) P1(3) P2(4) P1(4) P2(5) P1(5) The lower plot of Fig. 1 shows the narrow wavelength range between 727 and 745 nm to demonstrate the good spectral resolution. The plotted range includes the full Q branch and the P branch up to lines with N = 5 of OH(8-3), a band of intermediate strength. The plot clearly shows the splitting of each rotational state by spin–orbit coupling. The Q and P lines 1 1 2 2 (F = 1) related to the electronic substate X are well separated from the fainter Q and P lines (F = 2) of X . For 3=2 2 2 1=2 00 0 the visible lines, the value of F does not change during the transition, i.e. F = F . The intercombination lines which show a change of F are much fainter and are therefore neglected in this study. The spectral resolving power of 42,000 is sufﬁciently high for seeing doubling. The separation of both components can already be found for Q and P lines with relatively low 1 1 0 0 N . In the lower plot of Fig. 1, the largest separation is visible for Q (N = 4). It amounts to 55 pm (Brooke et al., 2016), i.e. this doublet is fully resolved. Separations of more than 200 pm are measurable for P lines with N 11. The faintest marked doublets Q (2) and Q (4) have intensities between 1 and 2 R. They can easily be measured in the UVES mean 2 1 spectrum, which allows one to also detect lines that are more than one order of magnitude fainter (Sect. 3.2). For a general overview of lines (not only OH) that can be accessed with UVES data, see the catalogue of Cosby et al. (2006). It is based on the night-sky atlas of Hanuschik (2003), which involves UVES observations with a total exposure time of 9 hours in the red and near-infrared wavelength range. 3.2 Line intensities As it is the most comprehensive list of calculated OH lines so far, we used the line wavelengths of Brooke et al. (2016) for the identiﬁcation of lines in the UVES mean spectrum (Sect. 3.1) and the derivation of their intensities. For this purpose, the calculated vacuum wavelengths were converted into air wavelengths by means of the formula of Edlén (1966) for standard air, which works well for the UVES data. The default line integration range was set to a width of about 2 resolution elements of the spectrograph (Sect. 2), which is 40 pm at 860 nm, plus the separation of the two -doublet components. If the latter was wider than the 2 resolution elements, the components were measured independently. For an optimal continuum subtraction, the two continuum points for a linear interpolation across the line were deﬁned manually as this approach can better handle contaminations by nearby emissions and absorptions of other lines than an automatic procedure, which was also tested. The wavelengths of the selected continuum points were also used as limits of the integration range. In particular, range modiﬁ- cations were necessary for lines at very long wavelengths around 1 m, where the spectrograph causes extended line wings. The resulting line intensities representing the zenith (Sect. 3.1) were also corrected for molecular absorption in the lower at- mosphere. The complex procedure involving high-resolution radiative transfer calculations and water vapour measurements in the astronomical target spectra is described by Noll et al. (2017). Further details are given by Noll et al. (2015). For the correction of the measured intensities, the derived line transmission values for the 533 individual spectra were averaged. The resulting mean absorption of the measured doublets was 3% and only 6% of the doublets was attenuated by more than 10%. Hence, the related intensity uncertainty after the correction, which can reduce the absorption by up to an order of magnitude, is negligible for most lines. As illustrated in Fig. 1, the mean spectrum is composed of UVES spectra of two different set-ups centred on 760 and 860 nm. The wavelength shift between both set-ups causes changes in the data properties depending on wavelength. In order 6 to minimise the impact of these changes on the measured line intensities, we investigated and corrected two effects: long-term variations in the OH line intensity and ﬂux calibration errors. The former are important since the two UVES set-ups cover very different parts of the sample-related period from May 2000 to July 2014. Before December 2004, there were only observations with the 860 nm set-up (Noll et al., 2017). On the other hand, spectra of this set-up are only present in the selected sample until May 2010. This results in mean 10.7 cm solar radio ﬂuxes (Tapping, 2013) for an averaging period of 27 days of 102 solar ﬂux units (sfu) for the 760 nm set-up and 140 sfu for the 860 nm set-up. According to Noll et al. (2017) (also based on UVES data), the mean solar cycle effect for v between 5 and 9 is 16:1 1:9 % per 100 sfu. There is no signiﬁcant change with v . We took this mean percentage, the set-up-speciﬁc mean solar radio ﬂuxes, and the wavelength-dependent fraction of 760 and 860 nm spectra to correct the OH line intensities to be representative of the mean solar radio ﬂux of the full sample of 123 sfu. The intensity corrections were up to a few per cent with line-dependent differences characterised by a standard deviation of 2.2%. We did not consider the impact of a possible linear long-term trend as it is not signiﬁcant for the UVES data (Noll et al., 2017). In order to test the ﬂux calibration, we calculated mean spectra for each set-up and also measured line intensities. The latter was performed automatically by using the same wavelengths for the line integration as in the case of the mean spectrum of the full sample. The continuum was measured in narrow intervals (about 0.25 resolution elements, i.e. 5 pm wide at 860 nm) around these limiting positions. The wavelength-dependent intensity ratios for the two set-ups were then used to derive correction factors depending on set-up and chip. Taking the wavelength range between the set-up gaps around 760 and 860 nm as the reference, we found correction factors close to 1 with a standard deviation of 2.7%, which is consistent with the relative ﬂux calibration uncertainty of about 2% for the UVES data set reported by Noll et al. (2017). Combining the solar activity and ﬂux calibration correction, the resulting standard deviation is only 1.7% as both effects partly cancel out. The quality of the ﬁnal line intensities was indicated by a ﬂag consisting of a primary and a secondary classiﬁer. Each classiﬁer is represented by a digit between 0 and 3. A value of 3 corresponds to a reliable measurement of the entire doublet, which requires symmetric line emission and a featureless underlying continuum, 1 and 2 refer to reliable measurements only for the -doublet component with e or f parity in the upper state, and 0 marks uncertainties for both components. For unresolved doublets, only 0 and 3 are possible digits. The introduction of the secondary classiﬁer allows for a ﬁner classiﬁcation scheme. For example, the combined classes 30 and 03 can be used for ambiguous cases. Reasons for measurement uncertainties are obvious or possible blends with other emission lines, regions of signiﬁcant absorptions in the continuum (often combined with very low transmission at the position of the line), and insufﬁcient signal-to-noise ratio in the case of very weak lines. Figure 2 shows a histogram of the measured lines depending on the decadal logarithm of the intensity in rayleighs. In total, 723 doublets are included. This neglects 13 measurements with digit 0 for the primary and secondary classiﬁer. These lines were not further used in this study. Other potential lines could not be measured due to a blend with a stronger line, no detection, or a line wavelength within the order gaps in the near-infrared (Fig. 1). The measured intensities of the 723 doublets range from 0.01 to 600 R, i.e. they comprise almost 5 orders of magnitude. Intensities around 1 R are most abundant. The median is 1.7 R. There is a conspicuous drop in the occurrence frequency below about 0.2 R, which suggests a strongly increasing incompleteness of detections for fainter lines. The intensity of weaker doublets also tends to be more uncertain as the intensity distributions for the different quality classes show. For the classes 3, 1+2, and 0, the median intensities are 2.6, 7 both Λ lines one Λ line ncertain −2 −1 0 1 2 3 log(I [R]) Figure 2. Distribution of the decadal logarithm of the measured OH -doublet intensities in rayleighs. Three categories are indicated: reliable measurement of both -doublet components (detached or unresolved) in green (class 3), reliable measurement of only one component in orange (classes 1 and 2), and uncertain intensities in violet (class 0 without 00). H + O uncertain one line both lines 3 4 5 6 7 8 9 Figure 3. Distribution of the upper level energies E of the measured roto-vibrational transitions in inverse centimetres separated for each upper vibrational level v . The same categories as in Fig. 2 are marked by different colours and line lengths (see legend). The likely exother- micity limit of the hydrogen–ozone reaction is marked by the grey dotted line. 0.83, and 0.086 R. 546 doublets or 76% belong to class 3, where the two components were measured independently in 34% of the cases. There are 122 doublets (17%) with only one reliable component (1+2). Finally, there are 55 cases (8%) with class 0, 75% of them with resolved doublets. Considering that detached components require independent line measurements (350 cases), the total number of measurements for the data in Figure 2 amounts to 1,073. ′ 1 E [cm ] N mber of do blets 0 0 0 The 723 studied doublets probe 236 different upper states characterised by v , N , and F . Up to nine doublets contribute 0 0 to the population data for a certain level. The distribution of level energies E depending on v is shown in Fig. 3. The energies 1 0 0 range from 10,211 to 28,051 cm . Except for v = 3, where N only up to 9 could be measured (mainly due to the wavelength 0 0 limitations of the UVES data), wide ranges of E are covered by the data for the different v . A maximum energy range of 1 0 0 0 8,861 cm is achieved for v = 4. This is possible due to N up to 24. The energy ranges shrink for higher v due to a steeper 0 0 decrease of the line intensities with increasing N , which reduces the detectability of high-N lines. An important reason for this is certainly the closer exothermicity limit of the hydrogen–ozone reaction, which produces the excited OH. Nevertheless, there are nine levels above this limit if we assume 3.38 eV (Cosby and Slanger, 2007), i.e. about 27,260 cm (Noll et al., 2018b). This suggests that the kinetic energy involved in the reaction is also important to populate the OH roto-vibrational 0 0 0 1 levels. For the excitation of the highest level found (v = 9, N = 12, F = 1), about 800 cm of additional energy would be needed. Figure 3 also provides the primary quality classes for the lines related to the displayed states. 79% of the levels are covered by at least one line with class 3. An exclusive class 0 contribution is found for only 16 states. However, excluding 0 0 0 0 uncertain lines can reduce the E range for a given v . In particular, the maximum N for v = 5 shows a decrease from 23 to 20 in this case. 3.3 Line positions The high resolving power of 42,000 of the UVES data used allows for a check of the quality of the input line positions, which were taken from Brooke et al. (2016) and were converted to standard air using the formula of Edlén (1966). The default positioning of the integration ranges for the line intensity measurements described in Sect. 3.2 worked well in most cases. However, signiﬁcant shifts were necessary for some high-N lines. For a systematic study of these offsets, we took the manually adapted integration windows to calculate the intensity-weighted centroid wavelength for each line. In Fig. 4, we show the difference between observed and model wavelengths in picometres as a function of upper state energy (neglecting the vibrational energy) for 406 individually measured -doublet components of the P branch. This selection rejects unresolved or only partly resolved doublets and lines with uncertain central wavelengths due to blending with other lines. 0 1 The remaining 66 Q-branch and 67 R-branch -doublet components are not plotted as they only probe E up to 5,000 cm and are essentially consistent with the P-branch lines, which tend to have higher signal-to-noise ratios. The plot shows for all 0 0 1 v a very good agreement of observed and modelled wavelengths in the case of low energies. For E lower than 2,000 cm , the mean value and standard deviation are 0:4 and 0:7 pm, respectively. The systematic offset is much less than the original pixel size in the UVES spectra (Sect. 3.1). Hence, it can be caused by uncertainties in the wavelength calibration. Moreover, the assumption of standard air conditions (1013 hPa, 288 K, and no H O) for the UVES instrument might cause a part of the offset. Beyond 3,000 cm , the displayed wavelength offsets show an increasing scatter. In part, this is caused by the higher measurement uncertainties for the fainter lines, but there are also clear trends depending on v . In general, the difference 0 0 between the measured and theoretical line wavelengths increases with E . This increase appears to be stronger for higher v . 0 1 0 0 While the change for v = 4 is only about 1 pm at around 8,000 cm , it is about 20 pm for v = 5. For higher v (at least for 6 and 7), the increase of the offsets appears to be even stronger. However, as the covered energy range decreases as well, the 9 Sampl : r sol+ d P-branch Λ doubl ts + = 3 + = 4 20 + = 5 + = 6 + = 7 + = 8 + = 9 0 2000 4000 6000 8000 ′ ′ ′ ′ ′ ′ −1 E (v , N , F ) − E (v , 1, 1) [cm ] Figure 4. Wavelength difference between the UVES-based measurement of isolated P-branch lines (resolved doublets) and the corre- sponding position provided by Brooke et al. (2016) (similar to HITRAN) in picometres as a function of the energy of the upper state of the 0 0 transition E with respect to the lowest energy of the corresponding vibrational level v in inverse centimetres. The coloured markers indicate v . 0 0 maximum offsets only amount to a few picometres. Hence, only the measured shifts for v = 5 and N of 22 and 23 are of the order of a spectral resolution element. For all other detected lines, the quality of the theoretical line positions is much better. The discussed results are for the line wavelengths published by Brooke et al. (2016). As the HITRAN line database (Gordon et al., 2017) is more frequently used, we also calculated the wavelength shifts for those data. We took the version HITRAN2012 (Rothman et al., 2013), which does not differ from the more recent version HITRAN2016 (Gordon et al., 2017) in terms of the OH data. The results are very similar to those in Fig. 4. Strong deviations above 10 pm are found for the same small sample of lines, although OH(5-1)P (23) is missing in the HITRAN database. The mean difference between the line wavelengths from Rothman et al. (2013) and Brooke et al. (2016) is 0.07 pm. The standard deviation only amounts to 0.40 pm. Based on the UVES mean spectrum of Hanuschik (2003) (Sect. 3.1), the accuracy of OH line wavelengths was already investigated by Cosby et al. (2006). Their theoretical line positions originate from Cosby et al. (2000) but should be similar to Goldman et al. (1998), the basis of the OH data in HITRAN, for low rotational levels. For higher rotational levels, the line wavelengths calculated by Cosby et al. (2000) should be more precise. Indeed, although the spectrum of Hanuschik (2003) is noisier, the critical OH(5-1) lines do not show clear systematic offsets. However, all OH data indicate a mean shift of the UVES-based wavelengths of about +1:0 pm. This offset is similar to those for atomic and molecular oxygen in the same wavelength range, which were also measured by Cosby et al. (2006). Thus, systematic errors in the wavelength calibration are the most likely explanation. In comparison, our measurements result in a negative mean offset. This could be caused by differences in the UVES sample, data processing, and analysis. − − − [pm] UVES B+16 4 Einstein-A coefﬁcients 4.1 Full OH level populations 0 00 The intensities I 0 00 derived in Sect. 3.2, where i and i are the upper and lower states of the roto-vibrational transition, can i i be converted into level populations by dividing Einstein coefﬁcients A 0 00 . For a visualisation, these populations are usually i i normalised by dividing the statistical weight (i.e. the degeneracy) of the upper state g = g and then logarithmised. Following Noll et al. (2015), we deﬁne I 0 00 i i y := ln ; (1) 0 00 A g i i where the intensity is given in rayleighs and the Einstein-A coefﬁcients are provided in inverse seconds, which is consistent 6 2 with population column densities in units of 10 cm (Noll et al., 2018b). For doublets, i is characterised by the vibrational level v, rotational level N , and electronic substate F . In the case of individual components, the parity p (i.e. e or f) is also a parameter and g is half as large as for the doublet. Apart from the uncertainties in the line intensities, the quality of the resulting populations also depends on the reliability of the Einstein-A coefﬁcients. As already mentioned in Sect. 1, the latter is not satisfactory as the available sets differ quite signiﬁcantly. With our large sample of energy levels, where the population of each state can be derived from up to nine different lines, we can carry out a comprehensive comparison of Einstein-A coefﬁcients. As a reference set, we take the coefﬁcients calculated by Brooke et al. (2016) (B+16), who provide the most recent and largest set of OH line parameters. The left upper panel of Fig. 5 shows the corresponding y for the 544 reliable doublets of class 3 (neglecting OH(7-1), i.e. two doublets) as a function of the energy of the upper state E . The distribution of populations displays the well-known pattern of steep population 0 0 decreases for low N and weaker population gradients for higher N (Pendleton et al., 1989, 1993; Cosby and Slanger, 2007; Oliva et al., 2015; Kalogerakis et al., 2018; Noll et al., 2018b). Moreover, it indicates the expected decrease of populations for 0 0 higher v with a remarkable exception for v = 8 (Cosby and Slanger, 2007; Noll et al., 2015). The latter is a signature of the nascent OH level population distribution, which mainly occupies v = 8 and 9. The population properties will be discussed in more detail in Sect. 5. It is now important to know how robust the observed pattern is with respect to changes in the set of Einstein-A coefﬁcients. For this purpose, we also consider the HITRAN database with the version from 2012 (Rothman et al., 2013) (see also Sect. 3.3), which is mainly based on the calculations of Goldman et al. (1998) for OH. Moreover, we use the coefﬁcients from van der Loo and Groenenboom (2008) (vdLG08), i.e. the corrected version of van der Loo and Groenenboom (2007), Turnbull and Lowe (1989) (TL89), Langhoff et al. (1986) (LWR86), and Mies (1974) (M74). We neglect the also still popular data from Nelson et al. (1990) as their line list only focuses on low v and low N . Except for B+16, our selection of sets agrees with Liu et al. (2015) and Hart (2019b), who studied the impact of Einstein-A coefﬁcients on the populations of low rotational levels. Other comparisons used a smaller number of sets (French et al., 2000; Cosby and Slanger, 2007; Noll et al., 2015; Parihar et al., 2017; Noll et al., 2018b). Note that the three oldest sets lack a signiﬁcant number of the measured 723 doublets. The set of Turnbull and Lowe (1989) only includes 634 doublets with maximum N between 13 (R branch) and 15 (Q and P 1 2 2 11 4 −2 −4 −6 B+16 HITRAN vdLG08 )2 )4 )6 TL89 LWR86 M≔4 10000 15000 20000 25000 10000 15000 20000 25000 10000 15000 20000 25000 ′ )1 ′ )1 ′ )1 E [cm ] E [cm ] E [cm ] Figure 5. Distribution of logarithmic OH level populations y, i.e. line intensity in rayleighs divided by Einstein-A coefﬁcient in inverse 0 0 seconds and the statistical weight of the upper level g , as a function of the upper level energy E in inverse centimetres. The six panels show the results for the reliable -doublet measurements (see Fig. 2) neglecting OH(7-1) (two lines) for different sets of Einstein-A coefﬁcients: Brooke et al. (2016) (B+16), HITRAN (Rothman et al., 2013), van der Loo and Groenenboom (2008) (vdLG08), Turnbull and Lowe (1989) (TL89), Langhoff et al. (1986) (LWR86), and Mies (1974) (M74). The number of plotted data points is smaller for older sets due to their limitations in the coverage of faint bands and high rotational levels. The populations derived from different branches (P, Q, and R) are highlighted by different symbols and colours (see legend in lower left panel). branches). In the case of LWR86 and M74, the N -related limits are higher by 1 but the bands with v = 6, i.e. essentially OH(8-2) and OH(9-3), are not covered. The number of doublets is therefore only 566 or 78% of the full sample. Figure 5 reveals clear discrepancies between the populations for the six investigated sets of Einstein-A coefﬁcients. The general structure of the distribution is similar but the y values are shifted. Taking 416 doublets with N 12 and v 5, which are present in all six sets, we ﬁnd mean y between 1:43 for TL89 and 0:02 for LWR86 (Table 1). This corresponds to ′ ′ ′ ′′ ′ ′′ ′ ′′ ′ ′′ ( ≔ ln(I /(A g )), i = (v, N, F) y ≔ ln(I /(A g )), i = (v, N, F) i i i i i i i i 00 Table 1. Comparison of populations y and population ratios y with respect to changes in branch and v for six sets of Einstein-A coefﬁ- cients a b c Set hyi hyi hyi hjyji hjyji hjyji hjyji d 00e 00 00 00 N N v v v v Q P R P P Q R B+16 0:18 0:25 0:08 0.11 0.08 0.17 0.10 HITRAN 0:69 0:31 0:28 0.19 0.12 0.27 0.27 vdLG08 0:13 0:25 0:06 0.17 0.13 0.23 0.19 TL89 1:43 0:37 0:51 0.38 0.28 0.44 0.53 LWR86 0:02 0:25 0:02 0.18 0.14 0.25 0.20 M74 0:33 0:27 0:24 0.41 0.42 0.36 0.43 N 416 82 96 127 65 30 32 sel mean logarithmic level population y mean difference of logarithmic level populations mean absolute difference of logarithmic level populations change in branch (Q P and R P) e 00 change in v for all and individual branches (P, Q, R) number of selected doublets (present in all sets of Einstein-A coefﬁcients) an unsatisfactorily large population ratio of about 4.1. Substituting the extreme TL89 y value by the next highest one of 0:69 for HITRAN, the ratio is still about 1.9. The coefﬁcients of M74, B+16, and vdLG08 result in intermediate mean y values of 0:33, 0:18, and 0:13, respectively. Any estimate of absolute OH level populations by means of OH line intensities will be highly uncertain with these results if the quality of the Einstein-A coefﬁcients used cannot be evaluated. Tests of the accuracy of the measured absolute populations require an alternative calculation which is less sensitive to the choice of the set of molecular parameters. Using a kinetic model for chemical OH production and excitation relaxation via collisions and radiative transitions is a solution, although various required rate coefﬁcients and molecular abundances (especially for atomic oxygen) are quite uncertain. Noll et al. (2018b) found that the higher populations related to B+16 tend to be more reliable than the lower ones based on HITRAN; this was for v = 9 based on OH line intensities from UVES and pressure, temperature, and molecular abundance proﬁles from the satellite- based Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument (Russell et al., 1999) and the empirical atmospheric NRLMISE-00 model (Picone et al., 2002). Hence, at least the very low populations related to TL89 appear to be quite unlikely. 13 B+16 HITRAN vdLG08 0.0 −0.5 −1.0 Q−P R−P −1.5 TL89 LWR86 M74 0.0 20.5 21.0 21.5 0 1000 2000 3000 0 1000 2000 3000 0 1000 2000 3000 ′ ′ ′ ′ ′ ′ 21 ′ ′ ′ ′ ′ ′ 21 ′ ′ ′ ′ ′ ′ 21 E (v , N , F ) 2 E (v , 1, 1) [c) E (v , N , F ) 2 E (v , 1, 1) [c) E (v , N , F ) 2 E (v , 1, 1) [c) Figure 6. Difference in logarithmic OH level populations y by the change of the branch of the measured transitions from P to Q or R as a 0 0 0 0 function of the upper state energy E relative to the v -speciﬁc zero point (corresponding to N = 1 and F = 1) in inverse centimetres for six sets of Einstein-A coefﬁcients. For more details, see legend and Fig. 5. 4.2 Detailed population comparisons Figure 5 marks the populations derived from lines of the R, Q, and P branches by different symbols. A good set of Einstein-A coefﬁcients should result in similar population distributions for the three branches. However, the TL89 y values for the R branch are distinctly below those of the P branch. A similar but weaker effect can also be seen for the data related to HITRAN and M74. In order to study these discrepancies in more detail, we calculated differences y for lines with the same upper state but different branches. We focus on Q versus P branch and R versus P branch. The corresponding results for the six sets of Einstein-A coefﬁcients and reliable doublets of class 3 are shown as a function of E relative to the lowest energy for a given v in Fig. 6. For B+16 Einstein-A coefﬁcients, 219 population ratios y are plotted. 10 fo, d(ffe,en. b,anche− 10ΔforΔdifferen.Δbranches ′ ′ Mean for 12 line pairs (N = 2, 3 for const. F ) P+R B+16 HITRAN vdLG08 TL89 LWR86 M74 Set of Einstein-A coefficients Figure 7. Rotational temperature T in kelvins for six sets of Einstein-A coefﬁcients (see Fig. 5). The plotted mean temperatures were rot 0 0 derived from 12 pairs of OH lines with ﬁxed electronic substate F originating from the second and third rotational upper level N . The branch was either P (circles) or R (stars). The diamonds show the resulting T for a combination of the y(E ) regression slopes for both rot branches (mean of inverse T ). rot 0 0 All sets indicate unsatisfactory ratios for the comparison of Q- and P-related populations, especially for high E or N , where y can be lower than 1. The mean y for a subsample with N 12 and v 5 (c.f. Sect. 4.1) are between 0:37 for TL89 and 0:25 for B+16, vdLG08, and LWR86 (Table 1), i.e. the different sets fail in a similar way. According to the theoretical considerations of Pendleton and Taylor (2002) triggered by the OH(6-2) line intensity ratios measured by French et al. (2000), this can be explained by the general negligence of orbital angular momentum uncoupling, which is related 2 2 + to rotational–electronic mixing of the electronic ground state X and the ﬁrst excited state A , for the calculation of the available Einstein-A coefﬁcients. OH line measurements in near-infrared spectroscopic data from the Nordic Optical Telescope at La Palma (Spain) by Franzen et al. (2019) indicate that too low Q-branch populations or too high Einstein-A coefﬁcients (based on HITRAN) are also an issue for OH bands with v = 2 and 3 not covered by our study. The comparison of populations based on R- and P-branch lines reveals a more complex situation than for the Q-branch data. All sets of Einstein-A coefﬁcients show negative mean y, i.e. lower R-related populations on average (Table 1). However, the range is relatively wide with values between 0:51 for TL89 and 0:02 for LWR86. The latter is the only satisfactory set 0 0 for this comparison at low E as it was already found by French et al. (2000) for OH(6-2) low-N lines. For levels with high rotational energy, y tends to be positive. For the other sets, Fig. 6 shows a clear decrease of y with increasing E at least 1 1 below 1,000 cm . For example, the most recent set B+16 shows mean y of 0:04 and 0:14 below and above 400 cm , respectively. At high E , the negative trend appears to vanish for B+16, HITRAN, and vdLG08, the latter even indicating an increase at the highest energies as in the case of LWR86. The especially bad performance of the TL89, HITRAN, and M74 coefﬁcients is probably related to an underestimation of the vibration–rotation interaction (Pendleton and Taylor, 2002). T [K] rot 0 The differences in the dependence of the P- and R-branch-based populations on E for the investigated sets of Einstein-A coefﬁcients imply deviations in the related rotational temperatures T = (2) rot dy B 0 dE (Mies, 1974; Noll et al., 2018b), where k is the Boltzmann constant and dy=dE represents the slope of a regression line in a y(E ) plot like Fig. 5 for the included level populations. For a quantitative T comparison, we considered pairs of levels with rot 0 0 a difference in N of 1 where the populations were derived from reliable lines (class 3) of the same OH band, F , and branch. Only those pairs that are available for the P and R branches were selected. This resulted in 35 pairs that are covered by all sets of Einstein-A coefﬁcients up to N = 12. The sample is relatively small since R-branch lines are often blended. The highest 0 0 number of 12 pairs is found for the combination of the lowest N of 2 and 3 (N = 1 does not exist for the R branch). Mean results for these 12 pairs (which minimise the measurement uncertainties) are shown in Fig. 7. The T differences based on rot higher N agree qualitatively. For the P branch, Fig. 7 reveals a wide range of mean temperatures between 195 K for LWR86 and 208 K for TL89, i.e. the selection of the set of Einstein-A coefﬁcients strongly affects the derivation of absolute T . The situation is better if the rot extremely high value for TL89 is neglected. In this case, the maximum difference (now limited by the HITRAN-related result) is only 6 K instead of 13 K. Moreover, the T for the two most recent sets B+16 and vdLG08 agree well with the minimum rot related to LWR86. The temperature differences are consistent with those derived by Liu et al. (2015) for low-N P -branch lines of the OH bands (3-0), (5-1), (6-2), (8-3), and (9-4) based on observations with a Czerny–Turner spectrometer at Xinglong in China. The differences between the highest and lowest T related to TL89 and LWR86, respectively (B+16 was not published rot yet), were between 9 K for OH(3-0) and 17 K for OH(8-3) with the same mean of 13 K. The trend of decreasing T differences rot for OH bands with longer central wavelengths can also be observed in our data. Taking the differences between HITRAN and B+16 as an example, we ﬁnd between 2 K for OH(3-0)P and 9 K for OH(6-1)P for the 12 selected line combinations. In this 1 2 context, the result of Hart (2019b) for the P branch of OH(4-2) is interesting. Based on data from an astronomical spectrograph at Apache Point in the USA, he found a maximum difference of 3 K for the same ﬁve sets investigated by Liu et al. (2015). If the minimum related to LWR86 is excluded, the variation is only about 1 K with the lowest T related to TL89. rot Figure 7 also shows T based on R-branch lines, which were not used in the discussed studies. The set-dependent results rot are remarkable since they mirror those for the P-branch lines. Now, T ranges from 179 K for TL89 to 193 K for LWR86, i.e. rot the maximum difference of 14 K is very similar to the result for the P branch but the sign is reversed. Moreover, all T related rot to the R branch are lower than those related to the P branch. Hence, the T difference between P and R branch is between 2 K rot for LWR86 and 29 K for TL89. For individual double pairs of lines, R-branch-related T can also be higher than those for the rot P branch, i.e. LWR86 might not show the smallest differences. However, the large discrepancies for TL89 are obvious in any case. As the P- and R-branch T show an oppositional behaviour, we averaged the slopes dy=dE for both branches to derive rot more robust temperatures. As demonstrated by Fig. 7, this was achieved. The mean value for all sets is 193.3 K with a standard deviation of only 1.0 K. The latter represents less than 20% of the variation for the individual branches. Consequently, the 16 B+16 HITRAN .dLG08 0.5 0.0 −0.5 Ref.: (4-0), (5-1), (6-2), (7-3), (8-3), (9-4) −1.0 TL89 LWR86 M74 0.5 0.0 −0.5 −1.0 0.6 0.7 0.8 0.9 1.0 0.6 0.7 0.8 0.9 1.0 0.6 0.7 0.8 0.9 1.0 Wa.e(eng-h [2) Wa.eleng−hΔ[μm Wa.eleng−hΔ[μm Figure 8. Difference in logarithmic OH level populations y by the change of the lower vibrational level v of the measured transitions as a function of the mean wavelength of the doublet in micrometres for six sets of Einstein-A coefﬁcients. For the reference bands OH(4-0), OH(5-1), OH(6-2), OH(7-3), OH(8-3), and OH(9-4), the resulting y are zero and are therefore not shown. OH(3-0) is also neglected since it is the only band with v = 3. Finally, additional bands at short wavelengths are missing in the case of the limited sets of Langhoff et al. (1986) and Mies (1974). For more details, see legend and Fig. 5. combination of P- and R-branch data can signiﬁcantly reduce the impact of the choice of the Einstein-A coefﬁcients on the quality of the resulting T . However, in practice, this will be difﬁcult to apply due to the difﬁculties in measuring R-branch rot lines at moderate spectral resolution. Hence, it is more promising to improve the Einstein-A coefﬁcients by a better handling of the vibration–rotation interaction (Pendleton and Taylor, 2002), which appears to be the main reason for the set-dependent T discrepancies. Data as plotted in Fig. 7 can provide important constraints for this purpose. rot Another population-independent evaluation of Einstein-A coefﬁcients is possible for transitions with the same upper and lower levels except for a different v . For the comparison of the related y, it was necessary to deﬁne a reference OH band for ′′ ′′ 1yΔforΔdifferen−Δ. 1yΔforΔdifferentΔ. 0 each v between 4 and 9, where we have line measurements for two or more bands. We preferentially selected bands with good quality data in the middle of the covered wavelength range: OH(4-0), OH(5-1), OH(6-2), OH(7-3), OH(8-3), and OH(9-4). The resulting y are plotted in Fig. 8 as a function of line wavelength for the six sets of Einstein-A coefﬁcients. In the case of B+16, population ratios for 182 pairs of reliable doublets (class 3) are shown. For LWR86 and M74, this number is only 136 due to the limitations in N and v (Sect. 4.1). The plots indicate a complex behaviour where the y depend on band, branch, 0 0 and N in a different way for each set of Einstein-A coefﬁcients. The data points for the lowest N , which tend to cluster for each band, show a clear trend with wavelength (or v) for all sets except B+16. The data for HITRAN, vdLG08, and TL89 indicate an increase of y with wavelength, whereas the M74 data show a decrease. For LWR86, y is mainly negative, i.e. the reference bands in the middle of the wavelength range with y = 0 (not plotted) indicate the highest relative populations. The overall performance of each set can be evaluated by measuring the mean absolute y for line pairs where Einstein-A coefﬁcients are available in all sets. The corresponding results for 127 line pairs fulﬁlling v 5 and N 12 are provided in Table 1. The highest and hence worsthjyji were found for M74 (0.41) and TL89 (0.38). Lower but still unsatisfactory values of around 0.18 were obtained for HITRAN, vdLG08, and LWR86. B+16 clearly shows the best performance with a value of 0.11. Table 1 also containshjyji depending on branch. The best results are obtained for the P branch for all sets except M74, which is unsatisfactory for all branches. The large y and their trend with wavelength for M74 and TL89 shown in Fig. 8 were already found by Cosby and Slanger (2007). Also using UVES data, they compared the populations derived from the P (1) line of the accessible OH bands with v of 6, 8 and 9. Including the transition probabilities of M74, LWR86, TL89, and Goldman et al. (1998), their analysis favoured the latter, i.e. the main input source for HITRAN. Cosby and Slanger (2007) explained the bad performance of the TL89 coefﬁcients by the erroneous intensity calibration of data used for the applied empirical dipole moment function (DMF), which is the basis for the calculation of the transition probabilities. Population comparisons for OH lines from near-infrared bands with low v mostly not covered by UVES were performed by Oliva et al. (2013) based on observations between 0.95 and 2.4 m with the high-resolution echelle spectrograph GIANO at the Telescopio Nazionale Galileo at the La Palma Observatory in Spain. The results show clear discrepancies between populations derived from lines of bands with v = 2, 3 and 4 for the Einstein-A coefﬁcients from van der Loo and Groenenboom (2007). Interestingly, the corresponding trend with wavelength displayed in Fig. 8 seems to be reversed for bands at longer wavelengths. In general, it can be expected that the accuracy of 0 0 Einstein-A coefﬁcients for bands with high v in the optical tends to be worse than in the case of bands with low v in the near- infrared. Theoretical ab initio DMF calculations as used by Mies (1974) and van der Loo and Groenenboom (2007, 2008) are more uncertain for internuclear distances between the O and H atom that are far from the equilibrium. Moreover, the input data for empirical DMFs (Turnbull and Lowe, 1988, 1989; Nelson et al., 1990), theoretically extended empirical DMFs (Goldman et al., 1998), and modiﬁed ab initio DMFs (Langhoff et al., 1986; Brooke et al., 2016) were mainly restricted to low v or low v. As discussed in Sect. 3.2, about half of the measured doublets are resolved due to the high spectral resolving power of UVES. This allowed us to systematically study deviations between the Einstein-A coefﬁcients of the e and f components. The older sets of transition probabilities (M74, LWR86, and TL89) do not provide information on the individual components. The 18 1.5 A OH(6-2)P (4) P 1 0.8 1.0 0.6 0.5 0.4 0.0 0.2 0.8504 0.8505 0.8506 Wa.elength [μm] 0.0 B OH(6-2)P (12) f 2 0.10 190 K 40.2 0.0μ −0.4 −0.6 Brooke et al. (2016) 0.00 −10 −5 0 5 10 15 0.8981 0.8982 0.8983 ′ ′ ′ ′ −1 Wa.elength [3m] E (p = f) − E () = e) [cm ] Figure 9. Difference in logarithmic OH level populations y of the separately measured components of reliable doublets indicated by the upper level parity f and e as a function of the corresponding difference in the upper level energy E . The results for different branches (P, Q, and R) are marked by different symbols and colours (see legend). Examples for small and large population discrepancies are marked by the letters A and B. On the right-hand side, the corresponding -doublet spectra are plotted with indicated line identiﬁcations. The dashed line in the main plot indicates the effect of a thermal population with a temperature of 190 K on y. HITRAN database (Gordon et al., 2017) contains these components but the Einstein-A coefﬁcients were just set to the value of the corresponding doublet. Finally, vdLG08 and B+16 consider doubling but the differences between the coefﬁcients are very small. For B+16, the mean relative difference for our sample of 723 doublets is only 0.04%. The largest deviations are related to P- and Q-branch lines with high N . The maximum in our sample of 0.25% is linked to OH(5-1)Q (6). As the corresponding values for vdLG08 are almost identical, it is sufﬁcient to use only B+16 coefﬁcients for the comparison of the -doublet components. Figure 9 shows the results for 185 reliable (class 3) doublets with resolved components. The logarithmic population ratio y 0 0 for f minus e for the upper state parity p is plotted as a function of the corresponding difference in E . The latter is negative for F = 2 lines of the P and R branch according to the parity deﬁnition used by Brooke et al. (2016). If the theoretically predicted ΔyΔ orΔdi erentΔ-2 doublet parities p −1 Radiance [R )m ] Radiance [R )m ] equality of the transition probabilities was true, y should be close to 0. Small deviations in the populations are possible due to the small differences in E . Assuming a Boltzmann-like distribution for a typical kinetic temperature of 190 K at altitudes 0 1 of the OH emission layer at Cerro Paranal (Noll et al., 2016), there would be y = +0:08 for E = 10 cm and the same amount with negative sign for the corresponding positive energy difference. However, the true y are about 1 order of magnitude larger. The average absolute discrepancy in y is 0.22. In the case of large energy differences of at least 5 cm , it would even be 0.37, which corresponds to a ratio of 1.4. The clear differences in the strengths of the -doublet components are also illustrated by two example spectra. The weakly separated OH(6-2)P (4) doublet already shows slight differences in the intensity (y = 0:04). For the widely separated OH(6-2)P (12) pair, the f component is about 1.7 times brighter than the e component (y = 0:53). It is astonishing that it appears that these large effects have not been recognised, so far. They can only be explained by inadequate Einstein-A coefﬁcients since the P, Q, and R branches behave differently. The y related to P-branch lines could be ﬁtted by a non-LTE Boltzmann-like distribution with about 20 K. However, the y distributions for the Q and R branches are less clear. A convincing regression line cannot be drawn, and even if the ﬁtting was performed, the slopes would be very different. This rules out a signiﬁcant impact of possible non-LTE-inducing propensity differences for the population of the e and f states, even if the parity deﬁnitions are changed or more complex dependencies involving several level parameters are considered. 4.3 Correction of Brooke et al. coefﬁcients The discussion in the previous section has shown that the currently available sets of Einstein-A coefﬁcients are not satisfactory, especially with respect to Q-branch lines and -doublet components. Overall, the B+16 set is the most promising since it is the most complete in terms of the included lines, shows the smallest band-to-band variations for constant v , and is only slightly worse than LWR86 with respect to the population deviations between different branches. Hence, we focus on the B+16 set for the rest of this paper. However, as the remaining issues can still negatively affect the evaluation of OH level population distributions as shown in Fig. 5, we tried to improve the coefﬁcients to result in more consistent population ratios in the diagnostic plots discussed in Sect. 4.2. Our approach is fully empirical, i.e. it is based on regression lines and correction factors, and consists of three steps related to the correction of the discrepancies revealed by Figs. 6, 8, and 9. Complex ﬁtting approaches involving theoretical DMF-based calculations of the Einstein-A coefﬁcients are out of the scope of this study. We started with a correction of the discrepancies between the populations derived from different branches as shown in Fig. 6. For this purpose, we used 184 suitable, highly reliable doublets classiﬁed as 33, i.e. with a primary and secondary class of 3 (Sect. 3.2). It turned out that a linear ﬁt of y for Q versus P and R versus P shows the best performance for N as the independent variable (instead of the plotted E ) and a separate ﬁt for each OH band. The resulting slopes and intercepts and their uncertainties are provided in Fig. 10. There are no data points for OH(5-0), OH(9-5), and OH(4-1) due to relatively high uncertainties caused by an insufﬁcient number of reliable line measurements. For the remaining 12 bands, 5 to 16 (3 to 8) doublets could be used for the ﬁt of the difference between R and P branch (Q and P branch). The slopes for the Q branch are clearly more negative than those for the R branch. The mean values are 0:15 and 0:02, i.e. the apparent populations decrease by 14 and 2% compared to the P branch for an increase of N by 1. The corresponding mean intercepts are +0:11 and +0:01, 20 ' B(ooke et al. (2016) 0.0 −0.1 R-P (band)) R-P (fit) −0.2 Q-P (band)) Q-P (fit) −0.3 0.6 0.7 0.8 0.9 1.0 0.3 0.2 0.1 0.0 −0.1 0.6 0.7 0.8 0.9 1.0 Wavelengt of Q (1) [μm] Figure 10. Fit of the differences in logarithmic OH level populations for Q (crosses) and R branches (horizontal bars) relative to the P branch as a function of the upper rotational level N for each band where a reliable linear regression could be performed. The symbols with error bars show the resulting slope (upper panel) and intercept (lower panel) and their respective uncertainties. These values based on Einstein-A coefﬁcients of Brooke et al. (2016) were also ﬁtted by performing error-weighted linear ﬁts depending on the wavelength of Q (1) for each band. The ﬁt lines are displayed. i.e. nearly zero for R versus P. The data points are plotted as a function of the central wavelength of the Q (1) doublet of each band as especially the slope indicates a signiﬁcant trend with wavelength, which suggests that the branch-related errors in the B+16 Einstein-A coefﬁcients depend on the energy difference between the upper and lower state. The discrepancies between the branches appear to be smaller for longer wavelengths and might even vanish around 1 m. We used these correlations for a more robust correction approach involving error-weighted linear ﬁts of slope and intercept as a function of Q (1) wavelength. The error weighting allowed us to consider the strong dependence of the uncertainties on the band. The resulting ﬁts are also shown in Fig. 10. The slopes change with +0:45 0:05 (Q) and +0:14 0:03 (R) per micrometre. For the intercepts, we obtain changes of 0:36 0:08 (Q) and +0:04 0:11 (R) per micrometre. According to Eq. (1), the reciprocals of the population ratios y from our two-step ﬁtting procedure indicate the systematic deviations of the Einstein-A coefﬁcients. Thanks to the derived regression lines, they can also be predicted for missing lines and bands, at least in the covered wavelength range. For the correction, we need to deﬁne a reference. While the Q-branch coefﬁcients do not appear to be reliable in general, it would be arbitrary to choose the P or the R branch. However, Fig. 7 suggests that the errors in the Einstein-A coefﬁcients are minimised if both branches are combined with the same weight. Hence, we corrected the transition probabilities for the three branches by using half the ﬁtted deviation between P and R branch as the reference. Figure 11 (left) shows y for the 184 considered pairs of doublets before and after the correction. Intercept for fit of .y ′ ′′(N) Slope for fit of .y ′ ′′(N) v , v v , v BΔ16 B+16 B+16 < Δy > ,Δσ 0.02, 0.25 0.02, 0.11 0.5 0.0 0.0 0.0 −0.5 < Δ− > ,σ < Δ− >, σ −0.5 -0.13, 0.16 -0.02, 0.11 −0.5 -0.00, 0.08 -0.01, 0.05 0 1000 2000 3000 0.6 0.7 0.8 0.9 1.0 −10 0 10 ′ ′ ′ ′ ′ ′ −1 ′ ′ ′ ′ −1 Wa,elen th [1m] E (v , N , F ) − E (, , 1, 1) [cm ] E (p = f) − E (p = e) [cm ] Figure 11. Correction of Brooke et al. (2016) Einstein-A coefﬁcients. The measured (circles) and corrected (crosses) OH level population ratios y are given for changes in the branch (left; c.f. Fig. 6), lower vibrational level v (middle; c.f. Fig. 8), and the -doublet component (right; c.f. Fig. 9). Each plot lists the mean y and standard deviation for the measured and corrected data. While in the former case the mean value and standard deviation are 0:13 and 0:16, the corrected data reveal 0:00 and 0:08. Thus, the offsets vanished completely on average and the scatter was reduced by a factor of 2. In the next step, we corrected y offsets between lines differing only in v as shown in Fig. 8. For this purpose, we focused 0 0 on relatively bright lines with N 4 for the P and R branch and N = 1 for the Q branch. As these doublets indicate similar y (scatter of 0.03), differences in the selected line subsets do not critically affect the mean values. For each considered OH band, 8 to 14 reliable pairs of doublets of class 33 were available (89 in total). For the correction of y by changing the Einstein-A coefﬁcients, we used the same reference bands for each v as discussed in Sect. 4.2. The choice is motivated by the accessibility and quality of the line measurements with UVES. It does not necessarily include the bands with the most realistic transition probabilities. As discussed in Sect. 4.2, coefﬁcients of bands with low v and v tend to be more reliable as the DMF calculations are less challenging and the experimental data are more abundant. Our reference bands have v between 3 and 9 and v between 3 and 5. Hence, the most promising bands are beyond the UVES wavelength range. Nevertheless, there does not appear to be a strong quality gradient with v or wavelength (v cannot be tested) since the B+16 coefﬁcients do not show such a dependence of y for the covered bands in Fig. 8 (in contrast to the other investigated sets). Note that this is different from the situation for the branches illustrated in Fig. 10. In the end, we shifted the mean y for eight OH bands to zero by multiplying the Einstein-A coefﬁcients by factors between 0.85 for OH(7-2) and 1.06 for OH(8-4). This reduces the scatter in the measured populations for ﬁxed v in any case, even if the choice of the reference bands might not be optimal. Figure 11 (middle) shows the corresponding results for 143 doublets. The discussed corrections (also including the branch-related modiﬁcations) change the mean y and standard deviation from 0:02 and 0:11 to 0:01 and 0:05. Finally, we corrected the y between the -doublet components as shown in Fig. 9. This is necessary in order to also use doublets with only one reliable component for the study of the OH level populations discussed in Sect. 5. For the change Δ− fo( diffe(ent b(anche) ′′ Δ−ΔforΔdifferentΔv Δ− fo( diffe(ent 0-doublet parities p of the Einstein-A coefﬁcients, we assumed a natural population discrepancy between the two components consistent with a temperature of 190 K as illustrated in Fig. 9. We ﬁtted the remaining y for each branch using N as the independent variable due to a better performance with respect to linear regressions compared to E . This approach requires us to ﬂip the sign for the data points with negative E . This is reasonable since the amount of the deviations is very similar for doublets with F of 1 and 2. For the ﬁts related to P, Q, and R, we considered 99, 18, and 11 resolved doublets of class 33. The resulting slopes (and intercepts) are 0:0350:003 (+0:130:03), +0:0650:019 ( 0:190:08), and 0:0350:035 (+0:430:48), respectively. Hence, the effect for the Q branch seems to be twice as large as for the P branch and to also have a different sign. The R branch might behave similarly to the P branch but the uncertainties are high. Assuming that the e and f components equally contribute to the ﬁtted differences, we corrected the Einstein-A coefﬁcients for the reliable N range of the ﬁts, i.e. we neglected doublets with low N where both components are not sufﬁciently separated or the ﬁt crossed y = 0. Figure 11 (right) shows the resulting change in the mean y and scatter for the investigated 128 doublets. While the small mean value of +0:02 did not signiﬁcantly change, the standard deviation was clearly reduced from 0:25 to 0:11. 5 OH level populations 5.1 Mean populations and rotational temperatures With the correction of the B+16 Einstein-A coefﬁcients in Sect. 4.3, we minimise the scatter in the OH level populations for upper states with different measured lines. Moreover, the change of the populations with N appears to be more reliable due to the promising combination of P- and R-branch data for T estimates (Fig. 7). The resulting population distributions for v rot 0 0 between 4 and 9 are shown in Fig. 12. We neglect v = 3 due to the lack of high-N states in the UVES data (Fig. 3), which does not allow us to describe this population distribution in detail. As already brieﬂy discussed in Sects. 1 and 4.1, there is a characteristic pattern for increasing level energy with a steep population decrease for low N and a rather slow decrease for 0 0 0 high N . The difference between high and low N tends to increase with decreasing v . Moreover, the change between the two extremes does not appear to happen continuously with increasing N . Instead, the transition is mostly localised in a narrow E interval of a few hundred inverse centimetres. This known pattern (Cosby and Slanger, 2007; Oliva et al., 2015) suggests the deﬁnition of a cold and a hot population for each v , which can be described by corresponding T and population ratios rot for ﬁxed level energies (Oliva et al., 2015; Kalogerakis et al., 2018; Kalogerakis, 2019). We applied this concept by ﬁtting the natural logarithm of the sum of two exponential Boltzmann terms as a function of 0 0 0 0 0 E (E relative to the energy for N = 1 and F = 1) to the corrected y for each v . We considered the populations from all doublets with quality classes above 0. For classes 1 and 2, we derived the doublet-related populations from the reliable 0 0 components. An inspection of the change of the populations with increasing E resulted in the rejection of the highest N levels for v = 8 and 9 as the related populations cannot be reproduced satisfactorily by a two-component ﬁt. In the case of 0 0 0 1 v = 9, the seven measurements with N 10 were neglected. All corresponding E are between 250 and 790 cm above the exothermicity limit of the hydrogen–ozone reaction (Sect. 3.2), which can explain the rapid population decrease in this energy range, which could not clearly be constrained before due to a lack of data (Cosby and Slanger, 2007; Noll et al., 2018b). In 23 4 ′ ′ ′ v = 4 v = 5 v = 6 N = 83 N = 88 N = 124 sel sel sel T = 196.4±4.4 K T = 193.0±4.6 K T = 189.9±3.0 K cold cold cold T = 6636±424 K T = 5952±569 K T = 4544±201 K hot hot hot r = 0.33±0.03 % r = 0.37±0.04 % r = 0.51±0.03 % pop, 0 pop, 0 pop, 0 −2 −4 −6 4 ′ ′ ′ − = 7 − = 8 − = 9 N = 113 N = 124 N = 86 sel sel sel T = 187.203.4 K T = 190.502.7 K T = 189.6±5.7 K cold cold cold T = 2181054 K T = 1336047 K T = 686±116 K hot hot hot ) = 1.1300.06 % ) = 1.6800.11 % ) = 2.94±1.11 % pop, 0 pop, 0 pop, 0 H+O exotherm. selected rejected cold population hot population 2-component fit 0 2000 4000 6000 8000 0 2000 4000 6000 8000 0 2000 4000 6000 8000 ′ ′ ′ ′ ′ ′ 11 ′ ′ ′ ′ ′ ′ 11 ′ ′ ′ ′ ′ ′ 11 E (v , N , F ) − E (−, 1, 1) [cm ] E (v , N , F ) − E (−, 1, 1) [cm ] E (v , N , F ) − E (−, 1, 1) [cm ] Figure 12. Distribution of logarithmic OH level populations for the corrected Brooke et al. (2016) Einstein-A coefﬁcients. Except for OH(7-1), all doublets with at least one reliable component are considered. The populations are plotted in separate panels for each up- 0 0 per vibrational level v between 4 and 9. The given level energies in inverse centimetres are relative to the lowest v -related energy. The populations can be ﬁtted by means of a two-component rotational temperature ﬁt. All data points involved are marked by circles. The few exceptions (indicated by crosses) are close to or above the exothermicity limit of the hydrogen–ozone reaction, which is marked by vertical dotted lines. For these data, a more complex ﬁt would be necessary. The ﬁnal ﬁt curves are displayed by solid lines. The underlying cold and hot components are shown by dashed and dot-dashed lines, respectively. Each panel lists the v -speciﬁc number of selected data points N , sel the rotational temperatures of the cold and hot components T and T in kelvins, and the ratio r of the populations of both linear ﬁt cold hot pop;0 components for the lowest v -related energy in per cent. 0 0 0 the case of v = 8, a strong decrease of the populations is found for eight measurements related to N 14 with E between 1 1 660 cm below and 90 cm above the exothermicity limit. This is an interesting result as it provides valuable constraints on the nascent populations and the relaxation process from v = 9 to 8. The drop of the populations below the exothermicity limit y y co)) corr also seems to be present in the population distribution of Cosby and Slanger (2007), also based on UVES spectra (Hanuschik, 2003; Cosby et al., 2006). However, the authors do not discuss this phenomenon. Our v 7 data, which are related to energies of more than 1,200 cm below the limit, do not show a cut in the populations. The data of Cosby and Slanger (2007) are not conclusive here either. 0 0 0 The remaining population measurements for each v , which varied between 83 for v = 4 and 124 for v = 6 and 8, were ﬁtted with our two-component model by means of robust least-squares minimisation, which resulted in the same best ﬁts for a wide range of start values. Figure 12 shows the ﬁnal best ﬁts and also indicates the corresponding temperatures T and cold T as well as the ratio of the hot and cold populations for E = 0, r . Under consideration of the ﬁt uncertainties, the hot pop;0 best-ﬁt T are very similar and consistent with a temperature of 190 K, i.e. the typical ambient temperature at OH emission cold altitudes (Noll et al., 2016). Only v = 4 with 196 4 K might slightly be higher. This could point to the weak inﬂuence of an 0 0 intermediate population for low v . The second highest T value of 193 5 K for v = 5 would be in agreement with this cold interpretation. Note that ﬁxing the ﬁt to a T of 190 K did not signiﬁcantly change the other parameters. The differences cold were much smaller than the uncertainties. In contrast to T , T shows a strong trend with v . The temperatures increase cold hot 0 0 0 from about 700 K for v = 9 to about 7,000 K for v = 4. In parallel, r decreases from about 3% for v = 9 to about 0.3% pop;0 0 0 for v = 4, i.e. hot populations with higher T show lower contributions to the total population at low N . The strong change hot 0 0 in r appears to be mainly caused by the decrease of the cold population with increasing v since the E = 0 intercepts pop;0 of the lines describing the hot populations are located at similar y values of around 2 in Fig. 12. The ﬁts for v 7 (no rejection of states) are convincing with respect to the assumption of a homogeneous hot population, which can be described by a single temperature. Nevertheless, some ﬁne structure might exist as the comparison of the individual measurements and the ﬁt lines suggest, although the possible population deviations appear to be not larger than 30%, which is small compared to population changes of the order of a magnitude in the v -dependent energy ranges most contributing to T . Hence, the two- hot component ﬁts are quite robust as the listed errors show. The highest uncertainties are related to v = 9 since the hot population 1 0 is essentially constrained in an energy range of less than 300 cm , which includes 12 measurements with N of 8 and 9. Two-component ﬁts were previously performed by Oliva et al. (2015) based on a near-infrared GIANO spectrum with a resolving power of 32,000 taken during 2 hours with the spectrograph directly pointing to the night sky at the La Palma Observatory. The investigated lines belong to OH bands with low v and are complementary to those covered by our study. For the calculation of the populations, Oliva et al. (2015) used the Einstein-A coefﬁcients from van der Loo and Groenenboom (2007). For the ﬁts, T was ﬁxed at 200 K. The resulting T and r varied from about 1,300 K and 1.8% for v = 8 to cold hot pop;0 about 7,000 K and 0.23% for v = 4. Although errors were not reported, these values are in good agreement with our results provided in Fig. 12. The ﬁt parameters for v = 9 are highly uncertain. However, Oliva et al. (2015) succeeded in ﬁtting the 0 0 0 populations for v = 2 and 3, which show an extension of the trend found for the higher v . For v = 2, T and r resulted hot pop;0 in 12,000 K and 0.14%, respectively. The GIANO data were reﬁtted by Kalogerakis et al. (2018) with unconstrained T , cold which resulted in temperatures of about 190 K but with larger scatter than in our case. For T , the general trend was the same hot 0 0 but with a large step from 900 K for v = 8 to 4,000 K for v = 7, which disagrees with our ﬁndings. Population ratios were not provided by Kalogerakis et al. (2018). Noll et al. (2018b) already published populations related to v = 9 and P-branch 25 lines based on the UVES data used in this study and B+16 Einstein-A coefﬁcients. Kalogerakis (2019) ﬁtted these populations and found T and T of about 180 K and 500 K, respectively. Both temperatures are lower than our results, but less than cold hot 2 standard deviations. The ﬁt of Kalogerakis (2019) based on fewer data points seems to be related to a higher impact of the 0 0 hot population at low N . Our results for T allow for an interesting comparison to the T of the nascent populations of v hot rot between 7 and 9, which were derived by Llewellyn and Long (1978) using laboratory data from Charters et al. (1971). Our best-ﬁt T of 690 120, 1340 50, and 2180 50 K agree well with their 760 20, 1230 30, and 1940 200 K, which hot implies that the OH relaxation processes do not appear to signiﬁcantly affect the hot populations of the highest v . The previous discussion has shown that bimodality is a good concept for the description of the population distributions for 0 0 each v . Moreover, the derived T are close to the expected effective ambient temperatures for the v -dependent OH emission cold 0 0 layers (Noll et al., 2016). The increasing trend of T derived from the lines with the lowest N for increasing v (Cosby and rot Slanger, 2007; Noll et al., 2015, 2017) is not found in the best-ﬁt T . Hence, our ﬁts could be used to estimate the non-LTE cold contributions to such T , T , which are an issue for the use of T as indicators of the temperatures in the mesopause rot NLTE rot region. Kalogerakis et al. (2018) and Kalogerakis (2019) compared T ﬁts with T from linear regressions for levels with cold rot 0 1 0 E lower than 500 and 250 cm , respectively. The results indicate higher T than T at least for the highest v (order rot cold of 20 K). However, the uncertainties are large due to the strong impact of the line selection (Noll et al., 2015), uncertainties in the line intensities (unclear for the GIANO data), and the choice of the Einstein-A coefﬁcients (Fig. 7). Hence, we applied a different approach by directly taking the two-component ﬁt for the measurement of T . For this purpose, we derived the rot populations related to the ﬁrst three P -branch lines, which are often taken for T determinations (e.g. Schmidt et al., 2013; 1 rot 0 0 Noll et al., 2016), from the ﬁt curve at the corresponding v -dependent E . The related T were then calculated by a linear rot regression of the three y for each v . Finally, the resulting T is just the difference between T and T . This method NLTE rot cold is very robust as it is fully based on the two-component ﬁt, which relies on a high number of population measurements. Hence, uncertainties related to individual lines are negligible. Our T for v between 4 and 9 are shown in Fig. 13. They increase relatively slowly between 4 and 6 from 1:2 0:1 K NLTE 0 0 to 1:80:1 K and then faster to 6:01:6 K for v = 9. As the latter value is relatively uncertain, 4:60:3 K for v = 8 might also be the maximum deviation. The errors were derived by displacing the hot component ﬁt in both y directions according to the uncertainty in r and reﬁtting T as the only parameter to obtain modiﬁed T . This approach considers that r pop;0 hot NLTE pop;0 and T are anticorrelated and that the relative uncertainty in T is relatively small. Additional systematic uncertainties are hot cold caused by assuming only two components. It is required that the ﬁt line for the hot component can be linearly extrapolated to E = 0. The good quality of the ﬁts in the transition region between the dominance of the cold and hot components are promising. Nevertheless, contributions of additional components of intermediate temperature cannot be excluded (as e.g. for v = 4 due to a possibly elevated T ). Fits with such an additional component and a ﬁxed T of 190 K, where the best-ﬁt cold cold parameters of the intermediate and hot populations are not well constrained, showed possible T increases by 10 to 30%, NLTE i.e. the signiﬁcance of positive non-LTE effects for all v would remain high. It is hard to imagine situations where T NLTE could signiﬁcantly drop. A sharp cut of the hot population for low N would be inconsistent with a Boltzmann-like distribution as expected for relaxation processes. Another source of possible systematic errors are the Einstein-A coefﬁcients, especially 26 2-component fit Oliva+15 Noll+16 Noll+16 adapted 5 − 2 3 4 5 6 7 8 9 v' Figure 13. Estimate of non-LTE contributions to rotational temperatures for the three lowest rotational lines with F = 1 as a function of the upper vibrational level v . The differences of the T for the two-component ﬁts and the related T shown in Fig. 12 are provided rot cold by circles and solid lines. Error bars (which are only larger than the symbols for v = 9) are also given. Two-component ﬁts of OH level populations were also performed by Oliva et al. (2015) based on data from the near-infrared spectrograph GIANO at La Palma Observatory in Spain. We used their ﬁt parameters (no errors) to calculate T for v from 2 to 9 (plus signs and dotted lines). For comparison, NLTE the ﬁgure also indicates the estimates by Noll et al. (2016) based on X-shooter spectra from Cerro Paranal and SABER data for a similar area (diamonds and dashes). As these T were derived for Einstein-A coefﬁcients from HITRAN, we recalculated the non-LTE effects NLTE using our corrected Brooke et al. (2016) coefﬁcients (crosses and dot-dashed lines). with respect to their dependence on N . The latter was changed in Sect. 4.3 for the B+16 coefﬁcients by the branch-speciﬁc corrections. Hence, we tested what happens if we consider the P- or R-branch data as the standards instead of a combination of both. These modiﬁcations would change directly measured T to values as indicated in Fig. 7. However, the effect on the rot two-component approach is much smaller. It is of the order of the already small ﬁt errors for all v . As expected, non-LTE contributions related to R as the standard are lower than those for P. Furthermore, we investigated the inﬂuence of the choice of the energy levels on T by also simulating line sets consisting of the ﬁrst two and ﬁrst four P -branch lines. For v = 8 as NLTE 1 an example, these changes cause T of 3.3 and 6.9 K, which clearly deviate from the plotted 4.6 K. Hence, the non-LTE NLTE contributions are very sensitive to the selected energy levels, which is consistent with the results from Noll et al. (2015, 2018b). As shown in Fig. 13, we also calculated T from the two-component ﬁts of Oliva et al. (2015). Excluding their very NLTE uncertain ﬁt parameters for v = 9, there is a very good agreement with differences smaller than 0.4 K. The only exception is v = 7, where our non-LTE contributions are about 1.6 K lower. As the data basis and analysis were completely different (including different Einstein-A coefﬁcients), this convincing result demonstrates the robustness of the approach. The GIANO- 0 0 related data for v = 2 and 3 suggest that T decreases only very slowly with decreasing v . The drop in r seems to NLTE pop;0 be nearly compensated by the increase in T . hot ΔT Δ[K] NLTE 0 Mean T for v from 2 to 9 at Cerro Paranal were already derived by Noll et al. (2016) based on measurements NLTE of 25 OH bands and 2 O bands (where non-LTE effects are less important) in optical and near-infrared spectra from the echelle spectrograph X-shooter as well as from OH emission and kinetic temperature proﬁle measurements with SABER. The T from the complex analysis for the ﬁrst three P -branch lines (derived from different band-speciﬁc line sets) are NLTE 1 shown in Fig. 13. The values with a conspicuous maximum of 13:2 2:0 K at v = 8 are clearly higher than those from the two-component ﬁt. However, Noll et al. (2016) used HITRAN Einstein-A coefﬁcients, which signiﬁcantly deviate from our modiﬁed B+16 coefﬁcients. As demonstrated by Fig. 7, the impact on T can be large. Hence, we recalculated the T rot NLTE of Noll et al. (2016) with the modiﬁed B+16 transition probabilities for the lines considered in Sect. 4.3 and the original ones (which result in about 2 K higher T on average) in all other cases. Figure 13 indicates a clear reduction of T for the rot NLTE 0 0 relevant v 4, which better matches our results based on two-component population ﬁts. Between v of 4 and 7 the non-LTE contributions are almost constant with a mean of 2.5 K. However, the absolute uncertainties are larger. Only the maximum of 8:4 2:7 K at v = 8 seems to be signiﬁcant. It might also be present (but less pronounced) in the population ﬁtting results. The high absolute uncertainties from temperature comparisons (v -related differences are safer) are a critical drawback of that method and imply that two-component population ﬁts provide the best constraints for T , so far. The higher errors NLTE compared to the original Noll et al. (2016) data are partly related to the unavoidable mixture of corrected and uncorrected B+16 coefﬁcients. However, B+16 line parameters also appear to cause a larger scatter in T for different bands with the same v rot 0 0 compared to HITRAN data. The change in the T differences between adjacent v (especially around v = 7) is mainly NLTE caused by a different calculation of T for the reference line set consisting of the ﬁrst three P -branch lines. Instead of using rot 1 a constant temperature offset for the conversion from the reference line set of Noll et al. (2015) including all P-branch lines up to N = 3 as discussed by Noll et al. (2016), we directly corrected the band-speciﬁc T to be representative of the more rot recent reference line set. 5.2 Population variability The discussion of the roto-vibrational level populations of OH in Sect. 5.1 was only based on line intensity measurements in a single mean spectrum. We can learn more about these populations if we also consider variations in the emission layer properties. In order to keep the signal-to-noise ratios high, we split the sample into two parts based on a characteristic layer parameter, calculated the corresponding mean spectra, and derived level populations from the measured line intensities for a comparison. For the split, we selected the effective height of the OH emission layer h , i.e. the centroid altitude weighted by the volume emission rate, as it is positively correlated with the strength of the non-LTE effects (Noll et al., 2017, 2018a). There are fewer thermalising collisions without v change at higher altitudes due to lower air densities but higher atomic oxygen mixing ratios (Noll et al., 2018b). The impact of this effect is clearly reﬂected by the observed higher h for higher v (e.g., von Savigny et al., 2012), which are more affected by the hot nascent population and have lower effective lifetimes. According 0 0 to the population modelling of Noll et al. (2018b) for v = 9, h also increases for higher N . In order to study the change of the population distribution for the different v depending on the OH emission altitude, we need adequate space-based measurements of the emission proﬁles to be linked with our ground-based OH level population 28 ′ 0.10 v = 4 v = 5 v = 6 v = 7 v = 8 0.05 v = 9 0.00 .0.05 UVES sample split: SABER 2.06 μm median h −0.10 0 1000 2000 3000 ′ ′ ′ ′ ′ ′ −1 E (v , N , F ) − E (v , 1, 1) [cm ] Figure 14. Difference in logarithmic OH level populations for UVES mean spectra based on equally sized subsamples of high and low effective emission height (89.61 and 88.72 km on average). The altitudes were derived from volume emission proﬁles of the SABER OH channel at 2.06m for the Cerro Paranal region. The data points are shown for different upper vibrational levels v (indicated by the corresponding coloured markers) as a function of the level energies in inverse centimetres relative to the lowest v -related energy. The plot only shows population measurements for reliable doublets which are covered by both UVES set-ups (see Fig. 1 for the valid wavelength ranges). data. This was already achieved by Noll et al. (2017) based on limb-sounding data for the Cerro Paranal region from the OH- speciﬁc channels of the SABER radiometer (Russell et al., 1999). Here, we focus on the channel centred on 2.06 m, which covers OH(8-6) and OH(9-7). The effective v is about 8.3 (Noll et al., 2016). Noll et al. (2017) connected the resulting h for 4,496 proﬁles to each UVES spectrum by a weighting procedure which involved temporal differences in day of year and local time measured in a two-dimensional climatology. The approach also included the correction of differences in the solar activity, as measured by the solar radio ﬂux (Sect. 3.2), for each UVES observation compared with the corresponding weighted h . Finally, a most likely h was available for each UVES spectrum. We used these data to split our sample of 533 spectra (Sect. 2) at a median h of 89.2 km for the 2.06 m OH channel. The resulting subsamples show mean h of 88.7 and 89.6 km, i.e. e e the height difference is almost 1 km. The situation is very similar for the other OH channel at 1.64 m representing an effective v of about 4.6 (Noll et al., 2016), where the corresponding heights are 87.3 and 88.3 km. Then, we performed the entire data analysis starting with the calculation of the mean spectra up to the derivation of the ﬁnal populations. In order to minimise systematic effects in the line measurement, the same wavelengths for the line integration and continuum derivation as for the full sample spectrum were used (see also Sect. 3.2). Figure 14 shows the resulting population ratios for the high and low h cases. For such a comparison, the choice of the Einstein-A coefﬁcients does not matter. We only plot y related to the 257 best-quality doublets (class 33) that are covered by all UVES spectra, i.e. reliable population ratios which are representative of the given h . We can identify three energy level regimes with respect to the population change by a rise of the OH emission layer. y(h = 89.61 km) . y(h = 88.72 km) eff eff 1 0 0 Up to about 600 cm , there is a general decrease of the populations, which strongly depends on v . From v = 4 to 9, the decrease shrinks from about 8% ( 0:08) to about 1%, i.e. the populations for higher v are more stable. The largest difference in y is between the two lowest v . The decrease of the OH intensity for a rising OH emission layer is well known (Yee et al., 1997; Melo et al., 1999; Liu and Shepherd, 2006). It is accompanied by a lower width of the layer ( 0:5 km for our subsamples) due to especially low OH production rates at the bottom side, which are caused by a depletion of ozone. As band emissions with lower v peak at lower altitudes (von Savigny et al., 2012; Noll et al., 2016), they seem to be more affected by this lack of fuel for the OH production. 0 1 At E above 1,200 cm , Fig. 14 shows a completely different behaviour. There is a general increase of the populations with a mean of +0:04 and no signiﬁcant dependence on v . This ﬁnding implies that the contribution of hot populations to the total population increases with h . The impact of non-LTE effects grows by a less efﬁcient thermalisation process. The reduced contributions from lower altitudes to the total emission certainly play an important role here since (similar to v ) emission related to higher N peaks higher in the atmosphere (Noll et al., 2018b). There, collisional thermalisation of the rotational level population distributions is hampered by a relatively low density of nitrogen molecules and a relatively high volume mixing ratio of v -deactivating (or even OH-destroying) atomic oxygen radicals (Noll et al., 2018b). Note that the location of the zero line in Fig. 14 is uncertain with respect to the degree of thermalisation as the increase of 4% could also be caused by a change in the OH column density. If the hot populations deﬁne the zero line as they appear to be the most stable 0 0 ones (which might be supported by the lack of a v dependence), the low-N populations would further decrease. Figure 14 is another good argument for the bimodality of rotational level population distributions. The change of y for 0 0 low and high N with increasing E seems to be very small. This suggests that cold and hot populations are relatively homogeneous, which supports our two-temperature ﬁt approach. Moreover, there is a quick transition between both populations 0 1 in the relatively narrow E range between 600 and 1,200 cm . In Fig. 12, this is the region where both ﬁt components signiﬁcantly contribute. Hence, a rise of the OH emission layer there should have the strongest impact on the slope of the population distribution (i.e. T ) by a change of the relative contribution of the cold and hot populations. As r increases, rot pop;0 there is also an effect on T as shown in Fig. 13. Estimates of this quantity based on the populations for low and high OH NLTE layer are relatively uncertain due to the distinctly lower number of suitable lines and the high impact of increased line intensity errors on the analysis of very small population differences. Nevertheless, we calculated T changes of the order of a few NLTE tenths of a kelvin for the altitude difference of about 1 km. This is clearly smaller than about 1 K per kilometre, the order of magnitude from observational and modelling studies by Noll et al. (2017, 2018a, b), which might point to limitations in the study of T variations based on two-component population ﬁts. NLTE 6 Conclusions Based on averaged high-quality high-resolution spectra from the UVES echelle spectrograph at Cerro Paranal, we performed a detailed study of OH roto-vibrational level population distributions. The mean populations for 723 doublets with upper 0 0 vibrational levels v between 3 and 9 and upper rotational levels N up to 24 were investigated. In about half the cases, the 30 doublet components were measured separately. The line wavelengths from literature (Rothman et al., 2013; Brooke et al., 2016) 0 0 turned out to be sufﬁciently accurate in most cases. Only a small number of lines with high N and intermediate v (especially v = 5) showed deviations by more than a few picometres. The quality of population measurements is limited by uncertainties in the Einstein-A coefﬁcients. We investigated this issue with comparisons of populations from different transitions with the same upper state. We tested six sets of transition probabilities: Brooke et al. (2016), HITRAN (Rothman et al., 2013), van der Loo and Groenenboom (2008), Turnbull and Lowe (1989), Langhoff et al. (1986), and Mies (1974). All sets fail in the case of Q-branch lines and the -doublet components, where unexpectedly large intensity ratios are possible. The comparison of populations from P- and R-branch lines indicated relatively small errors for the coefﬁcients by Langhoff et al. (1986), van der Loo and Groenenboom (2008), and Brooke et al. (2016), whereas those from Turnbull and Lowe (1989) are clearly the worst. The comparison of OH bands with the same v showed a similar order of the different sets with respect to their quality. For this case, the coefﬁcients of Brooke et al. (2016) performed best. The widely used HITRAN data are only of intermediate and hence unsatisfactory quality. For the population analysis, we focused on the Einstein-A coefﬁcients from Brooke et al. (2016) due to their relatively good performance and the highest number of included lines. In order to minimise the scatter in the populations, we further improved these coefﬁcients by empirically correcting the found population discrepancies via regression lines related to N and wavelength as well as band-dependent correction factors. For the correction of the branch-related differences, we used P- and R-branch data combined with equal weights as the reference since this strongly reduced the deviations between the different sets of Einstein-A coefﬁcients with respect to rotational temperatures T , i.e. the change of the populations with increasing rot N . The whole correction procedure lowered the discrepancies in the coefﬁcients by more than a factor of 2 for the measured lines. Nevertheless, the development of an improved set for all lines would need a more sophisticated approach including modelling of the molecular parameters. The resulting v -dependent population distributions show clearly bimodal structures, which were convincingly reproduced 0 0 by two-temperature ﬁts only excluding steep population decreases for v = 8 and 9 at the highest N with energies slightly below and above the exothermicity limit of the OH-producing hydrogen–ozone reaction, respectively. The ﬁts show a cold population with nearly ambient temperature of about 190 K dominating at low N and a hot population with temperatures 0 0 0 between 700 K for v = 9 and 7,000 K for v = 4 at high N . In contrast, the ratio of the hot and cold populations at the level with the lowest energy of a given v changes from 3 to 0.3% mainly due to a decrease of the cold component. The signiﬁcant contribution of a hot population to low N causes deviations between T and ambient temperature, which we estimated by rot ﬁtting our two-component model for the energy levels related to the ﬁrst three P -branch lines. The results indicate non-LTE 0 0 0 contributions that increase from about 1 K for v = 4 to about 5 K for v = 8. The best-ﬁt value for v = 9 is even higher (about 6 K), but the ﬁt uncertainties are by far the highest. In general, the applied approach is much more robust than the previously used method based on comparisons of temperatures from different sources as it only weakly depends on uncertainties in the line intensities and Einstein-A coefﬁcients. Our approach is mostly limited by the reliability of the assumption of only two Boltzmann-like population distributions. There are hints of the existence of a more complex pattern but the impact of these additional components appears to be small. 31 This conclusion is supported by the change of the populations due to a rise of the OH emission layer, which we studied by the separation of the sample of spectra into two parts depending on the effective emission height as obtained from height- 1 1 resolved SABER OH volume emission rates. The energy regimes up to about 600 cm and above about 1,200 cm relative to the lowest energy for a given v show clearly distinct variability in agreement with the energy ranges dominated by the cold and hot components in the derived population distributions. While the cold populations show a decrease, which is stronger for lower N , the hot populations are relatively stable (or even increase) with increasing emission altitude. The largest measured effect is a 12% decrease of the cold population at v = 4 relative to the hot population for a height difference of almost 1 km. The success of the two-component model for OH rotational level population distributions has implications for the thermalisa- tion process of the highly non-thermal nascent populations. There are still high uncertainties with respect to the rate coefﬁcients 0 0 for collisions with and without change of v . In particular, the modiﬁcation of the rotational level population by v -changing 0 0 collisions is not known. Hence, the origin of the very hot populations at high N of low v is puzzling. Consequently, there is hope that the high-quality population data of this study can help to better understand relaxation processes in OH by detailed modelling. This will be important knowledge with respect to the use of OH as an indicator of mesopause temperatures and for retrievals of atomic abundances like those of oxygen. Data availability. This project made use of the ESO Science Archive Facility at http://archive.eso.org. UVES Phase 3 spectra from different observing programmes of the period from April 2000 to March 2015 were analysed. The v2.0 SABER data products used for this study were taken from http://saber.gats-inc.com. Author contributions. S. Noll has developed the project, processed the data, performed the analysis, produced the ﬁgures, and is the main author of the paper text, where all co-authors have made signiﬁcant contributions. H. Winkler and O. Goussev have also inﬂuenced the design of the study. In addition, H. Winkler has checked parts of the analysis and B. Proxauf has been involved in the post-processing of the UVES Phase 3 products. Competing interests. The authors declare that they have no conﬂict of interest. Acknowledgements. We thank reviewer E. Oliva and one anonymous referee for their positive and helpful reports. S. Noll is ﬁnanced by the project NO 1328/1-1 of the German Research Foundation (DFG). H. Winkler is funded by the DFG project NO 404/21-1. 32 References Adler-Golden, S.: Kinetic parameters for OH nightglow modeling consistent with recent laboratory measurements, J. Geophys. Res., 102, 19 969–19 976, https://doi.org/10.1029/97JA01622, 1997. Baker, D. J. and Stair, Jr., A. T.: Rocket measurements of the altitude distributions of the hydroxyl airglow, Phys. Scripta, 37, 611–622, https://doi.org/10.1088/0031-8949/37/4/021, 1988. 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. Beig, G., Keckhut, P., Lowe, R. P., Roble, R. G., Mlynczak, M. G., Scheer, J., Fomichev, V. I., Offermann, D., French, W. J. R., Shepherd, M. G., Semenov, A. I., Remsberg, E. E., She, C. Y., Lübken, F. J., Bremer, J., Clemesha, B. R., Stegman, J., Sigernes, F., and Fadnavis, S.: Review of mesospheric temperature trends, Rev. Geophys., 41, RG1015, https://doi.org/10.1029/2002RG000121, 2003. Brooke, J. S. A., Bernath, P. F., Western, C. M., Sneden, C., Afsar, ¸ M., Li, G., and Gordon, I. E.: Line strengths of rovibrational and rotational transitions in the X ground state of OH, J. Quant. Spectrosc. Radiat. Transf., 168, 142–157, https://doi.org/10.1016/j.jqsrt.2015.07.021, Charters, P. E., MacDonald, R. G., and Polanyi, J. C.: Formation of vibrationally excited OH by the reaction H + O3, Appl. Optics, 10, 1747–1754, https://doi.org/10.1364/AO.10.001747, 1971. 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. Cosby, P. C., Slanger, T. G., Huestis, D. L., and Osterbrock, D. E.: Term energies, line positions, and spectroscopic constants for the OH Meinel band system, in: 55th International Symposium on Molecular Spectroscopy, Ohio State University, Columbus, Ohio, USA, https://www.asc.ohio-state.edu/miller.104/molspect/symposium_55/symposium/Abstracts/p319.pdf, 2000. Cosby, P. C., Sharpee, B. D., Slanger, T. G., Huestis, D. L., and Hanuschik, R. W.: High-resolution terrestrial nightglow emission line atlas from UVES/VLT: Positions, intensities, and identiﬁcations for 2808 lines at 314-1043 nm, J. Geophys. Res., 111, A12307, https://doi.org/10.1029/2006JA012023, 2006. Dekker, H., D’Odorico, S., Kaufer, A., Delabre, B., and Kotzlowski, H.: Design, construction, and performance of UVES, the echelle spectrograph for the UT2 Kueyen Telescope at the ESO Paranal Observatory, in: Optical and IR Telescope Instrumentation and Detectors, edited by Iye, M. and Moorwood, A. F., vol. 4008 of SPIE Proc. Ser., pp. 534–545, https://doi.org/10.1117/12.395512, 2000. Dodd, J. A., Armstrong, P. S., Lipson, S. J., Lowell, J. R., Blumberg, W. A. M., Nadile, R. M., Adler-Golden, S. M., Marinelli, W. J., Holtzclaw, K. W., and Green, B. D.: Analysis of hydroxyl earthlimb airglow emissions: Kinetic model for state-to-state dynamics of OH(v,N), J. Geophys. Res., 99, 3559–3586, https://doi.org/10.1029/93JD03338, 1994. Edlén, B.: The Refractive Index of Air, Metrologia, 2, 71–80, https://doi.org/10.1088/0026-1394/2/2/002, 1966. Franzen, C., Espy, P. J., Hofmann, N., Hibbins, R. E., and Djupvik, A. A.: Airglow Derived Measurements of Q-Branch Transition Probabil- ities for Several Hydroxyl Meinel Bands, Atmosphere, 10, 637, https://doi.org/10.3390/atmos10100637, 2019. French, W. J. R., Burns, G. B., Finlayson, K., Greet, P. A., Lowe, R. P., and Williams, P. F. B.: Hydroxyl (6-2) airglow emission intensity ratios for rotational temperature determination, Ann. Geophys., 18, 1293–1303, https://doi.org/10.1007/s00585-000-1293-2, 2000. Goldman, A., Schoenfeld, W. G., Goorvitch, D., Chackerian, Jr., C., Dothe, H., Mélen, F., Abrams, M. C., and Selby, J. E. A.: Updated 2 2 0 00 line parameters for OH X II-X II (v ,v ) transitions., J. Quant. Spectrosc. Radiat. Transf., 59, 453–469, https://doi.org/10.1016/S0022- 4073(97)00112-X, 1998. 33 Gordon, I. E., Rothman, L. S., Hill, C., Kochanov, R. V., Tan, Y., Bernath, P. F., Birk, M., Boudon, V., Campargue, A., Chance, K. V., Drouin, B. J., Flaud, J. M., Gamache, R. R., Hodges, J. T., Jacquemart, D., Perevalov, V. I., Perrin, A., Shine, K. P., Smith, M. A. H., Tennyson, J., Toon, G. C., Tran, H., Tyuterev, V. G., Barbe, A., Császár, A. G., Devi, V. M., Furtenbacher, T., Harrison, J. J., Hartmann, J. M., Jolly, A., Johnson, T. J., Karman, T., Kleiner, I., Kyuberis, A. A., Loos, J., Lyulin, O. M., Massie, S. T., Mikhailenko, S. N., Moazzen-Ahmadi, N., Müller, H. S. P., Naumenko, O. V., Nikitin, A. V., Polyansky, O. L., Rey, M., Rotger, M., Sharpe, S. W., Sung, K., Starikova, E., Tashkun, S. A., Auwera, J. V., Wagner, G., Wilzewski, J., Wcisło, P., Yu, S., and Zak, E. J.: The HITRAN2016 molecular spectroscopic database, J. Quant. Spectrosc. Radiat. Transf., 203, 3–69, https://doi.org/10.1016/j.jqsrt.2017.06.038, 2017. Hanuschik, R. W.: A ﬂux-calibrated, high-resolution atlas of optical sky emission from UVES, Astron. Astrophys., 407, 1157–1164, https://doi.org/10.1051/0004-6361:20030885, 2003. Hart, M.: Long-term Spectroscopic Observations of the Atmospheric Airglow by the Sloan Digital Sky Survey, Publ. Astron. Soc. Pac., 131, 015 003, https://doi.org/10.1088/1538-3873/aae972, 2019a. Hart, M.: A Comparison of Einstein A Coefﬁcients for OH Rotational Temperature Measurements Using a Large Astronomical Data Set, Atmosphere, 10, 569, https://doi.org/10.3390/atmos10100569, 2019b. Kalogerakis, K. S.: Technical note: Bimodality in mesospheric OH rotational population distributions and implications for temperature measurements, Atmos. Chem. Phys., 19, 2629–2634, https://doi.org/10.5194/acp-19-2629-2019, 2019. 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. Khomich, V. Y., Semenov, A. I., and Shefov, N. N.: Airglow as an Indicator of Upper Atmospheric Structure and Dynamics, Springer, Berlin, Langhoff, S. R., Werner, H.-J., and Rosmus, P.: Theoretical transition probabilities for the OH meinel system, J. Mol. Spectrosc., 118, 507–529, https://doi.org/10.1016/0022-2852(86)90186-4, 1986. Liu, G. and Shepherd, G. G.: An empirical model for the altitude of the OH nightglow emission, Geophys. Res. Lett., 33, L09805, https://doi.org/10.1029/2005GL025297, 2006. 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 coefﬁcients, J. Geophys. Res. Space Phys., 120, 10,069–10,082, https://doi.org/10.1002/2015JA021886, 2015. Llewellyn, E. J. and Long, B. H.: The OH Meinel bands in the airglow - The radiative lifetime, Can. J. Phys., 56, 581–586, https://doi.org/10.1139/p78-076, 1978. Meinel, A. B.: OH Emission Bands in the Spectrum of the Night Sky. I, Astrophys. J., 111, 555–564, https://doi.org/10.1086/145296, 1950. Melo, S. M. L., Lowe, R. P., and Takahashi, H.: The nocturnal behavior of the hydroxyl airglow at the equatorial and low latitudes as observed by WINDII: Comparison with ground-based measurements, J. Geophys. Res., 104, 24 657–24 666, https://doi.org/10.1029/1999JA900291, Mies, F. H.: Calculated vibrational transition probabilities of OH( X ), J. Molec. Spectrosc., 53, 150–188, https://doi.org/10.1016/0022- 2852(74)90125-8, 1974. 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 thermo- 34 sphere derived from SABER: Algorithm theoretical basis and measurement uncertainty, J. Geophys. Res. Atmos., 118, 5724–5735, https://doi.org/10.1002/jgrd.50401, 2013. Nelson, Jr., D. D., Schiffman, A., Nesbitt, D. J., Orlando, J. J., and Burkholder, J. B.: H + O Fourier-transform infrared emission and laser absorption studies of OH(X ) radical: An experimental dipole moment function and state-to-state Einstein A coefﬁcients, J. Chem. Phys., 93, 7003–7019, https://doi.org/10.1063/1.459476, 1990. Noll, S., Kausch, W., Barden, M., Jones, A. M., Szyszka, C., Kimeswenger, S., and Vinther, J.: An atmospheric radiation model for Cerro Paranal. I. The optical spectral range, Astron. Astrophys., 543, A92, https://doi.org/10.1051/0004-6361/201219040, 2012. 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. Noll, S., Kausch, W., Kimeswenger, S., Unterguggenberger, S., and Jones, A. M.: Comparison of VLT/X-shooter OH and O rota- tional temperatures with consideration of TIMED/SABER emission and temperature proﬁles, Atmos. Chem. Phys., 16, 5021–5042, https://doi.org/10.5194/acp-16-5021-2016, 2016. Noll, S., Kimeswenger, S., Proxauf, B., Unterguggenberger, S., Kausch, W., and Jones, A. M.: 15 years of VLT/UVES OH intensities and tem- peratures in comparison with TIMED/SABER data, J. Atmos. Sol.-Terr. Phys., 163, 54–69, https://doi.org/10.1016/j.jastp.2017.05.012, Noll, S., Proxauf, B., Kausch, W., and Kimeswenger, S.: Mechanisms for varying non-LTE contributions to OH rotational temperatures from measurements and modelling. I. Climatology, J. Atmos. Sol.-Terr. Phys., 175, 87–99, https://doi.org/10.1016/j.jastp.2018.05.004, 2018a. Noll, S., Proxauf, B., Kausch, W., and Kimeswenger, S.: Mechanisms for varying non-LTE contributions to OH rotational temperatures from measurements and modelling. II. Kinetic model, J. Atmos. Sol.-Terr. Phys., 175, 100–119, https://doi.org/10.1016/j.jastp.2018.05.005, 2018b. Noll, S., Plane, J. M. C., Feng, W., Proxauf, B., Kimeswenger, S., and Kausch, W.: Observations and modeling of potassium emission in the terrestrial nightglow, J. Geophys. Res. Atmos., 124, 6612–6629, https://doi.org/10.1029/2018JD030044, 2019. Oliva, E., Origlia, L., Maiolino, R., Baffa, C., Biliotti, V., Bruno, P., Falcini, G., Gavriousev, V., Ghinassi, F., Giani, E., Gonzalez, M., Leone, F., Lodi, M., Massi, F., Montegriffo, P., Mochi, I., Pedani, M., Rossetti, E., Scuderi, S., Sozzi, M., Tozzi, A., and Valenti, E.: A GIANO-TNG high-resolution infrared spectrum of the airglow emission, Astron. Astrophys., 555, A78, https://doi.org/10.1051/0004- 6361/201321366, 2013. 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. Parihar, N., Singh, D., and Gurubaran, S.: A comparison of ground-based hydroxyl airglow temperatures with SABER/TIMED measurements over 23 N, India, Ann. Geophys., 35, 353–363, https://doi.org/10.5194/angeo-35-353-2017, 2017. Pendleton, Jr., W., Espy, P., Baker, D., Steed, A., and Fetrow, M.: Observation of OH Meinel (7,4) P(N = 13) transitions in the night airglow, J. Geophys. Res., 94, 505–510, https://doi.org/10.1029/JA094iA01p00505, 1989. Pendleton, Jr., W. R. and Taylor, M. J.: The impact of L-uncoupling on Einstein coefﬁcients for the OH Meinel (6,2) band: implications for Q-branch rotational temperatures, J. Atmos. Sol.-Terr. Phys., 64, 971–983, https://doi.org/10.1016/S1364-6826(02)00051-2, 2002. Pendleton, Jr., W. R., Espy, P. J., and Hammond, M. R.: Evidence for non-local-thermodynamic-equilibrium rotation in the OH nightglow, J. Geophys. Res., 98, 11, https://doi.org/10.1029/93JA00740, 1993. 35 Picone, J. M., Hedin, A. E., Drob, D. P., and Aikin, A. C.: NRLMSISE-00 empirical model of the atmosphere: Statistical comparisons and scientiﬁc issues, J. Geophys. Res., 107, 1468, https://doi.org/10.1029/2002JA009430, 2002. Reisin, E. R., Scheer, J., Dyrland, M. E., Sigernes, F., Deehr, C. S., Schmidt, C., Höppner, K., Bittner, M., Ammosov, P. P., Gavrilyeva, G. A., Stegman, J., Perminov, V. I., Semenov, A. I., Knieling, P., Koppmann, R., Shiokawa, K., Lowe, R. P., López-González, M. J., Rodríguez, E., Zhao, Y., Taylor, M. J., Buriti, R. A., Espy, P. J., French, W. J. R., Eichmann, K.-U., Burrows, J. P., and von Savigny, C.: Traveling planetary wave activity from mesopause region airglow temperatures determined by the Network for the Detection of Mesospheric Change (NDMC), J. Atmos. Sol.-Terr. Phys., 119, 71–82, https://doi.org/10.1016/j.jastp.2014.07.002, 2014. Rothman, L. S., Gordon, I. E., Babikov, Y., Barbe, A., Chris Benner, D., Bernath, P. F., Birk, M., Bizzocchi, L., Boudon, V., Brown, L. R., Campargue, A., Chance, K., Cohen, E. A., Coudert, L. H., Devi, V. M., Drouin, B. J., Fayt, A., Flaud, J.-M., Gamache, R. R., Harrison, J. J., Hartmann, J.-M., Hill, C., Hodges, J. T., Jacquemart, D., Jolly, A., Lamouroux, J., Le Roy, R. J., Li, G., Long, D. A., Lyulin, O. M., Mackie, C. J., Massie, S. T., Mikhailenko, S., Müller, H. S. P., Naumenko, O. V., Nikitin, A. V., Orphal, J., Perevalov, V., Perrin, A., Polovtseva, E. R., Richard, C., Smith, M. A. H., Starikova, E., Sung, K., Tashkun, S., Tennyson, J., Toon, G. C., Tyuterev, V. G., and Wagner, G.: The HITRAN2012 molecular spectroscopic database, J. Quant. Spectrosc. Radiat. Transf., 130, 4–50, https://doi.org/10.1016/j.jqsrt.2013.07.002, 2013. Rousselot, P., Lidman, C., Cuby, J.-G., Moreels, G., and Monnet, G.: Night-sky spectral atlas of OH emission lines in the near-infrared, Astron. Astrophys., 354, 1134–1150, 2000. Russell, III, J. M., Mlynczak, M. G., Gordley, L. L., Tansock, J., and Esplin, R.: Overview of the SABER experiment and preliminary calibration results, in: Optical Spectroscopic Techniques and Instrumentation for Atmospheric and Space Research III, edited by Larar, A. M., vol. 3756 of SPIE Proc. Ser., pp. 277–288, https://doi.org/10.1117/12.366382, 1999. Schmidt, C., Höppner, K., and Bittner, M.: A ground-based spectrometer equipped with an InGaAs array for routine observations of OH(3- 1) rotational temperatures in the mesopause region, J. Atmos. Sol.-Terr. Phys., 102, 125–139, https://doi.org/10.1016/j.jastp.2013.05.001, Sedlak, R., Hannawald, P., Schmidt, C., Wüst, S., and Bittner, M.: High-resolution observations of small-scale gravity waves and turbulence features in the OH airglow layer, Atmos. Meas. Tech., 9, 5955–5963, https://doi.org/10.5194/amt-9-5955-2016, 2016. Tapping, K. F.: The 10.7 cm solar radio ﬂux (F ), Space Weather, 11, 394–406, https://doi.org/10.1002/swe.20064, 2013. 10:7 Taylor, M. J., Pendleton, W. R., Clark, S., Takahashi, H., Gobbi, D., and Goldberg, R. A.: Image measurements of short-period gravity waves at equatorial latitudes, J. Geophys. Res., 102, 26,283–26,299, https://doi.org/10.1029/96JD03515, 1997. Turnbull, D. N. and Lowe, R. P.: An empirical determination of the dipole moment function of OH(X ), J. Chem. Phys., 89, 2763–2767, https://doi.org/10.1063/1.455028, 1988. Turnbull, D. N. and Lowe, R. P.: New hydroxyl transition probabilities and their importance in airglow studies, Planet. Space Sci., 37, 723–738, https://doi.org/10.1016/0032-0633(89)90042-1, 1989. van der Loo, M. P. J. and Groenenboom, G. C.: Theoretical transition probabilities for the OH Meinel system, J. Chem. Phys., 126, 114 314– 114 314, https://doi.org/10.1063/1.2646859, 2007. 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, 159 902–159 902, https://doi.org/10.1063/1.2899016, 2008. van Rhijn, P. J.: On the brightness of the sky at night and the total amount of starlight, Publ. Kapteyn Astron. Lab. Groningen, 31, 1–83, 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. 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., 117, D02301, https://doi.org/10.1029/2011JD016342, 2012. 1 1 Yee, J.-H., Crowley, G., Roble, R. G., Skinner, W. R., Burrage, M. D., and Hays, P. B.: Global simulations and observations of O( S), O ( ) and OH mesospheric nightglow emissions, J. Geophys. Res., 102, 19 949–19 968, https://doi.org/10.1029/96JA01833, 1997.
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.pngPhysicsarXiv (Cornell University)http://www.deepdyve.com/lp/arxiv-cornell-university/oh-level-populations-and-accuracies-of-einstein-a-coefficients-from-gd1Hr600fO
OH level populations and accuracies of Einstein-A coefficients from hundreds of measured lines